Showing 1–50 of 217 results for author: Law, C

ALMA high-resolution observations unveil planet formation shaping molecular emission in the PDS 70 disk

Authors: L. Rampinelli, S. Facchini, M. Leemker, J. Bae, M. Benisty, R. Teague, C. J. Law, K. I. Öberg, B. Portilla-Revelo, A. J. Cridland

Abstract: With two directly detected protoplanets, the PDS 70 system is a unique source in which to study the complex interplay between forming planets and their natal environment. The large dust cavity carved by the two giant planets can affect the disk chemistry, and therefore the molecular emission morphology. On the other hand, chemical properties of the gas component of the disk are expected to leave a… ▽ More With two directly detected protoplanets, the PDS 70 system is a unique source in which to study the complex interplay between forming planets and their natal environment. The large dust cavity carved by the two giant planets can affect the disk chemistry, and therefore the molecular emission morphology. On the other hand, chemical properties of the gas component of the disk are expected to leave an imprint on the planetary atmospheres. In this work, we reconstruct the emission morphology of a rich inventory of molecular tracers in the PDS 70 disk, and we look for possible chemical signatures of the two actively accreting protoplanets, PDS b and c. We leverage Atacama Large Millimeter/submillimeter Array (ALMA) band 6 high-angular-resolution and deep-sensitivity line emission observations, together with image and $uv$-plane techniques, to boost the detection of faint lines. We robustly detect ring-shaped emission from $^{12}$CO, $^{13}$CO, C$^{18}$O, H$^{13}$CN, HC$^{15}$N, DCN, H$_2$CO, CS, C$_2$H, and H$^{13}$CO$^{+}$ lines in unprecedented detail. Most of the molecular tracers show a peak of the emission inside the millimeter dust peak. We interpret this as the direct impact of the effective irradiation of the cavity wall, as a result of the planet formation process. Moreover, we have found evidence of an O-poor gas reservoir in the outer disk, which is supported by the observations of bright C-rich molecules, the non-detection of SO, and a lower limit on the $\mathrm{CS/SO}$ ratio of $\sim1$. Eventually, we provide the first detection of the c-C$_3$H$_2$ transitions at 218.73 GHz, and the marginal detection of an azimuthal asymmetry in the higher-energy H$_2$CO (3$_{2,1}$-2$_{2,0}$) line, which could be due to accretion heating near PDS 70b. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: 25 pages, 10 figures, 3 tables, 8 figures and one table in appendix. Accepted for publication in A&A

A Heavily Scattered Fast Radio Burst Is Viewed Through Multiple Galaxy Halos

Authors: Jakob T. Faber, Vikram Ravi, Stella Koch Ocker, Myles B. Sherman, Kritti Sharma, Liam Connor, Casey Law, Nikita Kosogorov, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, James W. Lamb, Paul Rasmussen, Jean J. Somalwar, Sander Weinreb, David P. Woody

Abstract: We present a multi-wavelength study of the apparently non-repeating, heavily scattered fast radio burst, FRB 20221219A, detected by the Deep Synoptic Array 110 (DSA-110). The burst exhibits a moderate dispersion measure (DM) of $706.7^{+0.6}_{-0.6}$ $\mathrm{pc}~\mathrm{cm}^{-3}$ and an unusually high scattering timescale of $τ_{\mathrm{obs}} = 19.2_{-2.7}^{+2.7}$ ms at 1.4 GHz. We associate the F… ▽ More We present a multi-wavelength study of the apparently non-repeating, heavily scattered fast radio burst, FRB 20221219A, detected by the Deep Synoptic Array 110 (DSA-110). The burst exhibits a moderate dispersion measure (DM) of $706.7^{+0.6}_{-0.6}$ $\mathrm{pc}~\mathrm{cm}^{-3}$ and an unusually high scattering timescale of $τ_{\mathrm{obs}} = 19.2_{-2.7}^{+2.7}$ ms at 1.4 GHz. We associate the FRB with a Milky Way-like host galaxy at $z_{\mathrm{host}} = 0.554$ of stellar mass $\mathrm{log}_{10}(M_{\star, \mathrm{host}}) = 10.20^{+0.04}_{-0.03} ~M_\odot$. We identify two intervening galaxy halos at redshifts $z_{\mathrm{igh1}} = 0.492$ and $z_{\mathrm{igh2}} = 0.438$, with low impact parameters, $b_{\mathrm{igh1}} = 43.0_{-11.3}^{+11.3}$ kpc and $b_{\mathrm{igh2}} = 36.1_{-11.3}^{+11.3}$ kpc, and intermediate stellar masses, $\mathrm{log}_{10}(M_{\star, \mathrm{igh1}}) = 10.01^{+0.02}_{-0.02} ~M_\odot$ and $\mathrm{log}_{10}(M_{\star, \mathrm{igh2}}) = 10.60^{+0.02}_{-0.02} ~M_\odot$. The presence of two such galaxies suggests that the sightline is significantly overcrowded compared to the median sightline to this redshift, as inferred from the halo mass function. We perform a detailed analysis of the sightline toward FRB 20221219A, constructing both DM and scattering budgets. Our results suggest that, unlike most well-localized sources, the host galaxy does not dominate the observed scattering. Instead, we posit that an intersection with a single partially ionized cloudlet in the circumgalactic medium of an intervening galaxy could account for the substantial scattering in FRB 20221219A and remain in agreement with typical electron densities inferred for extra-planar dense cloud-like structures in the Galactic and extragalactic halos (e.g., high-velocity clouds). △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 18 pages, 6 figures, submitted to ApJ, comments appreciated

Neutral hydrogen filaments in interstellar media: Are they physical?

Authors: Ka Ho Yuen, Ka Wai Ho, Chi Yan Law, Avi Chen

Abstract: The trending term "filament" is extensively used in the interstellar medium (ISM) and the star formation community, and is believed to be one of the most important objects that gauge molecular cloud and star formation. However, the physical definition of these ubiquitous, elongated, high contrast features is poorly defined and still actively debated. Despite the absence of a unified consensus, fil… ▽ More The trending term "filament" is extensively used in the interstellar medium (ISM) and the star formation community, and is believed to be one of the most important objects that gauge molecular cloud and star formation. However, the physical definition of these ubiquitous, elongated, high contrast features is poorly defined and still actively debated. Despite the absence of a unified consensus, filaments are believed to be involved in many important physical processes from galaxy structure formation to the emergence of protostellar objects. Therefore, understanding how filaments form, what constrains their growth, and their general physical properties, are extremely important for theorists and observers who study the dynamics of the ISM and consequent star formations. This review serves as a collection of the community's views and develops the concept of "filaments" in the context of the ISM and star-forming clouds. Observationally, filaments are seen across the entire sky and often carry an aspect ratio of the order of hundreds. In the context of the ISM, filaments are believed to form by stretching and tearing from magnetized ISM turbulence. ISM filaments are subjected to heating and cooling phases, and are likely to be magnetically aligned. Cold clouds are formed inside ISM due to turbulence instability. This review updates the understanding of ISM filaments in the community. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Journal ref: Rev. Mod. Plasma Phys. 8, 21 (2024)

Self-similarity of the magnetic field at different scales: the case of G31.41+0.31

Authors: M. T. Beltrán, M. Padovani, D. Galli, N. Áñez-López, J. M. Girart, R. Cesaroni, D. Dall'Olio, G. Anglada, C. Y. Law, A. Lorenzani, L. Moscadelli, Á. Sánchez-Monge, M. Osorio, Q. Zhang

Abstract: Context. Dust polarization observations of the massive protocluster G31.41+0.31 carried out at ~1'' (~3750 au) resolution with the SMA at 870 microm have revealed one of the clearest examples to date of an hourglass-shaped magnetic field morphology in the high-massregime. Additionally, ~0.24'' (~900 au) resolution observations with ALMA at 1.3 mm have confirmed these results. The next step is to i… ▽ More Context. Dust polarization observations of the massive protocluster G31.41+0.31 carried out at ~1'' (~3750 au) resolution with the SMA at 870 microm have revealed one of the clearest examples to date of an hourglass-shaped magnetic field morphology in the high-massregime. Additionally, ~0.24'' (~900 au) resolution observations with ALMA at 1.3 mm have confirmed these results. The next step is to investigate whether the magnetic field maintains its hourglass-shaped morphology down to circumstellar scales. Aims. To study the magnetic field morphology toward the four (proto)stars A, B, C, and D contained in G31.41+0.31 and examine whether the self-similarity observed at core scales (1'' and 0.24'' resolution) still holds at circumstellar scales, we carried out ALMA observations of the polarized dust continuum emission at 1.3 mm and 3.1 mm at an angular resolution of ~0.068'' (~250 au), sufficient to resolve the envelope emission of the embedded protostars. Methods. We used ALMA to perform full polarization observations at 233 GHz (Band 6) and 97.5 GHz (Band 3) with a synthesized beam of 0.072'' x 0.064''. We carried out polarization observations at two different wavelengths to confirm that the polarization traces magnetically aligned dust grains and is not due to dust self-scattering. Results. The polarized emission and the direction of the magnetic field obtained at the two wavelengths are basically the same, except for an area between the embedded sources C and B. In such an area, the emission at 1.3 mm could be optically thick and affected by dichroic extinction. ... △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: 15 pages, 15 figures, Accepted for publication in A&A

Anatomy of the Class I protostar L1489 IRS with NOEMA -- I. Disk, streamers, outflow(s) and bubbles at 3mm

Authors: M. Tanious, R. Le Gal, R. Neri, A. Faure, A. Gupta, C. J. Law, J. Huang, N. Cuello, J. P. Williams, F. Ménard

Abstract: Over the past few years, chemical studies have revealed multiple structures in the vicinity of young stellar objects (YSOs). It has become evident that specific physical conditions are associated with the emission of particular molecular lines, allowing us to use molecular probes of the YSO physics. Consequently, chemical surveys are now necessary to fully constrain the origin of the observed stru… ▽ More Over the past few years, chemical studies have revealed multiple structures in the vicinity of young stellar objects (YSOs). It has become evident that specific physical conditions are associated with the emission of particular molecular lines, allowing us to use molecular probes of the YSO physics. Consequently, chemical surveys are now necessary to fully constrain the origin of the observed structures. Several surveys have been conducted to explore the chemistry of YSOs, focusing on Class 0 and Class II objects. However, our knowledge of intermediate objects, that are Class I objects, remains limited. To bridge the gap and establish the relationship between observed structures and molecular line emission at the Class I evolutionary stage, we investigate the spatial distribution of key molecular gas species in the low-mass Class I protostar L1489 IRS (IRAS 04016+2610), a source part of the ChemYSO survey. We performed a 3mm line survey at high spatial and high spectral resolution using the NOEMA interferometer and the IRAM-30m telescope. We present here the ten brightest lines of our survey, in which we identified a new ~ 3 000 au long streamer in HC3N, C2H, and c-C3H2 emission, likely associated with more localized accretion shocks probed in SO. In addition, two ~ 10 000 au bubbles are seen with the dense molecular tracers HCO+, CS, and HCN around the YSO. Additionally, potential indicators of a second outflow appear in CS and HCN emission, but its nature remains to be confirmed. The late infall identified at large scales may originate from the nearby prestellar core L1489 and is likely responsible for the formation of an external warped disk in this system. The detection of a potential second outflow could be the direct evidence of a binary system. Finally, we hypothesize that the bubbles may result from the magnetic pressure as observed in numerical simulations. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: 16 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

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

The First Spatially-resolved Detection of $^{13}$CN in a Protoplanetary Disk and Evidence for Complex Carbon Isotope Fractionation

Authors: Tomohiro C. Yoshida, Hideko Nomura, Kenji Furuya, Richard Teague, Charles J. Law, Takashi Tsukagoshi, Seokho Lee, Christian Rab, Karin I. Öberg, Ryan A. Loomis

Abstract: Recent measurements of carbon isotope ratios in both protoplanetary disks and exoplanet atmospheres have suggested a possible transfer of significant carbon isotope fractionation from disks to planets. For a clearer understanding of the isotopic link between disks and planets, it is important to measure the carbon isotope ratios in various species. In this paper, we present a detection of the… ▽ More Recent measurements of carbon isotope ratios in both protoplanetary disks and exoplanet atmospheres have suggested a possible transfer of significant carbon isotope fractionation from disks to planets. For a clearer understanding of the isotopic link between disks and planets, it is important to measure the carbon isotope ratios in various species. In this paper, we present a detection of the $^{13}$CN $N=2-1$ hyperfine lines in the TW Hya disk with the Atacama Large Millimeter/submillimeter Array. This is the first spatially-resolved detection of $^{13}$CN in disks, which enables us to measure the spatially resolved $^{12}$CN/$^{13}$CN ratio for the first time. We conducted non-local thermal equilibrium modeling of the $^{13}$CN lines in conjunction with previously observed $^{12}$CN lines to derive the kinetic temperature, ${\rm H_2}$ volume density, and column densities of $^{12}$CN and $^{13}$CN. The ${\rm H_2}$ volume density is found to range between $ (4 - 10)\times10^7 \ {\rm cm^{-3}}$, suggesting that CN molecules mainly reside in the disk upper layer. The $^{12}$CN/$^{13}$CN ratio is measured to be $ 70^{+9}_{-6}$ at $30 < r < 80$ au from the central star, which is similar to the $\rm ^{12}C/^{13}C$ ratio in the interstellar medium. However, this value differs from the previously reported values found for other carbon-bearing molecules (CO and HCN) in the TW Hya disk. This could be self-consistently explained by different emission layer heights for different molecules combined with preferential sequestration of $\rm ^{12}C$ into the solid phase towards the disk midplane. This study reveals the complexity of the carbon isotope fractionation operating in disks. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 16 pages, 7 figures, accepted for publication in ApJ

Simultaneous multi-wavelength observations of the repeating fast radio burst FRB 20190520B with Swift and FAST

Authors: Zhen Yan, Wenfei Yu, Kim L. Page, Jie Lin, Di Li, Chenhui Niu, Casey Law, Bing Zhang, Shami Chatterjee, Xian Zhang, Reshma Anna-Thomas

Abstract: Fast radio bursts (FRBs) are bright, millisecond-duration radio bursts of cosmic origin. There have been several dozen FRBs found to repeat. Among them, those precisely localized provide the best opportunity to probe their multi-wavelength counterparts, local environment, and host galaxy that would reveal their origins. Here we report our X-ray, ultraviolet (UV) and optical observations with the… ▽ More Fast radio bursts (FRBs) are bright, millisecond-duration radio bursts of cosmic origin. There have been several dozen FRBs found to repeat. Among them, those precisely localized provide the best opportunity to probe their multi-wavelength counterparts, local environment, and host galaxy that would reveal their origins. Here we report our X-ray, ultraviolet (UV) and optical observations with the $Swift$ satellite that were performed simultaneously in the radio band with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) observations of the repeating FRB 20190520B, aiming at detection of possible multi-wavelength bursts in association with radio bursts and multi-wavelength counterpart of the persistent radio source (PRS). While a total of 10 radio bursts were detected by FAST at the same time of $Swift$ observations, we detected neither X-ray, UV or optical bursts in accompany of the radio bursts, nor persistent multi-wavelength counterpart of the PRS. We obtained the energy upper limits ($3σ$) on any multi-wavelength bursts as $5.03 \times 10^{47}$ erg in the hard X-ray band (15-150 keV), $7.98 \times 10^{45}$ erg in the soft X-ray band (0.3-10 keV), and $4.51 \times 10^{44}$ erg in the U band, respectively. The energy ratio between soft X-ray (0.3-10 keV) and radio emission of the bursts is constrained as $<6\times10^{7}$, and the ratio between optical (U band) and radio as $<1.19\times10^{6}$. The 3$σ$ luminosity upper limits at the position of PRS are 1.04$\times10^{47}$ (15-150 keV), 8.81$\times10^{42}$ (0.3-10 keV), 9.26$\times10^{42}$ (UVW1), and 2.54$\times10^{42}$ erg s$^{-1}$ (U), respectively. We show that the PRS is much more radio loud than representative pulsar wind nebulae, supernova remnants, extended jet of Galactic X-ray binaries and ultraluminous X-ray sources, suggestive of boosted radio emission of the PRS. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 12 pages, 7 figures, submitted to ApJ

An ALMA molecular inventory of warm Herbig Ae disks: II. Abundant complex organics and volatile sulphur in the IRS 48 disk

Authors: Alice S. Booth, Milou Temmink, Ewine F. van Dishoeck, Lucy Evans, John D. Ilee, Mihkel Kama, Luke Keyte, Charles J. Law, Margot Leemker, Nienke van der Marel, Hideko Nomura, Shota Notsu, Karin Öberg, Catherine Walsh

Abstract: The Atacama Large Millimeter/submillimeter Array (ALMA) can probe the molecular content of planet-forming disks with unprecedented sensitivity. These observations allow us to build up an inventory of the volatiles available for forming planets and comets. Herbig Ae transition disks are fruitful targets due to the thermal sublimation of complex organic molecule (COM) and likely H2O-rich ices in the… ▽ More The Atacama Large Millimeter/submillimeter Array (ALMA) can probe the molecular content of planet-forming disks with unprecedented sensitivity. These observations allow us to build up an inventory of the volatiles available for forming planets and comets. Herbig Ae transition disks are fruitful targets due to the thermal sublimation of complex organic molecule (COM) and likely H2O-rich ices in these disks. The IRS 48 disk shows a particularly rich chemistry that can be directly linked to its asymmetric dust trap. Here, we present ALMA observations of the IRS 48 disk where we detect 16 different molecules and make the first robust detections of H213CO, 34SO, 33SO and c-H2COCH2 (ethylene oxide) in a protoplanetary disk. All of the molecular emissions, aside from CO, are colocated with the dust trap and this includes newly detected simple molecules such as HCO+, HCN and CS. Interestingly, there are spatial offsets between different molecular families, including between the COMs and sulphur-bearing species, with the latter being more azimuthally extended and located radially further from the star. The abundances of the newly detected COMs relative to CH3OH are higher than the expected protostellar ratios, which implies some degree of chemical processing of the inherited ices during the disk lifetime. These data highlight IRS 48 as a unique astrochemical laboratory to unravel the full volatile reservoir at the epoch of planet and comet formation and the role of the disk in (re)setting chemical complexity. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: Accepted to AJ, 21 pages, 7 figures

An ALMA molecular inventory of warm Herbig Ae disks: I. Molecular rings, asymmetries and complexity in the HD 100546 disk

Authors: Alice S. Booth, Margot Leemker, Ewine F. van Dishoeck, Lucy Evans, John D. Ilee, Mihkel Kama, Luke Keyte, Charles J. Law, Nienke van der Marel, Hideko Nomura, Shota Notsu, Karin Öberg, Milou Temmink, Catherine Walsh

Abstract: Observations of disks with the Atacama Large Millimeter/submillimeter Array (ALMA) allow us to map the chemical makeup of nearby protoplanetary disks with unprecedented spatial resolution and sensitivity. The typical outer Class II disk observed with ALMA is one with an elevated C/O ratio and a lack of oxygen-bearing complex organic molecules, but there are now some interesting exceptions: three t… ▽ More Observations of disks with the Atacama Large Millimeter/submillimeter Array (ALMA) allow us to map the chemical makeup of nearby protoplanetary disks with unprecedented spatial resolution and sensitivity. The typical outer Class II disk observed with ALMA is one with an elevated C/O ratio and a lack of oxygen-bearing complex organic molecules, but there are now some interesting exceptions: three transition disks around Herbig Ae stars all show oxygen-rich gas traced via the unique detections of the molecules SO and CH3OH. We present the first results of an ALMA line survey at 337 to 357 GHz of such disks and focus this paper on the first Herbig Ae disk to exhibit this chemical signature - HD 100546. In these data, we detect 19 different molecules including NO, SO and CH3OCHO (methyl formate). We also make the first tentative detections of H213CO and 34SO in protoplanetary disks. Multiple molecular species are detected in rings, which are, surprisingly, all peaking just beyond the underlying millimeter continuum ring at 200 au. This result demonstrates a clear connection between the large dust distribution and the chemistry in this flat outer disk. We discuss the physical and/or chemical origin of these sub-structures in relation to ongoing planet formation in the HD 100546 disk. We also investigate how similar and/or different the molecular make up of this disk is to other chemically well-characterised Herbig Ae disks. The line-rich data we present motivates the need for more ALMA line surveys to probe the observable chemistry in Herbig Ae systems which offer unique insight into the composition of disk ices, including complex organic molecules. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: Accepted to AJ, 25 pages, 11 figures

JWST-MIRI Spectroscopy of Warm Molecular Emission and Variability in the AS 209 Disk

Authors: Carlos E. Muñoz-Romero, Karin I. Öberg, Andrea Banzatti, Klaus M. Pontoppidan, Sean M. Andrews, David J. Wilner, Edwin A. Bergin, Ian Czekala, Charles J. Law, Colette Salyk, Richard Teague, Chunhua Qi, Jennifer B. Bergner, Jane Huang, Catherine Walsh, Viviana V. Guzmán, L. Ilsedore Cleeves, Yuri Aikawa, Jaehan Bae, Alice S. Booth, Gianni Cataldi, John D. Ilee, Romane Le Gal, Feng Long, Ryan A. Loomis , et al. (2 additional authors not shown)

Abstract: We present MIRI MRS observations of the large, multi-gapped protoplanetary disk around the T-Tauri star AS 209. The observations reveal hundreds of water vapor lines from 4.9 to 25.5 $μ$m towards the inner $\sim1$ au in the disk, including the first detection of ro-vibrational water emission in this disk. The spectrum is dominated by hot ($\sim800$ K) water vapor and OH gas, with only marginal det… ▽ More We present MIRI MRS observations of the large, multi-gapped protoplanetary disk around the T-Tauri star AS 209. The observations reveal hundreds of water vapor lines from 4.9 to 25.5 $μ$m towards the inner $\sim1$ au in the disk, including the first detection of ro-vibrational water emission in this disk. The spectrum is dominated by hot ($\sim800$ K) water vapor and OH gas, with only marginal detections of CO$_2$, HCN, and a possible colder water vapor component. Using slab models with a detailed treatment of opacities and line overlap, we retrieve the column density, emitting area, and excitation temperature of water vapor and OH, and provide upper limits for the observable mass of other molecules. Compared to MIRI spectra of other T-Tauri disks, the inner disk of AS 209 does not appear to be atypically depleted in CO$_2$ nor HCN. Based on \textit{Spitzer IRS} observations, we further find evidence for molecular emission variability over a 10-year baseline. Water, OH, and CO$_2$ line luminosities have decreased by factors 2-4 in the new MIRI epoch, yet there are minimal continuum emission variations. The origin of this variability is yet to be understood. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJ

Filamentary Network and Magnetic Field Structures Revealed with BISTRO in the High-Mass Star-Forming Region NGC2264 : Global Properties and Local Magnetogravitational Configurations

Authors: Jia-Wei Wang, Patrick M. Koch, Seamus D. Clarke, Gary Fuller, Nicolas Peretto, Ya-Wen Tang, Hsi-Wei Yen, Shih-Ping Lai, Nagayoshi Ohashi, Doris Arzoumanian, Doug Johnstone, Ray Furuya, Shu-ichiro Inutsuka, Chang Won Lee, Derek Ward-Thompson, Valentin J. M. Le Gouellec, Hong-Li Liu, Lapo Fanciullo, Jihye Hwang, Kate Pattle, Frédérick Poidevin, Mehrnoosh Tahani, Takashi Onaka, Mark G. Rawlings, Eun Jung Chung , et al. (132 additional authors not shown)

Abstract: We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from… ▽ More We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the Astrophysical Journal. 43 pages, 32 figures, and 4 tables (including Appendix)

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

Mapping the Vertical Gas Structure of the Planet-hosting PDS 70 Disk

Authors: Charles J. Law, Myriam Benisty, Stefano Facchini, Richard Teague, Jaehan Bae, Andrea Isella, Inga Kamp, Karin I. Öberg, Bayron Portilla-Revelo, Luna Rampinelli

Abstract: PDS 70 hosts two massive, still-accreting planets and the inclined orientation of its protoplanetary disk presents a unique opportunity to directly probe the vertical gas structure of a planet-hosting disk. Here, we use high-spatial-resolution (${\approx}$0."1;10 au) observations in a set of CO isotopologue lines and HCO$^+$ J=4-3 to map the full 2D $(r,z)$ disk structure from the disk atmosphere,… ▽ More PDS 70 hosts two massive, still-accreting planets and the inclined orientation of its protoplanetary disk presents a unique opportunity to directly probe the vertical gas structure of a planet-hosting disk. Here, we use high-spatial-resolution (${\approx}$0."1;10 au) observations in a set of CO isotopologue lines and HCO$^+$ J=4-3 to map the full 2D $(r,z)$ disk structure from the disk atmosphere, as traced by $^{12}$CO, to closer to the midplane, as probed by less abundant isotopologues and HCO$^+$. In the PDS 70 disk, $^{12}$CO traces a height of $z/r\approx0.3$, $^{13}$CO is found at $z/r\approx0.1$, and C$^{18}$O originates at, or near, the midplane. The HCO$^+$ surface arises from $z/r\approx0.2$ and is one of the few non-CO emission surfaces constrained with high fidelity in disks to date. In the $^{12}$CO J=3-2 line, we resolve a vertical dip and steep rise in height at the cavity wall, making PDS 70 the first transition disk where this effect is directly seen in line emitting heights. In the outer disk, the CO emission heights of PDS 70 appear typical for its stellar mass and disk size and are not substantially altered by the two inner embedded planets. By combining CO isotopologue and HCO$^+$ lines, we derive the 2D gas temperature structure and estimate a midplane CO snowline of ${\approx}$56-85 au. This implies that both PDS 70b and 70c are located interior to the CO snowline and are likely accreting gas with a high C/O ratio of ${\approx}$1.0, which provides context for future planetary atmospheric measurements from, e.g., JWST, and for properly modeling their formation histories. △ Less

Submitted 5 January, 2024; originally announced January 2024.

Comments: 24 pages, 14 figures, accepted for publication in ApJ

An unidentified Fermi source emitting radio bursts in the Galactic bulge

Authors: Reshma Anna-Thomas, Sarah Burke-Spolaor, Casey J. Law, F. K. Schinzel, Kshitij Aggarwal, Geoffrey C. Bower, Liam Connor, Paul B. Demorest

Abstract: We report on the detection of radio bursts from the Galactic bulge using the real-time transient detection and localization system, realfast. The pulses were detected commensally on the Karl G. Jansky Very Large Array during a survey of unidentified Fermi $γ$-ray sources. The bursts were localized to subarcsecond precision using realfast fast-sampled imaging. Follow-up observations with the Green… ▽ More We report on the detection of radio bursts from the Galactic bulge using the real-time transient detection and localization system, realfast. The pulses were detected commensally on the Karl G. Jansky Very Large Array during a survey of unidentified Fermi $γ$-ray sources. The bursts were localized to subarcsecond precision using realfast fast-sampled imaging. Follow-up observations with the Green Bank Telescope detected additional bursts from the same source. The bursts do not exhibit periodicity in a search up to periods of 480s, assuming a duty cycle of < 20%. The pulses are nearly 100% linearly polarized, show circular polarization up to 12%, have a steep radio spectral index of -2.7, and exhibit variable scattering on timescales of months. The arcsecond-level realfast localization links the source confidently with the Fermi $γ$-ray source and places it nearby (though not coincident with) an XMM-Newton X-ray source. Based on the source's overall properties, we discuss various options for the nature of this object and propose that it could be a young pulsar, magnetar, or a binary pulsar system. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: Submitted to The Astrophysical Journal

C I Traces the Disk Atmosphere in the IM Lup Protoplanetary Disk

Authors: Charles J. Law, Felipe Alarcón, L. Ilsedore Cleeves, Karin I. Öberg, Teresa Paneque-Carreño

Abstract: The central star and its energetic radiation fields play a vital role in setting the vertical and radial chemical structure of planet-forming disks. We present observations that, for the first time, clearly reveal the UV-irradiated surface of a protoplanetary disk. Specifically, we spatially resolve the atomic-to-molecular (C I-to-CO) transition in the IM Lup disk with ALMA archival observations o… ▽ More The central star and its energetic radiation fields play a vital role in setting the vertical and radial chemical structure of planet-forming disks. We present observations that, for the first time, clearly reveal the UV-irradiated surface of a protoplanetary disk. Specifically, we spatially resolve the atomic-to-molecular (C I-to-CO) transition in the IM Lup disk with ALMA archival observations of [C I] $^3$P$_1$-$^3$P$_0$. We derive a C I emitting height of z/r $\gtrsim$ 0.5 with emission detected out to a radius of ${\approx}$600 au. Compared to other systems with C I heights inferred from unresolved observations or models, the C I layer in the IM Lup disk is at scale heights almost double that of other disks, confirming its highly flared nature. C I arises from a narrow, optically-thin layer that is substantially more elevated than that of $^{12}$CO (z/r $\approx$ 0.3-0.4), which allows us to directly constrain the physical gas conditions across the C I-to-CO transition zone. We also compute a radially-resolved C I column density profile and find a disk-averaged C I column density of 2$\times10^{16}$ cm$^{-2}$, which is ${\approx}$3-20$\times$ lower than that of other disks with spatially-resolved C I detections. We do not find evidence for vertical substructures or spatially-localized deviations in C I due, e.g., to either an embedded giant planet or a photoevaporative wind that have been proposed in the IM Lup disk, but emphasize that deeper observations are required for robust constraints. △ Less

Submitted 13 December, 2023; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: 12 pages, 6 figures, accepted for publication in ApJL; references added

The Future of Astronomical Data Infrastructure: Meeting Report

Authors: Michael R. Blanton, Janet D. Evans, Dara Norman, William O'Mullane, Adrian Price-Whelan, Luca Rizzi, Alberto Accomazzi, Megan Ansdell, Stephen Bailey, Paul Barrett, Steven Berukoff, Adam Bolton, Julian Borrill, Kelle Cruz, Julianne Dalcanton, Vandana Desai, Gregory P. Dubois-Felsmann, Frossie Economou, Henry Ferguson, Bryan Field, Dan Foreman-Mackey, Jaime Forero-Romero, Niall Gaffney, Kim Gillies, Matthew J. Graham , et al. (47 additional authors not shown)

Abstract: The astronomical community is grappling with the increasing volume and complexity of data produced by modern telescopes, due to difficulties in reducing, accessing, analyzing, and combining archives of data. To address this challenge, we propose the establishment of a coordinating body, an "entity," with the specific mission of enhancing the interoperability, archiving, distribution, and productio… ▽ More The astronomical community is grappling with the increasing volume and complexity of data produced by modern telescopes, due to difficulties in reducing, accessing, analyzing, and combining archives of data. To address this challenge, we propose the establishment of a coordinating body, an "entity," with the specific mission of enhancing the interoperability, archiving, distribution, and production of both astronomical data and software. This report is the culmination of a workshop held in February 2023 on the Future of Astronomical Data Infrastructure. Attended by 70 scientists and software professionals from ground-based and space-based missions and archives spanning the entire spectrum of astronomical research, the group deliberated on the prevailing state of software and data infrastructure in astronomy, identified pressing issues, and explored potential solutions. In this report, we describe the ecosystem of astronomical data, its existing flaws, and the many gaps, duplication, inconsistencies, barriers to access, drags on productivity, missed opportunities, and risks to the long-term integrity of essential data sets. We also highlight the successes and failures in a set of deep dives into several different illustrative components of the ecosystem, included as an appendix. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: 59 pages; please send comments and/or questions to foadi@googlegroups.com

ALMA Observations of Supernova Remnant N49 in the Large Magellanic Cloud. II. Non-LTE Analysis of Shock-heated Molecular Clouds

Authors: H. Sano, Y. Yamane, J. Th. van Loon, K. Furuya, Y. Fukui, R. Z. E. Alsaberi, A. Bamba, R. Enokiya, M. D. Filipović, R. Indebetouw, T. Inoue, A. Kawamura, M. Lakićević, C. J. Law, N. Mizuno, T. Murase, T. Onishi, S. Park, P. P. Plucinsky, J. Rho, A. M. S. Richards, G. Rowell, M. Sasaki, J. Seok, P. Sharda , et al. (6 additional authors not shown)

Abstract: We present the first compelling evidence of shock-heated molecular clouds associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC). Using $^{12}$CO($J$ = 2-1, 3-2) and $^{13}$CO($J$ = 2-1) line emission data taken with the Atacama Large Millimeter/Submillimeter Array, we derived the H$_2$ number density and kinetic temperature of eight $^{13}$CO-detected clouds using the… ▽ More We present the first compelling evidence of shock-heated molecular clouds associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC). Using $^{12}$CO($J$ = 2-1, 3-2) and $^{13}$CO($J$ = 2-1) line emission data taken with the Atacama Large Millimeter/Submillimeter Array, we derived the H$_2$ number density and kinetic temperature of eight $^{13}$CO-detected clouds using the large velocity gradient approximation at a resolution of 3.5$''$ (~0.8 pc at the LMC distance). The physical properties of the clouds are divided into two categories: three of them near the shock front show the highest temperatures of ~50 K with densities of ~500-700 cm$^{-3}$, while other clouds slightly distant from the SNR have moderate temperatures of ~20 K with densities of ~800-1300 cm$^{-3}$. The former clouds were heated by supernova shocks, but the latter were dominantly affected by the cosmic-ray heating. These findings are consistent with the efficient production of X-ray recombining plasma in N49 due to thermal conduction between the cold clouds and hot plasma. We also find that the gas pressure is roughly constant except for the three shock-engulfed clouds inside or on the SNR shell, suggesting that almost no clouds have evaporated within the short SNR age of ~4800 yr. This result is compatible with the shock-interaction model with dense and clumpy clouds inside a low-density wind bubble. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: 12 pages, 5 figures, 1 table, accepted for publication in The Astrophysical Journal (ApJ)

VLASS tidal disruption events with optical flares I: the sample and a comparison to optically-selected TDEs

Authors: Jean J. Somalwar, Vikram Ravi, Dillon Z. Dong, Erica Hammerstein, Gregg Hallinan, Casey Law, Jessie Miller, Steven T. Myers, Yuhan Yao, Richard Dekany, Matthew Graham, Steven L. Groom, Josiah Purdum, Avery Wold

Abstract: In this work, we use the Jansky VLA Sky Survey (VLASS) to compile the first sample of six radio-selected tidal disruption events (TDEs) with transient optical counterparts. While we still lack the statistics to do detailed population studies of radio-selected TDEs, we use these events to suggest trends in host galaxy and optical light curve properties that may correlate with the presence of radio… ▽ More In this work, we use the Jansky VLA Sky Survey (VLASS) to compile the first sample of six radio-selected tidal disruption events (TDEs) with transient optical counterparts. While we still lack the statistics to do detailed population studies of radio-selected TDEs, we use these events to suggest trends in host galaxy and optical light curve properties that may correlate with the presence of radio emission, and hence can inform optically-selected TDE radio follow-up campaigns. We find that radio-selected TDEs tend to have faint and cool optical flares, as well as host galaxies with low SMBH masses. Our radio-selected TDEs also tend to have more energetic, larger radio emitting regions than radio-detected, optically-selected TDEs. We consider possible explanations for these trends, including by invoking super-Eddington accretion and enhanced circumnuclear media. Finally, we constrain the radio-emitting TDE rate to be $\gtrsim 10$ Gpc$^{-3}$ yr$^{-1}$. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 26 pages, 5 tables, 11 figures, submitted to ApJ

A JWST inventory of protoplanetary disk ices: The edge-on protoplanetary disk HH 48 NE, seen with the Ice Age ERS program

Authors: J. A. Sturm, M. K. McClure, T. L. Beck, D. Harsono, J. B. Bergner, E. Dartois, A. C. A. Boogert, J. E. Chiar, M. A. Cordiner, M. N. Drozdovskaya, S. Ioppolo, C. J. Law, H. Linnartz, D. C. Lis, G. J. Melnick, B. A. McGuire, J. A. Noble, K. I. Öberg, M. E. Palumbo, Y. J. Pendleton, G. Perotti, K. M. Pontoppidan, D. Qasim, W. R. M. Rocha, H. Terada , et al. (2 additional authors not shown)

Abstract: Ices are the main carriers of volatiles in protoplanetary disks and are crucial to our understanding of the chemistry that ultimately sets the organic composition of planets. The ERS program Ice Age on the JWST follows the ice evolution through all stages of star and planet formation. JWST/NIRSpec observations of the edge-on Class II protoplanetary disk HH~48~NE reveal spatially resolved absorptio… ▽ More Ices are the main carriers of volatiles in protoplanetary disks and are crucial to our understanding of the chemistry that ultimately sets the organic composition of planets. The ERS program Ice Age on the JWST follows the ice evolution through all stages of star and planet formation. JWST/NIRSpec observations of the edge-on Class II protoplanetary disk HH~48~NE reveal spatially resolved absorption features of the major ice components H$_2$O, CO$_2$, CO, and multiple weaker signatures from less abundant ices NH$_3$, OCN$^-$, and OCS. Isotopologue $^{13}$CO$_2$ ice has been detected for the first time in a protoplanetary disk. Since multiple complex light paths contribute to the observed flux, the ice absorption features are filled in by ice-free scattered light. The $^{12}$CO$_2$/$^{13}$CO$_2$ ratio of 14 implies that the $^{12}$CO$_2$ feature is saturated, without the flux approaching 0, indicative of a very high CO$_2$ column density on the line of sight, and a corresponding abundance with respect to hydrogen that is higher than ISM values by a factor of at least a few. Observations of rare isotopologues are crucial, as we show that the $^{13}$CO$_2$ observation allows us to determine the column density of CO$_2$ to be at an order of magnitude higher than the lower limit directly inferred from the observed optical depth. Radial variations in ice abundance, e.g., snowlines, are significantly modified since all observed photons have passed through the full radial extent of the disk. CO ice is observed at perplexing heights in the disk, extending to the top of the CO-emitting gas layer. We argue that the most likely interpretation is that we observe some CO ice at high temperatures, trapped in less volatile ices like H$_2$O and CO$_2$. Future radiative transfer models will be required to constrain the implications on our current understanding of disk physics and chemistry. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 16 pages, 8 figures, accepted for publication in A&A

Improving Undergraduate Astronomy Students' Skills with Research Literature via Accessible Summaries: A Case Study with Astrobites-based Lesson Plans

Authors: Briley L. Lewis, Abygail R. Waggoner, Emma Clarke, Alison L. Crisp, Mark Dodici, Graham M. Doskoch, Michael M. Foley, Ryan Golant, Katya Gozman, Sahil Hegde, Macy J. Huston, Charles J. Law, Roel R. Lefever, Ishan Mishra, Mark Popinchalk, Sabina Sagynbayeva, Wei Yan, Kaitlin L. Ingraham Dixie, K. Supriya

Abstract: Undergraduate physics and astronomy students are expected to engage with scientific literature as they begin their research careers, but reading comprehension skills are rarely explicitly taught in major courses. We seek to determine the efficacy of lesson plans designed to improve undergraduate astronomy (or related) majors' perceived ability to engage with research literature by using accessible… ▽ More Undergraduate physics and astronomy students are expected to engage with scientific literature as they begin their research careers, but reading comprehension skills are rarely explicitly taught in major courses. We seek to determine the efficacy of lesson plans designed to improve undergraduate astronomy (or related) majors' perceived ability to engage with research literature by using accessible summaries of current research written by experts in the field. During the 2022-2023 academic year, twelve faculty members incorporated lesson plans using accessible summaries from Astrobites into their undergraduate astronomy major courses, surveyed their students before and after the activities, and participated in follow-up interviews with our research team. Quantitative and qualitative survey data clearly show that students' perceptions of their abilities with jargon, identifying main takeaways of a paper, conceptual understanding of physics and astronomy, and communicating scientific results all improved with use of the tested lesson plans. Additionally, students show evidence of increased confidence of their abilities within astronomy after exposure to these lessons, and instructors valued a ready-to-use resource to incorporate reading comprehension in their pedagogy. This case study with Astrobites-based lesson plans suggests that incorporating current research in the undergraduate classroom through accessible literature summaries may increase students' confidence and ability to engage with research literature, as well as their preparation for participation in research and applied careers. △ Less

Submitted 19 January, 2024; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: Submitted to PRPER

A severe challenge to the MOND phenomenology in our Galaxy

Authors: Man Ho Chan, Ka Chung Law

Abstract: Modified Newtonian Dynamics (MOND) is one of the most popular alternative theories of dark matter to explain the missing mass problem in galaxies. Although it remains controversial regarding MOND as a fundamental theory, MOND phenomenology has been shown to widely apply in different galaxies, which gives challenges to the standard $Λ$ cold dark matter model. In this article, we derive analytically… ▽ More Modified Newtonian Dynamics (MOND) is one of the most popular alternative theories of dark matter to explain the missing mass problem in galaxies. Although it remains controversial regarding MOND as a fundamental theory, MOND phenomenology has been shown to widely apply in different galaxies, which gives challenges to the standard $Λ$ cold dark matter model. In this article, we derive analytically the galactic rotation curve gradient in the MOND framework and present a rigorous analysis to examine the MOND phenomenology in our Galaxy. By assuming a benchmark baryonic disk density profile and two popular families of MOND interpolating functions, we show for the first time that the recent discovery of the declining Galactic rotation curve in the outer region ($R \approx 17-23$ kpc) can almost rule out the MOND phenomenology at more than $5σ$. This strongly supports some of the previous studies claiming that MOND is neither a fundamental theory nor a universal description of galactic properties. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: Accepted in ApJ

Journal ref: ApJ 957, 24 (2023)

MAPS: Constraining Serendipitous Time Variability in Protoplanetary Disk Molecular Ion Emission

Authors: Abygail R. Waggoner, L. Ilsedore Cleeves, Ryan A. Loomis, Yuri Aikawa, Jaehan Bae, Jennifer B. Bergner, Alice S. Booth, Jenny K. Calahan, Gianni Cataldi, Charles J. Law, Romane Le Gal, Feng Long, Karin I. Öberg, Richard Teague, David J. Wilner

Abstract: Theoretical models and observations suggest that the abundances of molecular ions in protoplanetary disks should be highly sensitive to the variable ionization conditions set by the young central star. We present a search for temporal flux variability of HCO+ J=1-0, which was observed as a part of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We split out and imaged t… ▽ More Theoretical models and observations suggest that the abundances of molecular ions in protoplanetary disks should be highly sensitive to the variable ionization conditions set by the young central star. We present a search for temporal flux variability of HCO+ J=1-0, which was observed as a part of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We split out and imaged the line and continuum data for each individual day the five sources were observed (HD 163296, AS 209, GM Aur, MWC 480, and IM Lup, with between 3 to 6 unique visits per source). Significant enhancement (>3σ) was not observed, but we find variations in the spectral profiles in all five disks. Variations in AS 209, GM Aur, and HD 163296 are tentatively attributed to variations in HCO+ flux, while variations in IM Lup and MWC 480 are most likely introduced by differences in the \textit{uv} coverage, which impact the amount of recovered flux during imaging. The tentative detections and low degree of variability are consistent with expectations of X-ray flare driven HCO+ variability, which requires relatively large flares to enhance the HCO+ rotational emission at significant (>20%) levels. These findings also demonstrate the need for dedicated monitoring campaigns with high signal to noise ratios to fully characterize X-ray flare driven chemistry. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: Accepted for publication in ApJ, 18 pages, 9 figures

Tracing snowlines and C/O ratio in a planet-hosting disk: ALMA molecular line observations towards the HD169142 disk

Authors: Alice S. Booth, Charles J. Law, Milou Temmink, Margot Leemker, Enrique Macias

Abstract: The composition of a forming planet is set by the material it accretes from its parent protoplanetary disk. Therefore, it is crucial to map the chemical make-up of the gas in disks to understand the chemical environment of planet formation. This paper presents molecular line observations taken with the Atacama Large Millimeter/submillimeter Array of the planet-hosting disk around the young star HD… ▽ More The composition of a forming planet is set by the material it accretes from its parent protoplanetary disk. Therefore, it is crucial to map the chemical make-up of the gas in disks to understand the chemical environment of planet formation. This paper presents molecular line observations taken with the Atacama Large Millimeter/submillimeter Array of the planet-hosting disk around the young star HD 169142. We detect N2H+, CH3OH, [CI], DCN, CS, C34S, 13CS, H2CS, H2CO, HC3N and c-C3H2 in this system for the first time. Combining these data with the recent detection of SO and previously published DCO+ data, we estimate the location of H2O and CO snowlines and investigate radial variations in the gas phase C/O ratio. We find that the HD 169142 disk has a relatively low N2H+ flux compared to the disks around Herbig stars HD 163296 and MWC 480 indicating less CO freeze-out and place the CO snowline beyond the millimetre disk at ~150 au. The detection of CH3OH from the inner disk is consistent with the H2O snowline being located at the edge of the central dust cavity at ~20 au. The radially varying CS/SO ratio across the proposed H2O snowline location is consistent with this interpretation. Additionally, the detection of CH3OH in such a warm disk adds to the growing evidence supporting the inheritance of complex ices in disks from the earlier, colder stages of star formation. Finally, we propose that the giant HD 169142 b located at 37 au is forming between the CO2 and H2O snowlines where the local elemental make of the gas is expected to have C/O=1.0. △ Less

Submitted 15 August, 2023; originally announced August 2023.

Comments: Accepted A&A 13th August 2023

Journal ref: A&A 678, A146 (2023)

Deep Synoptic Array Science: Implications of Faraday Rotation Measures of Localized Fast Radio Bursts

Authors: Myles B. Sherman, Liam Connor, Vikram Ravi, Casey Law, Ge Chen, Kritti Sharma, Morgan Catha, Jakob T. Faber, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, James W. Lamb, Paul Rasmussen, Jun Shi, Dana Simard, Jean Somalwar, Reynier Squillace, Sander Weinreb, David P. Woody, Nitika Yadlapalli

Abstract: Faraday rotation measures (RMs) of fast radio bursts (FRBs) offer the prospect of directly measuring extragalactic magnetic fields. We present an analysis of the RMs of ten as yet non-repeating FRBs detected and localized to host galaxies by the 110-antenna Deep Synoptic Array (DSA-110). We combine this sample with published RMs of 15 localized FRBs, nine of which are repeating sources. For each F… ▽ More Faraday rotation measures (RMs) of fast radio bursts (FRBs) offer the prospect of directly measuring extragalactic magnetic fields. We present an analysis of the RMs of ten as yet non-repeating FRBs detected and localized to host galaxies by the 110-antenna Deep Synoptic Array (DSA-110). We combine this sample with published RMs of 15 localized FRBs, nine of which are repeating sources. For each FRB in the combined sample, we estimate the host-galaxy dispersion measure (DM) contributions and extragalactic RM. We find compelling evidence that the extragalactic components of FRB RMs are often dominated by contributions from the host-galaxy interstellar medium (ISM). Specifically, we find that both repeating and as yet non-repeating FRBs show a correlation between the host-DM and host-RM in the rest frame, and we find an anti-correlation between extragalactic RM (in the observer frame) and redshift for non-repeaters, as expected if the magnetized plasma is in the host galaxy. Important exceptions to the ISM origin include a dense, magnetized circum-burst medium in some repeating FRBs, and the intra-cluster medium (ICM) of host or intervening galaxy clusters. We find that the estimated ISM magnetic-field strengths, $\bar{B}_{||}$, are characteristically larger than those inferred from Galactic radio pulsars. This suggests either increased ISM magnetization in FRB hosts in comparison with the Milky Way, or that FRBs preferentially reside in regions of increased magnetic-field strength within their hosts. △ Less

Submitted 7 September, 2023; v1 submitted 13 August, 2023; originally announced August 2023.

Comments: 13 pages, 4 figures, 1 table ; added references to section 1, 2, and 3; added a paragraph in Section 3 to discuss simulations addressing interplay of local and ISM magnetic field; added bar to B|| to indicate average over line of sight

Deep Synoptic Array Science: Polarimetry of 25 New Fast Radio Bursts Provides Insights into their Origins

Authors: Myles B. Sherman, Liam Connor, Vikram Ravi, Casey Law, Ge Chen, Morgan Catha, Jakob T. Faber, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, James W. Lamb, Paul Rasmussen, Kritti Sharma, Jun Shi, Dana Simard, Jean Somalwar, Reynier Squillace, Sander Weinreb, David P. Woody, Nitika Yadlapalli

Abstract: We report on a full-polarization analysis of the first 25 as yet non-repeating FRBs detected at 1.4 GHz by the 110-antenna Deep Synoptic Array (DSA-110) during commissioning observations. We present details of the data-reduction, calibration, and analysis procedures developed for this novel instrument. Faraday rotation measures (RMs) are searched between $\pm10^6$ rad m$^{-2}$ and detected for 20… ▽ More We report on a full-polarization analysis of the first 25 as yet non-repeating FRBs detected at 1.4 GHz by the 110-antenna Deep Synoptic Array (DSA-110) during commissioning observations. We present details of the data-reduction, calibration, and analysis procedures developed for this novel instrument. Faraday rotation measures (RMs) are searched between $\pm10^6$ rad m$^{-2}$ and detected for 20 FRBs with magnitudes ranging from $4-4670$ rad m$^{-2}$. $15/25$ FRBs are consistent with 100% polarization, 10 of which have high ($\ge70\%$) linear-polarization fractions and 2 of which have high ($\ge30\%$) circular-polarization fractions. Our results disfavor multipath RM scattering as a dominant depolarization mechanism. Polarization-state and possible RM variations are observed in the four FRBs with multiple sub-components. We combine the DSA-110 sample with polarimetry of previously published FRBs, and compare the polarization properties of FRB sub-populations and FRBs with Galactic pulsars. Although FRB polarization fractions are typically higher than those of Galactic pulsars, and cover a wider range than those of pulsar single pulses, they resemble those of the youngest (characteristic ages $<10^{5}$ yr) pulsars. Our results support a scenario wherein FRB emission is intrinsically highly linearly polarized, and propagation effects can result in conversion to circular polarization and depolarization. Young pulsar emission and magnetospheric-propagation geometries may form a useful analogy for the origin of FRB polarization. △ Less

Submitted 15 February, 2024; v1 submitted 13 August, 2023; originally announced August 2023.

Comments: 57 pages, 24 figures, 5 tables ; version accepted by ApJ ; revised classification to include polarized signal-to-noise; added appendices for PA error/variability, J1935 polarization ; corrected literature FRB sample and supplemented with omitted data; added discussion of biases, limitations, and selection criteria; revised sigmaRM analysis ; additional plots added ; revised FRB20220424E analysis

Temporal and Spectral Properties of the Persistent Radio Source Associated with FRB 20190520B with the VLA

Authors: Xian Zhang, Wenfei Yu, Casey Law, Di Li, Shami Chatterjee, Paul Demorest, Zhen Yan, Chenhui Niu, Kshitij Aggarwal, Reshma Anna-Thomas, Sarah Burke-Spolaor, Liam Connor, Chao-wei Tsai, Weiwei Zhu, Gan Luo

Abstract: Among more than 800 known fast radio bursts (FRBs), only two, namely FRB 20121102A and FRB 20190520B, are confirmed to be associated with a persistent radio sources (PRS). Here we report evidence of apparent temporal variability in the PRS associated with the bursting FRB 20190520B based on the Karl G. Jansky Very Large Array (VLA) observations taken in 2020 and 2021. Based on the analysis of epoc… ▽ More Among more than 800 known fast radio bursts (FRBs), only two, namely FRB 20121102A and FRB 20190520B, are confirmed to be associated with a persistent radio sources (PRS). Here we report evidence of apparent temporal variability in the PRS associated with the bursting FRB 20190520B based on the Karl G. Jansky Very Large Array (VLA) observations taken in 2020 and 2021. Based on the analysis of epoch-to-epoch variability of the PRS at L, S, C, and X band in 1-12 GHz, we detected not only overall marginal variability but also a likely radio flux decrease ($\sim$ 3.2 $σ$) between the observations taken in 2020 and 2021 at 3 GHz. Assuming no spectral variation in the PRS during these observations, we found the evidence for an overall broadband radio flux decrease by about 20 percent between the 2020 and the 2021 observations, suggesting that the PRS probably evolves on the yearly time scale. If we attribute the marginal variability at 3 GHz as intrinsic or due to scintillation, the size of potential variable component of the PRS is constrained to be sub-parsec. On the other hand, the size of the PRS can be also constrained to be larger than about 0.22 parsec from the averaged radio spectrum and the integrated radio luminosity in the 1-12 GHz band based on equipartition and self-absorption arguments. We discuss potential origins of the PRS and suggest that an accreting compact object origin might be able to explain the PRS's temporal and spectral properties. Confirmation of variability or flux decline of the PRS would be critical to our understanding of the PRS and its relation to the bursting source. △ Less

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

Comments: 12 pages, 3 figures, accepted for publication in ApJ

Deep Synoptic Array Science: First FRB and Host Galaxy Catalog

Authors: C. J. Law, K. Sharma, V. Ravi, G. Chen, M. Catha, L. Connor, J. T. Faber, G. Hallinan, C. Harnach, G. Hellbourg, R. Hobbs, D. Hodge, M. Hodges, J. W. Lamb, P. Rasmussen, M. B. Sherman, J. Shi, D. Simard, R. Squillace, S. Weinreb, D. P. Woody, N. Yadlapalli

Abstract: Fast Radio Bursts (FRBs) are a powerful and mysterious new class of transient that are luminous enough to be detected at cosmological distances. By associating FRBs to host galaxies, we can measure intrinsic and environmental properties that test FRB origin models, in addition to using them as precise probes of distant cosmic gas. The Deep Synoptic Array (DSA-110) is a radio interferometer built t… ▽ More Fast Radio Bursts (FRBs) are a powerful and mysterious new class of transient that are luminous enough to be detected at cosmological distances. By associating FRBs to host galaxies, we can measure intrinsic and environmental properties that test FRB origin models, in addition to using them as precise probes of distant cosmic gas. The Deep Synoptic Array (DSA-110) is a radio interferometer built to maximize the rate at which it can simultaneously detect and localize FRBs. Here, we present the first sample of FRBs and host galaxies discovered by the DSA-110. This sample of 11 FRBs is the largest, most uniform sample of localized FRBs to date, as it is selected based on association to host galaxies identified in optical imaging by Pan-STARRS1. These FRBs have not been observed to repeat and their radio properties (dispersion, temporal scattering, energy) are similar to that of the known non-repeating FRB population. Most host galaxies have ongoing star formation, as has been identified before for FRB hosts. Two hosts of the new sample are massive, quiescent galaxies. The distribution of star-formation history across this host-galaxy sample shows that the delay-time distribution is wide, with a powerlaw model that spans from $\sim100$\,Myr to $\gtrsim2$\,Gyr. This requires the existence of one or more progenitor formation channels associated with old stellar populations, such as the binary evolution of compact objects. △ Less

Submitted 24 January, 2024; v1 submitted 6 July, 2023; originally announced July 2023.

Comments: 21 pages. Submitted to AAS Journals. Includes changes based on referee comments and improved host galaxy analysis

Constraining the gas distribution in the PDS 70 disk as a method to assess the effect of planet-disk interactions

Authors: B. Portilla-Revelo, I. Kamp, S. Facchini, E. F. van Dishoeck, C. Law, Ch. Rab, J. Bae, M. Benisty, K. Öberg, R. Teague

Abstract: Embedded planets are potentially the cause of substructures like gaps and cavities observed in several protoplanetary disks. Thus, the substructures observed in the continuum and in line emission encode information about the presence of planets in the system and how they interact with the natal disk. The pre-transitional disk around the star PDS 70 is the first case of two young planets imaged wit… ▽ More Embedded planets are potentially the cause of substructures like gaps and cavities observed in several protoplanetary disks. Thus, the substructures observed in the continuum and in line emission encode information about the presence of planets in the system and how they interact with the natal disk. The pre-transitional disk around the star PDS 70 is the first case of two young planets imaged within a dust depleted gap that was likely carved by themselves. We aim to determine the spatial distribution of the gas and dust components in the PDS 70 disk. The axisymmetric substructures observed in the resulting profiles are interpreted in the context of planet-disk interactions. We develop a thermo-chemical forward model for an axisymmetric disk to explain a subset of the Atacama Large Millimeter/Submillimeter Array (ALMA) band 6 observations of three CO isotopologues plus the continuum towards PDS 70. Combining the inferred gas and dust distributions, the model results in a variable gas-to-dust ratio profile throughout the disk that spans two orders of magnitude within the first $130$ au and shows a step gradient towards the outer disk, which is consistent with the presence of a pressure maxima driven by planet-disk interactions. We find a gas density drop factor of ${\sim} 19$ at the location of the planet PDS 70 c with respect to the peak gas density at $75$ au. Combining this value with literature results on the hydrodynamics of planet-disk interactions, we find this gas gap depth to be consistent with independent planet mass estimates from infrared observations. Our findings point towards gas stirring processes taking place in the common gap due to the gravitational perturbation of both planets. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 677, A76 (2023)

SO and SiS Emission Tracing an Embedded Planet and Compact $^{12}$CO and $^{13}$CO Counterparts in the HD 169142 Disk

Authors: Charles J. Law, Alice S. Booth, Karin I. Öberg

Abstract: Planets form in dusty, gas-rich disks around young stars, while at the same time, the planet formation process alters the physical and chemical structure of the disk itself. Embedded planets will locally heat the disk and sublimate volatile-rich ices, or in extreme cases, result in shocks that sputter heavy atoms such as Si from dust grains. This should cause chemical asymmetries detectable in mol… ▽ More Planets form in dusty, gas-rich disks around young stars, while at the same time, the planet formation process alters the physical and chemical structure of the disk itself. Embedded planets will locally heat the disk and sublimate volatile-rich ices, or in extreme cases, result in shocks that sputter heavy atoms such as Si from dust grains. This should cause chemical asymmetries detectable in molecular gas observations. Using high-angular-resolution ALMA archival data of the HD 169142 disk, we identify compact SO J=8$_8$-7$_7$ and SiS J=19-18 emission coincident with the position of a ${\sim}$2 M$_{\rm{Jup}}$ planet seen as a localized, Keplerian NIR feature within a gas-depleted, annular dust gap at ${\approx}$38 au. The SiS emission is located along an azimuthal arc and has a similar morphology as a known $^{12}$CO kinematic excess. This is the first tentative detection of SiS emission in a protoplanetary disk and suggests that the planet is driving sufficiently strong shocks to produce gas-phase SiS. We also report the discovery of compact $^{12}$CO and $^{13}$CO J=3-2 emission coincident with the planet location. Taken together, a planet-driven outflow provides the best explanation for the properties of the observed chemical asymmetries. We also resolve a bright, azimuthally-asymmetric SO ring at ${\approx}$24 au. While most of this SO emission originates from ice sublimation, its asymmetric distribution implies azimuthal temperature variations driven by a misaligned inner disk or planet-disk interactions. Overall, the HD 169142 disk shows several distinct chemical signatures related to giant planet formation and presents a powerful template for future searches of planet-related chemical asymmetries in protoplanetary disks. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Comments: 22 pages, 12 figures, accepted for publication in ApJL

Detecting Fast Radio Bursts with Spectral Structure using the Continuous Forward Algorithm

Authors: Pravir Kumar, Barak Zackay, Casey J. Law

Abstract: Detecting fast radio bursts (FRBs) with frequency-dependent intensity remains a challenge, as existing search algorithms do not account for the spectral shape, potentially leading to non-detections. We propose a novel detection statistic, which we call the Kalman detector, that improves the sensitivity of FRB signal detection by incorporating spectral shape information. The detection statistic is… ▽ More Detecting fast radio bursts (FRBs) with frequency-dependent intensity remains a challenge, as existing search algorithms do not account for the spectral shape, potentially leading to non-detections. We propose a novel detection statistic, which we call the Kalman detector, that improves the sensitivity of FRB signal detection by incorporating spectral shape information. The detection statistic is based on an optimal matched filter, marginalizing over all possible intensity functions, weighted by a random walk probability distribution, considering some decorrelation bandwidth. Our analysis of previously detected FRBs demonstrates that the Kalman score provides a comparable yet independent source of information for bursts with significant spectral structure, and the sensitivity improvement is of the order 0%--200% with a median improvement of 20%. We also applied the Kalman detector to existing data from FRB 20201124A and detected two new repeat bursts that were previously missed. Furthermore, we suggest a practical implementation for real-time surveys by employing a low significance soft-trigger from initial flux integration-based detection algorithms. The Kalman detector has the potential to significantly enhance FRB detection capabilities and enable new insights into the spectral properties of these enigmatic astrophysical phenomena. △ Less

Submitted 8 December, 2023; v1 submitted 13 June, 2023; originally announced June 2023.

Comments: 16 pages, 8 figures; accepted for publication in ApJ

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)

The JCMT BISTRO Survey: Studying the Complex Magnetic Field of L43

Authors: Janik Karoly, Derek Ward-Thompson, Kate Pattle, David Berry, Anthony Whitworth, Jason Kirk, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar… ▽ More We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to $\sim$160$\pm$30 $μ$G in the main starless core and up to $\sim$90$\pm$40 $μ$G in the more diffuse, extended region. These field strengths give magnetically super- and sub-critical values respectively and both are found to be roughly trans-Alfvénic. We also present a new method of data reduction for these denser but fainter objects like starless cores. △ Less

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

Comments: Accepted for publication in ApJ. 23 pages, 9 figures (7 main text, 2 appendix)

The edge-on protoplanetary disk HH 48 NE II. Modeling ices and silicates

Authors: J. A. Sturm, M. K. McClure, J. B. Bergner, D. Harsono, E. Dartois, M. N. Drozdovskaya, S. Ioppolo, K. I. Öberg, C. J. Law, M. E. Palumbo, Y. J. Pendleton, W. R. M. Rocha, H. Terada, R. G. Urso

Abstract: The abundance and distribution of ice in protoplanetary disks (PPD) is critical to understand the linkage between the composition of circumstellar matter and the composition of exoplanets. Edge-on PPDs are a useful tool to constrain such ice composition and its location in the disk, as ice spectral signatures can be observed in absorption against the continuum emission arising from the warmer cent… ▽ More The abundance and distribution of ice in protoplanetary disks (PPD) is critical to understand the linkage between the composition of circumstellar matter and the composition of exoplanets. Edge-on PPDs are a useful tool to constrain such ice composition and its location in the disk, as ice spectral signatures can be observed in absorption against the continuum emission arising from the warmer central disk regions. The aim of this work is to model ice absorption features in PPDs and determine how well the abundance of the main ice species across the disk can be determined within the uncertainty of the physical parameter space. The edge-on PPD around HH 48 NE, a target of the JWST ERS program IceAge, is used as a reference system. We use RADMC-3D to raytrace the mid-infrared continuum. Using a constant parameterized ice abundance, ice opacities are added to the dust opacity in regions wherever the disk is cold enough for the main carbon, oxygen and nitrogen carriers to freeze out. The global abundance of the main ice carriers in HH 48 NE can be determined within a factor of 3, when taking the uncertainty of the physical parameters into account. Ice features in PPDs can be saturated at an optical depth <1, due to local saturation. Spatially observed ice optical depths cannot be directly related to column densities due to radiative transfer effects. Vertical snowlines will not be a clear transition due to the radially increasing height of the snowsurface, but their location may be constrained from observations using radiative transfer modeling. Radial snowlines are not really accesible. Not only the ice abundance, but also inclination, settling, grain size distribution and disk mass have strong impact on the observed ice absorption features in disks. Relative changes in ice abundance can be inferred from observations only if the source structure is well constrained △ Less

Submitted 8 May, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 18 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

The edge-on protoplanetary disk HH 48 NE I. Modeling the geometry and stellar parameters

Authors: J. A. Sturm, M. K. McClure, C. J. Law, D. Harsono, J. B. Bergner, E. Dartois, M. N. Drozdovskaya, S. Ioppolo, K. I. Öberg, M. E. Palumbo, Y. J. Pendleton, W. R. M. Rocha, H. Terada, R. G. Urso

Abstract: Context. Observations of edge-on disks are an important tool for constraining general protoplanetary disk properties that cannot be determined in any other way. However, most radiative transfer models cannot simultaneously reproduce the spectral energy distributions (SEDs) and resolved scattered light and submillimeter observations of these systems, due to the differences in geometry and dust prop… ▽ More Context. Observations of edge-on disks are an important tool for constraining general protoplanetary disk properties that cannot be determined in any other way. However, most radiative transfer models cannot simultaneously reproduce the spectral energy distributions (SEDs) and resolved scattered light and submillimeter observations of these systems, due to the differences in geometry and dust properties at different wavelengths. Aims. We simultaneously constrain the geometry of the edge-on protoplanetary disk HH 48 NE and the characteristics of the host star. HH 48 NE is part of the JWST early release science program Ice Age. This work serves as a stepping stone towards a better understanding of the disk physical structure and icy chemistry in this particular source. This kind of modeling lays the groundwork for studying other edge-on sources to be observed with the JWST. Methods. We fit a parameterized dust model to HH 48 NE by coupling the radiative transfer code RADMC-3D and an MCMC framework. The dust structure was fitted independently to a compiled SED, a scattered light image at 0.8 $μ$m and an ALMA dust continuum observation at 890 $μ$m. Results. We find that 90% of the dust mass in HH 48 NE is settled to the disk midplane, less than in average disks, and that the atmospheric layers of the disk contain exclusively large grains (0.3-10 $μ$m). The exclusion of small grains in the upper atmosphere likely has important consequences for the chemistry due to the deep penetration of high-energy photons. The addition of a relatively large cavity (ca. 50 au in radius) is necessary to explain the strong mid-infrared emission, and to fit the scattered light and continuum observations simultaneously. △ Less

Submitted 8 May, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 16 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 677, A17 (2023)

Molecules with ALMA at Planet-forming Scales (MAPS). Complex Kinematics in the AS 209 Disk Induced by a Forming Planet and Disk Winds

Authors: Maria Galloway-Sprietsma, Jaehan Bae, Richard Teague, Myriam Benisty, Stefano Facchini, Yuri Aikawa, Felipe Alarcón, Sean M. Andrews, Edwin Bergin, Gianni Cataldi, L. Ilsedore Cleeves, Ian Czekala, Viviana V. Guzmán, Jane Huang, Charles J. Law, Romane Le Gal, Yao Liu, Feng Long, François Ménard, Karin I. Öberg, Catherine Walsh, David J. Wilner

Abstract: We study the kinematics of the AS 209 disk using the J=2-1 transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O. We derive the radial, azimuthal, and vertical velocity of the gas, taking into account the lowered emission surface near the annular gap at ~1.7 (200 au) within which a candidate circumplanetary disk-hosting planet has been reported previously. In $^{12}$CO and $^{13}$CO, we find a cohere… ▽ More We study the kinematics of the AS 209 disk using the J=2-1 transitions of $^{12}$CO, $^{13}$CO, and C$^{18}$O. We derive the radial, azimuthal, and vertical velocity of the gas, taking into account the lowered emission surface near the annular gap at ~1.7 (200 au) within which a candidate circumplanetary disk-hosting planet has been reported previously. In $^{12}$CO and $^{13}$CO, we find a coherent upward flow arising from the gap. The upward gas flow is as fast as $150~{\rm m~s}^{-1}$ in the regions traced by $^{12}$CO emission, which corresponds to about 50% of the local sound speed or $6\%$ of the local Keplerian speed. Such an upward gas flow is difficult to reconcile with an embedded planet alone. Instead, we propose that magnetically driven winds via ambipolar diffusion are triggered by the low gas density within the planet-carved gap, dominating the kinematics of the gap region. We estimate the ambipolar Elsasser number, Am, using the HCO$^+$ column density as a proxy for ion density and find that Am is ~0.1 at the radial location of the upward flow. This value is broadly consistent with the value at which numerical simulations find ambipolar diffusion drives strong winds. We hypothesize the activation of magnetically-driven winds in a planet-carved gap can control the growth of the embedded planet. We provide a scaling relationship which describes the wind-regulated terminal mass: adopting parameters relevant to 100 au from a solar-mass star, we find the wind-regulated terminal mass is about one Jupiter mass, which may help explain the dearth of directly imaged super-Jovian-mass planets. △ Less

Submitted 12 May, 2023; v1 submitted 7 April, 2023; originally announced April 2023.

Comments: This paper has been accepted for publication in the Astrophysical Journal (ApJ)

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

An SMA Survey of Chemistry in Disks around Herbig AeBe Stars

Authors: Jamila Pegues, Karin I. Öberg, Chunhua Qi, Sean M. Andrews, Jane Huang, Charles J. Law, Romane Le Gal, Luca Matrà, David J. Wilner

Abstract: Protoplanetary disks around Herbig AeBe stars are exciting targets for studying the chemical environments where giant planets form. Save for a few disks, however, much of Herbig AeBe disk chemistry is an open frontier. We present a Submillimeter Array (SMA) $\sim$213-268 GHz pilot survey of mm continuum, CO isotopologues, and other small molecules in disks around five Herbig AeBe stars (HD 34282,… ▽ More Protoplanetary disks around Herbig AeBe stars are exciting targets for studying the chemical environments where giant planets form. Save for a few disks, however, much of Herbig AeBe disk chemistry is an open frontier. We present a Submillimeter Array (SMA) $\sim$213-268 GHz pilot survey of mm continuum, CO isotopologues, and other small molecules in disks around five Herbig AeBe stars (HD 34282, HD 36112, HD 38120, HD 142666, and HD 144432). We detect or tentatively detect $^{12}$CO 2--1 and $^{13}$CO 2--1 from four disks; C$^{18}$O 2--1 and HCO$^+$ 3--2 from three disks; HCN 3--2, CS 5--4, and DCO$^+$ 3--2 from two disks; and C$_2$H 3--2 and DCN 3--2 from one disk each. H$_2$CO 3--2 is undetected at the sensitivity of our observations. The mm continuum images of HD 34282 suggest a faint, unresolved source $\sim$5\farcs0 away, which could arise from a distant orbital companion or an extended spiral arm. We fold our sample into a compilation of T Tauri and Herbig AeBe/F disks from the literature. Altogether, most line fluxes generally increase with mm continuum flux. Line flux ratios between CO 2--1 isotopologues are nearest to unity for the Herbig AeBe/F disks. This may indicate emitting layers with relatively similar, warmer temperatures and more abundant CO relative to disk dust mass. Lower HCO$^+$ 3--2 flux ratios may reflect less ionization in Herbig AeBe/F disks. Smaller detection rates and flux ratios for DCO$^+$ 3--2, DCN 3--2, and H$_2$CO 3--2 suggest smaller regimes of cold chemistry around the luminous Herbig AeBe/F stars. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Comments: 43 pages (23 pages in the main document, 20 pages in the appendix), 16 figures, 16 tables. Accepted in ApJ (February 2023)

Deep Synoptic Array science: Two fast radio burst sources in massive galaxy clusters

Authors: Liam Connor, Vikram Ravi, Morgan Catha, Ge Chen, Jakob T. Faber, James W. Lamb, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, Casey Law, Paul Rasmussen, Jack Sayers, Kritti Sharma, Myles B. Sherman, Jun Shi, Dana Simard, Jean Somalwar, Reynier Squillace, Sander Weinreb, David P. Woody, Nitika Yadlapalli

Abstract: The hot gas that constitutes the intracluster medium (ICM) has been studied at X-ray and millimeter/sub-millimeter wavelengths (Sunyaev-Zeldovich effect) for decades. Fast radio bursts (FRBs) offer an additional method of directly measuring the ICM and gas surrounding clusters, via observables such as dispersion measure (DM) and Faraday rotation measure (RM). We report the discovery of two FRB sou… ▽ More The hot gas that constitutes the intracluster medium (ICM) has been studied at X-ray and millimeter/sub-millimeter wavelengths (Sunyaev-Zeldovich effect) for decades. Fast radio bursts (FRBs) offer an additional method of directly measuring the ICM and gas surrounding clusters, via observables such as dispersion measure (DM) and Faraday rotation measure (RM). We report the discovery of two FRB sources detected with the Deep Synoptic Array (DSA-110) whose host galaxies belong to massive galaxy clusters. In both cases, the FRBs exhibit excess extragalactic DM, some of which likely originates in the ICM of their respective clusters. FRB 20220914A resides in the galaxy cluster Abell 2310 at z=0.1125 with a projected offset from the cluster center of 520 kpc. The host of a second source, FRB 20220509G, is an elliptical galaxy at z=0.0894 that belongs to the galaxy cluster Abell 2311 at projected offset of 870 kpc. These sources represent the first time an FRB has been localized to a galaxy cluster. We combine our FRB data with archival X-ray, SZ, and optical observations of these clusters in order to infer properties of the ICM, including a measurement of gas temperature from DM and ySZ of 0.8-3.9 keV. We then compare our results to massive cluster halos from the IllustrisTNG simulation. Finally, we describe how large samples of localized FRBs from future surveys will constrain the ICM, particularly beyond the virial radius of clusters. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Deep Synoptic Array science: A massive elliptical host among two galaxy-cluster fast radio bursts

Authors: Kritti Sharma, Jean Somalwar, Casey Law, Vikram Ravi, Morgan Catha, Ge Chen, Liam Connor, Jakob T. Faber, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, James W. Lamb, Paul Rasmussen, Myles B. Sherman, Jun Shi, Dana Simard, Reynier Squillace, Sander Weinreb, David P. Woody, Nitika Yadlapalli

Abstract: The stellar population environments associated with fast radio burst (FRB) sources provide important insights for developing their progenitor theories. We expand the diversity of known FRB host environments by reporting two FRBs in massive galaxy clusters discovered by the Deep Synoptic Array (DSA-110) during its commissioning observations. FRB 20220914A has been localized to a star-forming, late-… ▽ More The stellar population environments associated with fast radio burst (FRB) sources provide important insights for developing their progenitor theories. We expand the diversity of known FRB host environments by reporting two FRBs in massive galaxy clusters discovered by the Deep Synoptic Array (DSA-110) during its commissioning observations. FRB 20220914A has been localized to a star-forming, late-type galaxy at a redshift of 0.1139 with multiple starbursts at lookback times less than $\sim$3.5 Gyr in the Abell 2310 galaxy cluster. Although the host galaxy of FRB 20220914A is similar to typical FRB hosts, the FRB 20220509G host stands out as a quiescent, early-type galaxy at a redshift of 0.0894 in the Abell 2311 galaxy cluster. The discovery of FRBs in both late and early-type galaxies adds to the body of evidence that the FRB sources have multiple formation channels. Therefore, even though FRB hosts are typically star-forming, there must exist formation channels consistent with old stellar population in galaxies. The varied star formation histories of the two FRB hosts we report indicate a wide delay-time distribution of FRB progenitors. Future work in constraining the FRB delay-time distribution, using methods we develop herein, will prove crucial in determining the evolutionary histories of FRB sources. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: 18 pages, 9 figures

First BISTRO observations of the dark cloud Taurus L1495A-B10: the role of the magnetic field in the earliest stages of low-mass star formation

Authors: Derek Ward-Thompson, Janik Karoly, Kate Pattle, Anthony Whitworth, Jason Kirk, David Berry, Pierre Bastien, Tao-Chung Ching, Simon Coude, Jihye Hwang, Woojin Kwon, Archana Soam, Jia-Wei Wang, Tetsuo Hasegawa, Shih-Ping Lai, Keping Qiu, Doris Arzoumanian, Tyler L. Bourke, Do-Young Byun, Huei-Ru Vivien Chen, Wen Ping Chen, Mike Chen, Zhiwei Chen, Jungyeon Cho, Minho Choi , et al. (133 additional authors not shown)

Abstract: We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to… ▽ More We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face-on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field-dominated to being matter-dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 years (Mestel 1965). △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 14 pages, 5 figures. ApJ accepted

HST Proper Motion Measurements of Supernova Remnant N132D: Center of Expansion and Age

Authors: John Banovetz, Dan Milisavljevic, Niharika Sravan, Kathryn E. Weil, Bhagya Subrayan, Robert A. Fesen, Daniel J. Patnaude, Paul P. Plucinsky, Charles J. Law, William P. Blair, Jon A. Morse

Abstract: We present proper motion measurements of oxygen-rich ejecta of the LMC supernova remnant N132D using two epochs of Hubble Space Telescope Advanced Camera for Surveys data spanning 16 years. The proper motions of 120 individual knots of oxygen-rich gas were measured and used to calculate a center of expansion (CoE) of $α$=05:25:01.71 and $δ$=-69:38:41.64 (J2000) with a 1-$σ$ uncertainty of 2.90 arc… ▽ More We present proper motion measurements of oxygen-rich ejecta of the LMC supernova remnant N132D using two epochs of Hubble Space Telescope Advanced Camera for Surveys data spanning 16 years. The proper motions of 120 individual knots of oxygen-rich gas were measured and used to calculate a center of expansion (CoE) of $α$=05:25:01.71 and $δ$=-69:38:41.64 (J2000) with a 1-$σ$ uncertainty of 2.90 arcseconds. This new CoE measurement is 9.2 and 10.8 arcseconds from two previous CoE estimates based on the geometry of the optically emitting ejecta. We also derive an explosion age of 2770 $\pm$ 500 yr, which is consistent with recent age estimates of $\approx 2500$ yr made from 3D ejecta reconstructions. We verify our estimates of the CoE and age using a new automated procedure that detected and tracked the proper motions of 137 knots, with 73 knots that overlap with the visually identified knots. We find the proper motions of ejecta are still ballistic, despite the remnant's age, and are consistent with the notion that the ejecta are expanding into an ISM cavity. Evidence for explosion asymmetry from the parent supernova is also observed. Using the visually measured proper motion measurements and corresponding center of expansion and age, we compare N132D to other supernova remnants with proper motion ejecta studies. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 26 pages, 16 figures, accepted to ApJ

Deep Synoptic Array science: a 50 Mpc fast radio burst constrains the mass of the Milky Way circumgalactic medium

Authors: Vikram Ravi, Morgan Catha, Ge Chen, Liam Connor, James M. Cordes, Jakob T. Faber, James W. Lamb, Gregg Hallinan, Charlie Harnach, Greg Hellbourg, Rick Hobbs, David Hodge, Mark Hodges, Casey Law, Paul Rasmussen, Kritti Sharma, Myles B. Sherman, Jun Shi, Dana Simard, Jean J. Somalwar, Reynier Squillace, Sander Weinreb, David P. Woody, Nitika Yadlapalli

Abstract: We present the Deep Synoptic Array (DSA-110) discovery and interferometric localization of the so far non-repeating FRB 20220319D. The FRB originates in a young, rapidly star-forming barred spiral galaxy, IRAS 02044$+$7048, at just 50 Mpc. Although the NE2001 and YMW16 models for the Galactic interstellar-medium (ISM) contribution to the DM of FRB 20220319D exceed its total observed DM, we show th… ▽ More We present the Deep Synoptic Array (DSA-110) discovery and interferometric localization of the so far non-repeating FRB 20220319D. The FRB originates in a young, rapidly star-forming barred spiral galaxy, IRAS 02044$+$7048, at just 50 Mpc. Although the NE2001 and YMW16 models for the Galactic interstellar-medium (ISM) contribution to the DM of FRB 20220319D exceed its total observed DM, we show that uncertainties in these models accommodate an extragalactic origin for the burst. We derive a conservative upper limit on the DM contributed by the circumgalactic medium (CGM) of the Milky Way: the limit is either 28.7 pc cm$^{-3}$ and 47.3 pc cm$^{-3}$, depending on which of two pulsars nearby on the sky to FRB 20220319D is used to estimate the ISM DM. These limits both imply that the total Galactic CGM mass is $<10^{11}M_{\odot}$, and that the baryonic mass of the Milky Way is $\lesssim60\%$ of the cosmological average given the total halo mass. More stringent albeit less conservative constraints are possible when the DMs of pulsars in the distant globular cluster M53 are additionally considered. Although our constraints are sensitive to possible anisotropy in the CGM and to the assumed form of the radial-density profile, they are not subject to uncertainties in the chemical and thermal properties of the CGM. Our results strongly support scenarios commonly predicted by galaxy-formation simulations wherein feedback processes expel baryonic matter from the halos of galaxies like the Milky Way. △ Less

Submitted 3 January, 2023; originally announced January 2023.

Comments: 17 pages, 8 figures, 3 tables, submitted to AAS Journals

JCMT BISTRO Observations: Magnetic Field Morphology of Bubbles Associated with NGC 6334

Authors: Mehrnoosh Tahani, Pierre Bastien, Ray S. Furuya, Kate Pattle, Doug Johnstone, Doris Arzoumanian, Yasuo Doi, Tetsuo Hasegawa, Shu-ichiro Inutsuka, Simon Coudé, Laura Fissel, Michael Chun-Yuan Chen, Frédérick Poidevin, Sarah Sadavoy, Rachel Friesen, Patrick M. Koch, James Di Francesco, Gerald H. Moriarty-Schieven, Zhiwei Chen, Eun Jung Chung, Chakali Eswaraiah, Lapo Fanciullo, Tim Gledhill, Valentin J. M. Le Gouellec, Thiem Hoang , et al. (120 additional authors not shown)

Abstract: We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles… ▽ More We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles indicate that the gas and magnetic field lines have been pushed away from the bubble, toward an almost tangential (to the bubble) magnetic field morphology. In the densest part of NGC 6334, where the magnetic field morphology is similar to an hourglass, the polarization observations do not exhibit observable impact from HII regions. We detect two nested radial polarization patterns in a bubble to the south of NGC 6334 that correspond to the previously observed bipolar structure in this bubble. Finally, using the results of this study, we present steps (incorporating computer vision; circular Hough Transform) that can be used in future studies to identify bubbles that have physically impacted magnetic field lines. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astrophysical Journal (ApJ)

Mapping Protoplanetary Disk Vertical Structure with CO Isotopologue Line Emission

Authors: Charles J. Law, Richard Teague, Karin I. Öberg, Evan A. Rich, Sean M. Andrews, Jaehan Bae, Myriam Benisty, Stefano Facchini, Kevin Flaherty, Andrea Isella, Sheng Jin, Jun Hashimoto, Jane Huang, Ryan A. Loomis, Feng Long, Carlos E. Muñoz-Romero, Teresa Paneque-Carreño, Laura M. Pérez, Chunhua Qi, Kamber R. Schwarz, Jochen Stadler, Takashi Tsukagoshi, David J. Wilner, Gerrit van der Plas

Abstract: High spatial resolution observations of CO isotopologue line emission in protoplanetary disks at mid-inclinations (${\approx}$30-75°) allow us to characterize the gas structure in detail, including radial and vertical substructures, emission surface heights and their dependencies on source characteristics, and disk temperature profiles. By combining observations of a suite of CO isotopologues, we… ▽ More High spatial resolution observations of CO isotopologue line emission in protoplanetary disks at mid-inclinations (${\approx}$30-75°) allow us to characterize the gas structure in detail, including radial and vertical substructures, emission surface heights and their dependencies on source characteristics, and disk temperature profiles. By combining observations of a suite of CO isotopologues, we can map the 2D (r, z) disk structure from the disk upper atmosphere, as traced by CO, to near the midplane, as probed by less abundant isotopologues. Here, we present high angular resolution (${\lesssim}$0."1 to ${\approx}$0."2; ${\approx}$15-30 au) observations of CO, $^{13}$CO, and C$^{18}$O in either or both J=2-1 and J=3-2 lines in the transition disks around DM Tau, Sz 91, LkCa 15, and HD 34282. We derived line emission surfaces in CO for all disks and in $^{13}$CO for the DM Tau and LkCa 15 disks. With these observations, we do not resolve the vertical structure of C$^{18}$O in any disk, which is instead consistent with C$^{18}$O emission originating from the midplane. Both the J=2-1 and J=3-2 lines show similar heights. Using the derived emission surfaces, we computed radial and vertical gas temperature distributions for each disk, including empirical temperature models for the DM Tau and LkCa 15 disks. After combining our sample with literature sources, we find that $^{13}$CO line emitting heights are also tentatively linked with source characteristics, e.g., stellar host mass, gas temperature, disk size, and show steeper trends than seen in CO emission surfaces. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Comments: 31 pages, 15 figures, accepted for publication in ApJ. Image cubes available at https://zenodo.org/record/7430257

Cold Deuterium Fractionation in the Nearest Planet-Forming Disk

Authors: Carlos E. Muñoz-Romero, Karin I. Öberg, Charles J. Law, Richard Teague, Yuri Aikawa, Jennifer B. Bergner, David J. Wilner, Jane Huang, Viviana V. Guzmán, L. Ilsedore Cleeves

Abstract: Deuterium fractionation provides a window to the thermal history of volatiles in the solar system and protoplanetary disks. While evidence of active molecular deuteration has been observed towards a handful of disks, it remains unclear whether this chemistry affects the composition of forming planetesimals due to limited observational constraints on the radial and vertical distribution of deuterat… ▽ More Deuterium fractionation provides a window to the thermal history of volatiles in the solar system and protoplanetary disks. While evidence of active molecular deuteration has been observed towards a handful of disks, it remains unclear whether this chemistry affects the composition of forming planetesimals due to limited observational constraints on the radial and vertical distribution of deuterated molecules. To shed light on this question, we introduce new ALMA observations of DCO$^+$ and DCN $J=2-1$ at an angular resolution of $0.5"$ (30 au) and combine them with archival data of higher energy transitions towards the protoplanetary disk around TW Hya. We carry out a radial excitation analysis assuming both LTE and non-LTE to localize the physical conditions traced by DCO$^+$ and DCN emission in the disk, thus assessing deuterium fractionation efficiencies and pathways at different disk locations. We find similar disk-averaged column densities of $1.9\times10^{12}$ and $9.8\times10^{11}$ cm$^{-2}$ for DCO$^{+}$ and DCN, with typical kinetic temperatures for both molecules of 20-30K, indicating a common origin near the comet- and planet-forming midplane. The observed DCO$^+$/DCN abundance ratio, combined with recent modeling results, provide tentative evidence of a gas phase C/O enhancement within $<40$ au. Observations of DCO$^+$ and DCN in other disks, as well as HCN and HCO$^+$, will be necessary to place the trends exhibited by TW Hya in context, and fully constrain the main deuteration mechanisms in disks. △ Less

Submitted 5 January, 2023; v1 submitted 13 December, 2022; originally announced December 2022.

Comments: Fixed author information

Clustered Formation of Massive Stars within an Ionized Rotating Disk

Authors: Roberto Galván-Madrid, Qizhou Zhang, Andrés Izquierdo, Charles J. Law, Thomas Peters, Eric Keto, Hauyu Baobab Liu, Paul T. P. Ho, Adam Ginsburg, Carlos Carrasco-González

Abstract: We present ALMA observations with a 800 au resolution and radiative-transfer modelling of the inner part ($r\approx6000$ au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultra-compact (UC) HII region G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich en… ▽ More We present ALMA observations with a 800 au resolution and radiative-transfer modelling of the inner part ($r\approx6000$ au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultra-compact (UC) HII region G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich envelope. The MCMC fits to the data give total stellar masses $M_\star$ from 120 to $200~M_\odot$, with much smaller ionized-gas masses $M_\mathrm{ion-gas} = 0.2$ to $0.25~M_\odot$. The stellar mass is distributed within the gravitational radius $R_g\approx 1000$ to 1500 au, where the ionized gas is bound. The viewing inclination angle from the face-on orientation is $i = 49$ to $56~°$. Radial motions at radii $r > R_g$ converge to $v_{r,0} \approx 8.7$ km/s, or about the speed of sound of ionized gas, indicating that this gas is marginally unbound at most. From additional constraints on the ionizing-photon rate and far-IR luminosity of the region, we conclude that the stellar cluster consists of a few massive stars with $M_\mathrm{star} = 32$ to $60~M_\odot$, or one star in this range of masses accompanied by a population of lower-mass stars. Any active accretion of ionized gas onto the massive (proto)stars is residual. The inferred cluster density is very large, comparable to that reported at similar scales in the Galactic Center. Stellar interactions are likely to occur within the next Myr. △ Less

Submitted 11 December, 2022; originally announced December 2022.

Comments: Accepted to ApJ Letters, December 7 2022

UV-driven Chemistry as a Signpost for Late-stage Planet Formation

Authors: Jenny K. Calahan, Edwin A. Bergin, Arthur D. Bosman, Evan Rich, Sean M. Andrews, Jennifer B. Bergner, L. Ilsedore Cleeves, Viviana V. Guzman, Jane Huang, John D. Ilee, Charles J. Law, Romane Le Gal, Karin I. Oberg, Richard Teague, Catherine Walsh, David J. Wilner, Ke Zhang

Abstract: The chemical reservoir within protoplanetary disks has a direct impact on planetary compositions and the potential for life. A long-lived carbon-and nitrogen-rich chemistry at cold temperatures (<=50K) is observed within cold and evolved planet-forming disks. This is evidenced by bright emission from small organic radicals in 1-10 Myr aged systems that would otherwise have frozen out onto grains w… ▽ More The chemical reservoir within protoplanetary disks has a direct impact on planetary compositions and the potential for life. A long-lived carbon-and nitrogen-rich chemistry at cold temperatures (<=50K) is observed within cold and evolved planet-forming disks. This is evidenced by bright emission from small organic radicals in 1-10 Myr aged systems that would otherwise have frozen out onto grains within 1 Myr. We explain how the chemistry of a planet-forming disk evolves from a cosmic-ray/X-ray-dominated regime to an ultraviolet-dominated chemical equilibrium. This, in turn, will bring about a temporal transition in the chemical reservoir from which planets will accrete. This photochemical dominated gas phase chemistry develops as dust evolves via growth, settling and drift, and the small grain population is depleted from the disk atmosphere. A higher gas-to-dust mass ratio allows for deeper penetration of ultraviolet photons is coupled with a carbon-rich gas (C/O > 1) to form carbon-bearing radicals and ions. This further results in gas phase formation of organic molecules, which then would be accreted by any actively forming planets present in the evolved disk. △ Less

Submitted 11 December, 2022; originally announced December 2022.

Comments: Accepted to Nature Astronomy, Published Dec 8th 2022

The JCMT BISTRO-2 Survey: Magnetic Fields of the Massive DR21 Filament

Authors: Tao-Chung Ching, Keping Qiu, Di Li, Zhiyuan Ren, Shih-Ping Lai, David Berry, Kate Pattle, Ray Furuya, Derek Ward-Thompson, Doug Johnstone, Patrick M. Koch, Chang Won Lee, Thiem Hoang, Tetsuo Hasegawa, Woojin Kwon, Pierre Bastien, Chakali Eswaraiah, Jia-Wei Wang, Kyoung Hee Kim, Jihye Hwang, Archana Soam, A-Ran Lyo, Junhao Liu, Valentin J. M. Le Gouellec, Doris Arzoumanian , et al. (132 additional authors not shown)

Abstract: We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parall… ▽ More We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1--10 pc scales. The magnetic fields revealed in the Planck data are well aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6--1.0 mG in the DR21 filament and $\sim$ 0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by JCMT. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and sub-filaments. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 26 pages, 13 figures, ApJ accepted

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