Showing 1–50 of 203 results for author: Booth, A

IoT Monitoring with Blockchain: Generating Smart Contracts from Service Level Agreements

Authors: Adam Booth, Awatif Alqahtani, Ellis Solaiman

Abstract: A Service Level Agreement (SLA) is a commitment between a client and provider that assures the quality of service (QoS) a client can expect to receive when purchasing a service. However, evidence of SLA violations in Internet of Things (IoT) service monitoring data can be manipulated by the provider or consumer, resulting in an issue of trust between contracted parties. The following research aims… ▽ More A Service Level Agreement (SLA) is a commitment between a client and provider that assures the quality of service (QoS) a client can expect to receive when purchasing a service. However, evidence of SLA violations in Internet of Things (IoT) service monitoring data can be manipulated by the provider or consumer, resulting in an issue of trust between contracted parties. The following research aims to explore the use of blockchain technology in monitoring IoT systems using smart contracts so that SLA violations captured are irrefutable amongst service providers and clients. The research focuses on the development of a Java library that is capable of generating a smart contract from a given SLA. A smart contract generated by this library is validated through a mock scenario presented in the form of a Remote Patient Monitoring IoT system. In this scenario, the findings demonstrate a 100 percent success rate in capturing all emulated violations. △ Less

Submitted 23 August, 2024; originally announced August 2024.

Journal ref: Managing Internet of Things Applications across Edge and Cloud Data Centres, IET, 2024

The UK Submillimetre and Millimetre Astronomy Roadmap 2024

Authors: K. Pattle, P. S. Barry, A. W. Blain, M. Booth, R. A. Booth, D. L. Clements, M. J. Currie, S. Doyle, D. Eden, G. A. Fuller, M. Griffin, P. G. Huggard, J. D. Ilee, J. Karoly, Z. A. Khan, N. Klimovich, E. Kontar, P. Klaassen, A. J. Rigby, P. Scicluna, S. Serjeant, B. -K. Tan, D. Ward-Thompson, T. G. Williams, T. A. Davis , et al. (9 additional authors not shown)

Abstract: In this Roadmap, we present a vision for the future of submillimetre and millimetre astronomy in the United Kingdom over the next decade and beyond. This Roadmap has been developed in response to the recommendation of the Astronomy Advisory Panel (AAP) of the STFC in the AAP Astronomy Roadmap 2022. In order to develop our stragetic priorities and recommendations, we surveyed the UK submillimetre a… ▽ More In this Roadmap, we present a vision for the future of submillimetre and millimetre astronomy in the United Kingdom over the next decade and beyond. This Roadmap has been developed in response to the recommendation of the Astronomy Advisory Panel (AAP) of the STFC in the AAP Astronomy Roadmap 2022. In order to develop our stragetic priorities and recommendations, we surveyed the UK submillimetre and millimetre community to determine their key priorities for both the near-term and long-term future of the field. We further performed detailed reviews of UK leadership in submillimetre/millimetre science and instrumentation. Our key strategic priorities are as follows: 1. The UK must be a key partner in the forthcoming AtLAST telescope, for which it is essential that the UK remains a key partner in the JCMT in the intermediate term. 2. The UK must maintain, and if possible enhance, access to ALMA and aim to lead parts of instrument development for ALMA2040. Our strategic priorities complement one another: AtLAST (a 50m single-dish telescope) and an upgraded ALMA (a large configurable interferometric array) would be in synergy, not competition, with one another. Both have identified and are working towards the same overarching science goals, and both are required in order to fully address these goals. △ Less

Submitted 3 September, 2024; v1 submitted 23 August, 2024; originally announced August 2024.

Comments: 91 pages plus cover, 38 figures. Submitted to the Science and Technology Facilities Council, August 2024. One figure corrected (v2); new appendix with STFC Q&A; corrected SMA access statement; updated references, acronyms & author list (v3)

DUNE Phase II: Scientific Opportunities, Detector Concepts, Technological Solutions

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, T. Alves, H. Amar, P. Amedo, J. Anderson, C. Andreopoulos, M. Andreotti , et al. (1347 additional authors not shown)

Abstract: The international collaboration designing and constructing the Deep Underground Neutrino Experiment (DUNE) at the Long-Baseline Neutrino Facility (LBNF) has developed a two-phase strategy toward the implementation of this leading-edge, large-scale science project. The 2023 report of the US Particle Physics Project Prioritization Panel (P5) reaffirmed this vision and strongly endorsed DUNE Phase I… ▽ More The international collaboration designing and constructing the Deep Underground Neutrino Experiment (DUNE) at the Long-Baseline Neutrino Facility (LBNF) has developed a two-phase strategy toward the implementation of this leading-edge, large-scale science project. The 2023 report of the US Particle Physics Project Prioritization Panel (P5) reaffirmed this vision and strongly endorsed DUNE Phase I and Phase II, as did the European Strategy for Particle Physics. While the construction of the DUNE Phase I is well underway, this White Paper focuses on DUNE Phase II planning. DUNE Phase-II consists of a third and fourth far detector (FD) module, an upgraded near detector complex, and an enhanced 2.1 MW beam. The fourth FD module is conceived as a "Module of Opportunity", aimed at expanding the physics opportunities, in addition to supporting the core DUNE science program, with more advanced technologies. This document highlights the increased science opportunities offered by the DUNE Phase II near and far detectors, including long-baseline neutrino oscillation physics, neutrino astrophysics, and physics beyond the standard model. It describes the DUNE Phase II near and far detector technologies and detector design concepts that are currently under consideration. A summary of key R&D goals and prototyping phases needed to realize the Phase II detector technical designs is also provided. DUNE's Phase II detectors, along with the increased beam power, will complete the full scope of DUNE, enabling a multi-decadal program of groundbreaking science with neutrinos. △ Less

Submitted 22 August, 2024; originally announced August 2024.

Report number: FERMILAB-TM-2833-LBNF

Monotonic warpings for additive and deep Gaussian processes

Authors: Steven D. Barnett, Lauren J. Beesley, Annie S. Booth, Robert B. Gramacy, Dave Osthus

Abstract: Gaussian processes (GPs) are canonical as surrogates for computer experiments because they enjoy a degree of analytic tractability. But that breaks when the response surface is constrained, say to be monotonic. Here, we provide a mono-GP construction for a single input that is highly efficient even though the calculations are non-analytic. Key ingredients include transformation of a reference proc… ▽ More Gaussian processes (GPs) are canonical as surrogates for computer experiments because they enjoy a degree of analytic tractability. But that breaks when the response surface is constrained, say to be monotonic. Here, we provide a mono-GP construction for a single input that is highly efficient even though the calculations are non-analytic. Key ingredients include transformation of a reference process and elliptical slice sampling. We then show how mono-GP may be deployed effectively in two ways. One is additive, extending monotonicity to more inputs; the other is as a prior on injective latent warping variables in a deep Gaussian process for (non-monotonic, multi-input) non-stationary surrogate modeling. We provide illustrative and benchmarking examples throughout, showing that our methods yield improved performance over the state-of-the-art on examples from those two classes of problems. △ Less

Submitted 2 August, 2024; originally announced August 2024.

First Measurement of the Total Inelastic Cross-Section of Positively-Charged Kaons on Argon at Energies Between 5.0 and 7.5 GeV

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, T. Alves, H. Amar, P. Amedo, J. Anderson, C. Andreopoulos, M. Andreotti , et al. (1341 additional authors not shown)

Abstract: ProtoDUNE Single-Phase (ProtoDUNE-SP) is a 770-ton liquid argon time projection chamber that operated in a hadron test beam at the CERN Neutrino Platform in 2018. We present a measurement of the total inelastic cross section of charged kaons on argon as a function of kaon energy using 6 and 7 GeV/$c$ beam momentum settings. The flux-weighted average of the extracted inelastic cross section at each… ▽ More ProtoDUNE Single-Phase (ProtoDUNE-SP) is a 770-ton liquid argon time projection chamber that operated in a hadron test beam at the CERN Neutrino Platform in 2018. We present a measurement of the total inelastic cross section of charged kaons on argon as a function of kaon energy using 6 and 7 GeV/$c$ beam momentum settings. The flux-weighted average of the extracted inelastic cross section at each beam momentum setting was measured to be 380$\pm$26 mbarns for the 6 GeV/$c$ setting and 379$\pm$35 mbarns for the 7 GeV/$c$ setting. △ Less

Submitted 1 August, 2024; originally announced August 2024.

Report number: CERN-EP-2024-211, FERMILAB-PUB-24-0216-V

Supernova Pointing Capabilities of DUNE

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, T. Alves, H. Amar, P. Amedo, J. Anderson, D. A. Andrade , et al. (1340 additional authors not shown)

Abstract: The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electr… ▽ More The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electron-neutrino charged-current absorption on $^{40}$Ar and elastic scattering of neutrinos on electrons. Procedures to reconstruct individual interactions, including a newly developed technique called ``brems flipping'', as well as the burst direction from an ensemble of interactions are described. Performance of the burst direction reconstruction is evaluated for supernovae happening at a distance of 10 kpc for a specific supernova burst flux model. The pointing resolution is found to be 3.4 degrees at 68% coverage for a perfect interaction-channel classification and a fiducial mass of 40 kton, and 6.6 degrees for a 10 kton fiducial mass respectively. Assuming a 4% rate of charged-current interactions being misidentified as elastic scattering, DUNE's burst pointing resolution is found to be 4.3 degrees (8.7 degrees) at 68% coverage. △ Less

Submitted 14 July, 2024; originally announced July 2024.

Comments: 25 pages, 16 figures

Report number: FERMILAB-PUB-24-0319-LBNF

Viscous circumbinary protoplanetary discs -- I. Structure of the inner cavity

Authors: Anna B. T. Penzlin, Richard A. Booth, Richard P. Nelson, Christoph M. Schäfer, Wilhelm Kley

Abstract: Many of the most intriguing features, including spirals and cavities, in the current disc observations are found in binary systems like GG Tau, HD 142527 or HD 100453. Such features are evidence of the dynamic interaction between binary stars and the viscous disc. Understanding these dynamic interactions and how they result in the structure and growth of asymmetric circumbinary discs is a difficul… ▽ More Many of the most intriguing features, including spirals and cavities, in the current disc observations are found in binary systems like GG Tau, HD 142527 or HD 100453. Such features are evidence of the dynamic interaction between binary stars and the viscous disc. Understanding these dynamic interactions and how they result in the structure and growth of asymmetric circumbinary discs is a difficult problem, for which there is no complete analytical solution, that predicts the shape of the observed disc accurately. We use numeric simulation to evolve circumbinary discs with varying disc viscosities and investigate the size and shape of the inner cavities in such protoplanetary discs. We have simulated over 140 locally isothermal 2D grid-based disc models for > 3e4 binary orbits each and mapped out the parameter space relevant for protoplanetary discs. With this, it becomes possible to create parametrised profiles for individual discs to compare to observation and find limits to their binary eccentricity or internal viscosity from the simulation data. In the long-term simulations larger cavity sizes than previously considered are possible within the parameter space (< 6 binary separations). As an example, we find that the eccentricity of the disc around HD 142527 suggests the impact of the binary dynamics on the disc. However, even considering the larger cavity sizes, the large size of the cavity in HD 142527 remains unexplained by the simulations considering the most recent orbital constraints. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 13 pages, 14 figures (appendix: 2 pages, 3 figures); accepted by mnras

BOWIE-ALIGN: How formation and migration histories of giant planets impact atmospheric compositions

Authors: Anna B. T. Penzlin, Richard A. Booth, James Kirk, James E. Owen, Eva-Maria Ahrer, Duncan A. Christie, Alastair B. Claringbold, Emma Esparza-Borges, M. López-Morales, N. J. Mayne, Mason McCormack, Annabella Meech, Vatsal Panwar, Diana Powell, Denis E. Sergeev, Jake Taylor, Peter J. Wheatley, Maria Zamyatina

Abstract: Hot Jupiters present a unique opportunity for measuring how planet formation history shapes present-day atmospheric composition. However, due to the myriad pathways influencing composition, a well-constructed sample of planets is needed to determine whether formation history can be accurately traced back from atmospheric composition. To this end, the BOWIE-ALIGN survey will compare the composition… ▽ More Hot Jupiters present a unique opportunity for measuring how planet formation history shapes present-day atmospheric composition. However, due to the myriad pathways influencing composition, a well-constructed sample of planets is needed to determine whether formation history can be accurately traced back from atmospheric composition. To this end, the BOWIE-ALIGN survey will compare the compositions of 8 hot Jupiters around F stars, 4 with orbits aligned with the stellar rotation axis and 4 misaligned. Using the alignment as an indicator for planets that underwent disc migration or high-eccentricity migration, one can determine whether migration history produces notable differences in composition between the two samples of planets. This paper describes the planet formation model that motivates our observing programme. Our model traces the accretion of chemical components from the gas and dust in the disc over a broad parameter space to create a full, unbiased model sample from which we can estimate the range of final atmospheric compositions. For high metallicity atmospheres (O/H > 10 times solar), the C/O ratios of aligned and misaligned planets diverge, with aligned planets having lower C/O (< 0.25) due to the accretion of oxygen-rich silicates from the inner disc. However, silicates may rain out instead of releasing their oxygen into the atmosphere. This would significantly increase the C/O of aligned planets (C/O > 0.6), inverting the trend between the aligned and misaligned planets. Nevertheless, by comparing statistically significant samples of aligned and misaligned planets, we expect atmospheric composition to constrain how planets form. △ Less

Submitted 4 July, 2024; v1 submitted 3 July, 2024; originally announced July 2024.

Comments: 11pages 10 figures, (appendix: 6 page, 4 figures), submitted to mnras

BOWIE-ALIGN: A JWST comparative survey of aligned vs misaligned hot Jupiters to test the dependence of atmospheric composition on migration history

Authors: James Kirk, Eva-Maria Ahrer, Anna B. T. Penzlin, James E. Owen, Richard A. Booth, Lili Alderson, Duncan A. Christie, Alastair B. Claringbold, Emma Esparza-Borges, Chloe E. Fisher, Mercedes López-Morales, N. J. Mayne, Mason McCormack, Annabella Meech, Vatsal Panwar, Diana Powell, Jake Taylor, Denis E. Sergeev, Daniel Valentine, Hannah R. Wakeford, Peter J. Wheatley, Maria Zamyatina

Abstract: A primary objective of exoplanet atmosphere characterisation is to learn about planet formation and evolution, however, this is challenged by degeneracies. To determine whether differences in atmospheric composition can be reliably traced to differences in evolution, we are undertaking a new survey with JWST to compare the compositions of a sample of hot Jupiters that orbit F stars above the Kraft… ▽ More A primary objective of exoplanet atmosphere characterisation is to learn about planet formation and evolution, however, this is challenged by degeneracies. To determine whether differences in atmospheric composition can be reliably traced to differences in evolution, we are undertaking a new survey with JWST to compare the compositions of a sample of hot Jupiters that orbit F stars above the Kraft break with different orbital alignments. Under the assumption that aligned planets migrate through the inner disc, while misaligned planets migrate after disc dispersal, the act of migrating through the inner disc should lead to a measurable difference in the C/O between aligned and misaligned planets. We expect the amplitude and sign of this difference to depend on the amount of planetesimal accretion and whether silicates accreted from the inner disc release their oxygen. Here, we identify all known exoplanets that are suitable for testing this hypothesis, describe our JWST survey, and use noise simulations and atmospheric retrievals to estimate our survey's sensitivity. With the selected sample of four aligned and four misaligned hot Jupiters, we will be sensitive to the predicted differences in C/O between aligned and misaligned hot Jupiters for a wide range of model scenarios. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 13 pages, 8 figures, submitted to RASTI

First Measurement of Deeply Virtual Compton Scattering on the Neutron with Detection of the Active Neutron

Authors: CLAS Collaboration, A. Hobart, S. Niccolai, M. Čuić, K. Kumerički, P. Achenbach, J. S. Alvarado, W. R. Armstrong, H. Atac, H. Avakian, L. Baashen, N. A. Baltzell, L. Barion, M. Bashkanov, M. Battaglieri, B. Benkel, F. Benmokhtar, A. Bianconi, A. S. Biselli, S. Boiarinov, M. Bondi, W. A. Booth, F. Bossù, K. -Th. Brinkmann, W. J. Briscoe , et al. (124 additional authors not shown)

Abstract: Measuring Deeply Virtual Compton Scattering on the neutron is one of the necessary steps to understand the structure of the nucleon in terms of Generalized Parton Distributions (GPDs). Neutron targets play a complementary role to transversely polarized proton targets in the determination of the GPD $E$. This poorly known and poorly constrained GPD is essential to obtain the contribution of the qua… ▽ More Measuring Deeply Virtual Compton Scattering on the neutron is one of the necessary steps to understand the structure of the nucleon in terms of Generalized Parton Distributions (GPDs). Neutron targets play a complementary role to transversely polarized proton targets in the determination of the GPD $E$. This poorly known and poorly constrained GPD is essential to obtain the contribution of the quarks' angular momentum to the spin of the nucleon. DVCS on the neutron was measured for the first time selecting the exclusive final state by detecting the neutron, using the Jefferson Lab longitudinally polarized electron beam, with energies up to 10.6 GeV, and the CLAS12 detector. The extracted beam-spin asymmetries, combined with DVCS observables measured on the proton, allow a clean quark-flavor separation of the imaginary parts of the GPDs $H$ and $E$. △ Less

Submitted 25 June, 2024; v1 submitted 21 June, 2024; originally announced June 2024.

Comments: 7 pages, 6 figures

Report number: JLAB-PHY-24-4089

C/O Ratios and the formation of wide separation exoplanets

Authors: Edwin A. Bergin, Richard A. Booth, Maria Jose Colmenares, John D. Ilee

Abstract: The gas and solid-state C/O ratios provide context to potentially link the atmospheric composition of planets to that of the natal disk. We provide a synthesis of extant estimates of the gaseous C/O and C/H ratios in planet-forming disks obtained primarily through analysis of Atacama Large Millimeter Array observations. These estimates are compared to atmospheric abundances of wide separation (> 1… ▽ More The gas and solid-state C/O ratios provide context to potentially link the atmospheric composition of planets to that of the natal disk. We provide a synthesis of extant estimates of the gaseous C/O and C/H ratios in planet-forming disks obtained primarily through analysis of Atacama Large Millimeter Array observations. These estimates are compared to atmospheric abundances of wide separation (> 10 au) gas giants. The resolved disk gas C/O ratios, from seven systems, generally exhibit C/O > 1 with sub-solar, or depleted, carbon content. In contrast, wide separation gas giants have atmospheric C/O ratios that cluster near or slightly above the presumed stellar value with a range of elemental C/H. From the existing disk composition, we infer that the solid-state mm/cm-sized pebbles have a total C/O ratio (solid cores and ices) that is solar (stellar) in content. We explore simple models that reconstruct the exoplanet atmospheric composition from the disk, while accounting for silicate cloud formation in the planet atmosphere. If wide separate planets formed via the core-accretion mechanism, they must acquire their metals from pebble or planetesimal accretion. Further, the dispersion in giant planet C/H content is best matched by a disk composition with modest and variable factors of carbon depletion. An origin of the wide separation gas giants via gravitational instability cannot be ruled out as stellar C/O ratios should natively form in this scenario. However, the variation in planet metallicity with a stellar C/O ratio potentially presents challenges to these models. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 17 pages, 6 figures, accepted by Astrophysical Journal Letters

Seeing the unseen: a method to detect unresolved rings in protoplanetary disks

Authors: Chiara E. Scardoni, Richard A. Booth, Cathie J. Clarke, Giovanni P. Rosotti, Alvaro Ribas

Abstract: While high resolution ALMA observations reveal a wealth of substructure in protoplanetary discs, they remain incapable of resolving the types of small scale dust structures predicted, for example, by numerical simulations of the streaming instability. In this Letter, we propose a method to find evidence for unresolved, optically thick dusty rings in protoplanetary disks. We demonstrate that, in pr… ▽ More While high resolution ALMA observations reveal a wealth of substructure in protoplanetary discs, they remain incapable of resolving the types of small scale dust structures predicted, for example, by numerical simulations of the streaming instability. In this Letter, we propose a method to find evidence for unresolved, optically thick dusty rings in protoplanetary disks. We demonstrate that, in presence of unresolved rings, the brightness of an inclined disc exhibits a distinctive emission peak at the minor axis. Furthermore, the azimuthal brightness depends on both the geometry of the rings and the dust optical properties; we can therefore use the azimuthal brightness variations to both detect unresolved rings and probe their properties. By analyzing the azimuthal brightness in the test-case of ring-like substructures formed by streaming instability, we show that the resulting peak is likely detectable by ALMA for typical disc parameters. Moreover, we present an analytic model that not only qualitatively but also quantitatively reproduces the peak found in the simulations, validating its applicability to infer the presence of unresolved rings in observations and characterize their optical properties and shape. This will contribute to the identification of disk regions where streaming instability (and thus planet formation) is occurring. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: Accepted by ApJ, 10 pages, 6 figures

Faraday tomography with CHIME: the `tadpole' feature G137+7

Authors: Nasser Mohammed, Anna Ordog, Rebecca A. Booth, Andrea Bracco, Jo-Anne C. Brown, Ettore Carretti, John M. Dickey, Simon Foreman, Mark Halpern, Marijke Haverkorn, Alex S. Hill, Gary Hinshaw, Joseph W Kania, Roland Kothes, T. L. Landecker, Joshua MacEachern, Kiyoshi W. Masui, Aimee Menard, Ryan R. Ransom, Wolfgang Reich, Patricia Reich, J. Richard Shaw, Seth R. Siegel, Mehrnoosh Tahani, Alec J. M. Thomson , et al. (5 additional authors not shown)

Abstract: A direct consequence of Faraday rotation is that the polarized radio sky does not resemble the total intensity sky at long wavelengths. We analyze G137+7, which is undetectable in total intensity but appears as a depolarization feature. We use the first polarization maps from the Canadian Hydrogen Intensity Mapping Experiment. Our $400-729$ MHz bandwidth and angular resolution, $17'$ to $30'$, all… ▽ More A direct consequence of Faraday rotation is that the polarized radio sky does not resemble the total intensity sky at long wavelengths. We analyze G137+7, which is undetectable in total intensity but appears as a depolarization feature. We use the first polarization maps from the Canadian Hydrogen Intensity Mapping Experiment. Our $400-729$ MHz bandwidth and angular resolution, $17'$ to $30'$, allow us to use Faraday synthesis to analyze the polarization structure. In polarized intensity and polarization angle maps, we find a "tail" extending $10^\circ$ from the "head" and designate the combined object the "tadpole". Similar polarization angles, distinct from the background, indicate that the head and tail are physically associated. The head appears as a depolarized ring in single channels, but wideband observations show that it is a Faraday rotation feature. Our investigations of H I and H$α$ find no connections to the tadpole. The tail suggests motion of either the gas or an ionizing star through the ISM; the B2(e) star HD 20336 is a candidate. While the head features a coherent, $\sim -8$ rad m$^2$ Faraday depth, Faraday synthesis also identifies multiple components in both the head and tail. We verify the locations of the components in the spectra using QU fitting. Our results show that $\sim$octave-bandwidth Faraday rotation observations at $\sim 600$ MHz are sensitive to low-density ionized or partially-ionized gas which is undetectable in other tracers. △ Less

Submitted 31 July, 2024; v1 submitted 24 May, 2024; originally announced May 2024.

Comments: ApJ in press. Replacement corrects typographical error in equation 6

Chemistry across dust and gas gaps in protoplanetary disks: modelling the co-spatial molecular rings in the HD 100546 disk

Authors: M. Leemker, A. S. Booth, E. F. van Dishoeck, L. Wölfer, B. Dent

Abstract: High-resolution observations show that typically both the dust and the gas in nearby extended protoplanetary disks are structured, possibly related to radial and azimuthal variations in the disk density and/or chemistry. The aim of this work is to identify the expected location and intensity of rings seen in molecular line emission of HCN, CN, C$_2$H, NO, [CI], and HCO$^+$ in gapped disks while ex… ▽ More High-resolution observations show that typically both the dust and the gas in nearby extended protoplanetary disks are structured, possibly related to radial and azimuthal variations in the disk density and/or chemistry. The aim of this work is to identify the expected location and intensity of rings seen in molecular line emission of HCN, CN, C$_2$H, NO, [CI], and HCO$^+$ in gapped disks while exploring a range of physical conditions across the gap. In particular, we model HD 100546 disk where molecular rings are co-spatial with the dust rings at $\sim$20 and $\sim$200 au, in contrast to most other gapped disks. The fiducial model of a gapped disk with a 15 au gas cavity, a 20 au dust cavity, and a shallow (a factor of $\lesssim10$) gas and deep dust gap at 40-175 au provides a good fit to the continuum, CO isotopologues, HCN, and HCO$^+$ in the HD 100546 disk. However, the predictions for [CI], CN, C$_2$H and NO do not match the intensity nor the morphology of the observations. An exploration of the parameter space shows that in general the molecular emission rings are only co-spatial with the dust rings if the gas gap between the dust rings is depleted by at least four orders of magnitude in gas or if the C/O ratio of the gas is varying as a function of radius. For shallower gaps the decrease in the UV field roughly balances the effect of a higher gas density for UV tracers such as CN, C$_2$H, and NO. Therefore, these radicals are not good tracers of the gas gap depth. The C/O ratio primarily effects the intensity of the lines without changing the morphology much. The co-spatial rings observed in the HD 100546 disk could be indicative of a radially varying C/O ratio in the HD 100546 disk with a C/O above 1 in a narrow region across the dust rings, together with a shallow gas gap that is depleted by a factor of $\sim$10 in gas, and a reduced background UV field. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: Accepted for publication in A&A

Introducing two improved methods for approximating radiative cooling in hydrodynamical simulations of accretion discs

Authors: Alison K. Young, Maggie Celeste, Richard A. Booth, Ken Rice, Adam Koval, Ethan Carter, Dimitris Stamatellos

Abstract: The evolution of many astrophysical systems depends strongly on the balance between heating and cooling, in particular star formation in giant molecular clouds and the evolution of young protostellar systems. Protostellar discs are susceptible to the gravitational instability, which can play a key role in their evolution and in planet formation. The strength of the instability depends on the rate… ▽ More The evolution of many astrophysical systems depends strongly on the balance between heating and cooling, in particular star formation in giant molecular clouds and the evolution of young protostellar systems. Protostellar discs are susceptible to the gravitational instability, which can play a key role in their evolution and in planet formation. The strength of the instability depends on the rate at which the system loses thermal energy. To study the evolution of these systems, we require radiative cooling approximations because full radiative transfer is generally too expensive to be coupled to hydrodynamical models. Here we present two new approximate methods for computing radiative cooling that make use of the polytropic cooling approximation. This approach invokes the assumption that each parcel of gas is located within a spherical pseudo-cloud which can then be used to approximate the optical depth. The first method combines the methods introduced by Stamatellos et al. and Lombardi et al. to overcome the limitations of each method at low and high optical depths respectively. The second, the "Modified Lombardi" method, is specifically tailored for self-gravitating discs. This modifies the scale height estimate from the method of Lombardi et al. using the analytical scale height for a self-gravitating disc. We show that the Modified Lombardi method provides an excellent approximation for the column density in a fragmenting disc, a regime in which the existing methods fail to recover the clumps and spiral structures. We therefore recommend this improved radiative cooling method for more realistic simulations of self-gravitating discs. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: Accepted to MNRAS. 10 pages

Searching for planet-driven dust spirals in ALMA visibilities

Authors: Edward T. Stevenson, Álvaro Ribas, Jessica Speedie, Richard A. Booth, Cathie J. Clarke

Abstract: ALMA (Atacama Large Millimetre/submillimetre Array) observations of the thermal emission from protoplanetary disc dust have revealed a wealth of substructures that could evidence embedded planets, but planet-driven spirals, one of the more compelling lines of evidence, remain relatively rare. Existing works have focused on detecting these spirals using methods that operate in image space. Here, we… ▽ More ALMA (Atacama Large Millimetre/submillimetre Array) observations of the thermal emission from protoplanetary disc dust have revealed a wealth of substructures that could evidence embedded planets, but planet-driven spirals, one of the more compelling lines of evidence, remain relatively rare. Existing works have focused on detecting these spirals using methods that operate in image space. Here, we explore the planet detection capabilities of fitting planet-driven spirals to disc observations directly in visibility space. We test our method on synthetic ALMA observations of planet-containing model discs for a range of disc/observational parameters, finding it significantly outperforms image residuals in identifying spirals in these observations and is able to identify spirals in regions of the parameter space in which no gaps are detected. These tests suggest that a visibility-space fitting approach warrants further investigation and may be able to find planet-driven spirals in observations that have not yet been found with existing approaches. We also test our method on six discs in the Taurus molecular cloud observed with ALMA at 1.33 mm, but find no evidence for planet-driven spirals. We find that the minimum planet masses necessary to drive detectable spirals range from $\approx$ 0.03 to 0.5 $M_{\text{Jup}}$ over orbital radii of 10 to 100 au, with planet masses below these thresholds potentially hiding in such disc observations. Conversely, we suggest that planets $\gtrsim$ 0.5 to 1 $M_{\text{Jup}}$ can likely be ruled out over orbital radii of $\approx$ 20 to 60 au on the grounds that we would have detected them if they were present. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 16 pages, 14 figures, 3 tables. Accepted 2024 April 18 for publication in MNRAS

Search for CP-violating Neutrino Non-Standard Interactions with the NOvA Experiment

Authors: NOvA Collaboration, M. A. Acero, B. Acharya, P. Adamson, L. Aliaga, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, N. Balashov, P. Baldi, B. A. Bambah, A. Bat, K. Bays, R. Bernstein, T. J. C. Bezerra, V. Bhatnagar, D. Bhattarai, B. Bhuyan, J. Bian, A. C. Booth, R. Bowles, B. Brahma , et al. (182 additional authors not shown)

Abstract: This Letter reports a search for charge-parity (CP) symmetry violating non-standard interactions (NSI) of neutrinos with matter using the NOvA Experiment, and examines their effects on the determination of the standard oscillation parameters. Data from $ν_μ(\barν_μ)\rightarrowν_μ(\barν_μ)$ and $ν_μ(\barν_μ)\rightarrowν_{e}(\barν_{e})$ oscillation channels are used to measure the effect of the NSI… ▽ More This Letter reports a search for charge-parity (CP) symmetry violating non-standard interactions (NSI) of neutrinos with matter using the NOvA Experiment, and examines their effects on the determination of the standard oscillation parameters. Data from $ν_μ(\barν_μ)\rightarrowν_μ(\barν_μ)$ and $ν_μ(\barν_μ)\rightarrowν_{e}(\barν_{e})$ oscillation channels are used to measure the effect of the NSI parameters $\varepsilon_{eμ}$ and $\varepsilon_{eτ}$. With 90% C.L. the magnitudes of the NSI couplings are constrained to be $|\varepsilon_{eμ}| \, \lesssim 0.3$ and $|\varepsilon_{eτ}| \, \lesssim 0.4$. A degeneracy at $|\varepsilon_{eτ}| \, \approx 1.8$ is reported, and we observe that the presence of NSI limits sensitivity to the standard CP phase $δ_{\tiny\text{CP}}$. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Report number: FERMILAB-PUB-24-0108-PPD

Performance of a modular ton-scale pixel-readout liquid argon time projection chamber

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, T. Alves, H. Amar, P. Amedo, J. Anderson, D. A. Andrade , et al. (1340 additional authors not shown)

Abstract: The Module-0 Demonstrator is a single-phase 600 kg liquid argon time projection chamber operated as a prototype for the DUNE liquid argon near detector. Based on the ArgonCube design concept, Module-0 features a novel 80k-channel pixelated charge readout and advanced high-coverage photon detection system. In this paper, we present an analysis of an eight-day data set consisting of 25 million cosmi… ▽ More The Module-0 Demonstrator is a single-phase 600 kg liquid argon time projection chamber operated as a prototype for the DUNE liquid argon near detector. Based on the ArgonCube design concept, Module-0 features a novel 80k-channel pixelated charge readout and advanced high-coverage photon detection system. In this paper, we present an analysis of an eight-day data set consisting of 25 million cosmic ray events collected in the spring of 2021. We use this sample to demonstrate the imaging performance of the charge and light readout systems as well as the signal correlations between the two. We also report argon purity and detector uniformity measurements, and provide comparisons to detector simulations. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 47 pages, 41 figures

Report number: FERMILAB-PUB-24-0073-LBNF

Introducing cuDisc: a 2D code for protoplanetary disc structure and evolution calculations

Authors: Alfie Robinson, Richard A. Booth, James E. Owen

Abstract: We present a new 2D axisymmetric code, cuDisc, for studying protoplanetary discs, focusing on the self-consistent calculation of dust dynamics, grain size distribution and disc temperature. Self-consistently studying these physical processes is essential for many disc problems, such as structure formation and dust removal, given that the processes heavily depend on one another. To follow the evolu… ▽ More We present a new 2D axisymmetric code, cuDisc, for studying protoplanetary discs, focusing on the self-consistent calculation of dust dynamics, grain size distribution and disc temperature. Self-consistently studying these physical processes is essential for many disc problems, such as structure formation and dust removal, given that the processes heavily depend on one another. To follow the evolution over substantial fractions of the disc lifetime, cuDisc uses the CUDA language and libraries to speed up the code through GPU acceleration. cuDisc employs a second-order finite-volume Godonuv solver for dust dynamics, solves the Smoluchowski equation for dust growth and calculates radiative transfer using a multi-frequency hybrid ray-tracing/flux-limited-diffusion method. We benchmark our code against current state-of-the-art codes. Through studying steady-state problems, we find that including 2D structure reveals that when collisions are important, the dust vertical structure appears to reach a diffusion-settling-coagulation equilibrium that can differ substantially from standard models that ignore coagulation. For low fragmentation velocities, we find an enhancement of intermediate-sized dust grains at heights of ~ 1 gas scale height due to the variation in collision rates with height, and for large fragmentation velocities, we find an enhancement of small grains around the disc mid-plane due to collisional ''sweeping'' of small grains by large grains. These results could be important for the analysis of disc SEDs or scattered light images, given these observables are sensitive to the vertical grain distribution. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Comments: 18 pages, accepted for publication in MNRAS. GitHub repository link: https://github.com/cuDisc/cuDisc

Voronoi Candidates for Bayesian Optimization

Authors: Nathan Wycoff, John W. Smith, Annie S. Booth, Robert B. Gramacy

Abstract: Bayesian optimization (BO) offers an elegant approach for efficiently optimizing black-box functions. However, acquisition criteria demand their own challenging inner-optimization, which can induce significant overhead. Many practical BO methods, particularly in high dimension, eschew a formal, continuous optimization of the acquisition function and instead search discretely over a finite set of s… ▽ More Bayesian optimization (BO) offers an elegant approach for efficiently optimizing black-box functions. However, acquisition criteria demand their own challenging inner-optimization, which can induce significant overhead. Many practical BO methods, particularly in high dimension, eschew a formal, continuous optimization of the acquisition function and instead search discretely over a finite set of space-filling candidates. Here, we propose to use candidates which lie on the boundary of the Voronoi tessellation of the current design points, so they are equidistant to two or more of them. We discuss strategies for efficient implementation by directly sampling the Voronoi boundary without explicitly generating the tessellation, thus accommodating large designs in high dimension. On a battery of test problems optimized via Gaussian processes with expected improvement, our proposed approach significantly improves the execution time of a multi-start continuous search without a loss in accuracy. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: comments very welcome

An ALMA molecular inventory of warm Herbig Ae disks: II. Abundant complex organics and volatile sulphur in the IRS 48 disk

Authors: Alice S. Booth, Milou Temmink, Ewine F. van Dishoeck, Lucy Evans, John D. Ilee, Mihkel Kama, Luke Keyte, Charles J. Law, Margot Leemker, Nienke van der Marel, Hideko Nomura, Shota Notsu, Karin Öberg, Catherine Walsh

Abstract: The Atacama Large Millimeter/submillimeter Array (ALMA) can probe the molecular content of planet-forming disks with unprecedented sensitivity. These observations allow us to build up an inventory of the volatiles available for forming planets and comets. Herbig Ae transition disks are fruitful targets due to the thermal sublimation of complex organic molecule (COM) and likely H2O-rich ices in the… ▽ More The Atacama Large Millimeter/submillimeter Array (ALMA) can probe the molecular content of planet-forming disks with unprecedented sensitivity. These observations allow us to build up an inventory of the volatiles available for forming planets and comets. Herbig Ae transition disks are fruitful targets due to the thermal sublimation of complex organic molecule (COM) and likely H2O-rich ices in these disks. The IRS 48 disk shows a particularly rich chemistry that can be directly linked to its asymmetric dust trap. Here, we present ALMA observations of the IRS 48 disk where we detect 16 different molecules and make the first robust detections of H213CO, 34SO, 33SO and c-H2COCH2 (ethylene oxide) in a protoplanetary disk. All of the molecular emissions, aside from CO, are colocated with the dust trap and this includes newly detected simple molecules such as HCO+, HCN and CS. Interestingly, there are spatial offsets between different molecular families, including between the COMs and sulphur-bearing species, with the latter being more azimuthally extended and located radially further from the star. The abundances of the newly detected COMs relative to CH3OH are higher than the expected protostellar ratios, which implies some degree of chemical processing of the inherited ices during the disk lifetime. These data highlight IRS 48 as a unique astrochemical laboratory to unravel the full volatile reservoir at the epoch of planet and comet formation and the role of the disk in (re)setting chemical complexity. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: Accepted to AJ, 21 pages, 7 figures

An ALMA molecular inventory of warm Herbig Ae disks: I. Molecular rings, asymmetries and complexity in the HD 100546 disk

Authors: Alice S. Booth, Margot Leemker, Ewine F. van Dishoeck, Lucy Evans, John D. Ilee, Mihkel Kama, Luke Keyte, Charles J. Law, Nienke van der Marel, Hideko Nomura, Shota Notsu, Karin Öberg, Milou Temmink, Catherine Walsh

Abstract: Observations of disks with the Atacama Large Millimeter/submillimeter Array (ALMA) allow us to map the chemical makeup of nearby protoplanetary disks with unprecedented spatial resolution and sensitivity. The typical outer Class II disk observed with ALMA is one with an elevated C/O ratio and a lack of oxygen-bearing complex organic molecules, but there are now some interesting exceptions: three t… ▽ More Observations of disks with the Atacama Large Millimeter/submillimeter Array (ALMA) allow us to map the chemical makeup of nearby protoplanetary disks with unprecedented spatial resolution and sensitivity. The typical outer Class II disk observed with ALMA is one with an elevated C/O ratio and a lack of oxygen-bearing complex organic molecules, but there are now some interesting exceptions: three transition disks around Herbig Ae stars all show oxygen-rich gas traced via the unique detections of the molecules SO and CH3OH. We present the first results of an ALMA line survey at 337 to 357 GHz of such disks and focus this paper on the first Herbig Ae disk to exhibit this chemical signature - HD 100546. In these data, we detect 19 different molecules including NO, SO and CH3OCHO (methyl formate). We also make the first tentative detections of H213CO and 34SO in protoplanetary disks. Multiple molecular species are detected in rings, which are, surprisingly, all peaking just beyond the underlying millimeter continuum ring at 200 au. This result demonstrates a clear connection between the large dust distribution and the chemistry in this flat outer disk. We discuss the physical and/or chemical origin of these sub-structures in relation to ongoing planet formation in the HD 100546 disk. We also investigate how similar and/or different the molecular make up of this disk is to other chemically well-characterised Herbig Ae disks. The line-rich data we present motivates the need for more ALMA line surveys to probe the observable chemistry in Herbig Ae systems which offer unique insight into the composition of disk ices, including complex organic molecules. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: Accepted to AJ, 25 pages, 11 figures

Doping Liquid Argon with Xenon in ProtoDUNE Single-Phase: Effects on Scintillation Light

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, H. Amar Es-sghir, P. Amedo, J. Anderson, D. A. Andrade, C. Andreopoulos , et al. (1297 additional authors not shown)

Abstract: Doping of liquid argon TPCs (LArTPCs) with a small concentration of xenon is a technique for light-shifting and facilitates the detection of the liquid argon scintillation light. In this paper, we present the results of the first doping test ever performed in a kiloton-scale LArTPC. From February to May 2020, we carried out this special run in the single-phase DUNE Far Detector prototype (ProtoDUN… ▽ More Doping of liquid argon TPCs (LArTPCs) with a small concentration of xenon is a technique for light-shifting and facilitates the detection of the liquid argon scintillation light. In this paper, we present the results of the first doping test ever performed in a kiloton-scale LArTPC. From February to May 2020, we carried out this special run in the single-phase DUNE Far Detector prototype (ProtoDUNE-SP) at CERN, featuring 720 t of total liquid argon mass with 410 t of fiducial mass. A 5.4 ppm nitrogen contamination was present during the xenon doping campaign. The goal of the run was to measure the light and charge response of the detector to the addition of xenon, up to a concentration of 18.8 ppm. The main purpose was to test the possibility for reduction of non-uniformities in light collection, caused by deployment of photon detectors only within the anode planes. Light collection was analysed as a function of the xenon concentration, by using the pre-existing photon detection system (PDS) of ProtoDUNE-SP and an additional smaller set-up installed specifically for this run. In this paper we first summarize our current understanding of the argon-xenon energy transfer process and the impact of the presence of nitrogen in argon with and without xenon dopant. We then describe the key elements of ProtoDUNE-SP and the injection method deployed. Two dedicated photon detectors were able to collect the light produced by xenon and the total light. The ratio of these components was measured to be about 0.65 as 18.8 ppm of xenon were injected. We performed studies of the collection efficiency as a function of the distance between tracks and light detectors, demonstrating enhanced uniformity of response for the anode-mounted PDS. We also show that xenon doping can substantially recover light losses due to contamination of the liquid argon by nitrogen. △ Less

Submitted 2 August, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: 36 pages, 20 figures. Corrected author list; corrected typos across paper and polished text

Report number: CERN-EP-2024-024; FERMILAB-PUB-23-0819-LBNF

JWST-MIRI Spectroscopy of Warm Molecular Emission and Variability in the AS 209 Disk

Authors: Carlos E. Muñoz-Romero, Karin I. Öberg, Andrea Banzatti, Klaus M. Pontoppidan, Sean M. Andrews, David J. Wilner, Edwin A. Bergin, Ian Czekala, Charles J. Law, Colette Salyk, Richard Teague, Chunhua Qi, Jennifer B. Bergner, Jane Huang, Catherine Walsh, Viviana V. Guzmán, L. Ilsedore Cleeves, Yuri Aikawa, Jaehan Bae, Alice S. Booth, Gianni Cataldi, John D. Ilee, Romane Le Gal, Feng Long, Ryan A. Loomis , et al. (2 additional authors not shown)

Abstract: We present MIRI MRS observations of the large, multi-gapped protoplanetary disk around the T-Tauri star AS 209. The observations reveal hundreds of water vapor lines from 4.9 to 25.5 $μ$m towards the inner $\sim1$ au in the disk, including the first detection of ro-vibrational water emission in this disk. The spectrum is dominated by hot ($\sim800$ K) water vapor and OH gas, with only marginal det… ▽ More We present MIRI MRS observations of the large, multi-gapped protoplanetary disk around the T-Tauri star AS 209. The observations reveal hundreds of water vapor lines from 4.9 to 25.5 $μ$m towards the inner $\sim1$ au in the disk, including the first detection of ro-vibrational water emission in this disk. The spectrum is dominated by hot ($\sim800$ K) water vapor and OH gas, with only marginal detections of CO$_2$, HCN, and a possible colder water vapor component. Using slab models with a detailed treatment of opacities and line overlap, we retrieve the column density, emitting area, and excitation temperature of water vapor and OH, and provide upper limits for the observable mass of other molecules. Compared to MIRI spectra of other T-Tauri disks, the inner disk of AS 209 does not appear to be atypically depleted in CO$_2$ nor HCN. Based on \textit{Spitzer IRS} observations, we further find evidence for molecular emission variability over a 10-year baseline. Water, OH, and CO$_2$ line luminosities have decreased by factors 2-4 in the new MIRI epoch, yet there are minimal continuum emission variations. The origin of this variability is yet to be understood. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJ

The Role of Drag and Gravity on Dust Concentration in a Gravitationally Unstable Disc

Authors: Sahl Rowther, Rebecca Nealon, Farzana Meru, James Wurster, Hossam Aly, Richard Alexander, Ken Rice, Richard A. Booth

Abstract: We carry out three dimensional smoothed particle hydrodynamics simulations to study the role of gravitational and drag forces on the concentration of large dust grains (St > 1) in the spiral arms of gravitationally unstable protoplanetary discs, and the resulting implications for planet formation. We find that both drag and gravity play an important role in the evolution of large dust grains. If w… ▽ More We carry out three dimensional smoothed particle hydrodynamics simulations to study the role of gravitational and drag forces on the concentration of large dust grains (St > 1) in the spiral arms of gravitationally unstable protoplanetary discs, and the resulting implications for planet formation. We find that both drag and gravity play an important role in the evolution of large dust grains. If we include both, grains that would otherwise be partially decoupled will become well coupled and trace the spirals. For the dust grains most influenced by drag (with Stokes numbers near unity), the dust disc quickly becomes gravitationally unstable and rapidly forms clumps with masses between 0.15 - 6 Earth masses. A large fraction of clumps are below the threshold where runaway gas accretion can occur. However, if dust self-gravity is neglected, the dust is unable to form clumps, despite still becoming trapped in the gas spirals. When large dust grains are unable to feel either gas gravity or drag, the dust is unable to trace the gas spirals. Hence, full physics is needed to properly simulate dust in gravitationally unstable discs. Dust trapping of large grains in spiral arms of discs stable to gas fragmentation could explain planet formation in very young discs by a population of planetesimals formed due to the combined roles of drag and gravity in the earliest stages of a disc's evolution. Furthermore, it highlights that gravitationally unstable discs are not just important for forming gas giants quickly, it can also rapidly form Earth mass bodies. △ Less

Submitted 24 January, 2024; v1 submitted 17 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS. Updated references

Constraining the gas mass of Herbig disks using CO isotopologues

Authors: L. M. Stapper, M. R. Hogerheijde, E. F. van Dishoeck, L. Lin, A. Ahmadi, A. S. Booth, S. L. Grant, K. Immer, M. Leemker, A. F. Pérez-Sánchez

Abstract: The total disk mass sets the formation potential for exoplanets. Carbon-monoxide (CO) has been used as a gas mass tracer in T Tauri disks, but was found to be less abundant than expected due to freeze-out and chemical conversion of CO on the surfaces of cold dust grains. The disks around more massive intermediate mass pre-main sequence stars called Herbig disks are likely to be warmer, allowing fo… ▽ More The total disk mass sets the formation potential for exoplanets. Carbon-monoxide (CO) has been used as a gas mass tracer in T Tauri disks, but was found to be less abundant than expected due to freeze-out and chemical conversion of CO on the surfaces of cold dust grains. The disks around more massive intermediate mass pre-main sequence stars called Herbig disks are likely to be warmer, allowing for the possibility of using CO as a more effective total gas mass tracer. Using ALMA archival data and new NOEMA data of 12CO, 13CO, and C18O transitions of 35 Herbig disks within 450 pc, the masses are determined using the thermo-chemical code Dust And LInes (DALI). The majority of Herbig disks for which 13CO and C18O are detected are optically thick in both. Computing the gas mass using a simple optically thin relation between line flux and column density results in an underestimate of the gas mass of at least an order of magnitude compared to the masses obtained with DALI. The inferred gas masses with DALI are consistent with a gas-to-dust ratio of at least 100. These gas-to-dust ratios are two orders of magnitude higher compared to those found for T Tauri disks using similar techniques, even over multiple orders of magnitude in dust mass, illustrating the importance of chemical conversion of CO in colder T Tauri disks. Similar high gas-to-dust ratios are found for Herbig group I and II disks. Since group II disks have dust masses comparable to T Tauri disks, their higher CO gas masses illustrate the determining role of temperature. Compared to debris disks, Herbig disks have four orders of magnitude higher gas masses. At least one Herbig disk, HD 163296, has a detected molecular disk wind, but our investigation has not turned up other detections of the CO disk wind in spite of similar sensitivities. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy and Astrophysics. 25 pages, 19 figures, plus appendices

The DUNE Far Detector Vertical Drift Technology, Technical Design Report

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, H. Amar, P. Amedo, J. Anderson, D. A. Andrade, C. Andreopoulos , et al. (1304 additional authors not shown)

Abstract: DUNE is an international experiment dedicated to addressing some of the questions at the forefront of particle physics and astrophysics, including the mystifying preponderance of matter over antimatter in the early universe. The dual-site experiment will employ an intense neutrino beam focused on a near and a far detector as it aims to determine the neutrino mass hierarchy and to make high-precisi… ▽ More DUNE is an international experiment dedicated to addressing some of the questions at the forefront of particle physics and astrophysics, including the mystifying preponderance of matter over antimatter in the early universe. The dual-site experiment will employ an intense neutrino beam focused on a near and a far detector as it aims to determine the neutrino mass hierarchy and to make high-precision measurements of the PMNS matrix parameters, including the CP-violating phase. It will also stand ready to observe supernova neutrino bursts, and seeks to observe nucleon decay as a signature of a grand unified theory underlying the standard model. The DUNE far detector implements liquid argon time-projection chamber (LArTPC) technology, and combines the many tens-of-kiloton fiducial mass necessary for rare event searches with the sub-centimeter spatial resolution required to image those events with high precision. The addition of a photon detection system enhances physics capabilities for all DUNE physics drivers and opens prospects for further physics explorations. Given its size, the far detector will be implemented as a set of modules, with LArTPC designs that differ from one another as newer technologies arise. In the vertical drift LArTPC design, a horizontal cathode bisects the detector, creating two stacked drift volumes in which ionization charges drift towards anodes at either the top or bottom. The anodes are composed of perforated PCB layers with conductive strips, enabling reconstruction in 3D. Light-trap-style photon detection modules are placed both on the cryostat's side walls and on the central cathode where they are optically powered. This Technical Design Report describes in detail the technical implementations of each subsystem of this LArTPC that, together with the other far detector modules and the near detector, will enable DUNE to achieve its physics goals. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 425 pages; 281 figures Central editing team: A. Heavey, S. Kettell, A. Marchionni, S. Palestini, S. Rajogopalan, R. J. Wilson

Report number: Fermilab Report no: TM-2813-LBNF

Observing planetesimal formation under streaming instability in the rings of HD 163296

Authors: Francesco Zagaria, Cathie J. Clarke, Richard A. Booth, Stefano Facchini, Giovanni P. Rosotti

Abstract: We introduce a new technique to determine the gas turbulence and surface density in bright disc rings, under the assumption that dust growth is limited by turbulent fragmentation at the ring centre. We benchmark this prescription in HD 163296, showing that our measurements are consistent with available turbulence upper limits and agree with independent estimates of the gas surface density within a… ▽ More We introduce a new technique to determine the gas turbulence and surface density in bright disc rings, under the assumption that dust growth is limited by turbulent fragmentation at the ring centre. We benchmark this prescription in HD 163296, showing that our measurements are consistent with available turbulence upper limits and agree with independent estimates of the gas surface density within a factor of two. We combine our results with literature measurements of the dust surface density and grain size to determine the dust-to-gas ratio and Stokes number in the 67 au and 100 au rings. Our estimates suggest that particle clumping is taking place under the effect of streaming instability (SI) in the 100 au ring. Even though in the presence of external isotropic turbulence this process might be hindered, we provide evidence that turbulence is non-isotropic in both rings and likely originating from mechanisms (such as ambipolar diffusion) that could ease particle clumping under SI. Finally, we determine the mass accretion rate under the assumption that the disc is in steady state and turbulence regulates angular momentum transport. Our results are in tension with spectroscopic measurements and suggest that other mechanisms might be responsible for accretion, in qualitative agreement with the detection of a magneto-centrifugal wind in this system. Applying our method to larger samples can be used to statistically assess if SI is a viable mechanism to form planetesimals in bright rings. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 13 pages, 4 figures; accepted for publication on ApJL

Expanding neutrino oscillation parameter measurements in NOvA using a Bayesian approach

Authors: NOvA Collaboration, M. A. Acero, B. Acharya, P. Adamson, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, N. Balashov, P. Baldi, B. A. Bambah, A. Bat, K. Bays, R. Bernstein, T. J. C. Bezerra, V. Bhatnagar, D. Bhattarai, B. Bhuyan, J. Bian, A. C. Booth, R. Bowles, B. Brahma, C. Bromberg , et al. (174 additional authors not shown)

Abstract: NOvA is a long-baseline neutrino oscillation experiment that measures oscillations in charged-current $ν_μ \rightarrow ν_μ$ (disappearance) and $ν_μ \rightarrow ν_{e}$ (appearance) channels, and their antineutrino counterparts, using neutrinos of energies around 2 GeV over a distance of 810 km. In this work we reanalyze the dataset first examined in our previous paper [Phys. Rev. D 106, 032004 (20… ▽ More NOvA is a long-baseline neutrino oscillation experiment that measures oscillations in charged-current $ν_μ \rightarrow ν_μ$ (disappearance) and $ν_μ \rightarrow ν_{e}$ (appearance) channels, and their antineutrino counterparts, using neutrinos of energies around 2 GeV over a distance of 810 km. In this work we reanalyze the dataset first examined in our previous paper [Phys. Rev. D 106, 032004 (2022)] using an alternative statistical approach based on Bayesian Markov Chain Monte Carlo. We measure oscillation parameters consistent with the previous results. We also extend our inferences to include the first NOvA measurements of the reactor mixing angle $θ_{13}$ and the Jarlskog invariant. We use these results to quantify the strength of our inferences about CP violation, as well as to examine the effects of constraints from short-baseline measurements of $θ_{13}$ using antineutrinos from nuclear reactors when making NOvA measurements of $θ_{23}$. Our long-baseline measurement of $θ_{13}$ is also shown to be consistent with the reactor measurements, supporting the general applicability and robustness of the PMNS framework for neutrino oscillations. △ Less

Submitted 27 May, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Comments: 20 pages, 17 figures; version accepted by Phys. Rev. D. Data associated with this paper is available at https://doi.org/10.15484/2349444

Report number: FERMILAB-PUB-23-667-AD-CSAID-ND

Journal ref: Phys.Rev.D 110 (2024) 1, 012005

Formation of misaligned second-generation discs through flyby encounters

Authors: Jeremy L. Smallwood, Rebecca Nealon, Nicolás Cuello, Ruobing Dong, Richard A. Booth

Abstract: Observations reveal protoplanetary discs being perturbed by flyby candidates. We simulate a scenario where an unbound perturber, i.e., a flyby, undergoes an inclined grazing encounter, capturing material and forming a second-generation protoplanetary disc. We run $N$--body and three-dimensional hydrodynamical simulations of a parabolic flyby grazing a particle disc and a gas-rich protoplanetary di… ▽ More Observations reveal protoplanetary discs being perturbed by flyby candidates. We simulate a scenario where an unbound perturber, i.e., a flyby, undergoes an inclined grazing encounter, capturing material and forming a second-generation protoplanetary disc. We run $N$--body and three-dimensional hydrodynamical simulations of a parabolic flyby grazing a particle disc and a gas-rich protoplanetary disc, respectively. In both our $N$--body and hydrodynamic simulations, we find that the captured, second-generation disc forms at a tilt twice the initial flyby tilt. This relationship is robust to variations in the flyby's tilt, position angle, periastron, and mass. We extend this concept by also simulating the case where the flyby has a disc of material prior to the encounter but we do not find the same trend. An inclined disc with respect to the primary disc around a misaligned flyby is tilted by a few degrees, remaining close to its initial disc tilt. Therefore, if a disc is present around the flyby before the encounter, the disc may not tilt up to twice the perturber tilt depending on the balance between the angular momentum of the circumsecondary disc and captured particles. In the case where the perturber has no initial disc, analyzing the orientation of these second-generation discs can give information about the orbital properties of the flyby encounter. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 17 pages, 15 figures, accepted to MNRAS

Informative path planning for scalar dynamic reconstruction using coregionalized Gaussian processes and a spatiotemporal kernel

Authors: Lorenzo Booth, Stefano Carpin

Abstract: The proliferation of unmanned vehicles offers many opportunities for solving environmental sampling tasks with applications in resource monitoring and precision agriculture. Informative path planning (IPP) includes a family of methods which offer improvements over traditional surveying techniques for suggesting locations for observation collection. In this work, we present a novel solution to the… ▽ More The proliferation of unmanned vehicles offers many opportunities for solving environmental sampling tasks with applications in resource monitoring and precision agriculture. Informative path planning (IPP) includes a family of methods which offer improvements over traditional surveying techniques for suggesting locations for observation collection. In this work, we present a novel solution to the IPP problem by using a coregionalized Gaussian processes to estimate a dynamic scalar field that varies in space and time. Our method improves previous approaches by using a composite kernel accounting for spatiotemporal correlations and at the same time, can be readily incorporated in existing IPP algorithms. Through extensive simulations, we show that our novel modeling approach leads to more accurate estimations when compared with formerly proposed methods that do not account for the temporal dimension. △ Less

Submitted 13 September, 2023; originally announced September 2023.

Comments: Accepted to IROS 2023

MAPS: Constraining Serendipitous Time Variability in Protoplanetary Disk Molecular Ion Emission

Authors: Abygail R. Waggoner, L. Ilsedore Cleeves, Ryan A. Loomis, Yuri Aikawa, Jaehan Bae, Jennifer B. Bergner, Alice S. Booth, Jenny K. Calahan, Gianni Cataldi, Charles J. Law, Romane Le Gal, Feng Long, Karin I. Öberg, Richard Teague, David J. Wilner

Abstract: Theoretical models and observations suggest that the abundances of molecular ions in protoplanetary disks should be highly sensitive to the variable ionization conditions set by the young central star. We present a search for temporal flux variability of HCO+ J=1-0, which was observed as a part of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We split out and imaged t… ▽ More Theoretical models and observations suggest that the abundances of molecular ions in protoplanetary disks should be highly sensitive to the variable ionization conditions set by the young central star. We present a search for temporal flux variability of HCO+ J=1-0, which was observed as a part of the Molecules with ALMA at Planet-forming Scales (MAPS) ALMA Large Program. We split out and imaged the line and continuum data for each individual day the five sources were observed (HD 163296, AS 209, GM Aur, MWC 480, and IM Lup, with between 3 to 6 unique visits per source). Significant enhancement (>3σ) was not observed, but we find variations in the spectral profiles in all five disks. Variations in AS 209, GM Aur, and HD 163296 are tentatively attributed to variations in HCO+ flux, while variations in IM Lup and MWC 480 are most likely introduced by differences in the \textit{uv} coverage, which impact the amount of recovered flux during imaging. The tentative detections and low degree of variability are consistent with expectations of X-ray flare driven HCO+ variability, which requires relatively large flares to enhance the HCO+ rotational emission at significant (>20%) levels. These findings also demonstrate the need for dedicated monitoring campaigns with high signal to noise ratios to fully characterize X-ray flare driven chemistry. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: Accepted for publication in ApJ, 18 pages, 9 figures

Tracing snowlines and C/O ratio in a planet-hosting disk: ALMA molecular line observations towards the HD169142 disk

Authors: Alice S. Booth, Charles J. Law, Milou Temmink, Margot Leemker, Enrique Macias

Abstract: The composition of a forming planet is set by the material it accretes from its parent protoplanetary disk. Therefore, it is crucial to map the chemical make-up of the gas in disks to understand the chemical environment of planet formation. This paper presents molecular line observations taken with the Atacama Large Millimeter/submillimeter Array of the planet-hosting disk around the young star HD… ▽ More The composition of a forming planet is set by the material it accretes from its parent protoplanetary disk. Therefore, it is crucial to map the chemical make-up of the gas in disks to understand the chemical environment of planet formation. This paper presents molecular line observations taken with the Atacama Large Millimeter/submillimeter Array of the planet-hosting disk around the young star HD 169142. We detect N2H+, CH3OH, [CI], DCN, CS, C34S, 13CS, H2CS, H2CO, HC3N and c-C3H2 in this system for the first time. Combining these data with the recent detection of SO and previously published DCO+ data, we estimate the location of H2O and CO snowlines and investigate radial variations in the gas phase C/O ratio. We find that the HD 169142 disk has a relatively low N2H+ flux compared to the disks around Herbig stars HD 163296 and MWC 480 indicating less CO freeze-out and place the CO snowline beyond the millimetre disk at ~150 au. The detection of CH3OH from the inner disk is consistent with the H2O snowline being located at the edge of the central dust cavity at ~20 au. The radially varying CS/SO ratio across the proposed H2O snowline location is consistent with this interpretation. Additionally, the detection of CH3OH in such a warm disk adds to the growing evidence supporting the inheritance of complex ices in disks from the earlier, colder stages of star formation. Finally, we propose that the giant HD 169142 b located at 37 au is forming between the CO2 and H2O snowlines where the local elemental make of the gas is expected to have C/O=1.0. △ Less

Submitted 15 August, 2023; originally announced August 2023.

Comments: Accepted A&A 13th August 2023

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

Contour Location for Reliability in Airfoil Simulation Experiments using Deep Gaussian Processes

Authors: Annie S. Booth, S. Ashwin Renganathan, Robert B. Gramacy

Abstract: Bayesian deep Gaussian processes (DGPs) outperform ordinary GPs as surrogate models of complex computer experiments when response surface dynamics are non-stationary, which is especially prevalent in aerospace simulations. Yet DGP surrogates have not been deployed for the canonical downstream task in that setting: reliability analysis through contour location (CL). In that context, we are motivate… ▽ More Bayesian deep Gaussian processes (DGPs) outperform ordinary GPs as surrogate models of complex computer experiments when response surface dynamics are non-stationary, which is especially prevalent in aerospace simulations. Yet DGP surrogates have not been deployed for the canonical downstream task in that setting: reliability analysis through contour location (CL). In that context, we are motivated by a simulation of an RAE-2822 transonic airfoil which demarcates efficient and inefficient flight conditions. Level sets separating passable versus failable operating conditions are best learned through strategic sequential designs. There are two limitations to modern CL methodology which hinder DGP integration in this setting. First, derivative-based optimization underlying acquisition functions is thwarted by sampling-based Bayesian (i.e., MCMC) inference, which is essential for DGP posterior integration. Second, canonical acquisition criteria, such as entropy, are famously myopic to the extent that optimization may even be undesirable. Here we tackle both of these limitations at once, proposing a hybrid criterion that explores along the Pareto front of entropy and (predictive) uncertainty, requiring evaluation only at strategically located "triangulation" candidates. We showcase DGP CL performance in several synthetic benchmark exercises and on the RAE-2822 airfoil. △ Less

Submitted 9 August, 2024; v1 submitted 8 August, 2023; originally announced August 2023.

Comments: 19 pages, 11 figures

SO and SiS Emission Tracing an Embedded Planet and Compact $^{12}$CO and $^{13}$CO Counterparts in the HD 169142 Disk

Authors: Charles J. Law, Alice S. Booth, Karin I. Öberg

Abstract: Planets form in dusty, gas-rich disks around young stars, while at the same time, the planet formation process alters the physical and chemical structure of the disk itself. Embedded planets will locally heat the disk and sublimate volatile-rich ices, or in extreme cases, result in shocks that sputter heavy atoms such as Si from dust grains. This should cause chemical asymmetries detectable in mol… ▽ More Planets form in dusty, gas-rich disks around young stars, while at the same time, the planet formation process alters the physical and chemical structure of the disk itself. Embedded planets will locally heat the disk and sublimate volatile-rich ices, or in extreme cases, result in shocks that sputter heavy atoms such as Si from dust grains. This should cause chemical asymmetries detectable in molecular gas observations. Using high-angular-resolution ALMA archival data of the HD 169142 disk, we identify compact SO J=8$_8$-7$_7$ and SiS J=19-18 emission coincident with the position of a ${\sim}$2 M$_{\rm{Jup}}$ planet seen as a localized, Keplerian NIR feature within a gas-depleted, annular dust gap at ${\approx}$38 au. The SiS emission is located along an azimuthal arc and has a similar morphology as a known $^{12}$CO kinematic excess. This is the first tentative detection of SiS emission in a protoplanetary disk and suggests that the planet is driving sufficiently strong shocks to produce gas-phase SiS. We also report the discovery of compact $^{12}$CO and $^{13}$CO J=3-2 emission coincident with the planet location. Taken together, a planet-driven outflow provides the best explanation for the properties of the observed chemical asymmetries. We also resolve a bright, azimuthally-asymmetric SO ring at ${\approx}$24 au. While most of this SO emission originates from ice sublimation, its asymmetric distribution implies azimuthal temperature variations driven by a misaligned inner disk or planet-disk interactions. Overall, the HD 169142 disk shows several distinct chemical signatures related to giant planet formation and presents a powerful template for future searches of planet-related chemical asymmetries in protoplanetary disks. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Comments: 22 pages, 12 figures, accepted for publication in ApJL

Deprojecting and constraining the vertical thickness of exoKuiper belts

Authors: James Terrill, Sebastian Marino, Richard A. Booth, Yinuo Han, Jeff Jennings, Mark C. Wyatt

Abstract: Constraining the vertical and radial structure of debris discs is crucial to understanding their formation, evolution and dynamics. To measure both the radial and vertical structure, a disc must be sufficiently inclined. However, if a disc is too close to edge-on, deprojecting its emission becomes non-trivial. In this paper we show how Frankenstein, a non-parametric tool to extract the radial brig… ▽ More Constraining the vertical and radial structure of debris discs is crucial to understanding their formation, evolution and dynamics. To measure both the radial and vertical structure, a disc must be sufficiently inclined. However, if a disc is too close to edge-on, deprojecting its emission becomes non-trivial. In this paper we show how Frankenstein, a non-parametric tool to extract the radial brightness profile of circumstellar discs, can be used to deproject their emission at any inclination as long as they are optically thin and axisymmetric. Furthermore, we extend Frankenstein to account for the vertical thickness of an optically thin disc ($H(r)$) and show how it can be constrained by sampling its posterior probability distribution and assuming a functional form (e.g. constant $h=H/r$), while fitting the radial profile non-parametrically. We use this new method to determine the radial and vertical structure of 16 highly inclined debris discs observed by ALMA. We find a wide range of vertical aspect ratios, $h$, ranging from $0.020\pm0.002$ (AU Mic) to $0.20\pm0.03$ (HD 110058), which are consistent with parametric models. We find a tentative correlation between $h$ and the disc fractional width, as expected if wide discs were more stirred. Assuming discs are self-stirred, the thinnest discs would require the presence of at least 500 km-sized planetesimals. The thickest discs would likely require the presence of planets. We also recover previously inferred and new radial structures, including a potential gap in the radial distribution of HD 61005. Finally, our new extension of Frankenstein also allows constraining how $h$ varies as a function of radius, which we test on 49 Ceti, finding that $h$ is consistent with being constant. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS. 17 pages. 16 figures

Investigating the asymmetric chemistry in the disk around the young star HD 142527

Authors: Milou Temmink, Alice S. Booth, Nienke van der Marel, Ewine F. van Dishoeck

Abstract: The atmospheric composition of planets is determined by the chemistry of the disks in which they form. Studying the gas-phase molecular composition of disks thus allows us to infer what the atmospheric composition of forming planets might be. Recent observations of the IRS 48 disk have shown that (asymmetric) dust traps can directly impact the observable chemistry, through radial and vertical tran… ▽ More The atmospheric composition of planets is determined by the chemistry of the disks in which they form. Studying the gas-phase molecular composition of disks thus allows us to infer what the atmospheric composition of forming planets might be. Recent observations of the IRS 48 disk have shown that (asymmetric) dust traps can directly impact the observable chemistry, through radial and vertical transport, and the sublimation of ices. The asymmetric HD 142527 disk provides another good opportunity to investigate the role of dust traps in setting the disk's chemical composition. In this work, we use archival ALMA observations of the HD 142527 disk to obtain an as large as possible molecular inventory, which allows us to investigate the possible influence of the asymmetric dust trap on the disk's chemistry. We present the first ALMA detections of [C I], 13C18O, DCO+, H2CO and additional transition of HCO+ and CS in this disk. In addition, we have acquired upper limits for non-detected species such as SO and CH3OH. For the majority of the observed molecules, a decrement in the emission at the location of the dust trap is found. For the main CO isotopologues continuum over-subtraction likely causes the observed asymmetry, while for CS and HCN we propose that the observed asymmetries are likely due to shadows cast by the misaligned inner disk. As the emission of the observed molecules is not co-spatial with the dust trap and no SO or CH3OH are found, thermal sublimation of icy mantles does not appear to play a major role in changing the gas-phase composition of the outer disk in HD 142527 disk. Using our observations of 13C18O and DCO+ and a RADMC-3D model, we determine the CO snowline to be located beyond the dust traps, favouring cold gas-phase formation of H2CO, rather than the hydrogenation of CO-ice and subsequent sublimation. △ Less

Submitted 13 April, 2023; originally announced April 2023.

Comments: Accepted for publication in A&A on 12/04/2023

Impact of cross-section uncertainties on supernova neutrino spectral parameter fitting in the Deep Underground Neutrino Experiment

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, P. Amedo, J. Anderson, D. A. Andrade , et al. (1294 additional authors not shown)

Abstract: A primary goal of the upcoming Deep Underground Neutrino Experiment (DUNE) is to measure the $\mathcal{O}(10)$ MeV neutrinos produced by a Galactic core-collapse supernova if one should occur during the lifetime of the experiment. The liquid-argon-based detectors planned for DUNE are expected to be uniquely sensitive to the $ν_e$ component of the supernova flux, enabling a wide variety of physics… ▽ More A primary goal of the upcoming Deep Underground Neutrino Experiment (DUNE) is to measure the $\mathcal{O}(10)$ MeV neutrinos produced by a Galactic core-collapse supernova if one should occur during the lifetime of the experiment. The liquid-argon-based detectors planned for DUNE are expected to be uniquely sensitive to the $ν_e$ component of the supernova flux, enabling a wide variety of physics and astrophysics measurements. A key requirement for a correct interpretation of these measurements is a good understanding of the energy-dependent total cross section $σ(E_ν)$ for charged-current $ν_e$ absorption on argon. In the context of a simulated extraction of supernova $ν_e$ spectral parameters from a toy analysis, we investigate the impact of $σ(E_ν)$ modeling uncertainties on DUNE's supernova neutrino physics sensitivity for the first time. We find that the currently large theoretical uncertainties on $σ(E_ν)$ must be substantially reduced before the $ν_e$ flux parameters can be extracted reliably: in the absence of external constraints, a measurement of the integrated neutrino luminosity with less than 10\% bias with DUNE requires $σ(E_ν)$ to be known to about 5%. The neutrino spectral shape parameters can be known to better than 10% for a 20% uncertainty on the cross-section scale, although they will be sensitive to uncertainties on the shape of $σ(E_ν)$. A direct measurement of low-energy $ν_e$-argon scattering would be invaluable for improving the theoretical precision to the needed level. △ Less

Submitted 7 July, 2023; v1 submitted 29 March, 2023; originally announced March 2023.

Comments: 25 pages, 21 figures

Report number: FERMILAB-PUB-23-132-CSAID-LBNF-ND-T

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

Azimuthal C/O Variations in a Planet-Forming Disk

Authors: Luke Keyte, Mihkel Kama, Alice S. Booth, Edwin A. Bergin, L. Ilsedore Cleeves, Ewine F. van Dishoeck, Maria N. Drozdovskaya, Kenji Furuya, Jonathan Rawlings, Oliver Shorttle, Catherine Walsh

Abstract: The elemental carbon-to-oxygen ratio (C/O) in the atmosphere of a giant planet is a promising diagnostic of that planet's formation history in a protoplanetary disk. Alongside efforts in the exoplanet community to measure C/O in planetary atmospheres, observational and theoretical studies of disks are increasingly focused on understanding how the gas-phase C/O varies both with radial location and… ▽ More The elemental carbon-to-oxygen ratio (C/O) in the atmosphere of a giant planet is a promising diagnostic of that planet's formation history in a protoplanetary disk. Alongside efforts in the exoplanet community to measure C/O in planetary atmospheres, observational and theoretical studies of disks are increasingly focused on understanding how the gas-phase C/O varies both with radial location and between disks. This is mostly tied to the icelines of major volatile carriers such as CO and H2O. Using ALMA observations of CS and SO, we have unearthed evidence for an entirely novel type of C/O variation in the protoplanetary disk around HD 100546: an azimuthal variation from a typical, oxygen-dominated ratio (C/O=0.5) to a carbon-dominated ratio (C/O>1.0). We show that the spatial distribution and peculiar line kinematics of both CS and SO molecules can be well-explained by azimuthal variations in the C/O ratio. We propose a shadowing mechanism that could lead to such a chemical dichotomy. Our results imply that tracing the formation history of giant exoplanets using their atmospheric C/O ratios will need to take into account time-dependent azimuthal C/O variations in a planet's accretion zone. △ Less

Submitted 15 March, 2023; originally announced March 2023.

Comments: Accepted in Nature Astronomy

A major asymmetric ice trap in a planet-forming disk IV. Nitric oxide gas and a lack of CN tracing sublimating ices and a C/O ratio $< 1$

Authors: M. Leemker, A. S. Booth, E. F. van Dishoeck, N. van der Marel, B. Tabone, N. F. W. Ligterink, N. G. C. Brunken, M. R. Hogerheijde

Abstract: [Abridged] Most well-resolved disks observed with ALMA show signs of dust traps. These dust traps set the chemical composition of the planet forming material in these disks, as the dust grains with their icy mantles are trapped at specific radii and could deplete the gas and dust of volatiles at smaller radii. In this work we analyse the first detection of nitric oxide (NO) in a protoplanetary dis… ▽ More [Abridged] Most well-resolved disks observed with ALMA show signs of dust traps. These dust traps set the chemical composition of the planet forming material in these disks, as the dust grains with their icy mantles are trapped at specific radii and could deplete the gas and dust of volatiles at smaller radii. In this work we analyse the first detection of nitric oxide (NO) in a protoplanetary disk. We aim to constrain the nitrogen chemistry and the gas-phase C/O ratio in the highly asymmetric dust trap in the Oph-IRS 48 disk. We use ALMA observations of NO, CN, C$_2$H, and related molecules and model the effect of the dust trap on the physical and chemical structure using the thermochemical code DALI. Furthermore, we explore how ice sublimation contributes to the observed emission lines. NO is only observed at the location of the dust trap but CN and C$_2$H are not detected in the Oph-IRS 48 disk. This results in an CN/NO column density ratio of $< 0.05$ and thus a low C/O ratio at the location of the dust trap. The main gas-phase formation pathways to NO through OH and NH in the fiducial model predict NO emission that is an order of magnitude lower than is observed. The gaseous NO column density can be increased by factors ranging from 2.8 to 10 when the H$_2$O and NH$_3$ gas abundances are significantly boosted by ice sublimation. However, these models are inconsistent with the upper limits on the H$_2$O and OH column densities derived from observations. We propose that the NO emission in the Oph-IRS 48 disk is closely related to the nitrogen containing ices sublimating in the dust trap. The non-detection of CN constrains the C/O ratio both inside and outside the dust trap to be $< 1$ if all nitrogen initially starts as N$_2$ and $\leq 0.6$, consistent with the Solar value, if (part of) the nitrogen initially starts as N or NH$_3$. △ Less

Submitted 1 March, 2023; originally announced March 2023.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 673, A7 (2023)

Highly-parallelized simulation of a pixelated LArTPC on a GPU

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, C. Alt, A. Alton, R. Alvarez, P. Amedo, J. Anderson , et al. (1282 additional authors not shown)

Abstract: The rapid development of general-purpose computing on graphics processing units (GPGPU) is allowing the implementation of highly-parallelized Monte Carlo simulation chains for particle physics experiments. This technique is particularly suitable for the simulation of a pixelated charge readout for time projection chambers, given the large number of channels that this technology employs. Here we pr… ▽ More The rapid development of general-purpose computing on graphics processing units (GPGPU) is allowing the implementation of highly-parallelized Monte Carlo simulation chains for particle physics experiments. This technique is particularly suitable for the simulation of a pixelated charge readout for time projection chambers, given the large number of channels that this technology employs. Here we present the first implementation of a full microphysical simulator of a liquid argon time projection chamber (LArTPC) equipped with light readout and pixelated charge readout, developed for the DUNE Near Detector. The software is implemented with an end-to-end set of GPU-optimized algorithms. The algorithms have been written in Python and translated into CUDA kernels using Numba, a just-in-time compiler for a subset of Python and NumPy instructions. The GPU implementation achieves a speed up of four orders of magnitude compared with the equivalent CPU version. The simulation of the current induced on $10^3$ pixels takes around 1 ms on the GPU, compared with approximately 10 s on the CPU. The results of the simulation are compared against data from a pixel-readout LArTPC prototype. △ Less

Submitted 28 February, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 26 pages, 15 figures

Report number: FERMILAB-PUB-22-926-LBNF

First CLAS12 measurement of DVCS beam-spin asymmetries in the extended valence region

Authors: CLAS Collaboration, G. Christiaens, M. Defurne, D. Sokhan, P. Achenbach, Z. Akbar, M. J. Amaryan, H. Atac, H. Avakian, C. Ayerbe Gayoso, L. Baashen, N. A. Baltzell, L. Barion, M. Bashkanov, M. Battaglieri, I. Bedlinskiy, B. Benkel, F. Benmokhtar, A. Bianconi, A. S. Biselli, M. Bondi, W. A. Booth, F. Bossù, S. Boiarinov, K. -Th. Brinkmann , et al. (146 additional authors not shown)

Abstract: Deeply virtual Compton scattering (DVCS) allows one to probe Generalized Parton Distributions (GPDs) describing the 3D structure of the nucleon. We report the first measurement of the DVCS beam-spin asymmetry using the CLAS12 spectrometer with a 10.2 and 10.6 GeV electron beam scattering from unpolarised protons. The results greatly extend the $Q^2$ and Bjorken-$x$ phase space beyond the existing… ▽ More Deeply virtual Compton scattering (DVCS) allows one to probe Generalized Parton Distributions (GPDs) describing the 3D structure of the nucleon. We report the first measurement of the DVCS beam-spin asymmetry using the CLAS12 spectrometer with a 10.2 and 10.6 GeV electron beam scattering from unpolarised protons. The results greatly extend the $Q^2$ and Bjorken-$x$ phase space beyond the existing data in the valence region and provide over 2000 new data points measured with unprecedented statistical uncertainty, setting new, tight constraints for future phenomenological studies. △ Less

Submitted 2 December, 2022; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: Revised Figure 4 and discussion around the number of effective ANNs after Bayesian reweighting

Report number: JLAB-PHY-22-3752

Identification and reconstruction of low-energy electrons in the ProtoDUNE-SP detector

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, C. Alt, A. Alton, R. Alvarez, P. Amedo, J. Anderson , et al. (1235 additional authors not shown)

Abstract: Measurements of electrons from $ν_e$ interactions are crucial for the Deep Underground Neutrino Experiment (DUNE) neutrino oscillation program, as well as searches for physics beyond the standard model, supernova neutrino detection, and solar neutrino measurements. This article describes the selection and reconstruction of low-energy (Michel) electrons in the ProtoDUNE-SP detector. ProtoDUNE-SP is… ▽ More Measurements of electrons from $ν_e$ interactions are crucial for the Deep Underground Neutrino Experiment (DUNE) neutrino oscillation program, as well as searches for physics beyond the standard model, supernova neutrino detection, and solar neutrino measurements. This article describes the selection and reconstruction of low-energy (Michel) electrons in the ProtoDUNE-SP detector. ProtoDUNE-SP is one of the prototypes for the DUNE far detector, built and operated at CERN as a charged particle test beam experiment. A sample of low-energy electrons produced by the decay of cosmic muons is selected with a purity of 95%. This sample is used to calibrate the low-energy electron energy scale with two techniques. An electron energy calibration based on a cosmic ray muon sample uses calibration constants derived from measured and simulated cosmic ray muon events. Another calibration technique makes use of the theoretically well-understood Michel electron energy spectrum to convert reconstructed charge to electron energy. In addition, the effects of detector response to low-energy electron energy scale and its resolution including readout electronics threshold effects are quantified. Finally, the relation between the theoretical and reconstructed low-energy electron energy spectrum is derived and the energy resolution is characterized. The low-energy electron selection presented here accounts for about 75% of the total electron deposited energy. After the addition of lost energy using a Monte Carlo simulation, the energy resolution improves from about 40% to 25% at 50~MeV. These results are used to validate the expected capabilities of the DUNE far detector to reconstruct low-energy electrons. △ Less

Submitted 31 May, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 19 pages, 10 figures

Report number: FERMILAB-PUB-22-784, CERN-EP-DRAFT-MISC-2022-008

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

Sulphur monoxide emission tracing an embedded planet in the HD 100546 protoplanetary disk

Authors: Alice S. Booth, John D. Ilee, Catherine Walsh, Mihkel Kama, Luke Keyte, Ewine F. van Dishoeck, Hideko Nomura

Abstract: Molecular line observations are powerful tracers of the physical and chemical conditions across the different evolutionary stages of star, disk and planet formation. Using the high angular resolution and unprecedented sensitivity of the Atacama Large Millimeter Array (ALMA) there is now a drive to detect small-scale gas structures in protoplanetary disks that can be attributed directly to forming… ▽ More Molecular line observations are powerful tracers of the physical and chemical conditions across the different evolutionary stages of star, disk and planet formation. Using the high angular resolution and unprecedented sensitivity of the Atacama Large Millimeter Array (ALMA) there is now a drive to detect small-scale gas structures in protoplanetary disks that can be attributed directly to forming planets. We report high angular resolution ALMA Band 7 observations of sulphur monoxide (SO) in the nearby planet-hosting disk around Herbig star HD 100546. SO is rarely detected in evolved protoplanetary disks but in other environments, it is most often utilised as a tracer of shocks. The SO emission from the HD 100546 disk is primarily originating from gas within the approx. 20 au mm-dust cavity and shows a clear azimuthal brightness asymmetry of a factor of 2. In addition, we see a significant difference in the line profile shape when comparing these new Cycle 7 data to Cycle 0 data of the same SO transitions. We discuss the different physical/chemical mechanisms that could be responsible for this asymmetry and time variability including disk winds, disk warps, and a shock triggered by a (forming) planet. We propose that the SO is enhanced in the cavity due to the presence of a giant planet. The SO asymmetry complements evidence for hot circumplanetary material around the giant planet HD 100546 c traced via CO ro-vibrational emission. This work sets the stage for further observational and modelling efforts to detect and understand the chemical imprint of a forming planet on its parent disk. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Comments: Accepted A&A 24th October

Journal ref: A&A 669, A53 (2023)

A multidimensional study of the structure function ratio $σ_{LT'}/σ_{0}$ from hard exclusive $π^+$ electro-production off protons in the GPD regime

Authors: S. Diehl, A. Kim, K. Joo, P. Achenbach, Z. Akbar, M. J. Amaryan, H. Atac, H. Avagyan, C. Ayerbe Gayoso, L. Baashen, L. Barion, M. Bashkanov, M. Battaglieri, I. Bedlinskiy, B. Benkel, F. Benmokhtar, A. Bianconi, A. S. Biselli, M. Bondi, W. A. Booth, F. Bossu, S. Boiarinov, K. -Th. Brinkmann, W. J. Briscoe, S. Bueltmann , et al. (129 additional authors not shown)

Abstract: A multidimensional extraction of the structure function ratio $σ_{LT'}/σ_{0}$ from the hard exclusive $\vec{e} p \to e^\prime n π^+$ reaction above the resonance region has been performed. The study was done based on beam-spin asymmetry measurements using a 10.6 GeV incident electron beam on a liquid-hydrogen target and the CLAS12 spectrometer at Jefferson Lab. The measurements focus on the very f… ▽ More A multidimensional extraction of the structure function ratio $σ_{LT'}/σ_{0}$ from the hard exclusive $\vec{e} p \to e^\prime n π^+$ reaction above the resonance region has been performed. The study was done based on beam-spin asymmetry measurements using a 10.6 GeV incident electron beam on a liquid-hydrogen target and the CLAS12 spectrometer at Jefferson Lab. The measurements focus on the very forward regime ($t/Q^{2}$ $\ll$ 1) with a wide kinematic range of $x_{B}$ in the valence regime (0.17 $<$ $x_{B}$ $<$ 0.55), and virtualities $Q^{2}$ ranging from 1.5 GeV$^{2}$ up to 6 GeV$^{2}$. The results and their comparison to theoretical models based on Generalized Parton Distributions demonstrate the sensitivity to chiral-odd GPDs and the directly related tensor charge of the nucleon. In addition, the data is compared to an extension of a Regge formalism at high photon virtualities. It was found that the Regge model provides a better description at low $Q^{2}$, while the GPD model is more appropriate at high $Q^{2}$. △ Less

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

Comments: arXiv admin note: text overlap with arXiv:2007.15677

First Measurement of $Λ$ Electroproduction off Nuclei in the Current and Target Fragmentation Regions

Authors: T. Chetry, L. El Fassi, W. K. Brooks, R. Dupré, A. El Alaoui, K. Hafidi, P. Achenbach, K. P. Adhikari, Z. Akbar, W. R. Armstrong, M. Arratia, H. Atac, H. Avakian, L. Baashen, N. A. Baltzell, L. Barion, M. Bashkanov, M. Battaglieri, I. Bedlinskiy, B. Benkel, F. Benmokhtar, A. Bianconi, A. S. Biselli, M. Bondi, W. A. Booth , et al. (129 additional authors not shown)

Abstract: We report results of $Λ$ hyperon production in semi-inclusive deep-inelastic scattering off deuterium, carbon, iron, and lead targets obtained with the CLAS detector and the Continuous Electron Beam Accelerator Facility 5.014~GeV electron beam. These results represent the first measurements of the $Λ$ multiplicity ratio and transverse momentum broadening as a function of the energy fraction~($z$)… ▽ More We report results of $Λ$ hyperon production in semi-inclusive deep-inelastic scattering off deuterium, carbon, iron, and lead targets obtained with the CLAS detector and the Continuous Electron Beam Accelerator Facility 5.014~GeV electron beam. These results represent the first measurements of the $Λ$ multiplicity ratio and transverse momentum broadening as a function of the energy fraction~($z$) in the current and target fragmentation regions. The multiplicity ratio exhibits a strong suppression at high~$z$~and~an enhancement at~low~$z$. The measured transverse momentum broadening is an order of magnitude greater than that seen for light mesons. This indicates that the propagating entity interacts very strongly with the nuclear medium, which suggests that propagation of diquark configurations in the nuclear medium takes place at least part of the time, even at high~$z$. The trends of these results are qualitatively described by the Giessen Boltzmann-Uehling-Uhlenbeck transport model, particularly for the multiplicity ratios. These observations will potentially open a new era of studies of the structure of the nucleon as well as of strange baryons. △ Less

Submitted 1 April, 2023; v1 submitted 24 October, 2022; originally announced October 2022.

Comments: 14 pages, 6 figures

A New Distance to the Supernova Remnant DA 530 Based on HI Absorption of Polarized Emission

Authors: Rebecca A. Booth, Roland Kothes, Tom Landecker, Jo-Anne Brown, Andrew Gray, Tyler Foster, Eric Greisen

Abstract: Supernova remnants (SNRs) are significant contributors of matter and energy to the interstellar medium. Understanding the impact and the mechanism of this contribution requires knowledge of the physical size, energy, and expansion rate of individual SNRs, which can only come if reliable distances can be obtained. We aim to determine the distance to the SNR DA 530 (G93.3+6.9), an object of low surf… ▽ More Supernova remnants (SNRs) are significant contributors of matter and energy to the interstellar medium. Understanding the impact and the mechanism of this contribution requires knowledge of the physical size, energy, and expansion rate of individual SNRs, which can only come if reliable distances can be obtained. We aim to determine the distance to the SNR DA 530 (G93.3+6.9), an object of low surface brightness. To achieve this, we used the Dominion Radio Astrophysical Observatory Synthesis Telescope and the National Radio Astronomy Observatory Very Large Array to observe the absorption by intervening HI of the polarized emission from DA 530. Significant absorption was detected at velocities $-28$ and -67 km/s (relative to the local standard of rest), corresponding to distances of 4.4 and 8.3 kpc, respectively. Based on the radio and X-ray characteristics of DA 530, we conclude that the minimum distance is 4.4$^{+0.4}_{-0.2}$ kpc. At this minimum distance, the diameter of the SNR is 34$^{+4}_{-1}$ pc, and the elevation above the Galactic plane is 537$^{+40}_{-32}$ pc. The $-67$ km/s absorption likely occurs in gas whose velocity is not determined by Galactic rotation. We present a new data processing method for combining Stokes $Q$ and $U$ observations of the emission from an SNR into a single HI absorption spectrum, which avoids the difficulties of the noise-bias subtraction required for the calculation of polarized intensity. The polarized absorption technique can be applied to determine distances to many more SNRs. △ Less

Submitted 21 October, 2022; originally announced October 2022.

A portable coding strategy to exploit vectorization on combustion simulations

Authors: Fabio Banchelli, Guillermo Oyarzun, Marta Garcia-Gasulla, Filippo Mantovani, Ambrus Both, Guillaume Houzeaux, Daniel Mira

Abstract: The complexity of combustion simulations demands the latest high-performance computing tools to accelerate its time-to-solution results. A current trend on HPC systems is the utilization of CPUs with SIMD or vector extensions to exploit data parallelism. Our work proposes a strategy to improve the automatic vectorization of finite element-based scientific codes. The approach applies a parametric c… ▽ More The complexity of combustion simulations demands the latest high-performance computing tools to accelerate its time-to-solution results. A current trend on HPC systems is the utilization of CPUs with SIMD or vector extensions to exploit data parallelism. Our work proposes a strategy to improve the automatic vectorization of finite element-based scientific codes. The approach applies a parametric configuration to the data structures to help the compiler detect the block of codes that can take advantage of vector computation while maintaining the code portable. A detailed analysis of the computational impact of this methodology on the different stages of a CFD solver is studied on the PRECCINSTA burner simulation. Our parametric implementation has proven to help the compiler generate more vector instructions in the assembly operation: this results in a reduction of up to 9.3 times of the total executed instruction maintaining constant the Instructions Per Cycle and the CPU frequency. The proposed strategy improves the performance of the CFD case under study up to 4.67 times on the MareNostrum 4 supercomputer. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Disentangling protoplanetary disk gas mass and carbon depletion through combined atomic and molecular tracers

Authors: J. A. Sturm, A. S. Booth, M. K. McClure, M. Leemker, E. F. van Dishoeck

Abstract: The total disk gas mass and elemental C, N, O composition of protoplanetary disks are crucial ingredients for our understanding of planet formation. Measuring the gas mass is complicated, since H$_2$ cannot be detected in the cold bulk of the disk and the elemental abundances with respect to hydrogen are degenerate with gas mass in all disk models. We present new NOEMA observations of CO, $^{13}$C… ▽ More The total disk gas mass and elemental C, N, O composition of protoplanetary disks are crucial ingredients for our understanding of planet formation. Measuring the gas mass is complicated, since H$_2$ cannot be detected in the cold bulk of the disk and the elemental abundances with respect to hydrogen are degenerate with gas mass in all disk models. We present new NOEMA observations of CO, $^{13}$CO, C$^{18}$O and optically thin C$^{17}$O $J$=2-1 lines, and use additional high angular resolution Atacama Large Millimeter Array millimeter continuum and CO data to construct a representative model of LkCa 15. The transitions that constrain the gas mass and carbon abundance most are C$^{17}$O 2-1, N${_2}$H$^+$ 3-2 and HD 1-0. Using these three molecules we find that the gas mass in the LkCa 15 disk is $M_\mathrm{g}=0.01 ^{+0.01}_{-0.004} M_{\odot}$, a factor of six lower than estimated before. The carbon abundance is C/H = ($3 \pm 1.5) \times10^{-5}$, implying a moderate depletion of elemental carbon by a factor of 3-9. All other analyzed transitions also agree with these numbers, within a modeling uncertainty of a factor of two. Using the resolved \ce{C2H} image we find a C/O ratio of $\sim$1, which is consistent with literature values of H$_2$O depletion in this disk. The lack of severe carbon depletion in the LkCa 15 disk is consistent with the young age of the disk, but contrasts with the higher depletions seen in older cold transition disks. Combining optically thin CO isotopologue lines with N$_2$H$^+$ is promising to break the degeneracy between gas mass and CO abundance. The moderate level of depletion for this source with a cold, but young disk, suggests that long carbon transformation timescales contribute to the evolutionary trend seen in the level of carbon depletion among disk populations, rather than evolving temperature effects and presence of dust traps alone. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 16 pages, 8 figures. Accepted for publication in A&A

The Profiled Feldman-Cousins technique for confidence interval construction in the presence of nuisance parameters

Authors: M. A. Acero, B. Acharya, P. Adamson, L. Aliaga, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, C. Backhouse, M. Baird, N. Balashov, P. Baldi, B. A. Bambah, S. Bashar, A. Bat, K. Bays, R. Bernstein, V. Bhatnagar, D. Bhattarai, B. Bhuyan, J. Bian, A. C. Booth, R. Bowles , et al. (196 additional authors not shown)

Abstract: Measuring observables to constrain models using maximum-likelihood estimation is fundamental to many physics experiments. The Profiled Feldman-Cousins method described here is a potential solution to common challenges faced in constructing accurate confidence intervals: small datasets, bounded parameters, and the need to properly handle nuisance parameters. This method achieves more accurate frequ… ▽ More Measuring observables to constrain models using maximum-likelihood estimation is fundamental to many physics experiments. The Profiled Feldman-Cousins method described here is a potential solution to common challenges faced in constructing accurate confidence intervals: small datasets, bounded parameters, and the need to properly handle nuisance parameters. This method achieves more accurate frequentist coverage than other methods in use, and is generally applicable to the problem of parameter estimation in neutrino oscillations and similar measurements. We describe an implementation of this method in the context of the NOvA experiment. △ Less

Submitted 1 August, 2022; v1 submitted 28 July, 2022; originally announced July 2022.

Comments: 19 pages, 12 figures

Report number: FERMILAB-PUB-22-476-ND