Showing 1–50 of 257 results for author: Clarke, J

Geophysical Observations of the 24 September 2023 OSIRIS-REx Sample Return Capsule Re-Entry

Authors: Elizabeth A. Silber, Daniel C. Bowman, Chris G. Carr, David P. Eisenberg, Brian R. Elbing, Benjamin Fernando, Milton A. Garcés, Robert Haaser, Siddharth Krishnamoorthy, Charles A. Langston, Yasuhiro Nishikawa, Jeremy Webster, Jacob F. Anderson, Stephen Arrowsmith, Sonia Bazargan, Luke Beardslee, Brant Beck, Jordan W. Bishop, Philip Blom, Grant Bracht, David L. Chichester, Anthony Christe, Kenneth Cummins, James Cutts, Lisa Danielson , et al. (57 additional authors not shown)

Abstract: Sample Return Capsules (SRCs) entering Earth's atmosphere at hypervelocity from interplanetary space are a valuable resource for studying meteor phenomena. The 24 September 2023 arrival of the OSIRIS-REx (Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer) SRC provided an unprecedented chance for geophysical observations of a well-characterized source with kn… ▽ More Sample Return Capsules (SRCs) entering Earth's atmosphere at hypervelocity from interplanetary space are a valuable resource for studying meteor phenomena. The 24 September 2023 arrival of the OSIRIS-REx (Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer) SRC provided an unprecedented chance for geophysical observations of a well-characterized source with known parameters, including timing and trajectory. A collaborative effort involving researchers from 16 institutions executed a carefully planned geophysical observational campaign at strategically chosen locations, deploying over 400 ground-based sensors encompassing infrasound, seismic, distributed acoustic sensing (DAS), and GPS technologies. Additionally, balloons equipped with infrasound sensors were launched to capture signals at higher altitudes. This campaign (the largest of its kind so far) yielded a wealth of invaluable data anticipated to fuel scientific inquiry for years to come. The success of the observational campaign is evidenced by the near-universal detection of signals across instruments, both proximal and distal. This paper presents a comprehensive overview of the collective scientific effort, field deployment, and preliminary findings. The early findings have the potential to inform future space missions and terrestrial campaigns, contributing to our understanding of meteoroid interactions with planetary atmospheres. Furthermore, the dataset collected during this campaign will improve entry and propagation models as well as augment the study of atmospheric dynamics and shock phenomena generated by meteoroids and similar sources. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 87 pages, 14 figures

Inner walls or vortices? Crescent-shaped asymmetries in ALMA observations of protoplanetary discs

Authors: Álvaro Ribas, Cathie J. Clarke, Francesco Zagaria

Abstract: Crescent-shaped asymmetries are common in millimetre observations of protoplanetary discs and are usually attributed to vortices or dust overdensities. However, they often appear on a single side of the major axis and roughly symmetric about the minor axis, suggesting a geometric origin. In this work, we interpret such asymmetries as emission from the exposed inner cavity walls of inclined discs a… ▽ More Crescent-shaped asymmetries are common in millimetre observations of protoplanetary discs and are usually attributed to vortices or dust overdensities. However, they often appear on a single side of the major axis and roughly symmetric about the minor axis, suggesting a geometric origin. In this work, we interpret such asymmetries as emission from the exposed inner cavity walls of inclined discs and use them to characterise their vertical extent. Here we focus on the discs around CIDA 9 and RY Tau, first modelling their observations in visibility space with a simple geometric prescription for the walls, and then exploring more detailed radiative transfer models. Accounting for the wall emission yields significantly better residuals than purely axisymmetric models, and we estimate the dust scale height of these systems to be 0.4 au at 37 au for CIDA 9 and 0.2 au at 12 au for RY Tau. Finally, we identify crescent-shaped asymmetries in twelve discs, nine of which have constraints on their orientation - in all cases, the asymmetry appears on the far-side of the disc, lending support to the hypothesis that they are due to their inner rims. Modelling this effect in larger samples of discs will help to build a statistical view of their vertical structure. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: Accepted for publication in MNRAS

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

Angular momentum transport via gravitational instability in the Elias 2-27 disc

Authors: Cristiano Longarini, Giuseppe Lodato, Cathie J. Clarke, Jessica Speedie, Teresa Paneque-Carreno, Edoardo Arrigoni, Pietro Curone, Claudia Toci, Cassandra Hall

Abstract: Gravitational instability is thought to be one of the main drivers of angular momentum transport in young protoplanetary discs. The disc around Elias 2-27 offers a unique example of gravitational instability at work. It is young and massive, displaying two prominent spiral arms in dust continuum emission and global non-axisymmetric kinematic signatures in molecular line data. In this work, we used… ▽ More Gravitational instability is thought to be one of the main drivers of angular momentum transport in young protoplanetary discs. The disc around Elias 2-27 offers a unique example of gravitational instability at work. It is young and massive, displaying two prominent spiral arms in dust continuum emission and global non-axisymmetric kinematic signatures in molecular line data. In this work, we used archival ALMA observations of $^{13}$CO line emission to measure the efficiency of angular momentum transport in the Elias 2-27 system through the kinematic signatures generated by gravitational instability, known as 'GI wiggles'. Assuming the angular momentum is transported by the observed spiral structure and leveraging previously-derived dynamical disc mass measurements, the amount of angular momentum transport we found corresponds to an $α-$viscosity of $α=0.038\pm0.018$. This value implies an accretion rate onto the central star of $\log_{10}\dot{M}_\star=-6.99\pm0.17\text{M}_\odot/\text{yr, which}$ reproduces the one observed value of $\log_{10}\dot{M}_{\star,\text{obs}}=-7.2\pm0.5\text{M}_\odot/\text{yr }$ very well. The excellent agreement we have found serves as further proof that gravitational instability is the main driver of angular momentum transport acting in this system. △ Less

Submitted 11 June, 2024; v1 submitted 9 June, 2024; originally announced June 2024.

Comments: Accepted for publication in A&A letters, 6 pages, 3 images

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

[CII] luminosity models and large-scale image cubes based on COSMOS 2020 and ALPINE-ALMA [CII] data back to the epoch of reionisation

Authors: J. Clarke, C. Karoumpis, D. Riechers, B. Magnelli, Y. Okada, A. Dev, T. Nikola, F. Bertoldi

Abstract: We have implemented a novel method to create simulated [CII] emission line intensity mapping (LIM) data cubes using COSMOS 2020 galaxy catalogue data. It allows us to provide solid lower limits for previous simulation-based model predictions and the expected signal strength of upcoming surveys. We applied [CII]158$μ$m luminosity models to COSMOS 2020 to create LIM cubes covering a 1.2$\times$1.2 d… ▽ More We have implemented a novel method to create simulated [CII] emission line intensity mapping (LIM) data cubes using COSMOS 2020 galaxy catalogue data. It allows us to provide solid lower limits for previous simulation-based model predictions and the expected signal strength of upcoming surveys. We applied [CII]158$μ$m luminosity models to COSMOS 2020 to create LIM cubes covering a 1.2$\times$1.2 deg$^2$ sky area. These models were derived using galaxy bulk property data from the ALPINE-ALMA survey over the redshift range $4.4<z<5.9$, while additional models were taken from the literature. The cubes cover $3.42<z<3.87$, $4.14<z<4.76$, $5.34<z<6.31$, and $6.75<z<8.27$, matched to planned observations from the EoR-Spec module of the Prime-Cam instrument in the Fred Young Submillimeter Telescope (FYST). We also created predictions including additional galaxies below current detection limits by `extrapolating' from the faint end of the COSMOS 2020 luminosity function, comparing these to predictions from the literature. In addition, we computed signal-to-noise (S/N) ratios for the power spectra, using parameters from the planned FYST survey with predicted instrumental noise levels. We find lower limits for the expected power spectrum using the likely incomplete empirical data: when normalised by $2π^2$, the amplitudes at $k$ = 1 Mpc$^{-1}$ are 3.06$\times$10$^7$, 1.43$\times$10$^7$, 9.80$\times$10$^5$, 2.77$\times$10$^5$ (Jy sr$^{-1}$)$^2$ for the aforementioned redshift ranges. For the extrapolated sample, the power spectra are consistent with prior predictions, indicating that extrapolation is a viable method for creating mock LIM cubes. In this case, we expect a result of S/N>1 when using FYST parameters. However, our high-redshift results remain inconclusive because of the poor completeness of COSMOS 2020 at $z$>6.3. These predictions will be improved on the basis of future JWST data. △ Less

Submitted 5 July, 2024; v1 submitted 8 April, 2024; originally announced April 2024.

Comments: 28 pages, 22 figures, accepted by A&A

Modeling JWST MIRI-MRS Observations of T Cha: Mid-IR Noble Gas Emission Tracing a Dense Disk Wind

Authors: Andrew D. Sellek, Naman S. Bajaj, Ilaria Pascucci, Cathie J. Clarke, Richard Alexander, Chengyan Xie, Giulia Ballabio, Dingshan Deng, Uma Gorti, Andras Gaspar, Jane Morrison

Abstract: [Ne II] 12.81 $μ\mathrm{m}$ emission is a well-used tracer of protoplanetary disk winds due to its blueshifted line profile. MIRI-MRS recently observed T Cha, detecting this line along with lines of [Ne III], [Ar II] and [Ar III], with the [Ne II] and [Ne III] lines found to be extended while the [Ar II] was not. In this complementary work, we use these lines to address long-debated questions abou… ▽ More [Ne II] 12.81 $μ\mathrm{m}$ emission is a well-used tracer of protoplanetary disk winds due to its blueshifted line profile. MIRI-MRS recently observed T Cha, detecting this line along with lines of [Ne III], [Ar II] and [Ar III], with the [Ne II] and [Ne III] lines found to be extended while the [Ar II] was not. In this complementary work, we use these lines to address long-debated questions about protoplanetary disk winds regarding their mass-loss rate, the origin of their ionization, and the role of magnetically-driven winds as opposed to photoevaporation. To this end, we perform photoionization radiative transfer on simple hydrodynamic wind models to map the line emission. We compare the integrated model luminosities to those observed with MIRI-MRS to identify which models most closely reproduce the data and produce synthetic images from these to understand what information is captured by measurements of the line extents. Along with the low degree of ionization implied by the line ratios, the relative compactness of [Ar II] compared to [Ne II] is particularly constraining. This requires Ne II production by hard X-rays and Ar II production by soft X-rays (and/or EUV) in an extended ($\gtrsim 10$ au) wind that is shielded from soft X-rays - necessitating a dense wind with material launched on scales down to ~1 au. Such conditions could be produced by photoevaporation, whereas an extended MHD wind producing equal shielding would likely underpredict the line fluxes. However, a tenuous inner MHD wind may still contribute to shielding the extended wind. This picture is consistent with constraints from spectrally-resolved line profiles. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: 32 pages, 16 figures, Accepted 14/03/24 to the Astronomical Journal. Complementary modeling to Bajaj et al. 2024 (arXiv:2403.01060)

Axion Dark Matter eXperiment: Run 1A Analysis Details

Authors: C. Boutan, B. H. LaRoque, E. Lentz, N. S. Oblath, M. S. Taubman, J. Tedeschi, J. Yang, A. M. Jones, T. Braine, N. Crisosto, L. J Rosenberg, G. Rybka, D. Will, D. Zhang, S. Kimes, R. Ottens, C. Bartram, D. Bowring, R. Cervantes, A. S. Chou, S. Knirck, D. V. Mitchell, A. Sonnenschein, W. Wester, R. Khatiwada , et al. (28 additional authors not shown)

Abstract: The ADMX collaboration gathered data for its Run 1A axion dark matter search from January to June 2017, scanning with an axion haloscope over the frequency range 645-680 MHz (2.66-2.81 ueV in axion mass) at DFSZ sensitivity. The resulting axion search found no axion-like signals comprising all the dark matter in the form of a virialized galactic halo over the entire frequency range, implying lower… ▽ More The ADMX collaboration gathered data for its Run 1A axion dark matter search from January to June 2017, scanning with an axion haloscope over the frequency range 645-680 MHz (2.66-2.81 ueV in axion mass) at DFSZ sensitivity. The resulting axion search found no axion-like signals comprising all the dark matter in the form of a virialized galactic halo over the entire frequency range, implying lower bound exclusion limits at or below DFSZ coupling at the 90% confidence level. This paper presents expanded details of the axion search analysis of Run 1A, including review of relevant experimental systems, data-taking operations, preparation and interpretation of raw data, axion search methodology, candidate handling, and final axion limits. △ Less

Submitted 27 December, 2023; originally announced December 2023.

Comments: 27 pages, 19 figures, accepted for publication in PRD

EWOCS-I: The catalog of X-ray sources in Westerlund 1 from the Extended Westerlund 1 and 2 Open Clusters Survey

Authors: M. G. Guarcello, E. Flaccomio, J. F. Albacete-Colombo, V. Almendros-Abad, K. Anastasopoulou, M. Andersen, C. Argiroffi, A. Bayo, E. S. Bartlett, N. Bastian, M. De Becker, W. Best, R. Bonito, A. Borghese, D. Calzetti, R. Castellanos, C. Cecchi-Pestellini, S. Clark, C. J. Clarke, F. Coti Zelati, F. Damiani, J. J. Drake, M. Gennaro, A. Ginsburg, E. K. Grebel , et al. (26 additional authors not shown)

Abstract: Context. With a mass exceeding several 10^4 solar masses and a rich and dense population of massive stars, supermassive young star clusters represent the most massive star-forming environment that is dominated by the feedback from massive stars and gravitational interactions among stars. Aims. In this paper we present the "Extended Westerlund 1 and 2 Open Clusters Survey" (EWOCS) project, which ai… ▽ More Context. With a mass exceeding several 10^4 solar masses and a rich and dense population of massive stars, supermassive young star clusters represent the most massive star-forming environment that is dominated by the feedback from massive stars and gravitational interactions among stars. Aims. In this paper we present the "Extended Westerlund 1 and 2 Open Clusters Survey" (EWOCS) project, which aims to investigate the influence of the starburst environment on the formation of stars and planets, and on the evolution of both low and high mass stars. The primary targets of this project are Westerlund 1 and 2, the closest supermassive star clusters to the Sun. Methods. The project is based primarily on recent observations conducted with the Chandra and JWST observatories. Specifically, the Chandra survey of Westerlund 1 consists of 36 new ACIS-I observations, nearly co-pointed, for a total exposure time of 1 Msec. Additionally, we included 8 archival Chandra/ACIS-S observations. This paper presents the resulting catalog of X-ray sources within and around Westerlund 1. Sources were detected by combining various existing methods, and photon extraction and source validation were carried out using the ACIS-Extract software. Results. The EWOCS X-ray catalog comprises 5963 validated sources out of the 9420 initially provided to ACIS-Extract, reaching a photon flux threshold of approximately 2x10^-8 photons/cm^2/s. The X-ray sources exhibit a highly concentrated spatial distribution, with 1075 sources located within the central 1 arcminute. We have successfully detected X-ray emissions from 126 out of the 166 known massive stars of the cluster, and we have collected over 71000 photons from the magnetar CXO J164710.20-455217 △ Less

Submitted 15 December, 2023; v1 submitted 14 December, 2023; originally announced December 2023.

Comments: The paper has been accepted for publication by Astronomy and Astrophysics

The enigmatic abundance of atomic hydrogen in Saturn's upper atmosphere

Authors: Lotfi Ben-Jaffel, Julie Moses, Robert A. West, M-K. aye, Eric T. Bradley, John T. Clarke, Jay B. Holber, Gilda E. Ballester

Abstract: A planet's Lyman-α (Lyα) emission is sensitive to its thermospheric structure. Here, we report joint Hubble Space Telescope (HST) and Cassini cross-calibration observations of the Saturn Lyα emission made two weeks before the Cassini grand finale. To investigate the long-term Saturn Lyα airglow observed by different ultraviolet instruments, we cross-correlate their calibration, finding that while… ▽ More A planet's Lyman-α (Lyα) emission is sensitive to its thermospheric structure. Here, we report joint Hubble Space Telescope (HST) and Cassini cross-calibration observations of the Saturn Lyα emission made two weeks before the Cassini grand finale. To investigate the long-term Saturn Lyα airglow observed by different ultraviolet instruments, we cross-correlate their calibration, finding that while the official Cassini/UVIS sensitivity should be lowered by ~75%, the Voyager 1/UVS sensitivities should be enhanced by ~20% at the Lyα channels. This comparison also allowed us to discover a permanent feature of the Saturn disk Lyα brightness that appears at all longitudes as a brightness excess (Lyα bulge) of ~30% (~12σ) extending over the latitude range ~5-35N compared to the regions at equator and ~60N. This feature is confirmed by three distinct instruments between 1980 & 2017 in the Saturn north hemisphere. To analyze the Lyα observations, we use a radiation transfer (RT) model of resonant scattering of solar and interplanetary Lyα photons, and a latitude-dependent photochemistry model of the upper atmosphere constrained by occultation and remote-sensing observations. For each latitude, we show that the Lyα observations are sensitive to the temperature profile in the upper stratosphere and lower thermosphere, thus providing useful information in a region of the atmosphere that is difficult to probe by other means. In the Saturn Lyα bulge region, at latitudes between ~5 to ~35°, the observed brightening and line broadening support seasonal effects, variation of the temperature vertical profile, and potential superthermal atoms that require confirmation. △ Less

Submitted 23 November, 2023; originally announced November 2023.

Comments: Published in Planetary Science Journal. Added reference to grants

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

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

Non-Virialized Axion Search Sensitive to Doppler Effects in the Milky Way Halo

Authors: C. Bartram, T. Braine, R. Cervantes, N. Crisosto, N. Du, C. Goodman, M. Guzzetti, C. Hanretty, S. Lee, G. Leum, L. J. Rosenberg, G. Rybka, J. Sinnis, D. Zhang, M. H. Awida, D. Bowring, A. S. Chou, M. Hollister, S. Knirck, A. Sonnenschein, W. Wester, R. Khatiwada, J. Brodsky, G. Carosi, L. D. Duffy , et al. (31 additional authors not shown)

Abstract: The Axion Dark Matter eXperiment (ADMX) has previously excluded Dine-Fischler-Srednicki-Zhitnisky (DFSZ) axions between 680-790 MHz under the assumption that the dark matter is described by the isothermal halo model. However, the precise nature of the velocity distribution of dark matter is still unknown, and alternative models have been proposed. We report the results of a non-virialized axion se… ▽ More The Axion Dark Matter eXperiment (ADMX) has previously excluded Dine-Fischler-Srednicki-Zhitnisky (DFSZ) axions between 680-790 MHz under the assumption that the dark matter is described by the isothermal halo model. However, the precise nature of the velocity distribution of dark matter is still unknown, and alternative models have been proposed. We report the results of a non-virialized axion search over the mass range 2.81-3.31 μeV, corresponding to the frequency range 680-800 MHz. This analysis marks the most sensitive search for non-virialized axions sensitive to Doppler effects in the Milky Way Halo to date. Accounting for frequency shifts due to the detector's motion through the Galaxy, we exclude cold flow relic axions with a velocity dispersion of order 10^-7 c with 95% confidence. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Radio multiwavelength analysis of the compact disk CX Tau: Presence of strong free-free variability or anomalous microwave emission

Authors: Pietro Curone, Leonardo Testi, Enrique Macias, Marco Tazzari, Stefano Facchini, Jonathan P. Williams, Cathie J. Clarke, Antonella Natta, Giovanni Rosotti, Claudia Toci, Giuseppe Lodato

Abstract: Protoplanetary disks emit radiation across a broad range of wavelengths, requiring a multiwavelength approach to fully understand their physical mechanisms and how they form planets. Observations at sub-millimeter to centimeter wavelengths can provide insights into the thermal emission from dust, free-free emission from ionized gas, and possible gyro-synchrotron emission from the stellar magnetosp… ▽ More Protoplanetary disks emit radiation across a broad range of wavelengths, requiring a multiwavelength approach to fully understand their physical mechanisms and how they form planets. Observations at sub-millimeter to centimeter wavelengths can provide insights into the thermal emission from dust, free-free emission from ionized gas, and possible gyro-synchrotron emission from the stellar magnetosphere. This work is focused on CX Tau, a ${\sim}0.4\,M_\odot$ star with an extended gas emission and a compact and apparently structureless dust disk, with an average millimeter flux when compared to Class II sources in Taurus. We present observations from the Karl G. Jansky Very Large Array (VLA) across four bands (between 9.0 mm and 6.0 cm) and combine them with archival data from the Atacama Large Millimeter/submillimeter Array (ALMA), the Submillimeter Array (SMA) and the Plateau de Bure Interferometer (PdBI). This multiwavelength approach allows us to separate the dust continuum from other emissions. After isolating the dust thermal emission, we derived an upper limit of the dust disk extent at 1.3 cm which is consistent with theoretical predictions of a radial drift-dominated disk. Centimeter data show a peculiar behavior: deep observations at 6.0 cm did not detect the source, while at 1.3 cm the flux density is anomalously higher than adjacent bands. Intraband spectral indices suggest a dominant contribution from free-free emission, whereas gyro-synchrotron emission is excluded. To explain these observations, we propose a strong variability among the free-free emission with timescales shorter than a month. Another possible interpretation is the presence of anomalous microwave emission from spinning dust grains. △ Less

Submitted 22 August, 2023; v1 submitted 20 July, 2023; originally announced July 2023.

Comments: Accepted for publication in A&A; 13 pages, 13 figures, 1 table

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

Testing protoplanetary disc evolution with CO fluxes. A proof of concept in Lupus and Upper Sco

Authors: Francesco Zagaria, Stefano Facchini, Anna Miotello, Carlo F. Manara, Claudia Toci, Cathie J. Clarke

Abstract: The Atacama Large Millimeter/submillimeter Array (ALMA) revolutionised our understanding of protoplanetary discs. However, the available data have not given conclusive answers yet on the underlying disc evolution mechanisms (viscosity or MHD winds). Improving upon the current results, mostly based on the analysis of disc sizes, is difficult because larger, deeper and higher angular resolution surv… ▽ More The Atacama Large Millimeter/submillimeter Array (ALMA) revolutionised our understanding of protoplanetary discs. However, the available data have not given conclusive answers yet on the underlying disc evolution mechanisms (viscosity or MHD winds). Improving upon the current results, mostly based on the analysis of disc sizes, is difficult because larger, deeper and higher angular resolution surveys would be required, which could be prohibitive even for ALMA. In this Letter, we introduce an alternative method to study disc evolution based on $^{12}$CO fluxes. In fact, fluxes can be readily collected using less time-consuming, lower resolution observations, while tracing the same disc physico-chemical processes as sizes: assuming that $^{12}$CO is optically thick, fluxes scale with the disc surface area. We developed a semi-analytical model to compute $^{12}$CO fluxes and benchmarked it against the results of DALI thermochemical models, recovering an agreement within a factor of three. As a proof of concept we compared our models with Lupus and Upper Sco data, taking advantage of the increased samples, by a factor 1.3 (Lupus) and 3.6 (Upper Sco), when studying fluxes instead of sizes. Models and data agree well only if CO depletion is considered. However, the uncertainties on the initial conditions limited our interpretation of the observations. Our new method can be used to design future ad hoc observational strategies to collect better data and give conclusive answers on disc evolution. △ Less

Submitted 7 April, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Comments: 10 pages, 9 figures; accepted for publication in A&A Letters. Code and data available at https://github.com/fzagaria/COpops.git

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

Accretion of sub-stellar companions as the origin of chemical abundance inhomogeneities in globular clusters

Authors: Andrew J. Winter, Cathie J. Clarke

Abstract: Globular clusters exhibit abundance variations, defining `multiple populations', which have prompted a protracted search for their origin. Properties requiring explanation include: the high fraction of polluted stars ($\sim 40{-}90$~percent, correlated with cluster mass), the absence of pollution in young clusters and the lower pollution rate with binarity and distance from the cluster centre. We… ▽ More Globular clusters exhibit abundance variations, defining `multiple populations', which have prompted a protracted search for their origin. Properties requiring explanation include: the high fraction of polluted stars ($\sim 40{-}90$~percent, correlated with cluster mass), the absence of pollution in young clusters and the lower pollution rate with binarity and distance from the cluster centre. We present a novel mechanism for late delivery of pollutants into stars via accretion of sub-stellar companions. In this scenario, stars move through a medium polluted with AGB and massive star ejecta, accreting material to produce companions with typical mass ratio $q\sim 0.1$. These companions undergo eccentricity excitation due to dynamical perturbations by passing stars, culminating in a merger with their host star. The accretion of the companion alters surface abundances via injected pollutant. Alongside other self-enrichment models, the companion accretion model can explain the dilution of pollutant and correlation with intra-cluster location. The model also explains the ubiquity and discreteness of the populations and correlations of enrichment rates with cluster mass, cluster age and stellar binarity. Abundance variations in some clusters can be broadly reproduced using AGB and massive binary ejecta abundances from the literature. In other clusters, some high companion mass ratios ($q\gtrsim 1$) are required. In these cases, the available mass budget necessitates a variable degree of mixing of the polluted material with the primary star, deviations from model ejecta abundances or mixing of internal burning products. We highlight the avenues of further investigation which are required to explore some of the key processes invoked in this model. △ Less

Submitted 2 March, 2023; v1 submitted 26 January, 2023; originally announced January 2023.

Comments: 29 pages, 20 figures, accepted for publication in MNRAS

A study of convective core overshooting as a function of stellar mass based on two-dimensional hydrodynamical simulations

Authors: I. Baraffe, J. Clarke, A. Morison, D. G. Vlaykov, T. Constantino, T. Goffrey, T. Guillet, A. Le Saux, J. Pratt

Abstract: We perform two-dimensional numerical simulations of core convection for zero-age-main-sequence stars covering a mass range from 3 $M_\odot$ to 20 $M_\odot$. The simulations are performed with the fully compressible time-implicit code MUSIC. We study the efficiency of overshooting, which describes the ballistic process of convective flows crossing a convective boundary, as a function of stellar mas… ▽ More We perform two-dimensional numerical simulations of core convection for zero-age-main-sequence stars covering a mass range from 3 $M_\odot$ to 20 $M_\odot$. The simulations are performed with the fully compressible time-implicit code MUSIC. We study the efficiency of overshooting, which describes the ballistic process of convective flows crossing a convective boundary, as a function of stellar mass and luminosity. We also study the impact of artificially increasing the stellar luminosity for 3 $M_\odot$ models. The simulations cover hundreds to thousands of convective turnover timescales. Applying the framework of extreme plume events previously developed for convective envelopes, we derive overshooting lengths as a function of stellar masses. We find that the overshooting distance ($d_{\rm ov}$) scales with the stellar luminosity ($L$) and the convective core radius ($r_{\rm conv}$). We derive a scaling law $d_{\rm ov} \propto L^{1/3} r_{\rm conv}^{1/2}$ which is implemented in a 1D stellar evolution code and the resulting stellar models are compared to observations. The scaling predicts values for the overshooting distance that significantly increase with stellar mass, in qualitative agreement with observations. Quantitatively, however, the predicted values are underestimated for masses $\gtrsim 10 M_\odot$. Our 2D simulations show the formation of a nearly-adiabatic layer just above the Schwarzschild boundary of the convective core, as exhibited in recent 3D simulations of convection. The most luminous models show a growth in size with time of the nearly-adiabatic layer. This growth seems to slow down as the upper edge of the nearly-adiabatic layer gets closer to the maximum overshooting length and as the simulation time exceeds the typical thermal diffusive timescale in the overshooting layer. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: 13 pages, accepted for publication in MNRAS

Dust dynamics in planet-forming discs in binary systems

Authors: Francesco Zagaria, Giovanni P. Rosotti, Richard D. Alexander, Cathie J. Clarke

Abstract: In multiple stellar systems interactions among the companion stars and their discs affect planet formation. In the circumstellar case tidal truncation makes protoplanetary discs smaller, fainter and less long-lived than those evolving in isolation, thereby reducing the amount of material (gas and dust) available to assemble planetary embryos. On the contrary, in the circumbinary case the reduced a… ▽ More In multiple stellar systems interactions among the companion stars and their discs affect planet formation. In the circumstellar case tidal truncation makes protoplanetary discs smaller, fainter and less long-lived than those evolving in isolation, thereby reducing the amount of material (gas and dust) available to assemble planetary embryos. On the contrary, in the circumbinary case the reduced accretion can increase the disc lifetime, with beneficial effects on planet formation. In this chapter we review the main observational results on discs in multiple stellar systems and discuss their possible explanations, focusing on recent numerical simulations, mainly dealing with dust dynamics and disc evolution. Finally, some open issues and future research directions are examined. △ Less

Submitted 15 December, 2022; originally announced December 2022.

Comments: Accepted for publication in EPJ Plus Focus Point on Environmental and Multiplicity Effects on Planet Formation G. Lodato and C.F. Manara (Guest editors)

Forming short period sub-stellar companions in 47 Tucanae -- II. Analytic expressions for the orbital evolution of planets in dense environments

Authors: Andrew J. Winter, Cathie J. Clarke, Giovanni Rosotti, Mirek Giersz

Abstract: Short period, massive planets, known as hot Jupiters (HJs), have been discovered around $\sim 1$ percent of local field stars. The inward migration necessary to produce HJs may be `low eccentricity', due to torques in the primordial disc, or `high eccentricity' (HEM). The latter involves exciting high orbital eccentricity, allowing sufficiently close passages with the host star to raise circularis… ▽ More Short period, massive planets, known as hot Jupiters (HJs), have been discovered around $\sim 1$ percent of local field stars. The inward migration necessary to produce HJs may be `low eccentricity', due to torques in the primordial disc, or `high eccentricity' (HEM). The latter involves exciting high orbital eccentricity, allowing sufficiently close passages with the host star to raise circularising tides in the planet. We present an analytic framework for quantifying the role of dynamical encounters in high density environments during HEM. We show that encounters can enhance or suppress HEM, depending on the local stellar density and the initial semi-major axis $a_0$. For moderate densities, external perturbations can excite large eccentricities that allow a planet to circularise over the stellar lifetime. At extremely high densities, these perturbations can instead result in tidal disruption of the planet, thus yielding no HJ. This may explain the apparent excess of HJs in M67 compared with their local field star abundance versus their apparent deficit in 47 Tuc. Applying our analytic framework, we demonstrate that for an initial massive planet population similar to the field, the expected HJ occurrence rate in 47 Tuc is $f_\mathrm{HJ}=2.2\times 10^{-3}$, which remains consistent with present constraints. Future large (sample sizes $\gtrsim 10^5$) or sensitive transit surveys of stars in globular clusters are required to refute the hypothesis that the initial planet population is similar to the solar neighbourhood average. Non-detection in such surveys would have broad consequences for planet formation theory, implying planet formation rates in globular clusters must be suppressed across a wide range of $a_0$. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS - 28 pages, 17 figures

Super-resolution trends in the ALMA Taurus survey: Structured inner discs and compact discs

Authors: Jeff Jennings, Marco Tazzari, Cathie J. Clarke, Richard A. Booth, Giovanni P. Rosotti

Abstract: The 1.33 mm survey of protoplanetary discs in the Taurus molecular cloud found annular gaps and rings to be common in extended sources (>~55 au), when their 1D visibility distributions were fit parametrically. We first demonstrate the advantages and limitations of nonparametric visibility fits for data at the survey's 0.12" resolution. Then we use the nonparametric model in Frankenstein ('frank')… ▽ More The 1.33 mm survey of protoplanetary discs in the Taurus molecular cloud found annular gaps and rings to be common in extended sources (>~55 au), when their 1D visibility distributions were fit parametrically. We first demonstrate the advantages and limitations of nonparametric visibility fits for data at the survey's 0.12" resolution. Then we use the nonparametric model in Frankenstein ('frank') to identify new substructure in three compact and seven extended sources. Among the new features we identify three trends: a higher occurrence rate of substructure in the survey's compact discs than previously seen, underresolved (potentially azimuthally asymmetric) substructure in the innermost disc of extended sources, and a 'shoulder' on the trailing edge of a ring in discs with strong depletion at small radii. Noting the shoulder morphology is present in multiple discs observed at higher resolution, we postulate it is tracing a common physical mechanism. We further demonstrate how a super-resolution frank brightness profile is useful in motivating an accurate parametric model, using the highly structured source DL Tau in which frank finds two new rings. Finally we show that sparse (u, v) plane sampling may be masking the presence of substructure in several additional compact survey sources. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: Accepted to MNRAS

Inward and outward migration of massive planets: moving towards a stalling radius

Authors: C. E. Scardoni, C. J. Clarke, G. P. Rosotti, R. A. Booth, R. D. Alexander, E. Ragusa

Abstract: Recent studies on the planet-dominated regime of Type II migration showed that, contrary to the conventional wisdom, massive planets can migrate outwards. Using `fixed-planet' simulations these studies found a correlation between the sign of the torques acting on the planet and the parameter $K'$ (which describes the depth of the gap carved by the planet in the disc). We perform `live-planet' simu… ▽ More Recent studies on the planet-dominated regime of Type II migration showed that, contrary to the conventional wisdom, massive planets can migrate outwards. Using `fixed-planet' simulations these studies found a correlation between the sign of the torques acting on the planet and the parameter $K'$ (which describes the depth of the gap carved by the planet in the disc). We perform `live-planet' simulations exploring a range of $K'$ and disc mass values to test and extend these results. The excitation of planet eccentricity in live-planet simulations breaks the direct dependence of migration rate (rate of change of semi-major axis) on the torques imposed, an effect that `fixed-planet' simulations cannot treat. By disentangling the contribution to the torque due to the semi-major axis evolution from that due to the eccentricity evolution, we recover the relation between the magnitude and sign of migration and $K'$ and argue that this relation may be better expressed in terms of the related gap depth parameter $K$. We present a toy model in which the sign of planetary migration changes at a limiting value of $K$, through which we explore planets' migration in viscously evolving discs. The existence of the torque reversal shapes the planetary system's architecture by accumulating planets either at the stalling radius or in a band around it (defined by the interplay between the planet migration and the disc evolution). In either case, planets pile up in the area $1-10$ au, disfavouring hot Jupiter formation through Type II migration in the planet-dominated regime. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 16 pages, 14 figures, accepted for publication in MNRAS

Modelling the secular evolution of proto-planetary disc dust sizes -- A comparison between the viscous and magnetic wind case

Authors: Francesco Zagaria, Giovanni P. Rosotti, Cathie J. Clarke, Benoît Tabone

Abstract: For many years proto-planetary discs have been thought to evolve viscously: angular momentum redistribution leads to accretion and outward disc spreading. Recently, the hypothesis that accretion is due, instead, to angular momentum removal by magnetic winds gained new popularity: no disc spreading is expected in this case. In this paper, we run several one-dimensional gas and dust simulations to m… ▽ More For many years proto-planetary discs have been thought to evolve viscously: angular momentum redistribution leads to accretion and outward disc spreading. Recently, the hypothesis that accretion is due, instead, to angular momentum removal by magnetic winds gained new popularity: no disc spreading is expected in this case. In this paper, we run several one-dimensional gas and dust simulations to make predictions on the time-evolution of disc sizes \textit{in the dust} and to assess whether they can be used to understand how discs evolve. We show that viscous and magnetic wind models have very different dust disc radii. In particular, MHD wind models are compact and their sizes either remain constant or decrease with time. On the contrary, discs become larger with time in the viscous case (when $α\gtrsim10^{-3}$). Although current observations lack enough sensitivity to discriminate between these two scenarios, higher-sensitivity surveys could be fruitful to this goal on a $1\,{\rm to}\,10\,{\rm Myr}$ age range. When compared with the available ALMA Band~7 data, both viscous and magnetic wind models are compatible with the observationally-inferred dust radii in Lupus, Chamaeleon~I and Upper Sco. Furthermore, in the drift-dominated regime, the size-luminosity correlation is reproduced in Lupus, both in Band~7 and 3, while in Upper Sco a different slope than in the data is predicted. Sub-structures (potentially undetected) can explain several outliers with large observed sizes. Higher-angular-resolution observations will be helpful to test our predictions in the case of more compact discs, expected in both frameworks, particularly at the age of Upper Sco. △ Less

Submitted 22 May, 2022; originally announced May 2022.

Comments: 19 pages, 10 figures; accepted for publication in MNRAS

The stellar mass distribution of the Milky Way's bar: an analytic model

Authors: Mattia C. Sormani, Ortwin Gerhard, Matthieu Portail, Eugene Vasiliev, Jonathan Clarke

Abstract: We present an analytic model of the stellar mass distribution of the Milky Way bar. The model is obtained by fitting a multi-component parametric density distribution to a made-to-measure N-body model of Portail et al., constructed to match a variety of density and kinematics observational data. The analytic model reproduces in detail the 3D density distribution of the N-body bar including the X-s… ▽ More We present an analytic model of the stellar mass distribution of the Milky Way bar. The model is obtained by fitting a multi-component parametric density distribution to a made-to-measure N-body model of Portail et al., constructed to match a variety of density and kinematics observational data. The analytic model reproduces in detail the 3D density distribution of the N-body bar including the X-shape. The model and the gravitational potential it generates are available as part of the software package AGAMA for galactic dynamics, and can be readily used for orbit integrations, hydrodynamical simulations or other applications. △ Less

Submitted 27 April, 2022; originally announced April 2022.

Comments: Accepted in MNRAS letters

The importance of X-ray frequency in driving photoevaporative winds

Authors: Andrew D. Sellek, Cathie J. Clarke, Barbara Ercolano

Abstract: Photoevaporative winds are a promising mechanism for dispersing protoplanetary discs, but so far theoretical models have been unable to agree on the relative roles that the X-ray, Extreme Ultraviolet or Far Ultraviolet play in driving the winds. This has been attributed to a variety of methodological differences between studies, including their approach to radiative transfer and thermal balance, t… ▽ More Photoevaporative winds are a promising mechanism for dispersing protoplanetary discs, but so far theoretical models have been unable to agree on the relative roles that the X-ray, Extreme Ultraviolet or Far Ultraviolet play in driving the winds. This has been attributed to a variety of methodological differences between studies, including their approach to radiative transfer and thermal balance, the choice of irradiating spectrum employed, and the processes available to cool the gas. We use the \textsc{mocassin} radiative transfer code to simulate wind heating for a variety of spectra on a static density grid taken from simulations of an EUV-driven wind. We explore the impact of choosing a single representative X-ray frequency on their ability to drive a wind by measuring the maximum heated column as a function of photon energy. We demonstrate that for reasonable luminosities and spectra, the most effective energies are at a few $100~\mathrm{eV}$, firmly in the softer regions of the X-ray spectrum, while X-rays with energies $\sim1000~\mathrm{eV}$ interact too weakly with disc gas to provide sufficient heating to drive a wind. We develop a simple model to explain these findings. We argue that further increases in the cooling above our models - for example due to molecular rovibrational lines - may further restrict the heating to the softer energies but are unlikely to prevent X-ray heated winds from launching entirely; increasing the X-ray luminosity has the opposite effect. The various results of photoevaporative wind models should therefore be understood in terms of the choice of irradiating spectrum. △ Less

Submitted 20 April, 2022; originally announced April 2022.

Comments: 20 pages, 11 figures; This article has been accepted for publication in MNRAS Published by Oxford University Press on behalf of the Royal Astronomical Society after minor revisions

An APEX search for carbon emission from NGC 1977 proplyds

Authors: Thomas J. Haworth, Jinyoung S. Kim, Lin Qiao, Andrew J. Winter, Jonathan P. Williams, Cathie J. Clarke, James E. Owen, Stefano Facchini, Megan Ansdell, Mikhel Kama, Giulia Ballabio

Abstract: We used the Atacama Pathfinder Experiment (APEX) telescope to search for CI 1-0 (492.16GHz) emission towards 8 proplyds in NGC 1977, which is an FUV radiation environment two orders of magnitude weaker than that irradiating the Orion Nebular Cluster (ONC) proplyds. CI is expected to enable us to probe the wind launching region of externally photoevaporating discs. Of the 8 targets observed, no 3… ▽ More We used the Atacama Pathfinder Experiment (APEX) telescope to search for CI 1-0 (492.16GHz) emission towards 8 proplyds in NGC 1977, which is an FUV radiation environment two orders of magnitude weaker than that irradiating the Orion Nebular Cluster (ONC) proplyds. CI is expected to enable us to probe the wind launching region of externally photoevaporating discs. Of the 8 targets observed, no 3$σ$ detections of the CI line were made despite reaching sensitivities deeper than the anticipated requirement for detection from prior APEX CI observations of nearby discs and models of external photoevaporation of quite massive discs. By comparing both the proplyd mass loss rates and CI flux constraints with a large grid of external photoevaporation simulations, we determine that the non-detections are in fact fully consistent with the models if the proplyd discs are very low mass. Deeper observations in CI and probes of the disc mass with other tracers (e.g. in the continuum and CO) can test this. If such a test finds higher masses, this would imply carbon depletion in the outer disc, as has been proposed for other discs with surprisingly low CI fluxes, though more massive discs would also be incompatible with models that can explain the observed mass loss rates and CI non-detections. The expected remaining lifetimes of the proplyds are estimated to be similar to those of proplyds in the ONC at 0.1Myr. Rapid destruction of discs is therefore also a feature of common, intermediate UV environments. △ Less

Submitted 8 March, 2022; originally announced March 2022.

Comments: 10 pages. Accepted for publication in MNRAS

Stellar multiplicity affects the correlation between proto-planetary disc masses and accretion rates: binaries explain high-accretors in Upper Sco

Authors: Francesco Zagaria, Cathie J. Clarke, Giovanni P. Rosotti, Carlo F. Manara

Abstract: In recent years a correlation between mass accretion rates onto new-born stars and their proto-planetary disc masses was detected in nearby young star-forming regions. Although such a correlation can be interpreted as due to viscous-diffusion processes in the disc, highly-accreting sources with low disc masses in more evolved regions remain puzzling. In this paper, we hypothesise that the presence… ▽ More In recent years a correlation between mass accretion rates onto new-born stars and their proto-planetary disc masses was detected in nearby young star-forming regions. Although such a correlation can be interpreted as due to viscous-diffusion processes in the disc, highly-accreting sources with low disc masses in more evolved regions remain puzzling. In this paper, we hypothesise that the presence of a stellar companion truncating the disc can explain these outliers. Firstly, we searched the literature for information on stellar multiplicity in Lupus, Chamaeleon~I and Upper Sco, finding that roughly 20 per cent of the discs involved in the correlation are in binaries or higher-order multiple stellar systems. We prove with high statistical significance that at any disc mass these sources have systematically higher accretion rates than those in single-stars, with the bulk of the binary population being clustered around $M_\mathrm{disc}/\dot{M}_\mathrm{acc}\approx0.1\,\mathrm{Myr}$. We then run coupled gas and dust one-dimensional evolutionary models of tidally truncated discs to be compared with the data. We find that these models are able to reproduce well most of the population of observed discs in Lupus and Upper Sco, even though the unknown eccentricity of each binary prevents an object by object comparison. In the latter region, the agreement improves if the grain coagulation efficiency is reduced, as may be expected in discs around close binaries. Finally, we mention that thermal winds and sub-structures can be important in explaining few outlying sources. △ Less

Submitted 3 March, 2022; originally announced March 2022.

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

Observation and origin of non-thermal hard X-rays from Jupiter

Authors: Kaya Mori, Charles Hailey, Gabriel Bridges, Shifra Mandel, Amani Garvin, Brian Grefenstette, William Dunn, Benjamin J. Hord, Graziella Branduardi-Raymont, John Clarke, Caitriona Jackman, Melania Nynka, Licia Ray

Abstract: Electrons accelerated on Earth by a rich variety of wave scattering or stochastic processes generate hard non-thermal X-ray bremsstrahlung up to >~ 1 MeV and power Earth's various types of aurorae. Although Jupiter's magnetic field is an order of magnitude larger than Earth's, space-based telescopes have previously detected X-rays only up to ~7 keV. On the basis of theoretical models of the Jovian… ▽ More Electrons accelerated on Earth by a rich variety of wave scattering or stochastic processes generate hard non-thermal X-ray bremsstrahlung up to >~ 1 MeV and power Earth's various types of aurorae. Although Jupiter's magnetic field is an order of magnitude larger than Earth's, space-based telescopes have previously detected X-rays only up to ~7 keV. On the basis of theoretical models of the Jovian auroral X-ray production, X-ray emission in the ~2-7 keV band has been interpreted as thermal (arising from electrons characterized by a Maxwell-Boltzmann distribution) bremsstrahlung. Here we report the observation of hard X-rays in the 8-20 keV band from the Jovian aurorae, obtained with the NuSTAR X-ray observatory. The X-rays fit to a flat power-law model with slope 0.60+/-0.22 - a spectral signature of non-thermal, hard X-ray bremsstrahlung. We determine the electron flux and spectral shape in the keV to MeV energy range using coeval in situ measurements by the Juno spacecraft's JADE and JEDI instruments. Jovian electron spectra of the form we observe have previously been interpreted to arise in stochastic acceleration, rather than coherent acceleration by electric fields. We reproduce the X-ray spectral shape and approximate flux observed by NuSTAR, and explain the non-detection of hard X-rays by Ulysses, by simulating the non-thermal population of electrons undergoing precipitating electron energy loss, secondary electron generation and bremsstrahlung emission in a model Jovian atmosphere. The results highlight the similarities between the processes generating hard X-ray auroras on Earth and Jupiter, which may be occurring on Saturn, too. △ Less

Submitted 16 February, 2022; originally announced February 2022.

Comments: 39 pages, 6 figures, 2 tables. Published in Nature Astronomy (https://www.nature.com/articles/s41550-021-01594-8)

Dark Matter Axion Search Using a Josephson Traveling Wave Parametric Amplifier

Authors: C. Bartram, T. Braine, R. Cervantes, N. Crisosto, N. Du, G. Leum, P. Mohapatra, T. Nitta, L. J Rosenberg, G. Rybka, J. Yang, John Clarke, I. Siddiqi, A. Agrawal, A. V. Dixit, M. H. Awida, A. S. Chou, M. Hollister, S. Knirck, A. Sonnenschein, W. Wester, J. R. Gleason, A. T. Hipp, S. Jois, P. Sikivie , et al. (26 additional authors not shown)

Abstract: We present a new exclusion bound of axion-like particle dark matter with axion-photon couplings above $\mathrm{10^{-13}}$ $\mathrm{GeV^{-1}}$ over the frequency range 4796.7--4799.5 MHz, corresponding to a narrow range of axion masses centered around 19.84 $μ$eV. This measurement represents the first implementation of a Josephson Traveling Wave Parametric Amplifier (JTWPA) in a dark matter search.… ▽ More We present a new exclusion bound of axion-like particle dark matter with axion-photon couplings above $\mathrm{10^{-13}}$ $\mathrm{GeV^{-1}}$ over the frequency range 4796.7--4799.5 MHz, corresponding to a narrow range of axion masses centered around 19.84 $μ$eV. This measurement represents the first implementation of a Josephson Traveling Wave Parametric Amplifier (JTWPA) in a dark matter search. The JTWPA was operated in the insert of the Axion Dark Matter eXperiment (ADMX) as part of an independent receiver chain that was attached to a 0.588-liter cavity. The ability of the JTWPA to deliver high gain over a wide (3 GHz) bandwidth has engendered interest from those aiming to perform broadband axion searches, a longstanding goal in this field. △ Less

Submitted 15 October, 2021; originally announced October 2021.

Journal ref: Rev. Sci. Instrum., 94, 044703, 2023

Search for "Invisible" Axion Dark Matter in the $3.3\text{-}4.2~μ$eV Mass Range

Authors: ADMX Collaboration, C. Bartram, T. Braine, E. Burns, R. Cervantes, N. Crisosto, N. Du, H. Korandla, G. Leum, P. Mohapatra, T. Nitta, L. J Rosenberg, G. Rybka, J. Yang, John Clarke, I. Siddiqi, A. Agrawal, A. V. Dixit, M. H. Awida, A. S. Chou, M. Hollister, S. Knirck, A. Sonnenschein, W. Wester, J. R. Gleason , et al. (27 additional authors not shown)

Abstract: We report the results from a haloscope search for axion dark matter in the $3.3\text{-}4.2~μ$eV mass range. This search excludes the axion-photon coupling predicted by one of the benchmark models of "invisible" axion dark matter, the KSVZ model. This sensitivity is achieved using a large-volume cavity, a superconducting magnet, an ultra low noise Josephson parametric amplifier, and sub-Kelvin temp… ▽ More We report the results from a haloscope search for axion dark matter in the $3.3\text{-}4.2~μ$eV mass range. This search excludes the axion-photon coupling predicted by one of the benchmark models of "invisible" axion dark matter, the KSVZ model. This sensitivity is achieved using a large-volume cavity, a superconducting magnet, an ultra low noise Josephson parametric amplifier, and sub-Kelvin temperatures. The validity of our detection procedure is ensured by injecting and detecting blind synthetic axion signals. △ Less

Submitted 29 December, 2021; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: 6 pages, 5 figures

Journal ref: Phys. Rev. Lett. 127, 261803 (2021)

The Milky Way's middle-aged inner ring

Authors: Shola M. Wylie, Jonathan P. Clarke, Ortwin E. Gerhard

Abstract: We investigate the metallicity, age, and orbital anatomy of the inner Milky Way, specifically focussing on the outer bar region. We integrated a sample of APOGEE DR16 inner Galaxy stars in a state of the art bar-bulge potential with a slow pattern speed and investigated the link between the resulting orbits and their [Fe/H] and ages. By superimposing the orbits, we built density, [Fe/H], and age m… ▽ More We investigate the metallicity, age, and orbital anatomy of the inner Milky Way, specifically focussing on the outer bar region. We integrated a sample of APOGEE DR16 inner Galaxy stars in a state of the art bar-bulge potential with a slow pattern speed and investigated the link between the resulting orbits and their [Fe/H] and ages. By superimposing the orbits, we built density, [Fe/H], and age maps of the inner Milky Way, which we divided further using the orbital parameters eccentricity, |Xmax|, and |Zmax|. We find that at low heights from the Galactic plane, the Galactic bar gradually transitions into a radially thick, vertically thin, elongated inner ring with average solar [Fe/H]. This inner ring is mainly composed of stars with AstroNN ages between 4 and 9 Gyr with a peak in age between 6 and 8 Gyr, making the average age of the ring ~6 Gyr. The vertical thickness of the ring decreases markedly towards younger ages. We also find very large L4 Lagrange orbits that have average solar to super-solar metallicities and intermediate ages. Lastly, we confirm a clear X-shape in the [Fe/H] and density distributions at large Galactic heights. The orbital structure obtained for the APOGEE stars reveals that the Milky Way hosts an inner ring-like structure between the planar bar and corotation. This structure is on average metal rich, intermediately aged, and enhances the horizontal metallicity gradient along the bar's major axis. △ Less

Submitted 21 January, 2022; v1 submitted 7 October, 2021; originally announced October 2021.

Comments: 8 pages, 6 figures; Accepted for publication in A&A

Journal ref: A&A 659, A80 (2022)

The Pattern Speed of the Milky Way Bar/Bulge from VIRAC & Gaia

Authors: Jonathan Clarke, Ortwin Gerhard

Abstract: We compare distance resolved, absolute proper motions in the Milky Way bar/bulge region to a grid of made-to-measure dynamical models with well defined pattern speeds. The data are obtained by combining the relative VVV Infrared Astrometric Catalog v1 proper motions with the Gaia DR2 absolute reference frame. We undertake a comprehensive analysis of the various errors in our comparison, from both… ▽ More We compare distance resolved, absolute proper motions in the Milky Way bar/bulge region to a grid of made-to-measure dynamical models with well defined pattern speeds. The data are obtained by combining the relative VVV Infrared Astrometric Catalog v1 proper motions with the Gaia DR2 absolute reference frame. We undertake a comprehensive analysis of the various errors in our comparison, from both the data and the models, and allow for additional, unknown, contributions by using an outlier-tolerant likelihood function to evaluate the best fitting model. We quantify systematic effects such as the region of data included in the comparison, the possible overlap from spiral arms, and the choice of synthetic luminosity function and bar angle used to predict the data from the models. Resulting variations in the best-fit parameters are included in their final errors. We thus measure the bar pattern speed to be $Ω_b=33.29\pm1.81 \mathrm{km \, s^{-1} \, kpc^{-1}}$ and the azimuthal solar velocity to be $V_{φ,\,\odot} = 251.31 \pm 1.95 \mathrm{km \, s^{-1}}$. These values, when combined with recent measurements of the Galactic rotation curve, yield the distance of corotation, 6.5 < $R_{CR}$ [kpc] < 7.5, the outer Lindblad resonance (OLR), 10.7 < $R_{OLR}$ [kpc] < 12.4, and the higher order, m = 4, OLR, 8.7 < $R_{{OLR}_4}$ [kpc] < 10.0. The measured pattern speed provides strong evidence for the "long-slow" bar scenario. △ Less

Submitted 28 April, 2022; v1 submitted 22 July, 2021; originally announced July 2021.

Comments: Published in MNRAS, 19 pages, 17 figures Note: Model grid slightly altered from original preprint. Measurement of the pattern speed adjusted down slightly

A2A: 21,000 bulge stars from the ARGOS survey with stellar parameters on the APOGEE scale

Authors: S. M. Wylie, O. E. Gerhard, M. K. Ness, J. P. Clarke, K. C. Freeman, J. Bland-Hawthorn

Abstract: We use the data-driven method, The Cannon, to bring 21,000 stars from the ARGOS bulge survey, including 10,000 red clump stars, onto the parameter and abundance scales of the cross-Galactic survey, APOGEE, obtaining rms precisions of 0.10 dex, 0.07 dex, 74 K, and 0.18 dex for [Fe/H], [Mg/Fe], Teff, and log(g), respectively. The re-calibrated ARGOS survey - which we refer to as the A2A survey - is… ▽ More We use the data-driven method, The Cannon, to bring 21,000 stars from the ARGOS bulge survey, including 10,000 red clump stars, onto the parameter and abundance scales of the cross-Galactic survey, APOGEE, obtaining rms precisions of 0.10 dex, 0.07 dex, 74 K, and 0.18 dex for [Fe/H], [Mg/Fe], Teff, and log(g), respectively. The re-calibrated ARGOS survey - which we refer to as the A2A survey - is combined with the APOGEE survey to investigate the abundance structure of the Galactic bulge. We find X-shaped [Fe/H] and [Mg/Fe] distributions in the bulge that are more pinched than the bulge density, a signature of its disk origin. The mean abundance along the major axis of the bar varies such that the stars are more [Fe/H]-poor and [Mg/Fe]-rich near the Galactic center than in the long bar/outer bulge region. The vertical [Fe/H] and [Mg/Fe] gradients vary between the inner bulge and long bar with the inner bulge showing a flattening near the plane that is absent in the long bar. The [Fe/H]-[Mg/Fe] distribution shows two main maxima, an ``[Fe/H]-poor [Mg/Fe]- rich'' maximum and an ``[Fe/H]-rich [Mg/Fe]-poor'' maximum, that vary in strength with position in the bulge. In particular, the outer long bar close to the Galactic plane is dominated by super-solar [Fe/H], [Mg/Fe]-normal stars. Stars composing the [Fe/H]-rich maximum show little kinematic dependence on [Fe/H], but for lower [Fe/H] the rotation and dispersion of the bulge increase slowly. Stars with [Fe/H]<-1 dex have a very different kinematic structure than stars with higher [Fe/H]. Comparing with recent models for the Galactic boxy-peanut bulge, the abundance gradients and distribution, and the relation between [Fe/H] and kinematics suggest that the stars comprising each maximum have separate disk origins with the ``[Fe/H]-poor [Mg/Fe]-rich'' stars originating from a thicker disk than the ``[Fe/H]-rich [Mg/Fe]-poor'' stars. △ Less

Submitted 27 June, 2021; originally announced June 2021.

Comments: 26 pages, 29 figures, revised version as submitted to A&A following referee report

Journal ref: A&A 653, A143 (2021)

The general applicability of self-similar solutions for thermal disc winds

Authors: Andrew D. Sellek, Cathie J. Clarke, Richard A. Booth

Abstract: Thermal disc winds occur in many contexts and may be particularly important to the secular evolution and dispersal of protoplanetary discs heated by high energy radiation from their central star. In this paper we generalise previous models of self-similar thermal winds - which have self-consistent morphology and variation of flow variables - to the case of launch from an elevated base and to non-i… ▽ More Thermal disc winds occur in many contexts and may be particularly important to the secular evolution and dispersal of protoplanetary discs heated by high energy radiation from their central star. In this paper we generalise previous models of self-similar thermal winds - which have self-consistent morphology and variation of flow variables - to the case of launch from an elevated base and to non-isothermal conditions. These solutions are well-reproduced by hydrodynamic simulations, in which, as in the case of isothermal winds launched from the mid-plane, we find winds launch at the maximum Mach number for which the streamline solutions extend to infinity without encountering a singularity. We explain this behaviour based on the fact that lower Mach number solutions do not fill the spatial domain. We also show that hydrodynamic simulations reflect the corresponding self-similar models across a range of conditions appropriate to photoevaporating protoplanetary discs, even when gravity, centrifugal forces, or changes in the density gradient mean the problem is not inherently scale free. Of all the parameters varied, the elevation of the wind base affected the launch velocity and flow morphology most strongly, with temperature gradients causing only minor differences. We explore how launching from an elevated base affects Ne II line profiles from winds, finding it increases (reduces) the full width at half maximum (FWHM) of the line at low (high) inclination to the line of sight compared with models launched from the disc mid-plane and thus weakens the dependence of the FWHM on inclination. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Comments: 21 pages, 16 figures, accepted after revision for publication in MNRAS

ALMA imaging of the M-dwarf Fomalhaut C's debris disc

Authors: Patrick F. Cronin-Coltsmann, Grant M. Kennedy, Paul Kalas, Julien Milli, Cathie J. Clarke, Gaspard Duchêne, Jane Greaves, Samantha M. Lawler, Jean-François Lestrade, Brenda C. Matthews, Andrew Shannon, Mark C. Wyatt

Abstract: Fomalhaut C (LP 876-10) is a low mass M4V star in the intriguing Fomalhaut triple system and, like Fomalhaut A, possesses a debris disc. It is one of very few nearby M-dwarfs known to host a debris disc and of these has by far the lowest stellar mass. We present new resolved observations of the debris disc around Fomalhaut C with the Atacama Large Millimetre Array which allow us to model its prope… ▽ More Fomalhaut C (LP 876-10) is a low mass M4V star in the intriguing Fomalhaut triple system and, like Fomalhaut A, possesses a debris disc. It is one of very few nearby M-dwarfs known to host a debris disc and of these has by far the lowest stellar mass. We present new resolved observations of the debris disc around Fomalhaut C with the Atacama Large Millimetre Array which allow us to model its properties and investigate the system's unique history. The ring has a radius of 26 au and a narrow full width at half maximum of at most 4.2 au. We find a 3$σ$ upper limit on the eccentricity of 0.14, neither confirming nor ruling out previous dynamic interactions with Fomalhaut A that could have affected Fomalhaut C's disc. We detect no $^{12}$CO J=3-2 emission in the system and do not detect the disc in scattered light with HST/STIS or VLT/SPHERE. We find the original Herschel detection to be consistent with our ALMA model's radial size. We place the disc in the context of the wider debris disc population and find that its radius is as expected from previous disc radius-host luminosity trends. Higher signal-to-noise observations of the system would be required to further constrain the disc properties and provide further insight to the history of the Fomalhaut triple system as a whole. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: 15 pages, 10 figures, 2 tables, Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 504, Issue 3, July 2021, Pages 4497-4510

The effect of the streaming instability on protoplanetary disc emission at millimetre wavelengths

Authors: Chiara E. Scardoni, Richard A. Booth, Cathie J. Clarke

Abstract: In this paper, we investigate whether overdensity formation via streaming instability is consistent with recent multi-wavelength ALMA observations in the Lupus star forming region. We simulate the local action of streaming instability in 2D using the code ATHENA, and examine the radiative properties at mm wavelengths of the resulting clumpy dust distribution by focusing on two observable quantitie… ▽ More In this paper, we investigate whether overdensity formation via streaming instability is consistent with recent multi-wavelength ALMA observations in the Lupus star forming region. We simulate the local action of streaming instability in 2D using the code ATHENA, and examine the radiative properties at mm wavelengths of the resulting clumpy dust distribution by focusing on two observable quantities: the optically thick fraction $ff$ (in ALMA band 6) and the spectral index $α$ (in bands 3-7). By comparing the simulated distribution in the $ff-α$ plane before and after the action of streaming instability, we observe that clump formation causes $ff$ to drop, because of the suppression of emission from grains that end up in optically thick clumps. $α$, instead, can either increase or decline after the action of streaming instability; we use a simple toy model to demonstrate that this behaviour depends on the sizes of the grains whose emission is suppressed by being incorporated in optically thick clumps. In particular, the sign of evolution of $α$ depends on whether grains near the opacity maximum at a few tenths of a mm end up in clumps. By comparing the simulation distributions before/after clump formation to the data distribution, we note that the action of streaming instability drives simulations towards the area of the plane where the data are located. We furthermore demonstrate that this behaviour is replicated in integrated disc models provided that the instability is operative over a region of the disc that contributes significantly to the total mm flux. △ Less

Submitted 24 March, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

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

Gas Dynamics in the Galaxy: Total Mass Distribution and the Bar Pattern Speed

Authors: Zhi Li, Juntai Shen, Ortwin Gerhard, Jonathan P. Clarke

Abstract: Gas morphology and kinematics in the Milky Way contain key information for understanding the formation and evolution of our Galaxy. We present a high resolution hydrodynamical simulation based on a realistic barred Milky Way potential constrained by recent observations. Our model can reproduce most features in the observed longitude-velocity diagram, including the Central Molecular Zone, the Near… ▽ More Gas morphology and kinematics in the Milky Way contain key information for understanding the formation and evolution of our Galaxy. We present a high resolution hydrodynamical simulation based on a realistic barred Milky Way potential constrained by recent observations. Our model can reproduce most features in the observed longitude-velocity diagram, including the Central Molecular Zone, the Near and Far 3-kpc arms, the Molecular Ring, and the spiral arm tangents. It can also explain the non-circular motions of masers obtained by the recent BeSSeL2 survey. The central gas kinematics are consistent with a mass of $6.9\times10^8\; {\rm M}_{\odot}$ in the Nuclear Stellar Disk. Our model predicts the formation of an elliptical gaseous ring surrounding the bar, which is composed of the 3-kpc arms, Norma arm, and the bar-spiral interfaces. This ring is similar to those "inner" rings in some Milky Way analogs with a boxy/peanut-shaped bulge. The kinematics of gas near the solar neighbourhood are governed by the Local arm, which is induced by the four major stellar spiral arms. The bar pattern speed constrained by our gas model is $37.5-40\; {\rm km}\;{\rm s}^{-1}\;{\rm kpc}^{-1}$, corresponding to a corotation radius of $R_{\rm CR}=6.0-6.4\;{\rm kpc}$. The rotation curve of our model rises gently within the central $\sim5\;{\rm kpc}$, which is significantly less steep than those predicted by modern zoom-in cosmological simulations such as Auriga. △ Less

Submitted 10 November, 2021; v1 submitted 18 March, 2021; originally announced March 2021.

Comments: 23 pages, 14 figures. ApJ accepted, comments are welcome

A super-resolution analysis of the DSHARP survey: Substructure is common in the inner 30 au

Authors: Jeff Jennings, Richard A. Booth, Marco Tazzari, Cathie J. Clarke, Giovanni P. Rosotti

Abstract: The DSHARP survey evidenced the ubiquity of substructure in the mm dust distribution of large, bright protoplanetary discs. Intriguingly, these datasets have yet higher resolution information that is not recovered in a CLEAN image. We first show that the intrinsic performance of the CLEAN algorithm is resolution-limited. Then analyzing all 20 DSHARP sources using the 1D, super-resolution code Fran… ▽ More The DSHARP survey evidenced the ubiquity of substructure in the mm dust distribution of large, bright protoplanetary discs. Intriguingly, these datasets have yet higher resolution information that is not recovered in a CLEAN image. We first show that the intrinsic performance of the CLEAN algorithm is resolution-limited. Then analyzing all 20 DSHARP sources using the 1D, super-resolution code Frankenstein (frank), we accurately fit the 1D visibilities to a mean factor of 4.3 longer baseline than the Fourier transform of the CLEAN images and a factor of 3.0 longer baseline than the transform of the CLEAN component models. This yields a higher resolution brightness profile for each source, identifying new substructure interior to 30 au in multiple discs; resolving known gaps to be deeper, wider, and more structured; and known rings to be narrower and brighter. Across the survey, high contrast gaps are an average 14% wider and 44% deeper in the frank profiles relative to CLEAN, and high contrast rings are an average 26% narrower. Categorizing the frank brightness profiles into trends, we find that the relative scarcity of features interior to 30 au in the survey's CLEAN images is an artifact of resolving power, rather than an intrinsic rarity of inner disc (or compact disc) substructure. Finally the rings in the frank profiles are narrower than the previously inferred deconvolved widths, indicating smaller alpha / St ratios in the local gas disc. △ Less

Submitted 1 November, 2021; v1 submitted 3 March, 2021; originally announced March 2021.

Comments: Revised version, accepted to MNRAS. Changes: added CLEAN .model profiles to Sec.3 and 4, quantified gaps/rings in frank profiles in Sec.5 and Table 3, expanded analysis of alpha/St in Sec.5.3.1. All frank fits are now available at https://zenodo.org/record/5587841

Modeling the delivery of dust from discs to ionized winds

Authors: Richard A. Booth, Cathie J. Clarke

Abstract: A necessary first step for dust removal in protoplanetary disc winds is the delivery of dust from the disc to the wind. In the case of ionized winds, the disc and wind are sharply delineated by a narrow ionization front where the gas density and temperature vary by more than an order of magnitude. Using a novel method that is able to model the transport of dust across the ionization front in the p… ▽ More A necessary first step for dust removal in protoplanetary disc winds is the delivery of dust from the disc to the wind. In the case of ionized winds, the disc and wind are sharply delineated by a narrow ionization front where the gas density and temperature vary by more than an order of magnitude. Using a novel method that is able to model the transport of dust across the ionization front in the presence of disc turbulence, we revisit the problem of dust delivery. Our results show that the delivery of dust to the wind is determined by the vertical gas flow through the disc induced by the mass loss, rather than turbulent diffusion (unless the turbulence is strong, i.e. $α\gtrsim 0.01$). Using these results we provide a simple relation between the maximum size of particle that can be delivered to the wind and the local mass-loss rate per unit area from the wind. This relation is independent of the physical origin of the wind and predicts typical sizes in the 0.01 -- $1\,μm$ range for EUV or X-ray driven winds. These values are a factor $\sim 10$ smaller than those obtained when considering only whether the wind is able to carry away the grains. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: Accepted for publication in MNRAS

Proplyds in the Flame Nebula NGC 2024

Authors: Thomas J. Haworth, Jinyoung S. Kim, Andrew J. Winter, Dean C. Hines, Cathie J. Clarke, Andrew D. Sellek, Giulia Ballabio, Karl R. Stapelfeldt

Abstract: A recent survey of the inner $0.35\times0.35$pc of the NGC 2024 star forming region revealed two distinct millimetre continuum disc populations that appear to be spatially segregated by the boundary of a dense cloud. The eastern (and more embedded) population is $\sim0.2-0.5$Myr old, with an ALMA mm continuum disc detection rate of about $45\,$per cent. However this drops to only $\sim15$per cent… ▽ More A recent survey of the inner $0.35\times0.35$pc of the NGC 2024 star forming region revealed two distinct millimetre continuum disc populations that appear to be spatially segregated by the boundary of a dense cloud. The eastern (and more embedded) population is $\sim0.2-0.5$Myr old, with an ALMA mm continuum disc detection rate of about $45\,$per cent. However this drops to only $\sim15$per cent in the 1Myr western population. When presenting this result, van Terwisga et al. (2020) suggested that the two main UV sources, IRS 1 (a B0.5V star in the western region) and IRS 2b (an O8V star in the eastern region, but embedded) have both been evaporating the discs in the depleted western population. In this paper we report the firm discovery in archival HST data of 4 proplyds and 4 further candidate proplyds in NGC 2024, confirming that external photoevaporation of discs is occurring. However, the locations of these proplyds changes the picture. Only three of them are in the depleted western population and their evaporation is dominated by IRS 1, with no obvious impact from IRS 2b. The other 5 proplyds are in the younger eastern region and being evaporated by IRS 2b. We propose that both populations are subject to significant external photoevaporation, which happens throughout the region wherever discs are not sufficiently shielded by the interstellar medium. The external photoevaporation and severe depletion of mm grains in the 0.2-0.5Myr eastern part of NGC 2024 would be in competition even with very early planet formation. △ Less

Submitted 16 December, 2020; originally announced December 2020.

Comments: 14 pages, Accepted for publication in MNRAS

Planet formation in intermediate-separation binary systems

Authors: O. Panić, T. J. Haworth, M. G. Petr-Gotzens, J. Miley, M. van den Ancker, M. Vioque, L. Siess, R. Parker, C. J. Clarke, I. Kamp, G. Kennedy, R. D. Oudmaijer, I. Pascucci, A. M. S. Richards, T. Ratzka, C. Qi

Abstract: We report the first characterisation of the individual discs in the intermediate separation binary systems KK Oph and HD 144668 at millimetre wavelengths. In both systems the circum-primary and the circum-secondary discs are detected in the millimetre continuum emission, but not in $^{13}$CO nor C$^{18}$O lines. Even though the disc structure is only marginally resolved, we find indications of lar… ▽ More We report the first characterisation of the individual discs in the intermediate separation binary systems KK Oph and HD 144668 at millimetre wavelengths. In both systems the circum-primary and the circum-secondary discs are detected in the millimetre continuum emission, but not in $^{13}$CO nor C$^{18}$O lines. Even though the disc structure is only marginally resolved, we find indications of large-scale asymmetries in the outer regions of the primary discs, most likely due to perturbation by the companion. The derived dust masses are firmly above debris disc level for all stars. The primaries have about three times more dust in their discs than the secondaries. In the case of HD 144668 the opacity spectral index of the primary and secondary differ by the large margin of 0.69 which may be a consequence of the secondary disc being more compact. Upper limits on the gas masses imply less than 0.1 M$_{\textrm{jup}}$ in any of these discs, meaning that giant planets can no longer form in them. Considering that there have been no massive gas discs identified to date in intermediate separation binaries (i.e., binaries at a few hundred au separation), this opens space for speculation whether their binarity causes the removal of gas, with tidal interaction truncating the discs and hence shortening the accretion timescale. More systematic studies in this respect are sorely needed. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: 12 pages. Accepted for publication in MNRAS

Dust delivery and entrainment in photoevaporative winds

Authors: Mark A. Hutchison, Cathie J. Clarke

Abstract: We model the gas and dust dynamics in a turbulent protoplanetary disc undergoing extreme-UV photoevaporation in order to better characterise the dust properties in thermal winds (e.g. size distribution, flux rate, trajectories). Our semi-analytic approach allows us to rapidly calculate these dust properties without resorting to expensive hydrodynamic simulations. We find that photoevaporation crea… ▽ More We model the gas and dust dynamics in a turbulent protoplanetary disc undergoing extreme-UV photoevaporation in order to better characterise the dust properties in thermal winds (e.g. size distribution, flux rate, trajectories). Our semi-analytic approach allows us to rapidly calculate these dust properties without resorting to expensive hydrodynamic simulations. We find that photoevaporation creates a vertical gas flow within the disc that assists turbulence in supplying dust to the ionisation front. We examine both the delivery of dust to the ionisation front and its subsequent entrainment in the overlying wind. We derive a simple analytic criterion for the maximum grain size that can be entrained and show that this is in good agreement with the results of previous simulations where photoevaporation is driven by a range of radiation types. We show that, in contrast to the case for magnetically driven winds, we do not expect large scale dust transport within the disc to be effected by photoevaporation. We also show that the maximum size of grains that can be entrained in the wind ($s_{\rm max}$) is around an order of magnitude larger than the maximum size of grains that can be delivered to the front by advection alone ($s_{\rm crit} \lesssim 1$ $μ{\rm m}$ for Herbig Ae/Be stars and $\lesssim 0.01$ $μ{\rm m}$ for T Tauri stars). We further investigate how larger grains, up to a limiting size $s_{\rm limit}$, can be delivered to the front by turbulent diffusion alone. In all cases, we find $s_{\rm max} \gtrsim s_{\rm limit}$ so that we expect that any dust that is delivered to the front can be entrained in the wind and that most entrained dust follows trajectories close to that of the gas. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: 16 pages, 6 figures, Accepted for publication in MNRAS

Axion Dark Matter eXperiment: Run 1B Analysis Details

Authors: ADMX Collaboration, C. Bartram, T. Braine, R. Cervantes, N. Crisosto, N. Du, G. Leum, L. J Rosenberg, G. Rybka, J. Yang, D. Bowring, A. S. Chou, R. Khatiwada, A. Sonnenschein, W. Wester, G. Carosi, N. Woollett, L. D. Duffy, M. Goryachev, B. McAllister, M. E. Tobar, C. Boutan, M. Jones, B. H. Laroque, N. S. Oblath , et al. (23 additional authors not shown)

Abstract: Searching for axion dark matter, the ADMX collaboration acquired data from January to October 2018, over the mass range 2.81--3.31 $μ$eV, corresponding to the frequency range 680--790 MHz. Using an axion haloscope consisting of a microwave cavity in a strong magnetic field, the ADMX experiment excluded Dine-Fischler-Srednicki-Zhitnisky (DFSZ) axions at 100% dark matter density over this entire fre… ▽ More Searching for axion dark matter, the ADMX collaboration acquired data from January to October 2018, over the mass range 2.81--3.31 $μ$eV, corresponding to the frequency range 680--790 MHz. Using an axion haloscope consisting of a microwave cavity in a strong magnetic field, the ADMX experiment excluded Dine-Fischler-Srednicki-Zhitnisky (DFSZ) axions at 100% dark matter density over this entire frequency range, except for a few gaps due to mode crossings. This paper explains the full ADMX analysis for Run 1B, motivating analysis choices informed by details specific to this run. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Comments: 17 pages, 17 figures

Journal ref: Phys. Rev. D 103, 032002 (2021)

Multi-wavelength continuum sizes of protoplanetary discs: scaling relations and implications for grain growth and radial drift

Authors: Marco Tazzari, Cathie J. Clarke, Leonardo Testi, Jonathan P. Williams, Stefano Facchini, Carlo F. Manara, Antonella Natta, Giovanni Rosotti

Abstract: We analyse spatially resolved ALMA observations at 0.9, 1.3, and 3.1 mm for the 26 brightest protoplanetary discs in the Lupus star-forming region. We characterise the discs multi-wavelength brightness profiles by fitting the interferometric visibilities in a homogeneous way, obtaining effective disc sizes at the three wavelengths, spectral index profiles and optical depth estimates. We report thr… ▽ More We analyse spatially resolved ALMA observations at 0.9, 1.3, and 3.1 mm for the 26 brightest protoplanetary discs in the Lupus star-forming region. We characterise the discs multi-wavelength brightness profiles by fitting the interferometric visibilities in a homogeneous way, obtaining effective disc sizes at the three wavelengths, spectral index profiles and optical depth estimates. We report three fundamental discoveries: first, the millimeter continuum size - luminosity relation already observed at 0.9 mm is also present at 1.3 mm with an identical slope, and at 3.1 mm with a steeper slope, confirming that emission at longer wavelengths becomes increasingly optically thin. Second, when observed at 3.1 mm the discs appear to be only 9% smaller than when observed at 0.9 mm, in tension with models of dust evolution which predict a starker difference. Third, by forward modelling the sample of measurements with a simple parametric disc model, we find that the presence of large grains ($a_\mathrm{max}>1 $mm) throughout the discs is the most favoured explanation for all discs as it reproduces simultaneously their spectral indices, optical depth, luminosity, and radial extent in the 0.9-1.3 mm wavelength range. We also find that the observations can be alternatively interpreted with the discs being dominated by optically thick, unresolved, substructures made of mm-sized grains with a high scattering albedo. △ Less

Submitted 23 June, 2021; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: MNRAS Accepted. 20 pages, 12 figures, 5 tables. Detailed fit results in Appendix E (27 pages, 78 figures), available as supplementary material. Machine-readable version of Table 1, Table 2, and Table 3 available at https://zenodo.org/record/4756381

Axion Dark Matter eXperiment: Detailed Design and Operations

Authors: R. Khatiwada, D. Bowring, A. S. Chou, A. Sonnenschein, W. Wester, D. V. Mitchell, T. Braine, C. Bartram, R. Cervantes, N. Crisosto, N. Du, S. Kimes, L. J Rosenberg, G. Rybka, J. Yang, D. Will, G. Carosi, N. Woollett, S. Durham, L. D. Duffy, R. Bradley, C. Boutan, M. Jones, B. H. LaRoque, N. S. Oblath , et al. (26 additional authors not shown)

Abstract: Axion Dark Matter eXperiment (ADMX) ultra low noise haloscope technology has enabled the successful completion of two science runs (1A and 1B) that looked for dark matter axions in the $2.66$ to $3.1$ $μ$eV mass range with Dine-Fischler-Srednicki-Zhitnisky (DFSZ) sensitivity Ref. [1,2]. Therefore, it is the most sensitive axion search experiment to date in this mass range. We discuss the technolog… ▽ More Axion Dark Matter eXperiment (ADMX) ultra low noise haloscope technology has enabled the successful completion of two science runs (1A and 1B) that looked for dark matter axions in the $2.66$ to $3.1$ $μ$eV mass range with Dine-Fischler-Srednicki-Zhitnisky (DFSZ) sensitivity Ref. [1,2]. Therefore, it is the most sensitive axion search experiment to date in this mass range. We discuss the technological advances made in the last several years to achieve this sensitivity, which includes the implementation of components, such as state-of-the-art quantum limited amplifiers and a dilution refrigerator. Furthermore, we demonstrate the use of a frequency tunable Microstrip Superconducting Quantum Interference Device (SQUID) Amplifier (MSA), in Run 1A, and a Josephson Parametric Amplifier (JPA), in Run 1B, along with novel analysis tools that characterize the system noise temperature. △ Less

Submitted 30 September, 2020; originally announced October 2020.

Comments: 23 pages, 28 figures

Report number: FERMILAB-PUB-20-331-AD-E-QIS

The science enabled by a dedicated solar system space telescope

Authors: Cindy L. Young, Michael H. Wong, Kunio M. Sayanagi, Shannon Curry, Kandis L. Jessup, Tracy Becker, Amanda Hendrix, Nancy Chanover, Stephanie Milam, Bryan J. Holler, Gregory Holsclaw, Javier Peralta, John Clarke, John Spencer, Michael S. P. Kelley, Janet Luhmann, David MacDonnell, Ronald J. Vervack Jr., Kurt Retherford, Leigh N. Fletcher, Imke de Pater, Faith Vilas, Lori Feaga, Oswald Siegmund, Jim Bell , et al. (13 additional authors not shown)

Abstract: The National Academy Committee on Astrobiology and Planetary Science (CAPS) made a recommendation to study a large/medium-class dedicated space telescope for planetary science, going beyond the Discovery-class dedicated planetary space telescope endorsed in Visions and Voyages. Such a telescope would observe targets across the entire solar system, engaging a broad spectrum of the science community… ▽ More The National Academy Committee on Astrobiology and Planetary Science (CAPS) made a recommendation to study a large/medium-class dedicated space telescope for planetary science, going beyond the Discovery-class dedicated planetary space telescope endorsed in Visions and Voyages. Such a telescope would observe targets across the entire solar system, engaging a broad spectrum of the science community. It would ensure that the high-resolution, high-sensitivity observations of the solar system in visible and UV wavelengths revolutionized by the Hubble Space Telescope (HST) could be extended. A dedicated telescope for solar system science would: (a) transform our understanding of time-dependent phenomena in our solar system that cannot be studied currently under programs to observe and visit new targets and (b) enable a comprehensive survey and spectral characterization of minor bodies across the solar system, which requires a large time allocation not supported by existing facilities. The time-domain phenomena to be explored are critically reliant on high spatial resolution UV-visible observations. This paper presents science themes and key questions that require a long-lasting space telescope dedicated to planetary science that can capture high-quality, consistent data at the required cadences that are free from effects of the terrestrial atmosphere and differences across observing facilities. Such a telescope would have excellent synergy with astrophysical facilities by placing planetary discoveries made by astrophysics assets in temporal context, as well as triggering detailed follow-up observations using larger telescopes. The telescope would support future missions to the Ice Giants, Ocean Worlds, and minor bodies across the solar system by placing the results of such targeted missions in the context of longer records of temporal activities and larger sample populations. △ Less

Submitted 18 August, 2020; originally announced August 2020.

Comments: A whitepaper submitted to the Planetary Science Decadal Survey

A dusty origin for the correlation between protoplanetary disc accretion rates and dust masses

Authors: Andrew D. Sellek, Richard A. Booth, Cathie J. Clarke

Abstract: Recent observations have uncovered a correlation between the accretion rates (measured from the UV continuum excess) of protoplanetary discs and their masses inferred from observations of the sub-mm continuum. While viscous evolution models predict such a correlation, the predicted values are in tension with data obtained from the Lupus and Upper Scorpius star forming regions; for example, they un… ▽ More Recent observations have uncovered a correlation between the accretion rates (measured from the UV continuum excess) of protoplanetary discs and their masses inferred from observations of the sub-mm continuum. While viscous evolution models predict such a correlation, the predicted values are in tension with data obtained from the Lupus and Upper Scorpius star forming regions; for example, they underpredict the scatter in accretion rates, particularly in older regions. Here we argue that since the sub-mm observations trace the discs' dust, by explicitly modelling the dust grain growth, evolution, and emission, we can better understand the correlation. We show that for turbulent viscosities with $α\lesssim 10^{-3}$, the depletion of dust from the disc due to radial drift means we can reproduce the range of masses and accretion rates seen in the Lupus and Upper Sco datasets. One consequence of this model is that the upper locus of accretion rates at a given dust mass does not evolve with the age of the region. Moreover, we find that internal photoevaporation is necessary to produce the lowest accretion rates observed. In order to replicate the correct dust masses at the time of disc dispersal, we favour relatively low photoevaporation rates $\lesssim 10^{-9}~M_{\odot}~\mathrm{yr^{-1}}$ for most sources but cannot discriminate between EUV or X-ray driven winds. A limited number of sources, particularly in Lupus, are shown to have higher masses than predicted by our models which may be evidence for variations in the properties of the dust or dust trapping induced in substructures. △ Less

Submitted 17 August, 2020; originally announced August 2020.

Comments: Accepted to MNRAS

Architectures and Technologies for a Space Telescope for Solar System Science

Authors: Kunio M. Sayanagi, Cindy L. Young, Lynn Bowman, Joseph Pitman, Bo Naasz, Bonnie Meinke, Tracy Becker, Jim Bell, Richard Cartwright, Nancy Chanover, John Clarke, Joshua Colwell, Shannon Curry, Imke de Pater, Gregory Delory, Lori Feaga, Leigh N. Fletcher, Thomas Greathouse, Amanda Hendrix, Bryan J. Holler, Gregory Holsclaw, Kandis L. Jessup, Michael S. P. Kelley, Robert Lillis, Rosaly M. C. Lopes , et al. (15 additional authors not shown)

Abstract: We advocate for a mission concept study for a space telescope dedicated to solar system science in Earth orbit. Such a study was recommended by the Committee on Astrobiology and Planetary Science (CAPS) report "Getting Ready for the Next Planetary Science Decadal Survey." The Mid-Decadal Review also recommended NASA to assess the role and value of space telescopes for planetary science. The need f… ▽ More We advocate for a mission concept study for a space telescope dedicated to solar system science in Earth orbit. Such a study was recommended by the Committee on Astrobiology and Planetary Science (CAPS) report "Getting Ready for the Next Planetary Science Decadal Survey." The Mid-Decadal Review also recommended NASA to assess the role and value of space telescopes for planetary science. The need for high-resolution, UV-Visible capabilities is especially acute for planetary science with the impending end of the Hubble Space Telescope (HST); however, NASA has not funded a planetary telescope concept study, and the need to assess its value remains. Here, we present potential design options that should be explored to inform the decadal survey. △ Less

Submitted 15 August, 2020; originally announced August 2020.

Comments: Whitepaper submitted to Planetary Science and Astrobiology Decadal Survey

Ultraviolet-Based Science in the Solar System: Advances and Next Steps

Authors: Amanda R. Hendrix, Tracy M. Becker, Dennis Bodewits, E. Todd Bradley, Shawn Brooks, Ben Byron, Josh Cahill, John Clarke, Lori Feaga, Paul Feldman, G. Randall Gladstone, Candice J. Hansen, Charles Hibbitts, Tommi T. Koskinen, Lizeth Magana, Philippa Molyneux, Shouleh Nikzad, John Noonan, Wayne Pryor, Ujjwal Raut, Kurt D. Retherford, Lorenz Roth, Emilie Royer, Ella Sciamma-O'Brien, Alan Stern , et al. (3 additional authors not shown)

Abstract: We review the importance of recent UV observations of solar system targets and discuss the need for further measurements, instrumentation and laboratory work in the coming decade. In the past decade, numerous important advances have been made in solar system science using ultraviolet (UV) spectroscopic techniques. Formerly used nearly exclusively for studies of giant planet atmospheres, planetar… ▽ More We review the importance of recent UV observations of solar system targets and discuss the need for further measurements, instrumentation and laboratory work in the coming decade. In the past decade, numerous important advances have been made in solar system science using ultraviolet (UV) spectroscopic techniques. Formerly used nearly exclusively for studies of giant planet atmospheres, planetary exospheres and cometary emissions, UV imaging spectroscopy has recently been more widely applied. The geyser-like plume at Saturn's moon Enceladus was discovered in part as a result of UV stellar occultation observations, and this technique was used to characterize the plume and jets during the entire Cassini mission. Evidence for a similar style of activity has been found at Jupiter's moon Europa using Hubble Space Telescope (HST) UV emission and absorption imaging. At other moons and small bodies throughout the solar system, UV spectroscopy has been utilized to search for activity, probe surface composition, and delineate space weathering effects; UV photometric studies have been used to uncover regolith structure. Insights from UV imaging spectroscopy of solar system surfaces have been gained largely in the last 1-2 decades, including studies of surface composition, space weathering effects (e.g. radiolytic products) and volatiles on asteroids (e.g. [2][39][48][76][84]), the Moon (e.g. [30][46][49]), comet nuclei (e.g. [85]) and icy satellites (e.g. [38][41-44][45][47][65]). The UV is sensitive to some species, minor contaminants and grain sizes often not detected in other spectral regimes. In the coming decade, HST observations will likely come to an end. New infrastructure to bolster future UV studies is critically needed. These needs include both developmental work to help improve future UV observations and laboratory work to help interpret spacecraft data. UV instrumentation will be a critical tool on missions to a variety of targets in the coming decade, especially for the rapidly expanding application of UV reflectance investigations of atmosphereless bodies. △ Less

Submitted 28 July, 2020; originally announced July 2020.

Forbidden line diagnostics of photoevaporative disc winds

Authors: G. Ballabio, R. D. Alexander, C. J. Clarke

Abstract: Photoevaporation driven by high energy radiation from the central star plays an important role in the evolution of protoplanetary discs. Photoevaporative winds have been unambiguously detected through blue-shifted emission lines, but their detailed properties remain uncertain. Here we present a new empirical approach to make observational predictions of these thermal winds, seeking to fill the gap… ▽ More Photoevaporation driven by high energy radiation from the central star plays an important role in the evolution of protoplanetary discs. Photoevaporative winds have been unambiguously detected through blue-shifted emission lines, but their detailed properties remain uncertain. Here we present a new empirical approach to make observational predictions of these thermal winds, seeking to fill the gap between theory and observations. We use a self-similar model of an isothermal wind to compute line profiles of several characteristic emission lines (the [Ne${\rm{\scriptsize II}}$] line at 12.81 $μ$m, and optical forbidden lines such as [O${\rm{\scriptsize I}}$] 6300 $\mathring{A}$ and [S${\rm{\scriptsize II}}$] 4068/4076 $\mathring{A}$), studying how the lines are affected by parameters such as the gas temperature, disc inclinations, and density profile. Our model successfully reproduces blue-shifted lines with $v_{\rm peak} \lesssim 10$ km/s, which decrease with increasing disc inclination. The line widths increase with increasing disc inclinations and range from $Δv \sim 15-30$ km/s. The predicted blue-shifts are mostly sensitive to the gas sound speed. The observed [Ne${\rm{\scriptsize II}}$] line profiles are consistent with a thermal wind and point towards a relatively high sound speed, as expected for EUV photoevaporation. However, the observed [O${\rm{\scriptsize I}}$] line profiles require lower temperatures, as expected in X-ray photoevaporation, and show a wider scatter that is difficult to reconcile with a single wind model; it seems likely that these lines trace different components of a multi-phase wind. We also note that the spectral resolution of current observations remains an important limiting factor in these studies, and that higher resolution spectra are required if emission lines are to further our understanding of protoplanetary disc winds. △ Less

Submitted 17 June, 2020; originally announced June 2020.

Comments: 14 pages, 10 figures. Accepted for publication in MNRAS

Frankenstein: Protoplanetary disc brightness profile reconstruction at sub-beam resolution with a rapid Gaussian process

Authors: Jeff Jennings, Richard A. Booth, Marco Tazzari, Giovanni P. Rosotti, Cathie J. Clarke

Abstract: Interferometric observations of the mm dust distribution in protoplanetary discs are now showing a ubiquity of annular gap and ring substructures. Their identification and accurate characterization is critical to probing the physical processes responsible. We present Frankenstein (frank), an open source code that recovers axisymmetric disc structures at sub-beam resolution. By fitting the visibili… ▽ More Interferometric observations of the mm dust distribution in protoplanetary discs are now showing a ubiquity of annular gap and ring substructures. Their identification and accurate characterization is critical to probing the physical processes responsible. We present Frankenstein (frank), an open source code that recovers axisymmetric disc structures at sub-beam resolution. By fitting the visibilities directly, the model reconstructs a disc's 1D radial brightness profile nonparametrically using a fast (<~1 min) Gaussian process. The code avoids limitations of current methods that obtain the radial brightness profile by either extracting it from the disc image via nonlinear deconvolution at the cost of reduced fit resolution, or by assumptions placed on the functional forms of disc structures to fit the visibilities parametrically. We use mock ALMA observations to quantify the method's intrinsic capability and its performance as a function of baseline-dependent signal-to-noise. Comparing the technique to profile extraction from a CLEAN image, we motivate how our fits accurately recover disc structures at a sub-beam resolution. Demonstrating the model's utility in fitting real high and moderate resolution observations, we conclude by proposing applications to address open questions on protoplanetary disc structure and processes. △ Less

Submitted 28 May, 2020; v1 submitted 15 May, 2020; originally announced May 2020.

Comments: 27 pages, 17 figures, accepted in MNRAS

Constraints on primordial gravitational waves from the Cosmic Microwave Background

Authors: Thomas J. Clarke, Edmund J. Copeland, Adam Moss

Abstract: Searches for primordial gravitational waves have resulted in constraints in a large frequency range from a variety of sources. The standard Cosmic Microwave Background (CMB) technique is to parameterise the tensor power spectrum in terms of the tensor-to-scalar ratio, $r$, and spectral index, $n_{\rm t}$, and constrain these using measurements of the temperature and polarization power spectra. Ano… ▽ More Searches for primordial gravitational waves have resulted in constraints in a large frequency range from a variety of sources. The standard Cosmic Microwave Background (CMB) technique is to parameterise the tensor power spectrum in terms of the tensor-to-scalar ratio, $r$, and spectral index, $n_{\rm t}$, and constrain these using measurements of the temperature and polarization power spectra. Another method, applicable to modes well inside the cosmological horizon at recombination, uses the shortwave approximation, under which gravitational waves behave as an effective neutrino species. In this paper we give model-independent CMB constraints on the energy density of gravitational waves, $Ω_\textrm{gw} h^2$, for the entire range of observable frequencies. On large scales, $f \lesssim 10^{-16}\, \text{Hz}$, we reconstruct the initial tensor power spectrum in logarithmic frequency bins, finding maximal sensitivity for scales close to the horizon size at recombination. On small scales, $f \gtrsim10^{-15}\,\mbox{Hz}$, we use the shortwave approximation, finding $Ω_\textrm{gw} h^2 < 1.7 \times10^{-6}$ for adiabatic initial conditions and $Ω_\textrm{gw} h^2 < 2.9 \times10^{-7}$ for homogeneous initial conditions (both $2σ$ upper limits). For scales close to the horizon size at recombination, we use second-order perturbation theory to calculate the back-reaction from gravitational waves, finding $Ω_\textrm{gw} h^2 < 8.4 \times10^{-7}$, in the absence of neutrino anisotropic stress and $Ω_\textrm{gw} h^2 < 8.6 \times10^{-7}$ when including neutrino anisotropic stress. These constraints are valid for $ 10^{-15}\, \text{Hz} \gtrsim f \gtrsim 3 \times 10^{-16}\, \text{Hz}$. △ Less

Submitted 13 July, 2020; v1 submitted 23 April, 2020; originally announced April 2020.

Comments: 32 pages, 9 figures, 1 appendix (with 1 figure). Prepared for submission to JCAP