-
An ALCHEMI inspection of sulphur-bearing species towards the central molecular zone of NGC 253
Authors:
M. Bouvier,
S. Viti,
E. Behrens,
J. Butterworth,
K. -Y. Huang,
J. G. Mangum,
N. Harada,
S. Martín,
V. M. Rivilla,
S. Muller,
K. Sakamoto,
Y. Yoshimura,
K. Tanaka,
K. Nakanishi,
R. Herrero-Illana,
L. Colzi,
M. D. Gorski,
C. Henkel,
P. K. Humire,
D. S. Meier,
P. P. van der Werf,
Y. T. Yan
Abstract:
Sulphur-bearing species are detected in various environments within Galactic star-forming regions and are particularly abundant in the gas phase of outflows and shocks, and photo-dissociation regions. In this work, we aim to investigate the nature of the emission from the most common sulphur-bearing species observable at millimetre wavelengths towards the nuclear starburst of the galaxy NGC 253. W…
▽ More
Sulphur-bearing species are detected in various environments within Galactic star-forming regions and are particularly abundant in the gas phase of outflows and shocks, and photo-dissociation regions. In this work, we aim to investigate the nature of the emission from the most common sulphur-bearing species observable at millimetre wavelengths towards the nuclear starburst of the galaxy NGC 253. We intend to understand which type of regions are probed by sulphur-bearing species and which process(es) dominate(s) the release of sulphur into the gas phase. We used the high-angular resolution (1.6" or 27 pc) observations from the ALCHEMI ALMA Large Program to image several sulphur-bearing species towards the central molecular zone (CMZ) of NGC 253. We performed local thermodynamic equilibrium (LTE) and non-LTE large velocity gradient (LVG) analyses to derive the physical conditions of the gas in which S-bearing species are emitted, and their abundance ratios across the CMZ. Finally, we compared our results with previous ALCHEMI studies and a few selected Galactic environments. We found that not all sulphur-bearing species trace the same type of gas: strong evidence indicates that H2S and part of the emission of OCS, H2CS, and SO, are tracing shocks whilst part of SO and CS emission rather trace the dense molecular gas. For some species, such as CCS and SO2, we could not firmly conclude on their origin of emission. The present analysis indicates that the emission from most sulphur-bearing species throughout the CMZ is likely dominated by shocks associated with ongoing star formation. In the inner part of the CMZ where the presence of super star clusters was previously indicated, we could not distinguish between shocks or thermal evaporation as the main process releasing the S-bearing species.
△ Less
Submitted 14 May, 2024;
originally announced May 2024.
-
Dark Matter from Anomaly Cancellation at the LHC
Authors:
Jon Butterworth,
Hridoy Debnath,
Pavel Fileviez Perez,
Yoran Yeh
Abstract:
We discuss a class of theories that predict a fermionic dark matter candidate from gauge anomaly cancellation. As an explicit example, we study the predictions in theories where the global symmetry associated with baryon number is promoted to a local gauge symmetry. In this context the symmetry-breaking scale has to be below the multi-TeV scale in order to be in agreement with the cosmological con…
▽ More
We discuss a class of theories that predict a fermionic dark matter candidate from gauge anomaly cancellation. As an explicit example, we study the predictions in theories where the global symmetry associated with baryon number is promoted to a local gauge symmetry. In this context the symmetry-breaking scale has to be below the multi-TeV scale in order to be in agreement with the cosmological constraints on the dark matter relic density. The new physical "Cucuyo" Higgs boson in the theory has very interesting properties, decaying mainly into two photons in the low mass region, and mainly into dark matter in the intermediate mass region. We study the most important signatures at the Large Hadron Collider, evaluating the experimental bounds. We discuss the correlation between the dark matter relic density, direct detection and collider constraints. We find that these theories are still viable, and are susceptible to being probed in current, and future high-luminosity, running.
△ Less
Submitted 6 May, 2024;
originally announced May 2024.
-
A Comparative Study of the Ground State Transitions of CO and [C I] as Molecular Gas Tracers at High Redshift
Authors:
Marta Frias Castillo,
Matus Rybak,
Jacqueline A. Hodge,
Paul Van der Werk,
Ian Smail,
Joshua Butterworth,
Jasper Jansen,
Theodoros Topkaras,
Chian-Chou Chen,
Scott C. Chapman,
Axel Weiss,
Hiddo Algera,
Jack E. Birkin,
Elisabete da Cunha,
Jianhang Chen,
Helmut Dannerbauer,
E. F. Jiménez-Andrade,
Soh Ikarashi,
Cheng-Lin Liao,
Eric J. Murphy,
A. M. Swinbank,
Fabian Walter,
Gabriela Calistro Rivera,
R. J. Ivison,
Claudia del P. Lagos
Abstract:
The CO(1--0) and [\ion{C}{1}](1--0) emission lines are well-established tracers of cold molecular gas mass in local galaxies. At high redshift, where the interstellar medium (ISM) is likely to be denser, there have been limited direct comparisons of both ground state transitions. Here we present a study of CO(1--0) and [\ion{C}{1}](1--0) emission in a sample of 20 unlensed dusty, star-forming gala…
▽ More
The CO(1--0) and [\ion{C}{1}](1--0) emission lines are well-established tracers of cold molecular gas mass in local galaxies. At high redshift, where the interstellar medium (ISM) is likely to be denser, there have been limited direct comparisons of both ground state transitions. Here we present a study of CO(1--0) and [\ion{C}{1}](1--0) emission in a sample of 20 unlensed dusty, star-forming galaxies at $z=2-5$. The CO(1--0)/[\ion{C}{1}](1--0) ratio is constant up to at least $z=5$, supporting the use of [CI](1-0) as a gas mass tracer. PDR modelling of the available data indicates a median H$_2$ density of log$(n~[$cm$^{-3}])=4.7\pm0.2$, and UV radiation field log$(G_{\mathrm{UV}} [G$_0$])=3.2\pm0.2$. We use the CO(1--0), [\ion{C}{1}](1--0) and 3mm dust continuum measurements to cross--calibrate the respective gas mass conversion factors, finding no dependence of these factors on either redshift or infrared luminosity. Assuming a variable CO conversion factor then implies [\ion{C}{1}] and dust conversion factors that differ from canonically assumed values but are consistent with the solar/super-solar metallicities expected for our sources. Radiative transfer modelling shows that the warmer CMB at high redshift can significantly affect the [\ion{C}{1}] as well as CO emission, which can change the derived molecular gas masses by up to 70\% for the coldest kinetic gas temperatures expected. Nevertheless, we show that the magnitude of the effect on the ratio of the tracers is within the known scatter of the $L'_\mathrm{CO}-L'_\mathrm{[CI]}$ relation. Further determining the absolute decrease of individual line intensities will require well-sampled spectral line energy distributions (SLEDs) to model the gas excitation conditions in more detail.
△ Less
Submitted 8 April, 2024;
originally announced April 2024.
-
Molecular isotopologue measurements toward super star clusters and the relation to their ages in NGC253 with ALCHEMI
Authors:
J. Butterworth,
S. Viti,
P. P. Van der Werf,
J. G. Mangum,
S. Martín,
N. Harada,
K. L. Emig,
S. Muller,
K. Sakamoto,
Y. Yoshimura,
K. Tanaka,
R. Herrero-Illana,
L. Colzi,
V. M. Rivilla,
K. Y. Huang,
M. Bouvier,
E. Behrens,
C. Henkel,
Y. T. Yan,
D. S. Meier,
D. Zhou
Abstract:
Determining the evolution of the CNO isotopes in the interstellar medium (ISM) of starburst galaxies can yield important constraints on the ages of superstar clusters (SSCs), or on other aspects and contributing factors of their evolution. Due to the time-dependent nature of the abundances of isotopes within the ISM as they are supplied from processes such as nucleosynthesis or chemical fractionat…
▽ More
Determining the evolution of the CNO isotopes in the interstellar medium (ISM) of starburst galaxies can yield important constraints on the ages of superstar clusters (SSCs), or on other aspects and contributing factors of their evolution. Due to the time-dependent nature of the abundances of isotopes within the ISM as they are supplied from processes such as nucleosynthesis or chemical fractionation, this provides the possible opportunity to probe the ability of isotopes ratios to trace the ages of high star forming regions, such as SSCs. The goal of this study is to investigate whether the isotopic variations in SSC regions within NGC253 are correlated with their different ages as derived from stellar population modelling. We have measured abundance ratios of CO, HCN and HCO$^+$ isotopologues in six regions containing SSCs within NGC253 using high spatial resolution (1.6",$\sim 28$pc) data from the ALCHEMI (ALma Comprehensive High-resolution Extragalactic Molecular Inventory) ALMA Large program. We have then analysed these ratios using RADEX radiative transfer modelling, with the parameter space sampled using the nested sampling Monte Carlo algorithm MLFriends. These abundance ratios were then compared to ages predicted in each region via the fitting of observed star formation tracers (such as Br$γ$) to starburst stellar population evolution models. We do not find any significant trend with age for the CO and HCN isotopologue ratios on the timescales for the ages of the SSC* regions observed. The driving factors of these ratios within SSCs could be the Initial Mass Function as well as possibly fractionation effects. To further probe these effects in SSCs over time a larger sample of SSCs must be observed spanning a larger age range.
△ Less
Submitted 16 February, 2024;
originally announced February 2024.
-
Volume density structure of the NGC253 CMZ through ALCHEMI excitation analysis
Authors:
Kunihiko Tanaka,
Jeffrey G. Mangum,
Serena Viti,
Sergio Martin,
Nanase Harada,
Kazushi Sakamoto,
Sebastien Muller,
Yuki Yoshimura,
Kouichiro Nakanishi,
Ruben Herrero Illana,
Kimberly L. Emig,
S. Muhle,
Hiroyuki Kaneko,
Tomoka Tosaki,
Erica Behrens,
Victor M. Rivilla,
Laura Colzi,
Yuri Nishimura,
P. K. Humire,
Mathilde Bouvier,
Ko-Yun Huang,
Joshua Butterworth,
David S. Meier,
Paul P. van der Werf
Abstract:
We present a spatially-resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the ALMA Large program ALCHEMI, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way's Galactic Center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3-7 transitions from nine molecular spec…
▽ More
We present a spatially-resolved excitation analysis for the central molecular zone (CMZ) of the starburst galaxy NGC 253 using the data from the ALMA Large program ALCHEMI, whereby we explore parameters distinguishing NGC 253 from the quiescent Milky Way's Galactic Center (GC). Non-LTE analyses employing a hierarchical Bayesian framework are applied to Band 3-7 transitions from nine molecular species to delineate the position-position-velocity distributions of column density ($N_\mathrm{H_2}$), volume density ($n_\mathrm{H_2}$), and temperature ($T_\mathrm{kin}$) at 27 pc resolution. Two distinct components are detected: a low-density component with $(n_\mathrm{H_2},\ T_\mathrm{kin})\sim(10^{3.3}\ \mathrm{cm}^{-3}, 85 K)$ and a high-density component with $(n_\mathrm{H_2},\ T_\mathrm{kin})\sim (10^{4.4}\ \mathrm{cm}^{-3}, 110\ \mathrm{K})$, separated at $n_\mathrm{H_2}\sim10^{3.8}\ \mathrm{cm}^{-3}$. NGC 253 has $\sim10$ times the high-density gas mass and $\sim3$ times the dense-gas mass fraction of the GC. These properties are consistent with their HCN/CO ratio but cannot alone explain the factor of $\sim30$ difference in their star formation efficiencies (SFEs), contradicting the dense-gas mass to star formation rate scaling law. The $n_\mathrm{H_2}$ histogram toward NGC 253 exhibits a shallow declining slope up to $n_\mathrm{H_2}\sim10^6\ \mathrm{cm}^{-3}$, while that of the GC steeply drops in $n_\mathrm{H_2}\gtrsim10^{4.5}\ \mathrm{cm}^{-3}$ and vanishes at $10^5\ \mathrm{cm}^{-3}$. Their dense-gas mass fraction ratio becomes consistent with their SFEs when the threshold $n_\mathrm{H_2}$ for the dense gas is taken at $\sim 10^{4.2\mbox{-}4.6}\ \mathrm{cm}^{-3}$. The rich abundance of gas above this density range in the NGC 253 CMZ, or its scarcity in the GC, is likely to be the critical difference characterizing the contrasting star formation in the centers of the two galaxies.
△ Less
Submitted 20 November, 2023;
originally announced November 2023.
-
Understanding Molecular Abundances in Star-Forming Regions Using Interpretable Machine Learning
Authors:
Johannes Heyl,
Joshua Butterworth,
Serena Viti
Abstract:
Astrochemical modelling of the interstellar medium typically makes use of complex computational codes with parameters whose values can be varied. It is not always clear what the exact nature of the relationship is between these input parameters and the output molecular abundances. In this work, a feature importance analysis is conducted using SHapley Additive exPlanations (SHAP), an interpretable…
▽ More
Astrochemical modelling of the interstellar medium typically makes use of complex computational codes with parameters whose values can be varied. It is not always clear what the exact nature of the relationship is between these input parameters and the output molecular abundances. In this work, a feature importance analysis is conducted using SHapley Additive exPlanations (SHAP), an interpretable machine learning technique, to identify the most important physical parameters as well as their relationship with each output. The outputs are the abundances of species and ratios of abundances. In order to reduce the time taken for this process, a neural network emulator is trained to model each species' output abundance and this emulator is used to perform the interpretable machine learning. SHAP is then used to further explore the relationship between the physical features and the abundances for the various species and ratios we considered. \ce{H2O} and CO's gas phase abundances are found to strongly depend on the metallicity. \ce{NH3} has a strong temperature dependence, with there being two temperature regimes (< 100 K and > 100K). By analysing the chemical network, we relate this to the chemical reactions in our network and find the increased temperature results in increased efficiency of destruction pathways. We investigate the HCN/HNC ratio and show that it can be used as a cosmic thermometer, agreeing with the literature. This ratio is also found to be correlated with the metallicity. The HCN/CS ratio serves as a density tracer, but also has three separate temperature-dependence regimes, which are linked to the chemistry of the two molecules.
△ Less
Submitted 13 September, 2023;
originally announced September 2023.
-
Understanding if molecular ratios can be used as diagnostics of AGN and starburst activity: The case of NGC 1068
Authors:
J. Butterworth,
J. Holdship,
S. Viti,
S. García-Burillo
Abstract:
Molecular line ratios, such as HCN(1-0)/HCO$^+$(1-0) and HCN(4-3)/CS(7-6), are routinely used to identify active galactic nuclei (AGN) activity in galaxies. Such ratios are, however, hard to interpret as they are highly dependent on the physics and energetics of the gas, and hence can seldom be used as a unique, unambiguous diagnostic. We used the composite galaxy NGC 1068 as a `laboratory' to inv…
▽ More
Molecular line ratios, such as HCN(1-0)/HCO$^+$(1-0) and HCN(4-3)/CS(7-6), are routinely used to identify active galactic nuclei (AGN) activity in galaxies. Such ratios are, however, hard to interpret as they are highly dependent on the physics and energetics of the gas, and hence can seldom be used as a unique, unambiguous diagnostic. We used the composite galaxy NGC 1068 as a `laboratory' to investigate whether molecular line ratios between HCN, HCO$^+$, and CS are useful tracers of AGN-dominated gas and determine the origin of the differences in such ratios across different types of gas. Such a determination will enable a more rigorous use of such ratios. First, we empirically examined the aforementioned ratios at different angular resolutions to quantify correlations. We then used local thermodynamic equilibrium (LTE) and non-LTE analyses coupled with Markov chain Monte Carlo (MCMC) sampling in order to determine the origin of the underlying differences in ratios. We propose that at high spatial resolution (< 50 pc) the HCN(4-3)/CS(2-1) is a reliable tracer of AGN activity. We also find that the variations in ratios are not a consequence of different densities or temperature but of different fractional abundances, yielding to the important result that it is essential to consider the chemical processes at play when drawing conclusions from radiative transfer calculations. From analyses at varying spatial scales, we find that previously proposed molecular line ratios, as well as a new one, have varying levels of consistency. We also determine from an investigation of radiative transfer modelling of our data that it is essential to consider the chemistry of the species when reaching conclusions from radiative transfer calculations.
△ Less
Submitted 21 September, 2022; v1 submitted 13 September, 2022;
originally announced September 2022.
-
The evolution of the UV luminosity and stellar mass functions of Lyman-alpha emitters from z~2 to z~6
Authors:
Sérgio Santos,
David Sobral,
Josh Butterworth,
Ana Paulino-Afonso,
Bruno Ribeiro,
Elisabete da Cunha,
João Calhau,
Ali Ahmad Khostovan,
Jorryt Matthee,
Pablo Arrabal Haro
Abstract:
We measure the evolution of the rest-frame UV luminosity function (LF) and the stellar mass function (SMF) of Lyman-alpha (Lya) emitters (LAEs) from z~2 to z~6 by exploring ~4000 LAEs from the SC4K sample. We find a correlation between Lya luminosity (LLya) and rest-frame UV (M_UV), with best-fit M_UV=-1.6+-0.2 log10(LLya/erg/s)+47+-12 and a shallower relation between LLya and stellar mass (Mstar)…
▽ More
We measure the evolution of the rest-frame UV luminosity function (LF) and the stellar mass function (SMF) of Lyman-alpha (Lya) emitters (LAEs) from z~2 to z~6 by exploring ~4000 LAEs from the SC4K sample. We find a correlation between Lya luminosity (LLya) and rest-frame UV (M_UV), with best-fit M_UV=-1.6+-0.2 log10(LLya/erg/s)+47+-12 and a shallower relation between LLya and stellar mass (Mstar), with best-fit log10( Mstar/Msun)=0.9+-0.1 log10(LLya/erg/s)-28+-4.0. An increasing LLya cut predominantly lowers the number density of faint M_UV and low Mstar LAEs. We estimate a proxy for the full UV LFs and SMFs of LAEs with simple assumptions of the faint end slope. For the UV LF, we find a brightening of the characteristic UV luminosity (M_UV*) with increasing redshift and a decrease of the characteristic number density (Phi*). For the SMF, we measure a characteristic stellar mass (Mstar*/Msun) increase with increasing redshift, and a Phi* decline. However, if we apply a uniform luminosity cut of log10 (LLya/erg/s) >= 43.0, we find much milder to no evolution in the UV and SMF of LAEs. The UV luminosity density (rho_UV) of the full sample of LAEs shows moderate evolution and the stellar mass density (rho_M) decreases, with both being always lower than the total rho_UV and rho_M of more typical galaxies but slowly approaching them with increasing redshift. Overall, our results indicate that both rho_UV and rho_M of LAEs slowly approach the measurements of continuum-selected galaxies at z>6, which suggests a key role of LAEs in the epoch of reionisation.
△ Less
Submitted 30 April, 2021;
originally announced May 2021.
-
Highlights of EPS HEP 2019
Authors:
Jon Butterworth
Abstract:
An opinionated and informal recap of highlights from the EPS HEP 2019 conference in Ghent, including some aspects of flavour physics, neutrinos, high-density QCD, astrophysics and energy frontier collider physics, and some thoughts about the future.
An opinionated and informal recap of highlights from the EPS HEP 2019 conference in Ghent, including some aspects of flavour physics, neutrinos, high-density QCD, astrophysics and energy frontier collider physics, and some thoughts about the future.
△ Less
Submitted 28 October, 2019;
originally announced October 2019.
-
The evolution of rest-frame UV properties, Lya EWs and the SFR-Stellar mass relation at z~2-6 for SC4K LAEs
Authors:
Sérgio Santos,
David Sobral,
Jorryt Matthee,
João Calhau,
Elisabete da Cunha,
Bruno Ribeiro,
Ana Paulino-Afonso,
Pablo Arrabal Haro,
Josh Butterworth
Abstract:
We explore deep rest-frame UV to FIR data in the COSMOS field to measure the individual spectral energy distributions (SED) of the ~4000 SC4K (Sobral et al. 2018) Lyman-alpha (Lya) emitters (LAEs) at z~2-6. We find typical stellar masses of 10$^{9.3\pm0.6}$ M$_{\odot}$ and star formation rates (SFR) of SFR$_{SED}=4.4^{+10.5}_{-2.4}$ M$_{\odot}$/yr and SFR$_{Lya}=5.9^{+6.3}_{-2.6}$ M$_{\odot}$/yr,…
▽ More
We explore deep rest-frame UV to FIR data in the COSMOS field to measure the individual spectral energy distributions (SED) of the ~4000 SC4K (Sobral et al. 2018) Lyman-alpha (Lya) emitters (LAEs) at z~2-6. We find typical stellar masses of 10$^{9.3\pm0.6}$ M$_{\odot}$ and star formation rates (SFR) of SFR$_{SED}=4.4^{+10.5}_{-2.4}$ M$_{\odot}$/yr and SFR$_{Lya}=5.9^{+6.3}_{-2.6}$ M$_{\odot}$/yr, combined with very blue UV slopes of beta=-2.1$^{+0.5}_{-0.4}$, but with significant variations within the population. M$_{UV}$ and beta are correlated in a similar way to UV-selected sources, but LAEs are consistently bluer. This suggests that LAEs are the youngest and/or most dust-poor subset of the UV-selected population. We also study the Lya rest-frame equivalent width (EW$_0$) and find 45 "extreme" LAEs with EW$_0>240$ A (3 $σ$), implying a low number density of $(7\pm1)\times10^{-7}$ Mpc$^{-3}$. Overall, we measure little to no evolution of the Lya EW$_0$ and scale length parameter ($w_0$) which are consistently high (EW$_0=140^{+280}_{-70}$ A, $w_0=129^{+11}_{-11}$ A) from z~6 to z~2 and below. However, $w_0$ is anti-correlated with M$_{UV}$ and stellar mass. Our results imply that sources selected as LAEs have a high Lya escape fraction (f$_{esc, Lya}$) irrespective of cosmic time, but f$_{esc, Lya}$ is still higher for UV-fainter and lower mass LAEs. The least massive LAEs ($<10^{9.5}$ M$_{\odot}$) are typically located above the star formation "Main Sequence" (MS), but the offset from the MS decreases towards z~6 and towards $10^{10}$ M$_{\odot}$. Our results imply a lack of evolution in the properties of LAEs across time and reveals the increasing overlap in properties of LAEs and UV-continuum selected galaxies as typical star-forming galaxies at high redshift effectively become LAEs.
△ Less
Submitted 10 January, 2020; v1 submitted 7 October, 2019;
originally announced October 2019.
-
Calibration and performance of the STAR Muon Telescope Detector using cosmic rays
Authors:
C. Yang,
X. J. Huang,
C. M. Du,
B. C. Huang,
Z. Ahammed,
A. Banerjee,
P. Bhattarari,
S. Biswas,
B. Bowen,
J. Butterworth,
M. Calderón de la Barca Sánchez,
H. Carson,
S. Chattopadhyay,
D. Cebra,
H. F. Chen,
J. P. Cheng,
M. Codrington,
G. Eppley,
C. Flores,
F. Geurts,
G. W. Hoffmann,
A. Jentsch,
A. Kesich,
C. Li,
Y. J. Li
, et al. (13 additional authors not shown)
Abstract:
We report the timing and spatial resolution from the Muon Telescope Detector (MTD) installed in the STAR experiment at RHIC. Cosmic ray muons traversing the STAR detector have an average transverse momentum of 6 GeV/c. Due to their very small multiple scattering, these cosmic muons provide an ideal tool to calibrate the detectors and measure their timing and spatial resolution. The values obtained…
▽ More
We report the timing and spatial resolution from the Muon Telescope Detector (MTD) installed in the STAR experiment at RHIC. Cosmic ray muons traversing the STAR detector have an average transverse momentum of 6 GeV/c. Due to their very small multiple scattering, these cosmic muons provide an ideal tool to calibrate the detectors and measure their timing and spatial resolution. The values obtained were ~100 ps and ~1-2 cm, respectively. These values are comparable to those obtained from cosmic-ray bench tests and test beams.
△ Less
Submitted 5 February, 2014;
originally announced February 2014.