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

Investigating Mutual Coupling in the Hydrogen Epoch of Reionization Array and Mitigating its Effects on the 21-cm Power Spectrum

Authors: E. Rath, R. Pascua, A. T. Josaitis, A. Ewall-Wice, N. Fagnoni, E. de Lera Acedo, Z. E. Martinot, Z. Abdurashidova, T. Adams, J. E. Aguirre, R. Baartman, A. P. Beardsley, L. M. Berkhout, G. Bernardi, T. S. Billings, J. D. Bowman, P. Bull, J. Burba, R. Byrne, S. Carey, K. -F. Chen, S. Choudhuri, T. Cox, D. R. DeBoer, M. Dexter , et al. (56 additional authors not shown)

Abstract: Interferometric experiments designed to detect the highly redshifted 21-cm signal from neutral hydrogen are producing increasingly stringent constraints on the 21-cm power spectrum, but some k-modes remain systematics-dominated. Mutual coupling is a major systematic that must be overcome in order to detect the 21-cm signal, and simulations that reproduce effects seen in the data can guide strategi… ▽ More Interferometric experiments designed to detect the highly redshifted 21-cm signal from neutral hydrogen are producing increasingly stringent constraints on the 21-cm power spectrum, but some k-modes remain systematics-dominated. Mutual coupling is a major systematic that must be overcome in order to detect the 21-cm signal, and simulations that reproduce effects seen in the data can guide strategies for mitigating mutual coupling. In this paper, we analyse 12 nights of data from the Hydrogen Epoch of Reionization Array and compare the data against simulations that include a computationally efficient and physically motivated semi-analytic treatment of mutual coupling. We find that simulated coupling features qualitatively agree with coupling features in the data; however, coupling features in the data are brighter than the simulated features, indicating the presence of additional coupling mechanisms not captured by our model. We explore the use of fringe-rate filters as mutual coupling mitigation tools and use our simulations to investigate the effects of mutual coupling on a simulated cosmological 21-cm power spectrum in a "worst case" scenario where the foregrounds are particularly bright. We find that mutual coupling contaminates a large portion of the "EoR Window", and the contamination is several orders-of-magnitude larger than our simulated cosmic signal across a wide range of cosmological Fourier modes. While our fiducial fringe-rate filtering strategy reduces mutual coupling by roughly a factor of 100 in power, a non-negligible amount of coupling cannot be excised with fringe-rate filters, so more sophisticated mitigation strategies are required. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 19 pages, 12 figures, submitted to MNRAS

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 a brief overview of the program and preliminary results regarding 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 galaxy. Two further objects are clear AGN candidates at z=2.0 and 3.2, based on morphology and/or spectroscopy, in particular 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 3x10^{-4} cMpc^{-3}. After accounting for variability and luminosity incompleteness, we estimate n_SMBH \gtrsim 8x10^{-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 24 March, 2024; originally announced March 2024.

Comments: Submitted to ApJ

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

A demonstration of the effect of fringe-rate filtering in the Hydrogen Epoch of Reionization Array delay power spectrum pipeline

Authors: Hugh Garsden, Philip Bull, Mike Wilensky, Zuhra Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Lindsay M. Berkhout, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Steven Carey, Chris L. Carilli, Kai-Feng Chen, Carina Cheng, Samir Choudhuri, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter , et al. (72 additional authors not shown)

Abstract: Radio interferometers targeting the 21cm brightness temperature fluctuations at high redshift are subject to systematic effects that operate over a range of different timescales. These can be isolated by designing appropriate Fourier filters that operate in fringe-rate (FR) space, the Fourier pair of local sidereal time (LST). Applications of FR filtering include separating effects that are correl… ▽ More Radio interferometers targeting the 21cm brightness temperature fluctuations at high redshift are subject to systematic effects that operate over a range of different timescales. These can be isolated by designing appropriate Fourier filters that operate in fringe-rate (FR) space, the Fourier pair of local sidereal time (LST). Applications of FR filtering include separating effects that are correlated with the rotating sky vs. those relative to the ground, down-weighting emission in the primary beam sidelobes, and suppressing noise. FR filtering causes the noise contributions to the visibility data to become correlated in time however, making interpretation of subsequent averaging and error estimation steps more subtle. In this paper, we describe fringe rate filters that are implemented using discrete prolate spheroidal sequences, and designed for two different purposes -- beam sidelobe/horizon suppression (the `mainlobe' filter), and ground-locked systematics removal (the `notch' filter). We apply these to simulated data, and study how their properties affect visibilities and power spectra generated from the simulations. Included is an introduction to fringe-rate filtering and a demonstration of fringe-rate filters applied to simple situations to aid understanding. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 21 pages, 18 figures, submitted to Monthly Notices of the Royal Astronomical Society

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

Hydrogen Epoch of Reionization Array (HERA) Phase II Deployment and Commissioning

Authors: Lindsay M. Berkhout, Daniel C. Jacobs, Zuhra Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Philip Bull, Jacob Burba, Steven Carey, Chris L. Carilli, Kai-Feng Chen, Carina Cheng, Samir Choudhuri, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Joshua S. Dillon , et al. (71 additional authors not shown)

Abstract: This paper presents the design and deployment of the Hydrogen Epoch of Reionization Array (HERA) phase II system. HERA is designed as a staged experiment targeting 21 cm emission measurements of the Epoch of Reionization. First results from the phase I array are published as of early 2022, and deployment of the phase II system is nearing completion. We describe the design of the phase II system an… ▽ More This paper presents the design and deployment of the Hydrogen Epoch of Reionization Array (HERA) phase II system. HERA is designed as a staged experiment targeting 21 cm emission measurements of the Epoch of Reionization. First results from the phase I array are published as of early 2022, and deployment of the phase II system is nearing completion. We describe the design of the phase II system and discuss progress on commissioning and future upgrades. As HERA is a designated Square Kilometer Array (SKA) pathfinder instrument, we also show a number of "case studies" that investigate systematics seen while commissioning the phase II system, which may be of use in the design and operation of future arrays. Common pathologies are likely to manifest in similar ways across instruments, and many of these sources of contamination can be mitigated once the source is identified. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Journal ref: PASP 2024 136 045002

matvis: A matrix-based visibility simulator for fast forward modelling of many-element 21 cm arrays

Authors: Piyanat Kittiwisit, Steven G. Murray, Hugh Garsden, Philip Bull, Christopher Cain, Aaron R. Parsons, Jackson Sipple, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Adam P. Beardsley, Lindsay M. Berkhout, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Jacob Burba, Steven Carey, Chris L. Carilli, Kai-Feng Chen, Carina Cheng , et al. (73 additional authors not shown)

Abstract: Detection of the faint 21 cm line emission from the Cosmic Dawn and Epoch of Reionisation will require not only exquisite control over instrumental calibration and systematics to achieve the necessary dynamic range of observations but also validation of analysis techniques to demonstrate their statistical properties and signal loss characteristics. A key ingredient in achieving this is the ability… ▽ More Detection of the faint 21 cm line emission from the Cosmic Dawn and Epoch of Reionisation will require not only exquisite control over instrumental calibration and systematics to achieve the necessary dynamic range of observations but also validation of analysis techniques to demonstrate their statistical properties and signal loss characteristics. A key ingredient in achieving this is the ability to perform high-fidelity simulations of the kinds of data that are produced by the large, many-element, radio interferometric arrays that have been purpose-built for these studies. The large scale of these arrays presents a computational challenge, as one must simulate a detailed sky and instrumental model across many hundreds of frequency channels, thousands of time samples, and tens of thousands of baselines for arrays with hundreds of antennas. In this paper, we present a fast matrix-based method for simulating radio interferometric measurements (visibilities) at the necessary scale. We achieve this through judicious use of primary beam interpolation, fast approximations for coordinate transforms, and a vectorised outer product to expand per-antenna quantities to per-baseline visibilities, coupled with standard parallelisation techniques. We validate the results of this method, implemented in the publicly-available matvis code, against a high-precision reference simulator, and explore its computational scaling on a variety of problems. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 25 pages, 20 figures, submitted to RAS Techniques and Instruments, matvis is publicly available at https://github.com/HERA-Team/matvis

Constraints on Brans-Dicke gravity from neutron star-black hole merger events using higher harmonics

Authors: Jing Tan, Baoxiang Wang

Abstract: In this paper, we derive a 90\% credible lower bound on the modified parameter of scalar-tensor theories as $\varphi_{-2}>-7.94\times10^{-4}$ by using dominant-mode correction. Specific to BD theory, we have the constraint $ω_{\rm BD}>4.75$. Asymmetric binary systems usually have a significant mass ratio; in such cases, higher harmonic modes cannot be neglected. Our work considers higher harmonic… ▽ More In this paper, we derive a 90\% credible lower bound on the modified parameter of scalar-tensor theories as $\varphi_{-2}>-7.94\times10^{-4}$ by using dominant-mode correction. Specific to BD theory, we have the constraint $ω_{\rm BD}>4.75$. Asymmetric binary systems usually have a significant mass ratio; in such cases, higher harmonic modes cannot be neglected. Our work considers higher harmonic corrections from scalar-tensor theories and provides a tighter constraint of $\varphi_{-2}>-7.59\times10^{-4}$. Transitioning to the BD theory, the constraint is $ω_{\rm BD}>5.06$, with a 6.5\% improvement. We also consider a plausible NSBH event, GW190814, which is a highly unequal mass ratio source and exhibits strong evidence for higher-order multipoles. We obtain poorly converged results when using the dominant mode while getting a constraint of $\varphi_{-2}>-6.60\times10^{-4}$ on scalar-tensor theories when including the higher harmonic modes. This suggests that the difference between the dominant mode and higher modes has a significant impact on our analysis. Furthermore, treating this suspected event as an NSBH event, we find $ω_{\rm BD}>6.12$ when including the higher harmonic modes. Combining GW200115 and GW190814 and including higher modes, the constraint is improved to $ω_{\rm BD}>110.55$. This is currently the strongest constraint utilizing GWs, contingent upon GW190814 being an NSBH event. Additionally, we take into account a BD-like theory, known as screened modified gravity (SMG), and investigate the coupling constant constraints, both with and without higher-mode corrections, by using data from both GW200115 and GW190814. △ Less

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

Comments: 8 pages, 2 figures, 2 tables

Journal ref: Phys.Rev.D 109 (2024) 8, 084036

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

Ejecta Evolution Following a Planned Impact into an Asteroid: The First Five Weeks

Authors: Theodore Kareta, Cristina Thomas, Jian-Yang Li, Matthew M. Knight, Nicholas Moskovitz, Agata Rozek, Michele T. Bannister, Simone Ieva, Colin Snodgrass, Petr Pravec, Eileen V. Ryan, William H. Ryan, Eugene G. Fahnestock, Andrew S. Rivkin, Nancy Chabot, Alan Fitzsimmons, David Osip, Tim Lister, Gal Sarid, Masatoshi Hirabayashi, Tony Farnham, Gonzalo Tancredi, Patrick Michel, Richard Wainscoat, Rob Weryk , et al. (63 additional authors not shown)

Abstract: The impact of the DART spacecraft into Dimorphos, moon of the asteroid Didymos, changed Dimorphos' orbit substantially, largely from the ejection of material. We present results from twelve Earth-based facilities involved in a world-wide campaign to monitor the brightness and morphology of the ejecta in the first 35 days after impact. After an initial brightening of ~1.4 magnitudes, we find consis… ▽ More The impact of the DART spacecraft into Dimorphos, moon of the asteroid Didymos, changed Dimorphos' orbit substantially, largely from the ejection of material. We present results from twelve Earth-based facilities involved in a world-wide campaign to monitor the brightness and morphology of the ejecta in the first 35 days after impact. After an initial brightening of ~1.4 magnitudes, we find consistent dimming rates of 0.11-0.12 magnitudes/day in the first week, and 0.08-0.09 magnitudes/day over the entire study period. The system returned to its pre-impact brightness 24.3-25.3 days after impact through the primary ejecta tail remained. The dimming paused briefly eight days after impact, near in time to the appearance of the second tail. This was likely due to a secondary release of material after re-impact of a boulder released in the initial impact, through movement of the primary ejecta through the aperture likely played a role. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 16 pages, 5 Figures, accepted in the Astrophysical Journal Letters (ApJL) on October 16, 2023

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)

The First GECAM Observation Results on Terrestrial Gamma-ray Flashes and Terrestrial Electron Beams

Authors: Y. Zhao, J. C. Liu, S. L. Xiong, W. C. Xue, Q. B. Yi, G. P. Lu, W. Xu, F. C. Lyu, J. C. Sun, W. X. Peng, C. Zheng, Y. Q. Zhang, C. Cai, S. Xiao, S. L. Xie, C. W. Wang, W. J. Tan, Z. H. An, G. Chen, Y. Q. Du, Y. Huang, M. Gao, K. Gong, D. Y. Guo, J. J. He , et al. (37 additional authors not shown)

Abstract: Gravitational-wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) is a space-borne instrument dedicated to monitoring high-energy transients, including Terrestrial Gamma-ray Flashes (TGFs) and Terrestrial Electron Beams (TEBs). We implemented a TGF/TEB search algorithm for GECAM, with which 147 bright TGFs, 2 typical TEBs and 2 special TEB-like events are identified during an effe… ▽ More Gravitational-wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) is a space-borne instrument dedicated to monitoring high-energy transients, including Terrestrial Gamma-ray Flashes (TGFs) and Terrestrial Electron Beams (TEBs). We implemented a TGF/TEB search algorithm for GECAM, with which 147 bright TGFs, 2 typical TEBs and 2 special TEB-like events are identified during an effective observation time of $\sim$9 months. We show that, with gamma-ray and charged particle detectors, GECAM can effectively identify and distinguish TGFs and TEBs, and measure their temporal and spectral properties in detail. A very high TGF-lightning association rate of $\sim$80\% is obtained between GECAM and GLD360 in east Asia region. △ Less

Submitted 17 June, 2023; originally announced June 2023.

Comments: The paper was accepted by Geophysical Research Letters on June 16th, 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)

Double-dipping to refine stellar rotation periods

Authors: Joanne Tan, Gibor Basri

Abstract: We present a refined analysis of 15038 Kepler main sequence light curves to determine the stellar rotation periods. The initial period estimates come from an autocorrelation function, as has been done before. We then measure the duration of every intensity dip in the light curve, expressed as fractions of the initial rotation period estimate. These dip duration distributions are subdivided into se… ▽ More We present a refined analysis of 15038 Kepler main sequence light curves to determine the stellar rotation periods. The initial period estimates come from an autocorrelation function, as has been done before. We then measure the duration of every intensity dip in the light curve, expressed as fractions of the initial rotation period estimate. These dip duration distributions are subdivided into several regions whose relation to each other helps determine which harmonic of the initial rotation period is most physically plausible. We compare our final rotation periods to those from McQuillan, Mazeh, & Aigrain (2014) and find that the great majority agree, but about 10% of their periods are doubtful (usually twice as long as is most plausible). We are still refining our method, and will later extend it to more stars to substantially increase the sample of reliable stellar rotation periods. △ Less

Submitted 4 June, 2023; originally announced June 2023.

Comments: Published in AN on 3 August 2020. 6 pages, 6 figures, preprint version

Journal ref: Volume 341, Issue 5, 2020, Pages 513-518

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

As a matter of colon: I am NOT digging cheeky titles (no, but actually yes :>)

Authors: Joanne Tan, Tie Sien Suk

Abstract: What's in a name, a poet once asked. To which we present this work, where we investigate the importance of a paper title in ensuring its best outcome. We queried astronomy papers using NASA ADS and ranked 6000 of them in terms of cheekiness level. We investigate the correlation between citation counts and (i) the presence of a colon, and (ii) cheekiness ranking. We conclude that colon matters in t… ▽ More What's in a name, a poet once asked. To which we present this work, where we investigate the importance of a paper title in ensuring its best outcome. We queried astronomy papers using NASA ADS and ranked 6000 of them in terms of cheekiness level. We investigate the correlation between citation counts and (i) the presence of a colon, and (ii) cheekiness ranking. We conclude that colon matters in the anatomy of a paper title. So does trying to be cheeky, but we find that too much cheekiness can lead to cringefests. Striking the right balance is therefore crucial. May we recommend aiming for a level 4 cheekiness on a scale of 1-5. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Comments: 4 pages, 4 figures, 1 meme. Submitted to Acta Prima Aprila

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

Improved Constraints on the 21 cm EoR Power Spectrum and the X-Ray Heating of the IGM with HERA Phase I Observations

Authors: The HERA Collaboration, Zara Abdurashidova, Tyrone Adams, James E. Aguirre, Paul Alexander, Zaki S. Ali, Rushelle Baartman, Yanga Balfour, Rennan Barkana, Adam P. Beardsley, Gianni Bernardi, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley, Daniela Breitman, Philip Bull, Jacob Burba, Steve Carey, Chris L. Carilli, Carina Cheng, Samir Choudhuri, David R. DeBoer, Eloy de Lera Acedo, Matt Dexter, Joshua S. Dillon , et al. (70 additional authors not shown)

Abstract: We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits (HERA Collaboration 2022a), we find at 95% confidence that $Δ^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that… ▽ More We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits (HERA Collaboration 2022a), we find at 95% confidence that $Δ^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that $Δ^2 (k = 0.36$ $h$ Mpc$^{-1}) \leq 3,496$ mK$^2$ at $z = 10.4$, an improvement by a factor of 2.1 and 2.6 respectively. These limits are mostly consistent with thermal noise over a wide range of $k$ after our data quality cuts, despite performing a relatively conservative analysis designed to minimize signal loss. Our results are validated with both statistical tests on the data and end-to-end pipeline simulations. We also report updated constraints on the astrophysics of reionization and the cosmic dawn. Using multiple independent modeling and inference techniques previously employed by HERA Collaboration (2022b), we find that the intergalactic medium must have been heated above the adiabatic cooling limit at least as early as $z = 10.4$, ruling out a broad set of so-called "cold reionization" scenarios. If this heating is due to high-mass X-ray binaries during the cosmic dawn, as is generally believed, our result's 99% credible interval excludes the local relationship between soft X-ray luminosity and star formation and thus requires heating driven by evolved low-metallicity stars. △ Less

Submitted 19 January, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 57 pages, 37 figures. Updated to match the accepted ApJ version. Corresponding author: Joshua S. Dillon

Journal ref: 2023 ApJ 945 124

Utilizing a global network of telescopes to update the ephemeris for the highly eccentric planet HD 80606 b and to ensure the efficient scheduling of JWST

Authors: Kyle A. Pearson, Chas Beichman, Benjamin J. Fulton, Thomas M. Esposito, Robert T. Zellem, David R. Ciardi, Jonah Rolfness, John Engelke, Tamim Fatahi, Rachel Zimmerman-Brachman, Arin Avsar, Dr. Varun Bhalerao, Pat Boyce, Marc Bretton, Alexandra D. Burnett, Jennifer Burt, Martin Fowler, Daniel Gallego, Edward Gomez, Bruno Guillet, Jerry Hilburn, Yves Jongen, Tiffany Kataria, Anastasia Kokori, Harsh Kumar , et al. (15 additional authors not shown)

Abstract: The transiting planet HD80606b undergoes a 1000-fold increase in insolation during its 111-day orbit due to it being highly eccentric (e=0.93). The planet's effective temperature increases from 400K to over 1400K in a few hours as it makes a rapid passage to within 0.03AU of its host star during periapsis. Spectroscopic observations during the eclipse (which is conveniently oriented a few hours be… ▽ More The transiting planet HD80606b undergoes a 1000-fold increase in insolation during its 111-day orbit due to it being highly eccentric (e=0.93). The planet's effective temperature increases from 400K to over 1400K in a few hours as it makes a rapid passage to within 0.03AU of its host star during periapsis. Spectroscopic observations during the eclipse (which is conveniently oriented a few hours before periapsis) of HD80606b with the James Webb Space Telescope (JWST) are poised to exploit this highly variable environment to study a wide variety of atmospheric properties, including composition, chemical and dynamical timescales, and large scale atmospheric motions. Critical to planning and interpreting these observations is an accurate knowledge of the planet's orbit. We report on observations of two full-transit events: 7 February 2020 as observed by the TESS spacecraft and 7--8 December 2021 as observed with a worldwide network of small telescopes. We also report new radial velocity observations which when analyzed with a coupled model to the transits greatly improve the planet's orbital ephemeris. Our new orbit solution reduces the uncertainty in the transit and eclipse timing of the JWST era from tens of minutes to a few minutes. When combined with the planned JWST observations, this new precision may be adequate to look for non-Keplerian effects in the orbit of HD80606b. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: Accepted for publication in AJ; in press

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

Col-OSSOS: Evidence for a compositional gradient inherited from the protoplanetary disk?

Authors: Michael Marsset, Wesley C. Fraser, Megan E. Schwamb, Laura E. Buchanan, Rosemary E. Pike, Nuno Peixinho, Susan Benecchi, Michele T. Bannister, Nicole J. Tan, J. J. Kavelaars

Abstract: In the present-day Kuiper Belt, the number of compositional classes and the orbital distributions of these classes hold important cosmogonic implications for the Solar System. In a companion paper by Fraser et al., we demonstrate that the observed color distribution of small (H>6) Trans-Neptunian Objects (TNOs) can be accounted for by the existence of only two composition classes, named brightIR a… ▽ More In the present-day Kuiper Belt, the number of compositional classes and the orbital distributions of these classes hold important cosmogonic implications for the Solar System. In a companion paper by Fraser et al., we demonstrate that the observed color distribution of small (H>6) Trans-Neptunian Objects (TNOs) can be accounted for by the existence of only two composition classes, named brightIR and faintIR, where the range of colors in each class is governed by a mixture of two material end members. Here, we investigate the orbital distribution of the two color classes identified by Fraser et al. and find that the orbital inclinations of the brightIR class objects are correlated with their optical colors. Using the output of numerical simulations investigating the orbital evolution of TNOs during their scattering phase with Neptune, we show that this correlation could reflect a composition gradient in the early protoplanetary disk, in the range of heliocentric distances over which TNOs from the brightIR class accreted. However, tensions between this interpretation and the existence of blue contaminants among cold classical TNOs, and possible alternative origins for the detected correlation, currently bear uncertainty on our proposed interpretation. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: 19 pages, 4 figures, 5 tables. Submitted, comments welcome!

Col-OSSOS: The Two Types of Kuiper Belt Surfaces

Authors: Wesley C. Fraser, Rosemary E. Pike Michael Marsset, Megan E. Schwamb, Michele T. Bannister, Laura Buchanan, JJ Kavelaars, Susan D. Benecchi, Nicole J. Tan, Nuno Peixinho, Stephen D. J. Gwyn, Ying-Tung Chen, Brett Gladman, Kathryn Volk

Abstract: The Colours of the Outer Solar System Origins Survey (Col-OSSOS) has gathered high quality, near-simultaneous (g-r) and (r-J) colours of 92 Kuiper Belt Objects (KBOs) with (u-g) and (r-z) gathered for some. We present the current state of the survey and data analysis. Recognizing that the optical colours of most icy bodies broadly follow the reddening curve, we present a new projection of the opti… ▽ More The Colours of the Outer Solar System Origins Survey (Col-OSSOS) has gathered high quality, near-simultaneous (g-r) and (r-J) colours of 92 Kuiper Belt Objects (KBOs) with (u-g) and (r-z) gathered for some. We present the current state of the survey and data analysis. Recognizing that the optical colours of most icy bodies broadly follow the reddening curve, we present a new projection of the optical-NIR colours, which rectifies the main non-linear features in the optical-NIR along the ordinates. We find evidence for a bifurcation in the projected colours which presents itself as a diagonal empty region in the optical-NIR. A reanalysis of past colour surveys reveals the same bifurcation. We interpret this as evidence for two separate surface classes: the BrightIR class spans the full range of optical colours and broadly follows the reddening curve, while the FaintIR objects are limited in optical colour, and are less bright in the NIR than the BrightIR objects. We present a two class model. Objects in each class consist of a mix of separate blue and red materials, and span a broad range in colour. Spectra are modelled as linear optical and NIR spectra with different slopes, that intersect at some transition wavelength. The underlying spectral properties of the two classes fully reproduce the observed structures in the UV-optical-NIR colour space ($0.4\lesssimλ\lesssim1.4 \mbox{ $μ$m}$), including the bifurcation observed in the Col-OSSOS and H/WTSOSS datasets, the tendency for cold classical KBOs to have lower (r-z) colours than excited objects, and the well known bimodal optical colour distribution. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: 13 figures, 5 tables, submitted

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

Gas inflows in the polar ring of NGC 4111: the birth of an AGN

Authors: Gabriel R. Hauschild-Roier, Thaisa Storchi-Bergmann, Richard M. McDermid, Jonelle L. Walsh, Joanne Tan, Jonathan Cohn, Davor Krajnović, Jenny Greene, Monica Valluri, Kayhan Gültekin, Sabine Thater, Glenn van de Ven, Karl Gebhardt, Nora Lützgendorf, Benjamin D. Boizelle, Chung-Pei Ma, Aaron J. Barth

Abstract: We have used Hubble Space Telescope (HST) images, SAURON Integral Field Spectroscopy (IFS) and adaptative optics assisted Gemini NIFS near-infrared K-band IFS to map the stellar and gas distribution, excitation and kinematics of the inner few kpc of the nearby edge-on S0 galaxy NGC 4111. The HST images map its $\approx$ 450 pc diameter dusty polar ring, with an estimated gas mass $\ge10^7$ M… ▽ More We have used Hubble Space Telescope (HST) images, SAURON Integral Field Spectroscopy (IFS) and adaptative optics assisted Gemini NIFS near-infrared K-band IFS to map the stellar and gas distribution, excitation and kinematics of the inner few kpc of the nearby edge-on S0 galaxy NGC 4111. The HST images map its $\approx$ 450 pc diameter dusty polar ring, with an estimated gas mass $\ge10^7$ M$_\odot$. The NIFS datacube maps the inner 110 pc radius at $\approx$ 7 pc spatial resolution revealing a $\approx$ 220 pc diameter polar ring in hot ($2267\pm166$ K) molecular H$_2$ 1-0 S(1) gas embedded in the polar ring. The stellar velocity field shows disk-dominated kinematics along the galaxy plane both in the SAURON large-scale and in the NIFS nuclear-scale data. The large-scale [O III] $\lambda5007$ Åvelocity field shows a superposition of two disk kinematics: one similar to that of the stars and another along the polar ring, showing non-circular motions that seem to connect with the velocity field of the nuclear H$_2$ ring, whose kinematics indicate accelerated inflow to the nucleus. The estimated mass inflow rate is enough not only to feed an Active Galactic Nucleus (AGN) but also to trigger circumnuclear star formation in the near future. We propose a scenario in which gas from the polar ring, which probably originated from the capture of a dwarf galaxy, is moving inwards and triggering an AGN, as supported by the local X-ray emission, which seems to be the source of the H$_2$ 1-0 S(1) excitation. The fact that we see neither near-UV nor Br$γ$ emission suggests that the nascent AGN is still deeply buried under the optically thick dust of the polar ring. △ Less

Submitted 18 July, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: 18 pages, 21 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 512, Issue 2, May 2022, Pages 2556-2572

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

A Stacking Survey of Gamma-ray Pulsars

Authors: Yuzhe Song, Timothy A. D. Paglione, Joshua Tan, Charles Lee-Georgescu, Danisbel Herrera

Abstract: We report on a likelihood stacking search for gamma-ray pulsars at 362 high-latitude locations that coincide with known radio pulsar positions. We observe a stacked signal conservatively 2.5$σ$ over the background. Stacking their likelihood profiles in spectral parameter space implies a pulsar-like spectral index and a characteristic flux a factor of 2 below the Fermi-LAT point source sensitivity,… ▽ More We report on a likelihood stacking search for gamma-ray pulsars at 362 high-latitude locations that coincide with known radio pulsar positions. We observe a stacked signal conservatively 2.5$σ$ over the background. Stacking their likelihood profiles in spectral parameter space implies a pulsar-like spectral index and a characteristic flux a factor of 2 below the Fermi-LAT point source sensitivity, assuming isotropic/unbeamed emission from all sample pulsars. The same procedures performed on empty control fields indicate that the pulsars as a population can be distinguished from the background with a $Δ$(TS) = 28 at the peak location (or 4.8$σ$), and the stacked spectra of the control fields are distinctly softer than those of the pulsars. This study also probes a unique region of parameter space populated by low \Edot\ pulsars, most of which have no \gray\ ephemeris available, and is sensitive to high duty cycles. We also discusses the possible \gray\ emission mechanism from such pulsars. △ Less

Submitted 26 July, 2023; v1 submitted 18 December, 2021; originally announced December 2021.

Comments: Accepted for publication to MNRAS

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)