Showing 1–48 of 48 results for author: Chen, M C

Discovery of optically emitting circumgalactic nebulae around the majority of UV-luminous quasars at intermediate redshift

Authors: Sean D. Johnson, Zhuoqi Will Liu, Jennifer I. Li, Joop Schaye, Jenny E. Greene, Sebastiano Cantalupo, Gwen C. Rudie, Zhijie Qu, Hsiao-Wen Chen, Marc Rafelski, Sowgat Muzahid, Mandy C. Chen, Thierry Contini, Wolfram Kollatschny, Nishant Mishra, Michael Rauch, Patrick Petitjean, Fakhri S. Zahedy

Abstract: We report the discovery of large ionized, [O II] emitting circumgalactic nebulae around the majority of thirty UV luminous quasars at $z=0.4-1.4$ observed with deep, wide-field integral field spectroscopy (IFS) with the Multi-Unit Spectroscopy Explorer (MUSE) by the Cosmic Ultraviolet Baryon Survey (CUBS) and MUSE Quasar Blind Emitters Survey (MUSEQuBES). Among the 30 quasars, seven (23%) exhibit… ▽ More We report the discovery of large ionized, [O II] emitting circumgalactic nebulae around the majority of thirty UV luminous quasars at $z=0.4-1.4$ observed with deep, wide-field integral field spectroscopy (IFS) with the Multi-Unit Spectroscopy Explorer (MUSE) by the Cosmic Ultraviolet Baryon Survey (CUBS) and MUSE Quasar Blind Emitters Survey (MUSEQuBES). Among the 30 quasars, seven (23%) exhibit [O II] emitting nebulae with major axis sizes greater than 100 kpc, twenty greater than 50 kpc (67%), and 27 (90%) greater than 20 kpc. Such large, optically emitting nebulae indicate that cool, dense, and metal-enriched circumgalactic gas is common in the halos of luminous quasars at intermediate redshift. Several of the largest nebulae exhibit morphologies that suggest interaction-related origins. We detect no correlation between the sizes and cosmological dimming corrected surface brightnesses of the nebulae and quasar redshift, luminosity, black hole mass, or radio-loudness, but find a tentative correlation between the nebulae and rest-frame [O II] equivalent width in the quasar spectra. This potential trend suggests a relationship between ISM content and gas reservoirs on CGM scales. The [O II]-emitting nebulae around the $z\approx1$ quasars are smaller and less common than Ly$α$ nebulae around $z\approx3$ quasars. These smaller sizes can be explained if the outer regions of the Ly$α$ halos arise from scattering in more neutral gas, by evolution in the cool CGM content of quasar host halos, by lower-than-expected metallicities on $\gtrsim50$ kpc scales around $z\approx1$ quasars, or by changes in quasar episodic lifetimes between $z=3$ and $1$. △ Less

Submitted 3 April, 2024; v1 submitted 29 March, 2024; originally announced April 2024.

Comments: 18 pages, 5 figures, 2 tables. Accepted for publication in the Astrophysical Journal

The Cosmic Ultraviolet Baryon Survey (CUBS) VIII: Group Environment of the Most Luminous Quasars at $z\approx1$

Authors: Jennifer I. Li, Sean D. Johnson, Erin Boettcher, Sebastiano Cantalupo, Hsiao-Wen Chen, Mandy C. Chen, David R. DePalma, Zhuoqi, Liu, Nishant Mishra, Patrick Petitjean, Zhijie Qu, Gwen C. Rudie, Joop Schaye, Fakhri S. Zahedy

Abstract: We investigate the group-scale environment of 15 luminous quasars (luminosity $L_{\rm 3000}>10^{46}$ erg s$^{-1}$) from the Cosmic Ultraviolet Baryon Survey (CUBS) at redshift $z\approx1$. Using the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph on the Very Large Telescope (VLT), we conduct a deep galaxy redshift survey in the CUBS quasar fields to identify group members and… ▽ More We investigate the group-scale environment of 15 luminous quasars (luminosity $L_{\rm 3000}>10^{46}$ erg s$^{-1}$) from the Cosmic Ultraviolet Baryon Survey (CUBS) at redshift $z\approx1$. Using the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph on the Very Large Telescope (VLT), we conduct a deep galaxy redshift survey in the CUBS quasar fields to identify group members and measure the physical properties of individual galaxies and galaxy groups. We find that the CUBS quasars reside in diverse environments. The majority (11 out of 15) of the CUBS quasars reside in overdense environments with typical halo masses exceeding $10^{13}{\rm M}_{\odot}$, while the remaining quasars reside in moderate-size galaxy groups. No correlation is observed between overdensity and redshift, black hole (BH) mass, or luminosity. Radio-loud quasars (5 out of 15 CUBS quasars) are more likely to be in overdense environments than their radio-quiet counterparts in the sample, consistent with the mean trends from previous statistical observations and clustering analyses. Nonetheless, we also observe radio-loud quasars in moderate groups and radio-quiet quasars in overdense environments, indicating a large scatter in the connection between radio properties and environment. We find that the most UV luminous quasars might be outliers in the stellar mass-to-halo mass relations or may represent departures from the standard single-epoch BH relations. △ Less

Submitted 6 March, 2024; originally announced March 2024.

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

The Cosmic Ultraviolet Baryon Survey (CUBS) VII: on the warm-hot circumgalactic medium probed by O VI and Ne VIII at 0.4 $\lesssim$ z $\lesssim$ 0.7

Authors: Zhijie Qu, Hsiao-Wen Chen, Sean D. Johnson, Gwen C. Rudie, Fakhri S. Zahedy, David DePalma, Joop Schaye, Erin T. Boettcher, Sebastiano Cantalupo, Mandy C. Chen, Claude-André Faucher-Giguère, Jennifer I-Hsiu Li, John S. Mulchaey, Patrick Petitjean, Marc Rafelski

Abstract: This paper presents a newly established sample of 103 unique galaxies or galaxy groups at $0.4\lesssim z\lesssim 0.7$ from the Cosmic Ultraviolet Baryon Survey (CUBS) for studying the warm-hot circumgalactic medium (CGM) probed by both O VI and Ne VIII absorption. The galaxies and associated neighbors are identified at $< 1$ physical Mpc from the sightlines toward 15 CUBS QSOs at… ▽ More This paper presents a newly established sample of 103 unique galaxies or galaxy groups at $0.4\lesssim z\lesssim 0.7$ from the Cosmic Ultraviolet Baryon Survey (CUBS) for studying the warm-hot circumgalactic medium (CGM) probed by both O VI and Ne VIII absorption. The galaxies and associated neighbors are identified at $< 1$ physical Mpc from the sightlines toward 15 CUBS QSOs at $z_{\rm QSO}\gtrsim 0.8$. A total of 30 galaxies or galaxy groups exhibit associated O VI $λλ$ 1031, 1037 doublet absorption within a line-of-sight velocity interval of $\pm250$ km/s, while the rest show no trace of O VI to a detection limit of $\log N_{\rm OVI}/{\rm cm^{-2}}\approx13.7$. Meanwhile, only five galaxies or galaxy groups exhibit the Ne VIII $λλ$ 770,780 doublet absorption, down to a limiting column density of $\log N_{\rm NeVIII}/{\rm cm^{-2}}\approx14.0$. These O VI- and Ne VIII-bearing halos reside in different galaxy environments with stellar masses ranging from $\log M_{\rm star}/M_\odot \approx 8$ to $\approx11.5$. The warm-hot CGM around galaxies of different stellar masses and star formation rates exhibits different spatial profiles and kinematics. In particular, star-forming galaxies with $\log M_{\rm star}/M_\odot\approx9-11$ show a significant concentration of metal-enriched warm-hot CGM within the virial radius, while massive quiescent galaxies exhibit flatter radial profiles of both column densities and covering fractions. In addition, the velocity dispersion of O VI absorption is broad with $σ_v > 40$ km/s for galaxies of $\log M_{\rm star}/M_\odot>9$ within the virial radius, suggesting a more dynamic warm-hot halo around these galaxies. Finally, the warm-hot CGM probed by O VI and Ne VIII is suggested to be the dominant phase in sub-$L^*$ galaxies with $\log M_{\rm star}/M_\odot\approx9-10$ based on their high ionization fractions in the CGM. △ Less

Submitted 14 February, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: Submitted to ApJ after addressing the referee's comments; 28 pages, 16 figures

Infall of material onto the filaments in Barnard 5

Authors: Spandan Choudhury, Jaime E. Pineda, Paola Caselli, Michael Chun-Yuan Chen, Stella S. R. Offner, Maria Teresa Valdivia-Mena

Abstract: Aims. We aim to study the structure and kinematics of the two filaments inside the subsonic core Barnard 5 in Perseus using high-resolution ($\approx$ 2400 au) NH3 data and a multi-component fit analysis. Methods. We used observations of NH3 (1,1) and (2,2) inversion transitions using the Very Large Array (VLA) and the Green Bank Telescope (GBT). We smoothed the data to a beam of 8'' to reliably… ▽ More Aims. We aim to study the structure and kinematics of the two filaments inside the subsonic core Barnard 5 in Perseus using high-resolution ($\approx$ 2400 au) NH3 data and a multi-component fit analysis. Methods. We used observations of NH3 (1,1) and (2,2) inversion transitions using the Very Large Array (VLA) and the Green Bank Telescope (GBT). We smoothed the data to a beam of 8'' to reliably fit multiple velocity components towards the two filamentary structures identified in B5. Results. Along with the core and cloud components, which dominate the flux in the line of sight, we detected two components towards the two filaments showing signs of infall. We also detected two additional components that can possibly trace new material falling into the subsonic core of B5. Conclusions. Following comparison with previous simulations of filament formation scenarios in planar geometry, we conclude that either the formation of the B5 filaments is likely to be rather cylindrically symmetrical or the filaments are magnetically supported. We also estimate infall rates of $1.6\times10^{-4}\,M_\odot\,yr^{-1}$ and $1.8\times10^{-4}\,M_\odot\,yr^{-1}$ (upper limits) for the material being accreted onto the two filaments. At these rates, the filament masses can change significantly during the core lifetime. We also estimate an upper limit of $3.5\times10^{-5}\,M_\odot\,yr^{-1}$ for the rate of possible infall onto the core itself. Accretion of new material onto cores indicates the need for a significant update to current core evolution models, where cores are assumed to evolve in isolation. △ Less

Submitted 28 December, 2023; originally announced December 2023.

Comments: 15 pages, 13 figures, accepted for publication on Astronomy and Astrophysics. The scripts used for the analysis presented in this paper can be found at https://github.com/SpandanCh/Barnard5_filaments

Astronomy as a Field: A Guide for Aspiring Astrophysicists

Authors: Ava Polzin, Yasmeen Asali, Sanah Bhimani, Madison Brady, Mandy C. Chen, Lindsay DeMarchi, Michelle Gurevich, Emily Lichko, Emma Louden, Julie Malewicz, Samantha Pagan, Malena Rice, Zili Shen, Emily Simon, Candice Stauffer, J. Luna Zagorac, Katie Auchettl, Katelyn Breivik, Hsiao-Wen Chen, Deanne Coppejans, Sthabile Kolwa, Raffaella Margutti, Priyamvada Natarajan, Erica Nelson, Kim L. Page , et al. (3 additional authors not shown)

Abstract: This book was created as part of the SIRIUS B VERGE program to orient students to astrophysics as a broad field. The 2023-2024 VERGE program and the printing of this book is funded by the Women and Girls in Astronomy Program via the International Astronomical Union's North American Regional Office of Astronomy for Development and the Heising-Simons Foundation; as a result, this document is written… ▽ More This book was created as part of the SIRIUS B VERGE program to orient students to astrophysics as a broad field. The 2023-2024 VERGE program and the printing of this book is funded by the Women and Girls in Astronomy Program via the International Astronomical Union's North American Regional Office of Astronomy for Development and the Heising-Simons Foundation; as a result, this document is written by women in astronomy for girls who are looking to pursue the field. However, given its universal nature, the material covered in this guide is useful for anyone interested in pursuing astrophysics professionally. △ Less

Submitted 26 December, 2023; v1 submitted 7 December, 2023; originally announced December 2023.

Comments: Introductory guide for students interested in pursuing astrophysics; to be submitted to BAAS

An ensemble study of turbulence in extended QSO nebulae at $z\approx0.5$--1

Authors: Mandy C. Chen, Hsiao-Wen Chen, Michael Rauch, Zhijie Qu, Sean D. Johnson, Joop Schaye, Gwen C. Rudie, Jennifer I-Hsiu Li, Zhuoqi, Liu, Fakhri S. Zahedy, Sebastiano Cantalupo, Erin Boettcher

Abstract: Turbulent motions in the circumgalactic medium (CGM) play a critical role in regulating the evolution of galaxies, yet their detailed characterization remains elusive. Using two-dimensional velocity maps constructed from spatially-extended [OII] and [OIII] emission, Chen et al. (2023b) measured the velocity structure functions (VSFs) of four quasar nebulae at $z\approx\!0.5$--1.1. One of these exh… ▽ More Turbulent motions in the circumgalactic medium (CGM) play a critical role in regulating the evolution of galaxies, yet their detailed characterization remains elusive. Using two-dimensional velocity maps constructed from spatially-extended [OII] and [OIII] emission, Chen et al. (2023b) measured the velocity structure functions (VSFs) of four quasar nebulae at $z\approx\!0.5$--1.1. One of these exhibits a spectacular Kolmogorov relation. Here we carry out an ensemble study using an expanded sample incorporating four new nebulae from three additional QSO fields. The VSFs measured for all eight nebulae are best explained by subsonic turbulence revealed by the line-emitting gas, which in turn strongly suggests that the cool gas ($T\!\sim\!10^4$ K) is dynamically coupled to the hot ambient medium. Previous work demonstrates that the largest nebulae in our sample reside in group environments with clear signs of tidal interactions, suggesting that environmental effects are vital in seeding and enhancing turbulence within the gaseous halos, ultimately promoting the formation of the extended nebulae. No discernible differences are observed in the VSF properties between radio-loud and radio-quiet QSO fields. We estimate the turbulent heating rate per unit volume, $Q_{\rm turb}$, in the QSO nebulae to be $\sim 10^{-26}$--$10^{-22}$ erg cm$^{-3}$ s$^{-1}$ for the cool phase and $\sim 10^{-28}$--$10^{-25}$ erg cm$^{-3}$ s$^{-1}$ for the hot phase. This range aligns with measurements in the intracluster medium and star-forming molecular clouds but is $\sim10^3$ times higher than the $Q_{\rm turb}$ observed inside cool gas clumps on scales $\lesssim1$ kpc using absorption-line techniques. We discuss the prospect of bridging the gap between emission and absorption studies by pushing the emission-based VSF measurements to below $\approx\!10$ kpc. △ Less

Submitted 12 January, 2024; v1 submitted 27 October, 2023; originally announced October 2023.

Comments: 23 pages; 7 figures, and 4 tables in main text; 9 figures in Appendix; accepted by ApJ. Comments welcome

The Cosmic Ultraviolet Baryon Survey: Empirical Characterization of Turbulence in the Cool Circumgalactic Medium

Authors: Hsiao-Wen Chen, Zhijie Qu, Michael Rauch, Mandy C. Chen, Fakhri S. Zahedy, Sean D. Johnson, Joop Schaye, Gwen C. Rudie, Erin Boettcher, Sebastiano Cantalupo, Claude-Andre Faucher-Giguere, Jenny E. Greene, Sebastian Lopez, Robert A. Simcoe

Abstract: This paper reports the first measurement of the relationship between turbulent velocity and cloud size in the diffuse circumgalactic medium (CGM) in typical galaxy halos at redshift z~0.4-1. Through spectrally-resolved absorption profiles of a suite of ionic transitions paired with careful ionization analyses of individual components, cool clumps of size as small as l_cl~1 pc and density lower tha… ▽ More This paper reports the first measurement of the relationship between turbulent velocity and cloud size in the diffuse circumgalactic medium (CGM) in typical galaxy halos at redshift z~0.4-1. Through spectrally-resolved absorption profiles of a suite of ionic transitions paired with careful ionization analyses of individual components, cool clumps of size as small as l_cl~1 pc and density lower than nH = 0.001 cm^-3 are identified in galaxy halos. In addition, comparing the line widths between different elements for kinematically matched components provides robust empirical constraints on the thermal temperature T and the non-thermal motions bNT, independent of the ionization models. On average, bNT is found to increase with l_cl following bNT \propto l_cl^0.3 over three decades in spatial scale from l_cl~1 pc to l_cl~1 kpc. Attributing the observed bNT to turbulent motions internal to the clumps, the best-fit bNT-l_cl relation shows that the turbulence is consistent with Kolmogorov at <1 kpc with a roughly constant energy transfer rate per unit mass of epsilon~0.003 cm^2 s^-3 and a dissipation time scale of <~ 100 Myr. No significant difference is found between massive quiescent and star-forming halos in the sample on scales less than 1 kpc. While the inferred epsilon is comparable to what is found in CIV absorbers at high redshift, it is considerably smaller than observed in star-forming gas or in extended line-emitting nebulae around distant quasars. A brief discussion of possible sources to drive the observed turbulence in the cool CGM is presented. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: Accepted for publication in ApJ Letters; 10 pages, 4 figures

The first comprehensive study of a giant nebula around a radio-quiet quasar in the $z < 1$ Universe

Authors: Zhuoqi Will Liu, Sean D. Johnson, Jennifer I-Hsiu Li, Gwen C. Rudie, Joop Schaye, Hsiao-Wen Chen, Jarle Brinchmann, Sebastiano Cantalupo, Mandy C. Chen, Wolfram Kollatschny, Michael V. Maseda, Nishant Mishra, Sowgat Muzahid

Abstract: We present the first comprehensive study of a giant, $\approx \! \! 70$ kpc-scale nebula around a radio-quiet quasar at $z<1$. The analysis is based on deep integral field spectroscopy with MUSE of the field of HE$\,$0238$-$1904, a luminous quasar at $z=0.6282$. The nebula emits strongly in $\mathrm{[O \, II]}$, $\rm H β$, and $\mathrm{[O \, III]}$, and the quasar resides in an unusually overdense… ▽ More We present the first comprehensive study of a giant, $\approx \! \! 70$ kpc-scale nebula around a radio-quiet quasar at $z<1$. The analysis is based on deep integral field spectroscopy with MUSE of the field of HE$\,$0238$-$1904, a luminous quasar at $z=0.6282$. The nebula emits strongly in $\mathrm{[O \, II]}$, $\rm H β$, and $\mathrm{[O \, III]}$, and the quasar resides in an unusually overdense environment for a radio-quiet system. The environment likely consists of two groups which may be merging, and in total have an estimated dynamical mass of $M_{\rm dyn}\approx 4\times 10^{13}$ to $10^{14}\ {\rm M_\odot}$. The nebula exhibits largely quiescent kinematics and irregular morphology. The nebula may arise primarily through interaction-related stripping of circumgalactic and interstellar medium (CGM/ISM) of group members, with some potential contributions from quasar outflows. The simultaneous presence of the giant nebula and a radio-quiet quasar in a rich environment suggests a correlation between such circum-quasar nebulae and environmental effects. This possibility can be tested with larger samples. The upper limits on the electron number density implied by the $\mathrm{[O \, II]}$ doublet ratio range from $\log(n_{\rm e, \, [O \, II]} / \mathrm{cm^{-3}}) < 1.2$ to $2.8$. However, assuming a constant quasar luminosity and negligible projection effects, the densities implied from the measured line ratios between different ions (e.g., $\mathrm{[O\,II]}$, $\mathrm{[O\,III]}$, and $\mathrm{[Ne\,V]}$) and photoionization simulations are often $10{-}400$ times larger. This large discrepancy can be explained by quasar variability on a timescale of $\approx 10^4{-}10^5$ years. △ Less

Submitted 30 October, 2023; v1 submitted 31 August, 2023; originally announced September 2023.

Comments: 19 pages, 9 figures, 3 tables; Accepted by MNRAS

Alignment of dense molecular core morphology and velocity gradients with ambient magnetic fields

Authors: A. Pandhi, R. K. Friesen, L. Fissel, J. E. Pineda, P. Caselli, M. C-Y. Chen, J. Di Francesco, A. Ginsburg, H. Kirk, P. C. Myers, S. S. R. Offner, A. Punanova, F. Quan, E. Redaelli, E. Rosolowsky, S. Scibelli, Y. M. Seo, Y. Shirley

Abstract: Studies of dense core morphologies and their orientations with respect to gas flows and the local magnetic field have been limited to only a small sample of cores with spectroscopic data. Leveraging the Green Bank Ammonia Survey alongside existing sub-millimeter continuum observations and Planck dust polarization, we produce a cross-matched catalogue of 399 dense cores with estimates of core morph… ▽ More Studies of dense core morphologies and their orientations with respect to gas flows and the local magnetic field have been limited to only a small sample of cores with spectroscopic data. Leveraging the Green Bank Ammonia Survey alongside existing sub-millimeter continuum observations and Planck dust polarization, we produce a cross-matched catalogue of 399 dense cores with estimates of core morphology, size, mass, specific angular momentum, and magnetic field orientation. Of the 399 cores, 329 exhibit 2D $\mathrm{v}_\mathrm{LSR}$ maps that are well fit with a linear gradient, consistent with rotation projected on the sky. We find a best-fit specific angular momentum and core size relationship of $J/M \propto R^{1.82 \pm 0.10}$, suggesting that core velocity gradients originate from a combination of solid body rotation and turbulent motions. Most cores have no preferred orientation between the axis of core elongation, velocity gradient direction, and the ambient magnetic field orientation, favouring a triaxial and weakly magnetized origin. We find, however, strong evidence for a preferred anti-alignment between the core elongation axis and magnetic field for protostellar cores, revealing a change in orientation from starless and prestellar populations that may result from gravitational contraction in a magnetically-regulated (but not dominant) environment. We also find marginal evidence for anti-alignment between the core velocity gradient and magnetic field orientation in the L1228 and L1251 regions of Cepheus, suggesting a preferred orientation with respect to magnetic fields may be more prevalent in regions with locally ordered fields. △ Less

Submitted 24 July, 2023; originally announced July 2023.

Comments: 33 pages, 28 figures, accepted to MNRAS

The Cosmic Ultraviolet Baryon Survey (CUBS) VI: Connecting Physical Properties of the Cool Circumgalactic Medium to Galaxies at $z\approx 1$

Authors: Zhijie Qu, Hsiao-Wen Chen, Gwen C. Rudie, Sean D. Johnson, Fakhri S. Zahedy, David DePalma, Erin Boettcher, Sebastiano Cantalupo, Mandy C. Chen, Kathy L. Cooksey, Claude-André Faucher-Giguère, Jennifer I-Hsiu Li, Sebastian Lopez, Joop Schaye, Robert A. Simcoe

Abstract: This paper presents a new sample of 19 unique galaxies and galaxy groups at $z\approx1$ from the CUBS program, which is designated as the CUBSz1 sample. In this CUBSz1 sample, nine galaxies or galaxy groups show absorption features, while ten systems do not have detectable absorption with 2-$σ$ upper limits of log$N$(HeI)/cm$^{-2}\lesssim 13.5$ and log$N$(OV)/cm$^{-2}\lesssim 13.3$. Environmental… ▽ More This paper presents a new sample of 19 unique galaxies and galaxy groups at $z\approx1$ from the CUBS program, which is designated as the CUBSz1 sample. In this CUBSz1 sample, nine galaxies or galaxy groups show absorption features, while ten systems do not have detectable absorption with 2-$σ$ upper limits of log$N$(HeI)/cm$^{-2}\lesssim 13.5$ and log$N$(OV)/cm$^{-2}\lesssim 13.3$. Environmental properties of the galaxies, including galaxy overdensities, the total stellar mass and gravitational potential summed over all nearby neighbors, and the presence of local ionizing sources, are found to have a significant impact on the observed CGM absorption properties. Specifically, massive galaxies and galaxies in overdense regions exhibit a higher rate of incidence of absorption. At the same time, the observed CGM absorption properties in galaxy groups appear to be driven by the galaxy closest to the QSO sightline, rather than by the most massive galaxy or by mass-weighted properties. We introduce a total projected gravitational potential $ψ$, defined as $-ψ/G =\sum M_{\rm halo}/d_{\rm proj}$ summed over all group members, to characterize the overall galaxy environment. This projected gravitational potential correlates linearly with the maximum density detected in each sightline, consistent with higher-pressure gas being confined in deeper gravitational potential wells. In addition, we find that the radial profile of cool gas density exhibits a general decline from the inner regions to the outskirts, being in pressure balance with the hot halo. Finally, we note that the ionizing flux from nearby galaxies can generate an elevated $N$(HI)/$N$(HeI) ratio, which in turn provides a unique diagnostic of possible local sources contributing to the ionizing radiation field. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: 17 pages, 11 figures and 1 table in the main text; 40 figures and 1 table in the appendix; accepted for publication in MNRAS

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)

Velocity-Coherent Substructure in TMC-1: Inflow and Fragmentation

Authors: Simon E. T. Smith, Rachel Friesen, Antoine Marchal, Jaime E. Pineda, Paola Caselli, Michael Chun-Yuan Chen, Spandan Choudhury, James Di Francesco, Adam Ginsburg, Helen Kirk, Chris Matzner, Anna Punanova, Samantha Scibelli, Yancy Shirley

Abstract: Filamentary structures have been found nearly ubiquitously in molecular clouds and yet their formation and evolution is still poorly understood. We examine a segment of Taurus Molecular Cloud 1 (TMC-1) that appears as a single, narrow filament in continuum emission from dust. We use the Regularized Optimization for Hyper-Spectral Analysis (ROHSA), a Gaussian decomposition algorithm which enforces… ▽ More Filamentary structures have been found nearly ubiquitously in molecular clouds and yet their formation and evolution is still poorly understood. We examine a segment of Taurus Molecular Cloud 1 (TMC-1) that appears as a single, narrow filament in continuum emission from dust. We use the Regularized Optimization for Hyper-Spectral Analysis (ROHSA), a Gaussian decomposition algorithm which enforces spatial coherence when fitting multiple velocity components simultaneously over a data cube. We analyze HC$_5$N (9-8) line emission as part of the Green Bank Ammonia Survey (GAS) and identify three velocity-coherent components with ROHSA. The two brightest components extend the length of the filament, while the third component is fainter and clumpier. The brightest component has a prominent transverse velocity gradient of $2.7 \pm 0.1$ km s$^{-1}$ pc$^{-1}$ that we show to be indicative of gravitationally induced inflow. In the second component, we identify regularly spaced emission peaks along its length. We show that the local minima between pairs of adjacent HC$_5$N peaks line up closely with submillimetre continuum emission peaks, which we argue is evidence for fragmentation along the spine of TMC-1. While coherent velocity components have been described as separate physical structures in other star-forming filaments, we argue that the two bright components identified in HC$_5$N emission in TMC-1 are tracing two layers in one filament: a lower density outer layer whose material is flowing under gravity towards the higher density inner layer of the filament. △ Less

Submitted 6 February, 2023; v1 submitted 18 November, 2022; originally announced November 2022.

Comments: 17 pages, 8 figures; Accepted for publication to MNRAS

Indirect Measurements of Gas Velocities in Galaxy Clusters: Effects of Ellipticity and Cluster Dynamic State

Authors: Irina Zhuravleva, Mandy C. Chen, Eugene Churazov, Alexander A. Schekochihin, Congyao Zhang, Daisuke Nagai

Abstract: While awaiting direct velocity measurement of gas motions in the hot intracluster medium, we rely on indirect probes, including gas perturbations in galaxy clusters. Using a sample of $\sim 80$ clusters in different dynamic states from Omega500 cosmological simulations, we examine scaling relations between the fluctuation amplitudes of gas density, $δρ/ρ$, pressure, $δP/P$, X-ray surface brightnes… ▽ More While awaiting direct velocity measurement of gas motions in the hot intracluster medium, we rely on indirect probes, including gas perturbations in galaxy clusters. Using a sample of $\sim 80$ clusters in different dynamic states from Omega500 cosmological simulations, we examine scaling relations between the fluctuation amplitudes of gas density, $δρ/ρ$, pressure, $δP/P$, X-ray surface brightness, Sunyaev-Zeldovich (SZ) y-parameter, and the characteristic Mach number of gas motions, $M_{\rm 1d}$. In relaxed clusters, accounting for halo ellipticities reduces $δρ/ρ$ or $δP/P$ by a factor of up to 2 within $r_{500c}$. We confirm a strong linear correlation between $δρ/ρ$ (or $δP/P$) and $M_{\rm 1d}$ in relaxed clusters, with the proportionality coefficient $η\approx 1$. For unrelaxed clusters, the correlation is less strong and has a larger $η\approx 1.3\pm 0.5$ ($1.5\pm0.5$) for $δρ/ρ$ ($δP/P$). Examination of the power-law scaling of $M_{\rm 1d}$ with $δρ/ρ$ shows that it is almost linear for relaxed clusters, while for the unrelaxed ones, it is closer to $δρ/ρ\propto M_{\rm 1d}^2$, supporting an increasing role of non-linear terms and compressive modes. In agreement with previous studies, we observe a strong correlation of $M_{\rm 1d}$ with radius. Correcting for these correlations leaves a residual scatter in $M_{\rm 1d}$ of $\sim 4 (7)$ per cent for relaxed (perturbed) clusters. Hydrostatic mass bias correlates with $M_{\rm 1d}$ as strongly as with $δρ/ρ$ in relaxed clusters. The residual scatters after correcting for derived trends is $\sim 6-7$ per cent. These predictions can be verified with existing X-ray and SZ observations of galaxy clusters combined with forthcoming velocity measurements with X-ray microcalorimeters. △ Less

Submitted 20 October, 2022; originally announced October 2022.

Comments: 15 pages, 11 figures, submitted

Empirical constraints on the turbulence in QSO host nebulae from velocity structure function measurements

Authors: Mandy C. Chen, Hsiao-Wen Chen, Michael Rauch, Zhijie Qu, Sean D. Johnson, Jennifer I-Hsiu Li, Joop Schaye, Gwen C. Rudie, Fakhri S. Zahedy, Erin Boettcher, Kathy L. Cooksey, Sebastiano Cantalupo

Abstract: We present the first empirical constraints on the turbulent velocity field of the diffuse circumgalactic medium around four luminous QSOs at $z\!\approx\!0.5$--1.1. Spatially extended nebulae of $\approx\!50$--100 physical kpc in diameter centered on the QSOs are revealed in [OII]$λλ\,3727,3729$ and/or [OIII]$λ\,5008$ emission lines in integral field spectroscopic observations obtained using MUSE… ▽ More We present the first empirical constraints on the turbulent velocity field of the diffuse circumgalactic medium around four luminous QSOs at $z\!\approx\!0.5$--1.1. Spatially extended nebulae of $\approx\!50$--100 physical kpc in diameter centered on the QSOs are revealed in [OII]$λλ\,3727,3729$ and/or [OIII]$λ\,5008$ emission lines in integral field spectroscopic observations obtained using MUSE on the VLT. We measure the second- and third-order velocity structure functions (VSFs) over a range of scales, from $\lesssim\!5$ kpc to $\approx\!20$--50 kpc, to quantify the turbulent energy transfer between different scales in these nebulae. While no constraints on the energy injection and dissipation scales can be obtained from the current data, we show that robust constraints on the power-law slope of the VSFs can be determined after accounting for the effects of atmospheric seeing, spatial smoothing, and large-scale bulk flows. Out of the four QSO nebulae studied, one exhibits VSFs in spectacular agreement with the Kolmogorov law, expected for isotropic, homogeneous, and incompressible turbulent flows. The other three fields exhibit a shallower decline in the VSFs from large to small scales. However, with a limited dynamic range in the spatial scales in seeing-limited data, no constraints can be obtained for the VSF slopes of these three nebulae. For the QSO nebula consistent with the Kolmogorov law, we determine a turbulence energy cascade rate of $\approx\!0.2$ cm$^{2}$ s$^{-3}$. We discuss the implication of the observed VSFs in the context of QSO feeding and feedback in the circumgalactic medium. △ Less

Submitted 3 November, 2022; v1 submitted 9 September, 2022; originally announced September 2022.

Comments: 20 pages, 6 figures in main text + 9 figures in appendix, 4 tables; accepted by MNRAS

Directly tracing cool filamentary accretion over >100 kpc into the interstellar medium of a quasar host at z=1

Authors: Sean D. Johnson, Joop Schaye, Gregory L. Walth, Jennifer I-Hsiu Li, Gwen C. Rudie, Hsiao-Wen Chen, Mandy C. Chen, Benoît Epinat, Massimo Gaspari, Sebastiano Cantalupo, Wolfram Kollatschny, Zhuoqi, Liu, Sowgat Muzahid

Abstract: We report the discovery of giant (50-100 kpc) [O II] emitting nebulae with the Multi-Unit Spectroscopic Explorer (MUSE) in the field of TXS 0206-048, a luminous quasar at z=1.13. Down-the-barrel UV spectra of the quasar show absorption at velocities coincident with those of the extended nebulae, enabling new insights into inflows and outflows around the quasar host. One nebula exhibits a filamenta… ▽ More We report the discovery of giant (50-100 kpc) [O II] emitting nebulae with the Multi-Unit Spectroscopic Explorer (MUSE) in the field of TXS 0206-048, a luminous quasar at z=1.13. Down-the-barrel UV spectra of the quasar show absorption at velocities coincident with those of the extended nebulae, enabling new insights into inflows and outflows around the quasar host. One nebula exhibits a filamentary morphology extending over 120 kpc from the halo toward the quasar and intersecting with another nebula surrounding the quasar host with a radius of 50 kpc. This is the longest cool filament observed to-date and arises at higher redshift and in a less massive system than those in cool-core clusters. The filamentary nebula has line-of-sight velocities >300 km/s from nearby galaxies but matches that of the nebula surrounding the quasar host where they intersect, consistent with accretion of cool inter- or circum-galactic medium or cooling hot halo gas. The kinematics of the nebulae surrounding the quasar host are unusual and complex, with redshifted and blueshifted spiral-like structures. The emission velocities at 5-10 kpc from the quasar match those of inflowing absorbing gas observed in UV spectra of the quasar. Together, the extended nebulae and associated redshifted absorption represent a compelling case of cool, filamentary gas accretion from halo scales into the extended interstellar medium and toward the nucleus of a massive quasar host. The inflow rate implied by the combined emission and absorption constraints is well below levels required to sustain the quasar's radiative luminosity, suggesting anisotropic or variable accretion. △ Less

Submitted 11 November, 2022; v1 submitted 9 September, 2022; originally announced September 2022.

Comments: 10 pages (main text), 3 figures, 2 tables, and one appendix figure.. Accepted to ApJ Letters

The Cosmic Ultraviolet Baryon Survey (CUBS) V: On the Thermodynamic Properties of the Cool Circumgalactic Medium at $z < 1$

Authors: Zhijie Qu, Hsiao-Wen Chen, Gwen C. Rudie, Fakhri S. Zahedy, Sean D. Johnson, Erin Boettcher, Sebastiano Cantalupo, Mandy C. Chen, Kathy L. Cooksey, David DePalma, Claude-André Faucher-Giguère, Michael Rauch, Joop Schaye, Robert A. Simcoe

Abstract: This paper presents a systematic study of the photoionization and thermodynamic properties of the cool circumgalactic medium (CGM) as traced by rest-frame ultraviolet absorption lines around 26 galaxies at redshift $z\lesssim1$. The study utilizes both high-quality far-ultraviolet and optical spectra of background QSOs and deep galaxy redshift surveys to characterize the gas density, temperature,… ▽ More This paper presents a systematic study of the photoionization and thermodynamic properties of the cool circumgalactic medium (CGM) as traced by rest-frame ultraviolet absorption lines around 26 galaxies at redshift $z\lesssim1$. The study utilizes both high-quality far-ultraviolet and optical spectra of background QSOs and deep galaxy redshift surveys to characterize the gas density, temperature, and pressure of individual absorbing components and to resolve their internal non-thermal motions. The derived gas density spans more than three decades, from $\log (n_{\rm H}/{\rm cm^{-3}}) \approx -4$ to $-1$, while the temperature of the gas is confined in a narrow range of $\log (T/{\rm K})\approx 4.3\pm 0.3$. In addition, a weak anti-correlation between gas density and temperature is observed, consistent with the expectation of the gas being in photoionization equilibrium. Furthermore, decomposing the observed line widths into thermal and non-thermal contributions reveals that more than 30% of the components at $z\lesssim 1$ exhibit line widths driven by non-thermal motions, in comparison to $<20$% found at $z\approx 2$-3. Attributing the observed non-thermal line widths to intra-clump turbulence, we find that massive quenched galaxies on average exhibit higher non-thermal broadening/turbulent energy in their CGM compared to star-forming galaxies at $z\lesssim 1$. Finally, strong absorption features from multiple ions covering a wide range of ionization energy (e.g., from Mg II to O IV) can be present simultaneously in a single absorption system with kinematically aligned component structure, but the inferred pressure in different phases may differ by a factor of $\approx 10$. △ Less

Submitted 2 September, 2022; originally announced September 2022.

Comments: 16 pages, 4 figures, 1 table, accepted for publication in MNRAS

Turbulence and Accretion: a High-resolution Study of the B5 Filaments

Authors: Michael Chun-Yuan Chen, James Di Francesco, Jaime E. Pineda, Stella S. Offner, Rachel K. Friesen

Abstract: High-resolution observations of the Perseus B5 "core" have previously revealed that this subsonic region actually consists of several filaments that are likely in the process of forming a quadruple stellar system. Since subsonic filaments are thought to be produced at the $\sim 0.1$ pc sonic scale by turbulent compression, a detailed kinematic study is crucial to test such a scenario in the contex… ▽ More High-resolution observations of the Perseus B5 "core" have previously revealed that this subsonic region actually consists of several filaments that are likely in the process of forming a quadruple stellar system. Since subsonic filaments are thought to be produced at the $\sim 0.1$ pc sonic scale by turbulent compression, a detailed kinematic study is crucial to test such a scenario in the context of core and star formation. Here we present a detailed kinematic follow-up study of the B5 filaments at a 0.009 pc resolution using the VLA and GBT combined observations fitted with multi-component spectral models. Using precisely identified filament spines, we find a remarkable resemblance between the averaged width profiles of each filament and Plummer-like functions, with filaments possessing FWHM widths of $\sim 0.03$ pc. The velocity dispersion profiles of the filaments also show decreasing trends towards the filament spines. Moreover, the velocity gradient field in B5 appears to be locally well ordered ($\sim 0.04$ pc) but globally complex, with kinematic behaviors suggestive of inhomogeneous turbulent accretion onto filaments and longitudinal flows towards a local overdensity along one of the filaments. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: 27 pages, 12 figures, published in The Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 935, Issue 1, id.57, 16 pp (2022)

RELICS: Strong Lens Model of SMACSJ0723.3-7327

Authors: Keren Sharon, Mandy C. Chen, Guillaume Mahler, Dan Coe, the RELICS Collaboration

Abstract: We present the details of a strong lens model of SMACS J0723.3-7327, which was made public as part of the data and high level science products (HLSP) release of the RELICS HST treasury program (Reionization Lensing Cluster Survey; GO-14096, PI: Coe). The model products were made available on the Mikulski Archive for Space Telescopes (MAST) via 10.17909/T9SP45 in 2017. Here, we provide the list of… ▽ More We present the details of a strong lens model of SMACS J0723.3-7327, which was made public as part of the data and high level science products (HLSP) release of the RELICS HST treasury program (Reionization Lensing Cluster Survey; GO-14096, PI: Coe). The model products were made available on the Mikulski Archive for Space Telescopes (MAST) via 10.17909/T9SP45 in 2017. Here, we provide the list of constraints that were used in the HST-based RELICS lens model, as well as other information related to our modeling choices, which were not published with the data and HLSP release. This model was computed with Lenstool, used multiple images of 8 sources, with no spectroscopic redshifts. The image plane RMS was 0".58. △ Less

Submitted 17 August, 2022; originally announced August 2022.

Comments: 4 pages, 1 figure, 1 table. Submitted to ApJS

The Davis-Chandrasekhar-Fermi Method Revisited

Authors: Che-Yu Chen, Zhi-Yun Li, Renato R. Mazzei, Jinsoo Park, Laura M. Fissel, Michael C. -Y. Chen, Richard I. Klein, Pak Shing Li

Abstract: Despite the rich observational results on interstellar magnetic fields in star-forming regions, it is still unclear how dynamically significant the magnetic fields are at varying physical scales, because direct measurement of the field strength is observationally difficult. The Davis-Chandrasekhar-Fermi (DCF) method has been the most commonly used method to estimate the magnetic field strength fro… ▽ More Despite the rich observational results on interstellar magnetic fields in star-forming regions, it is still unclear how dynamically significant the magnetic fields are at varying physical scales, because direct measurement of the field strength is observationally difficult. The Davis-Chandrasekhar-Fermi (DCF) method has been the most commonly used method to estimate the magnetic field strength from polarization data. It is based on the assumption that gas turbulent motion is the driving source of field distortion via linear Alfvén waves. In this work, using MHD simulations of star-forming clouds, we test the validity of the assumption underlying the DCF method by examining its accuracy in the real 3D space. Our results suggest that the DCF relation between turbulent kinetic energy and magnetic energy fluctuation should be treated as a statistical result instead of a local property. We then develop and investigate several modifications to the original DCF method using synthetic observations, and propose new recipes to improve the accuracy of DCF-derived magnetic field strength. We further note that the biggest uncertainty in the DCF analysis may come from the linewidth measurement instead of the polarization observation, especially since the line-of-sight gas velocity can be used to estimate the gas volume density, another critical parameter in the DCF method. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: 20 pages, 17 figures, accepted for publication in MNRAS

Discovery of a damped Ly$α$ absorber originating in a spectacular interacting dwarf galaxy pair at $z = 0.026$

Authors: Erin Boettcher, Neeraj Gupta, Hsiao-Wen Chen, Mandy C. Chen, Gyula I. G. Józsa, Gwen C. Rudie, Sebastiano Cantalupo, Sean D. Johnson, S. A. Balashev, Françoise Combes, Kathy L. Cooksey, Claude-André Faucher-Giguère, Jens-Kristian Krogager, Sebastian Lopez, Emmanuel Momjian, Pasquier Noterdaeme, Patrick Petitjean, Marc Rafelski, Raghunathan Srianand, Gregory L. Walth, Fakhri S. Zahedy

Abstract: We present the discovery of neutral gas detected in both damped Ly$α$ absorption (DLA) and HI 21-cm emission outside of the stellar body of a galaxy, the first such detection in the literature. A joint analysis between the Cosmic Ultraviolet Baryon Survey and the MeerKAT Absorption Line Survey reveals an HI bridge connecting two interacting dwarf galaxies (log… ▽ More We present the discovery of neutral gas detected in both damped Ly$α$ absorption (DLA) and HI 21-cm emission outside of the stellar body of a galaxy, the first such detection in the literature. A joint analysis between the Cosmic Ultraviolet Baryon Survey and the MeerKAT Absorption Line Survey reveals an HI bridge connecting two interacting dwarf galaxies (log$(M_{\text{star}}/\text{M}_{\odot}) = 8.5 \pm 0.2$) that host a $z = 0.026$ DLA with log[$N$(HI)/cm$^{-2}$]$ = 20.60 \pm 0.05$ toward the QSO J2339-5523 ($z_{\text{QSO}} = 1.35$). At impact parameters of $d = 6$ and $33$ kpc, the dwarf galaxies have no companions more luminous than $\approx 0.05L_{*}$ within at least $Δv = \pm 300$ km s$^{-1}$ and $d \approx 350$ kpc. HI 21-cm emission is spatially coincident with the DLA at the 2-3$σ$ level per spectral channel over several adjacent beams. However, HI 21-cm absorption is not detected against the radio-bright QSO; if the background UV and radio sources are spatially aligned, the gas is either warm or clumpy (with spin temperature to covering factor ratio $T_{s}/f_{c} > 1880$ K). VLT-MUSE observations demonstrate that the $α$-element abundance of the ionized ISM is consistent with the DLA ($\approx 10$% solar), suggesting that the neutral gas envelope is perturbed ISM gas. This study showcases the impact of dwarf-dwarf interactions on the physical and chemical state of neutral gas outside of star-forming regions. In the SKA era, joint UV and HI 21-cm analyses will be critical for connecting the cosmic neutral gas content to galaxy environments. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: Accepted for publication in ApJL. 16 pages, 4 figures, 2 tables

Velocity structure of the 50 pc-long NGC 6334 filamentary cloud: Hints of multiple compressions and their impact on the cloud properties?

Authors: Doris Arzoumanian, Delphine Russeil, Annie Zavagno, Michael Chun-Yuan Chen, Philippe André, Shu-ichiro Inutsuka, Yoshiaki Misugi, Álvaro Sánchez-Monge, Peter Schilke, Alexander Men'shchikov, Mikito Kohno

Abstract: [Abridged] The interstellar medium is observed to be organised in filamentary structures, as well as neutral (HI) and ionized (HII) bubbles. The expanding nature of these bubbles makes them shape their surroundings and possibly play a role in the formation and evolution of interstellar filaments. We present APEX $^{13}$CO and C$^{18}$O(2-1) observations of the NGC 6334 molecular cloud. We investig… ▽ More [Abridged] The interstellar medium is observed to be organised in filamentary structures, as well as neutral (HI) and ionized (HII) bubbles. The expanding nature of these bubbles makes them shape their surroundings and possibly play a role in the formation and evolution of interstellar filaments. We present APEX $^{13}$CO and C$^{18}$O(2-1) observations of the NGC 6334 molecular cloud. We investigate the gas velocity structure along and across the 50 pc-long cloud and towards 75 identified velocity-coherent-filaments (VCFs). We measure a wealth of velocity gradients along the VCFs. We derive the column density and velocity power spectra of the VCFs. These power spectra are well represented with power laws showing similar slopes for both quantities (with a mean of about -2), albeit some differ by up to a factor of two. The position velocity diagrams perpendicular to three VCFs show the V-shaped velocity pattern, corresponding to a bent structure in velocity space with the filament at the tip of the V surrounded by an extended structure connected to it with a velocity gradient. This velocity structure is qualitatively similar to that resulting from numerical simulations of filament formation from large-scale compression from propagating shock fronts. In addition, the radial profiles perpendicular to these VCFs hint to small-scale internal impacts from neighbouring HII bubbles. The observed opposite curvature in velocity space towards the VCFs points to various origins of large-scale external compressions from propagating HI bubbles. This suggests the plausible importance of multiple HI compressions, separated in space and time, in the formation and evolution of molecular clouds and their star formation history. These latter atomic compressions due to past and distant star formation events are complemented by the impact of HII bubbles from present time and local star formation activity. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 660, A56 (2022)

The Cosmic Ultraviolet Baryon Survey (CUBS) IV: The Complex Multiphase Circumgalactic Medium as Revealed by Partial Lyman Limit Systems

Authors: Thomas J. Cooper, Gwen C. Rudie, Hsiao-Wen Chen, Sean D. Johnson, Fakhri S. Zahedy, Mandy C. Chen, Erin Boettcher, Gregory L. Walth, Sebastiano Cantalupo, Kathy L. Cooksey, Claude-André Faucher-Giguère, Jenny E. Greene, Sebastian Lopez, John S. Mulchaey, Steven V. Penton, Patrick Petitjean, Mary E. Putman, Marc Rafelski, Michael Rauch, Joop Schaye, Robert A. Simcoe

Abstract: We present a detailed study of two partial Lyman limit systems (pLLSs) of neutral hydrogen column density $N_\mathrm{H\,I}\approx(1-3)\times10^{16}\,\mathrm{cm}^{-2}$ discovered at $z=0.5$ in the Cosmic Ultraviolet Baryon Survey (CUBS). Available far-ultraviolet spectra from the Hubble Space Telescope Cosmic Origins Spectrograph and optical echelle spectra from MIKE on the Magellan Telescopes enab… ▽ More We present a detailed study of two partial Lyman limit systems (pLLSs) of neutral hydrogen column density $N_\mathrm{H\,I}\approx(1-3)\times10^{16}\,\mathrm{cm}^{-2}$ discovered at $z=0.5$ in the Cosmic Ultraviolet Baryon Survey (CUBS). Available far-ultraviolet spectra from the Hubble Space Telescope Cosmic Origins Spectrograph and optical echelle spectra from MIKE on the Magellan Telescopes enable a comprehensive ionization analysis of diffuse circumgalactic gas based on resolved kinematics and abundance ratios of atomic species spanning five different ionization stages. These data provide unambiguous evidence of kinematically aligned multi-phase gas that masquerades as a single-phase structure and can only be resolved by simultaneous accounting of the full range of observed ionic species. Both systems are resolved into multiple components with inferred $α$-element abundance varying from $[α/\text{H}]\approx\!{-0.8}$ to near solar and densities spanning over two decades from $\log n_\mathrm{H}\mathrm{cm}^{-3}\approx\!-2.2$ to $<-4.3$. Available deep galaxy survey data from the CUBS program taken with VLT/MUSE, Magellan/LDSS3-C and Magellan/IMACS reveal that the $z=0.47$ system is located 55 kpc from a star-forming galaxy with prominent Balmer absorption of stellar mass $M_\star\approx2\times10^{10}M_\odot$, while the $z=0.54$ system resides in an over-dense environment of 11 galaxies within 750 kpc in projected distance, with the most massive being a luminous red galaxy of $M_\star\approx2\times10^{11}M_\odot$ at 375 kpc. The study of these two pLLSs adds to an emerging picture of the complex, multiphase circumgalactic gas that varies in chemical abundances and density on small spatial scales in diverse galaxy environments. The inhomogeneous nature of metal enrichment and density revealed in observations must be taken into account in theoretical models of diffuse halo gas. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: Accepted for publication in MNRAS, 27 pages, 21 figures

Are massive dense clumps truly sub-virial? A new analysis using Gould Belt ammonia data

Authors: Ayushi Singh, Christopher D. Matzner, Rachel K. Friesen, Peter G. Martin, Jaime E. Pineda, Erik W. Rosolowsky, Felipe Alves, Ana Chacón-Tanarro, Hope How-Huan Chen, Michael Chun-Yuan Chen, Spandan Choudhury, James Di Francesco, Jared Keown, Helen Kirk, Anna Punanova, Youngmin Seo, Yancy Shirley, Adam Ginsburg, Stella S. R. Offner, Héctor G. Arce, Paola Caselli, Alyssa A. Goodman, Philip C. Myers, Elena Redaelli

Abstract: Dynamical studies of dense structures within molecular clouds often conclude that the most massive clumps contain too little kinetic energy for virial equilibrium, unless they are magnetized to an unexpected degree. This raises questions about how such a state might arise, and how it might persist long enough to represent the population of massive clumps. In an effort to re-examine the origins of… ▽ More Dynamical studies of dense structures within molecular clouds often conclude that the most massive clumps contain too little kinetic energy for virial equilibrium, unless they are magnetized to an unexpected degree. This raises questions about how such a state might arise, and how it might persist long enough to represent the population of massive clumps. In an effort to re-examine the origins of this conclusion, we use ammonia line data from the Green Bank Ammonia Survey and Planck-calibrated dust emission data from Herschel to estimate the masses and kinetic and gravitational energies for dense clumps in the Gould Belt clouds. We show that several types of systematic error can enhance the appearance of low kinetic-to-gravitational energy ratios: insufficient removal of foreground and background material; ignoring the kinetic energy associated with velocity differences across a resolved cloud; and over-correcting for stratification when evaluating the gravitational energy. Using an analysis designed to avoid these errors, we find that the most massive Gould Belt clumps harbor virial motions, rather than sub-virial ones. As a byproduct, we present a catalog of masses, energies, and virial energy ratios for 85 Gould Belt clumps. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: Submitted to ApJ

The Cosmic Ultraviolet Baryon Survey (CUBS) -- III. Physical properties and elemental abundances of Lyman limit systems at $z<1$

Authors: Fakhri S. Zahedy, Hsiao-Wen Chen, Thomas M. Cooper, Erin T. Boettcher, Sean D. Johnson, Gwen C. Rudie, Mandy C. Chen, Sebastiano Cantalupo, Kathy L. Cooksey, Claude-André Faucher-Giguère, Jenny E. Greene, Sebastian Lopez, John S. Mulchaey, Steven V. Penton, Patrick Petitjean, Mary E. Putman, Marc Rafelski, Michael Rauch, Joop Schaye, Robert A. Simcoe, Gregory L. Walth

Abstract: (Abridged) We present a systematic investigation of physical conditions and elemental abundances in four optically thick Lyman-limit systems (LLSs) at $z=0.36-0.6$ discovered within the Cosmic Ultraviolet Baryon Survey (CUBS). CUBS LLSs exhibit multi-component kinematic structure and a complex mix of multiphase gas, with associated metal transitions from multiple ionization states that span severa… ▽ More (Abridged) We present a systematic investigation of physical conditions and elemental abundances in four optically thick Lyman-limit systems (LLSs) at $z=0.36-0.6$ discovered within the Cosmic Ultraviolet Baryon Survey (CUBS). CUBS LLSs exhibit multi-component kinematic structure and a complex mix of multiphase gas, with associated metal transitions from multiple ionization states that span several hundred km/s in line-of-sight velocity. Specifically, higher column density components (log N(HI)>16) in all four absorbers comprise dynamically cool gas with $\langle T \rangle =(2\pm1) \times10^4\,$K and modest non-thermal broadening of $5\pm3\,$ km/s. The high quality of the QSO absorption spectra allows us to infer the physical conditions of the gas, using a detailed ionization modeling that takes into account the resolved component structures of HI and metal transitions. The range of inferred gas densities indicates that these absorbers consist of spatially compact clouds with a median line-of-sight thickness of $160^{+140}_{-50}$ pc. While obtaining robust metallicity constraints for the low-density, highly ionized phase remains challenging due to the uncertain N(HI), we demonstrate that the cool-phase gas in LLSs has a median metallicity of $\mathrm{[α/H]_{1/2}}=-0.7^{+0.1}_{-0.2}$, with a 16-84 percentile range of $\mathrm{[α/H]}=(-1.3,-0.1)$. Furthermore, the wide range of inferred elemental abundance ratios ($\mathrm{[C/α]}$, $\mathrm{[N/α]}$, and $\mathrm{[Fe/α]}$) indicate a diversity of chemical enrichment histories. Combining the absorption data with deep galaxy survey data characterizing the galaxy environment of these absorbers, we discuss the physical connection between star-forming regions in galaxies and diffuse gas associated with optically thick absorption systems in the $z<1$ circumgalactic medium. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Comments: Accepted for publication in MNRAS after a minor revision. 18 pages, 11 figures, and a nine-page Appendix with four additional figures

Resolved galactic superwinds reconstructed around their host galaxies at z>3

Authors: Mandy C. Chen, Hsiao-Wen Chen, Max Gronke, Michael Rauch, Tom Broadhurst

Abstract: This paper presents a detailed analysis of two giant Lyman-alpha (Lya) arcs detected near known galaxies at z=3.038 and z=3.754 lensed by the massive cluster MACS 1206 (z=0.44). The Lya nebulae revealed in deep MUSE observations exhibit a double-peak profile with a dominant red peak that indicates expansion/outflowing motions. One of the arcs stretches over 1' around the Einstein radius of the clu… ▽ More This paper presents a detailed analysis of two giant Lyman-alpha (Lya) arcs detected near known galaxies at z=3.038 and z=3.754 lensed by the massive cluster MACS 1206 (z=0.44). The Lya nebulae revealed in deep MUSE observations exhibit a double-peak profile with a dominant red peak that indicates expansion/outflowing motions. One of the arcs stretches over 1' around the Einstein radius of the cluster, resolving the velocity field of the line-emitting gas on kpc scales around a group of three star-forming galaxies of 0.3-1.6L* at z=3.038. The second arc spans 15'' in size, roughly centered around a pair of low-mass Lya emitters of ~0.03L* at z=3.754. All three galaxies in the z=3.038 group exhibit prominent damped Lya absorption (DLA) and several metal absorption lines, in addition to nebular emission lines such as HeII1640 and CIII]1906,1908. Extended Lya emission appears to emerge from star-forming regions to larger distances with suppressed surface brightness at the center of each galaxy, suggesting the presence of dusty outflowing cones of size 1-5 kpc across. There are significant spatial variations in the Lya line profile, consistent with the presence of a steep negative velocity gradient in a continuous flow of high column density gas from star-forming regions into a low-density halo environment. While the observed UV nebular line ratios show no evidence of AGN activity in the galaxies, the observed Lya signals can be explained by a combination of resonant scattering and recombination radiation due to photoionization by ionizing photons escaping from the nearby star-forming regions. These observations provide the most detailed insights yet into the kinematics of galactic superwinds associated with star-forming galaxies thought to be responsible for the chemical enrichment in the intergalactic medium. △ Less

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

Comments: 29 pages, 12 figures, 6 tables; this is the journal accepted version

The Cosmic Ultraviolet Baryon Survey (CUBS) II: Discovery of an H$_{2}$-Bearing DLA in the Vicinity of an Early-Type Galaxy at z = 0.576

Authors: Erin Boettcher, Hsiao-Wen Chen, Fakhri S. Zahedy, Thomas J. Cooper, Sean D. Johnson, Gwen C. Rudie, Mandy C. Chen, Patrick Petitjean, Sebastiano Cantalupo, Kathy L. Cooksey, Claude-André Faucher-Giguère, Jenny E. Greene, Sebastian Lopez, John S. Mulchaey, Steven V. Penton, Mary E. Putman, Marc Rafelski, Michael Rauch, Joop Schaye, Robert A. Simcoe, Gregory L. Walth

Abstract: We report the serendipitous detection of an H$_{2}$-bearing damped Lyman-$α$ absorber at z = 0.576 in the spectrum of the QSO J0111-0316 in the Cosmic Ultraviolet Baryon Survey. Spectroscopic observations from HST-COS in the far-ultraviolet reveal a damped absorber with log[N(HI)/cm^-2] = 20.1 +/- 0.2 and log[N(H$_{2}$)/cm^-2] = 18.97 (-0.06, +0.05). The diffuse molecular gas is found in two veloc… ▽ More We report the serendipitous detection of an H$_{2}$-bearing damped Lyman-$α$ absorber at z = 0.576 in the spectrum of the QSO J0111-0316 in the Cosmic Ultraviolet Baryon Survey. Spectroscopic observations from HST-COS in the far-ultraviolet reveal a damped absorber with log[N(HI)/cm^-2] = 20.1 +/- 0.2 and log[N(H$_{2}$)/cm^-2] = 18.97 (-0.06, +0.05). The diffuse molecular gas is found in two velocity components separated by dv = 60 km/s, with >99.9% of the total H$_{2}$ column density concentrated in one component. At a metallicity of $\approx$ 50% of solar, there is evidence for Fe enhancement and dust depletion, with a dust-to-gas ratio $κ_{\text{O}} \approx$ 0.4. A galaxy redshift survey conducted with IMACS and LDSS-3C on Magellan reveals an overdensity of nine galaxies at projected distance d <= 600 proper kpc (pkpc) and line-of-sight velocity offset dv$_{g}$ <= 300 km/s from the absorber. The closest is a massive, early-type galaxy at d = 41 pkpc which contains $\approx$ 70% of the total stellar mass identified at d <= 310 pkpc of the H$_{2}$ absorber. The close proximity of the H$_{2}$-bearing gas to the quiescent galaxy and the Fe-enhanced chemical abundance pattern of the absorber suggest a physical connection, in contrast to a picture in which DLAs are primarily associated with gas-rich dwarfs. This case study illustrates that deep galaxy redshift surveys are needed to gain insight into the diverse environments that host dense and potentially star-forming gas. △ Less

Submitted 19 March, 2021; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: 19 pages, 6 figures. Accepted for publication in ApJ

Ubiquitous $\rm NH_3$ supersonic component in L1688 coherent cores

Authors: Spandan Choudhury, Jaime E. Pineda, Paola Caselli, Adam Ginsburg, Stella S. R. Offner, Erik Rosolowsky, Rachel K. Friesen, Felipe O. Alves, Ana Chacón-Tanarro, Anna Punanova, Elena Redaelli, Helen Kirk, Philip C. Myers, Peter G. Martin, Yancy Shirley, Michael Chun-Yuan Chen, Alyssa A. Goodman, James Di Francesco

Abstract: Context : Star formation takes place in cold dense cores in molecular clouds. Earlier observations have found that dense cores exhibit subsonic non-thermal velocity dispersions. In contrast, CO observations show that the ambient large-scale cloud is warmer and has supersonic velocity dispersions. Aims : We aim to study the ammonia ($\rm NH_3$) molecular line profiles with exquisite sensitivity tow… ▽ More Context : Star formation takes place in cold dense cores in molecular clouds. Earlier observations have found that dense cores exhibit subsonic non-thermal velocity dispersions. In contrast, CO observations show that the ambient large-scale cloud is warmer and has supersonic velocity dispersions. Aims : We aim to study the ammonia ($\rm NH_3$) molecular line profiles with exquisite sensitivity towards the coherent cores in L1688 in order to study their kinematical properties in unprecedented detail. Methods : We used $\rm NH_3$ (1,1) and (2,2) data from the first data release (DR1) in the Green Bank Ammonia Survey (GAS). We first smoothed the data to a larger beam of 1' to obtain substantially more extended maps of velocity dispersion and kinetic temperature, compared to the DR1 maps. We then identified the coherent cores in the cloud and analysed the averaged line profiles towards the cores. Results : For the first time, we detected a faint (mean $\rm NH_3$(1,1) peak brightness $<$0.25 K in $T_{MB}$), supersonic component towards all the coherent cores in L1688. We fitted two components, one broad and one narrow, and derived the kinetic temperature and velocity dispersion of each component. The broad components towards all cores have supersonic linewidths ($\mathcal{M}_S \ge 1$). This component biases the estimate of the narrow dense core component's velocity dispersion by $\approx$28% and the kinetic temperature by $\approx$10%, on average, as compared to the results from single-component fits. Conclusions : Neglecting this ubiquitous presence of a broad component towards all coherent cores causes the typical single-component fit to overestimate the temperature and velocity dispersion. This affects the derived detailed physical structure and stability of the cores estimated from $\rm NH_3$ observations. △ Less

Submitted 20 July, 2020; v1 submitted 14 July, 2020; originally announced July 2020.

Comments: Accepted for publication in Astronomy & Astrophysics on 06/07/2020. 15 pages, 16 figures, 1 table. Language edits from previous version

Journal ref: A&A 640, L6 (2020)

Geometric Support for Dark Matter by an Unaligned Einstein Ring in Abell 3827

Authors: Mandy C. Chen, Tom Broadhurst, Jeremy Lim, Sandor M. Molnar, Jose M. Diego, Masamune Oguri, Lilian L. Lee

Abstract: The non-detection of dark matter (DM) particles in increasingly stringent laboratory searches has encouraged alternative gravity theories where gravity is sourced only from visible matter. Here, we consider whether such theories can pass a two-dimensional test posed by gravitational lensing -- to reproduce a particularly detailed Einstein ring in the core of the galaxy cluster Abell 3827. We find… ▽ More The non-detection of dark matter (DM) particles in increasingly stringent laboratory searches has encouraged alternative gravity theories where gravity is sourced only from visible matter. Here, we consider whether such theories can pass a two-dimensional test posed by gravitational lensing -- to reproduce a particularly detailed Einstein ring in the core of the galaxy cluster Abell 3827. We find that when we require the lensing mass distribution to strictly follow the shape (ellipticity and position angle) of the light distribution of cluster member galaxies, intracluster stars, and the X-ray emitting intracluster medium, we cannot reproduce the Einstein ring, despite allowing the mass-to-light ratios of these visible components to freely vary with radius to mimic alternative gravity theories. Alternatively, we show that the detailed features of the Einstein ring are accurately reproduced by allowing a smooth, freely oriented DM halo in the lens model, with relatively small contributions from the visible components at a level consistent with their observed brightnesses. This dominant DM component is constrained to have the same orientation as the light from the intracluster stars, indicating that the intracluster stars trace the gravitational potential of this component. The Einstein ring of Abell 3827 therefore presents a new challenge for alternative gravity theories: not only must such theories find agreement between the total lensing mass and visible mass, but they must also find agreement between the projected sky distribution of the lensing mass and that of the visible matter, a more stringent test than has hitherto been posed by lensing data. △ Less

Submitted 10 July, 2020; originally announced July 2020.

Comments: 20 pages, 7 figures, 8 tables. Accepted for publication in ApJ

The JCMT BISTRO Survey: Magnetic Fields Associated with a Network of Filaments in NGC 1333

Authors: Yasuo Doi, Tetsuo Hasegawa, Ray S. Furuya, Simon Coudé, Charles L. H. Hull, Doris Arzoumanian, Pierre Bastien, Michael Chun-Yuan Chen, James di Francesco, Rachel Friesen, Martin Houde, Shu-ichiro Inutsuka, Steve Mairs, Masafumi Matsumura, Takashi Onaka, Sarah Sadavoy, Yoshito Shimajiri, Mehrnoosh Tahani, Kohji Tomisaka, Chakali Eswaraiah, Patrick M. Koch, Kate Pattle, Chang Won Lee, Motohide Tamura, David Berry , et al. (113 additional authors not shown)

Abstract: We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary… ▽ More We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multiscale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~1 pc and remains continuous from the scales of filaments (~0.1 pc) to that of protostellar envelopes (~0.005 pc or ~1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data toward the filament network. △ Less

Submitted 21 July, 2020; v1 submitted 30 June, 2020; originally announced July 2020.

Comments: 33 pages, 21 figures, accepted for publication in ApJ

The Cosmic Ultraviolet Baryon Survey (CUBS) I. Overview and the diverse environments of Lyman limit systems at z<1

Authors: Hsiao-Wen Chen, Fakhri S. Zahedy, Erin Boettcher, Thomas M. Cooper, Sean D. Johnson, Gwen C. Rudie, Mandy C. Chen, Gregory L. Walth, Sebastiano Cantalupo, Kathy L. Cooksey, Claude-Andre Faucher-Gigu`ere, Jenny E. Greene, Sebastian Lopez, John S. Mulchaey, Steven V. Penton, Patrick Petitjean, Mary E. Putman, Marc Rafelski, Michael Rauch, Joop Schaye, Robert A. Simcoe, Benjamin J. Weiner

Abstract: We present initial results from the Cosmic Ultraviolet Baryon Survey (CUBS). CUBS is designed to map diffuse baryonic structures at redshift z<~1 using absorption-line spectroscopy of 15 UV-bright QSOs with matching deep galaxy survey data. CUBS QSOs are selected based on their NUV brightness to avoid biases against the presence of intervening Lyman Limit Systems (LLSs) at zabs<1. We report five n… ▽ More We present initial results from the Cosmic Ultraviolet Baryon Survey (CUBS). CUBS is designed to map diffuse baryonic structures at redshift z<~1 using absorption-line spectroscopy of 15 UV-bright QSOs with matching deep galaxy survey data. CUBS QSOs are selected based on their NUV brightness to avoid biases against the presence of intervening Lyman Limit Systems (LLSs) at zabs<1. We report five new LLSs of log N(HI)/cm^-2 >~ 17.2 over a total redshift survey pathlength of dz=9.3, and a number density of n(z)=0.43 (-0.18, +0.26). Considering all absorbers with log N(HI)/cm^-2 > 16.5 leads to n(z)=1.08 (-0.25, +0.31) at z<1. All LLSs exhibit a multi-component structure and associated metal transitions from multiple ionization states such as CII, CIII, MgII, SiII, SiIII, and OVI absorption. Differential chemical enrichment levels as well as ionization states are directly observed across individual components in three LLSs. We present deep galaxy survey data obtained using the VLT-MUSE integral field spectrograph and the Magellan Telescopes, reaching sensitivities necessary for detecting galaxies fainter than 0.1L* at d<~300 physical kpc (pkpc) in all five fields. A diverse range of galaxy properties is seen around these LLSs, from a low-mass dwarf galaxy pair, a co-rotating gaseous halo/disk, a star-forming galaxy, a massive quiescent galaxy, to a galaxy group. The closest galaxies have projected distances ranging from d=15 to 72 pkpc and intrinsic luminosities from ~0.01L* to ~3L*. Our study shows that LLSs originate in a variety of galaxy environments and trace gaseous structures with a broad range of metallicities. △ Less

Submitted 16 June, 2020; v1 submitted 5 May, 2020; originally announced May 2020.

Comments: 26 pages, 14 figures, MNRAS in press

Relative Alignment between Dense Molecular Cores and Ambient Magnetic Field: The Synergy of Numerical Models and Observations

Authors: Che-Yu Chen, Erica A. Behrens, Jasmin E. Washington, Laura M. Fissel, Rachel K. Friesen, Zhi-Yun Li, Jaime E. Pineda, Adam Ginsburg, Helen Kirk, Samantha Scibelli, Felipe Alves, Elena Redaelli, Paola Caselli, Anna Punanova, James Di Francesco, Erik Rosolowsky, Stella S. R. Offner, Peter G. Martin, Ana Chacón-Tanarro, Hope H. -H. Chen, Michael C. -Y. Chen, Jared Keown, Youngmin Seo, Yancy Shirley, Hector G. Arce , et al. (4 additional authors not shown)

Abstract: The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular cl… ▽ More The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud. △ Less

Submitted 24 March, 2020; originally announced March 2020.

Comments: 18 pages, 11 figures, accepted for publication in MNRAS

Velocity-coherent Filaments in NGC 1333: Evidence for Accretion Flow?

Authors: Michael Chun-Yuan Chen, James Di Francesco, Erik Rosolowsky, Jared Keown, Jaime E. Pineda, Rachel K. Friesen, Paola Caselli, How-Huan Chen, Christopher D. Matzner, Stella S. Offner, Anna Punanova, Elena Redaelli, Samantha Scibelli, Yancy Shirley

Abstract: Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficie… ▽ More Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficient, robust, and automatic method to fit ammonia lines with multiple velocity components, generalizable to other molecular species. We also present CRISPy, a Python package to identify filament spines in 3D images (e.g., position-position-velocity cubes), along with a complementary technique to sort fitted velocity components into velocity-coherent filaments. In NGC 1333, we find a wealth of velocity gradient structures on a beam-resolved scale of ~0.05 pc. Interestingly, these local velocity gradients are not randomly oriented with respect to filament spines and their perpendicular, i.e., radial, component decreases in magnitude towards the spine for many filaments. Together with remarkably constant velocity gradients on larger scales along many filaments, these results suggest a scenario in which gas falling onto filaments is progressively damped and redirected to flow along these filaments. △ Less

Submitted 26 February, 2020; originally announced February 2020.

Comments: Accepted to ApJ

Journal ref: ApJ 891 (2020) 84

CLOVER: Convnet Line-fitting Of Velocities in Emission-line Regions

Authors: Jared Keown, James Di Francesco, Hossen Teimoorinia, Erik Rosolowsky, Michael Chun-Yuan Chen

Abstract: When multiple star-forming gas structures overlap along the line-of-sight and emit optically thin emission at significantly different radial velocities, the emission can become non-Gaussian and often exhibits two distinct peaks. Traditional line-fitting techniques can fail to account adequately for these double-peaked profiles, providing inaccurate cloud kinematics measurements. We present a new m… ▽ More When multiple star-forming gas structures overlap along the line-of-sight and emit optically thin emission at significantly different radial velocities, the emission can become non-Gaussian and often exhibits two distinct peaks. Traditional line-fitting techniques can fail to account adequately for these double-peaked profiles, providing inaccurate cloud kinematics measurements. We present a new method called Convnet Line-fitting Of Velocities in Emission-line Regions (CLOVER) for distinguishing between one-component, two-component, and noise-only emission lines using 1D convolutional neural networks trained with synthetic spectral cubes. CLOVER utilizes spatial information in spectral cubes by predicting on $3\times3$ pixel sub-cubes, using both the central pixel's spectrum and the average spectrum over the $3\times3$ grid as input. On an unseen set of 10,000 synthetic spectral cubes in each predicted class, CLOVER has classification accuracies of $\sim99\%$ for the one-component class and $\sim97\%$ for the two-component class. For the noise-only class, which is analogous to a signal-to-noise cutoff of four for traditional line-fitting methods, CLOVER has classification accuracy of $100\%$. CLOVER also has exceptional performance on real observations, correctly distinguishing between the three classes across a variety of star-forming regions. In addition, CLOVER quickly and accurately extracts kinematics directly from spectra identified as two-component class members. Moreover, we show that CLOVER is easily scalable to emission lines with hyperfine splitting, making it an attractive tool in the new era of large-scale NH$_3$ and N$_2$H$^+$ mapping surveys. △ Less

Submitted 18 September, 2019; originally announced September 2019.

Comments: Accepted for publication in ApJ

KFPA Examinations of Young STellar Object Natal Environments (KEYSTONE): Hierarchical Ammonia Structures in Galactic Giant Molecular Clouds

Authors: Jared Keown, James Di Francesco, Erik Rosolowsky, Ayushi Singh, Charles Figura, Helen Kirk, L. D. Anderson, Michael Chun-Yuan Chen, Davide Elia, Rachel Friesen, Adam Ginsburg, A. Marston, Stefano Pezzuto, Eugenio Schisano, Sylvain Bontemps, Paola Caselli, Hong-Li Liu, Steven Longmore, Frederique Motte, Philip C. Myers, Stella S. R. Offner, Patricio Sanhueza, Nicola Schneider, Ian Stephens, James Urquhart , et al. (1 additional authors not shown)

Abstract: We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includ… ▽ More We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includes the NH$_3$ (1,1) and (2,2) maps for each cloud, which are modeled to produce maps of kinetic temperature, centroid velocity, velocity dispersion, and ammonia column density. Median cloud kinetic temperatures range from $11.4\pm2.2$ K in the coldest cloud (MonR1) to $23.0\pm6.5$ K in the warmest cloud (M17). Using dendrograms on the NH$_3$ (1,1) integrated intensity maps, we identify 856 dense gas clumps across the eleven clouds. Depending on the cloud observed, $40-100\%$ of the clumps are aligned spatially with filaments identified in H$_2$ column density maps derived from SED-fitting of dust continuum emission. A virial analysis reveals that 523 of the 835 clumps ($\sim63\%$) with mass estimates are bound by gravity alone. We find no significant difference between the virial parameter distributions for clumps aligned with the dust-continuum filaments and those unaligned with filaments. In some clouds, however, hubs or ridges of dense gas with unusually high mass and low virial parameters are located within a single filament or at the intersection of multiple filaments. These hubs and ridges tend to host water maser emission, multiple 70$μ$m-detected protostars, and have masses and radii above an empirical threshold for forming massive stars. △ Less

Submitted 29 August, 2019; v1 submitted 27 August, 2019; originally announced August 2019.

Comments: Accepted for publication in ApJ

The JCMT BISTRO Survey: The Magnetic Field of the Barnard 1 Star-Forming Region

Authors: Simon Coudé, Pierre Bastien, Martin Houde, Sarah Sadavoy, Rachel Friesen, James Di Francesco, Doug Johnstone, Steve Mairs, Tetsuo Hasegawa, Woojin Kwon, Shih-Ping Lai, Keping Qiu, Derek Ward-Thompson, David Berry, Michael Chun-Yuan Chen, Jason Fiege, Erica Franzmann, Jennifer Hatchell, Kevin Lacaille, Brenda C. Matthews, Gerald H. Moriarty-Schieven, Andy Pon, Philippe André, Doris Arzoumanian, Yusuke Aso , et al. (96 additional authors not shown)

Abstract: We present the POL-2 850 $μ$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We th… ▽ More We present the POL-2 850 $μ$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $μ$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $μ$G and a criticality criterion $λ_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component. △ Less

Submitted 23 April, 2019; v1 submitted 15 April, 2019; originally announced April 2019.

Comments: 22 pages, 7 figures, accepted in the Astrophysical Journal

The JCMT BISTRO Survey: The Magnetic Field In The Starless Core $ρ$ Ophiuchus C

Authors: Junhao Liu, Keping Qiu, David Berry, James Di Francesco, Pierre Bastien, Patrick M. Koch, Ray S. Furuya, Kee-Tae Kim, Simon Coudé, Chang Won Lee, Archana Soam, Chakali Eswaraiah, Di Li, Jihye Hwang, A-Ran Lyo, Kate Pattle, Tetsuo Hasegawa, Woojin Kwon, Shih-Ping Lai, Derek Ward-Thompson, Tao-Chung Ching, Zhiwei Chen, Qilao Gu, Dalei Li, Hua-bai Li , et al. (106 additional authors not shown)

Abstract: We report 850~$μ$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $ρ$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic… ▽ More We report 850~$μ$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $ρ$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $μ$G, 136 $\pm$ 69 $μ$G, and 213 $\pm$ 115 $μ$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties. △ Less

Submitted 20 February, 2019; originally announced February 2019.

Comments: 22 pages, 12 figures

JCMT BISTRO survey: Magnetic Fields within the Hub-Filament Structure in IC 5146

Authors: Jia-Wei Wang, Shih-Ping Lai, Chakali Eswaraiah, Kate Pattle, James Di Francesco, Doug Johnstone, Patrick M. Koch, Tie Liu, Motohide Tamura, Ray S. Furuya, Takashi Onaka, Derek Ward-Thompson, Archana Soam, Kee-Tae Kim, Chang Won Lee, Chin-Fei Lee, Steve Mairs, Doris Arzoumanian, Gwanjeong Kim, Thiem Hoang, Jihye Hwang, Sheng-Yuan Liu, David Berry, Pierre Bastien, Tetsuo Hasegawa , et al. (108 additional authors not shown)

Abstract: We present the 850 $μ$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scal… ▽ More We present the 850 $μ$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfvénic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important. △ Less

Submitted 27 March, 2019; v1 submitted 14 December, 2018; originally announced December 2018.

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

Droplets I: Pressure-Dominated Sub-0.1 pc Coherent Structures in L1688 and B18

Authors: Hope How-Huan Chen, Jaime E. Pineda, Alyssa A. Goodman, Andreas Burkert, Stella S. R. Offner, Rachel K. Friesen, Philip C. Myers, Felipe Alves, Hector G. Arce, Paola Caselli, Ana Chacon-Tanarro, Michael Chun-Yuan Chen, James Di Francesco, Adam Ginsburg, Jared Keown, Helen Kirk, Peter G. Martin, Christopher Matzner, Anna Punanova, Elena Redaelli, Erik Rosolowsky, Samantha Scibelli, Young Min Seo, Yancy Shirley, Ayushi Singh

Abstract: We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transiti… ▽ More We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process. △ Less

Submitted 15 May, 2019; v1 submitted 26 September, 2018; originally announced September 2018.

Comments: Accepted by ApJ in April, 2019

Journal ref: 2019ApJ...877...93C

Magnetic fields towards Ophiuchus-B derived from SCUBA-2 polarization measurements

Authors: Archana Soam, Kate Pattle, Derek Ward-Thompson, Chang Won Lee, Sarah Sadavoy, Patrick M. Koch, Gwanjeong Kim, Jungmi Kwon, Woojin Kwon, Doris Arzoumanian, David Berry, Thiem Hoang, Motohide Tamura, Sang-Sung Lee, Tie Liu, Kee-Tae Kim, Doug Johnstone, Fumitaka Nakamura, A-Ran Lyo, Takashi Onaka, Jongsoo Kim, Ray S. Furuya, Tetsuo Hasegawa, Shih-Ping Lai, Pierre Bastien , et al. (99 additional authors not shown)

Abstract: We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fi… ▽ More We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $μ$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $λ$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical. △ Less

Submitted 7 June, 2018; v1 submitted 16 May, 2018; originally announced May 2018.

Comments: 17 pages, 10 figures, Accepted for publication in ApJ

A Likely Super Massive Black Hole Revealed by its Einstein Radius in Hubble Frontier Fields Images

Authors: Mandy C. Chen, Tom Broadhurst, Jeremy Lim, Jose M. Diego, Youichi Ohyama, Holland Ford, Narciso Benítez

Abstract: At cosmological distances, gravitational lensing can provide a direct measurement of supermassive black hole (SMBH) masses irrespective of their luminosities. Here, we directly estimate the mass of a SMBH in the brightest cluster galaxy (BCG) of MACS J1149+2223.5 at $z=0.54$ through one of the multiply-lensed images of a background spiral galaxy at $z=1.49$ projected close to the BCG. In this part… ▽ More At cosmological distances, gravitational lensing can provide a direct measurement of supermassive black hole (SMBH) masses irrespective of their luminosities. Here, we directly estimate the mass of a SMBH in the brightest cluster galaxy (BCG) of MACS J1149+2223.5 at $z=0.54$ through one of the multiply-lensed images of a background spiral galaxy at $z=1.49$ projected close to the BCG. In this particular image, an intrinsically compact region in one of the spiral arms is lensed into an arc that curves towards the BCG center. This arc has a radius of curvature of only $\sim$0."6, betraying the presence of a local compact deflector. Its curvature is most simply reproduced by a point-like object with a mass of $8.4^{+4.3}_{-1.8}\times10^{9}M_\odot$, similar to SMBH masses in local elliptical galaxies having comparable luminosities. The SMBH is noticeably offset by $4.4\pm0.3$ kpc from the BCG light centre, plausibly the result of a kick imparted $\sim2.0\times10^7$ years ago during the merger of two SMBHs, placing it just beyond the stellar core. A similar curvature can be produced by replacing the offset SMBH with a compact galaxy having a mass of $\sim2\times 10^{10}M_\odot$ within a cutoff radius of $<4$ kpc, and an unusually large $M/L>50(M/L)_\odot$ to make it undetectable in the deep Hubble Frontiers Fields image, at or close to the cluster redshift; such a lensing galaxy, however, perturbs the adjacent lensed images in an undesirable way. △ Less

Submitted 20 August, 2018; v1 submitted 14 May, 2018; originally announced May 2018.

Comments: Published version

Journal ref: The Astrophysical Journal, 863:135 (24pp), 2018 August 20

A First Look at BISTRO Observations of The $ρ$ Oph-A core

Authors: Jungmi Kwon, Yasuo Doi, Motohide Tamura, Masafumi Matsumura, Kate Pattle, David Berry, Sarah Sadavoy, Brenda C. Matthews, Derek Ward-Thompson, Tetsuo Hasegawa, Ray S. Furuya, Andy Pon, James Di Francesco, Doris Arzoumanian, Saeko S. Hayashi, Koji S. Kawabata, Takashi Onaka, Minho Choi, Miju Kang, Thiem Hoang, Chang Won Lee, Sang-Sung Lee, Hong-Li Liu, Tie Liu, Shu-Ichiro Inutsuka , et al. (97 additional authors not shown)

Abstract: We present 850 $μ$m imaging polarimetry data of the $ρ$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We iden… ▽ More We present 850 $μ$m imaging polarimetry data of the $ρ$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $ρ$ Oph A core; some of them can be part of a coherent magnetic field structure in the $ρ$ Oph region. The results are consistent with previous observations of the brightest regions of $ρ$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components. △ Less

Submitted 24 April, 2018; originally announced April 2018.

Comments: 22 pages, 11 figures, accepted for publication in ApJ

The Green Bank Ammonia Survey: Unveiling the Dynamics of the Barnard 59 star-forming Clump

Authors: E. Redaelli, F. O. Alves, P. Caselli, J. E. Pineda, R. K. Friesen, A. Chacón-Tanarro, C. D. Matzner, A. Ginsburg, E. Rosolowsky, J. Keown, S. S. R. Offner, J. Di Francesco, H. Kirk, P. C. Myers, A. Hacar, A. Cimatti, H. H. Chen, M. C. Chen, Y. M. Seo, K. I. Lee

Abstract: Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe N… ▽ More Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe Nebula, to achieve a comprehensive view of the kinematic properties of the source. These information were derived by simultaneously fitting ammonia inversion transition lines (1,1) and (2,2). Our analysis unveils the imprint of protostellar feedback, such as increasing line widths, temperature and turbulent motions in our molecular data. Combined with complementary observations of dust thermal emission, we estimate that the core is gravitationally bound following a virial analysis. If the core is not contracting, another source of internal pressure, most likely the magnetic field, is supporting it against gravitational collapse and limits its star formation efficiency. △ Less

Submitted 6 November, 2017; originally announced November 2017.

Comments: 18 pages, 18 figures

The Green Bank Ammonia Survey: Observations of Hierarchical Dense Gas Structures in Cepheus-L1251

Authors: Jared Keown, James Di Francesco, Helen Kirk, Rachel K. Friesen, Jaime E. Pineda, Erik Rosolowsky, Adam Ginsburg, Stella S. R. Offner, Paola Caselli, Felipe Alves, Ana Chacón-Tanarro, Anna Punanova, Elena Redaelli, Young Min Seo, Christopher D. Matzner, Michael Chun-Yuan Chen, Alyssa A. Goodman, How-Huan Chen, Yancy Shirley, Ayushi Singh, Hector G. Arce, Peter Martin, Philip C. Myers

Abstract: We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc… ▽ More We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores. △ Less

Submitted 12 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJ

The Green Bank Ammonia Survey: Dense Cores Under Pressure in Orion A

Authors: Helen Kirk, Rachel K. Friesen, Jaime E. Pineda, Erik Rosolowsky, Stella S. R. Offner, Christopher D. Matzner, Philip C. Myers, James Di Francesco, Paola Caselli, Felipe O. Alves, Ana Chacón-Tanarro, How-Huan Chen, Michael Chun-Yuan Chen, Jared Keown, Anna Punanova, Young Min Seo, Yancy Shirley, Adam Ginsburg, Christine Hall, Ayushi Singh, Héctor G. Arce, Alyssa A. Goodman, Peter Martin, Elena Redaelli

Abstract: We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity… ▽ More We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined. △ Less

Submitted 17 August, 2017; originally announced August 2017.

Comments: Accepted for publication in ApJ Aug 10, 2017

First results from BISTRO -- a SCUBA-2 polarimeter survey of the Gould Belt

Authors: Derek Ward-Thompson, Kate Pattle, Pierre Bastien, Ray S. Furuya, Woojin Kwon, Shih-Ping Lai, Keping Qiu, David Berry, Minho Choi, Simon Coudé, James Di Francesco, Thiem Hoang, Erica Franzmann, Per Friberg, Sarah F. Graves, Jane S. Greaves, Martin Houde, Doug Johnstone, Jason M. Kirk, Patrick M. Koch, Jungmi Kwon, Chang Won Lee, Di Li, Brenda C. Matthews, Joseph C. Mottram , et al. (89 additional authors not shown)

Abstract: We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will a… ▽ More We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments. △ Less

Submitted 27 April, 2017; originally announced April 2017.

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

The Green Bank Ammonia Survey (GAS): First Results of NH3 mapping the Gould Belt

Authors: Rachel K. Friesen, Jaime E. Pineda, Erik Rosolowsky, Felipe Alves, Ana Chacón-Tanarro, Hope How-Huan Chen, Michael Chun-Yuan Chen, James Di Francesco, Jared Keown, Helen Kirk, Anna Punanova, Youngmin Seo, Yancy Shirley, Adam Ginsburg, Christine Hall, Stella S. R. Offner, Ayushi Singh, Héctor G. Arce, Paola Caselli, Alyssa A. Goodman, Peter G. Martin, Christopher Matzner, Philip C. Myers, Elena Redaelli

Abstract: We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt… ▽ More We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH$_3$ emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH$_3$ is present in over 60\% of lines-of-sight with $A_V \gtrsim 7$ mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH$_3$ is detected toward ~ 40\% of lines-of-sight with $A_V \gtrsim 7$ mag. Moreover, we find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, $σ_v$, and gas kinetic temperature, $T_K$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of $σ_v$ and $T_K$ with increasing star formation activity. The data presented in this paper are publicly available. △ Less

Submitted 20 April, 2017; originally announced April 2017.

Comments: 33 pages, 27 figures, accepted to ApJS. Datasets are publicly available: https://dataverse.harvard.edu/dataverse/GAS_DR1

The JCMT Gould Belt Survey: Evidence for Dust Grain Evolution in Perseus Star-forming Clumps

Authors: Michael Chun-Yuan Chen, J. Di Francesco, D. Johnstone, S. Sadavoy, J. Hatchell, J. C. Mottram, H. Kirk, J. Buckle, D. S. Berry, H. Broekhoven-Fiene, M. J. Currie, M. Fich, T. Jenness, D. Nutter, K. Pattle, J. E. Pineda, C. Quinn, C. Salji, S. Tisi, M. R. Hogerheijde, D. Ward-Thompson, P. Bastien, D. Bresnahan, H. Butner, A. Chrysostomou , et al. (34 additional authors not shown)

Abstract: The dust emissivity spectral index, $β$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $β$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $β$, dust temperature, and optical depth maps of… ▽ More The dust emissivity spectral index, $β$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $β$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $β$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $μ$m, 250 $μ$m, 350 $μ$m, 500 $μ$m, and 850 $μ$m bands. Most of the derived $β$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $β$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $β$ variations within individual clumps and find low $β$ regions correlate with local temperature peaks, hinting at the possible origins of low $β$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape. △ Less

Submitted 19 May, 2016; originally announced May 2016.

Comments: Accepted to ApJ, 15 figures, 3 tables

The dependence of the evolution of SN type Ia progenitors on the C burning rate uncertainty and parameters of convective boundary mixing

Authors: Michael C. Chen, Falk Herwig, Pavel A. Denissenkov, Bill Paxton

Abstract: Evolution of a supernova type Ia progenitor requires formation of a CO white dwarf, which implies a dependence on the C burning rate (CBR). It can also be affected by the recently identified possibility of C flame quenching by convective boundary mixing. We present first results of our study of the combined effect of these two potential sources of uncertainty on the SN Ia progenitor evolution. We… ▽ More Evolution of a supernova type Ia progenitor requires formation of a CO white dwarf, which implies a dependence on the C burning rate (CBR). It can also be affected by the recently identified possibility of C flame quenching by convective boundary mixing. We present first results of our study of the combined effect of these two potential sources of uncertainty on the SN Ia progenitor evolution. We consider the possibility that the CBR is higher than its currently recommended value by as much as a factor of 1000 if unidentified resonances are important, or that it is significantly lower because of the hindrance effect. For stellar models that assume the Schwarzschild boundary for convection, the maximum initial mass for the formation of CO WDs increases from M_i ~ 5.5 Msun for the CBR factor of 1000 to M_i > 7.0 Msun for the CBR factor of 0.01. For C-flame quenching models, hybrid C-O-Ne WDs form for a range of initial mass of Delta M_i ~ 1 Msun, which increases a fraction of stars that form WDs capable of igniting C in a thermonuclear runaway. The most extreme case is found for the CBR factor of 0.1 that is supported by the hindrance model. This nuclear physics assumption, combined with C flame quenching, leads to the formation of a hybrid C-O-Ne WD with a mass of 1.3 Msun. Such WDs do not need to accrete much mass to reach the Chandrasekhar limit. △ Less

Submitted 15 January, 2014; v1 submitted 7 October, 2013; originally announced October 2013.

Comments: 8 pages, 6 figures, accepted by MNRAS