Showing 1–50 of 230 results for author: Tan, J C

Asymmetric Kinematics in Young Clusters: The λ Ori Cluster

Authors: Joseph J. Armstrong, Jonathan C. Tan

Abstract: Context. Most stars form in clusters or associations but only a small number of these groups are expected to remain bound for longer than a few Myr. Once star formation has ended and the molecular gas around young stellar objects has been expelled via feedback processes, most initially bound young clusters lose the majority of their binding mass and begin to disperse into the Galactic field. Aims.… ▽ More Context. Most stars form in clusters or associations but only a small number of these groups are expected to remain bound for longer than a few Myr. Once star formation has ended and the molecular gas around young stellar objects has been expelled via feedback processes, most initially bound young clusters lose the majority of their binding mass and begin to disperse into the Galactic field. Aims. This process can be investigated by analysing the structure and kinematic trends in nearby young clusters, particularly expansion, the tell-tale sign that a cluster is no longer gravitationally bound but is dispersing into the field. Methods. We combine Gaia DR3 5-parameter astrometry with calibrated radial velocities for members of the nearby young cluster λ Ori (Collinder 69). Results. We characterise the plane-of-sky substructure of the cluster using the Q-parameter and Angular Dispersion parameter. We find evidence that the cluster contains significant substructure, but that this is preferentially located away from the central cluster core, which is smooth and likely remains bound. We find strong evidence for expansion in λ Ori in the plane-of-sky using a number of metrics, but also that the trends are asymmetric at the 5σ significance level. with the maximum rate of expansion being directed nearly parallel to the Galactic plane. We then invert the maximum rate of expansion of 0.144^{+0.003}_{-0.003} kms^{-1}pc^{-1} to give an expansion timescale of 6.944^{+0.148}_{-0.142} Myr, which is slightly larger than typical literature age estimates for the cluster. We also find asymmetry in the velocity dispersion, potential signatures of cluster rotation, and calculate kinematic ages for individual cluster members by tracing their motion back in time to their closest approach to the cluster center. △ Less

Submitted 16 July, 2024; originally announced July 2024.

Comments: 20 pages, 17 figures, submitted to A&A

The formation of supermassive black holes from Population III.1 seeds. III. Galaxy evolution and black hole growth from semi-analytic modelling

Authors: Vieri Cammelli, Pierluigi Monaco, Jonathan C. Tan, Jasbir Singh, Fabio Fontanot, Gabriella De Lucia, Michaela Hirschmann, Lizhi Xie

Abstract: We present an implementation of Pop III.1 seeding of supermassive black holes (SMBHs) in a theoretical model of galaxy formation and evolution to assess the growth the SMBH population and the properties of the host galaxies. The model of Pop III.1 seeding involves SMBH formation at redshifts $z\gtrsim 20$ in dark matter minihalos that are isolated from external radiative feedback, parameterized by… ▽ More We present an implementation of Pop III.1 seeding of supermassive black holes (SMBHs) in a theoretical model of galaxy formation and evolution to assess the growth the SMBH population and the properties of the host galaxies. The model of Pop III.1 seeding involves SMBH formation at redshifts $z\gtrsim 20$ in dark matter minihalos that are isolated from external radiative feedback, parameterized by isolation distance $d_{\rm iso}$. Within a standard $Λ$CDM cosmology, we generate dark matter halos using the code \textsc{pinocchio} and seed them according to the Pop III.1 scenario, exploring values of $d_{\rm iso}$ from 50 to 100~kpc (proper distance). We consider two alternative cases of SMBH seeding: a Halo Mass Threshold (HMT) model in which all halos $>7\times10^{10}\:M_\odot$ are seeded with $\sim 10^5\:M_\odot$ black holes; an All Light Seed (ALS) model in which all halos are seeded with low, stellar-mass black holes. We follow the redshift evolution of the halos, populating them with galaxies using the GAlaxy Evolution and Assembly theoretical model of galaxy formation, including accretion on SMBHs and related feedback processes. Here we present predictions for the properties of galaxy populations, focusing on stellar masses, star formation rates, and black hole masses. The local, $z\sim0$ metrics of occupation fraction as a function of the galaxy stellar mass, galaxy stellar mass function (GSMF), and black hole mass function (BHMF) all suggest a constraint of $d_{\rm iso}<75\:$kpc. We discuss the implications of this result for the Pop III.1 seeding mechanism. △ Less

Submitted 13 July, 2024; originally announced July 2024.

Comments: Submitted to MNRAS, comments welcome

Glimmers in the Cosmic Dawn: A Census of the Youngest Supermassive Black Holes by Photometric Variability

Authors: Matthew J. Hayes, Jonathan C. Tan, Richard S. Ellis, Alice R. Young, Vieri Cammelli, Jasbir Singh, Axel Runnholm, Aayush Saxena, Ragnhild Lunnan, Benjamin W. Keller, Pierluigi Monaco, Nicolas Laporte, Jens Melinder

Abstract: We report first results from a deep near infrared campaign with the Hubble Space Telescope to obtain late-epoch images of the Hubble Ultra-Deep Field (HUDF), 10-15 years after the first epoch data were obtained. The main objectives are to search for faint active galactic nuclei (AGN) at high redshifts by virtue of their photometric variability, and measure (or constrain) the comoving number densit… ▽ More We report first results from a deep near infrared campaign with the Hubble Space Telescope to obtain late-epoch images of the Hubble Ultra-Deep Field (HUDF), 10-15 years after the first epoch data were obtained. The main objectives are to search for faint active galactic nuclei (AGN) at high redshifts by virtue of their photometric variability, and measure (or constrain) the comoving number density of supermassive black holes (SMBHs), n_{SMBH}, at early times. In this Letter we present an overview of the program and preliminary results concerning eight objects. Three variables are supernovae, two of which are apparently hostless with indeterminable redshifts, although one has previously been recorded at a z\approx 6 object precisely because of its transient nature. Two further objects are clear AGN at z= 2.0 and 3.2, based on morphology and/or infrared spectroscopy from JWST. Three variable targets are identified at z = 6-7, which are also likely AGN candidates. These sources provide a first measure of n_{SMBH} in the reionization epoch by photometric variability, which places a firm lower limit of 3 \times 10^{-4} cMpc^{-3}. After accounting for variability and luminosity incompleteness, we estimate n_{SMBH} \gtrsim 8 \times 10{-3} cMpc{-3}, which is the largest value so far reported at these redshifts. This SMBH abundance is also strikingly similar to estimates of n_{SMBH} in the local Universe. We discuss how these results test various theories for SMBH formation. △ Less

Submitted 16 July, 2024; v1 submitted 24 March, 2024; originally announced March 2024.

Comments: In press at ApJL

The Core Mass Function Across Galactic Environments. IV. The Galactic Center

Authors: Alva V. I. Kinman, Maya A. Petkova, Jonathan C. Tan, Giuliana Cosentino, Yu Cheng

Abstract: The origin of the stellar Initial Mass Function (IMF) and how it may vary with galactic environment is a matter of debate. Certain star formation theories involve a close connection between the IMF and the Core Mass Function (CMF) and so it is important to measure this CMF in a range of locations in the Milky Way. Here we study the CMF of three Galactic Center clouds: G0.253+0.016 ("The Brick"), S… ▽ More The origin of the stellar Initial Mass Function (IMF) and how it may vary with galactic environment is a matter of debate. Certain star formation theories involve a close connection between the IMF and the Core Mass Function (CMF) and so it is important to measure this CMF in a range of locations in the Milky Way. Here we study the CMF of three Galactic Center clouds: G0.253+0.016 ("The Brick"), Sgr B2 (Deep South field) and Sgr C. We use ALMA 1 mm continuum images and identify cores as peaks in thermal dust emission via the dendrogram algorithm. We develop a completeness correction method via synthetic core insertion, where a realistic mass-dependent size distribution is used for the synthetic cores. After corrections, a power law of the form $\text{d}N/\text{d}\log M \propto M^{-α}$ is fit to the CMFs above 2 M$_\odot$. The three regions show disparate CMFs, with the Brick showing a Salpeter-like power law index $α=1.21\pm0.11$ and the other two regions showing shallower indices ($α=0.92\pm0.09$ for Sgr C and $α=0.66\pm0.05$ for Sgr B2-DS). Furthermore, we analyze the spatial distribution and mass segregation of cores in each region. Sgr C and Sgr B2-DS show signs of mass segregation, but the Brick does not. We compare our results to several other CMFs from different Galactic regions derived with the same methods. Finally, we discuss how these results may help define an evolutionary sequence of star cluster formation and can be used to test star formation theories. △ Less

Submitted 9 May, 2024; v1 submitted 6 March, 2024; originally announced March 2024.

Comments: Submitted to ApJ, comments welcome. 26 pages, 14 figures

Low-mass Runaways from the Orion Nebula Cluster -- Kinematic Age Constraints on Star Cluster Formation

Authors: Muhammad Fajrin, Joseph J. Armstrong, Jonathan C. Tan, Juan Farias, Laurent Eyer

Abstract: In their early, formative stages star clusters can undergo rapid dynamical evolution leading to strong gravitational interactions and ejection of ``runaway'' stars at high velocities. While O/B runaway stars have been well studied, lower-mass runaways are so far very poorly characterised, even though they are expected to be much more common. We carried out spectroscopic observations with MAG2-MIKE… ▽ More In their early, formative stages star clusters can undergo rapid dynamical evolution leading to strong gravitational interactions and ejection of ``runaway'' stars at high velocities. While O/B runaway stars have been well studied, lower-mass runaways are so far very poorly characterised, even though they are expected to be much more common. We carried out spectroscopic observations with MAG2-MIKE to follow-up 27 high priority candidate runaways consistent with having been ejected from the Orion Nebula Cluster (ONC) $>2.5$ Myr ago, based on Gaia astrometry. We derive spectroscopic youth indicators (Li \& H$α$) and radial velocities, enabling detection of bona fide runaway stars via signatures of youth and 3D traceback. We successfully confirmed 10 of the candidates as low-mass Young Stellar Objects (YSOs) on the basis of our spectroscopic criteria and derived radial velocities (RVs) with which we performed 3D traceback analysis. Three of these confirmed YSOs have kinematic ejection ages $>4\:$Myr, with the oldest being 4.7~Myr. This yields an estimate for the overall formation time of the ONC to be at least $\sim 5\:$Myr, i.e., about 10 free-fall times, and with a mean star formation efficiency per free-fall time of $\barε_{\rm ff}\lesssim0.05$. These results favor a scenario of slow, quasi-equilibrium star cluster formation, regulated by magnetic fields and/or protostellar outflow feedback. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: submitted to MNRAS, comments welcome

Polarized Light from Massive Protoclusters (POLIMAP). I. Dissecting the role of magnetic fields in the massive infrared dark cloud G28.37+0.07

Authors: C-Y Law, Jonathan C. Tan, Raphael Skalidis, Larry Morgan, Duo Xu, Felipe de Oliveira Alves, Ashley T. Barnes, Natalie Butterfield, Paola Caselli, Giuliana Cosentino, Francesco Fontani, Jonathan D. Henshaw, Izaskun Jimenez-Serra, Wanggi Lim

Abstract: Magnetic fields may play a crucial role in setting the initial conditions of massive star and star cluster formation. To investigate this, we report SOFIA-HAWC+ $214\:μ$m observations of polarized thermal dust emission and high-resolution GBT-Argus C$^{18}$O(1-0) observations toward the massive Infrared Dark Cloud (IRDC) G28.37+0.07. Considering the local dispersion of $B$-field orientations, we p… ▽ More Magnetic fields may play a crucial role in setting the initial conditions of massive star and star cluster formation. To investigate this, we report SOFIA-HAWC+ $214\:μ$m observations of polarized thermal dust emission and high-resolution GBT-Argus C$^{18}$O(1-0) observations toward the massive Infrared Dark Cloud (IRDC) G28.37+0.07. Considering the local dispersion of $B$-field orientations, we produce a map of $B$-field strength of the IRDC, which exhibits values between $\sim0.03 - 1\:$mG based on a refined Davis-Chandrasekhar-Fermi (r-DCF) method proposed by Skalidis \& Tassis. Comparing to a map of inferred density, the IRDC exhibits a $B-n$ relation with a power law index of $0.51\pm0.02$, which is consistent with a scenario of magnetically-regulated anisotropic collapse. Consideration of the mass-to-flux ratio map indicates that magnetic fields are dynamically important in most regions of the IRDC. A virial analysis of a sample of massive, dense cores in the IRDC, including evaluation of magnetic and kinetic internal and surface terms, indicates consistency with virial equilibrium, sub-Alfvénic conditions and a dominant role for $B-$fields in regulating collapse. A clear alignment of magnetic field morphology with direction of steepest column density gradient is also detected. However, there is no preferred orientation of protostellar outflow directions with the $B-$field. Overall, these results indicate that magnetic fields play a crucial role in regulating massive star and star cluster formation and so need to be accounted for in theoretical models of these processes. △ Less

Submitted 21 January, 2024; originally announced January 2024.

Comments: Submitted to ApJ, comments welcome

Near-Infrared Observations of Outflows and YSOs in the Massive Star-Forming Region AFGL 5180

Authors: S. Crowe, R. Fedriani, J. C. Tan, M. Whittle, Y. Zhang, A. Caratti o Garatti, J. P. Farias, A. Gautam, Z. Telkamp, B. Rothberg, M. Grudic, M. Andersen, G. Cosentino, R. Garcia-Lopez, V. Rosero, K. Tanaka, E. Pinna, F. Rossi, D. Miller, G. Agapito, C. Plantet, E. Ghose, J. Christou, J. Power, A. Puglisi , et al. (8 additional authors not shown)

Abstract: Methods: Broad- and narrow-band imaging of AFGL 5180 was made in the NIR with the LBT, in both seeing-limited ($\sim0.5\arcsec$) and high angular resolution ($\sim0.09\arcsec$) Adaptive Optics (AO) modes, as well as with HST. Archival ALMA continuum data was also utilized. Results: At least 40 jet knots were identified via NIR emission from H$_2$ and [FeII] tracing shocked gas. Bright jet knots… ▽ More Methods: Broad- and narrow-band imaging of AFGL 5180 was made in the NIR with the LBT, in both seeing-limited ($\sim0.5\arcsec$) and high angular resolution ($\sim0.09\arcsec$) Adaptive Optics (AO) modes, as well as with HST. Archival ALMA continuum data was also utilized. Results: At least 40 jet knots were identified via NIR emission from H$_2$ and [FeII] tracing shocked gas. Bright jet knots outflowing from the central most massive protostar, S4, are detected towards the east of the source and are resolved in fine detail with the AO imaging. Additional knots are distributed throughout the field, likely indicating the presence of multiple driving sources. Sub-millimeter sources detected by ALMA are shown to be grouped in two main complexes, AFGL 5180 M and a small cluster $\sim15\arcsec$ to the south, AFGL 5180 S. From our NIR continuum images we identify YSO candidates down to masses of $\sim 0.1\:M_\odot$. Combined with the sub-mm sources, this yields a surface number density of such YSOs of $N_* \sim 10^3 {\rm pc}^{-2}$ within a projected radius of about 0.1 pc. Such a value is similar to those predicted by models of both Core Accretion from a turbulent clump environment and Competitive Accretion. The radial profile of $N_*$ is relatively flat on scales out to 0.2~pc, with only modest enhancement around the massive protostar inside 0.05~pc. Conclusions: This study demonstrates the utility of high-resolution NIR imaging, in particular with AO, for detecting outflow activity and YSOs in distant regions. The presented images reveal the complex morphology of outflow-shocked gas within the large-scale bipolar flow of a massive protostar, as well as clear evidence for several other outflow driving sources in the region. Finally, this work presents a novel approach to compare the observed YSO surface number density from our study against different models of massive star formation. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: Accepted to Astronomy & Astrophysics (A&A)

The role of turbulence in high-mass star formation: Subsonic and transonic turbulence are ubiquitously found at early stages

Authors: Chao Wang, Ke Wang, Feng-Wei Xu, Patricio Sanhueza, Hauyu Baobab Liu, Qizhou Zhang, Xing Lu, F. Fontani, Paola Caselli, Gemma Busquet, Jonathan C. Tan, Di Li, J. M. Jackson, Thushara Pillai, Paul T. P. Ho, Andrés E. Guzmán, Nannan Yue

Abstract: Context. Traditionally, supersonic turbulence is considered to be one of the most likely mechanisms to slow down the gravitational collapse in dense clumps, thereby enabling the formation of massive stars. However, several recent studies have raised differing points of view based on observations carried out with sufficiently high spatial and spectral resolution. These studies call for a re-evaluat… ▽ More Context. Traditionally, supersonic turbulence is considered to be one of the most likely mechanisms to slow down the gravitational collapse in dense clumps, thereby enabling the formation of massive stars. However, several recent studies have raised differing points of view based on observations carried out with sufficiently high spatial and spectral resolution. These studies call for a re-evaluation of the role turbulence plays in massive star-forming regions. Aims. Our aim is to study the gas properties, especially the turbulence, in a sample of massive star-forming regions with sufficient spatial and spectral resolution, which can both resolve the core fragmentation and the thermal line width. Methods. We observed NH3 metastable lines with the Very Large Array (VLA) to assess the intrinsic turbulence. Results. Analysis of the turbulence distribution histogram for 32 identified NH3 cores reveals the presence of three distinct components. Furthermore, our results suggest that (1) sub- and transonic turbulence is a prevalent (21 of 32) feature of massive star-forming regions and those cold regions are at early evolutionary stage. This investigation indicates that turbulence alone is insufficient to provide the necessary internal pressure required for massive star formation, necessitating further exploration of alternative candidates; and (2) studies of seven multi-core systems indicate that the cores within each system mainly share similar gas properties and masses. However, two of the systems are characterized by the presence of exceptionally cold and dense cores that are situated at the spatial center of each system. Our findings support the hub-filament model as an explanation for this observed distribution △ Less

Submitted 7 February, 2024; v1 submitted 27 October, 2023; originally announced October 2023.

Comments: 34 pages, 15 figures, 4 tables. Accepted for publication on A&A

The JWST Galactic Center Survey -- A White Paper

Authors: Rainer Schoedel, Steve Longmore, Jonny Henshaw, Adam Ginsburg, John Bally, Anja Feldmeier, Matt Hosek, Francisco Nogueras Lara, Anna Ciurlo, Mélanie Chevance, J. M. Diederik Kruijssen, Ralf Klessen, Gabriele Ponti, Pau Amaro-Seoane, Konstantina Anastasopoulou, Jay Anderson, Maria Arias, Ashley T. Barnes, Cara Battersby, Giuseppe Bono, Lucía Bravo Ferres, Aaron Bryant, Miguel Cano Gonzáalez, Santi Cassisi, Leonardo Chaves-Velasquez , et al. (85 additional authors not shown)

Abstract: The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of… ▽ More The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of the most well-studied regions in astrophysics. Due to its proximity, we can study the center of our Galaxy on scales down to a few hundred AU, a hundred times better than in similar Local Group galaxies and thousands of times better than in the nearest active galaxies. The Galactic Center (GC) is therefore of outstanding astrophysical interest. However, in spite of intense observational work over the past decades, there are still fundamental things unknown about the GC. JWST has the unique capability to provide us with the necessary, game-changing data. In this White Paper, we advocate for a JWST NIRCam survey that aims at solving central questions, that we have identified as a community: i) the 3D structure and kinematics of gas and stars; ii) ancient star formation and its relation with the overall history of the Milky Way, as well as recent star formation and its implications for the overall energetics of our galaxy's nucleus; and iii) the (non-)universality of star formation and the stellar initial mass function. We advocate for a large-area, multi-epoch, multi-wavelength NIRCam survey of the inner 100\,pc of the Galaxy in the form of a Treasury GO JWST Large Program that is open to the community. We describe how this survey will derive the physical and kinematic properties of ~10,000,000 stars, how this will solve the key unknowns and provide a valuable resource for the community with long-lasting legacy value. △ Less

Submitted 14 March, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: This White Paper will be updated when required (e.g. new authors joining, editing of content). Most recent update: 24 Oct 2023

A HST Study of the Substellar Population of NGC 2024

Authors: Massimo Robberto, Mario Gennaro, Nicola Da Rio, Giovanni Maria Strampelli, Leonardo Ubeda, Elena Sabbi, Dana Koeppe, Jonathan C. Tan, David R. Soderblom

Abstract: We performed a HST/WFC3-IR imaging survey of the young stellar cluster NGC 2024 in three filters probing the 1.4~$μ$m H$_2$O absorption feature, characteristic of the population of low mass and sub-stellar mass objects down to a few Jupyter masses. We detect 812 point sources, 550 of them in all 3 filters with signal to noise greater than 5. Using a distance-independent two-color diagram we determ… ▽ More We performed a HST/WFC3-IR imaging survey of the young stellar cluster NGC 2024 in three filters probing the 1.4~$μ$m H$_2$O absorption feature, characteristic of the population of low mass and sub-stellar mass objects down to a few Jupyter masses. We detect 812 point sources, 550 of them in all 3 filters with signal to noise greater than 5. Using a distance-independent two-color diagram we determine extinction values as high as $A_V\simeq 40$. We also find that the change of effective wavelengths in our filters results in higher $A_V$ values as the reddening increases. Reconstructing a dereddened color-magnitude diagram we derive a luminosity histogram for both the full sample of candidate cluster members and for an extinction-limited sub-sample containing the 50% of sources with $A_V\lesssim 15$. Assuming a standard extinction law like Cardelli et al. (1989) with a nominal $R_V$=3.1 we produce a luminosity function in good agreement with the one resulting from a Salpeter-like Initial Mass Function for a 1~Myr isochrone. There is some evidence of an excess of luminous stars in the most embedded region. We posit that the correlation may be due to those sources being younger, and therefore overluminous than the more evolved and less extinct cluster's stars. We compare our classification scheme based on the depth of the 1.4$μ$m photometric feature with the results from the spectroscopic survey of Levine et al. (2006), and we report a few peculiar sources and morphological features typical of the rich phenomenology commonly encountered in young star-forming regions. △ Less

Submitted 7 October, 2023; originally announced October 2023.

Comments: 26 Pages, 23 Figures The Astrophysical Jornal, accepted

Disk Wind Feedback from High-mass Protostars. IV. Shock-Ionized Jets

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

Abstract: Massive protostars launch accretion-powered, magnetically-collimated outflows, which play crucial roles in the dynamics and diagnostics of the star formation process. Here we calculate the shock heating and resulting free-free radio emission in numerical models of outflows of massive star formation within the framework of the Turbulent Core Accretion model. We post-process 3D magneto-hydrodynamic… ▽ More Massive protostars launch accretion-powered, magnetically-collimated outflows, which play crucial roles in the dynamics and diagnostics of the star formation process. Here we calculate the shock heating and resulting free-free radio emission in numerical models of outflows of massive star formation within the framework of the Turbulent Core Accretion model. We post-process 3D magneto-hydrodynamic simulation snapshots of a protostellar disk wind interacts with an infalling core envelope, and calculate shock temperatures, ionization fractions, and radio free-free emission. We find heating up to ~10^7 K and near complete ionization in shocks at the interface between the outflow cavity and infalling envelope. However, line-of-sight averaged ionization fractions peak around ~10%, in agreement with values reported from observations of massive protostar G35.20-0.74N. By calculating radio continuum fluxes and spectra, we compare our models with observed samples of massive protostars. We find our fiducial models produce radio luminosities similar to those seen from low and intermediate-mass protostars that are thought to be powered by shock ionization. Comparing to more massive protostars, we find our model radio luminosities are ~10 to 100 times less luminous. We discuss how this apparent discrepancy either reflects aspects of our modeling related to the treatment of cooling of the post-shock gas or a dominant contribution in the observed systems from photoionization. Finally, our models exhibit 10-year radio flux variability of ~5%, especially in the inner 1000 au region, comparable to observed levels in some hyper-compact HII regions. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: Submitted to ApJ

Disk Wind Feedback from High-mass Protostars. III. Synthetic CO Line Emission

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

Abstract: To test theoretical models of massive star formation it is important to compare their predictions with observed systems. To this end, we conduct CO molecular line radiative transfer post-processing of 3D magneto-hydrodynamic (MHD) simulations of various stages in the evolutionary sequence of a massive protostellar core, including its infall envelope and disk wind outflow. Synthetic position-positi… ▽ More To test theoretical models of massive star formation it is important to compare their predictions with observed systems. To this end, we conduct CO molecular line radiative transfer post-processing of 3D magneto-hydrodynamic (MHD) simulations of various stages in the evolutionary sequence of a massive protostellar core, including its infall envelope and disk wind outflow. Synthetic position-position-velocity (PPV) cubes of various transitions of CO, 13CO, and C18O emission are generated. We also carry out simulated Atacama Large Millimeter/submillimeter Array (ALMA) observations of this emission. We compare the mass, momentum and kinetic energy estimates obtained from molecular lines to the true values, finding that the mass and momentum estimates can have uncertainties of up to a factor of four. However, the kinetic energy estimated from molecular lines is more significantly underestimated. Additionally, we compare the mass outflow rate and momentum outflow rate obtained from the synthetic spectra with the true values. Finally, we compare the synthetic spectra with real examples of ALMA-observed protostars and determine the best fitting protostellar masses and outflow inclination angles. We then calculate the mass outflow rate and momentum outflow rate for these sources, finding that both rates agree with theoretical protostellar evolutionary tracks. △ Less

Submitted 7 March, 2024; v1 submitted 7 September, 2023; originally announced September 2023.

Comments: ApJ Accepted

GMC Collisions As Triggers of Star Formation. IX. Chemical Evolution

Authors: Chia-Jung Hsu, Jonathan C. Tan, Jonathan Holdship, Duo, Xu, Serena Viti, Benjamin Wu, Brandt Gaches

Abstract: Collisions between giant molecular clouds (GMCs) have been proposed as a mechanism to trigger massive star and star cluster formation. To investigate the astrochemical signatures of such collisions, we carry out 3D magnetohydrodynamics simulations of colliding and non-colliding clouds exposed to a variety of cosmic ray ionization rates (CRIRs), $ζ$, following chemical evolution including gas and i… ▽ More Collisions between giant molecular clouds (GMCs) have been proposed as a mechanism to trigger massive star and star cluster formation. To investigate the astrochemical signatures of such collisions, we carry out 3D magnetohydrodynamics simulations of colliding and non-colliding clouds exposed to a variety of cosmic ray ionization rates (CRIRs), $ζ$, following chemical evolution including gas and ice-phase components. At the GMC scale, carbon starts mostly in $\rm{C^+}$, but then transitions into C, CO, followed by ice-phase CO and $\rm{CH_3OH}$ as dense, cooler filaments, clumps and cores form from the clouds. The oxygen budget is dominated by O, CO and water ice. In dense regions, we explore the gas phase CO depletion factor, $f_D$, that measures the extent of its freeze-out onto dust grains, including dependence on CRIR and observables of mass surface density and temperature. We also identify dense clumps and analyze their physical and chemical properties, including after synthetic line emission modeling, investigating metrics used in studies of infrared dark clouds (IRDCs), especially abundances of CO, $\rm HCO^+$ and $\rm N_2H^+$. For the colliding case, we find clumps have typical densities of $n_{\rm H}\sim10^5\:{\rm{cm}}^{-3}$ and temperatures of $\sim20\:$K, while those in non-colliding GMCs are cooler. Depending on $ζ$ and GMC dynamical history, we find CO depletion factors of up to $f_D\sim10$, and abundances of HCO$^+\sim 10^{-9}$ to $10^{-8}$ and $\rm{N_2H^+}\sim10^{-11}$ to $10^{-10}$. Comparison with observed IRDC clumps indicates a preference for low CRIRs ($\sim10^{-18}\:{\rm{s}}^{-1}$) and a more quiescent (non-colliding), cooler and evolved chemodynamical history. We discuss the general implications of our results and their caveats for interpretation of molecular cloud observations. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 32 pages, 30 figures, 2 tables, submitted to MNRAS, comments welcome

On the Lifetime of Molecular Clouds with the "Tuning-Fork" Analysis

Authors: Jin Koda, Jonathan C. Tan

Abstract: The "tuning-fork" (TF) analysis of CO and Halpha emission has been used to estimate the lifetimes of molecular clouds in nearby galaxies. With simple model calculations, we show that this analysis does not necessarily estimate cloud lifetimes, but instead captures a duration of the cloud evolutionary cycle, from dormant to star forming, and then back to a dormant phase. We adopt a hypothetical set… ▽ More The "tuning-fork" (TF) analysis of CO and Halpha emission has been used to estimate the lifetimes of molecular clouds in nearby galaxies. With simple model calculations, we show that this analysis does not necessarily estimate cloud lifetimes, but instead captures a duration of the cloud evolutionary cycle, from dormant to star forming, and then back to a dormant phase. We adopt a hypothetical setup in which molecular clouds (e.g., traced in CO) live forever and form stars (e.g., HII regions) at some frequency, which then drift away from the clouds. The TF analysis still returns a timescale for the immortal clouds. This model requires drifting motion to separate the newborn stars from the clouds, and we discuss its origin. We also discuss the physical origin of the characteristic spatial separation term in the TF analysis and a bias due to systematic error in the determination of the reference timescale. △ Less

Submitted 29 October, 2023; v1 submitted 22 August, 2023; originally announced August 2023.

Comments: ApJ accepted

Spitzer thermal phase curve of WASP-121 b

Authors: Giuseppe Morello, Quentin Changeat, Achrène Dyrek, Pierre-Olivier Lagage, Jonathan C. Tan

Abstract: Aims. We analyse unpublished Spitzer observations of the thermal phase-curve of WASP-121 b, a benchmark ultra-hot Jupiter. Methods. We adopted the wavelet pixel-independent component analysis technique to remove challenging instrumental systematic effects in these datasets and we fit them simultaneously with parametric light-curve models. We also performed phase-curve retrievals to better understa… ▽ More Aims. We analyse unpublished Spitzer observations of the thermal phase-curve of WASP-121 b, a benchmark ultra-hot Jupiter. Methods. We adopted the wavelet pixel-independent component analysis technique to remove challenging instrumental systematic effects in these datasets and we fit them simultaneously with parametric light-curve models. We also performed phase-curve retrievals to better understand the horizontal and vertical thermal structure of the planetary atmosphere. Results. We measured planetary brightness temperatures of $\sim$2700\,K (dayside) and $\sim$700--1100\,K (nightside), along with modest peak offsets of 5.9$^{\circ} \pm$1.6 (3.6\,$μ$m) and 5.0$^{\circ}$$_{-3.1}^{+3.4}$ (4.5\,$μ$m) after mid-eclipse. These results suggest inefficient heat redistribution in the atmosphere of WASP-121 b. The inferred atmospheric Bond albedo and circulation efficiency align well with observed trends for hot giant exoplanets. Interestingly, the measured peak offsets correspond to a westward hot spot, which has rarely been observed. We also report consistent transit depths at 3.6 and 4.5\,$μ$m, along with updated geometric and orbital parameters. Finally, we compared our Spitzer results with previous measurements, including recent JWST observations. Conclusions. We extracted new information on the thermal properties and dynamics of an exoplanet atmosphere from an especially problematic dataset. This study probes the reliability of exoplanet phase-curve parameters obtained from Spitzer observations when state-of-the-art pipelines are adopted to remove the instrumental systematic effects. It demonstrates that Spitzer phase-curve observations provide a useful baseline for comparison with JWST observations, and shows the increase in parameters precision achieved with the newer telescope. △ Less

Submitted 9 August, 2023; v1 submitted 2 July, 2023; originally announced July 2023.

Comments: 14 pages, 10 figures

Journal ref: A&A 676, A54 (2023)

The sharpest view on the high-mass star-forming region S255IR. Near-InfraRed Adaptive Optics Imaging on the Outbursting Source NIRS3

Authors: R. Fedriani, A. Caratti o Garatti, R. Cesaroni, J. C. Tan, B. Stecklum, L. Moscadelli, M. Koutoulaki, G. Cosentino, M. Whittle

Abstract: Massive stars have an impact on their surroundings from early in their formation until the end of their lives. However, very little is known about their formation. Episodic accretion may play a crucial role, but observations of these events have only been reported towards a handful of massive protostars. We aim to investigate the outburst event from the high-mass star-forming region S255IR where r… ▽ More Massive stars have an impact on their surroundings from early in their formation until the end of their lives. However, very little is known about their formation. Episodic accretion may play a crucial role, but observations of these events have only been reported towards a handful of massive protostars. We aim to investigate the outburst event from the high-mass star-forming region S255IR where recently the protostar NIRS3 underwent an accretion outburst. We follow the evolution of this source both in photometry and morphology of its surroundings. Methods: We perform near-infrared adaptive optics observations on the S255IR central region using the Large Binocular Telescope in the K$_{\rm s}$ broad-band and the H$_2$ and Br$γ$ narrow-band filters with an angular resolution of $\sim0\farcs06$, close to the diffraction limit. We discover a new near-infrared knot north-east from NIRS3 that we interpret as a jet knot that was ejected during the last accretion outburst and observed in the radio regime as part of a follow-up after the outburst. We measure a mean tangential velocity for this knot of $450\pm50\,\mathrm{km\,s^{-1}}$. We analyse the continuum-subtracted images from H$_2$ which traces jet shocked emission, and Br$γ$ which traces scattered light from a combination of accretion activity and UV radiation from the central massive protostar. We observe a significant decrease in flux at the location of NIRS3, with K=13.48\,mag being the absolute minimum in the historic series. Our observations strongly suggest a scenario where the episodic accretion is followed by an episodic ejection response in the near-infrared, as it was seen in the earlier radio follow-up. The 30 years of $\sim2\,μ{\rm m}$ photometry suggests that NIRS3 might have undergone another outburst in the late 1980s, being the first massive protostar with such evidence observed in the near-infrared. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: Accepted for publication in Astronomy and Astrophysics. 10 pages, 8 figures

Journal ref: A&A 676, A107 (2023)

Deuterium Fractionation across the Infrared Dark Cloud G034.77-00.55 interacting with the Supernova Remnant W44

Authors: G. Cosentino, J. C. Tan, I. Jiménez-Serra, F. Fontani, P. Caselli, J. D. Henshaw, A. T. Barnes, C. -Y. Law, S. Viti, R. Fedriani, C. -J. Hsu, P. Gorai, S. Zeng

Abstract: Supernova remnants (SNRs) may regulate star formation in galaxies. For example, SNR-driven shocks may form new molecular gas or compress pre-existing clouds and trigger the formation of new stars. To test this scenario, we measure the deuteration of $N_2H^+$, $D_{frac}^{N_2H^+}$, a well-studied tracer of pre-stellar cores, across the Infrared Dark Cloud (IRDC) G034.77-00.55, known to be experienci… ▽ More Supernova remnants (SNRs) may regulate star formation in galaxies. For example, SNR-driven shocks may form new molecular gas or compress pre-existing clouds and trigger the formation of new stars. To test this scenario, we measure the deuteration of $N_2H^+$, $D_{frac}^{N_2H^+}$, a well-studied tracer of pre-stellar cores, across the Infrared Dark Cloud (IRDC) G034.77-00.55, known to be experiencing a shock interaction with the SNR W44. We use N$_2$H$^+$ and N$_2$D$^+$ J=1-0 single pointing observations obtained with the 30m antenna at the Instituto de Radioastronomia Millimetrica to infer $D_{frac}^{N_2H^+}$ toward five positions across the cloud, namely a massive core, different regions across the shock front, a dense clump and ambient gas. We find $D_{frac}^{N_2H^+}$ in the range 0.03-0.1, several orders of magnitude larger than the cosmic D/H ratio ($\sim$10$^{-5}$). Across the shock front, $D_{frac}^{N_2H^+}$ is enhanced by more than a factor of 2 ($D_{frac}^{N_2H^+}\sim$0.05-0.07) with respect to the ambient gas ($\leq$0.03) and similar to that measured generally in pre-stellar cores. Indeed, in the massive core and dense clump regions of this IRDC we measure $D_{frac}^{N_2H^+}$}$\sim$0.1. We find enhanced deuteration of $N_2H^+$ across the region of the shock, at a level that is enhanced with respect to regions of unperturbed gas. It is possible that this has been induced by shock compression, which would then be indirect evidence that the shock is triggering conditions for future star formation. However, since unperturbed dense regions also show elevated levels of deuteration, further, higher-resolution studies are needed to better understand the structure and kinematics of the deuterated material in the shock region, e.g., if it still in relatively diffuse form or already organised in a population of low-mass pre-stellar cores. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: Accepted for publication on A&A; 8 pages, 5 figures

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

Chemical Differentiation around Five Massive Protostars Revealed by ALMA -Carbon-Chain Species, Oxygen-/Nitrogen-Bearing Complex Organic Molecules-

Authors: Kotomi Taniguchi, Liton Majumdar, Paola Caselli, Shigehisa Takakuwa, Tien-Hao Hsieh, Masao Saito, Zhi-Yun Li, Kazuhito Dobashi, Tomomi Shimoikura, Fumitaka Nakamura, Jonathan C. Tan, Eric Herbst

Abstract: We present Atacama Large Millimeter/submillimeter Array Band 3 data toward five massive young stellar objects (MYSOs), and investigate relationships between unsaturated carbon-chain species and saturated complex organic molecules (COMs). An HC$_{5}$N ($J=35-34$) line has been detected from three MYSOs, where nitrogen(N)-bearing COMs (CH$_{2}$CHCN and CH$_{3}$CH$_{2}$CN) have been detected. The HC… ▽ More We present Atacama Large Millimeter/submillimeter Array Band 3 data toward five massive young stellar objects (MYSOs), and investigate relationships between unsaturated carbon-chain species and saturated complex organic molecules (COMs). An HC$_{5}$N ($J=35-34$) line has been detected from three MYSOs, where nitrogen(N)-bearing COMs (CH$_{2}$CHCN and CH$_{3}$CH$_{2}$CN) have been detected. The HC$_{5}$N spatial distributions show compact features and match with a methanol (CH$_{3}$OH) line with an upper-state energy around 300 K, which should trace hot cores. The hot regions are more extended around the MYSOs where N-bearing COMs and HC$_{5}$N have been detected compared to two MYSOs without these molecular lines, while there are no clear differences in the bolometric luminosity and temperature. We run chemical simulations of hot-core models with a warm-up stage, and compare with the observational results. The observed abundances of HC$_{5}$N and COMs show good agreements with the model at the hot-core stage with temperatures above 160 K. These results indicate that carbon-chain chemistry around the MYSOs cannot be reproduced by warm carbon-chain chemistry, and a new type of carbon-chain chemistry occurs in hot regions around MYSOs. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: Accepted by the publication for The Astrophysical Journal Supplement Series, 32 pages,18 figures, 11 tables

Denoising Diffusion Probabilistic Models to Predict the Density of Molecular Clouds

Authors: Duo Xu, Jonathan C. Tan, Chia-Jung Hsu, Ye Zhu

Abstract: We introduce the state-of-the-art deep learning Denoising Diffusion Probabilistic Model (DDPM) as a method to infer the volume or number density of giant molecular clouds (GMCs) from projected mass surface density maps. We adopt magnetohydrodynamic simulations with different global magnetic field strengths and large-scale dynamics, i.e., noncolliding and colliding GMCs. We train a diffusion model… ▽ More We introduce the state-of-the-art deep learning Denoising Diffusion Probabilistic Model (DDPM) as a method to infer the volume or number density of giant molecular clouds (GMCs) from projected mass surface density maps. We adopt magnetohydrodynamic simulations with different global magnetic field strengths and large-scale dynamics, i.e., noncolliding and colliding GMCs. We train a diffusion model on both mass surface density maps and their corresponding mass-weighted number density maps from different viewing angles for all the simulations. We compare the diffusion model performance with a more traditional empirical two-component and three-component power-law fitting method and with a more traditional neural network machine learning approach (CASI-2D). We conclude that the diffusion model achieves an order of magnitude improvement on the accuracy of predicting number density compared to that by other methods. We apply the diffusion method to some example astronomical column density maps of Taurus and the Infrared Dark Clouds (IRDCs) G28.37+0.07 and G35.39-0.33 to produce maps of their mean volume densities. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: ApJ accepted

Carbon-Chain Chemistry in the Interstellar Medium

Authors: Kotomi Taniguchi, Prasanta Gorai, Jonathan C. Tan

Abstract: The presence of carbon-chain molecules in the interstellar medium (ISM) has been known since the early 1970s and $>130$ such species have been identified to date, making up $\sim 43$% of the total of detected ISM molecules. They are prevalent not only in star-forming regions in our Galaxy but also in other galaxies. These molecules provide important information on physical conditions, gas dynamics… ▽ More The presence of carbon-chain molecules in the interstellar medium (ISM) has been known since the early 1970s and $>130$ such species have been identified to date, making up $\sim 43$% of the total of detected ISM molecules. They are prevalent not only in star-forming regions in our Galaxy but also in other galaxies. These molecules provide important information on physical conditions, gas dynamics, and evolutionary stages of star-forming regions. Larger species of polycyclic aromatic hydrocarbons (PAHs) and fullerenes (C$_{60}$ and C$_{70}$), which may be related to the formation of the carbon-chain molecules, have been detected in circumstellar envelopes around carbon-rich Asymptotic Giant Branch (AGB) stars and planetary nebulae, while PAHs are also known to be a widespread component of the ISM in most galaxies. Recently, two line survey projects toward Taurus Molecular Cloud-1 with large single-dish telescopes have detected many new carbon-chain species, including molecules containing benzene rings. These new findings raise fresh questions about carbon-bearing species in the Universe. This article reviews various aspects of carbon-chain molecules, including observational studies, chemical simulations, quantum calculations, and laboratory experiments, and discusses open questions and how future facilities may answer them. △ Less

Submitted 12 April, 2024; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: A review article in Astrophysics and Space Science (https://link.springer.com/article/10.1007/s10509-024-04292-9)

Journal ref: Astrophysics and Space Science, 369, 34 (2024)

Mother of Dragons: A Massive, quiescent core in the dragon cloud (IRDC G028.37+00.07)

Authors: A. T. Barnes, J. Liu, Q. Zhang, J. C. Tan, F. Bigiel, P. Caselli, G. Cosentino, F. Fontani, J. D. Henshaw, I. Jiménez-Serra, D-S. Kalb, C. Y. Law, S. N. Longmore, R. J. Parker, J. E. Pineda, A. Sánchez-Monge, W. Lim, K. Wang

Abstract: Context: Core accretion models of massive star formation require the existence of massive, starless cores within molecular clouds. Yet, only a small number of candidates for such truly massive, monolithic cores are currently known. Aims: Here we analyse a massive core in the well-studied infrared-dark cloud (IRDC) called the 'dragon cloud' (also known as G028.37+00.07 or 'Cloud C'). This core (C2c… ▽ More Context: Core accretion models of massive star formation require the existence of massive, starless cores within molecular clouds. Yet, only a small number of candidates for such truly massive, monolithic cores are currently known. Aims: Here we analyse a massive core in the well-studied infrared-dark cloud (IRDC) called the 'dragon cloud' (also known as G028.37+00.07 or 'Cloud C'). This core (C2c1) sits at the end of a chain of a roughly equally spaced actively star-forming cores near the centre of the IRDC. Methods: We present new high-angular resolution 1 mm ALMA dust continuum and molecular line observations of the massive core. Results: The high-angular resolution observations show that this region fragments into two cores C2c1a and C2c1b, which retain significant background-subtracted masses of 23 Msun and 2 Msun (31 Msun and 6 Msun without background subtraction), respectively. The cores do not appear to fragment further on the scales of our highest angular resolution images (0.200 arcsec, 0.005 pc ~ 1000 AU). We find that these cores are very dense (nH2 > 10^6 cm-3) and have only trans-sonic non-thermal motions (Ms ~ 1). Together the mass, density and internal motions imply a virial parameter of < 1, which suggests the cores are gravitationally unstable, unless supported by strong magnetic fields with strengths of ~ 1 - 10 mG. From CO line observations, we find that there is tentative evidence for a weak molecular outflow towards the lower-mass core, and yet the more massive core remains devoid of any star formation indicators. Conclusions: We present evidence for the existence of a massive, pre-stellar core, which has implications for theories of massive star formation. This source warrants follow-up higher-angular-resolution observations to further assess its monolithic and pre-stellar nature. △ Less

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

Comments: 8+4 pages, 4+2 Figures, 2 Tables, accepted for publication in Astronomy & Astrophysics

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

The formation of supermassive black holes from Population III.1 seeds. II. Evolution to the local universe

Authors: Jasbir Singh, Pierluigi Monaco, Jonathan C. Tan

Abstract: We present predictions for cosmic evolution of populations of supermassive black holes (SMBHs) forming from Population III.1 seeds, i.e., early, metal-free dark matter minihalos forming far from other sources, parameterized by isolation distance, $d_{\rm{iso}}$. Extending previous work that explored this scenario to $z=10$, we follow evolution of a $(60\:{\rm{Mpc}})^3$ volume to $z=0$. We focus on… ▽ More We present predictions for cosmic evolution of populations of supermassive black holes (SMBHs) forming from Population III.1 seeds, i.e., early, metal-free dark matter minihalos forming far from other sources, parameterized by isolation distance, $d_{\rm{iso}}$. Extending previous work that explored this scenario to $z=10$, we follow evolution of a $(60\:{\rm{Mpc}})^3$ volume to $z=0$. We focus on evolution of SMBH comoving number densities, halo occupation fractions, angular clustering and 3D clustering, exploring a range of $d_{\rm{iso}}$ constrained by observed local number densities of SMBHs. We also compute synthetic projected observational fields, in particular a case comparable to the Hubble Ultra Deep Field. We compare Pop III.1 seeding to a simple halo mass threshold model, commonly adopted in cosmological simulations of galaxy formation. Major predictions of the Pop III.1 model include that all SMBHs form by $z\sim25$, after which their comoving number densities are near-constant, with low merger rates. Occupation fractions evolve to concentrate SMBHs in the most massive halos by $z=0$, but with rare cases in halos down to $\sim10^8\:M_\odot$. The $d_{\rm{iso}}$ scale at epoch of formation, e.g., $100\:$kpc-proper at $z\sim30$, i.e., $\sim3\:$Mpc-comoving, is imprinted in the SMBH two-point angular correlation function, remaining discernible as a low-amplitude feature to $z\sim1$. The SMBH 3D two-point correlation function at $z=0$ also shows lower amplitude compared to equivalently massive halos. We discuss prospects for testing these predictions with observational surveys of SMBH populations. △ Less

Submitted 31 July, 2023; v1 submitted 26 January, 2023; originally announced January 2023.

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

GMC Collisions As Triggers of Star Formation. VIII. The Core Mass Function

Authors: Chia-Jung Hsu, Jonathan C. Tan, Duncan Christie, Yu Cheng, Theo J. O'Neill

Abstract: Compression in giant molecular cloud (GMC) collisions is a promising mechanism to trigger formation of massive star clusters and OB associations. We simulate colliding and non-colliding magnetised GMCs and examine the properties of prestellar cores, selected from projected mass surface density maps, including after synthetic {\it ALMA} observations. We then examine core properties, including mass,… ▽ More Compression in giant molecular cloud (GMC) collisions is a promising mechanism to trigger formation of massive star clusters and OB associations. We simulate colliding and non-colliding magnetised GMCs and examine the properties of prestellar cores, selected from projected mass surface density maps, including after synthetic {\it ALMA} observations. We then examine core properties, including mass, size, density, velocity, velocity dispersion, temperature and magnetic field strength. After four Myr, $\sim1,000$ cores have formed in the GMC collision and the high-mass end of the core mass function (CMF) can be fit by a power law $dN/d{\rm{log}}M\propto{M}^{-α}$ with $α\simeq0.7$, i.e., relatively top-heavy compared to a Salpeter mass function. Depending on how cores are identified, a break in the power law can appear around a few $\times10\:M_\odot$. The non-colliding GMCs form fewer cores with a CMF with $α\simeq0.8$ to 1.2, i.e., closer to the Salpeter index. We compare the properties of these CMFs to those of several observed samples of cores. Considering other properties, cores formed from colliding clouds are typically warmer, have more disturbed internal kinematics and are more likely to be gravitational unbound, than cores formed from non-colliding GMCs. The dynamical state of the protocluster of cores formed in the GMC-GMC collision is intrinsically subvirial, but can appear to be supervirial if the total mass measurement is affected by observations that miss mass on large scales or at low densities. △ Less

Submitted 3 March, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: 22 pages, 22 figures, 2 tables, accepted to MNRAS

Star Cluster Formation from Turbulent Clumps. III. Across the mass spectrum

Authors: Juan P. Farias, Jonathan C. Tan

Abstract: We study the formation and early evolution of star clusters that have a wide range of masses and background cloud mass surface densities, $Σ_{\rm cloud}$, which help set the initial sizes, densities, and velocity dispersions of the natal gas clumps. Initial clump masses of 300, $3,000$ and $30,000$ $M_\odot$ are considered, from which star clusters are born with an assumed 50% overall star formati… ▽ More We study the formation and early evolution of star clusters that have a wide range of masses and background cloud mass surface densities, $Σ_{\rm cloud}$, which help set the initial sizes, densities, and velocity dispersions of the natal gas clumps. Initial clump masses of 300, $3,000$ and $30,000$ $M_\odot$ are considered, from which star clusters are born with an assumed 50% overall star formation efficiency and with 50% primordial binarity. This formation is gradual, i.e., with a range of star formation efficiencies per free-fall time from 1% to 100%, so that the formation time can range from 0.7 Myr for low-mass, high-$Σ_{\rm cloud}$ clumps to $\sim30$ Myr for high-mass, low-$Σ_{\rm cloud}$ clumps. Within this framework of the Turbulent Clump model, for a given $Σ_{\rm cloud}$, clumps of higher mass are of lower initial volume density, but their dynamical evolution leads to higher bound fractions and causes them to form much higher density cluster cores and maintain these densities for longer periods. This results in systematic differences in the evolution of binary properties, degrees of mass segregation and rates of creation of dynamically ejected runaways. We discuss the implications of these results for observed star clusters and stellar populations. △ Less

Submitted 7 June, 2023; v1 submitted 21 January, 2023; originally announced January 2023.

Comments: 28 pages, 20 Figures. Published by MNRAS

What Sets the Star Formation Rate of Molecular Clouds? The Density Distribution as a Fingerprint of Compression and Expansion Rates

Authors: Sabrina M. Appel, Blakesley Burkhart, Vadim A. Semenov, Christoph Federrath, Anna L. Rosen, Jonathan C. Tan

Abstract: We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycl… ▽ More We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycling between high and low density states. We find that the inclusion of protostellar jets produces rapidly expanding and compressing low-density gas. We find that the gas mass flux peaks at the transition between the lognormal and power-law forms of the density probability distribution function (PDF). This is consistent with the transition density tracking the post-shock density, which promotes an enhancement of mass at this density (i.e., shock compression and filament formation). At high densities, the gas dynamics are dominated by self-gravity: the compression rate in all of our runs matches the rate of the run with only gravity, suggesting that processes other than self-gravity have little effect at these densities. The net gas mass flux becomes constant at a density below the sink formation threshold, where it equals the star formation rate. The density at which the net gas mass flux equals the star formation rate is one order of magnitude lower than our sink threshold density, corresponds to the formation of the second power-law tail in the density PDF, and sets the overall star formation rates of these simulations. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 22 pages, 13 figures. Submitted to ApJ. Comments welcome

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

Application of Convolutional Neural Networks to Predict Magnetic Fields Directions in Turbulent Clouds

Authors: Duo Xu, Chi-Yan Law, Jonathan C. Tan

Abstract: We adopt the deep learning method CASI-3D (Convolutional Approach to Structure Identification-3D) to infer the orientation of magnetic fields in sub-/trans- Alfvenic turbulent clouds from molecular line emission. We carry out magnetohydrodynamic simulations with different magnetic field strengths and use these to generate synthetic observations. We apply the 3D radiation transfer code RADMC-3d to… ▽ More We adopt the deep learning method CASI-3D (Convolutional Approach to Structure Identification-3D) to infer the orientation of magnetic fields in sub-/trans- Alfvenic turbulent clouds from molecular line emission. We carry out magnetohydrodynamic simulations with different magnetic field strengths and use these to generate synthetic observations. We apply the 3D radiation transfer code RADMC-3d to model 12CO and 13CO (J = 1-0) line emission from the simulated clouds and then train a CASI-3D model on these line emission data cubes to predict magnetic field morphology at the pixel level. The trained CASI-3D model is able to infer magnetic field directions with low error (< 10deg for sub-Alfvenic samples and <30deg for trans-Alfvenic samples). We furthermore test the performance of CASI-3D on a real sub-/trans- Alfvenic region in Taurus. The CASI-3D prediction is consistent with the magnetic field direction inferred from Planck dust polarization measurements. We use our developed methods to produce a new magnetic field map of Taurus that has a three-times higher angular resolution than the Planck map. △ Less

Submitted 25 November, 2022; originally announced November 2022.

Comments: ApJ Accepted

A Census of Outflow to Magnetic Field Orientations in Nearby Molecular Clouds

Authors: Duo Xu, Stella S. R. Offner, Robert Gutermuth, Jonathan C. Tan

Abstract: We define a sample of 200 protostellar outflows showing blue and redshifted CO emission in the nearby molecular clouds Ophiuchus, Taurus, Perseus and Orion to investigate the correlation between outflow orientations and local, but relatively large-scale, magnetic field directions traced by Planck 353 GHz dust polarization. At high significance (p~1e-4), we exclude a random distribution of relative… ▽ More We define a sample of 200 protostellar outflows showing blue and redshifted CO emission in the nearby molecular clouds Ophiuchus, Taurus, Perseus and Orion to investigate the correlation between outflow orientations and local, but relatively large-scale, magnetic field directions traced by Planck 353 GHz dust polarization. At high significance (p~1e-4), we exclude a random distribution of relative orientations and find that there is a preference for alignment of projected plane of sky outflow axes with magnetic field directions. The distribution of relative position angles peaks at ~30deg and exhibits a broad dispersion of ~50deg. These results indicate that magnetic fields have dynamical influence in regulating the launching and/or propagation directions of outflows. However, the significant dispersion around perfect alignment orientation implies that there are large measurement uncertainties and/or a high degree of intrinsic variation caused by other physical processes, such as turbulence or strong stellar dynamical interactions. Outflow to magnetic field alignment is expected to lead to a correlation in the directions of nearby outflow pairs, depending on the degree of order of the field. Analyzing this effect we find limited correlation, except on relatively small scales < 0.5 pc. Furthermore, we train a convolutional neural network to infer the inclination angle of outflows with respect to the line of sight and apply it to our outflow sample to estimate their full 3D orientations. We find that the angles between outflow pairs in 3D space also show evidence of small-scale alignment. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Comments: ApJ Accepted

The Disk Population in a Distant Massive Protocluster

Authors: Yu Cheng, Jonathan C. Tan, John J. Tobin, Ruben Fedriani, Morten Andersen, Junfeng Wang

Abstract: The unprecedented angular resolution and sensitivity of ALMA makes it possible to unveil disk populations in distant ($>$2 kpc), embedded young cluster environments. We have conducted an observation towards the central region of the massive protocluster G286.21+0.16 at 1.3 mm. With a spatial resolution of 23 mas and a sensitivity of 15 $\rm μJy~beam^{-1}$, we detect a total of 38 protostellar disk… ▽ More The unprecedented angular resolution and sensitivity of ALMA makes it possible to unveil disk populations in distant ($>$2 kpc), embedded young cluster environments. We have conducted an observation towards the central region of the massive protocluster G286.21+0.16 at 1.3 mm. With a spatial resolution of 23 mas and a sensitivity of 15 $\rm μJy~beam^{-1}$, we detect a total of 38 protostellar disks. These disks have dust masses ranging from about 53 to 1825 $M_\oplus$, assuming a dust temperature of 20 K. This sample is not closely associated with previously identified dense cores, as would be expected for disks around Class 0 protostars. Thus, we expect our sample, being flux limited, to be mainly composed of Class I/flat-spectrum source disks, since these are typically more massive than Class II disks. Furthermore, we find that the distributions of disk masses and radii are statistically indistinguishable with those of the Class I/flat-spectrum objects in the Orion molecular cloud, indicating similar processes are operating in G286.21+0.16 to regulate disk formation and evolution. The cluster center appears to host a massive protostellar system composed of three sources within 1200 au, including a potential binary with 600 au projected separation. Relative to this center, there is no evidence for widespread mass segregation in the disk population. We do find a tentative trend of increasing disk radius versus distance from the cluster center, which may point to the influence of dynamical interactions being stronger in the central regions. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 23 pages, 7 Figures, accepted by ApJ

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 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

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

Inside-Out Planet Formation. VII. Astrochemical Models of Protoplanetary Disks and Implications for Planetary Compositions

Authors: Arturo Cevallos Soto, Jonathan C. Tan, Xiao Hu, Chia-Jung Hsu, Catherine Walsh

Abstract: Inside-Out Planet Formation (IOPF) proposes that the abundant systems of close-in Super-Earths and Mini-Neptunes form in situ at the pressure maximum associated with the Dead Zone Inner Boundary (DZIB). We present a model of physical and chemical evolution of protoplanetary disk midplanes that follows gas advection, radial drift of pebbles and gas-grain chemistry to predict abundances from 300~au… ▽ More Inside-Out Planet Formation (IOPF) proposes that the abundant systems of close-in Super-Earths and Mini-Neptunes form in situ at the pressure maximum associated with the Dead Zone Inner Boundary (DZIB). We present a model of physical and chemical evolution of protoplanetary disk midplanes that follows gas advection, radial drift of pebbles and gas-grain chemistry to predict abundances from 300~au down to the DZIB near 0.2 au. We consider typical disk properties relevant for IOPF, i.e., accretion rates 1E-9 < dM/dt / (Msun/yr) < 1E-8 and viscosity parameter alpha = 1E-4, and evolve for fiducial duration of t = 1E5 years. For outer, cool disk regions, we find that C and up to 90% of O nuclei start locked in CO and O2 ice, which keeps abundances of CO2 and H2O one order of magnitude lower. Radial drift of icy pebbles is influential, with gas-phase abundances of volatiles enhanced up to two orders of magnitude at ice-lines, while the outer disk becomes depleted of dust. Disks with decreasing accretion rates gradually cool, which draws in icelines closer to the star. At <~1 au, advective models yield water-rich gas with C/O ratios <~ 0.1, which may be inherited by atmospheres of planets forming here via IOPF. For planetary interiors built by pebble accretion, IOPF predicts volatile-poor compositions. However, advectively-enhanced volatile mass fractions of ~10% can occur at the water ice line. △ Less

Submitted 13 September, 2022; v1 submitted 4 February, 2022; originally announced February 2022.

Comments: Submitted to MNRAS; comments welcome; for summary of Inside-Out Planet Formation papers see http://www.cosmicorigins.space/iopf

Negative and Positive Feedback from a Supernova Remnant with SHREC: A detailed Study of the Shocked Gas in IC443

Authors: G. Cosentino, I. Jiménez-Serra, J. C. Tan, J. D. Henshaw, A. T. Barnes, C. -Y. Law, S. Zeng, F. Fontani, P. Caselli, S. Viti, S. Zahorecz, F. Rico-Villas, A. Megías, M. Miceli, S. Orlando, S. Ustamujic, E. Greco, G. Peres, F. Bocchino, R. Fedriani, P. Gorai, L. Testi, J. Martín-Pintado

Abstract: Supernova remnants (SNRs) contribute to regulate the star formation efficiency and evolution of galaxies. As they expand into the interstellar medium (ISM), they transfer vast amounts of energy and momentum that displace, compress and heat the surrounding material. Despite the extensive work in galaxy evolution models, it remains to be observationally validated to what extent the molecular ISM is… ▽ More Supernova remnants (SNRs) contribute to regulate the star formation efficiency and evolution of galaxies. As they expand into the interstellar medium (ISM), they transfer vast amounts of energy and momentum that displace, compress and heat the surrounding material. Despite the extensive work in galaxy evolution models, it remains to be observationally validated to what extent the molecular ISM is affected by the interaction with SNRs. We use the first results of the ESO-ARO Public Spectroscopic Survey SHREC, to investigate the shock interaction between the SNR IC443 and the nearby molecular clump G. We use high sensitivity SiO(2-1) and H$^{13}$CO$^+$(1-0) maps obtained by SHREC together with SiO(1-0) observations obtained with the 40m telescope at the Yebes Observatory. We find that the bulk of the SiO emission is arising from the ongoing shock interaction between IC443 and clump G. The shocked gas shows a well ordered kinematic structure, with velocities blue-shifted with respect to the central velocity of the SNR, similar to what observed toward other SNR-cloud interaction sites. The shock compression enhances the molecular gas density, n(H$_2$), up to $>$10$^5$ cm$^{-3}$, a factor of >10 higher than the ambient gas density and similar to values required to ignite star formation. Finally, we estimate that up to 50\% of the momentum injected by IC443 is transferred to the interacting molecular material. Therefore the molecular ISM may represent an important momentum carrier in sites of SNR-cloud interactions. △ Less

Submitted 9 January, 2022; originally announced January 2022.

Comments: Accepted for Publication in MNRAS; 11 pages, 7 figures

Isolated Massive Star Formation in G28.20-0.05

Authors: Chi-Yan Law, Jonathan C. Tan, Prasanta Gorai, Yichen Zhang, Rubén Fedriani, Daniel Tafoya, Kei Tanaka, Giuliana Cosentino, Yao-Lun Yang, Diego Mardones, Maria Teresa Beltrán, Guido Garay

Abstract: We report high-resolution 1.3~mm continuum and molecular line observations of the massive protostar G28.20-0.05 with ALMA. The continuum image reveals a ring-like structure with 2,000~au radius, similar to morphology seen in archival 1.3~cm VLA observations. Based on its spectral index and associated H$30α$ emission, this structure mainly traces ionised gas. However, there is evidence for… ▽ More We report high-resolution 1.3~mm continuum and molecular line observations of the massive protostar G28.20-0.05 with ALMA. The continuum image reveals a ring-like structure with 2,000~au radius, similar to morphology seen in archival 1.3~cm VLA observations. Based on its spectral index and associated H$30α$ emission, this structure mainly traces ionised gas. However, there is evidence for $\sim30$~M$_{\odot}$ of dusty gas near the main mm continuum peak on one side of the ring, as well as in adjacent regions within 3,000~au. A virial analysis on scales of $\sim$2,000~au from hot core line emission yields a dynamical mass of $\sim80\:M_\odot$. A strong velocity gradient in the H$30α$ emission is evidence for a rotating, ionized disk wind, which drives a larger-scale molecular outflow. An infrared SED analysis indicates a current protostellar mass of $m_*\sim40\:M_\odot$ forming from a core with initial mass $M_c\sim300\:M_\odot$ in a clump with mass surface density of $Σ_{\rm cl}\sim 0.8\:{\rm g\:cm}^{-2}$. Thus the SED and other properties of the system can be understood in the context of core accretion models. Structure-finding analysis on the larger-scale continuum image indicates G28.20-0.05 is forming in a relatively isolated environment, with no other concentrated sources, i.e., protostellar cores, above $\sim 1\:M_\odot$ found from $\sim$0.1 to 0.4~pc around the source. This implies that a massive star can form in relative isolation and the dearth of other protostellar companions within the $\sim1$~pc environs is a strong constraint on massive star formation theories that predict the presence of a surrounding protocluster. △ Less

Submitted 8 September, 2022; v1 submitted 4 January, 2022; originally announced January 2022.

Comments: Accepted to ApJ

NIR jets from a clustered region of massive star formation: Morphology and composition in the IRAS 18264-1152 region

Authors: A. R. Costa Silva, R. Fedriani, J. C. Tan, A. Caratti o Garatti, S. Ramsay, V. Rosero, G. Cosentino, P. Gorai, S. Leurini

Abstract: Massive stars form deeply embedded in their parental clouds, making it challenging to directly observe these stars and their immediate environments. It is known that accretion and ejection processes are intrinsically related, thus observing massive protostellar outflows can provide crucial information about the processes governing massive star formation close to the central engine. We aim to probe… ▽ More Massive stars form deeply embedded in their parental clouds, making it challenging to directly observe these stars and their immediate environments. It is known that accretion and ejection processes are intrinsically related, thus observing massive protostellar outflows can provide crucial information about the processes governing massive star formation close to the central engine. We aim to probe the IRAS 18264-1152 (G19.88-0.53) high-mass star-forming complex in the near infrared (NIR) through its molecular hydrogen (H2) jets to analyse the morphology and composition of the line emitting regions and to compare with other outflow tracers. We observed the H2 NIR jets via K-band (1.9-2.5um) observations obtained with the integral field units VLT/SINFONI and VLT/KMOS. SINFONI provides the highest NIR angular resolution achieved so far for the central region (~0.2''). We compared the geometry of the NIR outflows with that of the associated molecular outflow probed by CO (2-1) emission mapped with SMA. We identify nine point sources. Four of these display a rising continuum in the K-band and are BrG emitters, revealing that they are young, potentially jet-driving sources. The spectro-imaging analysis focusses on the H2 jets, for which we derived visual extinction, temperature, column density, area, and mass. The intensity, velocity, and excitation maps based on H2 emission strongly support the existence of a protostellar cluster, with at least two (and up to four) different large-scale outflows. The literature is in agreement with the outflow morphology found here. We derived a stellar density of ~4000 stars pc^-3. Our study reveals the presence of several outflows driven by young sources from a forming cluster of young, massive stars. The derived stellar number density together with the geometry of the outflows suggest that stars can form in a relatively ordered manner in this cluster. △ Less

Submitted 19 July, 2022; v1 submitted 8 December, 2021; originally announced December 2021.

Comments: 23 pages, 14 figures. Published in A&A. Updated to correct typos and match proofs

Journal ref: A&A 659, A23 (2022)

Astrochemical modelling of infrared dark clouds

Authors: Negar Entekhabi, Jonathan C. Tan, Giuliana Cosentino, Chia-Jung Hsu, Paola Caselli, Catherine Walsh, Wanggi Lim, Jonathan D. Henshaw, Ashley T. Barnes, Francesco Fontani, Izaskun Jiménez-Serra

Abstract: Infrared dark clouds (IRDCs) are cold, dense regions of the interstellar medium (ISM) that are likely to represent the initial conditions for massive star formation. It is thus important to study the physical and chemical conditions of IRDCs to provide constraints and inputs for theoretical models of these processes. We aim to determine the astrochemical conditions, especially cosmic ray ionisatio… ▽ More Infrared dark clouds (IRDCs) are cold, dense regions of the interstellar medium (ISM) that are likely to represent the initial conditions for massive star formation. It is thus important to study the physical and chemical conditions of IRDCs to provide constraints and inputs for theoretical models of these processes. We aim to determine the astrochemical conditions, especially cosmic ray ionisation rate (CRIR) and chemical age, in different regions of the massive IRDC G28.37+00.07 by comparing observed abundances of multiple molecules and molecular ions with the predictions of astrochemical models. We have computed a series of single-zone astrochemical models with a gas-grain network that systematically explores the parameter space of density, temperature, CRIR, and visual extinction. We have also investigated the effects of choices of CO ice binding energy and temperatures achieved in transient heating of grains when struck by cosmic rays. We selected 10 positions across the IRDC that are known to have a variety of star formation activity. We utilised mid-infrared (MIR) extinction maps and sub-mm emission maps to measure the mass surface densities of these regions, needed for abundance and volume density estimates. The sub-mm emission maps were also used to measure temperatures. We then used IRAM-30m observations of various tracers to estimate column densities and thus abundances. Using estimates of the abundances of CO, HCO$^+$ and N$_2$H$^+$ we find consistency with astrochemical models that have relatively low CRIRs of $ζ\sim10^{-18}$ to $\sim10^{-17}\:{\rm s}^{-1}$, with no evidence for systematic variation with the level of star formation activity. Astrochemical ages are found to be < 1 Myr. We discuss potential sources of systematic uncertainties in these results and the overall implications for IRDC evolutionary history and astrochemical models.(abridged for arXiv) △ Less

Submitted 23 March, 2022; v1 submitted 7 November, 2021; originally announced November 2021.

Comments: Accepted by A&A

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

MRI-active inner regions of protoplanetary discs. II. Dependence on dust, disc and stellar parameters

Authors: Marija R. Jankovic, Subhanjoy Mohanty, James E. Owen, Jonathan C. Tan

Abstract: Close-in super-Earths are the most abundant exoplanets known. It has been hypothesized that they form in the inner regions of protoplanetary discs, out of the dust that may accumulate at the boundary between the inner region susceptible to the magneto-rotational instability (MRI) and an MRI-dead zone further out. In Paper I we presented a model for the viscous inner disc which includes heating due… ▽ More Close-in super-Earths are the most abundant exoplanets known. It has been hypothesized that they form in the inner regions of protoplanetary discs, out of the dust that may accumulate at the boundary between the inner region susceptible to the magneto-rotational instability (MRI) and an MRI-dead zone further out. In Paper I we presented a model for the viscous inner disc which includes heating due to both irradiation and MRI-driven accretion; thermal and non-thermal ionization; dust opacities; and dust effects on ionization. Here we examine how the inner disc structure varies with stellar, disc and dust parameters. For high accretion rates and small dust grains, we find that: (1) the main sources of ionization are thermal ionization and thermionic and ion emission; (2) the disc features a hot, high-viscosity inner region, and a local gas pressure maximum at the outer edge of this region (in line with previous studies); and (3) an increase in the dust-to-gas ratio pushes the pressure maximum outwards. Consequently, dust can accumulate in such inner discs without suppressing the MRI, with the amount of accumulation depending on the viscosity in the MRI-dead regions. Conversely, for low accretion rates and large dust grains, there appears to be an additional steady-state solution in which: (1) stellar X-rays become the main source of ionization; (2) MRI-viscosity is high throughout the disc; and (3) the pressure maximum ceases to exist. Hence, if planets form in the inner disc, larger accretion rates (and thus younger disks) are favoured. △ Less

Submitted 26 November, 2021; v1 submitted 27 August, 2021; originally announced August 2021.

Comments: accepted to MNRAS

Is there any linkage between interstellar aldehyde and alcohol?

Authors: Suman Kumar Mondal, Prasanta Gorai, Milan Sil, Rana Ghosh, Emmanuel E. Etim, Sandip K Chakrabarti, Takashi Shimonishi, Naoki Nakatani, Kenji Furuya, Jonathan C. Tan, Ankan Das

Abstract: It is speculated that there might be some linkage between interstellar aldehydes and their corresponding alcohols. Here, an observational study and astrochemical modeling are coupled together to illustrate the connection between them. The ALMA Cycle 4 data of a hot molecular core, G10.47+0.03 is utilized for this study. Various aldehydes (acetaldehyde, propanal, and glycolaldehyde), alcohols (meth… ▽ More It is speculated that there might be some linkage between interstellar aldehydes and their corresponding alcohols. Here, an observational study and astrochemical modeling are coupled together to illustrate the connection between them. The ALMA Cycle 4 data of a hot molecular core, G10.47+0.03 is utilized for this study. Various aldehydes (acetaldehyde, propanal, and glycolaldehyde), alcohols (methanol and ethylene glycol), and a ketone (acetone) are identified in this source. The excitation temperatures and the column densities of these species were derived via the rotation diagram method assuming LTE conditions. An extensive investigation is carried out to understand the formation of these species. Six pairs of aldehyde-alcohol: i) methanal and methanol; ii) ethanal and ethanol; iii) propanal and 1-propanol; iv) propenal and allyl alcohol; v) propynal and propargyl alcohol; vi) glycolaldehyde and ethylene glycol; vii) along with one pair of ketone-alcohol (acetone and isopropanol) and viii) ketene-alcohol (ethenone and vinyl alcohol) are considered for this study. Two successive hydrogenation reactions in the ice phase are examined to form these alcohols from aldehydes, ketone, and ketene, respectively. Quantum chemical methods are extensively executed to review the ice phase formation route and the kinetics of these species. Based on the obtained kinetic data, astrochemical modeling is employed to derive the abundances of these aldehydes, alcohols, ketone, and ketene in this source. It is seen that our model could successfully explain the observed abundances of various species in this hot molecular core. △ Less

Submitted 13 August, 2021; originally announced August 2021.

Comments: 34 pages, 11 figures, Accepted for the publication in the Astrophysical Journal

The Core Mass Function Across Galactic Environments. III. Massive Protoclusters

Authors: Theo J. O'Neill, Giuliana Cosentino, Jonathan C. Tan, Yu Cheng, Mengyao Liu

Abstract: The stellar initial mass function (IMF) is fundamental for many areas of astrophysics, but its origin remains poorly understood. It may be inherited from the core mass function (CMF) or arise as a result of more chaotic, competitive accretion. Dense, gravitationally bound cores are seen in molecular clouds and some observations have suggested that the CMF is similar in shape to the IMF, though tra… ▽ More The stellar initial mass function (IMF) is fundamental for many areas of astrophysics, but its origin remains poorly understood. It may be inherited from the core mass function (CMF) or arise as a result of more chaotic, competitive accretion. Dense, gravitationally bound cores are seen in molecular clouds and some observations have suggested that the CMF is similar in shape to the IMF, though translated to higher masses by a factor of $\sim3$. Here we measure the CMF in 28 dense clumps within 3.5 kpc that are likely to be central regions of massive protoclusters, observed via $1.3\:{\rm{mm}}$ dust continuum emission by the ALMAGAL project. We identify 222 cores using the dendrogram algorithm with masses ranging from 0.04 to $252\:M_{\odot}$. We apply completeness corrections for flux and number recovery, estimated from core insertion and recovery experiments. At higher masses, the final derived CMF is well described by a single power law of the form $dN/d\:{\textrm{log}}\:M\propto\:M^{-α}$ with $α\simeq0.94\pm0.08$. However, we find evidence of a break in this power-law behavior between $\sim5$ and $15\:M_{\odot}$, which is, to our knowledge, the first time such a break has been found in distant ($\gtrsim 1$ kpc) regions by ALMA. We compare this massive protocluster CMF with those derived using the same methods in the G286 protocluster and a sample of Infrared Dark Clouds. The massive protocluster CMF is significantly different, i.e., containing more massive cores, which is a potential indication of the role of environment on the CMF and IMF. △ Less

Submitted 26 May, 2021; v1 submitted 18 April, 2021; originally announced April 2021.

Comments: Accepted for publication in ApJ, 25 pages, 10 figures, 4 tables

The Hi-GAL compact source catalogue -- II. The 360° catalogue of clump physical properties

Authors: D. Elia, M. Merello, S. Molinari, E. Schisano, A. Zavagno, D. Russeil, P. Mège, P. G. Martin, L. Olmi, M. Pestalozzi, R. Plume, S. E. Ragan, M. Benedettini, D. J. Eden, T. J. T. Moore, A. Noriega-Crespo, R. Paladini, P. Palmeirim, S. Pezzuto, G. L. Pilbratt, K. L. J. Rygl, P. Schilke, F. Strafella, J. C. Tan, A. Traficante , et al. (7 additional authors not shown)

Abstract: We present the $360^\circ$ catalogue of physical properties of Hi-GAL compact sources, detected between 70 and 500 $μ$m. This release not only completes the analogous catalogue previously produced by the Hi-GAL collaboration for $-71^\circ \lesssim \ell \lesssim 67^\circ$, but also meaningfully improves it thanks to a new set of heliocentric distances, 120808 in total. About a third of the 150223… ▽ More We present the $360^\circ$ catalogue of physical properties of Hi-GAL compact sources, detected between 70 and 500 $μ$m. This release not only completes the analogous catalogue previously produced by the Hi-GAL collaboration for $-71^\circ \lesssim \ell \lesssim 67^\circ$, but also meaningfully improves it thanks to a new set of heliocentric distances, 120808 in total. About a third of the 150223 entries are located in the newly added portion of the Galactic plane. A first classification based on detection at 70 $μ$m as a signature of ongoing star-forming activity distinguishes between protostellar sources (23~per cent of the total) and starless sources, with the latter further classified as gravitationally bound (pre-stellar) or unbound. The integral of the spectral energy distribution, including ancillary photometry from $λ=21$ to 1100 $μ$m, gives the source luminosity and other bolometric quantities, while a modified black body fitted to data for $λ\geq 160\, μ$m yields mass and temperature. All tabulated clump properties are then derived using photometry and heliocentric distance, where possible. Statistics of these quantities are discussed with respect to both source Galactic location and evolutionary stage. No strong differences in the distributions of evolutionary indicators are found between the inner and outer Galaxy. However, masses and densities in the inner Galaxy are on average significantly larger, resulting in a higher number of clumps that are candidates to host massive star formation. Median behaviour of distance-independent parameters tracing source evolutionary status is examined as a function of the Galactocentric radius, showing no clear evidence of correlation with spiral arm positions. △ Less

Submitted 10 April, 2021; originally announced April 2021.

Comments: accepted by MNRAS, april 2021

Inside-Out Planet Formation: VI. Oligarchic Coagulation of Planetesimals from a Pebble Ring?

Authors: Maxwell X. Cai, Jonathan C. Tan, Simon Portegies Zwart

Abstract: Inside-Out Planet Formation (IOPF) is a theory addressing the origin of Systems of Tightly-Packed Inner Planets (STIPs) via {\it in situ} formation and growth of the planets. It predicts that a pebble ring is established at the pressure maximum associated with the dead zone inner boundary (DZIB) with an inner disk magnetorotational instability (MRI)-active region. Using direct $N$-body simulations… ▽ More Inside-Out Planet Formation (IOPF) is a theory addressing the origin of Systems of Tightly-Packed Inner Planets (STIPs) via {\it in situ} formation and growth of the planets. It predicts that a pebble ring is established at the pressure maximum associated with the dead zone inner boundary (DZIB) with an inner disk magnetorotational instability (MRI)-active region. Using direct $N$-body simulations, we study the collisional evolution of planetesimals formed from such a pebble ring, in particular examining whether a single dominant planet emerges. We consider a variety of models, including some in which the planetesimals are continuing to grow via pebble accretion. We find that the planetesimal ring undergoes oligarchic evolution, and typically turns into 2 or 3 surviving oligarchs on nearly coplanar and circular orbits, independent of the explored initial conditions or form of pebble accretion. The most massive oligarchs typically consist of about $70\%$ of the total mass, with the building-up process typically finishing within $\sim 10^5$ years. However, a relatively massive secondary planet always remains with $\sim30-65\%$ of the mass of the primary. Such secondary planets have properties that are inconsistent with the observed properties of the innermost pairs of planets in STIPs. Thus, for IOPF to be a viable theory for STIP formation, it needs to be shown how oligarchic growth of a relatively massive secondary from the initial pebble ring can be avoided. We discuss some potential additional physical processes that should be included in the modeling and explored as next steps. △ Less

Submitted 13 December, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

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

ALMA-IRDC II. First high-angular resolution measurements of the 14N/15N ratio in a large sample of infrared-dark cloud cores

Authors: F. Fontani, A. T. Barnes, P. Caselli, J. D. Henshaw, G. Cosentino, I. Jiménez-Serra, J. C. Tan, J. E. Pineda, C. Y. Law

Abstract: The 14N/15N ratio in molecules exhibits a large variation in star-forming regions, especially when measured from N2H+ isotopologues. However, there are only a few studies performed at high-angular resolution. We present the first interferometric survey of the 14N/15N ratio in N2H+ obtained with the Atacama Large Millimeter Array towards four infrared-dark clouds harbouring 3~mm continuum cores ass… ▽ More The 14N/15N ratio in molecules exhibits a large variation in star-forming regions, especially when measured from N2H+ isotopologues. However, there are only a few studies performed at high-angular resolution. We present the first interferometric survey of the 14N/15N ratio in N2H+ obtained with the Atacama Large Millimeter Array towards four infrared-dark clouds harbouring 3~mm continuum cores associated with different physical properties. We detect N15NH+ (1-0) in about 20-40% of the cores, depending on the host cloud. The 14N/15N values measured towards the millimeter continuum cores range from a minimum of 80 up to a maximum of 400. The spread of values is narrower than that found in any previous single-dish survey of high-mass star-forming regions, and than that obtained using the total power data only. This suggests that the 14N/15N ratio is on average higher in the diffuse gaseous envelope of the cores, and stresses the need for high-angular resolution maps to measure correctly the 14N/15N ratio in dense cores embedded in IRDCs. The average 14N/15N ratio of 210 is also lower than the interstellar value at the Galactocentric distance of the clouds (300-330), although the sensitivity of our observations does not allow us to unveil 14N/15N ratios higher than 400. No clear trend is found between the 14N/15N ratio and the core physical properties. We find only a tentative positive trend between 14N/15N and the H2 column density. However, firmer conclusions can be drawn only with higher sensitivity measurements. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 19 pages, 7 figures, 6 tables, 1 appendix, accepted for publication in MNRAS

ALMA-IRDC: Dense gas mass distribution from cloud to core scales

Authors: A. T. Barnes, J. D. Henshaw, F. Fontani, J. E. Pineda, G. Cosentino, J. C. Tan, P. Caselli, I. Jiménez-Serra, C. Y. Law, A. Avison, F. Bigiel, S. Feng, S. Kong, S. N. Longmore, L. Moser, R. J. Parker, Á. Sánchez-Monge, K. Wang

Abstract: Infrared dark clouds (IRDCs) are potential hosts of the elusive early phases of high-mass star formation (HMSF). Here we conduct an in-depth analysis of the fragmentation properties of a sample of 10 IRDCs, which have been highlighted as some of the best candidates to study HMSF within the Milky Way. To do so, we have obtained a set of large mosaics covering these IRDCs with ALMA at band 3 (or 3mm… ▽ More Infrared dark clouds (IRDCs) are potential hosts of the elusive early phases of high-mass star formation (HMSF). Here we conduct an in-depth analysis of the fragmentation properties of a sample of 10 IRDCs, which have been highlighted as some of the best candidates to study HMSF within the Milky Way. To do so, we have obtained a set of large mosaics covering these IRDCs with ALMA at band 3 (or 3mm). These observations have a high angular resolution (~3arcsec or ~0.05pc), and high continuum and spectral line sensitivity (~0.15mJy/beam and ~0.2K per 0.1km/s channel at the N2H+(1-0) transition). From the dust continuum emission, we identify 96 cores ranging from low- to high-mass (M = 3.4 to 50.9Msun) that are gravitationally bound (alpha_vir = 0.3 to 1.3) and which would require magnetic field strengths of B = 0.3 to 1.0mG to be in virial equilibrium. We combine these results with a homogenised catalogue of literature cores to recover the hierarchical structure within these clouds over four orders of magnitude in spatial scale (0.01pc to 10pc). Using supplementary observations at an even higher angular resolution, we find that the smallest fragments (<0.02pc) within this hierarchy do not currently have the mass and/or the density required to form high-mass stars. Nonetheless, the new ALMA observations presented in this paper have facilitated the identification of 19 (6 quiescent and 13 star-forming) cores that retain >16Msun without further fragmentation. These high-mass cores contain trans-sonic non-thermal motions, are kinematically sub-virial, and require moderate magnetic field strengths for support against collapse. The identification of these potential sites of high-mass star formation represents a key step in allowing us to test the predictions from high-mass star and cluster formation theories. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 21(+6 in appendix) pages, 12(+3) figures, 3(+3) tables. Machine-readable versions of Table A1, A2 and A3 are available online. Accepted for publication in MNRAS

MRI-active inner regions of protoplanetary discs. I. A detailed model of disc structure

Authors: Marija R. Jankovic, James E. Owen, Subhanjoy Mohanty, Jonathan C. Tan

Abstract: Short-period super-Earth-sized planets are common. Explaining how they form near their present orbits requires understanding the structure of the inner regions of protoplanetary discs. Previous studies have argued that the hot inner protoplanetary disc is unstable to the magneto-rotational instability (MRI) due to thermal ionization of potassium, and that a local gas pressure maximum forms at the… ▽ More Short-period super-Earth-sized planets are common. Explaining how they form near their present orbits requires understanding the structure of the inner regions of protoplanetary discs. Previous studies have argued that the hot inner protoplanetary disc is unstable to the magneto-rotational instability (MRI) due to thermal ionization of potassium, and that a local gas pressure maximum forms at the outer edge of this MRI-active zone. Here we present a steady-state model for inner discs accreting viscously, primarily due to the MRI. The structure and MRI-viscosity of the inner disc are fully coupled in our model; moreover, we account for many processes omitted in previous such models, including disc heating by both accretion and stellar irradiation, vertical energy transport, realistic dust opacities, dust effects on disc ionization and non-thermal sources of ionization. For a disc around a solar-mass star with a standard gas accretion rate ($\dot{M}$$\sim$$10^{-8}$M$_\odot$yr$^{-1}$) and small dust grains, we find that the inner disc is optically thick, and the accretion heat is primarily released near the midplane. As a result, both the disc midplane temperature and the location of the pressure maximum are only marginally affected by stellar irradiation, and the inner disc is also convectively unstable. As previously suggested, the inner disc is primarily ionized through thermionic and potassium ion emission from dust grains, which, at high temperatures, counteract adsorption of free charges onto grains. Our results show that the location of the pressure maximum is determined by the threshold temperature above which thermionic and ion emission become efficient. △ Less

Submitted 1 April, 2021; v1 submitted 25 February, 2021; originally announced February 2021.

Comments: accepted for publication in MNRAS

Star Formation in a Strongly Magnetized Cloud

Authors: Yu Cheng, Jonathan C. Tan, Paola Caselli, Laura Fissel, Hector G. Arce, Francesco Fontani, Matthew D. Goodson, Mengyao Liu, Nicholas Galitzki

Abstract: We study star formation in the Center Ridge 1 (CR1) clump in the Vela C giant molecular cloud, selected as a high column density region that shows the lowest level of dust continuum polarization angle dispersion, likely indicating that the magnetic field is relatively strong. We observe the source with the ALMA 7m-array at 1.05~mm and 1.3~mm wavelengths, which enable measurements of dust temperatu… ▽ More We study star formation in the Center Ridge 1 (CR1) clump in the Vela C giant molecular cloud, selected as a high column density region that shows the lowest level of dust continuum polarization angle dispersion, likely indicating that the magnetic field is relatively strong. We observe the source with the ALMA 7m-array at 1.05~mm and 1.3~mm wavelengths, which enable measurements of dust temperature, core mass and astrochemical deuteration. A relatively modest number of eleven dense cores are identified via their dust continuum emission, with masses spanning from 0.17 to 6.7 Msun. Overall CR1 has a relatively low compact dense gas fraction compared with other typical clouds with similar column densities, which may be a result of the strong magnetic field and/or the very early evolutionary stage of this region. The deuteration ratios, Dfrac, of the cores, measured with N2H+(3-2) and N2D+(3-2) lines, span from 0.011 to 0.85, with the latter being one of the highest values yet detected. The level of deuteration appears to decrease with evolution from prestellar to protostellar phase. A linear filament, running approximately parallel with the large scale magnetic field orientation, is seen connecting the two most massive cores, each having CO bipolar outflows aligned orthogonally to the filament. The filament contains the most deuterated core, likely to be prestellar and located midway between the protostars. The observations permit measurement of the full deuteration structure of the filament along its length, which we present. We also discuss the kinematics and dynamics of this structure, as well as of the dense core population. △ Less

Submitted 25 May, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: 25 pages, 17 figures, accepted by ApJ

Carbon Chain Chemistry in Hot-Core Regions around Three Massive Young Stellar Objects Associated with 6.7 GHz Methanol Masers

Authors: Kotomi Taniguchi, Eric Herbst, Liton Majumdar, Paola Caselli, Jonathan C. Tan, Zhi-Yun Li, Tomomi Shimoikura, Kazuhito Dobashi, Fumitaka Nakamura, Masao Saito

Abstract: We have carried out observations of CCH ($N=1-0$), CH$_{3}$CN ($J=5-4$), and three $^{13}$C isotopologues of HC$_{3}$N ($J=10-9$) toward three massive young stellar objects (MYSOs), G12.89+0.49, G16.86--2.16, and G28.28--0.36, with the Nobeyama 45-m radio telescope. Combined with previous results on HC$_{5}$N, the column density ratios of $N$(CCH)/$N$(HC$_{5}$N), hereafter the CCH/HC$_{5}$N ratios… ▽ More We have carried out observations of CCH ($N=1-0$), CH$_{3}$CN ($J=5-4$), and three $^{13}$C isotopologues of HC$_{3}$N ($J=10-9$) toward three massive young stellar objects (MYSOs), G12.89+0.49, G16.86--2.16, and G28.28--0.36, with the Nobeyama 45-m radio telescope. Combined with previous results on HC$_{5}$N, the column density ratios of $N$(CCH)/$N$(HC$_{5}$N), hereafter the CCH/HC$_{5}$N ratios, in the MYSOs are derived to be $\sim 15$. This value is lower than that in a low-mass warm carbon chain chemistry (WCCC) source by more than one order of magnitude. We compare the observed CCH/HC$_{5}$N ratios with hot-core model calculations (Taniguchi et al. 2019). The observed ratios in the MYSOs can be best reproduced by models when the gas temperature is $\sim 85$ K, which is higher than in L1527, a low-mass WCCC source ($\sim 35$ K). These results suggest that carbon-chain molecules detected around the MYSOs exist at least partially in higher temperature regions than those in low-mass WCCC sources. There is no significant difference in column density among the three $^{13}$C isotopologues of HC$_{3}$N in G12.89+0.49 and G16.86-2.16, while HCC$^{13}$CN is more abundant than the others in G28.28--0.36. We discuss carbon-chain chemistry around the three MYSOs based on the CCH/HC$_{5}$N ratio and the $^{13}$C isotopic fractionation of HC$_{3}$N. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Comments: Accepted by The Astrophysical Journal, 17 pages, 6 figures, 6 tables

Photodissociation Region Diagnostics Across Galactic Environments

Authors: Thomas G. Bisbas, Jonathan C. Tan, Kei E. I. Tanaka

Abstract: We present three-dimensional astrochemical simulations and synthetic observations of magnetised, turbulent, self-gravitating molecular clouds. We explore various galactic interstellar medium environments, including cosmic-ray ionization rates in the range of $ζ_{\rm CR}=10^{-17}$-$10^{-14}\,{\rm s}^{-1}$, far-UV intensities in the range of $G_0=1$-$10^3$ and metallicities in the range of $Z=0.1$-… ▽ More We present three-dimensional astrochemical simulations and synthetic observations of magnetised, turbulent, self-gravitating molecular clouds. We explore various galactic interstellar medium environments, including cosmic-ray ionization rates in the range of $ζ_{\rm CR}=10^{-17}$-$10^{-14}\,{\rm s}^{-1}$, far-UV intensities in the range of $G_0=1$-$10^3$ and metallicities in the range of $Z=0.1$-$2\,{\rm Z}_{\odot}$. The simulations also probe a range of densities and levels of turbulence, including cases where the gas has undergone recent compression due to cloud-cloud collisions. We examine: i) the column densities of carbon species across the cycle of CII, CI and CO, along with OI, in relation to the HI-to-H$_2$ transition; ii) the velocity-integrated emission of [CII]~$158μ$m, [$^{13}$CII]~$158μ$m, [CI]~$609μ$m and $370μ$m, [OI]~$63μ$m and $146μ$m, and of the first ten $^{12}$CO rotational transitions; iii) the corresponding Spectral Line Energy Distributions; iv) the usage of [CII] and [OI]~$63μ$m to describe the dynamical state of the clouds; v) the behavior of the most commonly used ratios between transitions of CO and [CI]; and vi) the conversion factors for using CO and CI as H$_2$-gas tracers. We find that enhanced cosmic-ray energy densities enhance all aforementioned line intensities. At low metallicities, the emission of [CII] is well connected with the H$_2$ column, making it a promising new H$_2$ tracer in metal-poor environments. The conversion factors of $X_{\rm CO}$ and $X_{\rm CI}$ depend on metallicity and the cosmic-ray ionization rate, but not on FUV intensity. In the era of ALMA, SOFIA and the forthcoming CCAT-prime telescope, our results can be used to understand better the behaviour of systems in a wide range of galactic and extragalactic environments. △ Less

Submitted 12 December, 2020; originally announced December 2020.

Comments: 33 pages, 24 figures, 3 tables. Submitted to MNRAS. Comments welcome!

SiO Outflows as Tracers of Massive Star Formation in Infrared Dark Clouds

Authors: Mengyao Liu, Jonathan C. Tan, Joshua Marvil, Shuo Kong, Viviana Rosero, Paola Caselli, Giuliana Cosentino

Abstract: To study the early phases of massive star formation, we present ALMA observations of SiO(5-4) emission and VLA observations of 6 cm continuum emission towards 32 Infrared Dark Cloud (IRDC) clumps, spatially resolved down to $\lesssim 0.05$ pc. Out of the 32 clumps, we detect SiO emission in 20 clumps, and in 11 of them the SiO emission is relatively strong and likely tracing protostellar outflows.… ▽ More To study the early phases of massive star formation, we present ALMA observations of SiO(5-4) emission and VLA observations of 6 cm continuum emission towards 32 Infrared Dark Cloud (IRDC) clumps, spatially resolved down to $\lesssim 0.05$ pc. Out of the 32 clumps, we detect SiO emission in 20 clumps, and in 11 of them the SiO emission is relatively strong and likely tracing protostellar outflows. Some SiO outflows are collimated, while others are less ordered. For the six strongest SiO outflows, we estimate basic outflow properties. In our entire sample, where there is SiO emission, we find 1.3 mm continuum and infrared emission nearby, but not vice versa. We build the spectral energy distributions (SEDs) of cores with 1.3 mm continuum emission and fit them with radiative transfer (RT) models. The low luminosities and stellar masses returned by SED fitting suggest these are early stage protostars. We see a slight trend of increasing SiO line luminosity with bolometric luminosity, which suggests more powerful shocks in the vicinity of more massive YSOs. We do not see a clear relation between the SiO luminosity and the evolutionary stage indicated by $L/M$. We conclude that as a protostar approaches a bolometric luminosity of $\sim 10^2 \: L_{\odot}$, the shocks in the outflow are generally strong enough to form SiO emission. The VLA 6 cm observations toward the 15 clumps with the strongest SiO emission detect emission in four clumps, which is likely shock ionized jets associated with the more massive ones of these protostellar cores. △ Less

Submitted 6 June, 2021; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 42 pages, 18 figures, accepted to ApJ