Showing 1–50 of 289 results for author: Barnes, J

Asteroid (101955) Bennu in the Laboratory: Properties of the Sample Collected by OSIRIS-REx

Authors: Dante S. Lauretta, Harold C. Connolly, Jr., Joseph E. Aebersold, Conel M. O. D. Alexander, Ronald-L. Ballouz, Jessica J. Barnes, Helena C. Bates, Carina A. Bennett, Laurinne Blanche, Erika H. Blumenfeld, Simon J. Clemett, George D. Cody, Daniella N. DellaGiustina, Jason P. Dworkin, Scott A. Eckley, Dionysis I. Foustoukos, Ian A. Franchi, Daniel P. Glavin, Richard C. Greenwood, Pierre Haenecour, Victoria E. Hamilton, Dolores H. Hill, Takahiro Hiroi, Kana Ishimaru, Fred Jourdan , et al. (28 additional authors not shown)

Abstract: On 24 September 2023, the NASA OSIRIS-REx mission dropped a capsule to Earth containing approximately 120 g of pristine carbonaceous regolith from Bennu. We describe the delivery and initial allocation of this asteroid sample and introduce its bulk physical, chemical, and mineralogical properties from early analyses. The regolith is very dark overall, with higher-reflectance inclusions and particl… ▽ More On 24 September 2023, the NASA OSIRIS-REx mission dropped a capsule to Earth containing approximately 120 g of pristine carbonaceous regolith from Bennu. We describe the delivery and initial allocation of this asteroid sample and introduce its bulk physical, chemical, and mineralogical properties from early analyses. The regolith is very dark overall, with higher-reflectance inclusions and particles interspersed. Particle sizes range from sub-micron dust to a stone about 3.5 cm long. Millimeter-scale and larger stones typically have hummocky or angular morphologies. A subset of the stones appears mottled by brighter material that occurs as veins and crusts. Hummocky stones have the lowest densities and mottled stones have the highest. Remote sensing of the surface of Bennu detected hydrated phyllosilicates, magnetite, organic compounds, carbonates, and scarce anhydrous silicates, all of which the sample confirms. We also find sulfides, presolar grains, and, less expectedly, Na-rich phosphates, as well as other trace phases. The sample composition and mineralogy indicate substantial aqueous alteration and resemble those of Ryugu and the most chemically primitive, low-petrologic-type carbonaceous chondrites. Nevertheless, we find distinct hydrogen, nitrogen, and oxygen isotopic compositions, and some of the material we analyzed is enriched in fluid-mobile elements. Our findings underscore the value of sample return, especially for low-density material that may not readily survive atmospheric entry, and lay the groundwork for more comprehensive analyses. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 73 pages, 22 figures

The variable magnetic field of V889 Her and the challenge of detecting exoplanets around young Suns using Gaussian process regression

Authors: E. L. Brown, S. C. Marsden, S. V. Jeffers, A. Heitzmann, J. R. Barnes, C. P. Folsom

Abstract: Discovering exoplanets orbiting young Suns can provide insight into the formation and early evolution of our own solar system, but the extreme magnetic activity of young stars obfuscates exoplanet detection. Here we monitor the long-term magnetic field and chromospheric activity variability of the young solar analogue V889 Her, model the activity-induced radial velocity variations and evaluate the… ▽ More Discovering exoplanets orbiting young Suns can provide insight into the formation and early evolution of our own solar system, but the extreme magnetic activity of young stars obfuscates exoplanet detection. Here we monitor the long-term magnetic field and chromospheric activity variability of the young solar analogue V889 Her, model the activity-induced radial velocity variations and evaluate the impacts of extreme magnetism on exoplanet detection thresholds. We map the magnetic field and surface brightness for 14 epochs between 2004 and 2019. Our results show potential 3-4 yr variations of the magnetic field which evolves from weak and simple during chromospheric activity minima to strong and complex during activity maxima but without any polarity reversals. A persistent, temporally-varying polar spot coexists with weaker, short-lived lower-latitude spots. Due to their different decay time-scales, significant differential rotation and the limited temporal coverage of our legacy data, we were unable to reliably model the activity-induced radial velocity using Gaussian Process regression. Doppler Imaging can be a useful method for modelling the magnetic activity jitter of extremely active stars using data with large phase gaps. Given our data and using Doppler Imaging to filter activity jitter, we estimate that we could detect Jupiter-mass planets with orbital periods of $\sim$3 d. A longer baseline of continuous observations is the best observing strategy for the detection of exoplanets orbiting highly active stars. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Comments: Accepted by MNRAS

HRMOS White Paper: Science Motivation

Authors: Laura Magrini, Thomas Bensby, Anna Brucalassi, Sofia Randich, Robin Jeffries, Gayandhi de Silva, Asa Skuladottir, Rodolfo Smiljanic, Oscar Gonzalez, Vanessa Hill, Nadege Lagarde, Eline Tolstoy, Jose' Maria Arroyo-Polonio, Martina Baratella, John R. Barnes, Giuseppina Battaglia, Holger Baumgardt, Michele Bellazzini, Katia Biazzo, Angela Bragaglia, Bradley Carter, Giada Casali, Gabriele Cescutti, Camilla Danielski, Elisa Delgado Mena , et al. (30 additional authors not shown)

Abstract: The High-Resolution Multi-Object Spectrograph (HRMOS) is a facility instrument that we plan to propose for the Very Large Telescope (VLT) of the European Southern Observatory (ESO), following the initial presentation at the VLT 2030 workshop held at ESO in June 2019. HRMOS provides a combination of capabilities that are essential to carry out breakthrough science across a broad range of active res… ▽ More The High-Resolution Multi-Object Spectrograph (HRMOS) is a facility instrument that we plan to propose for the Very Large Telescope (VLT) of the European Southern Observatory (ESO), following the initial presentation at the VLT 2030 workshop held at ESO in June 2019. HRMOS provides a combination of capabilities that are essential to carry out breakthrough science across a broad range of active research areas from stellar astrophysics and exoplanet studies to Galactic and Local Group archaeology. HRMOS fills a gap in capabilities amongst the landscape of future instrumentation planned for the next decade. The key characteristics of HRMOS will be high spectral resolution (R = 60000 - 80000) combined with multi-object (20-100) capabilities and long term stability that will provide excellent radial velocity precision and accuracy (10m/s). Initial designs predict that a SNR~100 will be achievable in about one hour for a star with mag(AB) = 15, while with the same exposure time a SNR~ 30 will be reached for a star with mag(AB) = 17. The combination of high resolution and multiplexing with wavelength coverage extending to relatively blue wavelengths (down to 380\,nm), makes HRMOS a spectrograph that will push the boundaries of our knowledge and that is envisioned as a workhorse instrument in the future. The science cases presented in this White Paper include topics and ideas developed by the Core Science Team with the contributions from the astronomical community, also through the wide participation in the first HRMOS Workshop (https://indico.ict.inaf.it/event/1547/) that took place in Firenze (Italy) in October 2021. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: 88 pages, 39 figures. Comments and expressions of interest are welcome by contacting members of the Core Science Team

Orbit decay in encounters between anisotropic spherical galaxies of equal mass

Authors: Lucas Saleh, Joshua E. Barnes

Abstract: We investigate the effect of radial anisotropy on the rate of orbit decay in parabolic encounters of identical spherical galaxies. Our galaxy models have Hernquist density profiles and Osipkov--Merritt velocity distributions. We find that radially anisotropic models merge in as little as half the time of their isotropic counterparts. Anisotropic models are more susceptible to tidal deformation; th… ▽ More We investigate the effect of radial anisotropy on the rate of orbit decay in parabolic encounters of identical spherical galaxies. Our galaxy models have Hernquist density profiles and Osipkov--Merritt velocity distributions. We find that radially anisotropic models merge in as little as half the time of their isotropic counterparts. Anisotropic models are more susceptible to tidal deformation; this accelerates the transfer of orbital angular momentum to internal degrees of freedom. Even during the initial approach, the anisotropic models become more distorted, and arrive at pericentre already having lost substantial amounts of angular momentum. Our results may have implications for estimates of merger rates and persistence of tidal tails. △ Less

Submitted 30 November, 2023; originally announced December 2023.

Comments: 15 pages, 16 figures. Accepted for publication in MNRAS

Combing the Brown Dwarf Desert with Gaia DR3

Authors: Adam T. Stevenson, Carole A. Haswell, John R. Barnes, Joanna K. Barstow

Abstract: We have created an up-to-date catalogue of 214 brown dwarfs (BDs) in binaries with $P < 10^4$ d. This allows us to examine the population statistics of the brown dwarf desert. We searched $\textit{Gaia}$ DR3 NSS results for orbital inclinations of BD candidates, deriving 12 new masses. Three remain as desert BDs whereas nine candidates are found to be low-mass stars. We improved the RV solutions f… ▽ More We have created an up-to-date catalogue of 214 brown dwarfs (BDs) in binaries with $P < 10^4$ d. This allows us to examine the population statistics of the brown dwarf desert. We searched $\textit{Gaia}$ DR3 NSS results for orbital inclinations of BD candidates, deriving 12 new masses. Three remain as desert BDs whereas nine candidates are found to be low-mass stars. We improved the RV solutions for three previously studied BD candidates. A further 19 BD masses with periods less than $\sim$ 1200 d were identified in the DR3 $\texttt{binary_masses}$ database. We confirm a valley in the mass distribution with a minimum around 30-35 M$_{\textrm{jup}}$, and find that periods $<$ 100 d are still under-populated in comparison to longer periods. The updated mass and eccentricity distribution of BDs still shows a marginally statistically significant split into high- and low-mass BD populations. This hints at two different parent distributions, and two potential origins $-$ either akin to planetary formation, or stellar. There are no low eccentricity BDs at periods around 100 d. The mass-metallicity distribution of BDs indicates that core accretion is not the dominant formation mechanism for BDs as they do not follow the same trends that giant exoplanets do with metallicity. We identify a diagonal envelope bounding the Gaia BDs in the mass-period plane due to the detection thresholds of the currently available NSS solutions from 34 months of data. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: 16 pages, 11 figures. Accepted for publication in MNRAS

JWST Observations of the Extraordinary GRB 221009A Reveal an Ordinary Supernova Without Signs of $r$-Process Enrichment in a Low-Metallicity Galaxy

Authors: Peter K. Blanchard, V. Ashley Villar, Ryan Chornock, Tanmoy Laskar, Yijia Li, Joel Leja, Justin Pierel, Edo Berger, Raffaella Margutti, Kate D. Alexander, Jennifer Barnes, Yvette Cendes, Tarraneh Eftekhari, Daniel Kasen, Natalie LeBaron, Brian D. Metzger, James Muzerolle Page, Armin Rest, Huei Sears, Daniel M. Siegel, S. Karthik Yadavalli

Abstract: Identifying the astrophysical sites of the $r$-process, one of the primary mechanisms by which heavy elements are formed, is a key goal of modern astrophysics. The discovery of the brightest gamma-ray burst of all time, GRB 221009A, at a relatively nearby redshift, presented the first opportunity to spectroscopically test the idea that $r$-process elements are produced following the collapse of ra… ▽ More Identifying the astrophysical sites of the $r$-process, one of the primary mechanisms by which heavy elements are formed, is a key goal of modern astrophysics. The discovery of the brightest gamma-ray burst of all time, GRB 221009A, at a relatively nearby redshift, presented the first opportunity to spectroscopically test the idea that $r$-process elements are produced following the collapse of rapidly rotating massive stars. Here we present spectroscopic and photometric $\textit{James Webb Space Telescope}$ (JWST) observations of GRB 221009A obtained $+168$ and $+170$ rest-frame days after the initial gamma-ray trigger, and demonstrate they are well-described by a supernova (SN) and power-law afterglow, with no evidence for an additional component from $r$-process emission, and that the SN component strongly resembles the near-infrared spectra of previous SNe, including SN 1998bw. We further find that the SN associated with GRB 221009A is slightly fainter than the expected brightness of SN 1998bw at this phase, concluding that the SN is therefore not an unusual GRB-SN. We infer a nickel mass of $\approx0.09$ M$_{\odot}$, consistent with the lack of an obvious SN detection in the early-time data. We find that the host galaxy of GRB 221009A has a very low metallicity of $\approx0.12$ Z$_{\odot}$ and our resolved host spectrum shows that GRB 221009A occurred in a unique environment in its host characterized by strong H$_2$ emission lines consistent with recent star formation, which may hint at environmental factors being responsible for its extreme energetics. △ Less

Submitted 27 August, 2023; originally announced August 2023.

Comments: 32 pages, 14 Figures, Submitted to Nature Astronomy

A large topographic feature on the surface of the trans-Neptunian object (307261) 2002 MS$_4$ measured from stellar occultations

Authors: F. L. Rommel, F. Braga-Ribas, J. L. Ortiz, B. Sicardy, P. Santos-Sanz, J. Desmars, J. I. B. Camargo, R. Vieira-Martins, M. Assafin, B. E. Morgado, R. C. Boufleur, G. Benedetti-Rossi, A. R. Gomes-Júnior, E. Fernández-Valenzuela, B. J. Holler, D. Souami, R. Duffard, G. Margoti, M. Vara-Lubiano, J. Lecacheux, J. L. Plouvier, N. Morales, A. Maury, J. Fabrega, P. Ceravolo , et al. (179 additional authors not shown)

Abstract: This work aims at constraining the size, shape, and geometric albedo of the dwarf planet candidate 2002 MS4 through the analysis of nine stellar occultation events. Using multichord detection, we also studied the object's topography by analyzing the obtained limb and the residuals between observed chords and the best-fitted ellipse. We predicted and organized the observational campaigns of nine st… ▽ More This work aims at constraining the size, shape, and geometric albedo of the dwarf planet candidate 2002 MS4 through the analysis of nine stellar occultation events. Using multichord detection, we also studied the object's topography by analyzing the obtained limb and the residuals between observed chords and the best-fitted ellipse. We predicted and organized the observational campaigns of nine stellar occultations by 2002 MS4 between 2019 and 2022, resulting in two single-chord events, four double-chord detections, and three events with three to up to sixty-one positive chords. Using 13 selected chords from the 8 August 2020 event, we determined the global elliptical limb of 2002 MS4. The best-fitted ellipse, combined with the object's rotational information from the literature, constrains the object's size, shape, and albedo. Additionally, we developed a new method to characterize topography features on the object's limb. The global limb has a semi-major axis of 412 $\pm$ 10 km, a semi-minor axis of 385 $\pm$ 17 km, and the position angle of the minor axis is 121 $^\circ$ $\pm$ 16$^\circ$. From this instantaneous limb, we obtained 2002 MS4's geometric albedo and the projected area-equivalent diameter. Significant deviations from the fitted ellipse in the northernmost limb are detected from multiple sites highlighting three distinct topographic features: one 11 km depth depression followed by a 25$^{+4}_{-5}$ km height elevation next to a crater-like depression with an extension of 322 $\pm$ 39 km and 45.1 $\pm$ 1.5 km deep. Our results present an object that is $\approx$138 km smaller in diameter than derived from thermal data, possibly indicating the presence of a so-far unknown satellite. However, within the error bars, the geometric albedo in the V-band agrees with the results published in the literature, even with the radiometric-derived albedo. △ Less

Submitted 23 August, 2023; v1 submitted 15 August, 2023; originally announced August 2023.

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

DMPP-4: Candidate sub-Neptune mass planets orbiting a naked-eye star

Authors: J. R. Barnes, M. R. Standing, C. A. Haswell, D. Staab, J. P. J. Doherty, M. Waller-Bridge, L. Fossati, M. Soto, G. Anglada-Escudé, J. Llama, C. McCune, F. W. Lewis

Abstract: We present radial velocity measurements of the very bright ($V\sim5.7$) nearby F star, DMPP-4 (HD 184960). The anomalously low Ca II H&K emission suggests mass loss from planets orbiting a low activity host star. Periodic radial velocity variability with $\sim 10$ ms$^{-1}$ amplitude is found to persist over a $>4$ year timescale. Although the non-simultaneous photometric variability in four TESS… ▽ More We present radial velocity measurements of the very bright ($V\sim5.7$) nearby F star, DMPP-4 (HD 184960). The anomalously low Ca II H&K emission suggests mass loss from planets orbiting a low activity host star. Periodic radial velocity variability with $\sim 10$ ms$^{-1}$ amplitude is found to persist over a $>4$ year timescale. Although the non-simultaneous photometric variability in four TESS sectors supports the view of an inactive star, we identify periodic photometric signals and also find spectroscopic evidence for stellar activity. We used a posterior sampling algorithm that includes the number of Keplerian signals, $N_\textrm{p}$, as a free parameter to test and compare (1) purely Keplerian models (2) a Keplerian model with linear activity correlation and (3) Keplerian models with Gaussian processes. A preferred model, with one Keplerian and quasi-periodic Gaussian process indicates a planet with a period of $P_\textrm{b} = 3.4982^{+0.0015}_{-0.0027}$ d and corresponding minimum mass of $m_\textrm{b}\,\textrm{sin}\,i = 12.2^{+1.8}_{-1.9}$ M$_\oplus$. Without further high time resolution observations over a longer timescale, we cannot definitively rule out the purely Keplerian model with 2 candidates planets with $P_\textrm{b} = 2.4570^{+0.0026}_{-0.0462}$ d, minimum mass $m_\textrm{b}\,\textrm{sin}\,i = 8.0^{+1.1}_{-1.5}$ M$_\oplus$ and $P_\textrm{c} = 5.4196^{+0.6766}_{-0.0030}$ d and corresponding minimum mass of $m_\textrm{b}\,\textrm{sin}\,i = 12.2^{+1.4}_{-1.6}$ M$_\oplus$. The candidate planets lie in the region below the lower-envelope of the Neptune Desert. Continued mass loss may originate from the highly irradiated planets or from an as yet undetected body in the system. △ Less

Submitted 24 July, 2023; v1 submitted 13 July, 2023; originally announced July 2023.

Comments: 19 pages, 11 figures. Accepted for publication in MNRAS

Can the orbital distribution of Neptune's 3:2 mean motion resonance result from stability sculpting?

Authors: Sricharan Balaji, Nihaal Zaveri, Nanae Hayashi, Arcelia Hermosillo Ruiz, Jackson Barnes, Ruth Murray-Clay, Kathryn Volk, Jake Gerhardt, Zain Syed

Abstract: We explore a simplified model of the outcome of an early outer Solar System gravitational upheaval during which objects were captured into Neptune's 3:2 mean motion resonance via scattering rather than smooth planetary migration. We use N-body simulations containing the Sun, the four giant planets, and test particles in the 3:2 resonance to determine whether long-term stability sculpting over 4.5… ▽ More We explore a simplified model of the outcome of an early outer Solar System gravitational upheaval during which objects were captured into Neptune's 3:2 mean motion resonance via scattering rather than smooth planetary migration. We use N-body simulations containing the Sun, the four giant planets, and test particles in the 3:2 resonance to determine whether long-term stability sculpting over 4.5 Gyr can reproduce the observed 3:2 resonant population from an initially randomly scattered 3:2 population. After passing our simulated 3:2 resonant objects through a survey simulator, we find that the semimajor axis (a) and eccentricity (e) distributions are consistent with the observational data (assuming an absolute magnitude distribution constrained by prior studies), suggesting that these could be a result of stability sculpting. However, the inclination (i) distribution cannot be produced be stability sculpting and thus must result from a distinct process that excited the inclinations. Our simulations modestly under-predict the number of objects with high libration amplitudes (Aφ), possibly because we do not model transient sticking. Finally, our model under-populates the Kozai subresonance compared to both observations and to smooth migration models. Future work is needed to determine whether smooth migration occurring as Neptune's eccentricity damped to its current value can resolve this discrepancy. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: 15 pages, 13 figures

DMPP-3: confirmation of short-period S-type planet(s) in a compact eccentric binary star system, and warnings about long-period RV planet detections

Authors: Adam T. Stevenson, Carole A. Haswell, John R. Barnes, Joanna K. Barstow, Zachary O. B. Ross

Abstract: We present additional HARPS radial velocity observations of the highly eccentric ($e \sim 0.6$) binary system DMPP-3AB, which comprises a K0V primary and a low-mass companion at the hydrogen burning limit. The binary has a $507$ d orbital period and a $1.2$ au semi-major axis. The primary component harbours a known $2.2$ M$_{\oplus}$ planet, DMPP-3A b, with a $6.67$ day orbit. New HARPS measuremen… ▽ More We present additional HARPS radial velocity observations of the highly eccentric ($e \sim 0.6$) binary system DMPP-3AB, which comprises a K0V primary and a low-mass companion at the hydrogen burning limit. The binary has a $507$ d orbital period and a $1.2$ au semi-major axis. The primary component harbours a known $2.2$ M$_{\oplus}$ planet, DMPP-3A b, with a $6.67$ day orbit. New HARPS measurements constrain periastron passage for the binary orbit and add further integrity to previously derived solutions for both companion and planet orbits. Gaia astrometry independently confirms the binary orbit, and establishes the inclination of the binary is $63.89 \pm 0.78 ^{\circ}$. We performed dynamical simulations which establish that the previously identified $\sim800$ d RV signal cannot be attributed to an orbiting body. The additional observations, a deviation from strict periodicity, and our new analyses of activity indicators suggest the $\sim800$ d signal is caused by stellar activity. We conclude that there may be long period planet 'detections' in other systems which are similar misinterpreted stellar activity artefacts. Without the unusual eccentric binary companion to the planet-hosting star we could have accepted the $\sim800$ d signal as a probable planet. Further monitoring of DMPP-3 will reveal which signatures can be used to most efficiently identify these imposters. We also report a threshold detection (0.2 per cent FAP) of a $\sim2.26$ d periodicity in the RVs, potentially attributed to an Earth-mass S-type planet interior to DMPP-3A b. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: 20 pages, 15 figures. Accepted for publication in MNRAS

Hydrodynamic mixing of accretion disk outflows in collapsars: implications for r-process signatures

Authors: Jennifer Barnes, Paul C. Duffell

Abstract: The astrophysical environments capable of triggering heavy-element synthesis via rapid neutron capture (the r-process) remain uncertain. While binary neutron star mergers (NSMs) are known to forge r-process elements, certain rare supernovae (SNe) have been theorized to supplement, or even dominate, r-production by NSMs. However, the most direct evidence for such SNe, unusual reddening of the emiss… ▽ More The astrophysical environments capable of triggering heavy-element synthesis via rapid neutron capture (the r-process) remain uncertain. While binary neutron star mergers (NSMs) are known to forge r-process elements, certain rare supernovae (SNe) have been theorized to supplement, or even dominate, r-production by NSMs. However, the most direct evidence for such SNe, unusual reddening of the emission caused by the high opacities of r-process elements, has not been observed. Recent work identified the distribution of r-process material within the SN ejecta as a key predictor of the ease with which signals associated with r-process enrichment could be discerned. Though this distribution results from hydrodynamic processes at play during the SN explosion, thus far it has been treated only in a parameterized way. We use hydrodynamic simulations to model how disk winds, the alleged locus of r-production in rare SNe, mix with initially r-process-free ejecta. We study mixing as a function of the wind mass and duration and of the initial SN explosion energy, and find that it increases with the first two of these and decreases with the third. This suggests that SNe accompanying the longest long-duration gamma-ray bursts are promising places to search for signs of r-process enrichment. We use semianalytic radiation transport to connect hydrodynamics to electromagnetic observables, allowing us to assess the mixing level at which the presence of r-process material can be diagnosed from SN light curves. Analytic arguments constructed atop this foundation imply that a wind-driven r-process-enriched SN model is unlikely to explain standard energetic SNe. △ Less

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

Comments: V2 accepted to ApJ. Absolute banger. The authors encourage the widest possible adoption of the term "collapsaR-process.''

Superheavy Elements in Kilonovae

Authors: Erika M. Holmbeck, Jennifer Barnes, Kelsey A. Lund, Trevor M. Sprouse, G. C. McLaughlin, Matthew R. Mumpower

Abstract: As LIGO-Virgo-KAGRA enters its fourth observing run, a new opportunity to search for electromagnetic counterparts of compact object mergers will also begin. The light curves and spectra from the first "kilonova" associated with a binary neutron star binary (NSM) suggests that these sites are hosts of the rapid neutron capture ("$r$") process. However, it is unknown just how robust elemental produc… ▽ More As LIGO-Virgo-KAGRA enters its fourth observing run, a new opportunity to search for electromagnetic counterparts of compact object mergers will also begin. The light curves and spectra from the first "kilonova" associated with a binary neutron star binary (NSM) suggests that these sites are hosts of the rapid neutron capture ("$r$") process. However, it is unknown just how robust elemental production can be in mergers. Identifying signposts of the production of particular nuclei is critical for fully understanding merger-driven heavy-element synthesis. In this study, we investigate the properties of very neutron rich nuclei for which superheavy elements ($Z\geq 104$) can be produced in NSMs and whether they can similarly imprint a unique signature on kilonova light-curve evolution. A superheavy-element signature in kilonovae represents a route to establishing a lower limit on heavy-element production in NSMs as well as possibly being the first evidence of superheavy element synthesis in nature. Favorable NSMs conditions yield a mass fraction of superheavy elements is $X_{Z\geq 104}\approx 3\times 10^{-2}$ at 7.5 hours post-merger. With this mass fraction of superheavy elements, we find that the component of kilonova light curves possibly containing superheavy elements may appear similar to those arising from lanthanide-poor ejecta. Therefore, photometric characterizations of superheavy-element rich kilonova may possibly misidentify them as lanthanide-poor events. △ Less

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

Comments: 8 pages, 5 figures, 4 data tables

Collapsars as Sites of r-process Nucleosynthesis: Systematic Near-Infrared Follow-up of Type Ic-BL Supernovae

Authors: Shreya Anand, Jennifer Barnes, Sheng Yang, Mansi M. Kasliwal, Michael W. Coughlin, Jesper Sollerman, Kishalay De, Christoffer Fremling, Alessandra Corsi, Anna Y. Q. Ho, Arvind Balasubramanian, Conor Omand, Gokul P. Srinivasaragavan, S. Bradley Cenko, Tomas Ahumada, Igor Andreoni, Aishwarya Dahiwale, Kaustav Kashyap Das, Jacob Jencson, Viraj Karambelkar, Harsh Kumar, Brian D. Metzger, Daniel Perley, Nikhil Sarin, Tassilo Schweyer , et al. (19 additional authors not shown)

Abstract: One of the open questions following the discovery of GW170817 is whether neutron star mergers are the only astrophysical sites capable of producing $r$-process elements. Simulations have shown that 0.01-0.1M$_\odot$ of $r$-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both neutron star… ▽ More One of the open questions following the discovery of GW170817 is whether neutron star mergers are the only astrophysical sites capable of producing $r$-process elements. Simulations have shown that 0.01-0.1M$_\odot$ of $r$-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both neutron star mergers and collapsing massive stars associated with long-duration gamma-ray bursts (collapsars). The hallmark signature of $r$-process nucleosynthesis in the binary neutron star merger GW170817 was its long-lasting near-infrared emission, thus motivating a systematic photometric study of the light curves of broadlined stripped-envelope (Ic-BL) supernovae (SNe) associated with collapsars. We present the first systematic study of 25 SNe Ic-BL -- including 18 observed with the Zwicky Transient Facility and 7 from the literature -- in the optical/near-infrared bands to determine what quantity of $r$-process material, if any, is synthesized in these explosions. Using semi-analytic models designed to account for $r$-process production in SNe Ic-BL, we perform light curve fitting to derive constraints on the $r$-process mass for these SNe. We also perform independent light curve fits to models without $r$-process. We find that the $r$-process-free models are a better fit to the light curves of the objects in our sample. Thus we find no compelling evidence of $r$-process enrichment in any of our objects. Further high-cadence infrared photometric studies and nebular spectroscopic analysis would be sensitive to smaller quantities of $r$-process ejecta mass or indicate whether all collapsars are completely devoid of $r$-process nucleosynthesis. △ Less

Submitted 12 February, 2024; v1 submitted 17 February, 2023; originally announced February 2023.

Comments: 37 pages, 8 figures, 4 tables

Journal ref: ApJ, 962, 68 (2024)

A Machine Learning approach for correcting radial velocities using physical observables

Authors: M. Perger, G. Anglada-Escudé, D. Baroch, M. Lafarga, I. Ribas, J. C. Morales, E. Herrero, P. J. Amado, J. R. Barnes, J. A. Caballero, S. V. Jeffers, A. Quirrenbach, A. Reiners

Abstract: Precision radial velocity (RV) measurements continue to be a key tool to detect and characterise extrasolar planets. While instrumental precision keeps improving, stellar activity remains a barrier to obtain reliable measurements below 1-2 m/s accuracy. Using simulations and real data, we investigate the capabilities of a Deep Neural Network approach to produce activity free Doppler measurements o… ▽ More Precision radial velocity (RV) measurements continue to be a key tool to detect and characterise extrasolar planets. While instrumental precision keeps improving, stellar activity remains a barrier to obtain reliable measurements below 1-2 m/s accuracy. Using simulations and real data, we investigate the capabilities of a Deep Neural Network approach to produce activity free Doppler measurements of stars. As case studies we use observations of two known stars (Eps Eridani and AUMicroscopii), both with clear signals of activity induced RV variability. Synthetic data using the starsim code are generated for the observables (inputs) and the resulting RV signal (labels), and used to train a Deep Neural Network algorithm. We identify an architecture consisting of convolutional and fully connected layers that is adequate to the task. The indices investigated are mean line-profile parameters (width, bisector, contrast) and multi-band photometry. We demonstrate that the RV-independent approach can drastically reduce spurious Doppler variability from known physical effects such as spots, rotation and convective blueshift. We identify the combinations of activity indices with most predictive power. When applied to real observations, we observe a good match of the correction with the observed variability, but we also find that the noise reduction is not as good as in the simulations, probably due to the lack of detail in the simulated physics. We demonstrate that a model-driven machine learning approach is sufficient to clean Doppler signals from activity induced variability for well known physical effects. There are dozens of known activity related observables whose inversion power remains unexplored indicating that the use of additional indicators, more complete models, and more observations with optimised sampling strategies can lead to significant improvements in our detrending capabilities. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Journal ref: A&A 672, A118 (2023)

Science Priorities for the Extraction of the Solid MSR Samples from their Sample Tubes

Authors: N. Dauphas, S. S. Russell, D. Beaty, F. Thiessen, J. Barnes, L. Bonal, J. Bridges, T. Bristow, J. Eiler, L. Ferriere, T. Fornaro, J. Gattacceca, B. Hoffman, E. J. Javaux, T. Kleine, H. Y. McSween, M. Prasad, L. Rampe, M. Schmidt, B. Schoene, K. L. Siebach, J. Stern, N. Tosca

Abstract: Preservation of the chemical and structural integrity of samples that will be brought back from Mars is paramount to achieving the scientific objectives of MSR. Given our knowledge of the nature of the samples retrieved at Jezero by Perseverance, at least two options need to be tested for opening the sample tubes: (1) One or two radial cuts at the end of the tube to slide the sample out. (2) Two r… ▽ More Preservation of the chemical and structural integrity of samples that will be brought back from Mars is paramount to achieving the scientific objectives of MSR. Given our knowledge of the nature of the samples retrieved at Jezero by Perseverance, at least two options need to be tested for opening the sample tubes: (1) One or two radial cuts at the end of the tube to slide the sample out. (2) Two radial cuts at the ends of the tube and two longitudinal cuts to lift the upper half of the tube and access the sample. Strategy 1 will likely minimize contamination but incurs the risk of affecting the physical integrity of weakly consolidated samples. Strategy 2 will be optimal for preserving the physical integrity of the samples but increases the risk of contamination and mishandling of the sample as more manipulations and additional equipment will be needed. A flexible approach to opening the sample tubes is therefore required, and several options need to be available, depending on the nature of the rock samples returned. Both opening strategies 1 and 2 may need to be available when the samples are returned to handle different sample types (e.g., loosely bound sediments vs. indurated magmatic rocks). This question should be revisited after engineering tests are performed on analogue samples. The MSR sample tubes will have to be opened under stringent BSL4 conditions and this aspect needs to be integrated into the planning. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: 8 pages, 3 figures, 1 table, report NASA-ESA Mars Rock Team Report

SOFIA and ALMA Investigate Magnetic Fields and Gas Structures in Massive Star Formation: The Case of the Masquerading Monster in BYF 73

Authors: Peter J. Barnes, Stuart D. Ryder, Giles Novak, Richard M. Crutcher, Laura M. Fissel, Rebecca L. Pitts, William J. Schap III

Abstract: We present SOFIA+ALMA continuum and spectral-line polarisation data on the massive molecular cloud BYF 73, revealing important details about the magnetic field morphology, gas structures, and energetics in this unusual massive star formation laboratory. The 154$μ$m HAWC+ polarisation map finds a highly organised magnetic field in the densest, inner 0.55$\times$0.40 pc portion of the cloud, compare… ▽ More We present SOFIA+ALMA continuum and spectral-line polarisation data on the massive molecular cloud BYF 73, revealing important details about the magnetic field morphology, gas structures, and energetics in this unusual massive star formation laboratory. The 154$μ$m HAWC+ polarisation map finds a highly organised magnetic field in the densest, inner 0.55$\times$0.40 pc portion of the cloud, compared to an unremarkable morphology in the cloud's outer layers. The 3mm continuum ALMA polarisation data reveal several more structures in the inner domain, including a pc-long, $\sim$500 M$_{\odot}$ "Streamer" around the central massive protostellar object MIR 2, with magnetic fields mostly parallel to the east-west Streamer but oriented north-south across MIR 2. The magnetic field orientation changes from mostly parallel to the column density structures to mostly perpendicular, at thresholds $N_{\rm crit}$ = 6.6$\times$10$^{26}$ m$^{-2}$, $n_{\rm crit}$ = 2.5$\times$10$^{11}$ m$^{-3}$, and $B_{\rm crit}$ = 42$\pm$7 nT. ALMA also mapped Goldreich-Kylafis polarisation in $^{12}$CO across the cloud, which traces in both total intensity and polarised flux, a powerful bipolar outflow from MIR 2 that interacts strongly with the Streamer. The magnetic field is also strongly aligned along the outflow direction; energetically, it may dominate the outflow near MIR 2, comprising rare evidence for a magnetocentrifugal origin to such outflows. A portion of the Streamer may be in Keplerian rotation around MIR 2, implying a gravitating mass 1350$\pm$50 M$_{\odot}$ for the protostar+disk+envelope; alternatively, these kinematics can be explained by gas in free fall towards a 950$\pm$35 M$_{\odot}$ object. The high accretion rate onto MIR 2 apparently occurs through the Streamer/disk, and could account for $\sim$33% of MIR 2's total luminosity via gravitational energy release. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Comments: 33 pages, 32 figures, accepted by ApJ. Line-Integral Convolution (LIC) images and movie versions of Figures 3b, 7, and 29 are available at https://gemelli.spacescience.org/~pbarnes/research/champ/papers/

Bayesian model selection for GRB 211211A through multi-wavelength analyses

Authors: Nina Kunert, Sarah Antier, Vsevolod Nedora, Mattia Bulla, Peter T. H. Pang, Shreya Anand, Michael Coughlin, Ingo Tews, Jennifer Barnes, Thomas Hussenot-Desenonges, Brian Healy, Theophile Jegou du Laz, Meili Pilloix, Weizmann Kiendrebeogo, Tim Dietrich

Abstract: Although GRB 211211A is one of the closest gamma-ray bursts (GRBs), its classification is challenging because of its partially inconclusive electromagnetic signatures. In this paper, we investigate four different astrophysical scenarios as possible progenitors for GRB~211211A: a binary neutron-star merger, a black-hole--neutron-star merger, a core-collapse supernova, and an r-process enriched core… ▽ More Although GRB 211211A is one of the closest gamma-ray bursts (GRBs), its classification is challenging because of its partially inconclusive electromagnetic signatures. In this paper, we investigate four different astrophysical scenarios as possible progenitors for GRB~211211A: a binary neutron-star merger, a black-hole--neutron-star merger, a core-collapse supernova, and an r-process enriched core collapse of a rapidly rotating massive star (a collapsar). We perform a large set of Bayesian multi-wavelength analyses based on different models describing these scenarios and priors to investigate which astrophysical scenarios and processes might be related to GRB~211211A. Our analysis supports previous studies in which the presence of an additional component, likely related to $r$-process nucleosynthesis, is required to explain the observed light curves of GRB~211211A, as it can not solely be explained as a GRB afterglow. Fixing the distance to about $350~\rm Mpc$, namely the distance of the possible host galaxy SDSS J140910.47+275320.8, we find a statistical preference for a binary neutron-star merger scenario. △ Less

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

Comments: 14 pages, 6 figures

Report number: LA-UR-22-32364

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 2, January 2024, Pages 3900-3911

A collapsar origin for GRB 211211A is (just barely) possible

Authors: Jennifer Barnes, Brian D. Metzger

Abstract: Gamma-ray bursts (GRBs) have historically been divided into two classes. Short-duration GRBs are associated with binary neutron-star mergers (NSMs), while long-duration bursts are connected to a subset of core-collapse supernovae (SNe). GRB 211211A recently made headlines as the first long-duration burst purportedly generated by an NSM. The evidence for an NSM origin was excess optical and near-in… ▽ More Gamma-ray bursts (GRBs) have historically been divided into two classes. Short-duration GRBs are associated with binary neutron-star mergers (NSMs), while long-duration bursts are connected to a subset of core-collapse supernovae (SNe). GRB 211211A recently made headlines as the first long-duration burst purportedly generated by an NSM. The evidence for an NSM origin was excess optical and near-infrared emission consistent with the kilonova observed after the gravitational wave-detected NSM GW170817. Kilonovae derive their unique electromagnetic signatures from the properties of the heavy elements synthesized by rapid neutron capture (the r-process) following the merger. Recent simulations suggest that the "collapsar" SNe that trigger long GRBs may also produce r-process elements. While observations of GRB 211211A and its afterglow ruled out an SN typical of those that follow long GRBs, an unusual collapsar could explain both the duration of GRB 211211A and the r-process-powered excess in its afterglow. We use semianalytic radiation transport modeling to evaluate low-mass collapsars as the progenitors of GRB 211211A-like events. We compare a suite of collapsar models to the afterglow-subtracted emission that followed GRB 211211A, and find the best agreement for models with high kinetic energies and an unexpected pattern of Ni-56 enrichment. We discuss how core-collapse explosions could produce such ejecta, and how distinct our predictions are from those generated by more straightforward kilonova models. We also show that radio observations can distinguish between kilonovae and the more massive collapsar ejecta we consider here. △ Less

Submitted 3 January, 2023; originally announced January 2023.

Comments: Submitted to ApJ. Comments welcome

A Method to Achieve High Dynamic Range in a CMOS Image Sensor Using Interleaved Row Readout

Authors: Thomas Wocial, Konstantin D. Stefanov, William E. Martin, John R. Barnes, Hugh R. A. Jones

Abstract: We present a readout scheme for CMOS image sensors that can be used to achieve arbitrarily high dynamic range (HDR) in principle. The linear full well capacity (LFWC) in high signal regions was extended 50 times from 20 ke$^{-}$ to 984 ke$^{-}$ via an interlaced row-wise readout order, whilst the noise floor remained unchanged in low signal regions, resulting in a 34 dB increase in DR. The peak si… ▽ More We present a readout scheme for CMOS image sensors that can be used to achieve arbitrarily high dynamic range (HDR) in principle. The linear full well capacity (LFWC) in high signal regions was extended 50 times from 20 ke$^{-}$ to 984 ke$^{-}$ via an interlaced row-wise readout order, whilst the noise floor remained unchanged in low signal regions, resulting in a 34 dB increase in DR. The peak signal-to-noise ratio (PSNR) is increased in a continuous fashion from 43 dB to 60 dB. This was achieved by summing user-selected rows which were read out multiple times. Centroiding uncertainties were lowered when template-fitting a projected pattern, compared to the standard readout scheme. Example applications are aimed at scientific imaging due to the linearity and PSNR increase. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: 14 pages, 10 figures. Accepted to IEEE Sensors on 14/09/2022, available as a preprint as of 10/10/2022

Jupiter and Saturn as Spectral Analogs for Extrasolar Gas Giants and Brown Dwarfs

Authors: Daniel J. Coulter, Jason W. Barnes, Jonathan J. Fortney

Abstract: With the advent of direct imaging spectroscopy, the number of spectra from brown dwarfs and extrasolar gas giants is growing rapidly. Many brown dwarfs and extrasolar gas giants exhibit spectroscopic and photometric variability, which is likely the result of weather patterns. However, for the foreseeable future, point-source observations will be the only viable method to extract brown dwarf and ex… ▽ More With the advent of direct imaging spectroscopy, the number of spectra from brown dwarfs and extrasolar gas giants is growing rapidly. Many brown dwarfs and extrasolar gas giants exhibit spectroscopic and photometric variability, which is likely the result of weather patterns. However, for the foreseeable future, point-source observations will be the only viable method to extract brown dwarf and exoplanet spectra. Models have been able to reproduce the observed variability, but ground truth observations are required to verify their results. To that end, we provide visual and near-infrared spectra of Jupiter and Saturn obtained from the \emph{Cassini} VIMS instrument. We disk-integrate the VIMS spectral cubes to simulate the spectra of Jupiter and Saturn as if they were directly imaged exoplanets or brown dwarfs. We present six empirical disk-integrated spectra for both Jupiter and Saturn with phase coverage of $1.7^\circ$ to $133.5^\circ$ and $39.6^\circ$ to $110.2^\circ$, respectively. To understand the constituents of these disk-integrated spectra, we also provide end member (single feature) spectra for permutations of illumination and cloud density, as well as for Saturn's rings. In tandem, these disk-integrated and end member spectra provide the ground truth needed to analyze point source spectra from extrasolar gas giants and brown dwarfs. Lastly, we discuss the impact that icy rings, such as Saturn's, have on disk-integrated spectra and consider the feasibility of inferring the presence of rings from direct imaging spectra. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Detailed stellar activity analysis and modelling of GJ 832: Reassessment of the putative habitable zone planet GJ 832c

Authors: P. Gorrini, N. Astudillo-Defru, S. Dreizler, M. Damasso, R. F. Díaz, X. Bonfils, S. V. Jeffers, J. R. Barnes, F. Del Sordo, J. -M. Almenara, E. Artigau, F. Bouchy, D. Charbonneau, X. Delfosse, R. Doyon, P. Figueira, T. Forveille, C. A. Haswell, M. J. López-González, C. Melo, R. E. Mennickent, G. Gaisné, N. Morales, F. Murgas, F. Pepe , et al. (5 additional authors not shown)

Abstract: Context. Gliese 832 (GJ 832) is an M2V star hosting a massive planet on a decade-long orbit, GJ 832b, discovered by radial velocity (RV). Later, a super Earth or mini-Neptune orbiting within the stellar habitable zone was reported (GJ 832c). The recently determined stellar rotation period (45.7 $\pm$ 9.3 days) is close to the orbital period of putative planet c (35.68 $\pm$ 0.03 days). Aims. We… ▽ More Context. Gliese 832 (GJ 832) is an M2V star hosting a massive planet on a decade-long orbit, GJ 832b, discovered by radial velocity (RV). Later, a super Earth or mini-Neptune orbiting within the stellar habitable zone was reported (GJ 832c). The recently determined stellar rotation period (45.7 $\pm$ 9.3 days) is close to the orbital period of putative planet c (35.68 $\pm$ 0.03 days). Aims. We aim to confirm or dismiss the planetary nature of the RV signature attributed to GJ 832c, by adding 119 new RV data points, new photometric data, and an analysis of the spectroscopic stellar activity indicators. Additionally, we update the orbital parameters of the planetary system and search for additional signals. Methods. We performed a frequency content analysis of the RVs to search for periodic and stable signals. Radial velocity time series were modelled with Keplerians and Gaussian process (GP) regressions alongside activity indicators to subsequently compare them within a Bayesian framework. Results. We updated the stellar rotational period of GJ 832 from activity indicators, obtaining $37.5^{+1.4}_{-1.5}$ days, improving the precision by a factor of 6. The new photometric data are in agreement with this value. We detected an RV signal near 18 days (FAP < 4.6%), which is half of the stellar rotation period. Two Keplerians alone fail at modelling GJ 832b and a second planet with a 35-day orbital period. Moreover, the Bayesian evidence from the GP analysis of the RV data with simultaneous activity indices prefers a model without a second Keplerian, therefore negating the existence of planet c. △ Less

Submitted 2 August, 2022; v1 submitted 15 June, 2022; originally announced June 2022.

Comments: A&A in press

Journal ref: A&A 664, A64 (2022)

Signatures of R-process Enrichment in Supernovae from Collapsars

Authors: Jennifer Barnes, Brian D. Metzger

Abstract: Despite recent progress, the astrophysical channels responsible for rapid neutron capture (r-process) nucleosynthesis remain an unsettled question. Observations of kilonovae following gravitational wave-detected neutron star mergers established mergers as one site of the r-process, but additional sources may be needed to fully explain r-process enrichment in the Universe. One intriguing possibilit… ▽ More Despite recent progress, the astrophysical channels responsible for rapid neutron capture (r-process) nucleosynthesis remain an unsettled question. Observations of kilonovae following gravitational wave-detected neutron star mergers established mergers as one site of the r-process, but additional sources may be needed to fully explain r-process enrichment in the Universe. One intriguing possibility is that rapidly rotating massive stars undergoing core collapse launch r-process-rich outflows off the accretion disks formed from their infalling matter. In this scenario, r-process winds comprise one component of the supernova (SN) ejecta produced by "collapsar" explosions. We present the first systematic study of the effects of r-process enrichment on the emission from collapsar-generated SNe. We semi-analytically model r-process SN emission from explosion out to late times, and determine its distinguishing features. The ease with which r-process SNe can be identified depends on how effectively wind material mixes into the initially r-process-free outer layers of the ejecta. In many cases, enrichment produces a near infrared (NIR) excess that can be detected within ~75 days of explosion. We also discuss optimal targets and observing strategies for testing the r-process collapsar theory, and find that frequent monitoring of optical and NIR emission from high-velocity SNe in the first few months after explosion offers a reasonable chance of success while respecting finite observing resources. Such early identification of r-process collapsar candidates also lays the foundation for nebular-phase spectroscopic follow-up in the near- and mid-infrared, for example with the James Webb Space Telescope. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: Submitted to ApJ. Comments welcome

Horizons: Nuclear Astrophysics in the 2020s and Beyond

Authors: H. Schatz, A. D. Becerril Reyes, A. Best, E. F. Brown, K. Chatziioannou, K. A. Chipps, C. M. Deibel, R. Ezzeddine, D. K. Galloway, C. J. Hansen, F. Herwig, A. P. Ji, M. Lugaro, Z. Meisel, D. Norman, J. S. Read, L. F. Roberts, A. Spyrou, I. Tews, F. X. Timmes, C. Travaglio, N. Vassh, C. Abia, P. Adsley, S. Agarwal , et al. (140 additional authors not shown)

Abstract: Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit… ▽ More Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 96 pages. Submitted to Journal of Physics G

Report number: LA-UR-22-23997

The CARMENES search for exoplanets around M dwarfs: Benchmarking the impact of activity in high-precision radial velocity measurements

Authors: S. V. Jeffers, J. R. Barnes, P. Scheofer, A. Quirrenbach, M. Zechmeister, P. J. Amado, J. A. Caballero, M. Fernandez, E. Rodriguez, I. Ribas, A. Reiners, C. Cardona Guillen, C. Cifuentes, S. Czesla, A. P. Hatzes, M. Kurster, D. Montes, J. C. Morales, S. Pedraz, S. Sadegi

Abstract: Current exoplanet surveys using the radial velocity (RV) technique are targeting M dwarfs because any habitable zone terrestrial-mass planets will induce a high RV and orbit on shorter periods than for more massive stars. One of the main caveats is that M dwarfs show a wide range of activity levels from inactive to very active, which can induce an asymmetry in the line profiles and, consequently,… ▽ More Current exoplanet surveys using the radial velocity (RV) technique are targeting M dwarfs because any habitable zone terrestrial-mass planets will induce a high RV and orbit on shorter periods than for more massive stars. One of the main caveats is that M dwarfs show a wide range of activity levels from inactive to very active, which can induce an asymmetry in the line profiles and, consequently, a spurious RV measurement. We aim to benchmark the impact of stellar activity on high-precision RV measurements using regular-cadence CARMENES visible and near-infrared observations of the active M3.5 dwarf EV Lac. We used the newly developed technique of low-resolution Doppler imaging to determine the centre-of-light, or spot-induced RV component, for eight observational epochs. We confirm a statistically significant and strong correlation between the independently measured centre-of-light and the chromatic index, which is a measure of the amplitude variation with wavelength of the RVs. We also find circular closed-loop relations of several activity indices with RV for a subset of data that covers only several rotation periods. We also investigate the implications of large phase gaps in the periodograms of activity indicators. Finally, by removing the spot-induced RV component we improve the planet-mass sensitivity by a factor of at least three. We conclude that for active M stars, a regular-cadence observing strategy is the most efficient way to identify and eliminate sources of correlated noise. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: accepted by A&A

Journal ref: A&A 663, A27 (2022)

"Super-Kilonovae" from Massive Collapsars as Signatures of Black-Hole Birth in the Pair-instability Mass Gap

Authors: Daniel M. Siegel, Aman Agarwal, Jennifer Barnes, Brian D. Metzger, Mathieu Renzo, V. Ashley Villar

Abstract: The core collapse of rapidly rotating massive ~10 Msun stars ("collapsars"), and resulting formation of hyper-accreting black holes, are a leading model for the central engines of long-duration gamma-ray bursts (GRB) and promising sources of r-process nucleosynthesis. Here, we explore the signatures of collapsars from progenitors with extremely massive helium cores >130 Msun above the pair-instabi… ▽ More The core collapse of rapidly rotating massive ~10 Msun stars ("collapsars"), and resulting formation of hyper-accreting black holes, are a leading model for the central engines of long-duration gamma-ray bursts (GRB) and promising sources of r-process nucleosynthesis. Here, we explore the signatures of collapsars from progenitors with extremely massive helium cores >130 Msun above the pair-instability mass gap. While rapid collapse to a black hole likely precludes a prompt explosion in these systems, we demonstrate that disk outflows can generate a large quantity (up to >50 Msun) of ejecta, comprised of >5-10 Msun in r-process elements and ~0.1-1 Msun of $^{56}$Ni, expanding at velocities ~0.1c. Radioactive heating of the disk-wind ejecta powers an optical/infrared transient, with a characteristic luminosity $\sim 10^{42}$ erg s$^{-1}$ and spectral peak in the near-infrared (due to the high optical/UV opacities of lanthanide elements) similar to kilonovae from neutron star mergers, but with longer durations $\gtrsim$ 1 month. These "super-kilonovae" (superKNe) herald the birth of massive black holes >60 Msun, which, as a result of disk wind mass-loss, can populate the pair-instability mass gap 'from above' and could potentially create the binary components of GW190521. SuperKNe could be discovered via wide-field surveys such as those planned with the Roman Space Telescope or via late-time infrared follow-up observations of extremely energetic GRBs. Gravitational waves of frequency ~0.1-50 Hz from non-axisymmetric instabilities in self-gravitating massive collapsar disks are potentially detectable by proposed third-generation intermediate and high-frequency observatories at distances up to hundreds of Mpc; in contrast to the "chirp" from binary mergers, the collapsar gravitational-wave signal decreases in frequency as the disk radius grows ("sad trombone"). △ Less

Submitted 4 November, 2021; originally announced November 2021.

Comments: 36 pages, 21 figures, 2 tables

Science goals and new mission concepts for future exploration of Titan's atmosphere geology and habitability: Titan POlar Scout/orbitEr and In situ lake lander and DrONe explorer (POSEIDON)

Authors: Sébastien Rodriguez, Sandrine Vinatier, Daniel Cordier, Gabriel Tobie, Richard K. Achterberg, Carrie M. Anderson, Sarah V. Badman, Jason W. Barnes, Erika L. Barth, Bruno Bézard, Nathalie Carrasco, Benjamin Charnay, Roger N. Clark, Patrice Coll, Thomas Cornet, Athena Coustenis, Isabelle Couturier-Tamburelli, Michel Dobrijevic, F. Michael Flasar, Remco de Kok, Caroline Freissinet, Marina Galand, Thomas Gautier, Wolf D. Geppert, Caitlin A. Griffith , et al. (39 additional authors not shown)

Abstract: In response to ESA Voyage 2050 announcement of opportunity, we propose an ambitious L-class mission to explore one of the most exciting bodies in the Solar System, Saturn largest moon Titan. Titan, a "world with two oceans", is an organic-rich body with interior-surface-atmosphere interactions that are comparable in complexity to the Earth. Titan is also one of the few places in the Solar System w… ▽ More In response to ESA Voyage 2050 announcement of opportunity, we propose an ambitious L-class mission to explore one of the most exciting bodies in the Solar System, Saturn largest moon Titan. Titan, a "world with two oceans", is an organic-rich body with interior-surface-atmosphere interactions that are comparable in complexity to the Earth. Titan is also one of the few places in the Solar System with habitability potential. Titan remarkable nature was only partly revealed by the Cassini-Huygens mission and still holds mysteries requiring a complete exploration using a variety of vehicles and instruments. The proposed mission concept POSEIDON (Titan POlar Scout/orbitEr and In situ lake lander DrONe explorer) would perform joint orbital and in situ investigations of Titan. It is designed to build on and exceed the scope and scientific/technological accomplishments of Cassini-Huygens, exploring Titan in ways that were not previously possible, in particular through full close-up and in situ coverage over long periods of time. In the proposed mission architecture, POSEIDON consists of two major elements: a spacecraft with a large set of instruments that would orbit Titan, preferably in a low-eccentricity polar orbit, and a suite of in situ investigation components, i.e. a lake lander, a "heavy" drone (possibly amphibious) and/or a fleet of mini-drones, dedicated to the exploration of the polar regions. The ideal arrival time at Titan would be slightly before the next northern Spring equinox (2039), as equinoxes are the most active periods to monitor still largely unknown atmospheric and surface seasonal changes. The exploration of Titan northern latitudes with an orbiter and in situ element(s) would be highly complementary with the upcoming NASA New Frontiers Dragonfly mission that will provide in situ exploration of Titan equatorial regions in the mid-2030s. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: arXiv admin note: substantial text overlap with arXiv:1908.01374

Galaxy Velocity Bias in Cosmological Simulations: Towards Percent-level Calibration

Authors: Dhayaa Anbajagane, Han Aung, August E. Evrard, Arya Farahi, Daisuke Nagai, David J. Barnes, Weiguang Cui, Klaus Dolag, Ian G. McCarthy, Elena Rasia, Gustavo Yepes

Abstract: Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the rat… ▽ More Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the ratio of galaxy and DM velocity dispersions -- as a function of redshift, host halo mass, and galaxy stellar mass threshold (Mstarsat), for massive halos (M200c > 1e13.5 msun) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of b_v increasing with M200c and decreasing with redshift and Mstarsat. The ensemble-estimated theoretical uncertainty in b_v is 2-3% but becomes percent-level when considering only the three highest resolution simulations. We update the mass-richness normalization previously estimated by Farahi et al. (2016) for an SDSS redMaPPer cluster sample. The improved accuracy of our b_v estimates reduces the mass normalization uncertainty from 22% to 8%, demonstrating that dynamical estimation techniques can be competitive with weak lensing in calibrating population mean masses. We discuss necessary steps for further improving this precision. Our estimates for b_v(M200c, Mstarsat, z) are made publicly available. △ Less

Submitted 21 December, 2021; v1 submitted 4 October, 2021; originally announced October 2021.

Comments: [v1]: 8 figures, 14 pages in main text + 5 figures, 3 pages in Appendix. [v2]: Minor text edits. Accepted in MNRAS

Exploring tidal obliquity variations with SMERCURY-T

Authors: Steven M. Kreyche, Jason W. Barnes, Billy L. Quarles, John E. Chambers

Abstract: We introduce our new code, SMERCURY-T, which is based on existing codes SMERCURY (Lissauer et al. 2012) and Mercury-T (Bolmont et al. 2015). The result is a mixed-variable symplectic N-body integrator that can compute the orbital and spin evolution of a planet within a multi-planet system under the influence of tidal spin torques from its star. We validate our implementation by comparing our exper… ▽ More We introduce our new code, SMERCURY-T, which is based on existing codes SMERCURY (Lissauer et al. 2012) and Mercury-T (Bolmont et al. 2015). The result is a mixed-variable symplectic N-body integrator that can compute the orbital and spin evolution of a planet within a multi-planet system under the influence of tidal spin torques from its star. We validate our implementation by comparing our experimental results to that of a secular model. As we demonstrate in a series of experiments, SMERCURY-T allows for the study of secular spin-orbit resonance crossings and captures for planets within complex multi-planet systems. These processes can drive a planet's spin state to evolve along vastly different pathways on its road toward tidal equilibrium, as tidal spin torques dampen the planet's spin rate and evolve its obliquity. Additionally, we show the results of a scenario that exemplifies the crossing of a chaotic region that exists as the overlap of two spin-orbit resonances. The test planet experiences violent and chaotic swings in its obliquity until its eventual escape from resonance as it tidally evolves. All of these processes are and have been important over the obliquity evolution of many bodies within the Solar System and beyond, and have implications for planetary climate and habitability. SMERCURY-T is a powerful and versatile tool that allows for further study of these phenomena. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 21 pages, 9 figures, 3 tables; accepted by The Planetary Science Journal

Census of High- and Medium-mass Protostars V. CO Abundance and the Galactic $X_{\text{CO}}$ Factor

Authors: Rebecca L. Pitts, Peter J. Barnes

Abstract: We present the second dust continuum data release in the Census of High- and Medium-mass Protostars (CHaMP), expanding the methodology trialed in Pitts et al. 2019 to the entire CHaMP survey area ($280^{\circ}<l<300^{\circ}$, $-4^{\circ}<b<+2^{\circ}$). This release includes maps of dust temperature ($T_d$), H$_2$ column density ($N_{H_2}$), gas-phase CO abundance, and temperature-density plots fo… ▽ More We present the second dust continuum data release in the Census of High- and Medium-mass Protostars (CHaMP), expanding the methodology trialed in Pitts et al. 2019 to the entire CHaMP survey area ($280^{\circ}<l<300^{\circ}$, $-4^{\circ}<b<+2^{\circ}$). This release includes maps of dust temperature ($T_d$), H$_2$ column density ($N_{H_2}$), gas-phase CO abundance, and temperature-density plots for every prestellar clump with Herschel coverage, showing no evidence of internal heating for most clumps in our sample. We show that CO abundance is a strong function of $T_d$, and can be fit with a second-order polynomial in log-space, with a typical dispersion of a factor of 2--3. The CO abundance peaks at $20.0^{+0.4}_{-1.0}$ K with a value of $7.4^{+0.2}_{-0.3}\times10^{-5}$ per H$_2$; the low $T_d$ at which this maximal abundance occurs relative to laboratory results is likely due to interstellar UV bombardment in the largest survey fields. Finally, we show that, as predicted by theoretical literature and hinted at in previous studies of individual clouds, the conversion factor from integrated $^{12}$CO line intensity ($I_{^{12}CO}$) to $N_{H_2}$, the $X_{CO}$-factor, varies as a broken power-law in $I_{^{12}CO}$ with a transition zone between 70 and 90 K km$^{-1}$. The $X_{CO}$-function we propose has $N_{H_2}\propto I_{^{12}CO}^{0.51}$ for $I_{^{12}CO}\lesssim70$ K km$^{-1}$ and $N_{H_2}\propto I_{^{12}CO}^{2.3}$ for $I_{^{12}CO}\gtrsim90$ K km$^{-1}$. The high-$I_{^{12}CO}$ side should be generalizable with known adjustments for metallicity, but the influence of interstellar UV fields on the low-$I_{^{12}CO}$ side may be sample specific. We discuss how these results expand upon previous works in the CHaMP series, and help tie together observational, theoretical, and laboratory studies on CO over the past decade. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: 33 pages, 13 figures, 2 tables, 51 pages of appendices; accepted to ApJS (in print)

Proto-magnetar jets as central engines for broad-lined type Ic supernovae

Authors: Swapnil Shankar, Philipp Mösta, Jennifer Barnes, Paul C. Duffell, Daniel Kasen

Abstract: A subset of type Ic supernovae (SNe Ic), broad-lined SNe Ic (SNe Ic-bl), show unusually high kinetic energies ($\sim 10^{52}$ erg) which cannot be explained by the energy supplied by neutrinos alone. Many SNe Ic-bl have been observed in coincidence with long gamma-ray bursts (GRBs) which suggests a connection between SNe and GRBs. A small fraction of core-collapse supernovae (CCSNe) form a rapidly… ▽ More A subset of type Ic supernovae (SNe Ic), broad-lined SNe Ic (SNe Ic-bl), show unusually high kinetic energies ($\sim 10^{52}$ erg) which cannot be explained by the energy supplied by neutrinos alone. Many SNe Ic-bl have been observed in coincidence with long gamma-ray bursts (GRBs) which suggests a connection between SNe and GRBs. A small fraction of core-collapse supernovae (CCSNe) form a rapidly-rotating and strongly-magnetized protoneutron star (PNS), a proto-magnetar. Jets from such magnetars can provide the high kinetic energies observed in SNe Ic-bl and also provide the connection to GRBs. In this work we use the jetted outflow produced in a 3D CCSN simulation from a consistently formed proto-magnetar as the central engine for full-star explosion simulations. We extract a range of central engine parameters and find that the extracted engine energy is in the range of $6.231 \times 10^{51}-1.725 \times 10^{52}$ erg, the engine time-scale in the range of $0.479-1.159$ s and the engine half-opening angle in the range of $\sim 9-19^{\circ}$. Using these as central engines, we perform 2D special-relativistic (SR) hydrodynamic (HD) and radiation transfer simulations to calculate the corresponding light curves and spectra. We find that these central engine parameters successfully produce SNe Ic-bl which demonstrates that jets from proto-magnetars can be viable engines for SNe Ic-bl. We also find that only the central engines with smaller opening angles ($\sim 10^{\circ}$) form a GRB implying that GRB formation is likely associated with narrower jet outflows and Ic-bl's without GRBs may be associated with wider outflows. △ Less

Submitted 17 May, 2021; originally announced May 2021.

Comments: 13 pages, 12 figures

Ariel: Enabling planetary science across light-years

Authors: Giovanna Tinetti, Paul Eccleston, Carole Haswell, Pierre-Olivier Lagage, Jérémy Leconte, Theresa Lüftinger, Giusi Micela, Michel Min, Göran Pilbratt, Ludovic Puig, Mark Swain, Leonardo Testi, Diego Turrini, Bart Vandenbussche, Maria Rosa Zapatero Osorio, Anna Aret, Jean-Philippe Beaulieu, Lars Buchhave, Martin Ferus, Matt Griffin, Manuel Guedel, Paul Hartogh, Pedro Machado, Giuseppe Malaguti, Enric Pallé , et al. (293 additional authors not shown)

Abstract: Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths.… ▽ More Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution. △ Less

Submitted 10 April, 2021; originally announced April 2021.

Comments: Ariel Definition Study Report, 147 pages. Reviewed by ESA Science Advisory Structure in November 2020. Original document available at: https://www.cosmos.esa.int/documents/1783156/3267291/Ariel_RedBook_Nov2020.pdf/

Report number: ESA/SCI(2020)1

Titan: Earth-like on the Outside, Ocean World on the Inside

Authors: Shannon M. MacKenzie, Samuel P. D. Birch, Sarah Horst, Christophe Sotin, Erika Barth, Juan M. Lora, Melissa G. Trainer, Paul Corlies, Michael J. Malaska, Ella Sciamma-O'Brien, Alexander E. Thelen, Elizabeth P. Turtle, Jani Radebaugh, Jennifer Hanley, Anezina Solomonidou, Claire Newman, Leonardo Regoli, Sebastien Rodriguez, Benoit Seignovert, Alexander G. Hayes, Baptiste Journaux, Jordan Steckloff, Delphine Nna-Mvondo, Thomas Cornet, Maureen Palmer , et al. (8 additional authors not shown)

Abstract: Thanks to the Cassini-Huygens mission, Titan, the pale orange dot of Pioneer and Voyager encounters has been revealed to be a dynamic, hydrologically-shaped, organic-rich ocean world offering unparalleled opportunities to explore prebiotic chemistry. And while Cassini-Huygens revolutionized our understanding of each of the three layers of Titan--the atmosphere, the surface, and the interior--we ar… ▽ More Thanks to the Cassini-Huygens mission, Titan, the pale orange dot of Pioneer and Voyager encounters has been revealed to be a dynamic, hydrologically-shaped, organic-rich ocean world offering unparalleled opportunities to explore prebiotic chemistry. And while Cassini-Huygens revolutionized our understanding of each of the three layers of Titan--the atmosphere, the surface, and the interior--we are only beginning to hypothesize how these realms interact. In this paper, we summarize the current state of Titan knowledge and discuss how future exploration of Titan would address some of the next decade's most compelling planetary science questions. We also demonstrate why exploring Titan, both with and beyond the Dragonfly New Frontiers mission, is a necessary and complementary component of an Ocean Worlds Program that seeks to understand whether habitable environments exist elsewhere in our solar system. △ Less

Submitted 16 February, 2021; originally announced February 2021.

Comments: Submitted to the PSJ Focus Issue on Ocean World Exploration

Solving the Alhazen-Ptolemy Problem: Determining Specular Points on Spherical Surfaces for Radiative Transfer of Titan's Seas

Authors: William J. Miller, Jason W. Barnes, Shannon M. MacKenzie

Abstract: Given a light source, a spherical reflector, and an observer, where on the surface of the sphere will the light be directly reflected to the observer, i.e. where is the the specular point? This is known as the Alhazen-Ptolemy problem, and finding this specular point for spherical reflectors is useful in applications ranging from computer rendering to atmospheric modeling to GPS communications. Exi… ▽ More Given a light source, a spherical reflector, and an observer, where on the surface of the sphere will the light be directly reflected to the observer, i.e. where is the the specular point? This is known as the Alhazen-Ptolemy problem, and finding this specular point for spherical reflectors is useful in applications ranging from computer rendering to atmospheric modeling to GPS communications. Existing solutions rely upon finding the roots of a quartic equation and evaluating numerically which root provides the real specular point. We offer a formulation, and two solutions thereof, for which the correct root is predeterminable, thereby allowing the construction of the fully analytical solutions we present. Being faster to compute, our solutions should prove useful in cases which require repeated calculation of the specular point, such as Monte-Carlo radiative transfer, including reflections off of Titan's hydrocarbon seas. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Comments: 11 pages with 5 figures. Submitted to Planetary Science Journal

The SEDIGISM survey: Molecular clouds in the inner Galaxy

Authors: A. Duarte-Cabral, D. Colombo, J. S. Urquhart, A. Ginsburg, D. Russeil, F. Schuller, L. D. Anderson, P. J. Barnes, M. T. Beltran, H. Beuther, S. Bontemps, L. Bronfman, T. Csengeri, C. L. Dobbs, D. Eden, A. Giannetti, J. Kauffmann, M. Mattern, S. -N. X. Medina, K. M. Menten, M. -Y. Lee, A. R. Pettitt, M. Riener, A. J. Rigby, A. Trafficante , et al. (35 additional authors not shown)

Abstract: We use the 13CO(2-1) emission from the SEDIGISM high-resolution spectral-line survey of the inner Galaxy, to extract the molecular cloud population with a large dynamic range in spatial scales, using the SCIMES algorithm. This work compiles a cloud catalogue with a total of 10663 molecular clouds, 10300 of which we were able to assign distances and compute physical properties. We study some of the… ▽ More We use the 13CO(2-1) emission from the SEDIGISM high-resolution spectral-line survey of the inner Galaxy, to extract the molecular cloud population with a large dynamic range in spatial scales, using the SCIMES algorithm. This work compiles a cloud catalogue with a total of 10663 molecular clouds, 10300 of which we were able to assign distances and compute physical properties. We study some of the global properties of clouds using a science sample, consisting of 6664 well resolved sources and for which the distance estimates are reliable. In particular, we compare the scaling relations retrieved from SEDIGISM to those of other surveys, and we explore the properties of clouds with and without high-mass star formation. Our results suggest that there is no single global property of a cloud that determines its ability to form massive stars, although we find combined trends of increasing mass, size, surface density and velocity dispersion for the sub-sample of clouds with ongoing high-mass star formation. We then isolate the most extreme clouds in the SEDIGISM sample (i.e. clouds in the tails of the distributions) to look at their overall Galactic distribution, in search for hints of environmental effects. We find that, for most properties, the Galactic distribution of the most extreme clouds is only marginally different to that of the global cloud population. The Galactic distribution of the largest clouds, the turbulent clouds and the high-mass star-forming clouds are those that deviate most significantly from the global cloud population. We also find that the least dynamically active clouds (with low velocity dispersion or low virial parameter) are situated further afield, mostly in the least populated areas. However, we suspect that part of these trends may be affected by some observational biases, and thus require further follow up work in order to be confirmed. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: 25 pages (+ appendices, 15 pages), 26 figures, MNRAS

SEDIGISM-ATLASGAL: Dense Gas Fraction and Star Formation Efficiency Across the Galactic Disk

Authors: J. S. Urquhart, C. Figura, J. R. Cross, M. R. A. Wells, T. J. T. Moore, D. J. Eden, S. E. Ragan, A. R. Pettitt, A. Duarte-Cabral, D. Colombo, F. Schuller, T. Csengeri, M. Mattern, H. Beuther, K. M. Menten, F. Wyrowski, L. D. Anderson, P. J. Barnes, M. T. Beltrán, S. J. Billington, L. Bronfman, A. Giannetti, J. Kainulainen, J. Kauffmann, M. -Y. Lee , et al. (10 additional authors not shown)

Abstract: By combining two surveys covering a large fraction of the molecular material in the Galactic disk we investigate the role the spiral arms play in the star formation process. We have matched clumps identified by ATLASGAL with their parental GMCs as identified by SEDIGISM, and use these giant molecular cloud (GMC) masses, the bolometric luminosities, and integrated clump masses obtained in a concurr… ▽ More By combining two surveys covering a large fraction of the molecular material in the Galactic disk we investigate the role the spiral arms play in the star formation process. We have matched clumps identified by ATLASGAL with their parental GMCs as identified by SEDIGISM, and use these giant molecular cloud (GMC) masses, the bolometric luminosities, and integrated clump masses obtained in a concurrent paper to estimate the dense gas fractions (DGF$_{\rm gmc}=\sum M_{\rm clump}/M_{\rm gmc}$) and the instantaneous star forming efficiencies (i.e., SFE$_{\rm gmc} = \sum L_{\rm clump}/M_{\rm gmc}$). We find that the molecular material associated with ATLASGAL clumps is concentrated in the spiral arms ($\sim$60% found within $\pm$10 km s$^{-1}$ of an arm). We have searched for variations in the values of these physical parameters with respect to their proximity to the spiral arms, but find no evidence for any enhancement that might be attributable to the spiral arms. The combined results from a number of similar studies based on different surveys indicate that, while spiral-arm location plays a role in cloud formation and HI to H$_2$ conversion, the subsequent star formation processes appear to depend more on local environment effects. This leads us to conclude that the enhanced star formation activity seen towards the spiral arms is the result of source crowding rather than the consequence of a any physical process. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: 15 pages, 12 figures. Full list of affiliations can be found at the end of the paper

The splashback boundary of haloes in hydrodynamic simulations

Authors: Stephanie O'Neil, David J. Barnes, Mark Vogelsberger, Benedikt Diemer

Abstract: The splashback radius, $R_{\rm sp}$, is a physically motivated halo boundary that separates infalling and collapsed matter of haloes. We study $R_{\rm sp}$ in the hydrodynamic and dark matter only IllustrisTNG simulations. The most commonly adopted signature of $R_{\rm sp}$ is the radius at which the radial density profiles are steepest. Therefore, we explicitly optimise our density profile fit to… ▽ More The splashback radius, $R_{\rm sp}$, is a physically motivated halo boundary that separates infalling and collapsed matter of haloes. We study $R_{\rm sp}$ in the hydrodynamic and dark matter only IllustrisTNG simulations. The most commonly adopted signature of $R_{\rm sp}$ is the radius at which the radial density profiles are steepest. Therefore, we explicitly optimise our density profile fit to the profile slope and find that this leads to a $\sim5\%$ larger radius compared to other optimisations. We calculate $R_{\rm sp}$ for haloes with masses between $10^{13-15}{\rm M}_{\odot}$ as a function of halo mass, accretion rate and redshift. $R_{\rm sp}$ decreases with mass and with redshift for haloes of similar $M_{\rm200m}$ in agreement with previous work. We also find that $R_{\rm sp}/R_{\rm200m}$ decreases with halo accretion rate. We apply our analysis to dark matter, gas and satellite galaxies associated with haloes to investigate the observational potential of $R_{\rm sp}$. The radius of steepest slope in gas profiles is consistently smaller than the value calculated from dark matter profiles. The steepest slope in galaxy profiles, which are often used in observations, tends to agree with dark matter profiles but is lower for less massive haloes. We compare $R_{\rm sp}$ in hydrodynamic and N-body dark matter only simulations and do not find a significant difference caused by the addition of baryonic physics. Thus, results from dark matter only simulations should be applicable to realistic haloes. △ Less

Submitted 25 May, 2021; v1 submitted 30 November, 2020; originally announced December 2020.

Comments: 18 pages, 15 figures, accepted by MNRAS

SHAPing the Gas: Understanding Gas Shapes in Dark Matter Haloes with Interpretable Machine Learning

Authors: Luis Fernando Machado Poletti Valle, Camille Avestruz, David J. Barnes, Arya Farahi, Erwin T. Lau, Daisuke Nagai

Abstract: The non-spherical shapes of dark matter and gas distributions introduce systematic uncertainties that affect observable-mass relations and selection functions of galaxy groups and clusters. However, the triaxial gas distributions depend on the non-linear physical processes of halo formation histories and baryonic physics, which are challenging to model accurately. In this study we explore a machin… ▽ More The non-spherical shapes of dark matter and gas distributions introduce systematic uncertainties that affect observable-mass relations and selection functions of galaxy groups and clusters. However, the triaxial gas distributions depend on the non-linear physical processes of halo formation histories and baryonic physics, which are challenging to model accurately. In this study we explore a machine learning approach for modelling the dependence of gas shapes on dark matter and baryonic properties. With data from the IllustrisTNG hydrodynamical cosmological simulations, we develop a machine learning pipeline that applies \pkg{XGBoost}, an implementation of gradient boosted decision trees, to predict radial profiles of gas shapes from halo properties. We show that \pkg{XGBoost} models can accurately predict gas shape profiles in dark matter haloes. We also explore model interpretability with \pkg{SHAP}, a method that identifies the most predictive properties at different halo radii. We find that baryonic properties best predict gas shapes in halo cores, whereas dark matter shapes are the main predictors in the halo outskirts. This work demonstrates the power of interpretable machine learning in modelling observable properties of dark matter haloes in the era of multi-wavelength cosmological surveys. △ Less

Submitted 2 August, 2021; v1 submitted 25 November, 2020; originally announced November 2020.

Comments: 14 pages, 8 figures, Accepted to MNRAS

The Cluster HEritage project with XMM-Newton: Mass Assembly and Thermodynamics at the Endpoint of structure formation. I. Programme overview

Authors: The CHEX-MATE Collaboration, :, M. Arnaud, S. Ettori, G. W. Pratt, M. Rossetti, D. Eckert, F. Gastaldello, R. Gavazzi, S. T. Kay, L. Lovisari, B. J. Maughan, E. Pointecouteau, M. Sereno, I. Bartalucci, A. Bonafede, H. Bourdin, R. Cassano, R. T. Duffy, A. Iqbal, S. Maurogordato, E. Rasia, J. Sayers, F. Andrade-Santos, H. Aussel , et al. (45 additional authors not shown)

Abstract: The Cluster HEritage project with XMM-Newton - Mass Assembly and Thermodynamics at the Endpoint of structure formation (CHEX-MATE) is a three mega-second Multi-Year Heritage Programme to obtain X-ray observations of a minimally-biased, signal-to-noise limited sample of 118 galaxy clusters detected by Planck through the Sunyaev-Zeldovich effect. The programme, described in detail in this paper, aim… ▽ More The Cluster HEritage project with XMM-Newton - Mass Assembly and Thermodynamics at the Endpoint of structure formation (CHEX-MATE) is a three mega-second Multi-Year Heritage Programme to obtain X-ray observations of a minimally-biased, signal-to-noise limited sample of 118 galaxy clusters detected by Planck through the Sunyaev-Zeldovich effect. The programme, described in detail in this paper, aims to study the ultimate products of structure formation in time and mass. It is composed of a census of the most recent objects to have formed (Tier-1: 0.05 < z < 0.2; 2 x 10e14 M_sun < M_500 < 9 x 10e14 M_sun), together with a sample of the highest-mass objects in the Universe (Tier-2: z < 0.6; M_500 > 7.25 x 10e14 M_sun). The programme will yield an accurate vision of the statistical properties of the underlying population, measure how the gas properties are shaped by collapse into the dark matter halo, uncover the provenance of non-gravitational heating, and resolve the major uncertainties in mass determination that limit the use of clusters for cosmological parameter estimation. We will acquire X-ray exposures of uniform depth, designed to obtain individual mass measurements accurate to 15-20% under the hydrostatic assumption. We present the project motivations, describe the programme definition, and detail the ongoing multi-wavelength observational (lensing, SZ, radio) and theoretical effort that is being deployed in support of the project. △ Less

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

Comments: 27 pages, 11 figures; A&A, in press

Journal ref: A&A 650, A104 (2021)

Kilonovae across the nuclear physics landscape: The impact of nuclear physics uncertainties on r-process-powered emission

Authors: Jennifer Barnes, Y. L. Zhu, K. A. Lund, T. M. Sprouse, N. Vassh, G. C. McLaughlin, M. R. Mumpower, R. Surman

Abstract: Merging neutron stars produce "kilonovae"---electromagnetic transients powered by the decay of unstable nuclei synthesized via rapid neutron capture (the r-process) in material that is gravitationally unbound during inspiral and coalescence. Kilonova emission, if accurately interpreted, can be used to characterize the masses and compositions of merger-driven outflows, helping to resolve a long-sta… ▽ More Merging neutron stars produce "kilonovae"---electromagnetic transients powered by the decay of unstable nuclei synthesized via rapid neutron capture (the r-process) in material that is gravitationally unbound during inspiral and coalescence. Kilonova emission, if accurately interpreted, can be used to characterize the masses and compositions of merger-driven outflows, helping to resolve a long-standing debate about the origins of r-process material in the Universe. We explore how the uncertain properties of nuclei involved in the r-process complicate the inference of outflow properties from kilonova observations. Using r-process simulations, we show how nuclear physics uncertainties impact predictions of radioactive heating and element synthesis. For a set of models that span a large range in both predicted heating and final abundances, we carry out detailed numerical calculations of decay product thermalization and radiation transport in a kilonova ejecta with a fixed mass and density profile. The light curves associated with our models exhibit great diversity in their luminosities, with peak brightness varying by more than an order of magnitude. We also find variability in the shape of the kilonova light curves and their color, which in some cases runs counter to the expectation that increasing levels of lanthanide and/or actinide enrichment will be correlated with longer, dimmer, redder emission. △ Less

Submitted 21 October, 2020; originally announced October 2020.

Comments: Submitted to ApJ; comments welcome

Modeling Kilonova Light Curves: Dependence on Nuclear Inputs

Authors: Y. L. Zhu, K. Lund, J. Barnes, T. M. Sprouse, N. Vassh, G. C. McLaughlin, M. R. Mumpower, R. Surman

Abstract: The mergers of binary neutron stars, as well as black hole-neutron star systems, are expected to produce an electromagnetic counterpart that can be analyzed to infer the element synthesis that occurred in these events. We investigate one source of uncertainties pertinent to lanthanide-rich outflows: the nuclear inputs to rapid neutron capture nucleosynthesis calculations. We begin by examining thi… ▽ More The mergers of binary neutron stars, as well as black hole-neutron star systems, are expected to produce an electromagnetic counterpart that can be analyzed to infer the element synthesis that occurred in these events. We investigate one source of uncertainties pertinent to lanthanide-rich outflows: the nuclear inputs to rapid neutron capture nucleosynthesis calculations. We begin by examining thirty-two different combinations of nuclear inputs: eight mass models, two types of spontaneous fission rates, and two types of fission daughter product distributions. We find that such nuclear physics uncertainties typically generate at least one order of magnitude uncertainty in key quantities such as the nuclear heating (one and a half orders of magnitude at one day post-merger), the bolometric luminosity (one order of magnitude at five days post-merger), and the inferred mass of material from the bolometric luminosity (factor of eight when considering the eight to ten days region). Since particular nuclear processes are critical for determining the electromagnetic signal, we provide tables of key nuclei undergoing $β$-decay, $α$-decay, and spontaneous fission important for heating at different times, identifying decays that are common among the many nuclear input combinations. △ Less

Submitted 7 October, 2020; originally announced October 2020.

Comments: 29 pages, 17 figures, 4 tables

Report number: LA-UR-20-22630

A Search for In-Situ Field OB Star Formation in the Small Magellanic Cloud

Authors: Irene Vargas-Salazar, M. S. Oey, Jesse R. Barnes, Xinyi Chen, N. Castro, Kaitlin M. Kratter, Timothy A. Faerber

Abstract: Whether any OB stars form in isolation is a question central to theories of massive star formation. To address this, we search for tiny, sparse clusters around 210 field OB stars from the Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC (RIOTS4), using friends-of-friends (FOF) and nearest neighbors (NN) algorithms. We also stack the target fields to evaluate the presence of an agg… ▽ More Whether any OB stars form in isolation is a question central to theories of massive star formation. To address this, we search for tiny, sparse clusters around 210 field OB stars from the Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC (RIOTS4), using friends-of-friends (FOF) and nearest neighbors (NN) algorithms. We also stack the target fields to evaluate the presence of an aggregate density enhancement. Using several statistical tests, we compare these observations with three random-field datasets, and we also compare the known runaways to non-runaways. We find that the local environments of non-runaways show higher aggregate central densities than for runaways, implying the presence of some "tips-of-iceberg" (TIB) clusters. We find that the frequency of these tiny clusters is low, $\sim 4-5\%$ of our sample. This fraction is much lower than some previous estimates, but is consistent with field OB stars being almost entirely runaway and walkaway stars. The lack of TIB clusters implies that such objects either evaporate on short timescales, or do not form, implying a higher cluster lower-mass limit and consistent with a relationship between maximum stellar mass ($m_{\rm max}$) and the mass of the cluster ($M_{\rm cl}$). On the other hand, we also cannot rule out that some OB stars may form in highly isolated conditions. Our results set strong constraints on the formation of massive stars in relative isolation. △ Less

Submitted 25 September, 2020; originally announced September 2020.

Comments: 24 pages, 13 figures, Accepted to ApJ

The Broad-band Counterpart of the Short GRB 200522A at $z=0.5536$: A Luminous Kilonova or a Collimated Outflow with a Reverse Shock?

Authors: W. Fong, T. Laskar, J. Rastinejad, A. Rouco Escorial, G. Schroeder, J. Barnes, C. D. Kilpatrick, K. Paterson, E. Berger, B. D. Metzger, Y. Dong, A. E. Nugent, R. Strausbaugh, P. K. Blanchard, A. Goyal, A. Cucchiara, G. Terreran, K. D. Alexander, T. Eftekhari, C. Fryer, B. Margalit, R. Margutti, M. Nicholl

Abstract: We present the discovery of the radio afterglow and near-infrared (NIR) counterpart of the Swift short GRB 200522A, located at a small projected offset of $\approx 1$ kpc from the center of a young, star-forming host galaxy at $z=0.5536$. The radio and X-ray luminosities of the afterglow are consistent with those of on-axis cosmological short GRBs. The NIR counterpart, revealed by our HST observat… ▽ More We present the discovery of the radio afterglow and near-infrared (NIR) counterpart of the Swift short GRB 200522A, located at a small projected offset of $\approx 1$ kpc from the center of a young, star-forming host galaxy at $z=0.5536$. The radio and X-ray luminosities of the afterglow are consistent with those of on-axis cosmological short GRBs. The NIR counterpart, revealed by our HST observations at a rest-frame time of $\approx2.3$ days, has a luminosity of $\approx (1.3-1.7) \times 10^{42}$ erg s$^{-1}$. This is substantially lower than on-axis short GRB afterglow detections, but is a factor of $\approx 8$-$17$ more luminous than the kilonova of GW170817, and significantly more luminous than any kilonova candidate for which comparable observations exist. The combination of the counterpart's color ($i-y = -0.08\pm 0.21$; rest-frame) and luminosity cannot be explained by standard radioactive heating alone. We present two scenarios to interpret the broad-band behavior of GRB 200522A: a synchrotron forward shock with a luminous kilonova (potentially boosted by magnetar energy deposition), or forward and reverse shocks from a $\approx14^{\circ}$, relativistic ($Γ_0 \gtrsim 80$) jet. Models which include a combination of enhanced radioactive heating rates, low-lanthanide mass fractions, or additional sources of heating from late-time central engine activity may provide viable alternate explanations. If a stable magnetar was indeed produced in GRB 200522A, we predict that late-time radio emission will be detectable starting $\approx 0.3$-$6$ years after the burst for a deposited energy of $\approx 10^{53}$ erg. Counterparts of similar luminosity to GRB 200522A associated with gravitational wave events will be detectable with current optical searches to $\approx\!250$ Mpc. △ Less

Submitted 19 August, 2020; originally announced August 2020.

Comments: 33 pages, 13 figures, 5 tables. Submitted to AAS Journals

First Light And Reionisation Epoch Simulations (FLARES) II: The Photometric Properties of High-Redshift Galaxies

Authors: Aswin P. Vijayan, Christopher C. Lovell, Stephen M. Wilkins, Peter A. Thomas, David J. Barnes, Dimitrios Irodotou, Jussi Kuusisto, Will Roper

Abstract: We present the photometric properties of galaxies in the First Light and Reionisation Epoch Simulations (FLARES). The simulations trace the evolution of galaxies in a range of overdensities through the Epoch of Reionistion (EoR). With a novel weighting scheme we combine these overdensities, extending significantly the dynamic range of observed composite distribution functions compared to periodic… ▽ More We present the photometric properties of galaxies in the First Light and Reionisation Epoch Simulations (FLARES). The simulations trace the evolution of galaxies in a range of overdensities through the Epoch of Reionistion (EoR). With a novel weighting scheme we combine these overdensities, extending significantly the dynamic range of observed composite distribution functions compared to periodic simulation boxes. FLARES predicts a significantly larger number of intrinsically bright galaxies, which can be explained through a simple model linking dust-attenuation to the metal content of the interstellar medium, using a line-of-sight (LOS) extinction model. With this model we present the photometric properties of the FLARES galaxies for $z \in [5,10]$. We show that the ultraviolet (UV) luminosity function (LF) matches the observations at all redshifts. The function is fit by Schechter and double power-law forms, with the latter being favoured at these redshifts by the FLARES composite UV LF. We also present predictions for the UV continuum slope as well as the attenuation in the UV. The impact of environment on the UV LF is also explored, with the brightest galaxies forming in the densest environments. We then present the line luminosity and equivalent widths of some prominent nebular emission lines arising from the galaxies, finding rough agreement with available observations. We also look at the relative contribution of obscured and unobscured star formation, finding comparable contributions at these redshifts. △ Less

Submitted 14 December, 2020; v1 submitted 13 August, 2020; originally announced August 2020.

Comments: 21 pages, 22 figures. Data available at https://flaresimulations.github.io/data.html

The Science Case for a Titan Flagship-class Orbiter with Probes

Authors: Conor A. Nixon, James Abshire, Andrew Ashton, Jason W. Barnes, Nathalie Carrasco, Mathieu Choukroun, Athena Coustenis, Louis-Alexandre Couston, Niklas Edberg, Alexander Gagnon, Jason D. Hofgartner, Luciano Iess, Stéphane Le Mouélic, Rosaly Lopes, Juan Lora, Ralph D. Lorenz, Adrienn Luspay-Kuti, Michael Malaska, Kathleen Mandt, Marco Mastrogiuseppe, Erwan Mazarico, Marc Neveu, Taylor Perron, Jani Radebaugh, Sébastien Rodriguez , et al. (14 additional authors not shown)

Abstract: We outline a flagship-class mission concept focused on studying Titan as a global system, with particular emphasis on the polar regions. Investigating Titan from the unique standpoint of a polar orbit would enable comprehensive global maps to uncover the physics and chemistry of the atmosphere, and the topography and geophysical environment of the surface and subsurface. The mission includes two k… ▽ More We outline a flagship-class mission concept focused on studying Titan as a global system, with particular emphasis on the polar regions. Investigating Titan from the unique standpoint of a polar orbit would enable comprehensive global maps to uncover the physics and chemistry of the atmosphere, and the topography and geophysical environment of the surface and subsurface. The mission includes two key elements: (1) an orbiter spacecraft, which also acts as a data relay, and (2) one or more small probes to directly investigate Titan's seas and make the first direct measurements of their liquid composition and physical environment. The orbiter would carry a sophisticated remote sensing payload, including a novel topographic lidar, a long-wavelength surface-penetrating radar, a sub-millimeter sounder for winds and for mesospheric/thermospheric composition, and a camera and near-infrared spectrometer. An instrument suite to analyze particles and fields would include a mass spectrometer to focus on the interactions between Titan's escaping upper atmosphere and the solar wind and Saturnian magnetosphere. The orbiter would enter a stable polar orbit around 1500 to 1800 km, from which vantage point it would make global maps of the atmosphere and surface. One or more probes, released from the orbiter, would investigate Titan's seas in situ, including possible differences in composition between higher and lower latitude seas, as well as the atmosphere during the parachute descent. The number of probes, as well as the instrument complement on the orbiter and probe, remain to be finalized during a mission study that we recommend to NASA as part of the NRC Decadal Survey for Planetary Science now underway, with the goal of an overall mission cost in the "small flagship" category of ~$2 bn. International partnerships, similar to Cassini-Huygens, may also be included for consideration. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: 13 pages, white paper submitted to the NRC Decadal Survey for Planetary Science and Astrobiology

A deep learning view of the census of galaxy clusters in IllustrisTNG

Authors: Y. Su, Y. Zhang, G. Liang, J. A. ZuHone, D. J. Barnes, N. B. Jacobs, M. Ntampaka, W. R. Forman, P. E. J. Nulsen, R. P. Kraft, C. Jones

Abstract: The origin of the diverse population of galaxy clusters remains an unexplained aspect of large-scale structure formation and cluster evolution. We present a novel method of using X-ray images to identify cool core (CC), weak cool core (WCC), and non cool core (NCC) clusters of galaxies, that are defined by their central cooling times. We employ a convolutional neural network, ResNet-18, which is c… ▽ More The origin of the diverse population of galaxy clusters remains an unexplained aspect of large-scale structure formation and cluster evolution. We present a novel method of using X-ray images to identify cool core (CC), weak cool core (WCC), and non cool core (NCC) clusters of galaxies, that are defined by their central cooling times. We employ a convolutional neural network, ResNet-18, which is commonly used for image analysis, to classify clusters. We produce mock Chandra X-ray observations for a sample of 318 massive clusters drawn from the IllustrisTNG simulations. The network is trained and tested with low resolution mock Chandra images covering a central 1 Mpc square for the clusters in our sample. Without any spectral information, the deep learning algorithm is able to identify CC, WCC, and NCC clusters, achieving balanced accuracies (BAcc) of 92%, 81%, and 83%, respectively. The performance is superior to classification by conventional methods using central gas densities, with an average BAcc = 81%, or surface brightness concentrations, giving BAcc = 73%. We use Class Activation Mapping to localize discriminative regions for the classification decision. From this analysis, we observe that the network has utilized regions from cluster centers out to r~300 kpc and r~500 kpc to identify CC and NCC clusters, respectively. It may have recognized features in the intracluster medium that are associated with AGN feedback and disruptive major mergers. △ Less

Submitted 25 August, 2020; v1 submitted 9 July, 2020; originally announced July 2020.

Comments: Accepted for publication in MNRAS

Tidal Currents Detected in Kraken Mare Straits from Cassini VIMS Sun Glitter Observations

Authors: Michael F. Heslar, Jason W. Barnes, Jason M. Soderblom, Benoit Seignovert, Rajani D. Dhingra, Christophe Sotin

Abstract: We present Cassini VIMS observations of sun glitter -- wave-induced reflections from a liquid surface offset from a specular point -- on Kraken Mare. Sun glitter reveals rough sea surfaces around Kraken Mare, namely the coasts and narrow straits. The sun glitter observations indicate wave activity driven by the winds and tidal currents in Kraken Mare during northern summer. T104 Cassini VIMS obser… ▽ More We present Cassini VIMS observations of sun glitter -- wave-induced reflections from a liquid surface offset from a specular point -- on Kraken Mare. Sun glitter reveals rough sea surfaces around Kraken Mare, namely the coasts and narrow straits. The sun glitter observations indicate wave activity driven by the winds and tidal currents in Kraken Mare during northern summer. T104 Cassini VIMS observations show three sun glitter features in Bayta Fretum indicative of variegated wave fields. We cannot uniquely determine one source for the coastal Bayta waves, but we lean toward the interpretation of surface winds, because tidal currents should be too weak to generate capillary-gravity waves in Bayta Fretum. T105 and T110 observations reveal wave fields in the straits of Seldon Fretum, Lulworth Sinus, and Tunu Sinus that likely originate from the constriction of tidal currents. Coastlines of Bermoothes and Hufaidh Insulae adjoin rough sea surfaces, suggesting a complex interplay of wind-roughened seas and localized tidal currents. Bermoothes and Hufaidh Insulae may share characteristics of either the Torres Strait off Australia or the Aland region of Finland, summarized as an island-dense strait with shallow bathymetry that hosts complex surface circulation patterns. Hufaidh Insulae could host seafloor bedforms formed by tidal currents with an abundant sediment supply, similar to the Torres Strait. The coastlines of Hufaidh and Bermoothes Insulae likely host ria or flooded coastal inlets, suggesting the Insulae may be local peaks of primordial crust isolated by an episode of sea-level rise or tectonic uplift. △ Less

Submitted 1 July, 2020; originally announced July 2020.

Comments: Accepted to the Planetary Science Journal, 25 pages, 13 figures

A multiple planet system of super-Earths orbiting the brightest red dwarf star GJ887

Authors: S. V. Jeffers, S. Dreizler, J. R. Barnes, C. A. Haswell, R. P. Nelson, E. Rodríguez, M. J. Lopez-Gonzalez, N. Morales, R. Luque, M. Zechmeister, S. S. Vogt, J. S. Jenkins, E. Palle, Z. M. Berdinas, G. A. L. Coleman, M. R. Diaz, I. Ribas, H. R. A. Jones, R. P. Butler, C. G. Tinney, J. Bailey, B. D. Carter, S. ~O'Toole, R. A. Wittenmyer, J. D. Crane , et al. (7 additional authors not shown)

Abstract: The nearest exoplanets to the Sun are our best possibilities for detailed characterization. We report the discovery of a compact multi-planet system of super-Earths orbiting the nearby red dwarf GJ 887, using radial velocity measurements. The planets have orbital periods of 9.3 and 21.8~days. Assuming an Earth-like albedo, the equilibrium temperature of the 21.8 day planet is approx 350 K; which i… ▽ More The nearest exoplanets to the Sun are our best possibilities for detailed characterization. We report the discovery of a compact multi-planet system of super-Earths orbiting the nearby red dwarf GJ 887, using radial velocity measurements. The planets have orbital periods of 9.3 and 21.8~days. Assuming an Earth-like albedo, the equilibrium temperature of the 21.8 day planet is approx 350 K; which is interior, but close to the inner edge, of the liquid-water habitable zone. We also detect a further unconfirmed signal with a period of 50 days which could correspond to a third super-Earth in a more temperate orbit. GJ 887 is an unusually magnetically quiet red dwarf with a photometric variability below 500 parts-per-million, making its planets amenable to phase-resolved photometric characterization. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: Published by Science

Studying galaxy cluster morphological metrics with Mock-X

Authors: Kaili Cao, David J. Barnes, Mark Vogelsberger

Abstract: Dynamically relaxed galaxy clusters have long played a role in galaxy cluster studies because it is thought their properties can be reconstructed more precisely and with less systematics. As relaxed clusters are desirable, there exist a plethora of criteria for classifying a galaxy cluster as relaxed. In this work, we examine $9$ commonly used observational and theoretical morphological metrics ex… ▽ More Dynamically relaxed galaxy clusters have long played a role in galaxy cluster studies because it is thought their properties can be reconstructed more precisely and with less systematics. As relaxed clusters are desirable, there exist a plethora of criteria for classifying a galaxy cluster as relaxed. In this work, we examine $9$ commonly used observational and theoretical morphological metrics extracted from $54,000$ Mock-X synthetic X-ray images of galaxy clusters taken from the IllustrisTNG, BAHAMAS and MACSIS simulation suites. We find that the simulated criteria distributions are in reasonable agreement with the observed distributions. Many criteria distributions evolve as a function of redshift, cluster mass, numerical resolution and subgrid physics, limiting the effectiveness of a single relaxation threshold value. All criteria are positively correlated with each other, however, the strength of the correlation is sensitive to redshift, mass and numerical choices. Driven by the intrinsic scatter inherent to all morphological metrics and the arbitrary nature of relaxation threshold values, we find the consistency of relaxed subsets defined by the different metrics to be relatively poor. Therefore, the use of relaxed cluster subsets introduces significant selection effects that are non-trivial to resolve. △ Less

Submitted 18 June, 2020; originally announced June 2020.

Comments: 21 pages, 6 figures, submitted to MNRAS. Comments welcome

Variability in Short Gamma-ray Bursts: Gravitationally Unstable Tidal Tails

Authors: Eric R. Coughlin, C. J. Nixon, Jennifer Barnes, Brian D. Metzger, R. Margutti

Abstract: Short gamma-ray bursts are thought to result from the mergers of two neutron stars or a neutron star and stellar mass black hole. The final stages of the merger are generally accompanied by the production of one or more tidal "tails" of ejecta, which fall back onto the remnant-disc system at late times. Using the results of a linear stability analysis, we show that if the material comprising these… ▽ More Short gamma-ray bursts are thought to result from the mergers of two neutron stars or a neutron star and stellar mass black hole. The final stages of the merger are generally accompanied by the production of one or more tidal "tails" of ejecta, which fall back onto the remnant-disc system at late times. Using the results of a linear stability analysis, we show that if the material comprising these tails is modeled as adiabatic and the effective adiabatic index satisfies $γ\ge 5/3$, then the tails are gravitationally unstable and collapse to form small-scale knots. We analytically estimate the properties of these knots, including their spacing along the tidal tail and the total number produced, and their effect on the mass return rate to the merger remnant. We perform hydrodynamical simulations of the disruption of a polytropic (with the polytropic and adiabatic indices $γ$ equal), $γ=2$ neutron star by a black hole, and find agreement between the predictions of the linear stability analysis and the distribution of knots that collapse out of the instability. The return of these knots to the black hole induces variability in the fallback rate, which can manifest as variability in the lightcurve of the GRB and -- depending on how rapidly the instability operates -- the prompt emission. The late-time variability induced by the return of these knots is also consistent with the extended emission observed in some GRBs. △ Less

Submitted 23 June, 2020; v1 submitted 4 June, 2020; originally announced June 2020.

Comments: Small corrections, additional references included to reflect ApJL published version

Journal ref: ApJ, 896, L38 (2020)

Redshift evolution of the hot intracluster gas metallicity in the C-EAGLE cluster simulations

Authors: Francesca A. Pearce, Scott T. Kay, David J. Barnes, Yannick M. Bahe, Richard G. Bower

Abstract: The abundance and distribution of metals in galaxy clusters contains valuable information about their chemical history and evolution. By looking at how metallicity evolves with redshift, it is possible to constrain the different metal production channels. We use the C-EAGLE clusters, a sample of 30 high resolution ($m_{gas} \simeq 1.8\times 10^{6}$ M$_{\odot}$) cluster zoom simulations, to investi… ▽ More The abundance and distribution of metals in galaxy clusters contains valuable information about their chemical history and evolution. By looking at how metallicity evolves with redshift, it is possible to constrain the different metal production channels. We use the C-EAGLE clusters, a sample of 30 high resolution ($m_{gas} \simeq 1.8\times 10^{6}$ M$_{\odot}$) cluster zoom simulations, to investigate the redshift evolution of metallicity, with particular focus on the cluster outskirts. The early enrichment model, in which the majority of metals are produced in the core of cluster progenitors at high redshift, suggests that metals in cluster outskirts have not significantly evolved since $z=2$. With the C-EAGLE sample, we find reasonable agreement with the early enrichment model as there is very little scatter in the metallicity abundance at large radius across the whole sample, out to at least $z=2$. The exception is Fe for which the radial dependence of metallicity was found to evolve at low redshift as a result of being mainly produced by Type Ia supernovae, which are more likely to be formed at later times than core-collapse supernovae. We also found considerable redshift evolution of metal abundances in the cores of the C-EAGLE clusters which has not been seen in other simulations or observation based metallicity studies. Since we find this evolution to be driven by accretion of low metallicity gas, it suggests that the interaction between outflowing, AGN heated material and the surrounding gas is important for determining the core abundances in clusters. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: 18 pages, 16 figures, submitted to MNRAS