Showing 1–50 of 85 results for author: Jones, S

Beyond Point Masses. II. Non-Keplerian Shape Effects are Detectable in Several TNO Binaries

Authors: Benjamin C. N. Proudfoot, Darin A. Ragozzine, Meagan L. Thatcher, Will Grundy, Dallin J. Spencer, Tahina M. Alailima, Sawyer Allen, Penelope C. Bowden, Susanne Byrd, Conner D. Camacho, Gibson H. Campbell, Edison P. Carlisle, Jacob A. Christensen, Noah K. Christensen, Kaelyn Clement, Benjamin J. Derieg, Mara K. Dille, Cristian Dorrett, Abigail L. Ellefson, Taylor S. Fleming, N. J. Freeman, Ethan J. Gibson, William G. Giforos, Jacob A. Guerrette, Olivia Haddock , et al. (38 additional authors not shown)

Abstract: About 40 transneptunian binaries (TNBs) have fully determined orbits with about 10 others being solved except for breaking the mirror ambiguity. Despite decades of study almost all TNBs have only ever been analyzed with a model that assumes perfect Keplerian motion (e.g., two point masses). In reality, all TNB systems are non-Keplerian due to non-spherical shapes, possible presence of undetected s… ▽ More About 40 transneptunian binaries (TNBs) have fully determined orbits with about 10 others being solved except for breaking the mirror ambiguity. Despite decades of study almost all TNBs have only ever been analyzed with a model that assumes perfect Keplerian motion (e.g., two point masses). In reality, all TNB systems are non-Keplerian due to non-spherical shapes, possible presence of undetected system components, and/or solar perturbations. In this work, we focus on identifying candidates for detectable non-Keplerian motion based on sample of 45 well-characterized binaries. We use MultiMoon, a non-Keplerian Bayesian inference tool, to analyze published relative astrometry allowing for non-spherical shapes of each TNB system's primary. We first reproduce the results of previous Keplerian fitting efforts with MultiMoon, which serves as a comparison for the non-Keplerian fits and confirms that these fits are not biased by the assumption of a Keplerian orbit. We unambiguously detect non-Keplerian motion in 8 TNB systems across a range of primary radii, mutual orbit separations, and system masses. As a proof of concept for non-Keplerian fitting, we perform detailed fits for (66652) Borasisi-Pabu, possibly revealing a $J_2 \approx 0.44$, implying Borasisi (and/or Pabu) may be a contact binary or an unresolved compact binary. However, full confirmation of this result will require new observations. This work begins the next generation of TNB analyses that go beyond the point mass assumption to provide unique and valuable information on the physical properties of TNBs with implications for their formation and evolution. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: Accepted to AJ

Journal ref: AJ 167 144 (2024)

TOI-2015b: A Warm Neptune with Transit Timing Variations Orbiting an Active mid M Dwarf

Authors: Sinclaire E. Jones, Gudmundur Stefansson, Kento Masuda, Jessica E. Libby-Roberts, Cristilyn N. Gardner, Rae Holcomb, Corey Beard, Paul Robertson, Caleb I. Cañas, Suvrath Mahadevan, Shubham Kanodia, Andrea S. J. Lin, Henry A. Kobulnicky, Brock A. Parker, Chad F. Bender, William D. Cochran, Scott A. Diddams, Rachel B. Fernandes, Arvind F. Gupta, Samuel Halverson, Suzanne L. Hawley, Fred R. Hearty, Leslie Hebb, Adam Kowalski, Jack Lubin , et al. (7 additional authors not shown)

Abstract: We report the discovery of a close-in ($P_{\mathrm{orb}} = 3.349\:\mathrm{days}$) warm Neptune with clear transit timing variations (TTVs) orbiting the nearby ($d=47.3\:\mathrm{pc}$) active M4 star, TOI-2015. We characterize the planet's properties using TESS photometry, precise near-infrared radial velocities (RV) with the Habitable-zone Planet Finder (HP) Spectrograph, ground-based photometry, a… ▽ More We report the discovery of a close-in ($P_{\mathrm{orb}} = 3.349\:\mathrm{days}$) warm Neptune with clear transit timing variations (TTVs) orbiting the nearby ($d=47.3\:\mathrm{pc}$) active M4 star, TOI-2015. We characterize the planet's properties using TESS photometry, precise near-infrared radial velocities (RV) with the Habitable-zone Planet Finder (HP) Spectrograph, ground-based photometry, and high-contrast imaging. A joint photometry and RV fit yields a radius $R_p~=~3.37_{-0.20}^{+0.15} \:\mathrm{R_\oplus}$, mass $m_p~=~16.4_{-4.1}^{+4.1}\:\mathrm{M_\oplus}$, and density $ρ_p~=~2.32_{-0.37}^{+0.38} \:\mathrm{g cm^{-3}}$ for TOI-2015b, suggesting a likely volatile-rich planet. The young, active host star has a rotation period of $P_{\mathrm{rot}}~=~8.7 \pm~0.9~\mathrm{days}$ and associated rotation-based age estimate of $1.1~\pm~0.1\:\mathrm{Gyr}$. Though no other transiting planets are seen in the TESS data, the system shows clear TTVs of super period $P_{\mathrm{sup}}~\approx~430\:\mathrm{days}$ and amplitude $\sim$$100\:\mathrm{minutes}$. After considering multiple likely period ratio models, we show an outer planet candidate near a 2:1 resonance can explain the observed TTVs while offering a dynamically stable solution. However, other possible two-planet solutions -- including 3:2 and 4:3 resonance -- cannot be conclusively excluded without further observations. Assuming a 2:1 resonance in the joint TTV-RV modeling suggests a mass of $m_b~=~13.3_{-4.5}^{+4.7}\:\mathrm{M_\oplus}$ for TOI-2015b and $m_c~=~6.8_{-2.3}^{+3.5}\:\mathrm{M_\oplus}$ for the outer candidate. Additional transit and RV observations will be beneficial to explicitly identify the resonance and further characterize the properties of the system. △ Less

Submitted 9 May, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: 29 pages, 15 figures, 6 tables. Accepted for publication in The Astronomical Journal

Ensemble Forecasts of Solar Wind Connectivity to 1 Rs using ADAPT-WSA

Authors: D. E. da Silva, S. Wallace, C. N. Arge, S. Jones

Abstract: The solar wind which arrives at any location in the solar system is, in principle, relatable to the outflow of solar plasma from a single source location. This source location, itself usually being part of a larger coronal hole, is traceable to 1 Rs along the Sun's magnetic field, in which the entire path from 1 Rs to a location in the heliosphere is referred to as the solar wind connectivity. Whi… ▽ More The solar wind which arrives at any location in the solar system is, in principle, relatable to the outflow of solar plasma from a single source location. This source location, itself usually being part of a larger coronal hole, is traceable to 1 Rs along the Sun's magnetic field, in which the entire path from 1 Rs to a location in the heliosphere is referred to as the solar wind connectivity. While not directly measurable, the connectivity between the near-Earth solar wind is of particular importance to space weather. The solar wind solar source region can be obtained by leveraging near-sun magnetic field models and a model of the interplanetary solar wind. In this article we present a method for making an ensemble forecast of the connectivity presented as a probability distribution obtained from a weighted collection of individual forecasts from the combined Air Force Data Assimilative Photospheric Flux Transport - Wang Sheeley Arge (ADAPT-WSA) model. The ADAPT model derives the photospheric magnetic field from synchronic magnetogram data, using flux transport physics and ongoing data assimilation processes. The WSA model uses a coupled set of potential field type models to derive the coronal magnetic field, and an empirical relationship to derive the terminal solar wind speed observed at Earth. Our method produces an arbitrary 2D probability distribution capable of reflecting complex source configurations with minimal assumptions about the distribution structure, prepared in a computationally efficient manner. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Comments: Accepted to the journal "Space Weather"

Type Ia Supernova Nucleosynthesis: Metallicity-Dependent Yields

Authors: James D Keegans, Marco Pignatari, Richard J Stancliffe, Claudia Travaglio, Samuel Jones, Brad K Gibson, Dean M Townsley, Broxton J Miles, Ken J Shen, Gareth Few

Abstract: Type Ia supernova explosions (SNIa) are fundamental sources of elements for the chemical evolution of galaxies. They efficiently produce intermediate-mass (with Z between 11 and 20) and iron group elements - for example, about 70% of the solar iron is expected to be made by SNIa. In this work, we calculate complete abundance yields for 39 models of SNIa explosions, based on three progenitors - a 1… ▽ More Type Ia supernova explosions (SNIa) are fundamental sources of elements for the chemical evolution of galaxies. They efficiently produce intermediate-mass (with Z between 11 and 20) and iron group elements - for example, about 70% of the solar iron is expected to be made by SNIa. In this work, we calculate complete abundance yields for 39 models of SNIa explosions, based on three progenitors - a 1.4M deflagration detonation model, a 1.0 double detonation model and a 0.8 M double detonation model - and 13 metallicities, with 22Ne mass fractions of 0, 1x10-7, 1x10-6, 1x10-5, 1x10-4, 1x10-3, 2x10-3, 5x10-3, 1x10-2, 1.4x10-2, 5x10-2, and 0.1 respectively. Nucleosynthesis calculations are done using the NuGrid suite of codes, using a consistent nuclear reaction network between the models. Complete tables with yields and production factors are provided online at Zenodo: Yields. We discuss the main properties of our yields in the light of the present understanding of SNIa nucleosynthesis, depending on different progenitor mass and composition. Finally, we compare our results with a number of relevant models from the literature. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 42 pages, 21 figures. Accepted for publication in ApJS 21-06-23

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

TOI-3785 b: A Low-Density Neptune Orbiting an M2-Dwarf Star

Authors: Luke C. Powers, Jessica Libby-Roberts, Andrea S. J. Lin, Caleb I. Cañas, Shubham Kanodia, Suvrath Mahadevan, Joe P. Ninan, Guðmundur Stefánsson, Arvind F. Gupta, Sinclaire Jones, Henry A. Kobulnicky, Andrew Monson, Brock A. Parker, Tera N. Swaby, Chad F. Bender, William D. Cochran, Leslie Hebb, Andrew J. Metcalf, Paul Robertson, Christian Schwab, John Wisniewski, Jason T. Wright

Abstract: Using both ground-based transit photometry and high-precision radial velocity (RV) spectroscopy, we confirm the planetary nature of TOI-3785 b. This transiting Neptune orbits an M2-Dwarf star with a period of ~4.67 days, a planetary radius of 5.14 +/- 0.16 Earth Radii, a mass of 14.95 +4.10, -3.92 Earth Masses, and a density of 0.61 +0.18, -0.17 g/cm^3. TOI-3785 b belongs to a rare population of N… ▽ More Using both ground-based transit photometry and high-precision radial velocity (RV) spectroscopy, we confirm the planetary nature of TOI-3785 b. This transiting Neptune orbits an M2-Dwarf star with a period of ~4.67 days, a planetary radius of 5.14 +/- 0.16 Earth Radii, a mass of 14.95 +4.10, -3.92 Earth Masses, and a density of 0.61 +0.18, -0.17 g/cm^3. TOI-3785 b belongs to a rare population of Neptunes (4 Earth Radii < Rp < 7 Earth Radii) orbiting cooler, smaller M-dwarf host stars, of which only ~10 have been confirmed. By increasing the number of confirmed planets, TOI-3785 b offers an opportunity to compare similar planets across varying planetary and stellar parameter spaces. Moreover, with a high transmission spectroscopy metric (TSM) of ~150 combined with a relatively cool equilibrium temperature of 582 +/- 16 K and an inactive host star, TOI-3785 b is one of the more promising low-density M-dwarf Neptune targets for atmospheric follow-up. Future investigation into atmospheric mass loss rates of TOI-3785 b may yield new insights into the atmospheric evolution of these low-mass gas planets around M-dwarfs. △ Less

Submitted 12 July, 2023; v1 submitted 10 April, 2023; originally announced April 2023.

Comments: 22 pages, 6 figures, 6 tables, Published to AJ. arXiv admin note: text overlap with arXiv:2301.10837

Journal ref: The Astronomical Journal, Volume 166, Number 2, Year 2023

Prediction and Verification of Parker Solar Probe Solar Wind Sources at 13.3 R$_\odot$

Authors: Samuel T. Badman, Pete Riley, Shaela I. Jones, Tae K. Kim, Robert C. Allen, C. Nick Arge, Stuart D. Bale, Carl J. Henney, Justin C. Kasper, Parisa Mostafavi, Nikolai V. Pogorelov, Nour E. Raouafi, Michael L. Stevens, J. L. Verniero

Abstract: Drawing connections between heliospheric spacecraft and solar wind sources is a vital step in understanding the evolution of the solar corona into the solar wind and contextualizing \textit{in situ} timeseries. Furthermore, making advanced predictions of this linkage for ongoing heliospheric missions, such as Parker Solar Probe (PSP), is necessary for achieving useful coordinated remote observatio… ▽ More Drawing connections between heliospheric spacecraft and solar wind sources is a vital step in understanding the evolution of the solar corona into the solar wind and contextualizing \textit{in situ} timeseries. Furthermore, making advanced predictions of this linkage for ongoing heliospheric missions, such as Parker Solar Probe (PSP), is necessary for achieving useful coordinated remote observations and maximizing scientific return. The general procedure for estimating such connectivity is straightforward (i.e. magnetic field line tracing in a coronal model) but validating the resulting estimates difficult due to the lack of an independent ground truth and limited model constraints. In its most recent orbits, PSP has reached perihelia of 13.3$R_\odot$ and moreover travels extremely fast prograde relative to the solar surface, covering over 120 degrees longitude in three days. Here we present footpoint predictions and subsequent validation efforts for PSP Encounter 10, the first of the 13.3$R_\odot$ orbits, which occurred in November 2021. We show that the longitudinal dependence of \textit{in situ} plasma data from these novel orbits provides a powerful method of footpoint validation. With reference to other encounters, we also illustrate that the conditions under which source mapping is most accurate for near-ecliptic spacecraft (such as PSP) occur when solar activity is low, but also requires that the heliospheric current sheet is strongly warped by mid-latitude or equatorial coronal holes. Lastly, we comment on the large-scale coronal structure implied by the Encounter 10 mapping, highlighting an empirical equatorial cut of the Alfvèn surface consisting of localized protrusions above unipolar magnetic separatrices. △ Less

Submitted 29 March, 2023; v1 submitted 8 March, 2023; originally announced March 2023.

Comments: 33 Pages, 7 Figures, JGR Space Physics, Published Online 2023/3/28, In Final Production

TOI-3984 A b and TOI-5293 A b: two temperate gas giants transiting mid-M dwarfs in wide binary systems

Authors: Caleb I. Cañas, Shubham Kanodia, Jessica Libby-Roberts, Andrea S. J. Lin, Maria Schutte, Luke Powers, Sinclaire Jones, Andrew Monson, Songhu Wang, Guðmundur Stefánsson, William D. Cochran, Paul Robertson, Suvrath Mahadevan, Adam F. Kowalski, John Wisniewski, Brock A. Parker, Alexander Larsen, Franklin A. L. Chapman, Henry A. Kobulnicky, Arvind F. Gupta, Mark E. Everett, Bryan Edward Penprase, Gregory Zeimann, Corey Beard, Chad F. Bender , et al. (8 additional authors not shown)

Abstract: We confirm the planetary nature of two gas giants discovered by TESS to transit M dwarfs with stellar companions at wide separations. TOI-3984 A ($J=11.93$) is an M4 dwarf hosting a short-period ($4.353326 \pm 0.000005$ days) gas giant ($M_p=0.14\pm0.03~\mathrm{M_{J}}$ and $R_p=0.71\pm0.02~\mathrm{R_{J}}$) with a wide separation white dwarf companion. TOI-5293 A ($J=12.47$) is an M3 dwarf hosting… ▽ More We confirm the planetary nature of two gas giants discovered by TESS to transit M dwarfs with stellar companions at wide separations. TOI-3984 A ($J=11.93$) is an M4 dwarf hosting a short-period ($4.353326 \pm 0.000005$ days) gas giant ($M_p=0.14\pm0.03~\mathrm{M_{J}}$ and $R_p=0.71\pm0.02~\mathrm{R_{J}}$) with a wide separation white dwarf companion. TOI-5293 A ($J=12.47$) is an M3 dwarf hosting a short-period ($2.930289 \pm 0.000004$ days) gas giant ($M_p=0.54\pm0.07~\mathrm{M_{J}}$ and $R_p=1.06\pm0.04~\mathrm{R_{J}}$) with a wide separation M dwarf companion. We characterize both systems using a combination of ground-based and space-based photometry, speckle imaging, and high-precision radial velocities from the Habitable-zone Planet Finder and NEID spectrographs. TOI-3984 A b ($T_{eq}=563\pm15$ K and $\mathrm{TSM}=138_{-27}^{+29}$) and TOI-5293 A b ($T_{eq}=675_{-30}^{+42}$ K and $\mathrm{TSM}=92\pm14$) are two of the coolest gas giants among the population of hot Jupiter-sized gas planets orbiting M dwarfs and are favorable targets for atmospheric characterization of temperate gas giants and three-dimensional obliquity measurements to probe system architecture and migration scenarios. △ Less

Submitted 27 June, 2023; v1 submitted 15 February, 2023; originally announced February 2023.

Comments: Published in AJ, 46 pages, 16 figures, 6 tables, updated to reflect published version. arXiv admin note: substantial text overlap with arXiv:2201.09963

Journal ref: AJ, 166, 30 (2023)

Implications of Different Solar Photospheric Flux-Transport Models for Global Coronal and Heliospheric Modeling

Authors: Graham Barnes, Marc L. DeRosa, Shaela I. Jones, Charles N. Arge, Carl J. Henney, Mark C. M. Cheung

Abstract: The concept of surface-flux transport (SFT) is commonly used in evolving models of the large-scale solar surface magnetic field. These photospheric models are used to determine the large-scale structure of the overlying coronal magnetic field, as well as to make predictions about the fields and flows that structure the solar wind. We compare predictions from two SFT models for the solar wind, open… ▽ More The concept of surface-flux transport (SFT) is commonly used in evolving models of the large-scale solar surface magnetic field. These photospheric models are used to determine the large-scale structure of the overlying coronal magnetic field, as well as to make predictions about the fields and flows that structure the solar wind. We compare predictions from two SFT models for the solar wind, open magnetic field footpoints, and the presence of coronal magnetic null points throughout various phases of a solar activity cycle, focusing on the months of April in even-numbered years between 2012 and 2020, inclusive. We find that there is a solar cycle dependence to each of the metrics considered, but there is not a single phase of the cycle in which all the metrics indicate good agreement between the models. The metrics also reveal large, transient differences between the models when a new active region is rotating into the assimilation window. The evolution of the surface flux is governed by a combination of large scale flows and comparatively small scale motions associated with convection. Because the latter flows evolve rapidly, there are intervals during which their impact on the surface flux can only be characterized in a statistical sense, thus their impact is modeled by introducing a random evolution that reproduces the typical surface flux evolution. We find that the differences between the predicted properties are dominated by differences in the model assumptions and implementation, rather than selection of a particular realization of the random evolution. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: Accepted for publication in The Astrophysical Journal

GJ 3929: High Precision Photometric and Doppler Characterization of an Exo-Venus and its Hot, Mini-Neptune-mass Companion

Authors: Corey Beard, Paul Robertson, Shubham Kanodia, Jack Lubin, Caleb I. Cañas, Arvind F. Gupta, Rae Holcomb, Sinclaire Jones, Jessica E. Libby-Roberts, Andrea S. J. Lin, Suvrath Mahadevan, Guðmundur Stefánsson, Chad F. Bender, Cullen H. Blake, William D. Cochran, Michael Endl, Mark Everett, Eric B. Ford, Connor Fredrick, Samuel Halverson, Leslie Hebb, Dan Li, Sarah E. Logsdon, Jacob Luhn, Michael W. McElwain , et al. (9 additional authors not shown)

Abstract: We detail the follow up and characterization of a transiting exo-Venus identified by TESS, GJ 3929b, (TOI-2013b) and its non-transiting companion planet, GJ 3929c (TOI-2013c). GJ 3929b is an Earth-sized exoplanet in its star's Venus-zone (P$_{b}$ = 2.616272 $\pm$ 0.000005 days; S$_{b}$ = 17.3$^{+0.8}_{-0.7}$ S$_{\oplus}$) orbiting a nearby M dwarf. GJ 3929c is most likely a non-transiting sub-Nept… ▽ More We detail the follow up and characterization of a transiting exo-Venus identified by TESS, GJ 3929b, (TOI-2013b) and its non-transiting companion planet, GJ 3929c (TOI-2013c). GJ 3929b is an Earth-sized exoplanet in its star's Venus-zone (P$_{b}$ = 2.616272 $\pm$ 0.000005 days; S$_{b}$ = 17.3$^{+0.8}_{-0.7}$ S$_{\oplus}$) orbiting a nearby M dwarf. GJ 3929c is most likely a non-transiting sub-Neptune. Using the new, ultra-precise NEID spectrometer on the WIYN 3.5 m Telescope at Kitt Peak National Observatory, we are able to modify the mass constraints of planet b reported in previous works and consequently improve the significance of the mass measurement to almost 4$σ$ confidence (M$_{b}$ = 1.75 $\pm$ 0.45 M$_{\oplus}$). We further adjust the orbital period of planet c from its alias at 14.30 $\pm$ 0.03 days to the likely true period of 15.04 $\pm$ 0.03 days, and adjust its minimum mass to m$\sin i$ = 5.71 $\pm$ 0.92 M$_{\oplus}$. Using the diffuser-assisted ARCTIC imager on the ARC 3.5 m telescope at Apache Point Observatory, in addition to publicly available TESS and LCOGT photometry, we are able to constrain the radius of planet b to R$_{p}$ = 1.09 $\pm$ 0.04 R$_{\oplus}$. GJ 3929b is a top candidate for transmission spectroscopy in its size regime (TSM = 14 $\pm$ 4), and future atmospheric studies of GJ 3929b stand to shed light on the nature of small planets orbiting M dwarfs. △ Less

Submitted 30 July, 2022; v1 submitted 21 July, 2022; originally announced July 2022.

Comments: Accepted for publication in ApJ. arXiv admin note: text overlap with arXiv:2204.09063

New insights into the first two PSP solar encounters enabled by modeling analysis with ADAPT-WSA

Authors: Samantha Wallace, Shaela I. Jones, C. Nick Arge, Nicholeen Viall, Carl J. Henney

Abstract: Parker Solar Probes's (PSP)'s unique orbital path allows us to observe the solar wind closer to the Sun than ever before. Essential to advancing our knowledge of solar wind and energetic particle formation is identifying the sources of PSP observations. We report on results for the first two PSP solar encounters derived using the Wang-Sheeley-Arge (WSA) model driven by Air Force Data Assimilative… ▽ More Parker Solar Probes's (PSP)'s unique orbital path allows us to observe the solar wind closer to the Sun than ever before. Essential to advancing our knowledge of solar wind and energetic particle formation is identifying the sources of PSP observations. We report on results for the first two PSP solar encounters derived using the Wang-Sheeley-Arge (WSA) model driven by Air Force Data Assimilative Photospheric Flux Transport (ADAPT) model maps. We derive the coronal magnetic field and the 1 Rs source regions of the PSP-observed solar wind. We validate our results with the solar wind speed and magnetic polarity observed at PSP. When modeling results are very reliable, we derive time series of model-derived spacecraft separation from the heliospheric current sheet, magnetic expansion factor, coronal hole boundary distance, and photospheric field strength along the field lines estimated to be connected to the spacecraft. We present new results for Encounter 1, which show time evolution of the far-side mid-latitude coronal hole that PSP co-rotates with. We discuss how this evolution coincides with solar wind speed, density, and temperature observed at the spacecraft. During Encounter 2, a new active region emerges on the far-side, making it difficult to model. We show that ADAPT-WSA output agrees well with PSP observations once this active region rotates onto the near-side, allowing us to reliably estimate the solar wind sources retrospectively for most of the encounter. We close with ways in which coronal modeling enables scientific interpretation of these encounters that would otherwise not have been possible. △ Less

Submitted 25 May, 2022; originally announced May 2022.

Comments: 25 pages, 13 figures, accepted by ApJ

Solar Polar Flux Redistribution based on Observed Coronal Holes

Authors: Samuel J. Schonfeld, Carl J. Henney, Shaela I. Jones, Charles N. Arge

Abstract: We explore the use of observed polar coronal holes (CHs) to constrain the flux distribution within the polar regions of global solar magnetic field maps in the absence of reliable quality polar field observations. Global magnetic maps, generated by the Air Force Data Assimilative Photospheric flux Transport (ADAPT) model, are modified to enforce field unipolarity thresholds both within and outside… ▽ More We explore the use of observed polar coronal holes (CHs) to constrain the flux distribution within the polar regions of global solar magnetic field maps in the absence of reliable quality polar field observations. Global magnetic maps, generated by the Air Force Data Assimilative Photospheric flux Transport (ADAPT) model, are modified to enforce field unipolarity thresholds both within and outside observed CH boundaries. The polar modified and unmodified maps are used to drive Wang-Sheeley-Arge (WSA) models of the corona and solar wind (SW). The WSA predicted CHs are compared with the observations, and SW predictions at the WIND and Ulysses spacecraft are also used to provide context for the new polar modified maps. We find that modifications of the polar flux never worsen and typically improve both the CH and SW predictions. We also confirm the importance of the choice of the domain over which WSA generates the coronal magnetic field solution but find that solutions optimized for one location in the heliosphere can worsen predictions at other locations. Finally, we investigate the importance of low-latitude (i.e., active region) magnetic fields in setting the boundary of polar CHs, determining that they have at least as much impact as the polar fields themselves. △ Less

Submitted 28 April, 2022; originally announced April 2022.

Comments: 23 pages, 10 figures, 2 tables, 4 appendices, 1 data repository. Accepted for publication in ApJ

TOI-1696 and TOI-2136: Constraining the Masses of Two Mini-Neptunes with HPF

Authors: Corey Beard, Paul Robertson, Shubham Kanodia, Jessica Libby-Roberts, Caleb I. Canas, Arvind F. Gupta, Rae Holcomb, Sinclaire Jones, Henry A. Kobulnicky, Andrea S. J. Lin, Jack Lubin, Marissa Maney, Brock A. Parker, Gudmundur Stefansson, William D. Cochran, Michael Endl, Leslie Hebb, Suvrath Mahadevan, John Wisniewski, Chad F. Bender, Scott A. Diddams, Mark Everett, Connor Fredrick, Samuel Halverson, Fred Hearty , et al. (7 additional authors not shown)

Abstract: We present the validation of two planets orbiting M dwarfs, TOI-1696b and TOI-2136b. Both planets are mini-Neptunes orbiting nearby stars, making them promising prospects for atmospheric characterization with the James Webb Space Telescope. We validated the planetary nature of both candidates using high contrast imaging, ground-based photometry, and near-infrared radial velocities. Adaptive Optics… ▽ More We present the validation of two planets orbiting M dwarfs, TOI-1696b and TOI-2136b. Both planets are mini-Neptunes orbiting nearby stars, making them promising prospects for atmospheric characterization with the James Webb Space Telescope. We validated the planetary nature of both candidates using high contrast imaging, ground-based photometry, and near-infrared radial velocities. Adaptive Optics images were taken using the ShARCS camera on the 3 m Shane Telescope. Speckle images were taken using the NN-Explore Exoplanet Stellar Speckle Imager on the WIYN 3.5 m telescope. Radii and orbital ephemerides were refined using a combination of TESS, the diffuser-assisted ARCTIC imager on the 3.5m ARC telescope at Apache Point Observatory, and the 0.6 m telescope at Red Buttes Observatory. We obtained radial velocities using the Habitable-Zone Planet Finder on the 10 m Hobby-Eberly Telescope, which enabled us to place upper limits on the masses of both transiting planets. TOI-1696b (P = 2.5 days; R$_{p}$ = 3.24 R$_{\oplus}$; M$_{p}$ $<$ 56.6 M$_{\oplus}$) falls into a sparsely-populated region of parameter space considering its host star's temperature (T$_{\rm{eff}}$ = 3168 K, M4.5), as planets of its size are quite rare around mid to late M dwarfs. On the other hand, TOI-2136b (P = 7.85 days; R$_{p}$ = 2.09 R$_{\oplus}$; M$_{p}$ $<$ 15.0 M$_{\oplus}$) is an excellent candidate for atmospheric follow-up with JWST. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Comments: Accepted for Publication in the Astronomical Journal

Radioactive nuclei in the early Solar System: analysis of the 15 isotopes produced by core-collapse supernovae

Authors: Thomas V. Lawson, Marco Pignatari, Richard J. Stancliffe, Jacqueline den Hartogh, Sam Jones, Chris L. Fryer, Brad K. Gibson, Maria Lugaro

Abstract: Short-lived radioactive isotopes (SLRs) with half-lives between 0.1 to 100 Myr can be used to probe the origin of the Solar System. In this work, we examine the core-collapse supernovae production of the 15 SLRs produced: $^{26}$Al, $^{36}$Cl, $^{41}$Ca, $^{53}$Mn, $^{60}$Fe, $^{92}$Nb, $^{97}$Tc, $^{98}$Tc, $^{107}$Pd, $^{126}$Sn, $^{129}$I, $^{135}$Cs, $^{146}$Sm, $^{182}$Hf, and $^{205}$Pb. We… ▽ More Short-lived radioactive isotopes (SLRs) with half-lives between 0.1 to 100 Myr can be used to probe the origin of the Solar System. In this work, we examine the core-collapse supernovae production of the 15 SLRs produced: $^{26}$Al, $^{36}$Cl, $^{41}$Ca, $^{53}$Mn, $^{60}$Fe, $^{92}$Nb, $^{97}$Tc, $^{98}$Tc, $^{107}$Pd, $^{126}$Sn, $^{129}$I, $^{135}$Cs, $^{146}$Sm, $^{182}$Hf, and $^{205}$Pb. We probe the impact of the uncertainties of the core-collapse explosion mechanism by examining a collection of 62 core-collapse models with initial masses of 15, 20, and 25M$_{\odot}$, explosion energies between 3.4$\times$10$^{50}$ and 1.8$\times$10$^{52}$ ergs and compact remnant masses between 1.5M$_{\odot}$and 4.89M$_{\odot}$. We identify the impact of both explosion energy and remnant mass on the final yields of the SLRs. Isotopes produced within the innermost regions of the star, such as $^{92}$Nb and $^{97}$Tc, are the most affected by the remnant mass, $^{92}$Nb varying by five orders of magnitude. Isotopes synthesised primarily in explosive C-burning and explosive He-burning, such as $^{60}$Fe, are most affected by explosion energies. $^{60}$Fe increases by two orders of magnitude from the lowest to the highest explosion energy in the 15M$_{\odot}$model. The final yield of each examined SLR is used to compare to literature models. △ Less

Submitted 20 March, 2022; originally announced March 2022.

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 511, Issue 1, March 2022, Pages 886 902

White Paper on Light Sterile Neutrino Searches and Related Phenomenology

Authors: M. A. Acero, C. A. Argüelles, M. Hostert, D. Kalra, G. Karagiorgi, K. J. Kelly, B. Littlejohn, P. Machado, W. Pettus, M. Toups, M. Ross-Lonergan, A. Sousa, P. T. Surukuchi, Y. Y. Y. Wong, W. Abdallah, A. M. Abdullahi, R. Akutsu, L. Alvarez-Ruso, D. S. M. Alves, A. Aurisano, A. B. Balantekin, J. M. Berryman, T. Bertólez-Martínez, J. Brunner, M. Blennow , et al. (147 additional authors not shown)

Abstract: This white paper provides a comprehensive review of our present understanding of experimental neutrino anomalies that remain unresolved, charting the progress achieved over the last decade at the experimental and phenomenological level, and sets the stage for future programmatic prospects in addressing those anomalies. It is purposed to serve as a guiding and motivational "encyclopedic" reference,… ▽ More This white paper provides a comprehensive review of our present understanding of experimental neutrino anomalies that remain unresolved, charting the progress achieved over the last decade at the experimental and phenomenological level, and sets the stage for future programmatic prospects in addressing those anomalies. It is purposed to serve as a guiding and motivational "encyclopedic" reference, with emphasis on needs and options for future exploration that may lead to the ultimate resolution of the anomalies. We see the main experimental, analysis, and theory-driven thrusts that will be essential to achieving this goal being: 1) Cover all anomaly sectors -- given the unresolved nature of all four canonical anomalies, it is imperative to support all pillars of a diverse experimental portfolio, source, reactor, decay-at-rest, decay-in-flight, and other methods/sources, to provide complementary probes of and increased precision for new physics explanations; 2) Pursue diverse signatures -- it is imperative that experiments make design and analysis choices that maximize sensitivity to as broad an array of these potential new physics signatures as possible; 3) Deepen theoretical engagement -- priority in the theory community should be placed on development of standard and beyond standard models relevant to all four short-baseline anomalies and the development of tools for efficient tests of these models with existing and future experimental datasets; 4) Openly share data -- Fluid communication between the experimental and theory communities will be required, which implies that both experimental data releases and theoretical calculations should be publicly available; and 5) Apply robust analysis techniques -- Appropriate statistical treatment is crucial to assess the compatibility of data sets within the context of any given model. △ Less

Submitted 17 May, 2023; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021 by the NF02 Topical Group (Understanding Experimental Neutrino Anomalies). Submitted to J. Phys. G as a Major Report

TOI-3757 b: A low density gas giant orbiting a solar-metallicity M dwarf

Authors: Shubham Kanodia, Jessica Libby-Roberts, Caleb I. Canas, Joe P. Ninan, Suvrath Mahadevan, Gudmundur Stefansson, Andrea S. J. Lin, Sinclaire Jones, Andrew Monson, Brock A. Parker, Henry A. Kobulnicky, Tera N. Swaby, Luke Powers, Corey Beard, Chad F. Bender, Cullen H. Blake, William D. Cochran, Jiayin Dong, Scott A. Diddams, Connor Fredrick, Arvind F. Gupta, Samuel Halverson, Fred Hearty, Sarah E. Logsdon, Andrew J. Metcalf , et al. (10 additional authors not shown)

Abstract: We present the discovery of a new Jovian-sized planet, TOI-3757 b, the lowest density planet orbiting an M dwarf (M0V). It orbits a solar-metallicity M dwarf discovered using TESS photometry and confirmed with precise radial velocities (RV) from HPF and NEID. With a planetary radius of $12.0^{+0.4}_{-0.5}$ $R_{\oplus}$ and mass of $85.3^{+8.8}_{-8.7}$ $M_{\oplus}$, not only does this object add to… ▽ More We present the discovery of a new Jovian-sized planet, TOI-3757 b, the lowest density planet orbiting an M dwarf (M0V). It orbits a solar-metallicity M dwarf discovered using TESS photometry and confirmed with precise radial velocities (RV) from HPF and NEID. With a planetary radius of $12.0^{+0.4}_{-0.5}$ $R_{\oplus}$ and mass of $85.3^{+8.8}_{-8.7}$ $M_{\oplus}$, not only does this object add to the small sample of gas giants ($\sim 10$) around M dwarfs, but also, its low density ($ρ=$ $0.27^{+0.05}_{-0.04}$ $\textrm{g~cm}^{-3}$) provides an opportunity to test theories of planet formation. We present two hypotheses to explain its low density; first, we posit that the low metallicity of its stellar host ($\sim$ 0.3 dex lower than the median metallicity of M dwarfs hosting gas giants) could have played a role in the delayed formation of a solid core massive enough to initiate runaway accretion. Second, using the eccentricity estimate of $0.14 \pm 0.06$ we determine it is also plausible for tidal heating to at least partially be responsible for inflating the radius of TOI-3757b b. The low density and large scale height of TOI-3757 b makes it an excellent target for transmission spectroscopy studies of atmospheric escape and composition (TSM $\sim$ 190). We use HPF to perform transmission spectroscopy of TOI-3757 b using the helium 10830 Å~ line. Doing this, we place an upper limit of 6.9 \% (with 90\% confidence) on the maximum depth of the absorption from the metastable transition of He at $\sim$ 10830 Å, which can help constraint the atmospheric mass loss rate in this energy limited regime. △ Less

Submitted 5 August, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: AJ. arXiv admin note: text overlap with arXiv:2107.13670

Journal ref: The Astronomical Journal, Volume 164, Number 3 (2022 AJ 164 81)

Constraining Global Coronal Models with Multiple Independent Observables

Authors: Samuel T. Badman, David H. Brooks, Nicolas Poirier, Harry P. Warren, Gordon Petrie, Alexis P. Rouillard, C. Nick Arge, Stuart D. Bale, Diego de Pablos Aguero, Louise Harra, Shaela I. Jones, Athanasios Kouloumvakos, Pete Riley, Olga Panasenco, Marco Velli, Samantha Wallace

Abstract: Global coronal models seek to produce an accurate physical representation of the Sun's atmosphere which can be used, for example, to drive space weather models. Assessing their accuracy is a complex task and there are multiple observational pathways to provide constraints and tune model parameters. Here, we combine several such independent constraints, defining a model-agnostic framework for stand… ▽ More Global coronal models seek to produce an accurate physical representation of the Sun's atmosphere which can be used, for example, to drive space weather models. Assessing their accuracy is a complex task and there are multiple observational pathways to provide constraints and tune model parameters. Here, we combine several such independent constraints, defining a model-agnostic framework for standardized comparison. We require models to predict the distribution of coronal holes at the photosphere, and neutral line topology at the model outer boundary. We compare these predictions to extreme ultraviolet (EUV) observations of coronal hole locations, white-light Carrington maps of the streamer belt and the magnetic sector structure measured \textit{in situ} by Parker Solar Probe and 1AU spacecraft. We study these metrics for Potential Field Source Surface (PFSS) models as a function of source surface height and magnetogram choice, as well as comparing to the more physical Wang-Sheeley-Arge (WSA) and the Magnetohydrodynamics Algorithm outside a Sphere (MAS) models. We find that simultaneous optimization of PFSS models to all three metrics is not currently possible, implying a trade-off between the quality of representation of coronal holes and streamer belt topology. WSA and MAS results show the additional physics they include addresses this by flattening the streamer belt while maintaining coronal hole sizes, with MAS also improving coronal hole representation relative to WSA. We conclude that this framework is highly useful for inter- and intra-model comparisons. Integral to the framework is the standardization of observables required of each model, evaluating different model aspects. △ Less

Submitted 14 April, 2022; v1 submitted 27 January, 2022; originally announced January 2022.

Comments: Accepted to ApJ 4/9/2022

TOI-3714 b and TOI-3629 b: Two gas giants transiting M dwarfs confirmed with HPF and NEID

Authors: Caleb I. Cañas, Shubham Kanodia, Chad F. Bender, Suvrath Mahadevan, Guðmundur Stefánsson, William D. Cochran, Andrea S. J. Lin, Hsiang-Chih Hwang, Luke Powers, Andrew Monson, Elizabeth M. Green, Brock A. Parker, Tera N. Swaby, Henry A. Kobulnicky, John Wisniewski, Arvind F. Gupta, Mark E. Everett, Sinclaire Jones, Benjamin Anjakos, Corey Beard, Cullen H. Blake, Scott A. Diddams, Zehao Dong, Connor Fredrick, Elnaz Hakemiamjad , et al. (14 additional authors not shown)

Abstract: We confirm the planetary nature of two gas giants discovered by TESS to transit M dwarfs. TOI-3714 ($V=15.24,~J=11.74$) is an M2 dwarf hosting a hot Jupiter ($M_p=0.70 \pm 0.03~\mathrm{M_J}$ and $R_p=1.01 \pm 0.03~\mathrm{R_J}$) on an orbital period of $2.154849 \pm 0.000001$ days with a resolved white dwarf companion. TOI-3629 ($V=14.63,~J=11.42$) is an M1 dwarf hosting a hot Jupiter (… ▽ More We confirm the planetary nature of two gas giants discovered by TESS to transit M dwarfs. TOI-3714 ($V=15.24,~J=11.74$) is an M2 dwarf hosting a hot Jupiter ($M_p=0.70 \pm 0.03~\mathrm{M_J}$ and $R_p=1.01 \pm 0.03~\mathrm{R_J}$) on an orbital period of $2.154849 \pm 0.000001$ days with a resolved white dwarf companion. TOI-3629 ($V=14.63,~J=11.42$) is an M1 dwarf hosting a hot Jupiter ($M_p=0.26 \pm 0.02~\mathrm{M_J}$ and $R_p=0.74 \pm 0.02~\mathrm{R_J}$) on an orbital period of $3.936551_{-0.000006}^{+0.000005}$ days. We characterize each transiting companion using a combination of ground-based and space-based photometry, speckle imaging, and high-precision velocimetry from the Habitable-zone Planet Finder and the NEID spectrographs. With the discovery of these two systems, there are now nine M dwarfs known to host transiting hot Jupiters. Among this population, TOI-3714 b ($T_{eq}=750\pm20$ K and $\mathrm{TSM}=98\pm7$) and TOI-3629 b ($T_{eq}=690\pm20$ K and $\mathrm{TSM}=80\pm9$) are warm gas giants amenable to additional characterization with transmission spectroscopy to probe atmospheric chemistry and, for TOI-3714, obliquity measurements to probe formation scenarios. △ Less

Submitted 15 June, 2022; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: 39 pages, 13 figures, 4 tables. Accepted for publication in AJ

Journal ref: AJ, 164, 50 (2022)

An eccentric Brown Dwarf eclipsing an M dwarf

Authors: Caleb I. Cañas, Suvrath Mahadevan, Chad F. Bender, Noah Isaac Salazar Rivera, Andrew Monson, Corey Beard, Jack Lubin, Paul Robertson, Arvind F. Gupta, William D. Cochran, Connor Fredrick, Fred Hearty, Sinclaire Jones, Shubham Kanodia, Andrea S. J. Lin, Joe P. Ninan, Lawrence W. Ramsey, Christian Schwab, Guðmundur Stefánsson

Abstract: We report the discovery of a $M=67\pm2~\mathrm{M_J}$ brown dwarf transiting the early M dwarf TOI-2119 on an eccentric orbit ($e=0.3362 \pm 0.0005$) at an orbital period of $7.200861 \pm 0.000005$ days. We confirm the brown dwarf nature of the transiting companion using a combination of ground-based and space-based photometry and high-precision velocimetry from the Habitable-zone Planet Finder. De… ▽ More We report the discovery of a $M=67\pm2~\mathrm{M_J}$ brown dwarf transiting the early M dwarf TOI-2119 on an eccentric orbit ($e=0.3362 \pm 0.0005$) at an orbital period of $7.200861 \pm 0.000005$ days. We confirm the brown dwarf nature of the transiting companion using a combination of ground-based and space-based photometry and high-precision velocimetry from the Habitable-zone Planet Finder. Detection of the secondary eclipse with TESS photometry enables a precise determination of the eccentricity and reveals the brown dwarf has a brightness temperature of $2100\pm80$ K, a value which is consistent with an early L dwarf. TOI-2119 is one of the most eccentric known brown dwarfs with $P<10$ days, possibly due to the long circularization timescales for an object orbiting an M dwarf. We assess the prospects for determining the obliquity of the host star to probe formation scenarios and the possibility of additional companions in the system using Gaia EDR3 and our radial velocities. △ Less

Submitted 25 January, 2022; v1 submitted 7 December, 2021; originally announced December 2021.

Comments: 31 pages, 8 figures, 3 tables, published in AJ, minor revisions to match published version

Journal ref: AJ, 163, 89 (2022)

A hot Mars-sized exoplanet transiting an M dwarf

Authors: Caleb I. Cañas, Suvrath Mahadevan, William D. Cochran, Chad F. Bender, Eric D. Feigelson, C. E. Harman, Ravi Kumar Kopparapu, Gabriel A. Caceres, Scott A. Diddams, Michael Endl, Eric B. Ford, Samuel Halverson, Fred Hearty, Sinclaire Jones, Shubham Kanodia, Andrea S. J. Lin, Andrew J. Metcalf, Andrew Monson, Joe P. Ninan, Lawrence W. Ramsey, Paul Robertson, Arpita Roy, Christian Schwab, Guðmundur Stefánsson

Abstract: We validate the planetary nature of an ultra-short period planet orbiting the M dwarf KOI-4777. We use a combination of space-based photometry from Kepler, high-precision, near-infrared Doppler spectroscopy from the Habitable-zone Planet Finder, and adaptive optics imaging to characterize this system. KOI-4777.01 is a Mars-sized exoplanet ($\mathrm{R}_{p}=0.51 \pm 0.03R_{\oplus}$) orbiting the hos… ▽ More We validate the planetary nature of an ultra-short period planet orbiting the M dwarf KOI-4777. We use a combination of space-based photometry from Kepler, high-precision, near-infrared Doppler spectroscopy from the Habitable-zone Planet Finder, and adaptive optics imaging to characterize this system. KOI-4777.01 is a Mars-sized exoplanet ($\mathrm{R}_{p}=0.51 \pm 0.03R_{\oplus}$) orbiting the host star every 0.412-days ($\sim9.9$-hours). This is the smallest validated ultra-short period planet known and we see no evidence for additional massive companions using our HPF RVs. We constrain the upper $3σ$ mass to $M_{p}<0.34~\mathrm{M_\oplus}$ by assuming the planet is less dense than iron. Obtaining a mass measurement for KOI-4777.01 is beyond current instrumental capabilities. △ Less

Submitted 7 December, 2021; originally announced December 2021.

Comments: 25 pages, 8 figures, accepted for publication in AJ

Journal ref: AJ, 163, 3 (2022)

First models of s process in AGB stars of solar metallicity for the stellar evolutionary code ATON with a novel stable explicit numerical solver

Authors: A. Yagüe López, D. A. García-Hernández, P. Ventura, C. L. Doherty, J. W. den Hartogh, S. W. Jones, M. Lugaro

Abstract: Aims. We describe the first s-process post-processing models for asymptotic giant branch (AGB) stars of masses 3, 4 and 5 M at solar metallicity (Z=0.018) computed using the input from the stellar evolutionary code aton. Methods. The models are computed with the new code snuppat(S-process NUcleosynthesis Post-Processing code for aton), including an advective scheme for the convective overshoot tha… ▽ More Aims. We describe the first s-process post-processing models for asymptotic giant branch (AGB) stars of masses 3, 4 and 5 M at solar metallicity (Z=0.018) computed using the input from the stellar evolutionary code aton. Methods. The models are computed with the new code snuppat(S-process NUcleosynthesis Post-Processing code for aton), including an advective scheme for the convective overshoot that leads to the formation of the main neutron source, 13C. Each model is post-processed with 3 different values of the free overshoot parameter. Included in the code snuppat is the novel Patankar-Euler-Deflhard explicit numerical solver, that we use to solve the nuclear network system of differential equations. Results. The results are compared to those from other s-process nucleosynthesis codes (Monash,fruity, and NuGrid), as well as observations of s-process enhancement in AGB stars, planetary nebulae, and barium stars. This comparison shows that the relatively high abundance of12C in the He-rich intershell in aton results in as-process abundance pattern that favours the second over the first s-process peak for all the masses explored. Also, our choice of an advective as opposed to diffusive numerical scheme for the convective overshoot results in significants-process nucleosynthesis also for the 5 M models, which may be in contradiction with observations. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Comments: 19 pages, 11 figures, 2 appendices

Journal ref: A&A 657, A28 (2022)

TOI-532b: The Habitable-zone Planet Finder confirms a Large Super Neptune in the Neptune Desert orbiting a metal-rich M dwarf host

Authors: Shubham Kanodia, Gudmundur Stefansson, Caleb I. Canas, Marissa Maney, Andrea S. Lin, Joe P. Ninan, Sinclaire Jones, Andrew J. Monson, Brock A. Parker, Henry A. Kobulnicky, Jason Rothenberg, Corey Beard, Jack Lubin, Paul Robertson, Arvind F. Gupta, Suvrath Mahadevan, William D. Cochran, Chad F. Bender, Scott A. Diddams, Connor Fredrick, Samuel P. Halverson, Suzanne L. Hawley, Fred R. Hearty, Leslie Hebb, Ravi K. Kopparapu , et al. (8 additional authors not shown)

Abstract: We confirm the planetary nature of TOI-532b, using a combination of precise near-infrared radial velocities with the Habitable-zone Planet Finder, TESS light curves, ground based photometric follow-up, and high-contrast imaging. TOI-532 is a faint (J$\sim 11.5$) metal-rich M dwarf with Teff = $3957\pm69$ K and [Fe/H] = $0.38\pm0.04$; it hosts a transiting gaseous planet with a period of… ▽ More We confirm the planetary nature of TOI-532b, using a combination of precise near-infrared radial velocities with the Habitable-zone Planet Finder, TESS light curves, ground based photometric follow-up, and high-contrast imaging. TOI-532 is a faint (J$\sim 11.5$) metal-rich M dwarf with Teff = $3957\pm69$ K and [Fe/H] = $0.38\pm0.04$; it hosts a transiting gaseous planet with a period of $\sim 2.3$ days. Joint fitting of the radial velocities with the TESS and ground-based transits reveal a planet with radius of $5.82\pm0.19$ R$_{\oplus}$, and a mass of $61.5_{-9.3}^{+9.7}$ M$_{\oplus}$. TOI-532b is the largest and most massive super Neptune detected around an M dwarf with both mass and radius measurements, and it bridges the gap between the Neptune-sized planets and the heavier Jovian planets known to orbit M dwarfs. It also follows the previously noted trend between gas giants and host star metallicity for M dwarf planets. In addition, it is situated at the edge of the Neptune desert in the Radius--Insolation plane, helping place constraints on the mechanisms responsible for sculpting this region of planetary parameter space. △ Less

Submitted 9 September, 2021; v1 submitted 28 July, 2021; originally announced July 2021.

Comments: 19 pages, 9 figures, 4 tables. arXiv admin note: text overlap with arXiv:2006.14546

Journal ref: Astronomical Journal, 162, 135, 2021

Searching for solar KDAR with DUNE

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, D. Adams, M. Adinolfi, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Ali-Mohammadzadeh, T. Alion, K. Allison, S. Alonso Monsalve, M. Alrashed, C. Alt, A. Alton, P. Amedo, J. Anderson, C. Andreopoulos, M. Andreotti, M. P. Andrews , et al. (1157 additional authors not shown)

Abstract: The observation of 236 MeV muon neutrinos from kaon-decay-at-rest (KDAR) originating in the core of the Sun would provide a unique signature of dark matter annihilation. Since excellent angle and energy reconstruction are necessary to detect this monoenergetic, directional neutrino flux, DUNE with its vast volume and reconstruction capabilities, is a promising candidate for a KDAR neutrino search.… ▽ More The observation of 236 MeV muon neutrinos from kaon-decay-at-rest (KDAR) originating in the core of the Sun would provide a unique signature of dark matter annihilation. Since excellent angle and energy reconstruction are necessary to detect this monoenergetic, directional neutrino flux, DUNE with its vast volume and reconstruction capabilities, is a promising candidate for a KDAR neutrino search. In this work, we evaluate the proposed KDAR neutrino search strategies by realistically modeling both neutrino-nucleus interactions and the response of DUNE. We find that, although reconstruction of the neutrino energy and direction is difficult with current techniques in the relevant energy range, the superb energy resolution, angular resolution, and particle identification offered by DUNE can still permit great signal/background discrimination. Moreover, there are non-standard scenarios in which searches at DUNE for KDAR in the Sun can probe dark matter interactions. △ Less

Submitted 26 October, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 19 pages, 13 figures

Report number: FERMILAB-PUB-21-322-LBNF-ND

Journal ref: JCAP10(2021)065

Synthetic observables for electron-capture supernovae and low-mass core collapse supernovae

Authors: Alexandra Kozyreva, Petr Baklanov, Samuel Jones, Georg Stockinger, Hans-Thomas Janka

Abstract: Stars in the mass range from 8 to 10 solar masses are expected to produce one of two types of supernovae (SNe), either electron-capture supernovae (ECSNe) or core-collapse supernovae (CCSNe), depending on their previous evolution. Either of the associated progenitors retain extended and massive hydrogen-rich envelopes, the observables of these SNe are, therefore, expected to be similar. In this st… ▽ More Stars in the mass range from 8 to 10 solar masses are expected to produce one of two types of supernovae (SNe), either electron-capture supernovae (ECSNe) or core-collapse supernovae (CCSNe), depending on their previous evolution. Either of the associated progenitors retain extended and massive hydrogen-rich envelopes, the observables of these SNe are, therefore, expected to be similar. In this study we explore the differences in these two types of SNe. Specifically, we investigate three different progenitor models: a solar-metallicity ECSN progenitor with an initial mass of 8.8 solar masses, a zero-metallicity progenitor with 9.6 solar masses, and a solar-metallicity progenitor with 9 solar masses, carrying out radiative transfer simulations for these progenitors. We present the resulting light curves for these models. The models exhibit very low photospheric velocity variations of about 2000 km/s, therefore, this may serve as a convenient indicator of low-mass SNe. The ECSN has very unique light curves in broad bands, especially the U band, and does not resemble any currently observed SN. This ECSN progenitor being part of a binary will lose its envelope for which reason the light curve becomes short and undetectable. The SN from the 9.6 solar masses progenitor exhibits also quite an unusual light curve, explained by the absence of metals in the initial composition. The artificially iron polluted 9.6 solar masses model demonstrates light curves closer to normal SNe IIP. The SN from the 9 solar masses progenitor remains the best candidate for so-called low-luminosity SNe IIP like SN 1999br and SN 2005cs. △ Less

Submitted 4 February, 2021; originally announced February 2021.

Comments: 19 pages, 19 figures. Accepted by MNRAS

Supernova Neutrino Burst Detection with the Deep Underground Neutrino Experiment

Authors: DUNE collaboration, B. Abi, R. Acciarri, M. A. Acero, G. Adamov, D. Adams, M. Adinolfi, Z. Ahmad, J. Ahmed, T. Alion, S. Alonso Monsalve, C. Alt, J. Anderson, C. Andreopoulos, M. P. Andrews, F. Andrianala, S. Andringa, A. Ankowski, M. Antonova, S. Antusch, A. Aranda-Fernandez, A. Ariga, L. O. Arnold, M. A. Arroyave, J. Asaadi , et al. (949 additional authors not shown)

Abstract: The Deep Underground Neutrino Experiment (DUNE), a 40-kton underground liquid argon time projection chamber experiment, will be sensitive to the electron-neutrino flavor component of the burst of neutrinos expected from the next Galactic core-collapse supernova. Such an observation will bring unique insight into the astrophysics of core collapse as well as into the properties of neutrinos. The gen… ▽ More The Deep Underground Neutrino Experiment (DUNE), a 40-kton underground liquid argon time projection chamber experiment, will be sensitive to the electron-neutrino flavor component of the burst of neutrinos expected from the next Galactic core-collapse supernova. Such an observation will bring unique insight into the astrophysics of core collapse as well as into the properties of neutrinos. The general capabilities of DUNE for neutrino detection in the relevant few- to few-tens-of-MeV neutrino energy range will be described. As an example, DUNE's ability to constrain the $ν_e$ spectral parameters of the neutrino burst will be considered. △ Less

Submitted 29 May, 2021; v1 submitted 15 August, 2020; originally announced August 2020.

Comments: 29 pages, 17 figures; paper based on DUNE Technical Design Report. arXiv admin note: substantial text overlap with arXiv:2002.03005

Report number: FERMILAB-PUB-20-380-LBNF

Cold Molecular Gas and Free-Free Emission from Hot, Dust-Obscured Galaxies at $z$~3

Authors: J. I. Penney, A. W. Blain, R. J. Assef, T. Diaz-Santos, J. J. González-López, C. -W. Tsai, M. Aravena, P. R. M. Eisenhardt, S. F. Jones, H. D. Jun, M. Kim, D. Stern, J. Wu

Abstract: We report on observations of redshifted CO(1-0) line emission and observed-frame $\rm\sim$ 30GHz radio continuum emission from five ultra-luminous, mid-IR selected hot, Dust-Obscured Galaxies (Hot DOGs) at $z\rm\gtrsim$ 3 using the Karl G. Jansky Very Large Array. We detect CO(1-0) line emission in all five Hot DOGs, with one of them at high signal to noise. We analyse FIR-radio spectral energy di… ▽ More We report on observations of redshifted CO(1-0) line emission and observed-frame $\rm\sim$ 30GHz radio continuum emission from five ultra-luminous, mid-IR selected hot, Dust-Obscured Galaxies (Hot DOGs) at $z\rm\gtrsim$ 3 using the Karl G. Jansky Very Large Array. We detect CO(1-0) line emission in all five Hot DOGs, with one of them at high signal to noise. We analyse FIR-radio spectral energy distributions, including dust, free-free and synchrotron emission for the galaxies. We find that most of the 115 GHz rest-frame continuum is mostly due to synchrotron or free-free emission, with only a potentially small contribution from thermal emission. We see a deficit in the rest-frame 115 GHz continuum emission compared to dusty star-forming galaxies (DSFGs) and sub-millimetre galaxies (SMGs) at high redshift, suggesting that Hot DOGs do not have similar cold gas reserves compared with star-forming galaxies. One target, W2305-0039, is detected in the FIRST 1.4 GHz survey, and is likely to possess compact radio jets. We compare to the FIR-radio correlation, and find that at least half of the Hot DOGs in our sample are radio-quiet with respect to normal galaxies. These findings suggest that Hot DOGs have comparably less cold molecular gas than star-forming galaxies at lower, $z\rm\sim$ 2 redshifts, and are dominated by powerful, yet radio-quiet AGN. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: 15 pages, 12 figures, accepted for publication in MNRAS

The Nucleosynthetic Yields of Core-Collapse Supernovae, prospects for the Next Generation of Gamma-Ray Astronomy

Authors: Sydney Andrews, Chris L. Fryer, Samuel W. Jones, Wesley P. Even, Marco Pignatari

Abstract: Though the neutrino-driven convection model for the core-collapse explosion mechanism has received strong support in recent years, there are still many uncertainties in the explosion parameters -- such as explosion energy, remnant mass, and end-of-life stellar abundances as initial conditions. Using a broad set of spherically symmetric core-collapse simulations we examine the effects of these key… ▽ More Though the neutrino-driven convection model for the core-collapse explosion mechanism has received strong support in recent years, there are still many uncertainties in the explosion parameters -- such as explosion energy, remnant mass, and end-of-life stellar abundances as initial conditions. Using a broad set of spherically symmetric core-collapse simulations we examine the effects of these key parameters on explosive nucleosynthesis and final explosion yields. Post-bounce temperature and density evolution of ZAMS 15, 20, and 25 solar mass progenitors are post-processed through the Nucleosynthesis Grid (NuGrid) nuclear network to obtain detailed explosive yields. In particular, this study focuses on radio-isotopes that are of particular interest to the next generation of gamma-ray astronomical observations; 43K, 47Ca, 44Sc, 47Sc, 48V, 48Cr, 51Cr, 52Mn, 59Fe, 56Co, 57Co, 57Ni. These nuclides may be key in advancing our understanding of the inner workings of core-collapse supernovae by probing the parameters of the explosion engine. We find that the isotopes that are strong indicators of explosion energy are 43K, 47Ca, 44Sc, 47Sc, and 59Fe, those that are dependent on the progenitor structure are 48V, 51Cr, and 57Co, and those that probe neither are 48Cr, 52Mn, 57Ni, and 56Co. We discuss prospects of observing these radionuclides in supernova remnants. △ Less

Submitted 22 December, 2019; originally announced December 2019.

Comments: 17 pages, 19 figures, submitted to ApJ

A New Model for Electron-Capture Supernovae in Galactic Chemical Evolution

Authors: S. Jones, B. Côté, F. K. Roepke, S. Wanajo

Abstract: We examine the contribution of electron-capture supernovae (ECSNe), low-mass SNe from collapsing Fe cores (FeCCSNe), and rotating massive stars to the chemical composition of the Galaxy. Our model includes contributions to chemical evolution from both thermonuclear ECSNe (tECSNe) and gravitational collapse ECSNe (cECSNe). We show that if ECSNe are predominantly gravitational collapse SNe but about… ▽ More We examine the contribution of electron-capture supernovae (ECSNe), low-mass SNe from collapsing Fe cores (FeCCSNe), and rotating massive stars to the chemical composition of the Galaxy. Our model includes contributions to chemical evolution from both thermonuclear ECSNe (tECSNe) and gravitational collapse ECSNe (cECSNe). We show that if ECSNe are predominantly gravitational collapse SNe but about 15% are partial thermonuclear explosions, the model is able to reproduce the solar abundances of several important and problematic isotopes including $^{48}$Ca, $^{50}$Ti and $^{54}$Cr together with $^{58}$Fe, $^{64}$Ni, $^{82}$Se and $^{86}$Kr and several of the Zn--Zr isotopes. A model in which no cECSNe occur, only tECSNe with low-mass FeCCSNe or rotating massive stars, proves also very successful at reproducing the solar abundances for these isotopes. Despite the small mass range for the progenitors of ECSNe and low-mass FeCCSNe, the large production factors suffice for the solar inventory of the above isotopes. Our model is compelling because it introduces no new tensions with the solar abundance distribution for a Milky Way model -- only tending to improve the model predictions for several isotopes. The proposed astrophysical production model thus provides a natural and elegant way to explain one of the last uncharted territories on the periodic table of astrophysical element production. △ Less

Submitted 6 August, 2019; originally announced August 2019.

Comments: 10 pages, 4 figures; accepted for publication in the Astrophysical Journal

All-sky Medium Energy Gamma-ray Observatory: Exploring the Extreme Multimessenger Universe

Authors: Julie McEnery, Juan Abel Barrio, Ivan Agudo, Marco Ajello, José-Manuel Álvarez, Stefano Ansoldi, Sonia Anton, Natalia Auricchio, John B. Stephen, Luca Baldini, Cosimo Bambi, Matthew Baring, Ulisses Barres, Denis Bastieri, John Beacom, Volker Beckmann, Wlodek Bednarek, Denis Bernard, Elisabetta Bissaldi, Peter Bloser, Harsha Blumer, Markus Boettcher, Steven Boggs, Aleksey Bolotnikov, Eugenio Bottacini , et al. (160 additional authors not shown)

Abstract: The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger… ▽ More The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band. △ Less

Submitted 25 November, 2019; v1 submitted 17 July, 2019; originally announced July 2019.

Comments: Astro2020 APC White Paper Updated to make small change to author list in metadata

Chromium Nucleosynthesis and Silicon-Carbon Shell Mergers in Massive Stars

Authors: Benoit Côté, Samuel Jones, Falk Herwig, Marco Pignatari

Abstract: We analyze the production of the element Cr in galactic chemical evolution (GCE) models using the NuGrid nucleosynthesis yields set. We show that the unusually large [Cr/Fe] abundance at [Fe/H] $\approx 0$ reported by previous studies using those yields and predicted by our Milky Way model originates from the merging of convective Si-burning and C-burning shells in a 20 $M_\odot$ model at metallic… ▽ More We analyze the production of the element Cr in galactic chemical evolution (GCE) models using the NuGrid nucleosynthesis yields set. We show that the unusually large [Cr/Fe] abundance at [Fe/H] $\approx 0$ reported by previous studies using those yields and predicted by our Milky Way model originates from the merging of convective Si-burning and C-burning shells in a 20 $M_\odot$ model at metallicity $Z=0.01$, about an hour before the star explodes. This merger mixes the incomplete burning material in the Si shell, including $^{51}$V and $^{52}$Cr, out to the edge of the carbon/oxygen (CO) core. The adopted supernova model ejects the outer 2 $M_\odot$ of the CO core, which includes a significant fraction of the Cr-rich material. When including this 20 $M_\odot$ model at $Z=0.01$ in the yields interpolation scheme of our GCE model for stars in between 15 and 25 $M_\odot$, we overestimate [Cr/Fe] by an order of magnitude at [Fe/H] $\approx$ 0 relative to observations in the Galactic disk. This raises a number of questions regarding the occurrence of Si-C shell mergers in nature, the accuracy of different simulation approaches, and the impact of such mergers on the pre-supernova structure and explosion dynamics. According to the conditions in this 1D stellar model, the substantial penetration of C-shell material into the Si-shell could launch a convective-reactive global oscillation, if a merger does take place. In any case, GCE provides stringent constraints on the outcome of this stellar evolution phase. △ Less

Submitted 17 June, 2019; originally announced June 2019.

Comments: 8 pages, 5 figures, submitted to ApJL

Discovery of an Exceptionally Strong $β$-Decay Transition of $^{20}$F and Implications for the Fate of Intermediate-Mass Stars

Authors: O. S. Kirsebom, S. Jones, D. F. Strömberg, G. Martínez-Pinedo, K. Langanke, F. K. Roepke, B. A. Brown, T. Eronen, H. O. U. Fynbo, M. Hukkanen, A. Idini, A. Jokinen, A. Kankainen, J. Kostensalo, I. Moore, H. Möller, S. T. Ohlmann, H. Penttilä, K. Riisager, S. Rinta-Antila, P. C. Srivastava, J. Suhonen, W. H. Trzaska, J. Äystö

Abstract: A significant fraction of stars between 7-11 solar masses are thought to become supernovae, but the explosion mechanism is unclear. The answer depends critically on the rate of electron capture on $^{20}$Ne in the degenerate oxygen-neon stellar core. However, due to the unknown strength of the transition between the ground states of $^{20}$Ne and $^{20}$F, it has not previously been possible to fu… ▽ More A significant fraction of stars between 7-11 solar masses are thought to become supernovae, but the explosion mechanism is unclear. The answer depends critically on the rate of electron capture on $^{20}$Ne in the degenerate oxygen-neon stellar core. However, due to the unknown strength of the transition between the ground states of $^{20}$Ne and $^{20}$F, it has not previously been possible to fully constrain the rate. By measuring the transition, we have established that its strength is exceptionally large and enhances the capture rate by several orders of magnitude. This has a decisive impact on the evolution of the core, increasing the likelihood that the star is (partially) disrupted by a thermonuclear explosion rather than collapsing to form a neutron star. Importantly, our measurement resolves the last remaining nuclear physics uncertainty in the final evolution of degenerate oxygen-neon stellar cores, allowing future studies to address the critical role of convection, which at present is poorly understood. △ Less

Submitted 2 November, 2019; v1 submitted 22 May, 2019; originally announced May 2019.

Comments: 6 pages, 5 figures, supplemental material; submitted to Science on January 26, 2019; revised version submitted to Physical Review Letters on May 22, 2019; re-submitted to Physical Review Letters on August 15, 2019, after minor revisions; revised title

Journal ref: Phys. Rev. Lett. 123, 262701 (2019)

$^{60}$Fe in core-collapse supernovae and prospects for X-ray and $γ$-ray detection in supernova remnants

Authors: Samuel Jones, Heiko Moeller, Chris L. Fryer, Christopher J. Fontes, Reto Trappitsch, Wesley P. Even, Aaron Couture, Matthew R. Mumpower, Samar Safi-Harb

Abstract: We investigate $^{60}$Fe in massive stars and core-collapse supernovae focussing on uncertainties that influence its production in 15, 20 and 25 $M_\odot$ stars at solar metallicity. We find that the $^{60}$Fe yield is a monotonic increasing function of the uncertain $^{59}$Fe$(n,γ)^{60}$Fe cross section and that a factor of 10 reduction in the reaction rate results in a factor 8-10 reduction in t… ▽ More We investigate $^{60}$Fe in massive stars and core-collapse supernovae focussing on uncertainties that influence its production in 15, 20 and 25 $M_\odot$ stars at solar metallicity. We find that the $^{60}$Fe yield is a monotonic increasing function of the uncertain $^{59}$Fe$(n,γ)^{60}$Fe cross section and that a factor of 10 reduction in the reaction rate results in a factor 8-10 reduction in the $^{60}$Fe yield; while a factor of 10 increase in the rate increases the yield by a factor 4-7. We find that none of the 189 simulations we have performed are consistent with a core-collapse supernova triggering the formation of the Solar System, and that only models using $^{59}$Fe$(n,γ)^{60}$Fe cross section that is less than or equal to that from NON-SMOKER can reproduce the observed $^{60}$Fe/$^{26}$Al line flux ratio in the diffuse ISM. We examine the prospects of detecting old core-collapse supernova remnants (SNRs) in the Milky Way from their $γ$-ray emission from the decay of $^{60}$Fe, finding that the next generation of gamma-ray missions could be able to discover up to $\sim100$ such old SNRs as well as measure the $^{60}$Fe yields of a handful of known Galactic SNRs. We also predict the X-ray spectrum that is produced by atomic transitions in $^{60}$Co following its ionization by internal conversion and give theoretical X-ray line fluxes as a function of remnant age as well as the Doppler and fine-structure line broadening effects. The X-ray emission presents an interesting prospect for addressing the missing SNR problem with future X-ray missions. △ Less

Submitted 15 February, 2019; originally announced February 2019.

Comments: 25 pages, 25 figures, 1 table. Accepted for publication in MNRAS

Catching Element Formation In The Act

Authors: Chris L. Fryer, Frank Timmes, Aimee L. Hungerford, Aaron Couture, Fred Adams, Wako Aoki, Almudena Arcones, David Arnett, Katie Auchettl, Melina Avila, Carles Badenes, Eddie Baron, Andreas Bauswein, John Beacom, Jeff Blackmon, Stephane Blondin, Peter Bloser, Steve Boggs, Alan Boss, Terri Brandt, Eduardo Bravo, Ed Brown, Peter Brown, Steve Bruenn. Carl Budtz-Jorgensen, Eric Burns , et al. (194 additional authors not shown)

Abstract: Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-ray… ▽ More Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions. △ Less

Submitted 7 February, 2019; originally announced February 2019.

Comments: 14 pages including 3 figures

Report number: LA-UR-18-29748

Nucleosynthetic Yields from Neutron Stars Accreting in Binary Common Envelopes

Authors: James D. Keegans, Chris L. Fryer, Samuel W. Jones, Benoit Cote, Krzysztof Belczynski, Falk Herwig, Marco Pignatari, Alison M. Laird, Christian Aa. Diget

Abstract: Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario f… ▽ More Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario for neutron star binaries, the system might pass through a phase where a neutron star spirals into the envelope of its giant star companion. These phases lead to mass accretion onto the neutron star. Accretion onto these common-envelope-phase neutron stars can eject matter that has undergone burning near to the neutron star surface. This paper presents nucleosynthetic yields of this ejected matter, using population synthesis models to study the importance of these nucleosynthetic yields in a galactic chemical evolution context. Depending on the extreme conditions in temperature and density found in the accreted material, both proton-rich and neutron-rich nucleosynthesis can be obtained, with efficient production of neutron rich isotopes of low Z material at the most extreme conditions, and proton rich isotopes, again at low Z, in lower density models. Final yields are found to be extremely sensitive to the physical modeling of the accretion phase. We show that neutron stars accreting in binary common envelopes might be a new relevant site for galactic chemical evolution, and therefore more comprehensive studies are needed to better constrain nucleosynthesis in these objects. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Comments: 21 pages, 14 figures

Report number: LA-UR-18-29776

Remnants and ejecta of thermonuclear electron-capture supernovae: Constraining oxygen-neon deflagrations in high-density white dwarfs

Authors: S. Jones, F. K. Roepke, C. Fryer, A. J. Ruiter, I. R. Seitenzahl, L. R. Nittler, S. T. Ohlmann, R. Reifarth, M. Pignatari, K. Belczynski

Abstract: (Abridged) The explosion mechanism of electron-capture supernovae (ECSNe) remains equivocal. We attempt to constrain the explosion mechanism (neutron-star-forming implosion or thermonuclear explosion) and the frequency of occurrence of ECSNe using nucleosynthesis simulations of the latter scenario, population synthesis, the solar abundance distribution, pre-solar meteoritic oxide grain isotopic ra… ▽ More (Abridged) The explosion mechanism of electron-capture supernovae (ECSNe) remains equivocal. We attempt to constrain the explosion mechanism (neutron-star-forming implosion or thermonuclear explosion) and the frequency of occurrence of ECSNe using nucleosynthesis simulations of the latter scenario, population synthesis, the solar abundance distribution, pre-solar meteoritic oxide grain isotopic ratio measurements and the white dwarf mass-radius relation. Tracer particles from 3d hydrodynamic simulations were post-processed with a large nuclear reaction network in order to determine the complete compositional state of the bound ONeFe remnant and the ejecta, and population synthesis simulations were performed in order to estimate the ECSN rate with respect to the CCSN rate. The 3d deflagration simulations drastically overproduce the neutron-rich isotopes $^{48}$Ca, $^{50}$Ti, $^{54}$Cr, $^{60}$Fe and several of the Zn isotopes relative to their solar abundances. Using the solar abundance distribution as our constraint, we place an upper limit on the frequency of thermonuclear ECSNe as 1$-$3~\% the frequency at which core-collapse supernovae (FeCCSNe) occur. This is on par with or 1~dex lower than the estimates for ECSNe from single stars. The upper limit from the yields is also in relatively good agreement with the predictions from our population synthesis simulations. The $^{54}$Cr/$^{52}$Cr and $^{50}$Ti/$^{48}$Ti isotopic ratios in the ejecta are a near-perfect match with recent measurements of extreme pre-solar meteoritc oxide grains, and $^{53}$Cr/$^{52}$Cr can also be matched if the ejecta condenses before mixing with the interstellar medium. Theoretical mass-radius relations for the bound ONeFe WD remnants of these explosions are apparently consistent with several observational WD candidates. △ Less

Submitted 4 January, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: 23 pages, 13 figures, 7 tables; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 622, A74 (2019)

The Environments of Luminous Radio - WISE Selected Infrared Galaxies

Authors: J. I. Penney, A. W. Blain, D. Wylezalek, N. A. Hatch, C. Lonsdale, A. Kimball, R. J. Assef, J. J. Condon, P. R. M. Eisenhardt, S. F. Jones, M. Kim, M. Lacy, S. I. Muldrew, S. Petty, A. Sajina, A. Silva, D. Stern, T. Diaz-Santos, C-W. Tsai, J. Wu

Abstract: We have observed the environments of a population of 33 heavily dust obscured, ultra-luminous, high-redshift galaxies, selected using WISE and NVSS at $z>$1.3 with the Infra-Red Array Camera on the $Spitzer$ Space Telescope over $\rm5.12\,'\times5.12\,'$ fields. Colour selections are used to quantify any potential overdensities of companion galaxies in these fields. We find no significant excess o… ▽ More We have observed the environments of a population of 33 heavily dust obscured, ultra-luminous, high-redshift galaxies, selected using WISE and NVSS at $z>$1.3 with the Infra-Red Array Camera on the $Spitzer$ Space Telescope over $\rm5.12\,'\times5.12\,'$ fields. Colour selections are used to quantify any potential overdensities of companion galaxies in these fields. We find no significant excess of galaxies with the standard colour selection for IRAC colours of $\rm[3.6]-[4.5]>-0.1$ consistent with galaxies at $z>$1.3 across the whole fields with respect to wide-area $Spitzer$ comparison fields, but there is a $\rm>2σ$ statistical excess within $\rm0.25\,'$ of the central radio-WISE galaxy. Using a colour selection of $\rm[3.6]-[4.5]>0.4$, 0.5 magnitudes redder than the standard method of selecting galaxies at $z>$1.3, we find a significant overdensity, in which $\rm76\%$ ($\rm33\%$) of the 33 fields have a surface density greater than the $\rm3σ$ ($\rm5σ$) level. There is a statistical excess of these redder galaxies within $\rm0.5\,'$, rising to a central peak $\rm\sim2$--4 times the average density. This implies that these galaxies are statistically linked to the radio-WISE selected galaxy, indicating similar structures to those traced by red galaxies around radio-loud AGN. △ Less

Submitted 16 November, 2018; originally announced November 2018.

Comments: 17 pages, 16 figures, 2 tables

Super-Eddington Accretion in the WISE-Selected Extremely Luminous Infrared Galaxy W2246-0526

Authors: Chao-Wei Tsai, Peter Eisenhardt, Hyunsung Jun, Jingwen Wu, Roberto Assef, Andrew Blain, Tanio Diaz-Santos, Suzy Jones, Daniel Stern, Edward Wright, Sherry Yeh

Abstract: We use optical and near-infrared spectroscopy to observe rest-UV emission lines and estimate the black hole mass of WISEA J224607.56-052634.9 (W2246-0526) at z = 4.601, the most luminous hot dust-obscured galaxy yet discovered by WISE. From the broad component of the MgII-2799A emission line, we measure a black hole mass of log (M_BH/M_sun) = 9.6 +- 0.4. The broad CIV-1549A line is asymmetric and… ▽ More We use optical and near-infrared spectroscopy to observe rest-UV emission lines and estimate the black hole mass of WISEA J224607.56-052634.9 (W2246-0526) at z = 4.601, the most luminous hot dust-obscured galaxy yet discovered by WISE. From the broad component of the MgII-2799A emission line, we measure a black hole mass of log (M_BH/M_sun) = 9.6 +- 0.4. The broad CIV-1549A line is asymmetric and significantly blue-shifted. The derived M_BH from the blueshift-corrected broad CIV line width agrees with the MgII result. From direct measurement using a well-sampled SED, the bolometric luminosity is 3.6 * 10^14 L_sun. The corresponding Eddington ratio for W2246-0526 is lambda_Edd = L_AGN / L_Edd = 2.8. This high Eddington ratio may reach the level where the luminosity is saturating due to photon trapping in the accretion flow, and be insensitive to the mass accretion rate. In this case, the M_BH growth rate in W2246-0526 would exceed the apparent accretion rate derived from the observed luminosity. △ Less

Submitted 6 October, 2018; originally announced October 2018.

Comments: 9 pages in emulateapj format, including 5 figures and 1 table. Accepted by the Astrophysical Journal

The i-process yields of rapidly-accreting white dwarfs from multicycle He-shell flash stellar evolution models with mixing parameterizations from 3D hydrodynamics simulations

Authors: Pavel Denissenkov, Falk Herwig, Paul Woodward, Robert Andrassy, Marco Pignatari, Samuel Jones

Abstract: We have modelled the multicycle evolution of rapidly-accreting CO white dwarfs (RAWDs) with stable H burning intermittent with strong He-shell flashes on their surfaces for $0.7\leq M_\mathrm{RAWD}/M_\odot\leq 0.75$ and [Fe/H] ranging from $0$ to $-2.6$. We have also computed the i-process nucleosynthesis yields for these models. The i process occurs when convection driven by the He-shell flash in… ▽ More We have modelled the multicycle evolution of rapidly-accreting CO white dwarfs (RAWDs) with stable H burning intermittent with strong He-shell flashes on their surfaces for $0.7\leq M_\mathrm{RAWD}/M_\odot\leq 0.75$ and [Fe/H] ranging from $0$ to $-2.6$. We have also computed the i-process nucleosynthesis yields for these models. The i process occurs when convection driven by the He-shell flash ingests protons from the accreted H-rich surface layer, which results in maximum neutron densities $N_\mathrm{n,max}\approx 10^{13}$-$10^{15}\ \mathrm{cm}^{-3}$. The H-ingestion rate and the convective boundary mixing (CBM) parameter $f_\mathrm{top}$ adopted in the one-dimensional nucleosynthesis and stellar evolution models are constrained through 3D hydrodynamic simulations. The mass ingestion rate and, for the first time, the scaling laws for the CBM parameter $f_\mathrm{top}$ have been determined from 3D hydrodynamic simulations. We confirm our previous result that the high-metallicity RAWDs have a low mass retention efficiency ($η< 10\%$). A new result is that RAWDs with [Fe/H]$< -2$ have $η> 20\%$, therefore their masses may reach the Chandrasekhar limit and they may eventually explode as SNeIa. This result and the good fits of the i-process yields from the metal-poor RAWDs to the observed chemical composition of the CEMP-r/s stars suggest that some of the present-day CEMP-r/s stars could be former distant members of triple systems, orbiting close binary systems with RAWDs that may have later exploded as SNeIa. △ Less

Submitted 9 July, 2019; v1 submitted 10 September, 2018; originally announced September 2018.

Comments: 15 pages, 12 figures, 2 tables, accepted for publication in MNRAS

Cyberhubs: Virtual Research Environments for Astronomy

Authors: Falk Herwig, Robert Andrassy, Nic Annau, Ondrea Clarkson, Benoit Cote, Aaron D'Sa, Sam Jones, Belaid Moa, Jericho O'Connell, David Porter, Christian Ritter, Paul Woodward

Abstract: Collaborations in astronomy and astrophysics are faced with numerous cyber infrastructure challenges, such as large data sets, the need to combine heterogeneous data sets, and the challenge to effectively collaborate on those large, heterogeneous data sets with significant processing requirements and complex science software tools. The cyberhubs system is an easy-to-deploy package for small to med… ▽ More Collaborations in astronomy and astrophysics are faced with numerous cyber infrastructure challenges, such as large data sets, the need to combine heterogeneous data sets, and the challenge to effectively collaborate on those large, heterogeneous data sets with significant processing requirements and complex science software tools. The cyberhubs system is an easy-to-deploy package for small to medium-sized collaborations based on the Jupyter and Docker technology, that allows web-browser enabled, remote, interactive analytic access to shared data. It offers an initial step to address these challenges. The features and deployment steps of the system are described, as well as the requirements collection through an account of the different approaches to data structuring, handling and available analytic tools for the NuGrid and PPMstar collaborations. NuGrid is an international collaboration that creates stellar evolution and explosion physics and nucleosynthesis simulation data. The PPMstar collaboration performs large-scale 3D stellar hydrodynamics simulation of interior convection in the late phases of stellar evolution. Examples of science that is presently performed on cyberhubs, in the areas 3D stellar hydrodynamic simulations, stellar evolution and nucleosynthesis and Galactic chemical evolution, are presented. △ Less

Submitted 6 February, 2018; originally announced February 2018.

Comments: 19 pages, 9 figures, ApJ Supplement Special Issue on Data, accepted

Convective boundary mixing in a post-He core burning massive star model

Authors: Austin Davis, Sam Jones, Falk Herwig

Abstract: Convective boundary mixing (CBM) in the advanced evolutionary stages of massive stars is not well understood. Structural changes caused by convection have an impact on the evolution as well as the subsequent supernova, or lack thereof. The effects of convectively driven mixing across convective boundaries during the post He core burning evolution of $25\mathrm{M}_{\odot}$, solar-metallicity, non-r… ▽ More Convective boundary mixing (CBM) in the advanced evolutionary stages of massive stars is not well understood. Structural changes caused by convection have an impact on the evolution as well as the subsequent supernova, or lack thereof. The effects of convectively driven mixing across convective boundaries during the post He core burning evolution of $25\mathrm{M}_{\odot}$, solar-metallicity, non-rotating stellar models is studied using the MESA stellar evolution code. CBM is modelled using the exponentially decaying diffusion coefficient equation, the free parameter of which, $f_{\mathrm{cbm}}$, is varied systematically throughout the course of the stellar model's evolution with values of $(0.002, 0.012, 0.022, 0.032)$. The effect of varying this parameter produces mass ranges at collapse in the ONe, Si, Fe cores of $(1.82\mathrm{M}_{\odot}, 4.36\mathrm{M}_{\odot})$, $(1.67\mathrm{M}_{\odot}, 1.99\mathrm{M}_{\odot})$ and $(1.46\mathrm{M}_{\odot}, 1.70\mathrm{M}_{\odot})$ respectively, with percent differences from the model with minimal CBM as large as 86.3%. At the presupernova stage, the compactness of the stellar cores from O'Connor & Ott (2011), $ξ_M$, exhibit a range of $(0.120, 0.354)$, suggesting that the extent of CBM in the advanced burning stages of massive stars is an important consideration for the explodability and type of compact remnant. The nucleosynthetic yields from the models, most notably C, O, Ne, Mg and Si are also significantly affected by the CBM assumptions, showing non-linear trends with increased mixing. The simulations show that interactions between convective C, Ne and O shells produce significant non-linear changes in the evolution, whereas from the end of Si burning, the structural changes attributed to the CBM are dominated by the growth of the convective C shell. Progenitor structures for all the models are available from HERE (link and DOI to appear). △ Less

Submitted 30 November, 2017; originally announced December 2017.

Comments: 19 pages, 19 figures

Report number: LA-UR-17-30635

NuGrid Stellar Data Set. II. Stellar Yields from H to Bi for Stellar Models with Mzams = 1 to 25Msun and Z = 0.0001 to 0.02

Authors: C. Ritter, F. Herwig, S. Jones, M. Pignatari, C. Fryer, R. Hirschi

Abstract: We provide here a significant extension of the NuGrid Set 1 models in mass coverage and toward lower metallicity, adopting the same physics assumptions. The combined data set now includes the initial masses M/Msun = 1, 1.65, 2, 3, 4, 5, 6, 7, 12, 15, 20, 25 for Z = 0.02, 0.01, 0.006, 0.001, 0.0001 with alpha-enhanced composition for the lowest three metallicities. These models are computed with th… ▽ More We provide here a significant extension of the NuGrid Set 1 models in mass coverage and toward lower metallicity, adopting the same physics assumptions. The combined data set now includes the initial masses M/Msun = 1, 1.65, 2, 3, 4, 5, 6, 7, 12, 15, 20, 25 for Z = 0.02, 0.01, 0.006, 0.001, 0.0001 with alpha-enhanced composition for the lowest three metallicities. These models are computed with the MESA stellar evolution code and are evolved up to the AGB, the white dwarf stage, or until core collapse. The nucleosynthesis was calculated for all isotopes in post-processing with the NuGrid mppnp code. Explosive nucleosynthesis is based on semi-analytic 1D shock models. Metallicity-dependent mass loss, convective boundary mixing in low- and intermediate mass models and H and He core burning massive star models is included. Convective O-C shell mergers in some stellar models lead to the strong production of odd-Z elements P, Cl, K and Sc. In AGB models with hot dredge-up the convective boundary mixing efficiency is reduced to accommodate for its energetic feedback. In both low-mass and massive star models at the lowest metallicity H-ingestion events are observed and lead to i-process nucleosynthesis and substantial N-15 production. Complete yield data tables, derived data products and online analytic data access are provided. △ Less

Submitted 31 August, 2018; v1 submitted 25 September, 2017; originally announced September 2017.

Comments: Published. Detailed online data available at Canadian Astronomical Data Center, DOI:10.11570/18.0002 (http://www.canfar.phys.uvic.ca/vospace/nodes/AstroDataCitationDOI/CISTI.CANFAR/18.0002/18.0002.html?view=data). See Appendix for details

Journal ref: MNRAS 480, 538, 2018

The SCUBA-2 850 $μm$ follow-up of WISE-selected, luminous dust-obscured quasars

Authors: Lulu Fan, Suzy F. Jones, Yunkun Han, Kirsten K. Knudsen

Abstract: Hot dust-obscured galaxies (Hot DOGs) are a new population recently discovered in the \wise All-Sky survey. Multiwavelength follow-up observations suggest that they are luminous, dust-obscured quasars at high redshift. Here we present the JCMT SCUBA-2 850 $μm$ follow-up observations of 10 Hot DOGs. Four out of ten Hot DOGs have been detected at $>3σ$ level. Based on the IR SED decomposition approa… ▽ More Hot dust-obscured galaxies (Hot DOGs) are a new population recently discovered in the \wise All-Sky survey. Multiwavelength follow-up observations suggest that they are luminous, dust-obscured quasars at high redshift. Here we present the JCMT SCUBA-2 850 $μm$ follow-up observations of 10 Hot DOGs. Four out of ten Hot DOGs have been detected at $>3σ$ level. Based on the IR SED decomposition approach, we derive the IR luminosities of AGN torus and cold dust components. Hot DOGs in our sample are extremely luminous with most of them having $L_{\rm IR}^{\rm tot}>10^{14} L_\odot$. The torus emissions dominate the total IR energy output. However, the cold dust contribution is still non-negligible, with the fraction of the cold dust contribution to the total IR luminosity $(\sim 8-24\%)$ being dependent on the choice of torus model. The derived cold dust temperatures in Hot DOGs are comparable to those in UV bright quasars with similar IR luminosity, but much higher than those in SMGs. Higher dust temperatures in Hot DOGs may be due to the more intense radiation field caused by intense starburst and obscured AGN activities. Fourteen and five submillimeter serendipitous sources in the 10 SCUBA-2 fields around Hot DOGs have been detected at $>3σ$ and $>3.5σ$ levels, respectively. By estimating their cumulative number counts, we confirm the previous argument that Hot DOGs lie in dense environments. Our results support the scenario in which Hot DOGs are luminous, dust-obscured quasars lying in dense environments, and being in the transition phase between extreme starburst and UV-bright quasars. △ Less

Submitted 22 September, 2017; originally announced September 2017.

Comments: 26 pages, 7 figures, PASP accepted

Image-Optimized Coronal Magnetic Field Models

Authors: Shaela I. Jones, Joseph M. Davila, Vadim M. Uritsky

Abstract: We have reported previously on a new method we are developing for using image-based information to improve global coronal magnetic field models. In that work we presented early tests of the method which proved its capability to improve global models based on flawed synoptic magnetograms, given excellent constraints on the field in the model volume. In this follow-up paper we present the results of… ▽ More We have reported previously on a new method we are developing for using image-based information to improve global coronal magnetic field models. In that work we presented early tests of the method which proved its capability to improve global models based on flawed synoptic magnetograms, given excellent constraints on the field in the model volume. In this follow-up paper we present the results of similar tests given field constraints of a nature that could realistically be obtained from quality white-light coronagraph images of the lower corona. We pay particular attention to difficulties associated with the line-of-sight projection of features outside of the assumed coronagraph image plane, and the effect on the outcome of the optimization of errors in localization of constraints. We find that substantial improvement in the model field can be achieved with this type of constraints, even when magnetic features in the images are located outside of the image plane. △ Less

Submitted 22 June, 2017; originally announced June 2017.

The evolutionary roads leading to low effective spins, high black hole masses, and O1/O2 rates of LIGO/Virgo binary black holes

Authors: K. Belczynski, J. Klencki, C. E. Fields, A. Olejak, E. Berti, G. Meynet, C. L. Fryer, D. E. Holz, R. O'Shaughnessy, D. A. Brown, T. Bulik, S. C. Leung, K. Nomoto, P. Madau, R. Hirschi, E. Kaiser, S. Jones, S. Mondal, M. Chruslinska, P. Drozda, D. Gerosa, Z. Doctor, M. Giersz, S. Ekstrom, C. Georgy , et al. (9 additional authors not shown)

Abstract: All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH… ▽ More All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org. △ Less

Submitted 6 March, 2020; v1 submitted 21 June, 2017; originally announced June 2017.

Comments: A&A accepted: revised comparison with LIGO/Virgo (41 pages)

Journal ref: A&A 636, A104 (2020)

Overdensities of SMGs around WISE-selected, ultra-luminous, high-redshift AGN

Authors: Suzy F. Jones, Andrew W. Blain, Roberto J. Assef, Peter Eisenhardt, Carol Lonsdale, James Condon, Duncan Farrah, Chao-Wei Tsai, Carrie Bridge, Jingwen Wu, Edward L. Wright, Tom Jarrett

Abstract: We investigate extremely luminous dusty galaxies in the environments around WISE-selected hot dust obscured galaxies (Hot DOGs) and WISE/radio-selected active galactic nuclei (AGNs) at average redshifts of z = 2.7 and z = 1.7, respectively. Previous observations have detected overdensities of companion submillimetre-selected sources around 10 Hot DOGs and 30 WISE/radio AGNs, with overdensities of… ▽ More We investigate extremely luminous dusty galaxies in the environments around WISE-selected hot dust obscured galaxies (Hot DOGs) and WISE/radio-selected active galactic nuclei (AGNs) at average redshifts of z = 2.7 and z = 1.7, respectively. Previous observations have detected overdensities of companion submillimetre-selected sources around 10 Hot DOGs and 30 WISE/radio AGNs, with overdensities of ~ 2 - 3 and ~ 5 - 6 , respectively. We find that the space densities in both samples to be overdense compared to normal star-forming galaxies and submillimetre galaxies (SMGs) in the SCUBA-2 Cosmology Legacy Survey (S2CLS). Both samples of companion sources have consistent mid-IR colours and mid-IR to submm ratios as SMGs. The brighter population around WISE/radio AGNs could be responsible for the higher overdensity reported. We also find the star formation rate density (SFRDs) are higher than the field, but consistent with clusters of dusty galaxies. WISE-selected AGNs appear to be good signposts for protoclusters at high redshift on arcmin scales. The results reported here provide an upper limit to the strength of angular clustering using the two-point correlation function. Monte Carlo simulations show no angular correlation, which could indicate protoclusters on scales larger than the SCUBA-2 1.5arcmin scale maps. △ Less

Submitted 10 May, 2017; originally announced May 2017.

Comments: 10 pages, 7 figures, 6 tables, MNRAS accepted

The role of the most luminous, obscured AGN in galaxy assembly at z~2

Authors: Duncan Farrah, Sara Petty, Brian Connolly, Andrew Blain, Andreas Efstathiou, Mark Lacy, Daniel Stern, Sean Lake, Tom Jarrett, Carrie Bridge, Peter Eisenhardt, Dominic Benford, Suzy Jones, Chao-Wei Tsai, Roberto Assef, Jingwen Wu, Leonidas Moustakas

Abstract: We present HST WFC3 F160W imaging and infrared spectral energy distributions for twelve extremely luminous, obscured AGN at $1.8<z<2.7$, selected via "Hot, Dust Obscured" mid-infrared colors. Their infrared luminosities span $2-15\times10^{13}$L$_{\odot}$, making them among the most luminous objects in the Universe at $z\sim2$. In all cases the infrared emission is consistent with arising at least… ▽ More We present HST WFC3 F160W imaging and infrared spectral energy distributions for twelve extremely luminous, obscured AGN at $1.8<z<2.7$, selected via "Hot, Dust Obscured" mid-infrared colors. Their infrared luminosities span $2-15\times10^{13}$L$_{\odot}$, making them among the most luminous objects in the Universe at $z\sim2$. In all cases the infrared emission is consistent with arising at least in most part from AGN activity. The AGN fractional luminosities are higher than those in either sub-millimeter galaxies, or AGN selected via other mid-infrared criteria. Adopting the $G$, M$_{20}$ and $A$ morphological parameters, together with traditional classification boundaries, infers that three quarters of the sample as mergers. Our sample do not, however, show any correlation between the considered morphological parameters and either infrared luminosity or AGN fractional luminosity. Moreover, their asymmetries and effective radii are distributed identically to those of massive galaxies at $z\sim2$. We conclude that our sample is not preferentially associated with mergers, though a significant merger fraction is still plausible. Instead, we propose that our sample are examples of the massive galaxy population at $z\sim2$ that harbor a briefly luminous, "flickering" AGN, and in which the $G$ and M$_{20}$ values have been perturbed, due to either the AGN, and/or the earliest formation stages of a bulge in an inside-out manner. Furthermore, we find that the mass assembly of the central black holes in our sample leads the mass assembly of any bulge component. Finally, we speculate that our sample represent a small fraction of the immediate antecedents of compact star-forming galaxies at $z\sim2$. △ Less

Submitted 6 June, 2017; v1 submitted 7 May, 2017; originally announced May 2017.

Comments: ApJ, accepted. Updated to reflect the accepted version

Testing a one-dimensional prescription of dynamical shear mixing with a two-dimensional hydrodynamic simulation

Authors: Philipp V. F. Edelmann, Friedrich K. Roepke, Raphael Hirschi, Cyril Georgy, Samuel Jones

Abstract: The treatment of mixing processes is still one of the major uncertainties in 1D stellar evolution models. This is mostly due to the need to parametrize and approximate aspects of hydrodynamics in hydrostatic codes. In particular, the effect of hydrodynamic instabilities in rotating stars, for example, dynamical shear instability, evades consistent description. We intend to study the accuracy of th… ▽ More The treatment of mixing processes is still one of the major uncertainties in 1D stellar evolution models. This is mostly due to the need to parametrize and approximate aspects of hydrodynamics in hydrostatic codes. In particular, the effect of hydrodynamic instabilities in rotating stars, for example, dynamical shear instability, evades consistent description. We intend to study the accuracy of the diffusion approximation to dynamical shear in hydrostatic stellar evolution models by comparing 1D models to a first-principle hydrodynamics simulation starting from the same initial conditions. We chose an initial model calculated with the stellar evolution code GENEC that is just at the onset of a dynamical shear instability but does not show any other instabilities (e.g., convection). This was mapped to the hydrodynamics code SLH to perform a 2D simulation in the equatorial plane. We compare the resulting profiles in the two codes and compute an effective diffusion coefficient for the hydro simulation. Shear instabilities develop in the 2D simulation in the regions predicted by linear theory to become unstable in the 1D model. Angular velocity and chemical composition is redistributed in the unstable region, thereby creating new unstable regions. Eventually the 2D simulation settles in a symmetric, steady state, which is Richardson stable everywhere, whereas the instability remains for longer in the 1D model due to current limitations in the 1D code. A spatially resolved diffusion coefficient is extracted by comparing the initial and final profiles of mean atomic mass. The presented simulation gives a first insight on hydrodynamics of shear instabilities in a real stellar environment and even allows us to directly extract an effective diffusion coefficient. We see evidence for a critical Richardson number of 0.25 as regions above this value remain stable for the course of the simulation. △ Less

Submitted 20 April, 2017; originally announced April 2017.

Comments: 15 pages, 14 figures, accepted for publication by A&A, movie available at https://slh-code.org/papers/dynshear

Journal ref: A&A 604, A25 (2017)

Convective-reactive nucleosynthesis of K, Sc, Cl and p-process isotopes in O-C shell mergers

Authors: C. Ritter, R. Andrassy, B. Côté, F. Herwig, P. R. Woodward, M. Pignatari, S. Jones

Abstract: We address the deficiency of odd-Z elements P, Cl, K and Sc in galactic chemical evolution models through an investigation of the nucleosynthesis of interacting convective O- and C shells in massive stars. 3D hydrodynamic simulations of O-shell convection with moderate C-ingestion rates show no dramatic deviation from spherical symmetry. We derive a spherically averaged diffusion coefficient for 1… ▽ More We address the deficiency of odd-Z elements P, Cl, K and Sc in galactic chemical evolution models through an investigation of the nucleosynthesis of interacting convective O- and C shells in massive stars. 3D hydrodynamic simulations of O-shell convection with moderate C-ingestion rates show no dramatic deviation from spherical symmetry. We derive a spherically averaged diffusion coefficient for 1D nucleosynthesis simulations which show that such convective-reactive ingestion events can be a production site for P, Cl, K and Sc. An entrainment rate of $10^{-3}{{\rm M}_\odot}$/s features overproduction factors $OP_\mathrm{s} \approx 7$. Full O-C shell mergers in our 1D stellar evolution massive star models have overproduction factors $OP_\mathrm{m}>1 \mathrm{dex}$ but for such cases 3D hydrodynamic simulations suggest deviations from spherical symmetry. p-process species can be produced with overproduction factors of $OP_\mathrm{m}>1 \mathrm{dex}$, e.g. for $^{130,132}Ba$. Using the uncertain prediction of the $15{{\rm M}_\odot}$, $Z=0.02$ massive star model ($OP_\mathrm{m} \approx 15$) as representative for merger or entrainment convective-reactive events involving O- and C-burning shells, and assume that such events occur in more than 50% of all stars, our chemical evolution models reproduce the observed Galactic trends of the odd-Z elements. △ Less

Submitted 25 September, 2017; v1 submitted 19 April, 2017; originally announced April 2017.

Comments: 6 pages, 6 figures, accepted for publication in MNRAS

The impact of (n,$γ$) reaction rate uncertainties of unstable isotopes near $N=50$ on the i process nucleosynthesis in He-shell flash white dwarfs

Authors: Pavel Denissenkov, Georgios Perdikakis, Falk Herwig, Hendrik Schatz, Christian Ritter, Marco Pignatari, Samuel Jones, Stylianos Nikas, Artemis Spyrou

Abstract: The first peak s-process elements Rb, Sr, Y and Zr in the post-AGB star Sakurai's object (V4334 Sagittarii) have been proposed to be the result of i-process nucleosynthesis in a post-AGB very-late thermal pulse event. We estimate the nuclear physics uncertainties in the i-process model predictions to determine whether the remaining discrepancies with observations are significant and point to poten… ▽ More The first peak s-process elements Rb, Sr, Y and Zr in the post-AGB star Sakurai's object (V4334 Sagittarii) have been proposed to be the result of i-process nucleosynthesis in a post-AGB very-late thermal pulse event. We estimate the nuclear physics uncertainties in the i-process model predictions to determine whether the remaining discrepancies with observations are significant and point to potential issues with the underlying astrophysical model. We find that the dominant source in the nuclear physics uncertainties are predictions of neutron capture rates on unstable neutron rich nuclei, which can have uncertainties of more than a factor 20 in the band of the i-process. We use a Monte Carlo variation of 52 neutron capture rates and a 1D multi-zone post-processing model for the i-process in Sakurai's object to determine the cumulative effect of these uncertainties on the final elemental abundance predictions. We find that the nuclear physics uncertainties are large and comparable to observational errors. Within these uncertainties the model predictions are consistent with observations. A correlation analysis of the results of our MC simulations reveals that the strongest impact on the predicted abundances of Rb, Sr, Y and Zr is made by the uncertainties in the (n,$γ$) reaction rates of $^{85}$Br, $^{86}$Br, $^{87}$Kr, $^{88}$Kr, $^{89}$Kr, $^{89}$Rb, $^{89}$Sr, and $^{92}$Sr. (abridged) △ Less

Submitted 9 April, 2018; v1 submitted 3 November, 2016; originally announced November 2016.

Comments: 30 pages, 10 figures, 5 tables, accepted for publication by The Journal of Physics G

A comprehensive long term study of the radio and X-ray Variability of NGC 4051 Paper II

Authors: Sadie Jones, Ian McHardy, Thomas J. Maccarone

Abstract: The origin of the low luminosity radio emission in radio-quiet AGN, is unknown. The detection of a positive correlation between the radio and X-ray emission would imply a jet-like origin, similar to that seen in `hard state' X-ray binary systems. In our previous work, we found no believable radio variability in the well known X-ray bright Seyfert 1 galaxy NGC 4051, despite large amplitude X-ray va… ▽ More The origin of the low luminosity radio emission in radio-quiet AGN, is unknown. The detection of a positive correlation between the radio and X-ray emission would imply a jet-like origin, similar to that seen in `hard state' X-ray binary systems. In our previous work, we found no believable radio variability in the well known X-ray bright Seyfert 1 galaxy NGC 4051, despite large amplitude X-ray variability. In this study we have carefully re-analysed radio and X-ray observations using the same methods as our previous work, we again find no evidence for core radio variability. In direct contrast to our findings, another study claim significant radio variability and a distinctive anti-correlation between radio and X-ray data for the same source. The other study report only integral flux values and do not consider the effect of the changing array on the synthesised beam. In both our studies of NGC 4051 we have taken great care to account for the effect that the changing beam size has on the measured radio flux and as a result we are confident that our method gives more accurate values for the intrinsic core radio flux. However, the lack of radio variability we find is hard to reconcile because radio images of NGC 4051 do show jet-like structure. We suggest that the radio structures observed are likely the result of a previous period of higher radio activity and that the current level of radio emission from a compact nuclear jet is low. △ Less

Submitted 28 October, 2016; originally announced October 2016.

Comments: 13 pages, 8 Figures