Showing 1–50 of 452 results for author: Thompson, A

DarkSide-20k sensitivity to light dark matter particles

Authors: DarkSide-20k Collaboration, :, F. Acerbi, P. Adhikari, P. Agnes, I. Ahmad, S. Albergo, I. F. M. Albuquerque, T. Alexander, A. K. Alton, P. Amaudruz, M. Angiolilli, E. Aprile, R. Ardito, M. Atzori Corona, D. J. Auty, M. Ave, I. C. Avetisov, O. Azzolini, H. O. Back, Z. Balmforth, A. Barrado Olmedo, P. Barrillon, G. Batignani, P. Bhowmick , et al. (289 additional authors not shown)

Abstract: The dual-phase liquid argon time projection chamber is presently one of the leading technologies to search for dark matter particles with masses below 10 GeV/c$^2$. This was demonstrated by the DarkSide-50 experiment with approximately 50 kg of low-radioactivity liquid argon as target material. The next generation experiment DarkSide-20k, currently under construction, will use 1,000 times more arg… ▽ More The dual-phase liquid argon time projection chamber is presently one of the leading technologies to search for dark matter particles with masses below 10 GeV/c$^2$. This was demonstrated by the DarkSide-50 experiment with approximately 50 kg of low-radioactivity liquid argon as target material. The next generation experiment DarkSide-20k, currently under construction, will use 1,000 times more argon and is expected to start operation in 2027. Based on the DarkSide-50 experience, here we assess the DarkSide-20k sensitivity to models predicting light dark matter particles, including Weakly Interacting Massive Particles (WIMPs) and sub-GeV/c$^2$ particles interacting with electrons in argon atoms. With one year of data, a sensitivity improvement to dark matter interaction cross-sections by at least one order of magnitude with respect to DarkSide-50 is expected for all these models. A sensitivity to WIMP--nucleon interaction cross-sections below $1\times10^{-42}$ cm$^2$ is achievable for WIMP masses above 800 MeV/c$^2$. With 10 years exposure, the neutrino fog can be reached for WIMP masses around 5 GeV/c$^2$. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: submitted to Nature Communications

Theory and Observation of Winds from Star-Forming Galaxies

Authors: Todd A. Thompson, Timothy M. Heckman

Abstract: Galactic winds shape the stellar, gas, and metal content of galaxies. To quantify their impact, we must understand their physics. We review potential wind-driving mechanisms and observed wind properties, with a focus on the warm ionized and hot X-ray-emitting gas. Energy and momentum injection by supernovae (SNe), cosmic rays, radiation pressure, and magnetic fields are considered in the light of… ▽ More Galactic winds shape the stellar, gas, and metal content of galaxies. To quantify their impact, we must understand their physics. We review potential wind-driving mechanisms and observed wind properties, with a focus on the warm ionized and hot X-ray-emitting gas. Energy and momentum injection by supernovae (SNe), cosmic rays, radiation pressure, and magnetic fields are considered in the light of observations: (1) Emission and absorption line measurements of cool/warm gas provide our best physical diagnostics of galactic outflows. (2) The critical unsolved problem is how to accelerate cool gas to the high velocities observed. Although conclusive evidence for no one mechanism exists, the momentum, energy, and mass-loading budgets observed compare well with theory. (3) A model where star formation provides a force $\sim L/c$, where $L$ is the bolometric luminosity, and cool gas is pushed out of the galaxy's gravitational potential, compares well with available data. The wind power is $\sim0.1$ that provided by SNe. (4) The very hot X-ray emitting phase, may be a (or the) prime mover. Momentum and energy exchange between the hot and cooler phases is critical to the gas dynamics. (5) Gaps in our observational knowledge include the hot gas kinematics and the size and structure of the outflows probed with UV absorption lines. Simulations are needed to more fully understand mixing, cloud-radiation, cloud-cosmic ray, and cloud-hot wind interactions, the collective effects of star clusters, and both distributed and clustered SNe. Observational works should seek secondary correlations in the wind data that provide evidence for specific mechanisms and compare spectroscopy with the column density-velocity results from theory. △ Less

Submitted 20 June, 2024; v1 submitted 12 June, 2024; originally announced June 2024.

Comments: 68 pages, 11 figures, 2 tables. Authors' version. To appear in Annual Reviews of Astronomy and Astrophysics 2024, Vol 62. Comments welcome

A study of Galactic Plane Planck Galactic Cold Clumps observed by SCOPE and the JCMT Plane Survey

Authors: D. J. Eden, Tie Liu, T. J. T. Moore, J. Di Francesco, G. Fuller, Kee-Tae Kim, Di Li, S. -Y. Liu, R. Plume, Ken'ichi Tatematsu, M. A. Thompson, Y. Wu, L. Bronfman, H. M. Butner, M. J. Currie, G. Garay, P. F. Goldsmith, N. Hirano, D. Johnstone, M. Juvela, S. -P. Lai, C. W. Lee, E. E. Mannfors, F. Olguin, K. Pattle , et al. (10 additional authors not shown)

Abstract: We have investigated the physical properties of Planck Galactic Cold Clumps (PGCCs) located in the Galactic Plane, using the JCMT Plane Survey (JPS) and the SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. By utilising a suite of molecular-line surveys, velocities and distances were assigned to the compact sources within the PGCCs, placing them in a Galactic context. Th… ▽ More We have investigated the physical properties of Planck Galactic Cold Clumps (PGCCs) located in the Galactic Plane, using the JCMT Plane Survey (JPS) and the SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. By utilising a suite of molecular-line surveys, velocities and distances were assigned to the compact sources within the PGCCs, placing them in a Galactic context. The properties of these compact sources show no large-scale variations with Galactic environment. Investigating the star-forming content of the sample, we find that the luminosity-to-mass ratio (L/M) is an order of magnitude lower than in other Galactic studies, indicating that these objects are hosting lower levels of star formation. Finally, by comparing ATLASGAL sources that are associated or are not associated with PGCCs, we find that those associated with PGCCs are typically colder, denser, and have a lower L/M ratio, hinting that PGCCs are a distinct population of Galactic Plane sources. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 18 pages, 14 figures, 7 tables. Accepted for publication in MNRAS

Gamma-ray Emission from a Young Star Cluster in the Star-Forming Region RCW 38

Authors: Paarmita Pandey, Laura A. Lopez, Anna L. Rosen, Todd A. Thompson, Tim Linden, Ian Blackstone, Lachlan Lancaster

Abstract: We report the detection of gamma-ray emission near the young Milky Way star cluster ($\approx$ 0.5 Myr old) in the star-forming region RCW 38. Using 15 years of data from the Fermi-LAT, we find a significant ($ σ> 22$) detection coincident with the cluster, producing a total $γ$-ray luminosity of $L_γ = (2.66\pm 0.92) \times 10^{34}$ erg s$^{-1}$ adopting a power-law spectral model (… ▽ More We report the detection of gamma-ray emission near the young Milky Way star cluster ($\approx$ 0.5 Myr old) in the star-forming region RCW 38. Using 15 years of data from the Fermi-LAT, we find a significant ($ σ> 22$) detection coincident with the cluster, producing a total $γ$-ray luminosity of $L_γ = (2.66\pm 0.92) \times 10^{34}$ erg s$^{-1}$ adopting a power-law spectral model ($Γ= 2.34\pm0.04$) in the 0.1 $-$ 500 GeV band. We estimate the total wind power to be $7 \times 10^{36}$ erg s$^{-1}$, corresponding to a CR acceleration efficiency of $η_{\rm CR} \simeq 0.4$ for a diffusion coefficient consistent with the local interstellar medium of $D = 10^{28}$ cm$^{2}$ s$^{-1}$. Alternatively, the $γ$-ray luminosity could also account for a lower acceleration efficiency of 0.1 if the diffusion coefficient in the star-forming region is smaller, $D\simeq 2.5\times10^{27}\,{\rm cm^2\,\,s^{-1}}$. We compare the hot-gas pressure from Chandra X-ray analysis to the CR pressure and find the former is four orders of magnitude greater, suggesting that the CR pressure is not dynamically important relative to the stellar wind feedback. As RCW 38 is too young for supernovae to have occurred, the high CR acceleration efficiency in RCW 38 demonstrates that stellar winds may be an important source of Galactic cosmic rays. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ

A Population Analysis of 20 Exoplanets Observed from the Optical to the Near-infrared Wavelengths with HST: Evidence for Widespread Stellar Contamination

Authors: Arianna Saba, Alexandra Thompson, Kai Hou Yip, Sushuang Ma, Angelos Tsiaras, Ahmed Faris Al-Refaie, Giovanna Tinetti

Abstract: We present a population study of 20 exoplanets, ranging from Neptune-like to inflated hot-Jupiter planets, observed during transit with the STIS and WFC3 instruments aboard the Hubble Space Telescope. To obtain spectral information from the near-UV to the near-infrared, we reanalysed sixteen WFC3 and over fifty STIS archival data sets with our dedicated HST pipeline. We also include twenty-four WF… ▽ More We present a population study of 20 exoplanets, ranging from Neptune-like to inflated hot-Jupiter planets, observed during transit with the STIS and WFC3 instruments aboard the Hubble Space Telescope. To obtain spectral information from the near-UV to the near-infrared, we reanalysed sixteen WFC3 and over fifty STIS archival data sets with our dedicated HST pipeline. We also include twenty-four WFC3 data sets previously reduced with the same software. Across our target sample we observe significant divergence among multiple observations conducted with the same STIS grating at various epochs, whilst we do not detect variations in the WFC3 data sets. These results are suggestive of stellar contamination, which we have investigated further using known Bayesian tools and other tailored metrics, facilitating a more objective assessment of stellar activity intensity within each system. Our findings reveal that stellar activity contaminates up to half of the studied exoplanet atmospheres, albeit at varying extents. Accounting for stellar activity can significantly alter planetary atmospheric parameters like molecular abundances (up to 6 orders of magnitude) and temperature (up to 145 %), contrasting with the results of analyses that neglect activity. Our results emphasise the importance of considering the effects of stellar contamination in exoplanet transit studies; this issue is particularly true for data sets obtained with facilities that do not cover the optical and/or UV spectral range where the activity is expected to be more impactful but also more easily detectable. Our results also provide a catalogue of potentially active stars for further investigation and monitoring. △ Less

Submitted 25 April, 2024; v1 submitted 23 April, 2024; originally announced April 2024.

Comments: 68 pages, 18 figures, 25 tables, submitted to ApJS

H-alpha emission and HII regions at the locations of recent supernovae in nearby galaxies

Authors: Ness Mayker Chen, Adam K. Leroy, Sumit K. Sarbadhicary, Laura A. Lopez, Todd A. Thompson, Ashley T. Barnes, Eric Emsellem, Brent Groves, Rupali Chandar, Mélanie Chevance, Ryan Chown, Daniel A. Dale, Oleg V. Egorov, Simon C. O. Glover, Kathryn Grasha, Ralf S. Klessen, Kathryn Kreckel, Jing Li, J. Eduardo Méndez-Delgado, Eric J. Murphy, Debosmita Pathak, Eva Schinnerer, David A. Thilker, Leonardo Úbeda, Thomas G. Williams

Abstract: We present a statistical analysis of the local, approximately 50-100 pc scale, H-alpha emission at the locations of recent (less than 125 years) supernovae (SNe) in nearby star-forming galaxies. Our sample consists of 32 SNe in 10 galaxies that are targets of the PHANGS-MUSE survey. We find that 41% (13/32) of these SNe occur coincident with a previously identified HII region. For comparison, HII… ▽ More We present a statistical analysis of the local, approximately 50-100 pc scale, H-alpha emission at the locations of recent (less than 125 years) supernovae (SNe) in nearby star-forming galaxies. Our sample consists of 32 SNe in 10 galaxies that are targets of the PHANGS-MUSE survey. We find that 41% (13/32) of these SNe occur coincident with a previously identified HII region. For comparison, HII regions cover 32% of the area within 1 kpc of any recent SN. Contrasting this local covering fraction with the fraction of SNe coincident with HII regions, we find a statistical excess of 7.6% +/- 8.7% of all SNe to be associated with HII regions. This increases to an excess of 19.2% +/- 10.4% when considering only core-collapse SNe. These estimates appear to be in good agreement with qualitative results from new, higher resolution HST H-alpha imaging, which also suggest many CCSNe detonate near but not in HII regions. Our results appear consistent with the expectation that only a modest fraction of stars explode during the first 5 Myr of the life of a stellar population, when H-alpha emission is expected to be bright. Of the HII region associated SNe, 8% (11/13) also have associated detected CO(2-1) emission, indicating the presence of molecular gas. The HII region associated SNe have typical Av extinctions approximately equal to 1 mag, consistent with a significant amount of pre-clearing of gas from the region before the SNe explode. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ; 33 pages, 13 figures, 3 tables in two-column AASTEX63 format

Predictions for CO emission and the CO-to-H$_2$ conversion factor in galaxy simulations with non-equilibrium chemistry

Authors: Oliver A. Thompson, Alexander J. Richings, Brad K. Gibson, Claude-André Faucher-Giguère, Robert Feldmann, Christopher C. Hayward

Abstract: Our ability to trace the star-forming molecular gas is important to our understanding of the Universe. We can trace this gas using CO emission, converting the observed CO intensity into the H$_2$ gas mass of the region using the CO-to-H$_2$ conversion factor (Xco). In this paper, we use simulations to study the conversion factor and the molecular gas within galaxies. We analysed a suite of simulat… ▽ More Our ability to trace the star-forming molecular gas is important to our understanding of the Universe. We can trace this gas using CO emission, converting the observed CO intensity into the H$_2$ gas mass of the region using the CO-to-H$_2$ conversion factor (Xco). In this paper, we use simulations to study the conversion factor and the molecular gas within galaxies. We analysed a suite of simulations of isolated disc galaxies, ranging from dwarfs to Milky Way-mass galaxies, that were run using the FIRE-2 subgrid models coupled to the CHIMES non-equilibrium chemistry solver. We use the non-equilibrium abundances from the simulations, and we also compare to results using abundances assuming equilibrium, which we calculate from the simulation in post-processing. Our non-equilibrium simulations are able to reproduce the relation between CO and H$_2$ column densities, and the relation between Xco and metallicity, seen within observations of the Milky Way. We also compare to the xCOLD GASS survey, and find agreement with their data to our predicted CO luminosities at fixed star formation rate. We also find the multivariate function used by xCOLD GASS overpredicts the H$_2$ mass for our simulations, motivating us to suggest an alternative multivariate function of our fitting, though we caution that this fitting is uncertain due to the limited range of galaxy conditions covered by our simulations. We also find that the non-equilibrium chemistry has little effect on the conversion factor (<5\%) for our high-mass galaxies, though still affects the H$_2$ mass and Lco by $\approx$25\%. △ Less

Submitted 12 June, 2024; v1 submitted 21 March, 2024; originally announced March 2024.

Comments: 19 pages, 16 figures (including appendices). Accepted for publication in MNRAS

Atacama Large Aperture Submillimeter Telescope (AtLAST) Science: Our Galaxy

Authors: Pamela Klaassen, Alessio Traficante, Maria T. Beltrán, Kate Pattle, Mark Booth, Joshua B. Lovell, Jonathan P. Marshall, Alvaro Hacar, Brandt A. L. Gaches, Caroline Bot, Nicolas Peretto, Thomas Stanke, Doris Arzoumanian, Ana Duarte Cabral, Gaspard Duchêne, David J. Eden, Antonio Hales, Jens Kauffmann, Patricia Luppe, Sebastian Marino, Elena Redaelli, Andrew J. Rigby, Álvaro Sánchez-Monge, Eugenio Schisano, Dmitry A. Semenov , et al. (16 additional authors not shown)

Abstract: As we learn more about the multi-scale interstellar medium (ISM) of our Galaxy, we develop a greater understanding for the complex relationships between the large-scale diffuse gas and dust in Giant Molecular Clouds (GMCs), how it moves, how it is affected by the nearby massive stars, and which portions of those GMCs eventually collapse into star forming regions. The complex interactions of those… ▽ More As we learn more about the multi-scale interstellar medium (ISM) of our Galaxy, we develop a greater understanding for the complex relationships between the large-scale diffuse gas and dust in Giant Molecular Clouds (GMCs), how it moves, how it is affected by the nearby massive stars, and which portions of those GMCs eventually collapse into star forming regions. The complex interactions of those gas, dust and stellar populations form what has come to be known as the ecology of our Galaxy. Because we are deeply embedded in the plane of our Galaxy, it takes up a significant fraction of the sky, with complex dust lanes scattered throughout the optically recognisable bands of the Milky Way. These bands become bright at (sub-)millimetre wavelengths, where we can study dust thermal emission and the chemical and kinematic signatures of the gas. To properly study such large-scale environments, requires deep, large area surveys that are not possible with current facilities. Moreover, where stars form, so too do planetary systems, growing from the dust and gas in circumstellar discs, to planets and planetesimal belts. Understanding the evolution of these belts requires deep imaging capable of studying belts around young stellar objects to Kuiper belt analogues around the nearest stars. Here we present a plan for observing the Galactic Plane and circumstellar environments to quantify the physical structure, the magnetic fields, the dynamics, chemistry, star formation, and planetary system evolution of the galaxy in which we live with AtLAST; a concept for a new, 50m single-dish sub-mm telescope with a large field of view which is the only type of facility that will allow us to observe our Galaxy deeply and widely enough to make a leap forward in our understanding of our local ecology. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 27 pages, submitted to Open Research Europe as part of the AtLAST collection: https://open-research-europe.ec.europa.eu/collections/atlast/about

Journal ref: Open Res Europe 2024, 4:112 [version 1; peer review: awaiting peer review]

SDHDF: A new file format for spectral-domain radio astronomy data

Authors: L. J. Toomey, G. Hobbs, D. C. Price, J. R. Dawson, T. Wenger, D. Lagoy, L. Staveley-Smith, J. A. Green, E. Carretti, A. Hafner, M. Huynh, J. Kaczmarek, S. Mader, V. McIntyre, J. Reynolds, T. Robishaw, J. Sarkissian, A. Thompson, C. Tremblay, A. Zic

Abstract: Radio astronomy file formats are now required to store wide frequency bandwidths and multiple simultaneous receiver beams and must be able to account for versatile observing modes and numerous calibration strategies. The need to capture and archive high-time and high frequency-resolution data, along with the comprehensive metadata that fully describe the data, implies that a new data format and ne… ▽ More Radio astronomy file formats are now required to store wide frequency bandwidths and multiple simultaneous receiver beams and must be able to account for versatile observing modes and numerous calibration strategies. The need to capture and archive high-time and high frequency-resolution data, along with the comprehensive metadata that fully describe the data, implies that a new data format and new processing software are required. This requirement is suited to a well-defined, hierarchically-structured and flexible file format. In this paper we present the Spectral-Domain Hierarchical Data Format (`SDHDF') -- a new file format for radio astronomy data, in particular for single dish or beam-formed data streams. Since 2018, SDHDF has been the primary format for data products from the spectral-line and continuum observing modes at Murriyang, the CSIRO Parkes 64-m radio telescope, and we demonstrate that this data format can also be used to store observations of pulsars and fast radio bursts. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: Supplementary material (SDHDF definition): https://ars.els-cdn.com/content/image/1-s2.0-S2213133724000192-mmc1.pdf

Journal ref: Astronomy and Computing, Volume 47 (April 2024)

Favorable conditions for heavy element nucleosynthesis in rotating proto-magnetar winds

Authors: Tejas Prasanna, Matthew S. B. Coleman, Todd A. Thompson

Abstract: The neutrino-driven wind cooling phase of proto-neutron stars (PNSs) follows successful supernovae. Wind models without magnetic fields or rotation fail to achieve the necessary conditions for production of the third $r$-process peak, but robustly produce a weak $r$-process in neutron-rich winds. Using 2D magnetohydrodynamic simulations with magnetar-strength magnetic fields and rotation, we show… ▽ More The neutrino-driven wind cooling phase of proto-neutron stars (PNSs) follows successful supernovae. Wind models without magnetic fields or rotation fail to achieve the necessary conditions for production of the third $r$-process peak, but robustly produce a weak $r$-process in neutron-rich winds. Using 2D magnetohydrodynamic simulations with magnetar-strength magnetic fields and rotation, we show that the PNS rotation rate significantly affects the thermodynamic conditions of the wind. We show that high entropy material is quasi-periodically ejected from the closed zone of the PNS magnetosphere with the required thermodynamic conditions to produce heavy elements. We show that maximum entropy $S$ of the material ejected depends systematically on the magnetar spin period $P_{\star}$ and scales as $S \propto P_{\star}^{-5/6}$ for sufficiently rapid rotation. We present results from simulations at a constant neutrino luminosity representative of $\sim 1-2$ s after the onset of cooling for $P_{\star}$ ranging from 5 ms to 200 ms and a few simulations with evolving neutrino luminosity where we follow the evolution of the magnetar wind until $10-14$ s after the onset of cooling. We estimate at magnetar polar magnetic field strength $B_0=3\times 10^{15}$ G, $10^{15}$ G, and $5\times 10^{14}$ G that neutron-rich magnetar winds can respectively produce at least $\sim 2-5\times 10^{-5}$ M$_{\odot}$, $\sim 3-4\times 10^{-6}$ M$_{\odot}$, and $\sim 2.5\times 10^{-8}$ M$_{\odot}$ of material with the required parameters for synthesis of the third $r-$process peak, within $1-2$ s, 10 s, and 14 s in that order after the onset of cooling. We show that proton-rich magnetar winds can have favorable conditions for production of $p-$nuclei, even at a modest $B_0=5\times 10^{14}$ G. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 21 pages, 14 figures

JWST Observations of Starbursts: Polycyclic Aromatic Hydrocarbon Emission at the Base of the M 82 Galactic Wind

Authors: Alberto D. Bolatto, Rebecca C. Levy, Elizabeth Tarantino, Martha L. Boyer, Deanne B. Fisher, Adam K. Leroy, Serena A. Cronin, Ralf S. Klessen, J. D. Smith, Dannielle A. Berg, Torsten Boeker, Leindert A. Boogaard, Eve C. Ostriker, Todd A. Thompson, Juergen Ott, Laura Lenkic, Laura A. Lopez, Daniel A. Dale, Sylvain Veilleux, Paul P. van der Werf, Simon C. O. Glover, Karin M. Sandstrom, Evan D. Skillman, John Chisholm, Vicente Villanueva , et al. (15 additional authors not shown)

Abstract: We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find promine… ▽ More We present new observations of the central 1 kpc of the M 82 starburst obtained with the James Webb Space Telescope (JWST) near-infrared camera (NIRCam) instrument at a resolution ~0.05"-0.1" (~1-2 pc). The data comprises images in three mostly continuum filters (F140M, F250M, and F360M), and filters that contain [FeII] (F164N), H2 v=1-0 (F212N), and the 3.3 um PAH feature (F335M). We find prominent plumes of PAH emission extending outward from the central starburst region, together with a network of complex filamentary substructure and edge-brightened bubble-like features. The structure of the PAH emission closely resembles that of the ionized gas, as revealed in Paschen alpha and free-free radio emission. We discuss the origin of the structure, and suggest the PAHs are embedded in a combination of neutral, molecular, and photoionized gas. △ Less

Submitted 21 April, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted to The Astrophysical Journal

High mass function ellipsoidal variables in the Gaia Focused Product Release: searching for black hole candidates in the binary zoo

Authors: D. M. Rowan, Todd A. Thompson, T. Jayasinghe, C. S. Kochanek, K. Z. Stanek

Abstract: The recent Gaia Focused Product Release contains radial velocity time-series for more than 9,000 Gaia long-period photometric variables. Here we search for binary systems with large radial velocity amplitudes to identify candidates with massive, unseen companions. Eight targets have binary mass function $f(M)>1\ M_\odot$, three of which are eclipsing binaries. The remaining five show evidence of e… ▽ More The recent Gaia Focused Product Release contains radial velocity time-series for more than 9,000 Gaia long-period photometric variables. Here we search for binary systems with large radial velocity amplitudes to identify candidates with massive, unseen companions. Eight targets have binary mass function $f(M)>1\ M_\odot$, three of which are eclipsing binaries. The remaining five show evidence of ellipsoidal modulations. We fit spectroscopic orbit models to the Gaia radial velocities, and fit the spectral energy distributions of three targets. For the two systems most likely to host dark companions, J0946 and J1640, we use PHOEBE to fit the ASAS-SN light curves and Gaia radial velocities. The derived companion masses are $>3 M_\odot$, but the high Galactic dust extinctions towards these objects limit our ability to rule out main sequence companions or subgiants hotter than the photometric primaries. These systems are similar to other stellar-mass black hole impostors, notably the Unicorn (V723 Mon) and the Giraffe (2M04123153$+$6738486). While it is possible that J1640 and J0946 are similar examples of stripped giant star binaries, high-resolution spectra can be used to determine the nature of their companions. △ Less

Submitted 22 March, 2024; v1 submitted 17 January, 2024; originally announced January 2024.

Comments: 18 pages, 13 figures. Accepted to The Open Journal of Astrophysics

The dynamic centres of infrared-dark clouds and the formation of cores

Authors: Andrew J. Rigby, Nicolas Peretto, Michael Anderson, Sarah E. Ragan, Felix D. Priestley, Gary A. Fuller, Mark A. Thompson, Alessio Traficante, Elizabeth J. Watkins, Gwenllian M. Williams

Abstract: High-mass stars have an enormous influence on the evolution of the interstellar medium in galaxies, so it is important that we understand how they form. We examine the central clumps within a sample of seven infrared-dark clouds (IRDCs) with a range of masses and morphologies. We use 1 pc-scale observations from NOEMA and the IRAM 30-m telescope to trace dense cores with 2.8 mm continuum, and gas… ▽ More High-mass stars have an enormous influence on the evolution of the interstellar medium in galaxies, so it is important that we understand how they form. We examine the central clumps within a sample of seven infrared-dark clouds (IRDCs) with a range of masses and morphologies. We use 1 pc-scale observations from NOEMA and the IRAM 30-m telescope to trace dense cores with 2.8 mm continuum, and gas kinematics in C$^{18}$O, HCO$^+$, HNC, and N$_2$H$^+$ ($J$=1$-$0). We supplement our continuum sample with six IRDCs observed at 2.9 mm with ALMA, and examine the relationships between core- and clump-scale properties. We have developed a fully-automated multiple-velocity component hyperfine line-fitting code called mwydyn which we employ to trace the dense gas kinematics in N$_2$H$^+$ (1$-$0), revealing highly complex and dynamic clump interiors. We find that parsec-scale clump mass is the most important factor driving the evolution; more massive clumps are able to concentrate more mass into their most massive cores - with a log-normally distributed efficiency of around 9% - in addition to containing the most dynamic gas. Distributions of linewidths within the most massive cores are similar to the ambient gas, suggesting that they are not dynamically decoupled, but are similarly chaotic. A number of studies have previously suggested that clumps are globally collapsing; in such a scenario, the observed kinematics of clump centres would be the direct result of gravity-driven mass inflows that become ever more complex as the clumps evolve, which in turn leads to the chaotic mass growth of their core populations. △ Less

Submitted 31 January, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: Accepted by MNRAS

Hot Gas Outflow Properties of the Starburst Galaxy NGC 4945

Authors: Natalia Porraz Barrera, Sebastian Lopez, Laura A. Lopez, Adi Foord, Dustin D. Nguyen, Todd A. Thompson, Smita Mathur, Alberto D. Bolatto

Abstract: We analyze 330 ks of {\it Chandra} X-ray imaging and spectra of the nearby, edge-on starburst and Seyfert Type 2 galaxy NGC 4945 to measure the hot gas properties along the galactic outflows. We extract and model spectra from 15 regions extending from $-$0.55 kpc to $+$0.85 kpc above and below the galactic disk to determine the best-fit parameters and metal abundances. We find that the hot gas tem… ▽ More We analyze 330 ks of {\it Chandra} X-ray imaging and spectra of the nearby, edge-on starburst and Seyfert Type 2 galaxy NGC 4945 to measure the hot gas properties along the galactic outflows. We extract and model spectra from 15 regions extending from $-$0.55 kpc to $+$0.85 kpc above and below the galactic disk to determine the best-fit parameters and metal abundances. We find that the hot gas temperatures and number densities peak in the central regions and decrease along the outflows. These profiles are inconsistent with a spherical, adiabatically-expanding wind model, suggesting the need to include mass loading and/or a non-spherical outflow geometry. We estimate the mass outflow rate of the hot wind to be $1.6\:M_{\odot}~\rm{yr}^{-1}$. Emission from charge exchange is detected in the northern outflow, and we estimate it contributes 12\% to the emitted, broad-band ($0.5-7$ keV) X-ray flux. △ Less

Submitted 11 April, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

Comments: 11 pages, 7 figures; submitted (13 December 2023), accepted ApJ (11 April 2023)

The SARAO MeerKAT 1.3 GHz Galactic Plane Survey

Authors: S. Goedhart, W. D. Cotton, F. Camilo, M. A. Thompson, G. Umana, M. Bietenholz, P. A. Woudt, L. D. Anderson, C. Bordiu, D. A. H. Buckley, C. S. Buemi, F. Bufano, F. Cavallaro, H. Chen, J. O. Chibueze, D. Egbo, B. S. Frank, M. G. Hoare, A. Ingallinera, T. Irabor, R. C. Kraan-Korteweg, S. Kurapati, P. Leto, S. Loru, M. Mutale , et al. (105 additional authors not shown)

Abstract: We present the SARAO MeerKAT Galactic Plane Survey (SMGPS), a 1.3 GHz continuum survey of almost half of the Galactic Plane (251°$\le l \le$ 358°and 2°$\le l \le$ 61°at $|b| \le 1.5°$). SMGPS is the largest, most sensitive and highest angular resolution 1 GHz survey of the Plane yet carried out, with an angular resolution of 8" and a broadband RMS sensitivity of $\sim$10--20 $μ$ Jy/beam. Here we d… ▽ More We present the SARAO MeerKAT Galactic Plane Survey (SMGPS), a 1.3 GHz continuum survey of almost half of the Galactic Plane (251°$\le l \le$ 358°and 2°$\le l \le$ 61°at $|b| \le 1.5°$). SMGPS is the largest, most sensitive and highest angular resolution 1 GHz survey of the Plane yet carried out, with an angular resolution of 8" and a broadband RMS sensitivity of $\sim$10--20 $μ$ Jy/beam. Here we describe the first publicly available data release from SMGPS which comprises data cubes of frequency-resolved images over 908--1656 MHz, power law fits to the images, and broadband zeroth moment integrated intensity images. A thorough assessment of the data quality and guidance for future usage of the data products are given. Finally, we discuss the tremendous potential of SMGPS by showcasing highlights of the Galactic and extragalactic science that it permits. These highlights include the discovery of a new population of non-thermal radio filaments; identification of new candidate supernova remnants, pulsar wind nebulae and planetary nebulae; improved radio/mid-IR classification of rare Luminous Blue Variables and discovery of associated extended radio nebulae; new radio stars identified by Bayesian cross-matching techniques; the realisation that many of the largest radio-quiet WISE HII region candidates are not true HII regions; and a large sample of previously undiscovered background HI galaxies in the Zone of Avoidance. △ Less

Submitted 2 May, 2024; v1 submitted 12 December, 2023; originally announced December 2023.

Comments: Accepted for publication in MNRAS. The data release is live and links can be found in the Data Availability Statement in the paper

A Two-Component Probability Distribution Function Describes the mid-IR Emission from the Disks of Star-Forming Galaxies

Authors: Debosmita Pathak, Adam K. Leroy, Todd A. Thompson, Laura A. Lopez, Francesco Belfiore, Mederic Boquien, Daniel A. Dale, Simon C. O. Glover, Ralf S. Klessen, Eric W. Koch, Erik Rosolowsky, Karin M. Sandstrom, Eva Schinnerer, Rowan Smith, Jiayi Sun, Jessica Sutter, Thomas G. Williams, Frank Bigiel, Yixian Cao, Jeremy Chastenet, Melanie Chevance, Ryan Chown, Eric Emsellem, Christopher M. Faesi, Kirsten L. Larson , et al. (6 additional authors not shown)

Abstract: High-resolution JWST-MIRI images of nearby spiral galaxies reveal emission with complex substructures that trace dust heated both by massive young stars and the diffuse interstellar radiation field. We present high angular (0."85) and physical resolution (20-80 pc) measurements of the probability distribution function (PDF) of mid-infrared (mid-IR) emission (7.7-21 $μ$m) from 19 nearby star-formin… ▽ More High-resolution JWST-MIRI images of nearby spiral galaxies reveal emission with complex substructures that trace dust heated both by massive young stars and the diffuse interstellar radiation field. We present high angular (0."85) and physical resolution (20-80 pc) measurements of the probability distribution function (PDF) of mid-infrared (mid-IR) emission (7.7-21 $μ$m) from 19 nearby star-forming galaxies from the PHANGS-JWST Cycle-1 Treasury. The PDFs of mid-IR emission from the disks of all 19 galaxies consistently show two distinct components: an approximately log-normal distribution at lower intensities and a high-intensity power-law component. These two components only emerge once individual star-forming regions are resolved. Comparing with locations of HII regions identified from VLT/MUSE H$α$-mapping, we infer that the power-law component arises from star-forming regions and thus primarily traces dust heated by young stars. In the continuum-dominated 21 $μ$m band, the power-law is more prominent and contains roughly half of the total flux. At 7.7-11.3 $μ$m, the power-law is suppressed by the destruction of small grains (including PAHs) close to HII regions while the log-normal component tracing the dust column in diffuse regions appears more prominent. The width and shape of the log-normal diffuse emission PDFs in galactic disks remain consistent across our sample, implying a log-normal gas column density $N$(H)$\approx10^{21}$cm$^{-2}$ shaped by supersonic turbulence with typical (isothermal) turbulent Mach numbers $\approx5-15$. Finally, we describe how the PDFs of galactic disks are assembled from dusty HII regions and diffuse gas, and discuss how the measured PDF parameters correlate with global properties such as star-formation rate and gas surface density. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: 30 pages without appendix, 17 figures, (with appendix images of full sample: 56 pages, 39 figures), accepted in AJ

Dynamics of baryon ejection in magnetar giant flares: implications for radio afterglows, r-process nucleosynthesis, and fast radio bursts

Authors: Jakub Cehula, Todd A. Thompson, Brian D. Metzger

Abstract: We explore the impact of a magnetar giant flare (GF) on the neutron star (NS) crust, and the associated baryon mass ejection. We consider that sudden magnetic energy dissipation creates a thin high-pressure shell above a portion of the NS surface, which drives a relativistic shockwave into the crust, heating a fraction of these layers sufficiently to become unbound along directions unconfined by t… ▽ More We explore the impact of a magnetar giant flare (GF) on the neutron star (NS) crust, and the associated baryon mass ejection. We consider that sudden magnetic energy dissipation creates a thin high-pressure shell above a portion of the NS surface, which drives a relativistic shockwave into the crust, heating a fraction of these layers sufficiently to become unbound along directions unconfined by the magnetic field. We explore this process using spherically-symmetric relativistic hydrodynamical simulations. For an initial shell pressure $P_{\rm GF}$ we find the total unbound ejecta mass roughly obeys the relation $M_{\rm{ej}}\sim4-9\times10^{24}\:\rm{g}\:(P_{\rm GF}/10^{30}\:\rm{ergs}\:\rm{cm}^{-3})^{1.43}$. For $P_{\rm{GF}}\sim10^{30}-10^{31}\:\rm{ergs}\:\rm{cm}^{-3}$ corresponding to the dissipation of a magnetic field of strength $\sim10^{15.5}-10^{16}\:\rm{G}$, we find $M_{\rm{ej}}\sim10^{25}-10^{26}\:\rm{g}$ with asymptotic velocities $v_{\rm{ej}}/c\sim0.3-0.6$ compatible with the ejecta properties inferred from the afterglow of the December 2004 GF from SGR 1806-20. Because the flare excavates crustal material to a depth characterized by an electron fraction $Y_e\approx0.40-0.46$, and is ejected with high entropy and rapid expansion timescale, the conditions are met for heavy element $r$-process nucleosynthesis via the alpha-rich freeze-out mechanism. Given an energetic GF rate of roughly once per century in the Milky Way, we find that magnetar GFs could be an appreciable heavy $r$-process source that tracks star formation. We predict that GFs are accompanied by short $\sim$minutes long, luminous $\sim10^{39}\:\rm{ergs}\:\rm{s}^{-1}$ optical transients powered by $r$-process decay ("nova brevis"), akin to scaled-down kilonovae. Our findings also have implications for the synchrotron nebulae surrounding some repeating fast radio burst sources. △ Less

Submitted 1 February, 2024; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 24 pages, 12 figures, 2 tables. v2: matches the accepted version in MNRAS (only minor edits were made). Videos showing our Fiducial simulation are available at https://dx.doi.org/10.5281/zenodo.10593900 and at https://www.youtube.com/playlist?list=PLtfS1xXqUmtC2A48QhRvGXxbeMK1vkNnC

Neural ODEs as a discovery tool to characterize the structure of the hot galactic wind of M82

Authors: Dustin D. Nguyen, Yuan-Sen Ting, Todd A. Thompson, Sebastian Lopez, Laura A. Lopez

Abstract: Dynamic astrophysical phenomena are predominantly described by differential equations, yet our understanding of these systems is constrained by our incomplete grasp of non-linear physics and scarcity of comprehensive datasets. As such, advancing techniques in solving non-linear inverse problems becomes pivotal to addressing numerous outstanding questions in the field. In particular, modeling hot g… ▽ More Dynamic astrophysical phenomena are predominantly described by differential equations, yet our understanding of these systems is constrained by our incomplete grasp of non-linear physics and scarcity of comprehensive datasets. As such, advancing techniques in solving non-linear inverse problems becomes pivotal to addressing numerous outstanding questions in the field. In particular, modeling hot galactic winds is difficult because of unknown structure for various physical terms, and the lack of \textit{any} kinematic observational data. Additionally, the flow equations contain singularities that lead to numerical instability, making parameter sweeps non-trivial. We leverage differentiable programming, which enables neural networks to be embedded as individual terms within the governing coupled ordinary differential equations (ODEs), and show that this method can adeptly learn hidden physics. We robustly discern the structure of a mass-loading function which captures the physical effects of cloud destruction and entrainment into the hot superwind. Within a supervised learning framework, we formulate our loss function anchored on the astrophysical entropy ($K \propto P/ρ^{5/3}$). Our results demonstrate the efficacy of this approach, even in the absence of kinematic data $v$. We then apply these models to real Chandra X-Ray observations of starburst galaxy M82, providing the first systematic description of mass-loading within the superwind. This work further highlights neural ODEs as a useful discovery tool with mechanistic interpretability in non-linear inverse problems. We make our code public at this GitHub repository (https://github.com/dustindnguyen/2023_NeurIPS_NeuralODEs_M82). △ Less

Submitted 28 November, 2023; v1 submitted 3 November, 2023; originally announced November 2023.

Comments: 9 Pages, 2 Figures, Accepted at the NeurIPS 2023 workshop on Machine Learning and the Physical Sciences

Prospects for detecting proto-neutron star rotation and spindown using supernova neutrinos

Authors: Tejas Prasanna, Todd A. Thompson, Christopher Hirata

Abstract: After a successful supernova, a proto-neutron star (PNS) cools by emitting neutrinos on $\sim 1-100$ s timescales. Provided that there are neutrino emission `hot-spots' or `cold-spots' on the surface of the rotating PNS, we can expect a periodic modulation in the number of neutrinos observable by detectors. We show that Fourier transform techniques can be used to determine the PNS rotation rate fr… ▽ More After a successful supernova, a proto-neutron star (PNS) cools by emitting neutrinos on $\sim 1-100$ s timescales. Provided that there are neutrino emission `hot-spots' or `cold-spots' on the surface of the rotating PNS, we can expect a periodic modulation in the number of neutrinos observable by detectors. We show that Fourier transform techniques can be used to determine the PNS rotation rate from the neutrino arrival times. Provided there is no spindown, a 1-parameter Discrete Fourier Transform (DFT) is sufficient to determine the spin period of the PNS. If the PNS is born as a magnetar with polar magnetic field strength $B_0 \gtrsim 10^{15}$ G and is `slowly' rotating with an initial spin period $\gtrsim 100$ ms, then it can spindown to periods of the order of seconds during the cooling phase. We propose a modified DFT technique with three frequency parameters to detect spindown. Due to lack of neutrino data from a nearby supernova except the $\sim20$ neutrinos detected from SN1987A, we use toy models and one physically motivated modulating function to generate neutrino arrival times. We use the false alarm rate (FAR) to quantify the significance of the Fourier power spectrum peaks. We show that PNS rotation and spindown are detected with $\rm FAR<2\%$ ($2σ$) for periodic signal content $\rm M\gtrsim 13-15\%$ if $5\times 10^{3}$ neutrinos are detected in $\sim 3$ s and with $\rm FAR<1\%$ for $\rm M\geq 5\%$ if $5\times 10^{4}$ neutrinos are detected in $\sim 3$ s. Since we can expect $\sim 10^{4}-10^{5}$ neutrino detections from a supernova at 10 kpc, detection of PNS rotation and spindown is possible using the neutrinos from the next Galactic supernova. △ Less

Submitted 15 February, 2024; v1 submitted 20 October, 2023; originally announced October 2023.

Comments: 18 pages, 16 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 528, Issue 4, March 2024, Pages 5649-5666

CONGRuENTS (COsmic-ray, Neutrino, Gamma-ray and Radio Non-Thermal Spectra). II. Population-level correlations between galactic infrared, radio, and γ-ray emission

Authors: Matt A. Roth, Mark R. Krumholz, Roland M. Crocker, Todd A. Thompson

Abstract: Galaxies obey a number of empirical correlations between their radio, γ-ray, and infrared emission, but the physical origins of these correlations remain uncertain. Here we use the CONGRuENTS model for broadband non-thermal emission from star-forming galaxies, which self-consistently calculates energy-dependent transport and non-thermal emission from cosmic ray hadrons and leptons, to predict radi… ▽ More Galaxies obey a number of empirical correlations between their radio, γ-ray, and infrared emission, but the physical origins of these correlations remain uncertain. Here we use the CONGRuENTS model for broadband non-thermal emission from star-forming galaxies, which self-consistently calculates energy-dependent transport and non-thermal emission from cosmic ray hadrons and leptons, to predict radio and γ-ray emission for a synthetic galaxy population with properties drawn from a large deep-field survey. We show that our synthetic galaxies reproduce observed relations such as the FIR-radio correlation, the FIR-γ correlation, and the distribution of radio spectral indices, and we use the model to explain the physical origins of these relations. Our results show that the FIR-radio correlation arises because the amount of cosmic ray electron power ultimately radiated as synchrotron emission varies only weakly with galaxy star formation rate as a result of the constraints imposed on gas properties by hydrostatic balance and turbulent dynamo action; the same physics dictates the extent of proton calorimetry in different galaxies, and thus sets the FIR-γ-ray correlation. We further show that galactic radio spectral indices result primarily from competition between thermal free-free emission and energy-dependent loss of cosmic ray electrons to bremsstrahlung and escape into galactic halos, with shaping of the spectrum by inverse Compton, synchrotron, and ionisation processes typically playing a sub-dominant role. In addition to explaining existing observations, we use our analysis to predict a heretofore unseen correlation between the curvature of galaxies' radio spectra and their pion-driven γ-ray emission, a prediction that will be testable with upcoming facilities. △ Less

Submitted 15 June, 2024; v1 submitted 9 October, 2023; originally announced October 2023.

Comments: 17 pages, 14 figures

Journal ref: MNRAS, Volume 530, Issue 2, May 2024, Pages 1849-1865

Outgassing Composition of the Murchison Meteorite: Implications for Volatile Depletion of Planetesimals and Interior-atmosphere Connections for Terrestrial Exoplanets

Authors: Maggie A. Thompson, Myriam Telus, Graham Harper Edwards, Laura Schaefer, Jasmeet Dhaliwal, Brian Dreyer, Jonathan J. Fortney, Kyle Kim

Abstract: Outgassing is a central process during the formation and evolution of terrestrial planets and their atmospheres both within and beyond the solar system. Although terrestrial planets' early atmospheres likely form via outgassing during planetary accretion, the connection between a planet's bulk composition and its initial atmospheric properties is not well understood. One way to inform this connect… ▽ More Outgassing is a central process during the formation and evolution of terrestrial planets and their atmospheres both within and beyond the solar system. Although terrestrial planets' early atmospheres likely form via outgassing during planetary accretion, the connection between a planet's bulk composition and its initial atmospheric properties is not well understood. One way to inform this connection is to analyze the outgassing compositions of meteorites, and in particular carbonaceous chondrites, because they are some of the most volatile-rich, primitive materials (in terms of their bulk compositions) that are available for direct study. In addition, they may serve as compositional analogs for the building block materials of terrestrial planets in our solar system and around other Sun-like stars. This study builds upon previous outgassing experiments that monitored the abundances of volatile species (e.g., H2O, CO, and CO2) released from the Murchison meteorite. To gain a more complete understanding of Murchison's outgassing composition, we perform a series of heating experiments under atmospheric pressure (1 bar) and vacuum (1E-9 bar) conditions on samples of the Murchison meteorite and subsequent bulk element analysis to inform the outgassing trends of a suite of major elements in Murchison (e.g., Fe, Mg, Zn, and S). Under both pressure conditions, sulfur outgases significantly at the highest temperatures (800C - 1000C). For the samples heated under vacuum conditions, we also detect outgassing of zinc. Combined with prior outgassing experiments, this study provides important insights into the volatile depletion patterns of undifferentiated planetesimals and the early outgassing compositions of terrestrial exoplanets. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 23 pages, 11 figures, 15 tables

Journal ref: Planet. Sci. J. 4 185 (2023)

The Scale of Stellar Yields: Implications of the Measured Mean Iron Yield of Core Collapse Supernovae

Authors: David H. Weinberg, Emily J. Griffith, James W. Johnson, Todd A. Thompson

Abstract: The scale of alpha-element yields is difficult to predict from theory because of uncertainties in massive star evolution, supernova physics, and black hole formation, and it is difficult to constrain empirically because the impact of higher yields can be compensated by greater metal loss in galactic winds. We use a recent measurement of the mean iron yield of core collapse supernovae (CCSN) by Rod… ▽ More The scale of alpha-element yields is difficult to predict from theory because of uncertainties in massive star evolution, supernova physics, and black hole formation, and it is difficult to constrain empirically because the impact of higher yields can be compensated by greater metal loss in galactic winds. We use a recent measurement of the mean iron yield of core collapse supernovae (CCSN) by Rodriguez et al. (RMN23), $\bar{y}_{\rm Fe}^{\rm cc} =0.058 \pm 0.007 M_\odot$, to infer the scale of alpha-element yields by assuming that the plateau of [alpha/Fe] abundance ratios observed in low metallicity stars represents the yield ratio of CCSN. For a Kroupa IMF and a plateau at [alpha/Fe]=0.45, we find that the population-averaged yields of O and Mg per unit mass of star formation are about equal to the mass fractions of these elements in the sun. The inferred O and Fe yields agree with predictions of the Sukhbold et al. (2016) CCSN models assuming their Z9.6+N20 neutrino-driven engine, a scenario in which many progenitors with $M<40M_\odot$ implode to black holes rather than exploding. The yields are lower than assumed in some models of galactic chemical evolution (GCE) and the galaxy mass-metallicity relation, reducing the level of outflows needed to match observed abundances. For straightforward assumptions, we find that one-zone GCE models with mass-loading factor $η\approx 0.6$ evolve to solar metallicity at late times. By requiring that models reach [alpha/Fe]=0 at late times, and assuming a mean Fe yield of $0.7M_\odot$ per Type Ia supernova, we infer a Hubble-time integrated SNIa rate of $1.1\times 10^{-3} M_\odot^{-1}$, compatible with estimates from supernova surveys. The RMN23 measurement provides one of the few empirical anchors for the absolute scale of nucleosynthetic yields, with wide-ranging implications for stellar and galactic astrophysics. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 19 pp., 9 figs, submitted to AAS journals

Highly-mass-loaded hot galactic winds are unstable to cool filament formation

Authors: Dustin D. Nguyen, Todd A. Thompson, Evan E. Schneider, Ashley P. Tarrant

Abstract: When cool clouds are ram-pressure accelerated by a hot supersonic galactic wind, some of the clouds may be shredded by hydrodynamical instabilities and incorporated into the hot flow. Recent one-dimensional steady-state calculations show how cool cloud entrainment directly affects the bulk thermodynamics, kinematics, and observational characteristics of the hot gas. In particular, mass-loading dec… ▽ More When cool clouds are ram-pressure accelerated by a hot supersonic galactic wind, some of the clouds may be shredded by hydrodynamical instabilities and incorporated into the hot flow. Recent one-dimensional steady-state calculations show how cool cloud entrainment directly affects the bulk thermodynamics, kinematics, and observational characteristics of the hot gas. In particular, mass-loading decelerates the hot flow and changes its entropy. Here, we investigate the stability of planar and spherical mass-loaded hot supersonic flows using both perturbation analysis and three-dimensional time-dependent radiative hydrodynamical simulations. We show that mass-loading is stable over a broad range of parameters and that the 1D time-steady analytic solutions exactly reproduce the 3D time-dependent calculations, provided that the flow does not decelerate sufficiently to become subsonic. For higher values of the mass-loading, the flow develops a sonic point and becomes thermally unstable, rapidly cooling and forming elongated dense cometary filaments. We explore the mass-loading parameters required to reach a sonic point and the radiative formation of these filaments. For certain approximations, we can derive simple analytic criteria. In general a mass-loading rate similar to the initial mass outflow rate is required. In this sense, the destruction of small cool clouds by a hot flow may ultimately spontaneously generate fast cool filaments, as observed in starburst superwinds. Lastly, we find that the kinematics of filaments is sensitive to the slope of the mass-loading function. Filaments move faster than the surrounding wind if mass-loading is over long distances whereas filaments move slower than their surroundings if mass-loading is abrupt. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 12 pages, 15 figures, submitted to MNRAS (21 July 2023)

A hidden population of massive white dwarfs: two spotted K+WD binaries

Authors: D. M. Rowan, T. Jayasinghe, M. A. Tucker, C. Y. Lam, Todd A. Thompson, C. S. Kochanek, N. S. Abrams, B. J. Fulton, I. Ilyin, H. Isaacson, J. Lu, D. V. Martin, B. Nicholson

Abstract: The identification and characterization of massive ($\gtrsim 0.8~M_\odot$) white dwarfs is challenging in part due to their low luminosity. Here we present two candidate single-lined spectroscopic binaries, Gaia DR3 4014708864481651840 and 5811237403155163520, with K-dwarf primaries and optically dark companions. Both have orbital periods of $P\sim 0.45$ days and show rotational variability, ellip… ▽ More The identification and characterization of massive ($\gtrsim 0.8~M_\odot$) white dwarfs is challenging in part due to their low luminosity. Here we present two candidate single-lined spectroscopic binaries, Gaia DR3 4014708864481651840 and 5811237403155163520, with K-dwarf primaries and optically dark companions. Both have orbital periods of $P\sim 0.45$ days and show rotational variability, ellipsoidal modulations, and high-amplitude radial velocity variations. Using light curves from the Transiting Exoplanet Survey Satellite (TESS), radial velocities from ground-based spectrographs, and spectral energy distributions, we characterize these binaries to describe the nature of the unseen companion. We find that both systems are consistent with a massive white dwarf companion. Unlike simple ellipsoidal variables, star spots cause the light curve morphology to change between TESS sectors. We attempt to constrain the orbital inclination using PHOEBE binary light curve models, but degeneracies in the light curves of spotted stars prevent a precise determination. Finally, we search for similar objects using Gaia DR3 and TESS, and comment on these systems in the context of recently claimed compact object binaries. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 17 pages, 15 figures. Submitted to MNRAS

SCOTCH -- Search for Clandestine Optically Thick Compact HIIs

Authors: A. L. Patel, J. S. Urquhart, A. Y. Yang, T. J. T Moore, K. M. Menten, M. A. Thompson, M. G. Hoare, T. Irabor, S. L. Breen, M. D. Smith

Abstract: This study uses archival high frequency continuum data to expand the search for Hypercompact HII regions and determine the conditions at which they appear, as this stage high mass star formation is short-lived and rare. We use 23 GHz continuum data taken towards methanol masers, which are an excellent signpost for very young embedded high-mass protostars. We have searched for high-frequency, optic… ▽ More This study uses archival high frequency continuum data to expand the search for Hypercompact HII regions and determine the conditions at which they appear, as this stage high mass star formation is short-lived and rare. We use 23 GHz continuum data taken towards methanol masers, which are an excellent signpost for very young embedded high-mass protostars. We have searched for high-frequency, optically thick radio sources to identify HC HII region candidates. The data cover 128 fields that include 141 methanol masers identified by the Methanol Multibeam (MMB) survey. We have detected 68 high-frequency radio sources and conducted a multi-wavelength analysis to determine their nature. This has identified 49 HII regions, 47 of which are embedded in dense clumps fourteen of which do not have a 5 GHz radio counterpart. We have identified 13 methanol maser sites that are coincident with radio sources that have a steep positive spectral index. The majority of these are not detected in the mid-infrared and have been classified as protostellar or young stellar objects in the literature and we therefore consider to be good HC HII region candidates, however, further work and higher resolution data are needed to confirm these candidates. △ Less

Submitted 16 July, 2023; originally announced July 2023.

Characterising a World Within the Hot Neptune Desert: Transit Observations of LTT 9779 b with HST WFC3

Authors: Billy Edwards, Quentin Changeat, Angelos Tsiaras, Andrew Allan, Patrick Behr, Simone R. Hagey, Michael D. Himes, Sushuang Ma, Keivan G. Stassun, Luis Thomas, Alexandra Thompson, Aaron Boley, Luke Booth, Jeroen Bouwman, Kevin France, Nataliea Lowson, Annabella Meech, Caprice L. Phillips, Aline A. Vidotto, Kai Hou Yip, Michelle Bieger, Amelie Gressier, Estelle Janin, Ing-Guey Jiang, Pietro Leonardi , et al. (5 additional authors not shown)

Abstract: We present an atmospheric analysis of LTT 9779 b, a rare planet situated in the hot Neptune desert, that has been observed with HST WFC3 G102 and G141. The combined transmission spectrum, which covers 0.8 - 1.6 $μ$m, shows a gradual increase in transit depth with wavelength. Our preferred atmospheric model shows evidence for H$_{\rm 2}$O, CO$_{\rm 2}$ and FeH with a significance of 3.1 $σ$, 2.4… ▽ More We present an atmospheric analysis of LTT 9779 b, a rare planet situated in the hot Neptune desert, that has been observed with HST WFC3 G102 and G141. The combined transmission spectrum, which covers 0.8 - 1.6 $μ$m, shows a gradual increase in transit depth with wavelength. Our preferred atmospheric model shows evidence for H$_{\rm 2}$O, CO$_{\rm 2}$ and FeH with a significance of 3.1 $σ$, 2.4 $σ$ and 2.1 $σ$, respectively. In an attempt to constrain the rate of atmospheric escape for this planet, we search for the 1.083 $μ$m Helium line in the G102 data but find no evidence of excess absorption that would indicate an escaping atmosphere using this tracer. We refine the orbital ephemerides of LTT 9779 b using our HST data and observations from TESS, searching for evidence of orbital decay or apsidal precession, which is not found. The phase-curve observation of LTT 9779 b with JWST NIRISS should provide deeper insights into the atmosphere of this planet and the expected atmospheric escape might be detected with further observations concentrated on other tracers such as Lyman $α$. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Comments: Accepted for publication in AJ

Supernova Rates and Luminosity Functions from ASAS-SN I: 2014--2017 Type Ia SNe and Their Subtypes

Authors: D. D. Desai, C. S. Kochanek, B. J. Shappee, T. Jayasinghe, K. Z. Stanek, T. W. -S. Holoien, T. A. Thompson, C. Ashall, J. F. Beacom, A. Do, S. Dong, J. L. Prieto

Abstract: We present the volumetric rates and luminosity functions (LFs) of Type Ia supernovae (SNe Ia) from the $V$-band All-Sky Automated Survey for Supernovae (ASAS-SN) catalogues spanning discovery dates from UTC 2014-01-26 to UTC 2017-12-29. Our standard sample consists of 404 SNe Ia with $m_{V,\mathrm{peak}}<17$ mag and Galactic latitude $|b|>15^{\circ}$. Our results are both statistically more precis… ▽ More We present the volumetric rates and luminosity functions (LFs) of Type Ia supernovae (SNe Ia) from the $V$-band All-Sky Automated Survey for Supernovae (ASAS-SN) catalogues spanning discovery dates from UTC 2014-01-26 to UTC 2017-12-29. Our standard sample consists of 404 SNe Ia with $m_{V,\mathrm{peak}}<17$ mag and Galactic latitude $|b|>15^{\circ}$. Our results are both statistically more precise and systematically more robust than previous studies due to the large sample size and high spectroscopic completeness. We make completeness corrections based on both the apparent and absolute magnitudes by simulating the detection of SNe Ia in ASAS-SN light curves. We find a total volumetric rate for all subtypes of $R_{\mathrm{tot}}=2.28^{+0.20}_{-0.20}\,\times 10^{4}\,\mathrm{yr}^{-1}\,\mathrm{Gpc}^{-3}\,h^{3}_{70}$ for $M_{V,\mathrm{peak}}<-16.5$ mag ($R_{\mathrm{tot}}=1.91^{+0.12}_{-0.12}\,\times 10^{4}\,\mathrm{yr}^{-1}\,\mathrm{Gpc}^{-3}\,h^{3}_{70}$ for $M_{V,\mathrm{peak}}<-17.5$ mag) at the median redshift of our sample, $z_{\mathrm{med}}=0.024$. This is in agreement ($1σ$) with the local volumetric rates found by previous studies. We also compile luminosity functions (LFs) for the entire sample as well as for subtypes of SNe Ia for the first time. The major subtypes with more than one SN include Ia-91bg, Ia-91T, Ia-CSM, and Ia-03fg with total rates of $R_{\mathrm{Ia-91bg}}=1.4^{+0.5}_{-0.5} \times 10^{3}\,\mathrm{yr}^{-1}\,\mathrm{Gpc}^{-3}\,h^{3}_{70}$, $R_{\mathrm{Ia-91T}}=8.5^{+1.6}_{-1.7} \times 10^{2}\,\mathrm{yr}^{-1}\,\mathrm{Gpc}^{-3}\,h^{3}_{70}$, $R_{\mathrm{Ia-CSM}}=10^{+7}_{-7}\,\mathrm{yr}^{-1}\,\mathrm{Gpc}^{-3}\,h^{3}_{70}$, and $R_{\mathrm{Ia-03fg}}=30^{+20}_{-20}\,\mathrm{yr}^{-1}\,\mathrm{Gpc}^{-3}\,h^{3}_{70}$, respectively. We estimate a mean host extinction of $E(V-r)\approx 0.2$ mag based on the shift between our $V$-band and the ZTF $r$-band LFs. △ Less

Submitted 23 May, 2024; v1 submitted 19 June, 2023; originally announced June 2023.

Comments: 14 pages, 9 figures, 8 tables, accepted to MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 530, Issue 4, Pages 5016-5029, June 2024

The early evolution of magnetar rotation -- II. Rapidly rotating magnetars: Implications for Gamma-Ray Bursts and Super Luminous Supernovae

Authors: Tejas Prasanna, Matthew S. B. Coleman, Matthias J. Raives, Todd A. Thompson

Abstract: Rapidly rotating magnetars have been associated with gamma-ray bursts (GRBs) and super-luminous supernovae (SLSNe). Using a suite of 2D magnetohydrodynamic simulations at fixed neutrino luminosity and a couple of evolutionary models with evolving neutrino luminosity and magnetar spin period, we show that magnetars are viable central engines for powering GRBs and SLSNe. We also present analytic est… ▽ More Rapidly rotating magnetars have been associated with gamma-ray bursts (GRBs) and super-luminous supernovae (SLSNe). Using a suite of 2D magnetohydrodynamic simulations at fixed neutrino luminosity and a couple of evolutionary models with evolving neutrino luminosity and magnetar spin period, we show that magnetars are viable central engines for powering GRBs and SLSNe. We also present analytic estimates of the energy outflow rate from the proto-neutron star (PNS) as a function of polar magnetic field strength $B_0$, PNS angular velocity $Ω_{\star}$, PNS radius $R_{\star}$ and mass outflow rate $\dot{M}$. We show that rapidly rotating magnetars with spin periods $P_{\star}\lesssim 4$ ms and polar magnetic field strength $B_0\gtrsim 10^{15}$ G can release $10^{50}-5\times 10^{51}$ ergs of energy during the first $\sim2$ s of the cooling phase. Based on this result, it is plausible that sustained energy injection by magnetars through the relativistic wind phase can power GRBs. We also show that magnetars with moderate field strengths of $B_0\lesssim 5\times 10^{14}$ G do not release a large fraction of their rotational kinetic energy during the cooling phase and hence, are not likely to power GRBs. Although we cannot simulate to times greater than $\sim 3-5$ s after a supernova, we can hypothesize that moderate field strength magnetars can brighten the supernova light curves by releasing their rotational kinetic energy via magnetic dipole radiation on timescales of days to weeks, since these do not expend most of their rotational kinetic energy during the early cooling phase. △ Less

Submitted 21 October, 2023; v1 submitted 25 May, 2023; originally announced May 2023.

Comments: 15 pages, 13 Figures

Journal ref: Monthly Notices of the Royal Astronomical Society 2023, Volume 526, Issue 2, pp.3141-3155

ASAS-SN Sky Patrol V2.0

Authors: K. Hart, B. J. Shappee, D. Hey, C. S. Kochanek, K. Z. Stanek, L. Lim, S. Dobbs, M. Tucker, T. Jayasinghe, J. F. Beacom, T. Boright, T. Holoien, J. M. Joel Ong, J. L. Prieto, T. A. Thompson, D. Will

Abstract: The All-Sky Automated Survey for Supernovae (ASAS-SN) began observing in late-2011 and has been imaging the entire sky with nightly cadence since late 2017. A core goal of ASAS-SN is to release as much useful data as possible to the community. Working towards this goal, in 2017 the first ASAS-SN Sky Patrol was established as a tool for the community to obtain light curves from our data with no pre… ▽ More The All-Sky Automated Survey for Supernovae (ASAS-SN) began observing in late-2011 and has been imaging the entire sky with nightly cadence since late 2017. A core goal of ASAS-SN is to release as much useful data as possible to the community. Working towards this goal, in 2017 the first ASAS-SN Sky Patrol was established as a tool for the community to obtain light curves from our data with no preselection of targets. Then, in 2020 we released static V-band photometry from 2013--2018 for 61 million sources. Here we describe the next generation ASAS-SN Sky Patrol, Version 2.0, which represents a major progression of this effort. Sky Patrol 2.0 provides continuously updated light curves for 111 million targets derived from numerous external catalogs of stars, galaxies, and solar system objects. We are generally able to serve photometry data within an hour of observation. Moreover, with a novel database architecture, the catalogs and light curves can be queried at unparalleled speed, returning thousands of light curves within seconds. Light curves can be accessed through a web interface (http://asas-sn.ifa.hawaii.edu/skypatrol/) or a Python client (https://asas-sn.ifa.hawaii.edu/documentation). The Python client can be used to retrieve up to 1 million light curves, generally limited only by bandwidth. This paper gives an updated overview of our survey, introduces the new Sky Patrol, and describes its system architecture. These results provide significant new capabilities to the community for pursuing multi-messenger and time-domain astronomy. △ Less

Submitted 7 April, 2023; originally announced April 2023.

Comments: Light curves can be accessed through a web interface http://asas-sn.ifa.hawaii.edu/skypatrol, or a Python client at http://asas-sn.ifa.hawaii.edu/documentation

Dust Eddington Ratios for Star-Forming Galaxy Subregions

Authors: Ian Blackstone, Todd A. Thompson

Abstract: Radiation pressure on dust is an important feedback process around star clusters and may eject gas from bright sub-regions in star-forming galaxies. The Eddington ratio has previously been constructed for galaxy-averaged observations, individual star clusters, and Galactic HII regions. Here we assess the role of radiation pressure in thousands of sub-regions across two local star-forming galaxies,… ▽ More Radiation pressure on dust is an important feedback process around star clusters and may eject gas from bright sub-regions in star-forming galaxies. The Eddington ratio has previously been constructed for galaxy-averaged observations, individual star clusters, and Galactic HII regions. Here we assess the role of radiation pressure in thousands of sub-regions across two local star-forming galaxies, NGC 6946 and NGC 5194. Using a model for the spectral energy distribution from stellar population synthesis and realistic dust grain scattering and absorption, we compute flux- and radiation pressure-mean opacities and population-averaged optical depth $\langleτ_{\rm RP}\rangle$. Using Monte-Carlo calculations, we assess the momentum coupling through a dusty column to the stellar continuum. Optically-thin regions around young stellar populations are $30-50$ times super-Eddington. We calculate the Eddington ratio for the sub-regions including the local mass of young and old stars and HI and molecular gas. We compute the fraction of the total star formation that is currently super-Eddington, and provide an assessment of the role of radiation pressure in the dusty gas dynamics. Depending on the assumed height of the dusty gas and the age of the stellar population, we find that $\sim0-10$% of the sightlines are super-Eddington. These regions may be accelerated to $\sim5-15$ km/s by radiation pressure alone. Additionally, our results show that for beamed radiation the function $1-\exp(-\langleτ_{\rm RP}\rangle)$ is an excellent approximation to the momentum transfer. Opacities and optical depths are tabulated for SEDs of different stellar ages and for continuous star formation. △ Less

Submitted 20 February, 2023; originally announced February 2023.

Comments: 17 pages, 9 figures

Magnetized Rotating Isothermal Winds

Authors: Matthias J. Raives, Matthew S. B. Coleman, Todd A. Thompson

Abstract: We consider the general problem of a Parker-type non-relativistic isothermal wind from a rotating and magnetic star. Using the magnetohydrodynamics (MHD) code athena++, we construct an array of simulations in the stellar rotation rate $Ω_\ast$ and the isothermal sound speed $c_T$, and calculate the mass, angular momentum, and energy loss rates across this parameter space. We also briefly consider… ▽ More We consider the general problem of a Parker-type non-relativistic isothermal wind from a rotating and magnetic star. Using the magnetohydrodynamics (MHD) code athena++, we construct an array of simulations in the stellar rotation rate $Ω_\ast$ and the isothermal sound speed $c_T$, and calculate the mass, angular momentum, and energy loss rates across this parameter space. We also briefly consider the three dimensional case, with misaligned magnetic and rotation axes. We discuss applications of our results to the spindown of normal stars, highly-irradiated exoplanets, and to nascent highly-magnetic and rapidly-rotating neutron stars born in massive star core collapse.8 △ Less

Submitted 30 October, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

Comments: 19 pages, 13 figures; v2: published version after minor revisions

Journal ref: MNRAS, 526, 3, 4418-4434 (2023)

Correcting Exoplanet Transmission Spectra for Stellar Activity with an Optimised Retrieval Framework

Authors: Alexandra Thompson, Alfredo Biagini, Gianluca Cracchiolo, Antonino Petralia, Quentin Changeat, Arianna Saba, Giuseppe Morello, Mario Morvan, Giuseppina Micela, Giovanna Tinetti

Abstract: The chromatic contamination that arises from photospheric heterogeneities e.g. spots and faculae on the host star presents a significant noise source for exoplanet transmission spectra. If this contamination is not corrected for, it can introduce substantial bias in our analysis of the planetary atmosphere. We utilise two stellar models of differing complexity, StARPA and ASteRA, to explore the bi… ▽ More The chromatic contamination that arises from photospheric heterogeneities e.g. spots and faculae on the host star presents a significant noise source for exoplanet transmission spectra. If this contamination is not corrected for, it can introduce substantial bias in our analysis of the planetary atmosphere. We utilise two stellar models of differing complexity, StARPA and ASteRA, to explore the biases introduced by stellar contamination in retrieval under differing degrees of stellar activity. We use the retrieval framework TauREx3 and a grid of 27 synthetic, spot-contaminated transmission spectra to investigate potential biases and to determine how complex our stellar models must be in order to accurately extract the planetary parameters from transmission spectra. The input observation is generated using the more complex model (StARPA), in which the spot latitude is an additional, fixable parameter. This observation is then fed into a combined stellar-planetary retrieval which contains a simplified stellar model (ASteRA). Our results confirm that the inclusion of stellar activity parameters in retrieval minimises bias under all activity regimes considered. ASteRA performs very well under low to moderate activity conditions, retrieving the planetary parameters with a high degree of accuracy. For the most active cases, characterised by larger, higher temperature contrast spots, some minor residual bias remains due to ASteRA neglecting the interplay between the spot and the limb darkening effect. As a result of this, we find larger errors in retrieved planetary parameters for central spots (0 degrees) and those found close to the limb (60 degrees) than those at intermediate latitudes (30 degrees). △ Less

Submitted 26 October, 2023; v1 submitted 9 February, 2023; originally announced February 2023.

Comments: 34 pages, 20 figures, accepted for publication in ApJ

Cosmic-Ray Driven Galactic Winds from the Warm Interstellar Medium

Authors: Shaunak Modak, Eliot Quataert, Yan-Fei Jiang, Todd A. Thompson

Abstract: We study the properties of cosmic-ray (CR) driven galactic winds from the warm interstellar medium using idealized spherically symmetric time-dependent simulations. The key ingredients in the model are radiative cooling and CR-streaming-mediated heating of the gas. Cooling and CR heating balance near the base of the wind, but this equilibrium is thermally unstable, leading to a multiphase wind wit… ▽ More We study the properties of cosmic-ray (CR) driven galactic winds from the warm interstellar medium using idealized spherically symmetric time-dependent simulations. The key ingredients in the model are radiative cooling and CR-streaming-mediated heating of the gas. Cooling and CR heating balance near the base of the wind, but this equilibrium is thermally unstable, leading to a multiphase wind with large fluctuations in density and temperature. In most of our simulations, the heating eventually overwhelms cooling, leading to a rapid increase in temperature and a thermally-driven wind; the exception to this is in galaxies with the shallowest potentials, which produce nearly isothermal $T \approx 10^4$ K winds driven by CR pressure. Many of the time-averaged wind solutions found here have a remarkable critical point structure, with two critical points. Scaled to real galaxies, we find mass outflow rates $\dot M$ somewhat larger than the observed star formation rate in low mass galaxies, and an approximately "energy-like" scaling $\dot M \propto v_{\rm esc}^{-2}$. The winds accelerate slowly and reach asymptotic wind speeds of only $\sim 0.4 v_{\rm esc}$. The total wind power is $\sim 1\%$ of the power from supernovae, suggesting inefficient preventive CR feedback for the physical conditions modeled here. We predict significant spatially extended emission and absorption lines from $10^4 - 10^{5.5}$ K gas; this may correspond to extraplanar diffuse ionized gas seen in star-forming galaxies. △ Less

Submitted 9 August, 2023; v1 submitted 7 February, 2023; originally announced February 2023.

Comments: Updated to accepted version: 19 pages, 14 figures. Main simulation results in Table 1, Figure 3, and Figure 7. Comments still welcome!

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 524, Issue 4, October 2023, Pages 6374-6391

The Co-Ordinated Radio and Infrared Survey for High-Mass Star Formation. V. The CORNISH-South Survey and Catalogue

Authors: T. Irabor, M. G. Hoare, M. Burton, W. D. Cotton, P. Diamond, S. Dougherty, S. P. Ellingsen, R. Fender, G. A. Fuller, S. Garrington, P. F. Goldsmith, J. Green, A. G. Gunn, J. Jackson, S. Kurtz, S. L. Lumsden, J. Marti, I. McDonald, S. Molinari, T. J. Moore, M. Mutale, T. Muxlow, T. OBrien, R. D. Oudmaijer, R. Paladini , et al. (10 additional authors not shown)

Abstract: We present the first high spatial resolution radio continuum survey of the southern Galactic plane. The CORNISH project has mapped the region defined by $295^{\circ} < l < 350^{\circ}$; $|b| < 1^{\circ}$ at 5.5-GHz, with a resolution of 2.5$^{''}$ (FWHM). As with the CORNISH-North survey, this is designed to primarily provide matching radio data to the Spitzer GLIMPSE survey region. The CORNISH-So… ▽ More We present the first high spatial resolution radio continuum survey of the southern Galactic plane. The CORNISH project has mapped the region defined by $295^{\circ} < l < 350^{\circ}$; $|b| < 1^{\circ}$ at 5.5-GHz, with a resolution of 2.5$^{''}$ (FWHM). As with the CORNISH-North survey, this is designed to primarily provide matching radio data to the Spitzer GLIMPSE survey region. The CORNISH-South survey achieved a root mean square noise level of $\sim$ 0.11 mJy beam$^{-1}$, using the 6A configuration of the Australia Telescope Compact Array (ATCA). In this paper, we discuss the observations, data processing and measurements of the source properties. Above a 7$σ$ detection limit, 4701 sources were detected, and their ensemble properties show similar distributions with their northern counterparts. The catalogue is highly reliable and is complete to 90 per cent at a flux density level of 1.1 mJy. We developed a new way of measuring the integrated flux densities and angular sizes of non-Gaussian sources. The catalogue primarily provides positions, flux density measurements and angular sizes. All sources with IR counterparts at 8$μm$ have been visually classified, utilizing additional imaging data from optical, near-IR, mid-IR, far-IR and sub-millimetre galactic plane surveys. This has resulted in the detection of 524 H II regions of which 255 are ultra-compact H II regions, 287 planetary nebulae, 79 radio stars and 6 massive young stellar objects. The rest of the sources are likely to be extra-galactic. These data are particularly important in the characterization and population studies of compact ionized sources such as UCHII regions and PNe towards the Galactic mid-plane. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comparing the Locations of Supernovae to CO (2-1) Emission in their Host Galaxies

Authors: Ness Mayker Chen, Adam K. Leroy, Laura A. Lopez, Samantha Benincasa, Mélanie Chevance, Simon C. O. Glover, Annie Hughes, Kathryn Kreckel, Sumit Sarbadhicary, Jiayi Sun, Todd A. Thompson, Dyas Utomo, Frank Bigiel, Guillermo A. Blanc, Daniel A. Dale, Kathryn Grasha, J. M. Diederik Kruijssen, Hsi-An Pan, Miguel Querejeta, Eva Schinnerer, Elizabeth J. Watkins, Thomas G. Williams

Abstract: We measure the molecular gas environment near recent ($< 100$ yr old) supernovae (SNe) using $\sim1''$ or $\leq 150$pc resolution CO (2-1) maps from the PHANGS-ALMA survey of nearby star-forming galaxies. This is arguably the first such study to approach the scales of individual massive molecular clouds ($M_{\rm mol} \gtrsim 10^{5.3}$ M$_{\odot}$). Using the Open Supernova Catalog (OSC), we identi… ▽ More We measure the molecular gas environment near recent ($< 100$ yr old) supernovae (SNe) using $\sim1''$ or $\leq 150$pc resolution CO (2-1) maps from the PHANGS-ALMA survey of nearby star-forming galaxies. This is arguably the first such study to approach the scales of individual massive molecular clouds ($M_{\rm mol} \gtrsim 10^{5.3}$ M$_{\odot}$). Using the Open Supernova Catalog (OSC), we identify 63 SNe within the PHANGS-ALMA footprint. We detect CO (2-1) emission near $\sim60\%$ of the sample at 150pc resolution, compared to $\sim35\%$ of map pixels with CO (2-1) emission, and up to $\sim95\%$ of the SNe at 1kpc resolution compared to $\sim80\%$ of map pixels with CO (2-1) emission. We expect the $\sim60\%$ of SNe within the same 150pc beam as a GMC will likely interact with these clouds in the future, consistent with the observation of widespread SN-molecular gas interaction in the Milky Way, while the other $\sim40\%$ of SNe without strong CO (2-1) detections will deposit their energy in the diffuse interstellar medium (ISM), perhaps helping drive large-scale turbulence or galactic outflows. Broken down by type, we detect CO (2-1) emission at the sites of $\sim85\%$ of our 9 stripped-envelope SNe (SESNe), $\sim40\%$ of our 34 Type II SNe, and $\sim35\%$ of our 13 Type Ia SNe, indicating that SESNe are most closely associated with the brightest CO (2-1) emitting regions in our sample. Our results confirm that SN explosions are not restricted to only the densest gas, and instead exert feedback across a wide range of molecular gas densities. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: Accepted to ApJ; 38 pages, 12 figures, 7 tables in two-column AASTEX63 format

CONGRuENTS (COsmic-ray, Neutrino, Gamma-ray and Radio Non-Thermal Spectra). I. A predictive model for galactic non-thermal emission

Authors: Matt A. Roth, Mark R. Krumholz, Roland M. Crocker, Todd A. Thompson

Abstract: The total luminosity and spectral shape of the non-thermal emission produced by cosmic rays depends on their interstellar environment, a dependence that gives rise to correlations between galaxies' bulk properties -- star formation rate, stellar mass, and others -- and their non-thermal spectra. Understanding the physical mechanisms of cosmic ray transport, loss, and emission is key to understandi… ▽ More The total luminosity and spectral shape of the non-thermal emission produced by cosmic rays depends on their interstellar environment, a dependence that gives rise to correlations between galaxies' bulk properties -- star formation rate, stellar mass, and others -- and their non-thermal spectra. Understanding the physical mechanisms of cosmic ray transport, loss, and emission is key to understanding these correlations. Here, in the first paper of the series, we present a new method to compute the non-thermal spectra of star-forming galaxies, and describe an open-source software package -- COsmic-ray, Neutrino, Gamma-ray and Radio Non-Thermal Spectra (CONGRuENTS) -- that implements it. As a crucial innovation, our method requires as input only a galaxy's effective radius, star formation rate, stellar mass, and redshift, all quantities that are readily available for large samples of galaxies and do not require expensive, spatially resolved gas measurements. From these inputs we derive individual, galaxy-by-galaxy models for the background gas and radiation field through which cosmic rays propagate, from which we compute steady state cosmic ray spectra for hadronic and leptonic particles in both the galactic disc and halo by solving the full kinetic equation. We invoke modern models for cosmic ray transport and include all significant emission and loss mechanisms. In this paper we describe the model and validate it against non-thermal emission measured in nearby star-forming galaxies that span four orders of magnitude in star formation rate. △ Less

Submitted 16 May, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 23 pages, 14 figures, 1 table, accepted for publication in MNRAS

The Value-Added Catalog of ASAS-SN Eclipsing Binaries III: Masses and Radii of Gaia Spectroscopic Binaries

Authors: D. M. Rowan, T. Jayasinghe, K. Z. Stanek, C. S. Kochanek, Todd A. Thompson, B. J. Shappee, W. Giles

Abstract: Masses and radii of stars can be derived by combining eclipsing binary light curves with spectroscopic orbits. In our previous work, we modeled the All-Sky Automated Survey for Supernovae (ASAS-SN) light curves of more than 30,000 detached eclipsing binaries using PHOEBE. Here we combine our results with 128 double-lined spectroscopic orbits from Gaia Data Release 3. We visually inspect ASAS-SN li… ▽ More Masses and radii of stars can be derived by combining eclipsing binary light curves with spectroscopic orbits. In our previous work, we modeled the All-Sky Automated Survey for Supernovae (ASAS-SN) light curves of more than 30,000 detached eclipsing binaries using PHOEBE. Here we combine our results with 128 double-lined spectroscopic orbits from Gaia Data Release 3. We visually inspect ASAS-SN light curves of double-lined spectroscopic binaries on the lower main sequence and the giant branch, adding 11 binaries to our sample. We find that only 50% of systems have Gaia periods and eccentricities consistent with the ASAS-SN values. We use emcee and PHOEBE to determine masses and radii for a total of 122 stars with median fractional uncertainties of 7.9% and 6.3%, respectively. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: 11 pages, 9 figures. Submitted to MNRAS. Data available online at https://asas-sn.osu.edu/binaries

Cool, Luminous, and Highly Variable Stars in the Magellanic Clouds. II: Spectroscopic and Environmental Analysis of Thorne-Żytkow Object and Super-AGB Star Candidates

Authors: Anna J. G. O'Grady, Maria R. Drout, B. M. Gaensler, C. S. Kochanek, Kathryn F. Neugent, Carolyn L. Doherty, Joshua S. Speagle, B. J. Shappee, Michael Rauch, Ylva Götberg, Bethany Ludwig, Todd A. Thompson

Abstract: In previous work we identified a population of 38 cool and luminous variable stars in the Magellanic Clouds and examined 11 in detail in order to classify them as either Thorne-Żytkow Objects (TŻOs, red supergiants with a neutron star cores) or super-AGB stars (the most massive stars that will not undergo core collapse). This population includes HV\,2112, a peculiar star previously considered in o… ▽ More In previous work we identified a population of 38 cool and luminous variable stars in the Magellanic Clouds and examined 11 in detail in order to classify them as either Thorne-Żytkow Objects (TŻOs, red supergiants with a neutron star cores) or super-AGB stars (the most massive stars that will not undergo core collapse). This population includes HV\,2112, a peculiar star previously considered in other works to be either a TŻO or high-mass AGB star. Here we continue this investigation, using the kinematic and radio environments and local star formation history of these stars to place constraints on the age of the progenitor systems and the presence of past supernovae. These stars are not associated with regions of recent star formation, and we find no evidence of past supernovae at their locations. Finally, we also assess the presence of heavy elements and lithium in their spectra compared to red supergiants. We find strong absorption in Li and s-process elements compared to RSGs in most of the sample, consistent with super-AGB nucleosynthesis, while HV\,2112 shows additional strong lines associated with TŻO nucleosynthesis. Coupled with our previous mass estimates, the results are consistent with the stars being massive (~4-6.5M$_{\odot}$) or super-AGB (~6.5-12M$_{\odot}$) stars in the thermally pulsing phase, providing crucial observations of the transition between low- and high-mass stellar populations. HV\,2112 is more ambiguous; it could either be a maximally massive sAGB star, or a TŻO if the minimum mass for stability extends down to <13 M$_\odot$. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Comments: 32 pages, 10 figures, 7 tables, spectroscopic data available at https://zenodo.org/record/7058608, accepted to The Astrophysical Journal

Dynamics of hot galactic winds launched from spherically-stratified starburst cores

Authors: Dustin D. Nguyen, Todd A. Thompson, Evan E. Schneider, Sebastian Lopez, Laura A. Lopez

Abstract: The analytic galactic wind model derived by Chevalier and Clegg in 1985 (CC85) assumes $\textit{uniform}$ energy and mass-injection within the starburst galaxy nucleus. However, the structure of nuclear star clusters, bulges, and star-forming knots are non-uniform. We generalize to cases with spherically-symmetric energy/mass injection that scale as $r^{-Δ}$ within the starburst volume $R$, provid… ▽ More The analytic galactic wind model derived by Chevalier and Clegg in 1985 (CC85) assumes $\textit{uniform}$ energy and mass-injection within the starburst galaxy nucleus. However, the structure of nuclear star clusters, bulges, and star-forming knots are non-uniform. We generalize to cases with spherically-symmetric energy/mass injection that scale as $r^{-Δ}$ within the starburst volume $R$, providing solutions for $Δ= 0$, 1/2, 1, 3/2, and 2. In marked contrast with the CC85 model ($Δ=0$), which predicts zero velocity at the center, for a singular isothermal sphere profile ($Δ=2$), we find that the flow maintains a $\textit{constant}$ Mach number of $\mathcal{M}=\sqrt{3/5} \simeq 0.77$ throughout the volume. The fast interior flow can be written as $v_{r < R} = (\dot{E}_T/3\dot{M}_T)^{1/2} \simeq 0.41 \, v_\infty$, where $v_\infty$ is the asymptotic velocity, and $\dot{E}_T$ and $\dot{M}_T$ are the total energy and mass injection rates. For $v_\infty \simeq 2000 \, \mathrm{km \, s^{-1}}$, $v_{r<R} \simeq 820 \, \mathrm{km\, s^{-1}}$ throughout the wind-driving region. The temperature and density profiles of the non-uniform models may be important for interpreting spatially-resolved maps of starburst nuclei. We compute velocity resolved spectra to contrast the $Δ=0$ (CC85) and $Δ=2$ models. Next generation X-ray space telescopes such as XRISM may assess these kinematic predictions. △ Less

Submitted 8 November, 2022; v1 submitted 13 October, 2022; originally announced October 2022.

Comments: submitted (13 October 2022), accepted MNRAS Letters (08 November 2022)

The ASAS-SN Bright Supernova Catalog -- V. 2018-2020

Authors: K. D. Neumann, T. W. -S. Holoien, C. S. Kochanek, K. Z. Stanek, P. J. Vallely, B. J. Shappee, J. L. Prieto, T. Pessi, T. Jayasinghe, J. Brimacombe, D. Bersier, E. Aydi, C. Basinger, J. F. Beacom, S. Bose, J. S. Brown, P. Chen, A. Clocchiatti, D. D. Desai, Subo Dong, E. Falco, S. Holmbo, N. Morrell, J. V. Shields, K. V. Sokolovsky , et al. (33 additional authors not shown)

Abstract: We catalog the 443 bright supernovae discovered by the All-Sky Automated Survey for Supernovae (ASAS-SN) in $2018-2020$ along with the 519 supernovae recovered by ASAS-SN and 516 additional $m_{peak}\leq18$ mag supernovae missed by ASAS-SN. Our statistical analysis focuses primarily on the 984 supernovae discovered or recovered in ASAS-SN $g$-band observations. The complete sample of 2427 ASAS-SN… ▽ More We catalog the 443 bright supernovae discovered by the All-Sky Automated Survey for Supernovae (ASAS-SN) in $2018-2020$ along with the 519 supernovae recovered by ASAS-SN and 516 additional $m_{peak}\leq18$ mag supernovae missed by ASAS-SN. Our statistical analysis focuses primarily on the 984 supernovae discovered or recovered in ASAS-SN $g$-band observations. The complete sample of 2427 ASAS-SN supernovae includes earlier $V$-band samples and unrecovered supernovae. For each supernova, we identify the host galaxy, its UV to mid-IR photometry, and the offset of the supernova from the center of the host. Updated light curves, redshifts, classifications, and host galaxy identifications supersede earlier results. With the increase of the limiting magnitude to $g\leq18$ mag, the ASAS-SN sample is roughly complete up to $m_{peak}=16.7$ mag and is $90\%$ complete for $m_{peak}\leq17.0$ mag. This is an increase from the $V$-band sample where it was roughly complete up to $m_{peak}=16.2$ mag and $70\%$ complete for $m_{peak}\leq17.0$ mag. △ Less

Submitted 24 February, 2023; v1 submitted 12 October, 2022; originally announced October 2022.

Comments: 14 pages, 7 figures, 4 tables. Updated to reflect changes made in the published version. Tables containing the catalog data presented in this submission are included in machine-readable format as ancillary files

Journal ref: MNRAS, 520, 4356 (2023)

The Value-Added Catalog of ASAS-SN Eclipsing Binaries II: Properties of Extra-Physics Systems

Authors: D. M. Rowan, T. Jayasinghe, K. Z. Stanek, C. S. Kochanek, Todd A. Thompson, B. J. Shappee, T. W. -S. Holoien, J. L. Prieto, W. Giles

Abstract: Detached eclipsing binaries are the primary tool used to measure precise masses and radii of stars. In our previous paper estimating the parameters of more than 30,000 detached eclipsing binaries, we identified 766 eclipsing binaries with additional features in their All-Sky Automated Survey for Supernovae (ASAS-SN) and Transiting Exoplanet Survey Satellite (TESS) light curves. Here, we characteri… ▽ More Detached eclipsing binaries are the primary tool used to measure precise masses and radii of stars. In our previous paper estimating the parameters of more than 30,000 detached eclipsing binaries, we identified 766 eclipsing binaries with additional features in their All-Sky Automated Survey for Supernovae (ASAS-SN) and Transiting Exoplanet Survey Satellite (TESS) light curves. Here, we characterize these "extra-physics" systems, identifying eclipsing binaries with spotted stars, pulsating components, and candidate triple/quadruple systems. We use the Gaia, ATLAS, ZTF, and ASAS-SN variable star catalogs to consider possible blends. We use MIST isochrones and evolutionary tracks to identify systems with main sequence, subgiant, and giant primaries and highlight systems in sparsely populated regions of the color-magnitude diagram. We find that the orbital period distribution of spotted binaries is divided by evolutionary state and find 68 with X-ray detections. For the candidate triple/quadruples and pulsating systems, we calculate the extra orbital/pulsational period and identify systems with resonances. Finally, we highlight a number of exotic systems, including eclipsing CVs, subdwarfs, and binaries with disks. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: 19 pages, 14 figures. Submitted to MNRAS

On Odd Radio Circles as Supernova Remnants

Authors: Sumit K. Sarbadhicary, Todd A. Thompson, Laura A. Lopez, Smita Mathur

Abstract: The origin of arcmin-sized Odd Radio Circles (ORCs) found in modern all-sky radio surveys remain uncertain, with explanations ranging from starburst/AGN-driven shocks to supernova remnants (SNRs) in the low-density ambient medium. Using well-calibrated radio light curve models, we assess the possibility that ORCs are radio SNRs evolving in low ambient densities. Our models imply that ORCs 1-5 and… ▽ More The origin of arcmin-sized Odd Radio Circles (ORCs) found in modern all-sky radio surveys remain uncertain, with explanations ranging from starburst/AGN-driven shocks to supernova remnants (SNRs) in the low-density ambient medium. Using well-calibrated radio light curve models, we assess the possibility that ORCs are radio SNRs evolving in low ambient densities. Our models imply that ORCs 1-5 and J0624-6948 (near the LMC) as SNRs must be within 200 kpc and 100 kpc from the Sun respectively, given their observed flux densities and angular sizes. To be evolving in the circumgalactic medium of the Milky Way, our models require ORCs 1-5 to be ejecta-dominated SNRs within 50 kpc, evolving in ambient densities of $(0.2-1.2) \times 10^{-3}$ cm$^{-3}$. However, this is statistically unlikely because ORCs 1-5 would have ages $<640$ yrs, much smaller than their expected lifetimes of $\gtrsim$10$^5$ yrs at these densities, and because the low SN rate and steep profile of the stellar halo imply a negligible number of ORC-like SNRs within 50 kpc. The circumgalactic medium SNR scenario for J0624-6948 is more likely (though still low probability) compared to ORCs 1-5, as our models allow J0624-6948 to be $\lesssim$3000 yrs. On the other hand, the interpretation of J0624-6948 as a Sedov-Taylor SNR at 50 kpc (LMC) distance is possible for a wide range of ambient densities ($6 \times 10^{-4} - 0.5$ cm$^{-3}$) and ages $\sim$$(0.2-2.6) \times 10^4$ yr, while also being consistent with the local HI environment. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: 9 pages, 4 figures. Submitted to MNRAAS

X-ray Properties of NGC 253's Starburst-Driven Outflow

Authors: Sebastian Lopez, Laura A. Lopez, Dustin D. Nguyen, Todd A. Thompson, Smita Mathur, Alberto D. Bolatto, Neven Vulic, Amy Sardone

Abstract: We analyze image and spectral data from $\approx$365~ks of observations from the {\it Chandra} X-ray Observatory of the nearby, edge-on starburst galaxy NGC 253 to constrain properties of the hot phase of the outflow. We focus our analysis on the $-$1.1 to $+$0.63 kpc region of the outflow and define several regions for spectral extraction where we determine best-fit temperatures and metal abundan… ▽ More We analyze image and spectral data from $\approx$365~ks of observations from the {\it Chandra} X-ray Observatory of the nearby, edge-on starburst galaxy NGC 253 to constrain properties of the hot phase of the outflow. We focus our analysis on the $-$1.1 to $+$0.63 kpc region of the outflow and define several regions for spectral extraction where we determine best-fit temperatures and metal abundances. We find that the temperatures and electron densities peak in the central $\sim$250 pc region of the outflow and decrease with distance. These temperature and density profiles are in disagreement with an adiabatic spherically expanding starburst wind model and suggest the presence of additional physics such as mass loading and non-spherical outflow geometry. Our derived temperatures and densities yield few-Myr cooling times in the nuclear region, which may imply that the hot gas can undergo bulk radiative cooling as it escapes along the minor axis. Our metal abundances of O, Ne, Mg, Si, S, and Fe all peak in the central region and decrease with distance along the outflow, with the exception of Ne which maintains a flat distribution. The metal abundances indicate significant dilution outside of the starburst region. We also find estimates on the mass outflow rates which are $2.8\:M_{\odot}/\rm{yr}$ in the northern outflow and $3.2\:M_{\odot}/\rm{yr}$ in the southern outflow. Additionally, we detect emission from charge exchange and find it has a significant contribution ($20-42$\%) to the total broad-band ($0.5-7$ keV) X-ray emission in the central and southern regions of the outflow. △ Less

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

Comments: 15 pages, 9 figures; submitted (19 September 2022), accepted ApJ (22 November 2022)

ATLASGAL -- Star forming efficiencies and the Galactic star formation rate

Authors: M. R. A. Wells, J. S. Urquhart, T. J. T. Moore, K. E. Browning, S. E. Ragan, A. J. Rigby, D. J. Eden, M. A. Thompson

Abstract: The ATLASGAL survey has characterised the properties of approximately 1000 embedded HII regions and found an empirical relationship between the clump mass and bolometric luminosity that covers 3-4 orders of magnitude. Comparing this relation with simulated clusters drawn from an initial mass function and using different star formation efficiencies we find that a single value is unable to fit the o… ▽ More The ATLASGAL survey has characterised the properties of approximately 1000 embedded HII regions and found an empirical relationship between the clump mass and bolometric luminosity that covers 3-4 orders of magnitude. Comparing this relation with simulated clusters drawn from an initial mass function and using different star formation efficiencies we find that a single value is unable to fit the observed luminosity to mass ($L/M$) relation. We have used a Monte Carlo simulation to generate 200,000 clusters using the $L/M$-ratio as a constraint to investigate how the star formation efficiency changes as a function of clump mass. This has revealed that the star formation efficiency decreases with increasing clump mass with a value of 0.2 for clumps with masses of a few hundred solar masses and dropping to 0.08 for clumps with masses of a few thousand solar masses. We find good agreement between our results and star formation efficiencies determined from counts of embedded objects in nearby molecular clouds. Using the star formation efficiency relationship and the infrared excess time for embedded star formation of $2\pm1$, Myr we estimate the Galactic star formation rate to be approximately $0.9\pm0.45$ Msun yr$^{-1}$, which is in good agreement with previously reported values. This model has the advantage of providing a direct means of determining the star formation rate and avoids the difficulties encountered in converting infrared luminosities to stellar mass that affect previous galactic and extragalactic studies. △ Less

Submitted 7 September, 2022; v1 submitted 23 August, 2022; originally announced August 2022.

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

The Early Evolution of Magnetar Rotation I: Slowly Rotating "Normal" Magnetars

Authors: Tejas Prasanna, Matthew S. B. Coleman, Matthias J. Raives, Todd A. Thompson

Abstract: In the seconds following their formation in core-collapse supernovae, "proto"-magnetars drive neutrino-heated magneto-centrifugal winds. Using a suite of two-dimensional axisymmetric MHD simulations, we show that relatively slowly rotating magnetars with initial spin periods of $P_{\star0}=50-500$ ms spin down rapidly during the neutrino Kelvin-Helmholtz cooling epoch. These initial spin periods a… ▽ More In the seconds following their formation in core-collapse supernovae, "proto"-magnetars drive neutrino-heated magneto-centrifugal winds. Using a suite of two-dimensional axisymmetric MHD simulations, we show that relatively slowly rotating magnetars with initial spin periods of $P_{\star0}=50-500$ ms spin down rapidly during the neutrino Kelvin-Helmholtz cooling epoch. These initial spin periods are representative of those inferred for normal Galactic pulsars, and much slower than those invoked for gamma-ray bursts and super-luminous supernovae. Since the flow is non-relativistic at early times, and because the Alfvén radius is much larger than the proto-magnetar radius, spindown is millions of times more efficient than the typically-used dipole formula. Quasi-periodic plasmoid ejections from the closed zone enhance spindown. For polar magnetic field strengths $B_0\gtrsim5\times10^{14}$ G, the spindown timescale can be shorter than than the Kelvin-Helmholtz timescale. For $B_0\gtrsim10^{15}$ G, it is of order seconds in early phases. We compute the spin evolution for cooling proto-magnetars as a function of $B_0$, $P_{\star0}$, and mass ($M$). Proto-magnetars born with $B_0$ greater than $\simeq1.3\times10^{15}\,{\rm\,G}\,(P_{\star0}/{400\,\rm\,ms})^{-1.4}(M/1.4\,{\rm M}_\odot)^{2.2}$ spin down to periods $> 1$ s in just the first few seconds of evolution, well before the end of the cooling epoch and the onset of classic dipole spindown. Spindown is more efficient for lower $M$ and for larger $P_{\star0}$. We discuss the implications for observed magnetars, including the discrepancy between their characteristic ages and supernova remnant ages. Finally, we speculate on the origin of 1E 161348-5055 in the remnant RCW 103, and the potential for other ultra-slowly rotating magnetars. △ Less

Submitted 18 August, 2022; originally announced August 2022.

Comments: 16 pages, 10 figures

Journal ref: Monthly Notices of the Royal Astronomical Society 2022, Volume 517, Issue 2, pp.3008-3023

A search for compact object companions to high mass function single-lined spectroscopic binaries in Gaia DR3

Authors: T. Jayasinghe, D. M. Rowan, Todd A. Thompson, C. S. Kochanek, K. Z. Stanek

Abstract: Gaia DR3 provides >181, 000 radial velocity solutions for single-lined spectroscopic binaries (SB1s) which can be used to search for non-interacting compact object+star binary candidates by selecting systems with large mass functions. We selected 234 such systems and identified 114 systems with good RV solutions in DR3. We used light curves from ASAS-SN and TESS to identify and remove 31 eclipsing… ▽ More Gaia DR3 provides >181, 000 radial velocity solutions for single-lined spectroscopic binaries (SB1s) which can be used to search for non-interacting compact object+star binary candidates by selecting systems with large mass functions. We selected 234 such systems and identified 114 systems with good RV solutions in DR3. We used light curves from ASAS-SN and TESS to identify and remove 31 eclipsing binaries to produce a catalog of 80 compact object+star candidates, including 38 ellipsoidal variables. The positions of these candidates on Gaia and 2MASS CMDs suggest that many of these systems are binaries with luminous companions. We compared the periods and eccentricities of detached eclipsing binaries in Gaia DR3 and ASAS-SN, and found that ~11% and ~60% of the binaries had different periods and eccentricities. We also compared RV solutions for 311 binaries in both Gaia DR3 and the Ninth Catalog of Spectroscopic Binary Orbits (SB9), and found similar results. We do not identify any strong candidates for non-interacting compact object+star binaries. △ Less

Submitted 11 July, 2022; originally announced July 2022.

Comments: 15 pages, 10 figures. Submitted to MNRAS

The Galactic Nova Rate: Estimates from the ASAS-SN and Gaia Surveys

Authors: A. Kawash, L. Chomiuk, J. Strader, K. V. Sokolovsky, E. Aydi, C. S. Kochanek, K. Z. Stanek, Z. Kostrzewa-Rutkowska, S. T. Hodgkin, K. Mukai, B. Shappee, T. Jayasinghe, M. Rizzo Smith, T. W. -S. Holoien, J. L. Prieto, T. A. Thompson

Abstract: We present the first estimate of the Galactic nova rate based on optical transient surveys covering the entire sky. Using data from the All-Sky Automated Survey for Supernovae (ASAS-SN) and \textit{Gaia} -- the only two all-sky surveys to report classical nova candidates -- we find 39 confirmed Galactic novae and 7 additional unconfirmed candidates discovered from 2019--2021, yielding a nova disco… ▽ More We present the first estimate of the Galactic nova rate based on optical transient surveys covering the entire sky. Using data from the All-Sky Automated Survey for Supernovae (ASAS-SN) and \textit{Gaia} -- the only two all-sky surveys to report classical nova candidates -- we find 39 confirmed Galactic novae and 7 additional unconfirmed candidates discovered from 2019--2021, yielding a nova discovery rate of $\approx 14$ yr$^{-1}$. Using accurate Galactic stellar mass models, three-dimensional dust maps, and incorporating realistic nova light curves, we have built a sophisticated Galactic nova model that allows an estimate of the recovery fraction of Galactic novae from these surveys over this time period. The observing capabilities of each survey are distinct: the high cadence of ASAS-SN makes it sensitive to fast novae, while the broad observing filter and high spatial resolution of \textit{Gaia} make it more sensitive to highly reddened novae across the entire Galactic plane and bulge. Despite these differences, we find that ASAS-SN and \textit{Gaia} give consistent Galactic nova rates, with a final joint nova rate of $26 \pm 5$ yr$^{-1}$. This inferred nova rate is substantially lower than found by many other recent studies. Critically assessing the systematic uncertainties in the Galactic nova rate, we argue that the role of faint fast-fading novae has likely been overestimated, but that subtle details in the operation of transient alert pipelines can have large, sometimes unappreciated effects on transient recovery efficiency. Our predicted nova rate can be directly tested with forthcoming red/near-infrared transient surveys in the southern hemisphere. △ Less

Submitted 28 June, 2022; originally announced June 2022.

Comments: 24 pages, 9 figures

Towards Powerful Probes of Neutrino Self-Interactions in Supernovae

Authors: Po-Wen Chang, Ivan Esteban, John F. Beacom, Todd A. Thompson, Christopher M. Hirata

Abstract: Neutrinos remain mysterious. As an example, enhanced self-interactions ($ν$SI), which would have broad implications, are allowed. At the high neutrino densities within core-collapse supernovae, $ν$SI should be important, but robust observables have been lacking. We show that $ν$SI make neutrinos form a tightly coupled fluid that expands under relativistic hydrodynamics. The outflow becomes either… ▽ More Neutrinos remain mysterious. As an example, enhanced self-interactions ($ν$SI), which would have broad implications, are allowed. At the high neutrino densities within core-collapse supernovae, $ν$SI should be important, but robust observables have been lacking. We show that $ν$SI make neutrinos form a tightly coupled fluid that expands under relativistic hydrodynamics. The outflow becomes either a burst or a steady-state wind; which occurs here is uncertain. Though the diffusive environment where neutrinos are produced may make a wind more likely, further work is needed to determine when each case is realized. In the burst-outflow case, $ν$SI increase the duration of the neutrino signal, and even a simple analysis of SN 1987A data has powerful sensitivity. For the wind-outflow case, we outline several promising ideas that may lead to new observables. Combined, these results are important steps towards solving the 35-year-old puzzle of how $ν$SI affect supernovae. △ Less

Submitted 28 June, 2023; v1 submitted 24 June, 2022; originally announced June 2022.

Comments: 5+19 pages, 4+14 figures. Accepted for publication in Phys.Rev.Lett. Match the version to be published: minor revisions for clarity; moved discussion about supernova neutrino properties to Supplemental Material C; updated references

Journal ref: Phys. Rev. Lett. 131, 071002 (2023)

Chandra, HST/STIS, NICER, Swift, and TESS Detail the Flare Evolution of the Repeating Nuclear Transient ASASSN-14ko

Authors: Anna V. Payne, Katie Auchettl, Benjamin J. Shappee, Christopher S. Kochanek, Patricia T. Boyd, Thomas W. -S. Holoien, Michael M. Fausnaugh, Chris Ashall, Jason T. Hinkle, Patrick J. Vallely, K. Z. Stanek, Todd A. Thompson

Abstract: ASASSN-14ko is a nuclear transient at the center of the AGN ESO 253-G003 that undergoes periodic flares. Optical flares were first observed in 2014 by the All-Sky Automated Survey for Supernovae (ASAS-SN) and their peak times are well-modeled with a period of $115.2^{+1.3}_{-1.2}$ days and period derivative of $-0.0026 \pm 0.0006$. Here we present ASAS-SN, Chandra, HST/STIS, NICER, Swift, and TESS… ▽ More ASASSN-14ko is a nuclear transient at the center of the AGN ESO 253-G003 that undergoes periodic flares. Optical flares were first observed in 2014 by the All-Sky Automated Survey for Supernovae (ASAS-SN) and their peak times are well-modeled with a period of $115.2^{+1.3}_{-1.2}$ days and period derivative of $-0.0026 \pm 0.0006$. Here we present ASAS-SN, Chandra, HST/STIS, NICER, Swift, and TESS data for the flares that occurred in December 2020, April 2021, July 2021, and November 2021. The HST/STIS UV spectra evolve from blue shifted broad absorption features to red shifted broad emission features over $\sim$10 days. The Swift UV/optical light curves peaked as predicted by the timing model, but the peak UV luminosities varied between flares and the UV flux in July 2021 was roughly half the brightness of all other peaks. The X-ray luminosities consistently decreased and the spectra became harder during the UV/optical rise but apparently without changes in absorption. Finally, two high-cadence TESS light curves from December 2020 and November 2018 showed that the slopes during the rising and declining phases changed over time, which indicates some stochasticity in the flare's driving mechanism. ASASSN-14ko remains observationally consistent with a repeating partial tidal disruption event, but, these rich multi-wavelength data are in need of a detailed theoretical model. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: 25 pages, 14 figures, 4 tables; Submitted to ApJ, comments welcome

BSM Targets at a "Target-less DUNE"

Authors: Vedran Brdar, Bhaskar Dutta, Wooyoung Jang, Doojin Kim, Ian M. Shoemaker, Zahra Tabrizi, Adrian Thompson, Jaehoon Yu

Abstract: In this work we demonstrate that a future accelerator-based neutrino experiment such as DUNE can greatly increase its sensitivity to a variety of new physics scenarios by operating in a mode where the proton beam impinges on a beam dump. We consider two new physics scenarios, namely light dark matter (LDM) and axion-like particles (ALPs) and show that by utilizing a dump mode at a DUNE-like experi… ▽ More In this work we demonstrate that a future accelerator-based neutrino experiment such as DUNE can greatly increase its sensitivity to a variety of new physics scenarios by operating in a mode where the proton beam impinges on a beam dump. We consider two new physics scenarios, namely light dark matter (LDM) and axion-like particles (ALPs) and show that by utilizing a dump mode at a DUNE-like experiment, unexplored new regions of parameter space can be probed with an exposure of only 3 months with half of its expected initial beam power. Specifically, target-less configuration of future high intensity neutrino experiments will probe the parameter space for thermal relic DM as well as the QCD axion (DFSZ and KSVZ). The strength of such configuration in the context of new physics searches stems from the fact that the neutrino flux is significantly reduced compared to that of the target, resulting in much smaller backgrounds from neutrino interactions. We have verified this in detail by explicitly computing neutrino fluxes which we make publicly available in order to facilitate further studies with a target-less configuration. △ Less

Submitted 4 April, 2023; v1 submitted 13 June, 2022; originally announced June 2022.

Comments: Based on the published version on PRD

Journal ref: Phys.Rev.D 107 (2023) 5, 055043