Showing 1–50 of 359 results for author: Wilson, C

Geophysical Observations of the 24 September 2023 OSIRIS-REx Sample Return Capsule Re-Entry

Authors: Elizabeth A. Silber, Daniel C. Bowman, Chris G. Carr, David P. Eisenberg, Brian R. Elbing, Benjamin Fernando, Milton A. Garcés, Robert Haaser, Siddharth Krishnamoorthy, Charles A. Langston, Yasuhiro Nishikawa, Jeremy Webster, Jacob F. Anderson, Stephen Arrowsmith, Sonia Bazargan, Luke Beardslee, Brant Beck, Jordan W. Bishop, Philip Blom, Grant Bracht, David L. Chichester, Anthony Christe, Kenneth Cummins, James Cutts, Lisa Danielson , et al. (57 additional authors not shown)

Abstract: Sample Return Capsules (SRCs) entering Earth's atmosphere at hypervelocity from interplanetary space are a valuable resource for studying meteor phenomena. The 24 September 2023 arrival of the OSIRIS-REx (Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer) SRC provided an unprecedented chance for geophysical observations of a well-characterized source with kn… ▽ More Sample Return Capsules (SRCs) entering Earth's atmosphere at hypervelocity from interplanetary space are a valuable resource for studying meteor phenomena. The 24 September 2023 arrival of the OSIRIS-REx (Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer) SRC provided an unprecedented chance for geophysical observations of a well-characterized source with known parameters, including timing and trajectory. A collaborative effort involving researchers from 16 institutions executed a carefully planned geophysical observational campaign at strategically chosen locations, deploying over 400 ground-based sensors encompassing infrasound, seismic, distributed acoustic sensing (DAS), and GPS technologies. Additionally, balloons equipped with infrasound sensors were launched to capture signals at higher altitudes. This campaign (the largest of its kind so far) yielded a wealth of invaluable data anticipated to fuel scientific inquiry for years to come. The success of the observational campaign is evidenced by the near-universal detection of signals across instruments, both proximal and distal. This paper presents a comprehensive overview of the collective scientific effort, field deployment, and preliminary findings. The early findings have the potential to inform future space missions and terrestrial campaigns, contributing to our understanding of meteoroid interactions with planetary atmospheres. Furthermore, the dataset collected during this campaign will improve entry and propagation models as well as augment the study of atmospheric dynamics and shock phenomena generated by meteoroids and similar sources. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 87 pages, 14 figures

Caustics: A Python Package for Accelerated Strong Gravitational Lensing Simulations

Authors: Connor Stone, Alexandre Adam, Adam Coogan, M. J. Yantovski-Barth, Andreas Filipp, Landung Setiawan, Cordero Core, Ronan Legin, Charles Wilson, Gabriel Missael Barco, Yashar Hezaveh, Laurence Perreault-Levasseur

Abstract: Gravitational lensing is the deflection of light rays due to the gravity of intervening masses. This phenomenon is observed in a variety of scales and configurations, involving any non-uniform mass such as planets, stars, galaxies, clusters of galaxies, and even the large scale structure of the universe. Strong lensing occurs when the distortions are significant and multiple images of the backgrou… ▽ More Gravitational lensing is the deflection of light rays due to the gravity of intervening masses. This phenomenon is observed in a variety of scales and configurations, involving any non-uniform mass such as planets, stars, galaxies, clusters of galaxies, and even the large scale structure of the universe. Strong lensing occurs when the distortions are significant and multiple images of the background source are observed. The lens objects must align on the sky of order ~1 arcsecond for galaxy-galaxy lensing, or 10's of arcseonds for cluster-galaxy lensing. As the discovery of lens systems has grown to the low thousands, these systems have become pivotal for precision measurements and addressing critical questions in astrophysics. Notably, they facilitate the measurement of the Universe's expansion rate, dark matter, supernovae, quasars, and the first stars among other topics. With future surveys expected to discover hundreds of thousands of lensing systems, the modelling and simulation of such systems must occur at orders of magnitude larger scale then ever before. Here we present `caustics`, a Python package designed to handle the extensive computational demands of modeling such a vast number of lensing systems. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 13 pages, 3 figures, submitted to JOSS

Fisher's Mirage: Noise Tightening of Cosmological Constraints in Simulation-Based Inference

Authors: Christopher Wilson, Rachel bean

Abstract: We systematically analyze the implications of statistical noise within numerical derivatives on simulation-based Fisher forecasts for large scale structure surveys. Noisy numerical derivatives resulting from a finite number of simulations, $N_{sims}$, act to bias the associated Fisher forecast such that the resulting marginalized constraints can be significantly tighter than the noise-free limit.… ▽ More We systematically analyze the implications of statistical noise within numerical derivatives on simulation-based Fisher forecasts for large scale structure surveys. Noisy numerical derivatives resulting from a finite number of simulations, $N_{sims}$, act to bias the associated Fisher forecast such that the resulting marginalized constraints can be significantly tighter than the noise-free limit. We show the source of this effect can be traced to the influence of the noise on the marginalization process. Parameters such as the neutrino mass, $\M$, for which higher-order forward differentiation schemes are commonly used, are more prone to noise; the predicted constraints can be akin to those purely from a random instance of statistical noise even using $(1\mathrm{Gpc}/h)^{3}$ simulations with $N_{sims}=500$ realizations. We demonstrate how derivative noise can artificially reduce parameter degeneracies and seemingly null the effects of adding nuisance parameters to the forecast, such as HOD fitting parameters. We mathematically characterize these effects through a full statistical analysis, and demonstrate how confidence intervals for the true noise-free, $N_{sims} \rightarrow \infty$, Fisher constraints can be recovered even when noise comprises a consequential component of the measured signal. The findings and approaches developed here are important for ensuring simulation-based analyses can be accurately used to assess upcoming survey capabilities. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 28 pages, 10 figures

Light Curve Models of Convective Common Envelopes

Authors: N. Noughani, J. Nordhaus, M. Richmond, E. C. Wilson

Abstract: Common envelopes are thought to be the main method for producing tight binaries in the universe as the orbital period shrinks by several orders of magnitude during this phase. Despite their importance for various evolutionary channels, direct detections are rare, and thus observational constraints on common envelope physics are often inferred from post-CE populations. Population constraints sugges… ▽ More Common envelopes are thought to be the main method for producing tight binaries in the universe as the orbital period shrinks by several orders of magnitude during this phase. Despite their importance for various evolutionary channels, direct detections are rare, and thus observational constraints on common envelope physics are often inferred from post-CE populations. Population constraints suggest that the CE phase must be highly inefficient at using orbital energy to drive envelope ejection for low-mass systems and highly efficient for high-mass systems. Such a dichotomy has been explained by an interplay between convection, radiation and orbital decay. If convective transport to the surface occurs faster than the orbit decays, the CE self-regulates and radiatively cools. Once the orbit shrinks such that convective transport is slow compared to orbital decay, a burst occurs as the release of orbital energy can be far in excess of that required to unbind the envelope. With the anticipation of first light for the Rubin Observatory, we calculate light curve models for convective common envelopes and provide the time evolution of apparent magnitudes for the Rubin filters. Convection imparts a distinct signature in the light curves and lengthens the timescales during which they are observable. Given Rubin limiting magnitudes, convective CEs should be detectable out to distances of ~8 Mpc at a rate of ~0.3 per day and provide an intriguing observational test of common envelope physics. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: Submitted to New Astronomy: 14 pages, 6 figures, 6 tables

Cloud-Scale Molecular Gas Properties of the Antennae Merger: A Comparative Study with PHANGS-ALMA Galaxies and NGC 3256

Authors: Nathan Brunetti, Christine D. Wilson, Hao He, Jiayi Sun, Adam K. Leroy, Erik Rosolowsky, Ashley Bemis, Frank Bigiel, Brent Groves, Toshiki Saito, Eva Schinnerer

Abstract: We present observations of the central 9 kpc of the Antennae merger (NGC 4038/9) at 55 pc resolution in the CO 2-1 line obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). We use a pixel-based analysis to compare the gas properties in the Antennae to those in 70 nearby spiral galaxies from the PHANGS-ALMA survey, as well as the merger and nearest luminous infrared galaxy NGC 325… ▽ More We present observations of the central 9 kpc of the Antennae merger (NGC 4038/9) at 55 pc resolution in the CO 2-1 line obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). We use a pixel-based analysis to compare the gas properties in the Antennae to those in 70 nearby spiral galaxies from the PHANGS-ALMA survey, as well as the merger and nearest luminous infrared galaxy NGC 3256. Compared to PHANGS galaxies at matched spatial resolution, the molecular gas in the Antennae exhibits some of the highest surface densities, velocity dispersions, peak brightness temperatures, and turbulent pressures. However, the virial parameters in the Antennae are consistent with many of the PHANGS galaxies. NGC 3256 has similar gas surface densities but higher nuclear velocity dispersions than the Antennae, as well as higher system-wide peak brightness temperatures and virial parameters. NGC 3256 is at a later stage in the merging process than the Antennae, which may result in more intense merger-driven gas flows that could drive up the turbulence in the gas. The high virial parameters in NGC 3256 may indicate that this increased turbulence is suppressing future star formation as NGC 3256 moves out of the starburst phase. In comparison, the relatively normal virial parameters in the Antennae may imply that it is about to undergo a new burst of star formation. △ Less

Submitted 6 April, 2024; originally announced April 2024.

Comments: 16 pages, 8 figures, accepted to MNRAS

Unraveling the Mystery of the Low CO-to-H$_2$ Conversion Factor in Starburst Galaxies: RADEX Modeling of the Antennae

Authors: Hao He, Christine D. Wilson, Jiayi Sun, Yu-Hsuan Teng, Erik Rosolowsky, Ashley R. Bemis

Abstract: CO emission has been widely used as a tracer of molecular gas mass. However, it is a long-standing issue to accurately constrain the CO-to-H$_2$ conversion factor ($α_{\mathrm{CO}}$) that converts CO luminosity to molecular gas mass, especially in starburst galaxies. We present the first resolved $α_{\mathrm{CO}}$ modeling results with multiple ALMA CO and $^{13}$CO transition observations at both… ▽ More CO emission has been widely used as a tracer of molecular gas mass. However, it is a long-standing issue to accurately constrain the CO-to-H$_2$ conversion factor ($α_{\mathrm{CO}}$) that converts CO luminosity to molecular gas mass, especially in starburst galaxies. We present the first resolved $α_{\mathrm{CO}}$ modeling results with multiple ALMA CO and $^{13}$CO transition observations at both giant molecular cloud (GMC) scale at 150 pc and kpc scale for one of the closest starburst mergers, the Antennae. By combining our CO modeling results and measurements of 350 GHz dust continuum, we find that most GMCs in the Antennae have $α_{\mathrm{CO}}$ values $\sim$4 times smaller than the commonly adopted Milky Way value (4.3). We find $α_{\mathrm{CO}}$ at GMC scales shows a strong dependence on CO intensity, $^{13}$CO/CO ratio and GMC velocity dispersion, which is consistent with various theoretical and simulation predictions. Specifically, we suggest that the $^{13}$CO/CO line ratio and the velocity dispersion can be used to infer $α_{\mathrm{CO}}$ in starburst regions. By applying our modeled $α_{\mathrm{CO}}$ in GMC analyses, we find that GMCs in the Antennae are less gravitationally bound than in normal spiral galaxies, which is more consistent with what is predicted by merger simulations. At kpc scale, we find that our modeled $α_{\mathrm{CO}}$ values are smaller than the modeled $α_{\mathrm{CO}}$ at GMC scale by 40%, which can be due to inclusion of a diffuse gas component with lower $α_{\mathrm{CO}}$ values. We find a similar correlation of $α_{\mathrm{CO}}$ and CO intensity at kpc scales to that at GMC scales. △ Less

Submitted 9 June, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: 22 pages and 16 figures in the main text; accepted to ApJ

Hidden Gems on a Ring: Infant Massive Clusters and Their Formation Timeline Unveiled by ALMA, HST, and JWST in NGC 3351

Authors: Jiayi Sun, Hao He, Kyle Batschkun, Rebecca C. Levy, Kimberly Emig, M. Jimena Rodriguez, Hamid Hassani, Adam K. Leroy, Eva Schinnerer, Eve C. Ostriker, Christine D. Wilson, Alberto D. Bolatto, Elisabeth A. C. Mills, Erik Rosolowsky, Janice C. Lee, Daniel A. Dale, Kirsten L. Larson, David A. Thilker, Leonardo Ubeda, Bradley C. Whitmore, Thomas G. Williams, Ashley. T. Barnes, Frank Bigiel, Melanie Chevance, Simon C. O. Glover , et al. (16 additional authors not shown)

Abstract: We study young massive clusters (YMCs) in their embedded "infant" phase with $\sim0.\!^{\prime\prime}1$ ALMA, HST, and JWST observations targeting the central starburst ring in NGC 3351, a nearby Milky Way analog galaxy. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on… ▽ More We study young massive clusters (YMCs) in their embedded "infant" phase with $\sim0.\!^{\prime\prime}1$ ALMA, HST, and JWST observations targeting the central starburst ring in NGC 3351, a nearby Milky Way analog galaxy. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on the ALMA continuum and molecular line data, as well as ancillary measurements for the HST and JWST counterparts, we identify 14 sources as infant star clusters with high stellar and/or gas masses (${\sim}10^5\;\mathrm{M_\odot}$), small radii (${\lesssim}\,5\;\mathrm{pc}$), large escape velocities ($6{-}10\;\mathrm{km/s}$), and short free-fall times ($0.5{-}1\;\mathrm{Myr}$). Their multiwavelength properties motivate us to divide them into four categories, likely corresponding to four evolutionary stages from starless clumps to exposed HII region-cluster complexes. Leveraging age estimates for HST-identified clusters in the same region, we infer an evolutionary timeline going from $\sim$1-2 Myr before cluster formation as starless clumps, to $\sim$4-6 Myr after as exposed HII region-cluster complexes. Finally, we show that the YMCs make up a substantial fraction of recent star formation across the ring, exhibit an non-uniform azimuthal distribution without a very coherent evolutionary trend along the ring, and are capable of driving large-scale gas outflows. △ Less

Submitted 10 April, 2024; v1 submitted 25 January, 2024; originally announced January 2024.

Comments: 27 pages, 12 figures; ApJ accepted

SDSS-IV from 2014 to 2016: A Detailed Demographic Comparison over Three Years

Authors: Amy M. Jones, Rachael L. Beaton, Brian A. Cherinka, Karen L. Masters, Sara Lucatello, Aleksandar M. Diamond-Stanic, Sarah A. Bird, Michael R. Blanton, Katia Cunha, Emily E. Farr, Diane Feuillet, Peter M. Frinchaboy, Alex Hagen, Karen Kinemuchi, Britt Lundgren, Mariarosa L. Marinelli, Adam D. Myers, Alexandre Roman-Lopes, Ashley J. Ross, Jose R. Sanchez-Gallego, Sarah J. Schmidt, Jennifer Sobeck, Keivan G. Stassun, Jamie Tayar, Mariana Vargas-Magana , et al. (2 additional authors not shown)

Abstract: The Sloan Digital Sky Survey (SDSS) is one of the largest international astronomy organizations. We present demographic data based on surveys of its members from 2014, 2015 and 2016, during the fourth phase of SDSS (SDSS-IV). We find about half of SDSS-IV collaboration members were based in North America, a quarter in Europe, and the remainder in Asia and Central and South America. Overall, 26-36%… ▽ More The Sloan Digital Sky Survey (SDSS) is one of the largest international astronomy organizations. We present demographic data based on surveys of its members from 2014, 2015 and 2016, during the fourth phase of SDSS (SDSS-IV). We find about half of SDSS-IV collaboration members were based in North America, a quarter in Europe, and the remainder in Asia and Central and South America. Overall, 26-36% are women (from 2014 to 2016), up to 2% report non-binary genders. 11-14% report that they are racial or ethnic minorities where they live. The fraction of women drops with seniority, and is also lower among collaboration leadership. Men in SDSS-IV were more likely to report being in a leadership role, and for the role to be funded and formally recognized. SDSS-IV collaboration members are twice as likely to have a parent with a college degree, than the general population, and are ten times more likely to have a parent with a PhD. This trend is slightly enhanced for female collaboration members. Despite this, the fraction of first generation college students (FGCS) is significant (31%). This fraction increased among collaboration members who are racial or ethnic minorities (40-50%), and decreased among women (15-25%). SDSS-IV implemented many inclusive policies and established a dedicated committee, the Committee on INclusiveness in SDSS (COINS). More than 60% of the collaboration agree that the collaboration is inclusive; however, collaboration leadership more strongly agree with this than the general membership. In this paper, we explain these results in full, including the history of inclusive efforts in SDSS-IV. We conclude with a list of suggested recommendations based on our findings, which can be used to improve equity and inclusion in large astronomical collaborations, which we argue is not only moral, but will also optimize their scientific output. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 30 pages, 9 figures, accepted in PASP

Stored in the archives: Uncovering the CN/CO intensity ratio with ALMA in nearby U/LIRGs

Authors: Blake Ledger, Toshiki Saito, Daisuke Iono, Christine D. Wilson

Abstract: We present an archival Atacama Large Millimeter/submillimeter Array (ALMA) study of the CN N = 1 - 0 / CO J = 1 - 0 intensity ratio in nearby (z < 0.05) Ultra Luminous and Luminous Infrared Galaxies (U/LIRGs). We identify sixteen U/LIRGs that have been observed in both CN and CO lines at $\sim$ 500 pc resolution based on sixteen different ALMA projects. We measure the (CN bright)/CO and (CN bright… ▽ More We present an archival Atacama Large Millimeter/submillimeter Array (ALMA) study of the CN N = 1 - 0 / CO J = 1 - 0 intensity ratio in nearby (z < 0.05) Ultra Luminous and Luminous Infrared Galaxies (U/LIRGs). We identify sixteen U/LIRGs that have been observed in both CN and CO lines at $\sim$ 500 pc resolution based on sixteen different ALMA projects. We measure the (CN bright)/CO and (CN bright)/(CN faint) intensity ratios at an ensemble of molecular clouds scales (CN bright = CN N = 1 - 0, J = 3/2 - 1/2; CN faint = CN N = 1 - 0, J = 1/2 - 1/2 hyperfine groupings). Our global measured (CN bright)/CO ratios range from 0.02-0.15 in LIRGs and 0.08-0.17 in ULIRGs. We attribute the larger spread in LIRGs to the variety of galaxy environments included in our sample. Overall, we find that the (CN bright)/CO ratio is higher in nuclear regions, where the physical and excitation conditions favour increased CN emission relative to the disk regions. 10 out of 11 galaxies which contain well-documented active galactic nuclei show higher ratios in the nucleus compared to the disk. Finally, we measure the median resolved (CN bright)/(CN faint) ratio and use it to estimate the total integrated CN line optical depth in ULIRGs ($τ\sim$ 0.96) and LIRGs ($τ\sim$ 0.23). The optical depth difference is likely due to the higher molecular gas surface densities found in the more compact ULIRG systems. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: Accepted to MNRAS; 18 pages, 9 figures

VERTICO and IllustrisTNG: The spatially resolved effects of environment on galactic gas

Authors: Adam R. H. Stevens, Toby Brown, Benedikt Diemer, Annalisa Pillepich, Lars Hernquist, Dylan Nelson, Yannick M. Bahé, Alessandro Boselli, Timothy A. Davis, Pascal J. Elahi, Sara L. Ellison, María J. Jiménez-Donaire, Ian D. Roberts, Kristine Spekkens, Vicente Villanueva, Adam B. Watts, Christine D. Wilson, Nikki Zabel

Abstract: It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies w… ▽ More It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies with $m_* > \! 10^9\,{\rm M}_\odot$ and ${\rm SFR} \! > 0.05\,{\rm M}_\odot\,{\rm yr}^{-1}$. We select galaxies in the simulation that are analogous to those in the HERACLES and VERTICO surveys, and mock-observe them to the common specifications of the data. Although TNG50 does not quantitatively match the observed gas surface densities in the centers of galaxies, the simulation does qualitatively reproduce the trends of gas truncation and central density suppression seen in VERTICO in both HI and H$_2$. This result promises that modern cosmological hydrodynamic simulations can be used to reliably model the post-infall histories of cluster satellite galaxies. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 9 pages, 4 figures, accepted in ApJL

Exploring the Origin of the Distance Bimodality of Stars in the Periphery of the Small Magellanic Cloud with APOGEE and Gaia

Authors: Andres Almeida, Steven R. Majewski, David L. Nidever, Knut A. G. Olsen, Antonela Monachesi, Nitya Kallivayalil, Sten Hasselquist, Yumi Choi, Joshua T. Povick, John C. Wilson, Doug Geisler, Richard R. Lane, Christian Nitschelm, Jennifer S. Sobeck, Guy S. Stringfellow

Abstract: The Magellanic Cloud system represents a unique laboratory for study of both interacting dwarf galaxies and the ongoing process of the formation of the Milky Way and its halo. We focus on one aspect of this complex, 3 body interaction - the dynamical perturbation of the Small Magellanic Cloud (SMC) by the Large Magellanic Cloud (LMC), and specifically potential tidal effects on the SMC's eastern s… ▽ More The Magellanic Cloud system represents a unique laboratory for study of both interacting dwarf galaxies and the ongoing process of the formation of the Milky Way and its halo. We focus on one aspect of this complex, 3 body interaction - the dynamical perturbation of the Small Magellanic Cloud (SMC) by the Large Magellanic Cloud (LMC), and specifically potential tidal effects on the SMC's eastern side. Using Gaia astrometry and the precise radial velocities and multi-element chemical abundances from APOGEE-2 DR17, we explore the well-known distance bimodality on the eastern side of the SMC. Through estimated stellar distances, proper motions, and radial velocities, we characterize the kinematics of the two populations in the bimodality and compare their properties with those of SMC populations elsewhere. Moreover, while all regions explored by APOGEE seem to show a single chemical enrichment history, the metallicity distribution function (MDF), of the "far" stars on the eastern periphery of the SMC is found to resemble that for the more metal-poor fields of the western periphery, whereas the MDF for the "near" stars on the eastern periphery resembles that for stars in the SMC center. The closer eastern periphery stars also show radial velocities (corrected for SMC rotation and bulk motion) that are, on average, approaching us relative to all other SMC populations sampled. We interpret these trends as evidence that the near stars on the eastern side of the SMC represent material pulled out of the central SMC as part of its tidal interaction with the LMC. △ Less

Submitted 25 August, 2023; originally announced August 2023.

VERTICO VII: Environmental quenching caused by suppression of molecular gas content and star formation efficiency in Virgo Cluster galaxies

Authors: Toby Brown, Ian D. Roberts, Mallory Thorp, Sara L. Ellison, Nikki Zabel, Christine D. Wilson, Yannick M. Bahé, Dhruv Bisaria, Alberto D. Bolatto, Alessandro Boselli, Aeree Chung, Luca Cortese, Barbara Catinella, Timothy A. Davis, María J. Jiménez-Donaire, Claudia D. P. Lagos, Bumhyun Lee, Laura C. Parker, Rory Smith, Kristine Spekkens, Adam R. H. Stevens, Vicente Villanueva, Adam B. Watts

Abstract: We study how environment regulates the star formation cycle of 33 Virgo Cluster satellite galaxies on 720 parsec scales. We present the first resolved star-forming main sequence for cluster galaxies, dividing the sample based on their global HI properties and comparing to a control sample of field galaxies. HI-poor cluster galaxies have reduced star formation rate (SFR) surface densities with resp… ▽ More We study how environment regulates the star formation cycle of 33 Virgo Cluster satellite galaxies on 720 parsec scales. We present the first resolved star-forming main sequence for cluster galaxies, dividing the sample based on their global HI properties and comparing to a control sample of field galaxies. HI-poor cluster galaxies have reduced star formation rate (SFR) surface densities with respect to both HI-normal cluster and field galaxies (0.5 dex), suggesting that mechanisms regulating the global HI content are responsible for quenching local star formation. We demonstrate that the observed quenching in HI-poor galaxies is caused by environmental processes such as ram pressure stripping (RPS) simultaneously reducing molecular gas surface density and star formation efficiency (SFE), compared to regions in HI-normal systems (by 0.38 and 0.22 dex, respectively). We observe systematically elevated SFRs that are driven by increased molecular gas surface densities at fixed stellar mass surface density in the outskirts of early-stage RPS galaxies, while SFE remains unchanged with respect to the field sample. We quantify how RPS and starvation affect the star formation cycle of inner and outer galaxy discs as they are processed by the cluster. We show both are effective quenching mechanisms with the key difference being that RPS acts upon the galaxy outskirts while starvation regulates the star formation cycle throughout disc, including within the truncation radius. For both processes, the quenching is caused by a simultaneous reduction in molecular gas surface densities and SFE at fixed stellar mass surface density. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 17 pages, 1 table, 5 figures, accepted for publication in ApJ

A Tale of Two Disks: Mapping the Milky Way with the Final Data Release of APOGEE

Authors: Julie Imig, Cathryn Price, Jon A. Holtzman, Alexander Stone-Martinez, Steven R. Majewski, David H. Weinberg, Jennifer A. Johnson, Carlos Allende Prieto, Rachael L. Beaton, Timothy C. Beers, Dmitry Bizyaev, Michael R. Blanton, Joel R. Brownstein, Katia Cunha, José G. Fernández-Trincado, Diane K. Feuillet, Sten Hasselquist, Christian R. Hayes, Henrik Jönsson, Richard R. Lane, Jianhui Lian, Szabolcs Mészáros, David L. Nidever, Annie C. Robin, Matthew Shetrone , et al. (2 additional authors not shown)

Abstract: We present new maps of the Milky Way disk showing the distribution of metallicity ([Fe/H]), $α$-element abundances ([Mg/Fe]), and stellar age, using a sample of 66,496 red giant stars from the final data release (DR17) of the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We measure radial and vertical gradients, quantify the distribution functions for age and metallicity,… ▽ More We present new maps of the Milky Way disk showing the distribution of metallicity ([Fe/H]), $α$-element abundances ([Mg/Fe]), and stellar age, using a sample of 66,496 red giant stars from the final data release (DR17) of the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We measure radial and vertical gradients, quantify the distribution functions for age and metallicity, and explore chemical clock relations across the Milky Way for the low-$α$ disk, high-$α$ disk, and total population independently. The low-$α$ disk exhibits a negative radial metallicity gradient of $-0.06 \pm 0.001$ dex kpc$^{-1}$, which flattens with distance from the midplane. The high-$α$ disk shows a flat radial gradient in metallicity and age across nearly all locations of the disk. The age and metallicity distribution functions shift from negatively skewed in the inner Galaxy to positively skewed at large radius. Significant bimodality in the [Mg/Fe]-[Fe/H] plane and in the [Mg/Fe]-age relation persist across the entire disk. The age estimates have typical uncertainties of $\sim0.15$ in $\log$(age) and may be subject to additional systematic errors, which impose limitations on conclusions drawn from this sample. Nevertheless, these results act as critical constraints on galactic evolution models, constraining which physical processes played a dominant role in the formation of the Milky Way disk. We discuss how radial migration predicts many of the observed trends near the solar neighborhood and in the outer disk, but an additional more dramatic evolution history, such as the multi-infall model or a merger event, is needed to explain the chemical and age bimodality elsewhere in the Galaxy. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 41 pages, 32 figures, accepted to ApJ

Journal ref: ApJ 954 124 (2023)

VERTICO VI: Cold-gas asymmetries in Virgo cluster galaxies

Authors: Ian D. Roberts, Toby Brown, Nikki Zabel, Christine D. Wilson, Aeree Chung, Laura C. Parker, Dhruv Bisaria, Alessandro Boselli, Barbara Catinella, Ryan Chown, Luca Cortese, Timothy A. Davis, Sara Ellison, Maria Jesus Jimenez-Donaire, Bumhyun Lee, Rory Smith, Kristine Spekkens, Adam R. H. Stevens, Mallory Thorp, Vincente Villanueva, Adam B. Watts, Charlotte Welker, Hyein Yoon

Abstract: We analyze cold-gas distributions in Virgo cluster galaxies using resolved CO(2-1) (tracing molecular hydrogen, H2) and HI observations from the Virgo Environment Traced In CO (VERTICO) and the VLA Imaging of Virgo in Atomic Gas (VIVA) surveys. From a theoretical perspective, it is expected that environmental processes in clusters will have a stronger influence on diffuse atomic gas compared to th… ▽ More We analyze cold-gas distributions in Virgo cluster galaxies using resolved CO(2-1) (tracing molecular hydrogen, H2) and HI observations from the Virgo Environment Traced In CO (VERTICO) and the VLA Imaging of Virgo in Atomic Gas (VIVA) surveys. From a theoretical perspective, it is expected that environmental processes in clusters will have a stronger influence on diffuse atomic gas compared to the relatively dense molecular gas component, and that these environmental perturbations can compress the cold interstellar medium in cluster galaxies leading to elevated star formation. In this work we observationally test these predictions for star-forming satellite galaxies within the Virgo cluster. We divide our Virgo galaxy sample into HI-normal, HI-tailed, and HI-truncated classes and show, unsurprisingly, that the HI-tailed galaxies have the largest quantitative HI asymmetries. We also compare to a control sample of non-cluster galaxies and find that Virgo galaxies, on average, have HI asymmetries that are 40 +/- 10 per cent larger than the control. There is less separation between control, HI-normal, HI-tailed, and HI-truncated galaxies in terms of H2 asymmetries, and on average, Virgo galaxies have H2 asymmetries that are only marginally (20 +/- 10 per cent) larger than the control sample. We find a weak correlation between HI and H2 asymmetries over our entire sample, but a stronger correlation for those specific galaxies being strongly impacted by environmental perturbations. Finally, we divide the discs of the HI-tailed Virgo galaxies into a leading half and trailing half according to the observed tail direction. We find evidence for excess molecular gas mass on the leading halves of the disc. This excess molecular gas on the leading half is accompanied by an excess in star formation rate such that the depletion time is, on average, unchanged. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 15 pages, 8 figures, 1 table, accepted for publication in A&A

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

Coronal Heating as Determined by the Solar Flare Frequency Distribution Obtained by Aggregating Case Studies

Authors: James Paul Mason, Alexandra Werth, Colin G. West, Allison A. Youngblood, Donald L. Woodraska, Courtney Peck, Kevin Lacjak, Florian G. Frick, Moutamen Gabir, Reema A. Alsinan, Thomas Jacobsen, Mohammad Alrubaie, Kayla M. Chizmar, Benjamin P. Lau, Lizbeth Montoya Dominguez, David Price, Dylan R. Butler, Connor J. Biron, Nikita Feoktistov, Kai Dewey, N. E. Loomis, Michal Bodzianowski, Connor Kuybus, Henry Dietrick, Aubrey M. Wolfe , et al. (977 additional authors not shown)

Abstract: Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms th… ▽ More Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfvén waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, $α=2$ as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed $>$600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: pre-flare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that $α= 1.63 \pm 0.03$. This is below the critical threshold, suggesting that Alfvén waves are an important driver of coronal heating. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 1,002 authors, 14 pages, 4 figures, 3 tables, published by The Astrophysical Journal on 2023-05-09, volume 948, page 71

VERTICO V: The environmentally driven evolution of the inner cold gas discs of Virgo cluster galaxies

Authors: Adam B. Watts, Luca Cortese, Barbara Catinella, Toby Brown, Christine D. Wilson, Nikki Zabel, Ian D. Roberts, Timothy A. Davis, Mallory Thorp, Aeree Chung, Adam R. H. Stevens, Sara L. Ellison, Kristine Spekkens, Laura C. Parker, Yannick M. Bahé, Vicente Villanueva, María Jiménez-Donaire, Dhruv Bisaria, Alessandro Boselli, Alberto D. Bolatto, Bumhyun Lee

Abstract: The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virg… ▽ More The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virgo Environment traced in CO survey (VERTICO) and archival HI data, we study the HI and molecular gas within the optical discs of Virgo cluster galaxies on 1.2-kpc scales with spatially resolved scaling relations between stellar ($Σ_{\star}$), HI ($Σ_\mathrm{HI}$), and molecular gas ($Σ_\mathrm{mol}$) surface densities. Adopting HI deficiency as a measure of environmental impact, we find evidence that, in addition to removing the HI at large radii, the cluster processes also lower the average $Σ_\mathrm{HI}$ of the remaining gas even in the central 1.2 kpc. The impact on molecular gas is comparatively weaker than on the HI, and we show that the lower $Σ_\mathrm{mol}$ gas is removed first. In the most HI-deficient galaxies, however, we find evidence that environmental processes reduce the typical $Σ_\mathrm{mol}$ of the remaining gas by nearly a factor of 3. We find no evidence for environment-driven elevation of $Σ_\mathrm{HI}$ or $Σ_\mathrm{mol}$ in HI-deficient galaxies. Using the ratio of $Σ_\mathrm{mol}$-to-$Σ_\mathrm{HI}$ in individual regions, we show that changes in the ISM physical conditions, estimated using the total gas surface density and midplane hydrostatic pressure, cannot explain the observed reduction in molecular gas content. Instead, we suggest that direct stripping of the molecular gas is required to explain our results. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: Accepted for publication in PASA. Main text is 19 pages including 12 figures and 3 tables, plus 1 appendix. A 2.5 min, high-level summary can be found at https://youtu.be/7djMmVEpDVc

A nearly constant CN/HCN line ratio in nearby galaxies: CN as a new tracer of dense gas

Authors: Christine D. Wilson, Ashley Bemis, Blake Ledger, Osvald Klimi

Abstract: We investigate the relationship between CN N = 1 - 0 and HCN J = 1 - 0 emission on scales from 30 pc to 400 pc using ALMA archival data, for which CN is often observed simultaneously with the CO J = 1 - 0 line. In a sample of 9 nearby galaxies ranging from ultra-luminous infrared galaxies to normal spiral galaxies, we measure a remarkably constant CN/HCN line intensity ratio of 0.86 $\pm$ 0.07 (st… ▽ More We investigate the relationship between CN N = 1 - 0 and HCN J = 1 - 0 emission on scales from 30 pc to 400 pc using ALMA archival data, for which CN is often observed simultaneously with the CO J = 1 - 0 line. In a sample of 9 nearby galaxies ranging from ultra-luminous infrared galaxies to normal spiral galaxies, we measure a remarkably constant CN/HCN line intensity ratio of 0.86 $\pm$ 0.07 (standard deviation of 0.20). This relatively constant CN/HCN line ratio is rather unexpected, as models of photon dominated regions have suggested that HCN emission traces shielded regions with high column densities while CN should trace dense gas exposed to high ultraviolet radiation fields. We find that the CN/HCN line ratio shows no significant correlation with molecular gas surface density, but shows a mild trend (increase of ~ 1.3 per dex) with both star formation rate surface density and star formation efficiency (the inverse of the molecular gas depletion time). Some starburst and active galactic nuclei show small enhancements in their CN/HCN ratio, while other nuclei show no significant difference from their surrounding disks. The nearly constant CN/HCN line ratio implies that CN, like HCN, can be used as a tracer of dense gas mass and dense gas fraction in nearby galaxies. △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: Accepted to MNRAS; 21 pages, 12 figures

The JCMT Nearby Galaxies Legacy Survey: SCUBA-2 observations of nearby galaxies

Authors: Kate Pattle, Walter Gear, Christine D. Wilson

Abstract: We present 850$μ$m observations of a sample of 8 nearby spiral galaxies, made using the SCUBA-2 camera on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT Nearby Galaxies Legacy Survey (NGLS). We corrected our data for the presence of the $^{12}$CO $J=3\to 2$ line in the SCUBA-2 850$μ$m bandwidth using NGLS HARP data, finding a typical $^{12}$CO contribution of $\sim 20$%. We measured… ▽ More We present 850$μ$m observations of a sample of 8 nearby spiral galaxies, made using the SCUBA-2 camera on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT Nearby Galaxies Legacy Survey (NGLS). We corrected our data for the presence of the $^{12}$CO $J=3\to 2$ line in the SCUBA-2 850$μ$m bandwidth using NGLS HARP data, finding a typical $^{12}$CO contribution of $\sim 20$%. We measured dust column densities, temperatures and opacity indices by fitting spectral energy distributions constructed from SCUBA-2 and archival Herschel observations, and used archival GALEX and Spitzer data to make maps of surface density of star formation ($Σ_{\rm SFR}$). Typically, comparing SCUBA-2-derived H$_2$ surface densities ($Σ_{\rm H_2}$) to $Σ_{\rm SFR}$ gives shallow star formation law indices within galaxies, with SCUBA-2-derived values typically being sublinear and Herschel-derived values typically being broadly linear. This difference is likely due to the effects of atmospheric filtering on the SCUBA-2 data. Comparing the mean values of $Σ_{\rm H_2}$ and $Σ_{\rm SFR}$ of the galaxies in our sample returns a steeper star formation law index, broadly consistent with both the Kennicutt-Schmidt value of 1.4 and linearity. Our results show that a SCUBA-2 detection is a good predictor of star formation. We suggest that Herschel emission traces gas in regions which will form stars on timescales $\sim 5-100$ Myr, comparable to the star formation timescale traced by GALEX and Spitzer data, while SCUBA-2 preferentially traces the densest gas within these regions, which likely forms stars on shorter timescales. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: Accepted for publication in MNRAS. 15 pages, 12 figures, 8 tables (plus 15 pages of appendices, with 31 figures)

Star Formation Laws and Efficiencies across 80 Nearby Galaxies

Authors: Jiayi Sun, Adam K. Leroy, Eve C. Ostriker, Sharon Meidt, Erik Rosolowsky, Eva Schinnerer, Christine D. Wilson, Dyas Utomo, Francesco Belfiore, Guillermo A. Blanc, Eric Emsellem, Christopher Faesi, Brent Groves, Annie Hughes, Eric W. Koch, Kathryn Kreckel, Daizhong Liu, Hsi-An Pan, Jerome Pety, Miguel Querejeta, Alessandro Razza, Toshiki Saito, Amy Sardone, Antonio Usero, Thomas G. Williams , et al. (15 additional authors not shown)

Abstract: We measure empirical relationships between the local star formation rate (SFR) and properties of the star-forming molecular gas on 1.5 kpc scales across 80 nearby galaxies. These relationships, commonly referred to as "star formation laws," aim at predicting the local SFR surface density from various combinations of molecular gas surface density, galactic orbital time, molecular cloud free-fall ti… ▽ More We measure empirical relationships between the local star formation rate (SFR) and properties of the star-forming molecular gas on 1.5 kpc scales across 80 nearby galaxies. These relationships, commonly referred to as "star formation laws," aim at predicting the local SFR surface density from various combinations of molecular gas surface density, galactic orbital time, molecular cloud free-fall time, and the interstellar medium dynamical equilibrium pressure. Leveraging a multiwavelength database built for the PHANGS survey, we measure these quantities consistently across all galaxies and quantify systematic uncertainties stemming from choices of SFR calibrations and the CO-to-H$_2$ conversion factors. The star formation laws we examine show 0.3-0.4 dex of intrinsic scatter, among which the molecular Kennicutt-Schmidt relation shows a $\sim$10% larger scatter than the other three. The slope of this relation ranges $β\approx0.9{-}1.2$, implying that the molecular gas depletion time remains roughly constant across the environments probed in our sample. The other relations have shallower slopes ($β\approx0.6{-}1.0$), suggesting that the star formation efficiency (SFE) per orbital time, the SFE per free-fall time, and the pressure-to-SFR surface density ratio (i.e., the feedback yield) may vary systematically with local molecular gas and SFR surface densities. Last but not least, the shapes of the star formation laws depend sensitively on methodological choices. Different choices of SFR calibrations can introduce systematic uncertainties of at least 10-15% in the star formation law slopes and 0.15-0.25 dex in their normalization, while the CO-to-H$_2$ conversion factors can additionally produce uncertainties of 20-25% for the slope and 0.10-0.20 dex for the normalization. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 10 pages main text + 2 appendices. ApJL in press. Data products available at https://www.canfar.net/storage/list/phangs/RELEASES/Sun_etal_2023 . Slides summarizing key results can be found at https://www.dropbox.com/s/5gsegexeo9n0t05/Sun_et_PHANGS_2023.pptx?dl=0

Journal ref: ApJL, 945, L19 (2023)

Characterization of low-mass companions to $\textit{Kepler}$ objects of interest observed with APOGEE-N

Authors: Caleb I. Cañas, Chad F. Bender, Suvrath Mahadevan, Dmitry Bizyaev, Nathan De Lee, Scott W. Fleming, Fred Hearty, Steven R. Majewski, Christian Nitschelm, Donald P. Schneider, Javier Serna, Keivan G. Stassun, Guðmundur Stefánsson, Guy S. Stringfellow, John C. Wilson

Abstract: We report the characterization of 28 low-mass ($0.02\mathrm{~M_\odot}\le\mathrm{~M_{2}}\le0.25\mathrm{~M_\odot}$) companions to $\textit{Kepler}$ objects of interest (KOIs), eight of which were previously designated confirmed planets. These objects were detected as transiting companions to Sun-like stars (G and F dwarfs) by the $\textit{Kepler}$ mission and are confirmed as single-lined spectrosco… ▽ More We report the characterization of 28 low-mass ($0.02\mathrm{~M_\odot}\le\mathrm{~M_{2}}\le0.25\mathrm{~M_\odot}$) companions to $\textit{Kepler}$ objects of interest (KOIs), eight of which were previously designated confirmed planets. These objects were detected as transiting companions to Sun-like stars (G and F dwarfs) by the $\textit{Kepler}$ mission and are confirmed as single-lined spectroscopic binaries in the current work using the northern multiplexed Apache Point Observatory Galactic Evolution Experiment near-infrared spectrograph (APOGEE-N) as part of the third and fourth Sloan Digital Sky Surveys. We have observed hundreds of KOIs using APOGEE-N and collected a total of 43,175 spectra with a median of 19 visits and a median baseline of $\sim1.9$ years per target. We jointly model the $\textit{Kepler}$ photometry and APOGEE-N radial velocities to derive fundamental parameters for this subset of 28 transiting companions. The radii for most of these low-mass companions are over-inflated (by $\sim10\%$) when compared to theoretical models. Tidally locked M dwarfs on short period orbits show the largest amount of inflation, but inflation is also evident for companions that are well separated from the host star. We demonstrate that APOGEE-N data provides reliable radial velocities when compared to precise high-resolution spectrographs that enable detailed characterization of individual systems and the inference of orbital elements for faint ($H>12$) KOIs. The data from the entire APOGEE-KOI program is public and presents an opportunity to characterize an extensive subset of the binary population observed by $\textit{Kepler}$. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 98 pages (include 56 for the figure sets), 10 tables, 7 figures, 2 figure sets, accepted for publication in ApJS

Journal ref: ApJS (2023), 265, 50

Molecular gas and star formation in nearby starburst galaxy mergers

Authors: Hao He, Connor Bottrell, Christine Wilson, Jorge Moreno, Blakesley Burkhart, Christopher C. Hayward, Lars Hernquist, Angela Twum

Abstract: We employ the Feedback In Realistic Environments (FIRE-2) physics model to study how the properties of giant molecular clouds (GMCs) evolve during galaxy mergers. We conduct a pixel-by-pixel analysis of molecular gas properties in both the simulated control galaxies and galaxy major mergers. The simulated GMC-pixels in the control galaxies follow a similar trend in a diagram of velocity dispersion… ▽ More We employ the Feedback In Realistic Environments (FIRE-2) physics model to study how the properties of giant molecular clouds (GMCs) evolve during galaxy mergers. We conduct a pixel-by-pixel analysis of molecular gas properties in both the simulated control galaxies and galaxy major mergers. The simulated GMC-pixels in the control galaxies follow a similar trend in a diagram of velocity dispersion ($σ_v$) versus gas surface density ($Σ_{\mathrm{mol}}$) to the one observed in local spiral galaxies in the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) survey. For GMC-pixels in simulated mergers, we see a significant increase of factor of 5 - 10 in both $Σ_{\mathrm{mol}}$ and $σ_v$, which puts these pixels above the trend of PHANGS galaxies in the $σ_v$ vs $Σ_{\mathrm{mol}}$ diagram. This deviation may indicate that GMCs in the simulated mergers are much less gravitationally bound compared with simulated control galaxies with virial parameter ($α_{\mathrm{vir}}$) reaching 10 - 100. Furthermore, we find that the increase in $α_{\mathrm{vir}}$ happens at the same time as the increase in global star formation rate (SFR), which suggests stellar feedback is responsible for dispersing the gas. We also find that the gas depletion time is significantly lower for high $α_{\mathrm{vir}}$ GMCs during a starburst event. This is in contrast to the simple physical picture that low $α_{\mathrm{vir}}$ GMCs are easier to collapse and form stars on shorter depletion times. This might suggest that some other physical mechanisms besides self-gravity are helping the GMCs in starbursting mergers collapse and form stars. △ Less

Submitted 31 March, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: 22 pages, 11 figures. Accepted to ApJ. Link to animation updated

Does the HCN/CO ratio trace the star-forming fraction of gas? I. A comparison with analytical models of star formation

Authors: Ashley R. Bemis, Christine D. Wilson

Abstract: We use archival ALMA observations of the HCN and CO $J=1-0$ transitions, in addition to the radio continuum at 93 GHz, to assess the relationship between dense gas, star formation, and gas dynamics in ten, nearby (U)LIRGs and late-type galaxy centers. We frame our results in the context of turbulent and gravoturbulent models of star formation to assess if the HCN/CO ratio tracks the gravitationall… ▽ More We use archival ALMA observations of the HCN and CO $J=1-0$ transitions, in addition to the radio continuum at 93 GHz, to assess the relationship between dense gas, star formation, and gas dynamics in ten, nearby (U)LIRGs and late-type galaxy centers. We frame our results in the context of turbulent and gravoturbulent models of star formation to assess if the HCN/CO ratio tracks the gravitationally-bound, star-forming gas in molecular clouds ($f_\mathrm{grav}$) at sub-kpc scales in nearby galaxies. We confirm that the HCN/CO ratio is a tracer of gas above $n_\mathrm{SF}\approx10^{4.5}$ cm$^{-3}$, but the sub-kpc variations in HCN/CO do not universally track $f_\mathrm{grav}$. We find strong evidence for the use of varying star formation density threshold models, which are able to reproduce trends observed in $t_\mathrm{dep}$ and $ε_\mathrm{ff}$ that fixed threshold models cannot. Composite lognormal and powerlaw models outperform pure lognormal models in reproducing the observed trends, even when using a fixed powerlaw slope. The ability of the composite models to better reproduce star formation properties of the gas provides additional indirect evidence that the star formation efficiency per free-fall time is proportional to the fraction of gravitationally-bound gas. △ Less

Submitted 16 January, 2023; originally announced January 2023.

Comments: 23 pages, 11 figures, 1 appendix, accepted for publication in ApJ

PHANGS-JWST First Results: A Global and Moderately Resolved View of Mid-Infrared and CO Line Emission from Galaxies at the Start of the JWST Era

Authors: Adam K. Leroy, Alberto D. Bolatto, Karin Sandstrom, Erik Rosolowsky, Ashley. T. Barnes, F. Bigiel, Médéric Boquien, Jakob S. den Brok, Yixian Cao, Jérémy Chastenet, Mélanie Chevance, I-Da Chiang, Ryan Chown, Dario Colombo, Sara L. Ellison, Eric Emsellem, Kathryn Grasha, Jonathan D. Henshaw, Annie Hughes, Ralf S. Klessen, Eric W. Koch, Jaeyeon Kim, Kathryn Kreckel, J. M. Diederik Kruijssen, Kirsten L. Larson , et al. (19 additional authors not shown)

Abstract: We explore the relationship between mid-infrared (mid-IR) and CO rotational line emission from massive star-forming galaxies, which is one of the tightest scalings in the local universe. We assemble a large set of unresolved and moderately ($\sim 1$ kpc) spatially resolved measurements of CO (1-0) and CO (2-1) intensity, $I_{\rm CO}$, and mid-IR intensity, $I_{\rm MIR}$, at 8, 12, 22, and 24$μ$m.… ▽ More We explore the relationship between mid-infrared (mid-IR) and CO rotational line emission from massive star-forming galaxies, which is one of the tightest scalings in the local universe. We assemble a large set of unresolved and moderately ($\sim 1$ kpc) spatially resolved measurements of CO (1-0) and CO (2-1) intensity, $I_{\rm CO}$, and mid-IR intensity, $I_{\rm MIR}$, at 8, 12, 22, and 24$μ$m. The $I_{\rm CO}$ vs. $I_{\rm MIR}$ relationship is reasonably described by a power law with slopes $0.7{-}1.2$ and normalization $I_{\rm CO} \sim 1$ K km s$^{-1}$ at $I_{\rm MIR} \sim 1$ MJy sr$^{-1}$. Both the slopes and intercepts vary systematically with choice of line and band. The comparison between the relations measured for CO~(1-0) and CO (2-1) allow us to infer that $R_{21} \propto I_{\rm MIR}^{0.2}$, in good agreement with other work. The $8μ$m and $12μ$m bands, with strong PAH features, show steeper CO vs. mid-IR slopes than the $22μ$m and $24μ$m, consistent with PAH emission arising not just from CO-bright gas but also from atomic or CO-dark gas. The CO-to-mid-IR ratio correlates with global galaxy stellar mass ($M_\star$) and anti-correlates with SFR/$M_\star$. At $\sim 1$ kpc resolution, the first four PHANGS-JWST targets show CO to mid-IR relationships that are quantitatively similar to our larger literature sample, including showing the steep CO-to-mid-IR slopes for the JWST PAH-tracing bands, although we caution that these initial data have a small sample size and span a limited range of intensities. △ Less

Submitted 27 December, 2022; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 29 pages, 13 figures, key quantitative results in Table 3, Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

Challenges in Constraining Gravity with Cosmic Voids

Authors: Christopher Wilson, Rachel Bean

Abstract: We compare void size and clustering statistics for nDGP and $f(R)$ gravity models and GR using N-body simulations. We show how it is critical to consider the statistics derived from mock galaxy catalogs rather than the dark matter halos alone. Marked differences between the void size functions for GR and $f(R)$ models which present when voids are identified using dark matter halos are removed when… ▽ More We compare void size and clustering statistics for nDGP and $f(R)$ gravity models and GR using N-body simulations. We show how it is critical to consider the statistics derived from mock galaxy catalogs rather than the dark matter halos alone. Marked differences between the void size functions for GR and $f(R)$ models which present when voids are identified using dark matter halos are removed when voids are identified, more realistically, from mock galaxy tracers of the halos. The void radial velocities and velocity dispersions in the $f(R)$ and nDGP models are enhanced relative to GR in both halos and mock galaxy identified voids. Despite this, we find that the redshift space void quadrupole moments derived from the mock galaxy tracers are strikingly similar across the three gravity models. The Gaussian Streaming Model (GSM) is shown to accurately reconstruct $ξ_2$ in modified gravity models and we employ the GSM, using a functional derivative approach, to analyze the insensitivity of $ξ_2$ to the gravity model. Assuming linear theory, we show the void quadrupole to be an unbiased estimator of the redshift space growth rate parameter $β=f/b$ in the modified gravity theories. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: 20 pages, 13 figures

VERTICO III: The Kennicutt-Schmidt relation in Virgo cluster galaxies

Authors: M. J. Jiménez-Donaire, T. Brown, C. D. Wilson, I. D. Roberts, N. Zabel, S. L. Ellison, M. Thorp, V. Villanueva, R. Chown, D. Bisaria, A. D. Bolatto, A. Boselli, B. Catinella, A. Chung, L. Cortese, T. A. Davis, C. D. P. Lagos, B. Lee, L. C. Parker, K. Spekkens, A. R. H. Stevens, J. Sun

Abstract: In this VERTICO science paper we aim to study how the star formation process depends on galactic environment and gravitational interactions in the context of galaxy evolution. We explore the scaling relation between the star formation rate (SFR) surface density and the molecular gas surface density, also known as the Kennicutt-Schmidt (KS) relation, in a subsample of Virgo cluster spiral galaxies.… ▽ More In this VERTICO science paper we aim to study how the star formation process depends on galactic environment and gravitational interactions in the context of galaxy evolution. We explore the scaling relation between the star formation rate (SFR) surface density and the molecular gas surface density, also known as the Kennicutt-Schmidt (KS) relation, in a subsample of Virgo cluster spiral galaxies. We use new ACA and TP observations from the VERTICO-ALMA Large Program at 720pc resolution to resolve the molecular gas content, as traced by the 12CO(2-1) transition, across the disks of 37 spiral galaxies in the Virgo cluster. In combination with archival observations, we estimate the parameters of the KS relation for the entire ensemble of galaxies, and within individual galaxies. We find the KS slope for the entire population to be N=0.97+/-0.07, with a characteristic molecular gas depletion time of 1.86Gyr for our full sample, in agreement with previous work in isolated star-forming galaxies. In individual galaxies, we find KS slopes ranging between 0.69 and 1.40, and typical star formation efficiencies (SFE) that can vary from galaxy to galaxy by a factor of ~4. These galaxy-to-galaxy variations account for ~0.20dex in scatter in the ensemble KS relation, which is characterized by a 0.42dex scatter. We find that the HI-deficient galaxies in the Virgo cluster show a steeper resolved KS relation and lower molecular gas efficiencies than HI-normal cluster galaxies. While the molecular gas content in Virgo cluster galaxies appears to behave similarly to that in isolated galaxies, our VERTICO sample shows that cluster environments play a key role in regulating star formation. The environmental mechanisms affecting the HI galaxy content also have a direct impact in the SFE of molecular gas in cluster galaxies, leading to longer depletion times in HI-deficient members. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 671, A3 (2023)

Tracking ALMA System Temperature with Water Vapor Data at High Frequency

Authors: Hao He, William R. F. Dent, Christine Wilson

Abstract: The ALMA observatory is now putting more focus on high-frequency observations (frequencies from 275-950 GHz). However, high-frequency observations often suffer from rapid variations in atmospheric opacity that directly affect the system temperature $T_{sys}$. Current observations perform discrete atmospheric calibrations (Atm-cals) every few minutes, with typically 10-20 occurring per hour for hig… ▽ More The ALMA observatory is now putting more focus on high-frequency observations (frequencies from 275-950 GHz). However, high-frequency observations often suffer from rapid variations in atmospheric opacity that directly affect the system temperature $T_{sys}$. Current observations perform discrete atmospheric calibrations (Atm-cals) every few minutes, with typically 10-20 occurring per hour for high frequency observation and each taking 30-40 seconds. In order to obtain more accurate flux measurements and reduce the number of atmospheric calibrations (Atm-cals), a new method to monitor $T_{sys}$ continuously is proposed using existing data in the measurement set. In this work, we demonstrate the viability of using water vapor radiometer (WVR) data to track the $T_{sys}$ continuously. We find a tight linear correlation between $T_{sys}$ measured using the traditional method and $T_{sys}$ extrapolated based on WVR data with scatter of 0.5-3%. Although the exact form of the linear relation varies among different data sets and spectral windows, we can use a small number of discrete $T_{sys}$ measurements to fit the linear relation and use this heuristic relationship to derive $T_{sys}$ every 10 seconds. Furthermore, we successfully reproduce the observed correlation using atmospheric transmission at microwave (ATM) modeling and demonstrate the viability of a more general method to directly derive the $T_{sys}$ from the modeling. We apply the semi-continuous $T_{sys}$ from heuristic fitting on a few data sets from Band 7 to Band 10 and compare the flux measured using these methods. We find the discrete and continuous $T_{sys}$ methods give us consistent flux measurements with differences up to 5%. Furthermore, this method has significantly reduced the flux uncertainty due to $T_{sys}$ variability for one dataset, which has large precipitable water vapor (PWV) fluctuation, from 10% to 0.7%. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Comments: 24 pages, 18 figures, accepted to PASP

VERTICO IV: Environmental Effects on the Gas Distribution and Star Formation Efficiency of Virgo Cluster Spirals

Authors: Vicente Villanueva, Alberto D. Bolatto, Stuart Vogel, Tobias Brown, Christine D. Wilson, Nikki Zabel, Sara Ellison, Adam R. H. Stevens, Maria Jesus Jimenez Donaire, Kristine Spekkens, Mallory Thorp, Timothy A. Davis, Laura C. Parker, Ian D. Roberts, Dhruv Bisaria, Alessandro Boselli, Barbara Catinella, Aeree Chung, Luca Cortese, Bumhyun Lee, Adam Watts

Abstract: We measure the molecular-to-atomic gas ratio, $R_{\rm mol}$, and the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}$, in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We determine their scale-lengths for the molecular and stellar components and find a roughly 3:5 ratio between them compared to $\sim$1:1 in field galaxies, indicating th… ▽ More We measure the molecular-to-atomic gas ratio, $R_{\rm mol}$, and the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}$, in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We determine their scale-lengths for the molecular and stellar components and find a roughly 3:5 ratio between them compared to $\sim$1:1 in field galaxies, indicating that the CO emission is more centrally concentrated than the stars. We compute $R_{\rm mol}$ as a function of different physical quantities. While the spatially-resolved $R_{\rm mol}$ on average decreases with increasing radius, we find that the mean molecular-to-atomic gas ratio within the stellar effective radius $R_{\rm e}$, $R_{\rm mol}(r<R_{\rm e})$, shows a systematic increase with the level of H$_{\rm I}$, truncation and/or asymmetry (H$_{\rm I}$ perturbation). Analysis of the molecular- and the atomic-to-stellar mass ratios within $R_{\rm e}$, $R^{\rm mol}_{\star}(r<R_{\rm e})$ and $R^{\rm atom}_{\star}(r<R_{\rm e})$, shows that VERTICO galaxies have increasingly lower $R^{\rm atom}_{\star}(r<R_{\rm e})$ for larger levels of H$_{\rm I}$perturbation (compared to field galaxies matched in stellar mass), but no significant change in $R^{\rm mol}_{\star}(r<R_{\rm e})$. We also measure a clear systematic decrease of the SFE$_{\rm mol}$ within $R_{\rm e}$, SFE$_{\rm mol}(r<R_{\rm e})$, with increasingly perturbed H$_{\rm I}$. Therefore, compared to galaxies from the field, VERTICO galaxies are more compact in CO emission in relation to their stellar distribution, but increasingly perturbed atomic gas increases their $R_{\rm mol}$ and decreases the efficiency with which their molecular gas forms stars. (abridged) △ Less

Submitted 1 November, 2022; v1 submitted 11 October, 2022; originally announced October 2022.

Investigating 2M1155-7919B: a Nearby, Young, Low-Mass Star Actively Accreting from a Nearly Edge-on, Dusty Disk

Authors: D. Annie Dickson-Vandervelde, Joel H. Kastner, Jonathan Gagné, Adam C. Schneider, Jacqueline Faherty, Emily C. Wilson, Christophe Pinte, Francois Ménard

Abstract: We investigate the nature of an unusually faint member of the $ε$ Cha Association ($D\sim100$ pc, age $\sim5$ Myr), the nearest region of star formation of age $<$8 Myr. This object, 2MASS J11550336-7919147 (2M1155$-$79B), is a wide ($\sim$580 AU) separation, comoving companion to low-mass (M3) $ε$ Cha Association member 2MASS J11550485-7919108 (2M1155$-$79A). We present near-infrared spectra of b… ▽ More We investigate the nature of an unusually faint member of the $ε$ Cha Association ($D\sim100$ pc, age $\sim5$ Myr), the nearest region of star formation of age $<$8 Myr. This object, 2MASS J11550336-7919147 (2M1155$-$79B), is a wide ($\sim$580 AU) separation, comoving companion to low-mass (M3) $ε$ Cha Association member 2MASS J11550485-7919108 (2M1155$-$79A). We present near-infrared spectra of both components, along with analysis of photometry from Gaia EDR3, 2MASS, VHS, and WISE. The near-IR spectrum of 2M1155$-$79B displays strong He I 1.083 emission, a sign of active accretion and/or accretion-driven winds from a circumstellar disk. Analysis of WISE archival data reveals that the mid-infrared excess previously associated with 2M1155$-$79A instead originates from the disk surrounding 2M1155$-$79B. Based on these results, as well as radiative transfer modeling of its optical/IR spectral energy distribution, we conclude that 2M1155$-$79B is most likely a young, late-M, star that is partially obscured by, and actively accreting from, a nearly edge-on circumstellar disk. This would place 2M1155$-$79B among the rare group of nearby ($D\lesssim100$ pc), young (age $<$10 Myr) mid-M stars that are orbited by and accreting from highly inclined protoplanetary disks. Like these systems, the 2M1155$-$79B system is a particularly promising subject for studies of star and planet formation around low-mass stars. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 11 pages, 6 figures, accepted for publication in AJ

The cold gas and dust properties of red star-forming galaxies

Authors: Ryan Chown, Laura C. Parker, Christine D. Wilson, Toby Brown, Fraser A. Evans, Yang Gao, Ho Seong Hwang, Lihwai Lin, Amelie Saintonge, Mark Sargent, Matthew W. L. Smith, Ting Xiao

Abstract: We study the cold gas and dust properties for a sample of red star forming galaxies called "red misfits." We collect single-dish CO observations and HI observations from representative samples of low-redshift galaxies, as well as our own JCMT CO observations of red misfits. We also obtain SCUBA-2 850 um observations for a subset of these galaxies. With these data we compare the molecular gas, tota… ▽ More We study the cold gas and dust properties for a sample of red star forming galaxies called "red misfits." We collect single-dish CO observations and HI observations from representative samples of low-redshift galaxies, as well as our own JCMT CO observations of red misfits. We also obtain SCUBA-2 850 um observations for a subset of these galaxies. With these data we compare the molecular gas, total cold gas, and dust properties of red misfits against those of their blue counterparts ("blue actives") taking non-detections into account using a survival analysis technique. We compare these properties at fixed position in the log SFR-log M* plane, as well as versus offset from the star-forming main sequence. Compared to blue actives, red misfits have slightly longer molecular gas depletion times, similar total gas depletion times, significantly lower molecular- and total-gas mass fractions, lower dust-to-stellar mass ratios, similar dust-to-gas ratios, and a significantly flatter slope in the $\log M_\mathrm{mol}$-$\log M_\star$ plane. Our results suggest that red misfits as a population are likely quenching due to a shortage in gas supply. △ Less

Submitted 24 August, 2022; v1 submitted 7 August, 2022; originally announced August 2022.

Comments: 16 pages, 7 Figures, accepted to MNRAS

Extreme giant molecular clouds in the luminous infrared galaxy NGC 3256

Authors: Nathan Brunetti, Christine D. Wilson

Abstract: (Abridged) We present a cloud decomposition of $^{12}$CO (2--1) observations of the merger and nearest luminous infrared galaxy, NGC 3256. 185 spatially and spectrally resolved clouds are identified across the central $\approx$ 130 kpc$^{2}$ at 90 pc resolution and completeness is estimated. We compare our cloud catalogue from NGC 3256 to ten galaxies observed in the PHANGS-ALMA survey. Distributi… ▽ More (Abridged) We present a cloud decomposition of $^{12}$CO (2--1) observations of the merger and nearest luminous infrared galaxy, NGC 3256. 185 spatially and spectrally resolved clouds are identified across the central $\approx$ 130 kpc$^{2}$ at 90 pc resolution and completeness is estimated. We compare our cloud catalogue from NGC 3256 to ten galaxies observed in the PHANGS-ALMA survey. Distributions in NGC 3256 of cloud velocity dispersions, luminosities, CO-estimated masses, mass surface densities, virial masses, virial parameters, size-linewidth coefficients, and internal turbulent pressures are significantly higher than in the PHANGS-ALMA galaxies. Cloud radii are slightly larger in NGC 3256 and free-fall times are shorter. The distribution of cloud eccentricities in NGC 3256 is indistinguishable from many PHANGS-ALMA galaxies, possibly because the dynamical state of clouds in NGC 3256 is similar to that of nearby spiral galaxies. However, the narrower distribution of virial parameters in NGC 3256 may reflect a narrower range of dynamical states than in PHANGS-ALMA galaxies. No clear picture of cloud alignment is detected, despite the large eccentricities. Correlations between cloud properties point to high external pressures in NGC 3256 keeping clouds bound and collapsing given such high velocity dispersions and star-formation rates. A fit to the cloud mass function gives a high-mass power-law slope of $-2.75^{+0.07}_{-0.01}$, near the average from PHANGS-ALMA galaxies. We also compare our results to a pixel-based analysis of these observations and find molecular-gas properties agree qualitatively, though peak brightness temperatures are somewhat higher and virial parameters and free-fall times are somewhat lower in this cloud-based analysis. △ Less

Submitted 11 July, 2022; originally announced July 2022.

Comments: accepted to MNRAS; 25 pages, 17 figures

On the scale-height of the molecular gas disc in Milky Way-like galaxies

Authors: Sarah M. R. Jeffreson, Jiayi Sun, Christine D. Wilson

Abstract: We study the relationship between the scale-height of the molecular gas disc and the turbulent velocity dispersion of the molecular interstellar medium within a simulation of a Milky Way-like galaxy in the moving-mesh code Arepo. We find that the vertical distribution of molecular gas can be described by a Gaussian function with a uniform scale-height of ~50 pc. We investigate whether this scale-h… ▽ More We study the relationship between the scale-height of the molecular gas disc and the turbulent velocity dispersion of the molecular interstellar medium within a simulation of a Milky Way-like galaxy in the moving-mesh code Arepo. We find that the vertical distribution of molecular gas can be described by a Gaussian function with a uniform scale-height of ~50 pc. We investigate whether this scale-height is consistent with a state of hydrostatic balance between gravity and turbulent pressure. We find that the hydrostatic prediction using the total turbulent velocity dispersion (as one would measure from kpc-scale observations) gives an over-estimate of the true molecular disc scale-height. The hydrostatic prediction using the velocity dispersion between the centroids of discrete giant molecular clouds (cloud-cloud velocity dispersion) leads to more-accurate estimates. The velocity dispersion internal to molecular clouds is elevated by the locally-enhanced gravitational field. Our results suggest that observations of molecular gas need to reach the scale of individual molecular clouds in order to accurately determine the molecular disc scale-height. △ Less

Submitted 29 June, 2022; originally announced June 2022.

Comments: MNRAS accepted, comments welcome. 14 pages, 10 figures

VERTICO II: effects of HI-identified environmental mechanisms on molecular gas

Authors: Nikki Zabel, Toby Brown, Christine D. Wilson, Timothy A. Davis, Luca Cortese, Laura C. Parker, Alessandro Boselli, Barbara Catinella, Ryan Chown, Aeree Chung, Tirna Deb, Sara L. Ellison, María J. Jiménez-Donaire, Bumhyun Lee, Ian D. Roberts, Kristine Spekkens, Adam R. H. Stevens, Mallory Thorp, Stephanie Tonnesen, Vicente Villanueva

Abstract: In this VERTICO early science paper we explore in detail how environmental mechanisms, identified in HI, affect the resolved properties of molecular gas reservoirs in cluster galaxies. The molecular gas is probed using ALMA ACA (+TP) observations of 12CO(2-1) in 51 spiral galaxies in the Virgo cluster (of which 49 are detected), all of which are included in the VIVA HI survey. The sample spans a s… ▽ More In this VERTICO early science paper we explore in detail how environmental mechanisms, identified in HI, affect the resolved properties of molecular gas reservoirs in cluster galaxies. The molecular gas is probed using ALMA ACA (+TP) observations of 12CO(2-1) in 51 spiral galaxies in the Virgo cluster (of which 49 are detected), all of which are included in the VIVA HI survey. The sample spans a stellar mass range of 9 < log M*/Msol < 11. We study molecular gas radial profiles, isodensity radii, and surface densities as a function of galaxy HI deficiency and morphology. There is a weak correlation between global HI and H2 deficiencies, and resolved properties of molecular gas correlate with HI deficiency: galaxies that have large HI deficiencies have relatively steep and truncated molecular gas radial profiles, which is due to the removal of low-surface density molecular gas on the outskirts. Therefore, while the environmental mechanisms observed in HI also affect molecular gas reservoirs, there is only a moderate reduction of the total amount of molecular gas. △ Less

Submitted 12 January, 2023; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: Published in ApJ. 22 pages, 6 figures, 1 table, 1 appendix. Erratum accepted for publication in ApJ

A direct measurement of the distance to the Galactic center using the kinematics of bar stars

Authors: Henry W. Leung, Jo Bovy, J. Ted Mackereth, Jason A. S. Hunt, Richard R. Lane, John C. Wilson

Abstract: The distance to the Galactic center $R_0$ is a fundamental parameter for understanding the Milky Way, because all observations of our Galaxy are made from our heliocentric reference point. The uncertainty in $R_0$ limits our knowledge of many aspects of the Milky Way, including its total mass and the relative mass of its major components, and any orbital parameters of stars employed in chemo-dynam… ▽ More The distance to the Galactic center $R_0$ is a fundamental parameter for understanding the Milky Way, because all observations of our Galaxy are made from our heliocentric reference point. The uncertainty in $R_0$ limits our knowledge of many aspects of the Milky Way, including its total mass and the relative mass of its major components, and any orbital parameters of stars employed in chemo-dynamical analyses. While measurements of $R_0$ have been improving over a century, measurements in the past few years from a variety of methods still find a wide range of $R_0$ being somewhere within $8.0$ to $8.5\,\mathrm{kpc}$. The most precise measurements to date have to assume that Sgr A$^*$ is at rest at the Galactic center, which may not be the case. In this paper, we use maps of the kinematics of stars in the Galactic bar derived from APOGEE DR17 and Gaia EDR3 data augmented with spectro-photometric distances from the \texttt{astroNN} neural-network method. These maps clearly display the minimum in the rotational velocity $v_T$ and the quadrupolar signature in radial velocity $v_R$ expected for stars orbiting in a bar. From the minimum in $v_T$, we measure $R_0 = 8.23 \pm 0.12\,\mathrm{kpc}$. We validate our measurement using realistic $N$-body simulations of the Milky Way. We further measure the pattern speed of the bar to be $Ω_\mathrm{bar} = 40.08\pm1.78\,\mathrm{km\,s}^{-1}\mathrm{kpc}^{-1}$. Because the bar forms out of the disk, its center is manifestly the barycenter of the bar+disc system and our measurement is therefore the most robust and accurate measurement of $R_0$ to date. △ Less

Submitted 26 April, 2022; originally announced April 2022.

Convection Reconciles the Difference in Efficiencies Between Low-Mass and High-Mass Common Envelopes

Authors: E. C. Wilson, J. Nordhaus

Abstract: The formation pathways for gravitational-wave merger sources are predicted to include common envelope (CE) evolution. Observations of high-mass post-common envelope binaries suggest that energy transfer to the envelope during the CE phase must be highly efficient. In contrast, observations of low-mass post-CE binaries indicate energy transfer during the CE phase must be highly inefficient. Convect… ▽ More The formation pathways for gravitational-wave merger sources are predicted to include common envelope (CE) evolution. Observations of high-mass post-common envelope binaries suggest that energy transfer to the envelope during the CE phase must be highly efficient. In contrast, observations of low-mass post-CE binaries indicate energy transfer during the CE phase must be highly inefficient. Convection, a process present in low-mass and high-mass stars, naturally explains this dichotomy. Using observations of Wolf-Rayet binaries, we study the effects of convection and radiative losses on the predicted final separations of high-mass common envelopes. Despite robust convection in massive stars, the effect is minimal as the orbit decays well before convection can transport the liberated orbital energy to the surface. In low-mass systems, convective transport occurs faster then the orbit decays, allowing the system to radiatively cool thereby lowering the efficiency. The inclusion of convection reproduces observations of low-mass and high-mass binaries and remains a necessary ingredient for determining outcomes of common envelopes. △ Less

Submitted 30 August, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: Accepted to MNRAS

Embedded Young Massive Star Clusters in the Antennae Merger

Authors: Hao He, Christine Wilson, Nathan Brunetti, Molly Finn, Ashley Bemis, Kelsey Johnson

Abstract: The properties of young massive clusters (YMCs) are key to understanding the star formation mechanism in starburst systems, especially mergers. We present ALMA high-resolution ($\sim$10 pc) continuum (100 and 345 GHz) data of YMCs in the overlap region of the Antennae galaxy. We identify 6 sources in the overlap region, including two sources that lie in the same giant molecular cloud (GMC). These… ▽ More The properties of young massive clusters (YMCs) are key to understanding the star formation mechanism in starburst systems, especially mergers. We present ALMA high-resolution ($\sim$10 pc) continuum (100 and 345 GHz) data of YMCs in the overlap region of the Antennae galaxy. We identify 6 sources in the overlap region, including two sources that lie in the same giant molecular cloud (GMC). These YMCs correspond well with radio sources in lower resolution continuum (100 and 220 GHz) images at GMC scales ($\sim$60 pc). We find most of these YMCs are bound clusters through virial analysis. We estimate their ages to be $\sim$1 Myr and to be either embedded or just beginning to emerge from their parent cloud. We also compare each radio source with Pa$β$ source and find they have consistent total ionizing photon numbers, which indicates they are tracing the same physical source. By comparing the free-free emission at $\sim$10 pc scale and $\sim$60 pc scale, we find that $\sim$50% of the free-free emission in GMCs actually comes from these YMCs. This indicates that roughly half of the stars in massive GMCs are formed in bound clusters. We further explore the mass correlation between YMCs and GMCs in the Antennae and find it generally agrees with the predictions of the star cluster simulations. The most massive YMC has a stellar mass that is 1% - 5% of its host GMC mass. △ Less

Submitted 23 February, 2022; v1 submitted 16 February, 2022; originally announced February 2022.

Comments: 20 pages, 11 figures; accepted for publication in ApJ; added adjunct affliation

The Seventeenth Data Release of the Sloan Digital Sky Surveys: Complete Release of MaNGA, MaStar and APOGEE-2 Data

Authors: Abdurro'uf, Katherine Accetta, Conny Aerts, Victor Silva Aguirre, Romina Ahumada, Nikhil Ajgaonkar, N. Filiz Ak, Shadab Alam, Carlos Allende Prieto, Andres Almeida, Friedrich Anders, Scott F. Anderson, Brett H. Andrews, Borja Anguiano, Erik Aquino-Ortiz, Alfonso Aragon-Salamanca, Maria Argudo-Fernandez, Metin Ata, Marie Aubert, Vladimir Avila-Reese, Carles Badenes, Rodolfo H. Barba, Kat Barger, Jorge K. Barrera-Ballesteros, Rachael L. Beaton , et al. (316 additional authors not shown)

Abstract: This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies… ▽ More This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies this data, providing observations of almost 30,000 stars through the MaNGA instrument during bright time. DR17 also contains the complete release of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) survey which publicly releases infra-red spectra of over 650,000 stars. The main sample from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), as well as the sub-survey Time Domain Spectroscopic Survey (TDSS) data were fully released in DR16. New single-fiber optical spectroscopy released in DR17 is from the SPectroscipic IDentification of ERosita Survey (SPIDERS) sub-survey and the eBOSS-RM program. Along with the primary data sets, DR17 includes 25 new or updated Value Added Catalogs (VACs). This paper concludes the release of SDSS-IV survey data. SDSS continues into its fifth phase with observations already underway for the Milky Way Mapper (MWM), Local Volume Mapper (LVM) and Black Hole Mapper (BHM) surveys. △ Less

Submitted 13 January, 2022; v1 submitted 3 December, 2021; originally announced December 2021.

Comments: 40 pages, 8 figures, 6 tables. In press at ApJSS (arxiv v2 corrects some minor typos and updates references)

No detection of SO2, H2S, or OCS in the atmosphere of Mars from the first two Martian years of observations from TGO/ACS

Authors: Ashwin S. Braude, F. Montmessin, K. S. Olsen, A. Trokhimovskiy, O. I. Korablev, F. Lefèvre, A. A. Fedorova, J. Alday, L. Baggio, A. Irbah, G. Lacombe, F. Forget, E. Millour, C. F. Wilson, A. Patrakeev, A. Shakun

Abstract: The detection of sulphur species in the Martian atmosphere would be a strong indicator of volcanic outgassing from the surface of Mars. We wish to establish the presence of SO2, H2S, or OCS in the Martian atmosphere or determine upper limits on their concentration in the absence of a detection. We perform a comprehensive analysis of solar occultation data from the mid-infrared channel of the Atmos… ▽ More The detection of sulphur species in the Martian atmosphere would be a strong indicator of volcanic outgassing from the surface of Mars. We wish to establish the presence of SO2, H2S, or OCS in the Martian atmosphere or determine upper limits on their concentration in the absence of a detection. We perform a comprehensive analysis of solar occultation data from the mid-infrared channel of the Atmospheric Chemistry Suite instrument, on board the ExoMars Trace Gas Orbiter, obtained during Martian years 34 and 35. For the most optimal sensitivity conditions, we determine 1-sigma upper limits of SO2 at 20 ppbv, H2S at 15 ppbv, and OCS at 0.4 ppbv; the last value is lower than any previous upper limits imposed on OCS in the literature. We find no evidence of any of these species above a 3-sigma confidence threshold. We therefore infer that passive volcanic outgassing of SO2 must be below 2 ktons/day. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: 13 pages, 12 figures, received for production in A&A on the 19th of November 2021

Journal ref: A&A 658, A86 (2022)

VERTICO: The Virgo Environment Traced In CO Survey

Authors: Toby Brown, Christine D. Wilson, Nikki Zabel, Timothy A. Davis, Alessandro Boselli, Aeree Chung, Sara L. Ellison, Claudia D. P. Lagos, Adam R. H. Stevens, Luca Cortese, Yannick M. Bahé, Dhruv Bisaria, Alberto D. Bolatto, Claire R. Cashmore, Barbara Catinella, Ryan Chown, Benedikt Diemer, Pascal J. Elahi, Maan H. Hani, María J. Jiménez-Donaire, Bumhyun Lee, Katya Leidig, Angus Mok, Karen Pardos Olsen, Laura C. Parker , et al. (11 additional authors not shown)

Abstract: We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map $^{12}$CO($2-1$), $^{13}$CO($2-1$), and C$^{18}$O($2-1$) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array (ALMA). The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star forma… ▽ More We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map $^{12}$CO($2-1$), $^{13}$CO($2-1$), and C$^{18}$O($2-1$) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array (ALMA). The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star formation and galaxy evolution, in dense environments. This first paper contains an overview of VERTICO's design and sample selection, $^{12}$CO($2-1$) observations, and data reduction procedures. We characterize global $^{12}$CO($2-1$) fluxes and molecular gas masses for the 49 detected VERTICO galaxies, provide upper limits for the two non-detections, and produce resolved $^{12}$CO($2-1$) data products (median resolution $= 8^{\prime\prime} \approx 640~{\rm pc}$). Azimuthally averaged $^{12}$CO($2-1$) radial intensity profiles are presented along with derived molecular gas radii. We demonstrate the scientific power of VERTICO by comparing the molecular gas size--mass scaling relation for our galaxies with a control sample of field galaxies, highlighting the strong effect that radius definition has on this correlation. We discuss the drivers of the form and scatter in the size--mass relation and highlight areas for future work. VERTICO is an ideal resource for studying the fate of molecular gas in cluster galaxies and the physics of environment-driven processes that perturb the star formation cycle. Upon public release, the survey will provide a homogeneous legacy dataset for studying galaxy evolution in our closest cluster. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 68 pages, 13 Figures, 2 Figure Sets, Accepted for publication in ApJS, Online FITS versions of Tables 1, 2, and 3 are available with the journal publication

The HASHTAG project: The First Submillimeter Images of the Andromeda Galaxy from the Ground

Authors: Matthew W. L. Smith, Stephen A. Eales, Thomas G. Williams, Bumhyun Lee, Zongnan Li, Pauline Barmby, Martin Bureau, Scott Chapman, Brian S. Cho, Aeree Chung, Eun Jung Chung, Hui-Hsuan Chung, Christopher J. R. Clark, David L. Clements, Timothy A. Davis, Ilse De Looze, David J. Eden, Gayathri Athikkat-Eknath, George P. Ford, Yu Gao, Walter Gear, Haley L. Gomez, Richard de Grijs, Jinhua He, Luis C. Ho , et al. (24 additional authors not shown)

Abstract: Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing… ▽ More Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing the galaxy's large-scale structure. We have developed a technique for producing high-resolution submillimeter images of galaxies of large angular size by using the telescope on the ground to determine the small-scale structure (the large Fourier components) and a space telescope (Herschel or Planck) to determine the large-scale structure (the small Fourier components). Using this technique, we are carrying out the HARP and SCUBA-2 High Resolution Terahertz Andromeda Galaxy Survey (HASHTAG), an international Large Program on the James Clerk Maxwell Telescope, with one aim being to produce the first high-fidelity high-resolution submillimeter images of Andromeda. In this paper, we describe the survey, the method we have developed for combining the space-based and ground-based data, and present the first HASHTAG images of Andromeda at 450 and 850um. We also have created a method to predict the CO(J=3-2) line flux across M31, which contaminates the 850um band. We find that while normally the contamination is below our sensitivity limit, the contamination can be significant (up to 28%) in a few of the brightest regions of the 10 kpc ring. We therefore also provide images with the predicted line emission removed. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Comments: 26 pages, 19 figures. Submitted to ApJS June 2021, Accepted September 2021

Chemical Cartography with APOGEE: Mapping Disk Populations with a Two-Process Model and Residual Abundances

Authors: David H. Weinberg, Jon A. Holtzman, Jennifer A. Johnson, Christian Hayes, Sten Hasselquist, Matthew Shetrone, Yuan-Sen Ting, Rachael L. Beaton, Timothy C. Beers, Jonathan C. Bird, Dmitry Bizyaev, Michael R. Blanton, Katia Cunha, Jose G. Fernandez-Trincado, Peter M. Frinchaboy, D. A. Garcia-Hernandez, Emily Griffith, James W. Johnson, Henrik Jonsson, Richard R. Lane, Henry W. Leung, J. Ted Mackereth, Steven R. Majewski, Szabolcz Meszaros, Christian Nitschelm , et al. (11 additional authors not shown)

Abstract: We apply a novel statistical analysis to measurements of 16 elemental abundances in 34,410 Milky Way disk stars from the final data release (DR17) of APOGEE-2. Building on recent work, we fit median abundance ratio trends [X/Mg] vs. [Mg/H] with a 2-process model, which decomposes abundance patterns into a "prompt" component tracing core collapse supernovae and a "delayed" component tracing Type Ia… ▽ More We apply a novel statistical analysis to measurements of 16 elemental abundances in 34,410 Milky Way disk stars from the final data release (DR17) of APOGEE-2. Building on recent work, we fit median abundance ratio trends [X/Mg] vs. [Mg/H] with a 2-process model, which decomposes abundance patterns into a "prompt" component tracing core collapse supernovae and a "delayed" component tracing Type Ia supernovae. For each sample star, we fit the amplitudes of these two components, then compute the residuals Δ[X/H] from this two-parameter fit. The rms residuals range from ~0.01-0.03 dex for the most precisely measured APOGEE abundances to ~0.1 dex for Na, V, and Ce. The correlations of residuals reveal a complex underlying structure, including a correlated element group comprised of Ca, Na, Al, K, Cr, and Ce and a separate group comprised of Ni, V, Mn, and Co. Selecting stars poorly fit by the 2-process model reveals a rich variety of physical outliers and sometimes subtle measurement errors. Residual abundances allow comparison of populations controlled for differences in metallicity and [α/Fe]. Relative to the main disk (R=3-13 kpc, |Z|<2 kpc), we find nearly identical abundance patterns in the outer disk (R=15-17 kpc), 0.05-0.2 dex depressions of multiple elements in LMC and Gaia Sausage/Enceladus stars, and wild deviations (0.4-1 dex) of multiple elements in ωCen. Residual abundance analysis opens new opportunities for discovering chemically distinctive stars and stellar populations, for empirically constraining nucleosynthetic yields, and for testing chemical evolution models that include stochasticity in the production and redistribution of elements. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Comments: Submitted to AAS journals. Long paper, many figures; see end of Section 1 for a reader's guide

Observed CN and HCN intensity ratios exhibit subtle variations in extreme galaxy environments

Authors: B. Ledger, C. D. Wilson, T. Michiyama, D. Iono, S. Aalto, T. Saito, A. Bemis, R. Aladro

Abstract: We use both new and archival ALMA data of three energy lines each of CN and HCN to explore intensity ratios in dense gas in NGC 3256, NGC 7469, and IRAS 13120-5453. The HCN (3-2)/HCN (1-0) intensity ratio varies in NGC 3256 and NGC 7469, with superlinear trends of 1.53$\pm$0.07 and 1.55$\pm$0.05, respectively. We find an offset to higher HCN (3-2)/HCN (1-0) intensity ratios (~0.8) in IRAS 13120-54… ▽ More We use both new and archival ALMA data of three energy lines each of CN and HCN to explore intensity ratios in dense gas in NGC 3256, NGC 7469, and IRAS 13120-5453. The HCN (3-2)/HCN (1-0) intensity ratio varies in NGC 3256 and NGC 7469, with superlinear trends of 1.53$\pm$0.07 and 1.55$\pm$0.05, respectively. We find an offset to higher HCN (3-2)/HCN (1-0) intensity ratios (~0.8) in IRAS 13120-5453 compared to NGC 3256 (~0.3-0.4) and NGC 7469 (~0.3-0.5). The HCN (4-3)/HCN (3-2) intensity ratio in NGC 7469 has a slope of 1.34$\pm$0.05. We attribute the variation within NGC 3256 to excitation associated with the northern and southern nuclei. In NGC 7469, the variations are localized to the region surrounding the active galactic nucleus. At our resolution (~700 pc), IRAS 13120-5453 shows little variation in the HCN intensity ratios. Individual galaxies show nearly constant CN (2-1)/CN (1-0) intensity ratios. We find an offset to lower CN (2-1)/CN (1-0) intensity ratios (~0.5) in NGC 3256 compared to the other two galaxies (~0.8). For the CN (3-2)/CN (2-1) intensity ratio, NGC 7469 has a superlinear trend of 1.55$\pm$0.04, with the peak localized toward the active galactic nucleus. We find high (~1.7) CN (1-0)/HCN (1-0) intensity ratios in IRAS 13120-5453 and in the northern nucleus of NGC 3256, compared to a more constant ratio (~1.1) in NGC 7469 and non-starbursting regions of NGC 3256. △ Less

Submitted 16 April, 2021; originally announced April 2021.

Comments: 12 pages, 4 figures, 5 figures in Appendix, resubmitted to MNRAS after small revisions to address the first referee report

PHANGS-ALMA Data Processing and Pipeline

Authors: Adam K. Leroy, Annie Hughes, Daizhong Liu, Jerome Pety, Erik Rosolowsky, Toshiki Saito, Eva Schinnerer, Andreas Schruba, Antonio Usero, Christopher M. Faesi, Cinthya N. Herrera, Melanie Chevance, Alexander P. S. Hygate, Amanda A. Kepley, Eric W. Koch, Miguel Querejeta, Kazimierz Sliwa, David Will, Christine D. Wilson, Gagandeep S. Anand, Ashley Barnes, Francesco Belfiore, Ivana Beslic, Frank Bigiel, Guillermo A. Blanc , et al. (43 additional authors not shown)

Abstract: We describe the processing of the PHANGS-ALMA survey and present the PHANGS-ALMA pipeline, a public software package that processes calibrated interferometric and total power data into science-ready data products. PHANGS-ALMA is a large, high-resolution survey of CO J=2-1 emission from nearby galaxies. The observations combine ALMA's main 12-m array, the 7-m array, and total power observations and… ▽ More We describe the processing of the PHANGS-ALMA survey and present the PHANGS-ALMA pipeline, a public software package that processes calibrated interferometric and total power data into science-ready data products. PHANGS-ALMA is a large, high-resolution survey of CO J=2-1 emission from nearby galaxies. The observations combine ALMA's main 12-m array, the 7-m array, and total power observations and use mosaics of dozens to hundreds of individual pointings. We describe the processing of the u-v data, imaging and deconvolution, linear mosaicking, combining interferometer and total power data, noise estimation, masking, data product creation, and quality assurance. Our pipeline has a general design and can also be applied to VLA and ALMA observations of other spectral lines and continuum emission. We highlight our recipe for deconvolution of complex spectral line observations, which combines multiscale clean, single scale clean, and automatic mask generation in a way that appears robust and effective. We also emphasize our two-track approach to masking and data product creation. We construct one set of "broadly masked" data products, which have high completeness but significant contamination by noise, and another set of "strictly masked" data products, which have high confidence but exclude faint, low signal-to-noise emission. Our quality assurance tests, supported by simulations, demonstrate that 12-m+7-m deconvolved data recover a total flux that is significantly closer to the total power flux than the 7-m deconvolved data alone. In the appendices, we measure the stability of the ALMA total power calibration in PHANGS--ALMA and test the performance of popular short-spacing correction algorithms. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: Accepted for publication in the Astrophysical Journal Supplement series. 65 pages, 33 figures. Software available at https://github.com/akleroy/phangs_imaging_scripts . For a full resolution version see https://sites.google.com/view/phangs/publications

The Case for Thermalization as a Contributor to the [CII] Deficit

Authors: Jessica Sutter, Daniel A. Dale, Karin Sandstrom, J. D. T. Smith, Alberto Bolatto, Mederic Boquien, Daniela Calzetti, Kevin V. Croxall, Ilse De Looze, Maude Galametz, Brent A. Groves, George Helou, Rodrigo Herrera-Camus, Leslie K. Hunt, Robert C. Kennicutt, Eric W. Pelligrini, Christine Wilson, Mark G. Wolfire

Abstract: The [CII] deficit, which describes the observed decrease in the ratio of [CII] 158 micron emission to continuum infrared emission in galaxies with high star formation surface densities, places a significant challenge to the interpretation of [CII] detections from across the observable universe. In an attempt to further decode the cause of the [CII] deficit, the [CII] and dust continuum emission fr… ▽ More The [CII] deficit, which describes the observed decrease in the ratio of [CII] 158 micron emission to continuum infrared emission in galaxies with high star formation surface densities, places a significant challenge to the interpretation of [CII] detections from across the observable universe. In an attempt to further decode the cause of the [CII] deficit, the [CII] and dust continuum emission from 18 Local Volume galaxies has been split based on conditions within the interstellar medium where it originated. This is completed using the Key Insights in Nearby Galaxies: a Far-Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) surveys and the wide-range of wavelength information, from UV to far-infrared emission lines, available for a selection of star-forming regions within these samples. By comparing these subdivided [CII] emissions to isolated infrared emission and other properties, we find that the thermalization (collisional de-excitation) of the [CII] line in HII regions plays a significant role in the deficit observed in our sample. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: To be published in MNRAS, 8 pages, 6 figures

Testing $f(R)$ Gravity With Scale Dependent Cosmic Void Velocity Profiles

Authors: Christopher Wilson, Rachel Bean

Abstract: We study the impact of cosmological scale modifications to General Relativity on the dynamics of halos within voids by comparing N-body simulations incorporating Hu-Sawicki $f(R)$ gravity, with $|f_{R0}|=10^{-6}$ and $10^{-5}$, to those of $Λ$CDM. By examining the radial velocity statistics within voids classified based on their size and density-profile, as "rising" ($R$-type) or "shell" ($S$-type… ▽ More We study the impact of cosmological scale modifications to General Relativity on the dynamics of halos within voids by comparing N-body simulations incorporating Hu-Sawicki $f(R)$ gravity, with $|f_{R0}|=10^{-6}$ and $10^{-5}$, to those of $Λ$CDM. By examining the radial velocity statistics within voids classified based on their size and density-profile, as "rising" ($R$-type) or "shell" ($S$-type), we find that halo motions in small $R$-type voids, with effective radius $<15 Mpc/h$, reveal distinctive differences between $f(R)$ and $Λ$CDM cosmologies. To understand this observed effect, we study the linear and nonlinear fifth forces, and develop an iterative algorithm to accurately solve the non-linear fifth force equation. We use this to characterize the Chameleon screening mechanism in voids and contrast the behavior with that observed in gravitationally collapsed objects. The force analysis underscores how smaller $R$-type voids exhibit the highest ratios of fifth force to Newtonian force, which source distinguishable differences in the velocity profiles and thereby provide rich environments in which to constrain gravity. △ Less

Submitted 9 June, 2021; v1 submitted 10 December, 2020; originally announced December 2020.

Comments: 16 pages, 9 figures, void catalogs including density and radial velocity profiles and analysis codes including method to calculate the 5th force within voids in f(R) gravity can be downloaded from https://github.com/CornellCosmology/Void_fR_velocity_profiles. Updated to match published version in PRD

Journal ref: Phys. Rev. D 104, 023512 (2021)

Gaia-based Isochronal, Kinematic, and Spatial Analysis of the $ε$ Cha Association

Authors: D. Annie Dickson-Vandervelde, Emily C. Wilson, Joel H. Kastner

Abstract: The precise parallax, proper motion, and photometric measurements contained in Gaia Data Release 2 (DR2) offer the opportunity to re-examine the membership and ages of Nearby Young Moving Groups (NYMGs), i.e., loose groups of stars of age $\lesssim100$ Myr in the solar vicinity. Here, we analyze the available DR2 data for members and candidate members of the $ε$ Cha Association ($ε$CA) which, at a… ▽ More The precise parallax, proper motion, and photometric measurements contained in Gaia Data Release 2 (DR2) offer the opportunity to re-examine the membership and ages of Nearby Young Moving Groups (NYMGs), i.e., loose groups of stars of age $\lesssim100$ Myr in the solar vicinity. Here, we analyze the available DR2 data for members and candidate members of the $ε$ Cha Association ($ε$CA) which, at an estimated age of $\sim$3--5 Myr, has previously been identified as among the youngest NYMGs. The several dozen confirmed members of $ε$CA include MP Mus and T Cha, two of the nearest stars of roughly solar mass that are known to host primordial protoplanetary disks, and the Herbig Ae/Be star HD 104237A. We have used Gaia DR2 data to ascertain the Galactic positions and kinematics and color-magnitude diagram positions of $ε$CA members and candidates, so as to reassess their membership status and thereby refine estimates of the distance, age, multiplicity, and disk fraction of the group. Our analysis yields 36 \textit{bona fide} $ε$CA members, as well as 20 provisional members, including 3 new members identified here as comoving companions to previously known $ε$CA stars. We determine a mean distance to $ε$CA of 101.0$\pm$4.6 pc and confirm that, at an age of $5^{+3}_{-2}$ Myr, $ε$CA represents the youngest stellar group within $\sim$100 pc of Earth. We identify several new photometric binary candidates, bringing the overall multiplicity fraction (MF) of $ε$CA to 40\%, intermediate between the MFs of young T associations and the field. △ Less

Submitted 3 December, 2020; v1 submitted 12 November, 2020; originally announced November 2020.

Comments: 26 pages, 9 figures, to be published in AJ, replacement with updated minor revisions

Highly turbulent gas on GMC-scales in NGC 3256, the nearest luminous infrared galaxy

Authors: Nathan Brunetti, Christine D. Wilson, Kazimierz Sliwa, Eva Schinnerer, Susanne Aalto, Alison B. Peck

Abstract: We present the highest resolution CO (2-1) observations obtained to date (0.25") of NGC 3256 and use them to determine the detailed properties of the molecular interstellar medium in the central 6 kpc of this merger. Distributions of physical quantities are reported from pixel-by-pixel measurements at 55 and 120 pc scales and compared to disc galaxies observed by PHANGS-ALMA. Mass surface densitie… ▽ More We present the highest resolution CO (2-1) observations obtained to date (0.25") of NGC 3256 and use them to determine the detailed properties of the molecular interstellar medium in the central 6 kpc of this merger. Distributions of physical quantities are reported from pixel-by-pixel measurements at 55 and 120 pc scales and compared to disc galaxies observed by PHANGS-ALMA. Mass surface densities range from 8 to 5500 M$_{\odot}$ pc$^{-2}$ and velocity dispersions from 10 to 200 km s$^{-1}$. Peak brightness temperatures as large as 37 K are measured, indicating the gas in NGC 3256 may be hotter than all regions in nearby disc galaxies measured by PHANGS-ALMA. Brightness temperatures even surpass those in the overlap region of NGC 4038/9 at the same scales. The majority of the gas appears unbound with median virial parameters of 7 to 19, although external pressure may bind some of the gas. High internal turbulent pressures of 10$^{5}$ to 10$^{10}$ K cm$^{-3}$ are found. Given the lack of significant trends in surface density, brightness temperature, and velocity dispersion with physical scale we argue the molecular gas is made up of a smooth medium down to 55 pc scales, unlike the more structured medium found in the PHANGS-ALMA disc galaxies. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: accepted to MNRAS, 21 pages, 8 figures

No evidence of phosphine in the atmosphere of Venus by independent analyses

Authors: Geronimo Villanueva, Martin Cordiner, Patrick Irwin, Imke de Pater, Bryan Butler, Mark Gurwell, Stefanie Milam, Conor Nixon, Statia Luszcz-Cook, Colin Wilson, Vincent Kofman, Giuliano Liuzzi, Sara Faggi, Thomas Fauchez, Manuela Lippi, Richard Cosentino, Alexander Thelen, Arielle Moullet, Paul Hartogh, Edward Molter, Steve Charnley, Giada Arney, Avi Mandell, Nicolas Biver, Ann Vandaele , et al. (2 additional authors not shown)

Abstract: The detection of phosphine (PH3) in the atmosphere of Venus has been recently reported based on millimeter-wave radio observations (Greaves et al. 2020), and its re-analyses (Greaves et al. 2021a/b). In this Matters Arising we perform an independent reanalysis, identifying several issues in the interpretation of the spectroscopic data. As a result, we determine sensitive upper-limits for PH3 in Ve… ▽ More The detection of phosphine (PH3) in the atmosphere of Venus has been recently reported based on millimeter-wave radio observations (Greaves et al. 2020), and its re-analyses (Greaves et al. 2021a/b). In this Matters Arising we perform an independent reanalysis, identifying several issues in the interpretation of the spectroscopic data. As a result, we determine sensitive upper-limits for PH3 in Venus' atmosphere (>75 km, above the cloud decks) that are discrepant with the findings in G2020 and G2021a/b. The measurements target the fundamental first rotational transition of PH3 (J=1-0) at 266.944513 GHz, which was observed with the James Clerk Maxwell Telescope (JCMT) in June 2017 and with the Atacama Large Millimeter/submillimeter Array (ALMA) in March 2019. This line's center is near the SO2 (J=309,21-318,24) transition at 266.943329 GHz (only 1.3 km/s away from the PH3 line) which represents a potential source of contamination. The JCMT and ALMA data, as presented in G2020, are at spectral resolutions comparable to the frequency separation of the two lines. Moreover, the spectral features identified are several km/s in width, and therefore do not permit distinct spectroscopic separation of the candidate spectral lines of PH3 and SO2. We present the radiative transfer modelling we have performed and then discuss the ALMA and JCMT analyses in turn. △ Less

Submitted 21 July, 2021; v1 submitted 27 October, 2020; originally announced October 2020.

Journal ref: Nature Astronomy 5, 631-635 (2021)

Physical properties of the ambient medium and of dense cores in the Perseus star-forming region derived from Herschel Gould Belt Survey observations

Authors: S. Pezzuto, M. Benedettini, J. Di Francesco, P. Palmeirim, S. Sadavoy, E. Schisano, G. Li Causi, Ph. André, D. Arzoumanian, J. -Ph. Bernard, S. Bontemps, D. Elia, E. Fiorellino, J. M. Kirk, V. Könyves, B. Ladjelate, A. Menshchikov, F. Motte, L. Piccotti, N. Schneider, L. Spinoglio, D. Ward-Thompson, C. D. Wilson

Abstract: (Abridged) In this paper, we present analyses of images taken with the Herschel ESA satellite from 70mu to 500mu. We first constructed column density and dust temperature maps. Next, we identified compact cores in the maps, and characterize the cores using modified blackbody fits to their SEDs: we identified 684 starless cores, of which 199 are bound and potential prestellar cores, and 132 protost… ▽ More (Abridged) In this paper, we present analyses of images taken with the Herschel ESA satellite from 70mu to 500mu. We first constructed column density and dust temperature maps. Next, we identified compact cores in the maps, and characterize the cores using modified blackbody fits to their SEDs: we identified 684 starless cores, of which 199 are bound and potential prestellar cores, and 132 protostars. We also matched the Herschel-identified young stars with GAIA sources to model distance variations across the Perseus cloud. We measure a linear gradient function with right ascension and declination for the entire cloud. From the SED fits, mass and temperature of cores were derived. The core mass function can be modelled with a log-normal distribution that peaks at 0.82~$M_\sun$ suggesting a star formation efficiency of 0.30. The high-mass tail can be modelled with a power law of slope $\sim-2.32$, close to the Salpeter's value. We also identify the filamentary structure of Perseus, confirming that stars form preferentially in filaments. We find that the majority of filaments where star formation is ongoing are transcritical against their own internal gravity because their linear masses are below the critical limit of 16~$M_\sun$pc$^{-1}$ above which we expect filaments to collapse. We find a possible explanation for this result, showing that a filament with a linear mass as low as 8~$M_\sun$pc$^{-1}$ can be already unstable. We confirm a linear relation between star formation efficiency and slope of dust probability density function and a similar relation is also seen with the core formation efficiency. We derive a lifetime for the prestellar core phase of $1.69\pm0.52$~Myr for Perseus but different regions have a wide range in prestellar core fractions, hint that star-formation has started only recently in some clumps. We also derive a free-fall time for prestellar cores of 0.16~Myr. △ Less

Submitted 4 November, 2020; v1 submitted 30 September, 2020; originally announced October 2020.

Comments: Comparison with Zari et al's (2016) work improved once Eleonora Zari told us that, contrarily to what we wrote, their data are publicly available (this was not written in the first version of their paper put on arxiv, then our misunderstaning). The modification is small, a paragraph on page 11, and does not have any impact on the content of the paper. Modification allowed by A&A editor

Journal ref: A&A 645, A55 (2021)

A new estimator of resolved molecular gas in nearby galaxies

Authors: Ryan Chown, Cheng Li, Laura C. Parker, Christine D. Wilson, Niu Li, Yang Gao

Abstract: A relationship between dust-reprocessed light from recent star formation and the amount of star-forming gas in a galaxy produces a correlation between WISE 12 $μ$m emission and CO line emission. Here we explore this correlation on kiloparsec scales with CO(1-0) maps from EDGE-CALIFA matched in resolution to WISE 12 $μ$m images. We find strong CO-12 $μ$m correlations within each galaxy and we show… ▽ More A relationship between dust-reprocessed light from recent star formation and the amount of star-forming gas in a galaxy produces a correlation between WISE 12 $μ$m emission and CO line emission. Here we explore this correlation on kiloparsec scales with CO(1-0) maps from EDGE-CALIFA matched in resolution to WISE 12 $μ$m images. We find strong CO-12 $μ$m correlations within each galaxy and we show that the scatter in the global CO-12 $μ$m correlation is largely driven by differences from galaxy to galaxy. The correlation is stronger than that between star formation rate and H$_2$ surface densities ($Σ(\mathrm{H_2})$). We explore multi-variable regression to predict $Σ(\mathrm{H_2})$ in star-forming pixels using the WISE 12 $μ$m data combined with global and resolved galaxy properties, and provide the fit parameters for the best estimators. We find that $Σ(\mathrm{H_2})$ estimators that include $Σ(\mathrm{12\>μm})$ are able to predict $Σ(\mathrm{H_2})$ more accurately than estimators that include resolved optical properties instead of $Σ(\mathrm{12\>μm})$. These results suggest that 12 $μ$m emission and H$_2$ as traced by CO emission are physically connected at kiloparsec scales. This may be due to a connection between polycyclic aromatic hydrocarbon (PAH) emission and the presence of H$_2$. The best single-property estimator is $\log \frac{Σ(\mathrm{H_2})}{\mathrm{M_\odot\>pc^{-2}}} = (0.48 \pm 0.01) + (0.71 \pm 0.01)\log \frac{Σ(\mathrm{12\>μm})}{\mathrm{L_\odot\>pc^{-2}}}$. This correlation can be used to efficiently estimate $Σ(\mathrm{H_2})$ down to at least $1 \> M_\odot \> \mathrm{pc^{-2}}$ in star-forming regions within nearby galaxies. △ Less

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

Comments: 19 pages, 8 figures, 6 tables, accepted by MNRAS. Data tables and a script for making 2D cutouts for each galaxy are available as MNRAS supplementary material and at https://lig.astro.tsinghua.edu.cn/astrodata/publication/20/estimator-resolved-molecular/

Is this an Early Stage Merger? A Case Study on Molecular Gas and Star Formation Properties of Arp 240

Authors: Hao He, C. D. Wilson, Kazimierz Sliwa, Daisuke Iono, Toshiki Saito

Abstract: We present new high resolution $^{12}$CO $J$=1-0, $J$=2-1, and $^{13}$CO $J$=1-0 maps of the early stage merger Arp 240 (NGC5257/8) obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). Simulations in the literature suggest that the merger has just completed its first passage; however, we find that this system has a lower global gas fraction but a higher star formation efficiency… ▽ More We present new high resolution $^{12}$CO $J$=1-0, $J$=2-1, and $^{13}$CO $J$=1-0 maps of the early stage merger Arp 240 (NGC5257/8) obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). Simulations in the literature suggest that the merger has just completed its first passage; however, we find that this system has a lower global gas fraction but a higher star formation efficiency compared to typical close galaxy pairs, which suggests that this system may already be in an advanced merger stage. We combine the ALMA data with $^{12}$CO $J$=3-2 observations from the Submillimeter Array and carry out RADEX modeling on several different regions. Both the RADEX modeling and a local thermal equilibrium (LTE) analysis show that the regions are most likely to have a CO-to-H$_2$ conversion factor $α_{\mathrm{CO}}$ close to or perhaps even smaller than the typical value for (ultra-)luminous infrared galaxies. Using 33 GHz data from the Very Large Array to measure the star formation rate, we find that most star forming regions have molecular gas depletion times of less than 100 Myr. We calculated the star formation efficiency (SFE) per free-fall time for different regions and find some regions appear to have values greater than 100%. We find these regions generally show evidence for young massive clusters (YMCs). After exploring various factors, we argue that this is mainly due to the fact that radio continuum emission in those regions is dominated by that from YMCs, which results in an overestimate of the SFE per free-fall time. △ Less

Submitted 22 June, 2020; originally announced June 2020.

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