Showing 1–50 of 967 results for author: Wright, J

Searching for GEMS: TOI-5688 A b, a low-density giant orbiting a high-metallicity early M-dwarf

Authors: Varghese Reji, Shubham Kanodia, Joe Ninan, Caleb I. Cañas, Jessica Libby-Roberts, Andrea S. J. Lin, Arvind F Gupta, Tera N. Sewaby, Alexander Larsen, Henry A. Kobulnicky, Philip I. Choi, Nez Evans, Sage Santomenna, Isabelle Winnick, Larry Yu, Jaime A. Alvarado-Montes, Chad Bender, Lia Marta Bernabò, Cullen H. Blake, William D. Cochran, Scott A. Diddams, Samuel Halverson, Te Han, Fred Hearty, Sarah E. Logsdon , et al. (9 additional authors not shown)

Abstract: We present the discovery of a low-density planet transiting TOI-5688 A b, a high-metallicity M2V star. This planet was discovered as part of the search for transiting giant planets ($R \gtrsim8$ M$_\oplus$) through the Searching for GEMS (Giant Exoplanets around M-dwarf Stars) survey. The planet TOI-5688 A b was discovered with the Transiting Exoplanet Survey Satellite (TESS), and characterized wi… ▽ More We present the discovery of a low-density planet transiting TOI-5688 A b, a high-metallicity M2V star. This planet was discovered as part of the search for transiting giant planets ($R \gtrsim8$ M$_\oplus$) through the Searching for GEMS (Giant Exoplanets around M-dwarf Stars) survey. The planet TOI-5688 A b was discovered with the Transiting Exoplanet Survey Satellite (TESS), and characterized with ground-based transits from Red Buttes Observatory (RBO), the Table Mountain Observatory of Pomona College, and radial velocity (RV) measurements with the Habitable-Zone Planet Finder (HPF) on the 10 m Hobby Eberly Telescope (HET) and NEID on the WIYN 3.5 m telescope. From the joint fit of transit and RV data, the mass of the planet is $124\pm24$ M$_\oplus$ and the radius is $10.4\pm0.7$ R$_\oplus$. This planet has a density of $0.61^{+0.20}_{-0.15}$ g/cm${}^3$, and is on a $\sim2.95$ day orbit around its host star. The spectroscopic and photometric analysis of the host star TOI-5688 A shows that it is a high metallicity ([Fe/H] $ = 0.47\pm0.16$ dex) M2V star, favoring the core-accretion formation pathway as the likely formation scenario for this planet. In this paper, we analyze potential mechanisms of planet formation in the context of the formation of TOI-5688 A b. Additionally, observations with Gaia suggest the presence of a wide-separation binary companion, TOI-5688 B, which has a projected separation of $\sim5"$ (1110 AU) and is an M4V. This makes TOI-5688 A b part of a growing number of GEMS in wide-separation binary systems. △ Less

Submitted 4 September, 2024; v1 submitted 2 September, 2024; originally announced September 2024.

Comments: 20 pages, 7 figures, Submitted to AJ, Comments are welcome

Virtual EVE: a Deep Learning Model for Solar Irradiance Prediction

Authors: Manuel Indaco, Daniel Gass, William James Fawcett, Richard Galvez, Paul J. Wright, Andrés Muñoz-Jaramillo

Abstract: Understanding space weather is vital for the protection of our terrestrial and space infrastructure. In order to predict space weather accurately, large amounts of data are required, particularly in the extreme ultraviolet (EUV) spectrum. An exquisite source of information for such data is provided by the Solar Dynamic Observatory (SDO), which has been gathering solar measurements for the past 13… ▽ More Understanding space weather is vital for the protection of our terrestrial and space infrastructure. In order to predict space weather accurately, large amounts of data are required, particularly in the extreme ultraviolet (EUV) spectrum. An exquisite source of information for such data is provided by the Solar Dynamic Observatory (SDO), which has been gathering solar measurements for the past 13 years. However, after a malfunction in 2014 affecting the onboard Multiple EUV Grating Spectrograph A (MEGS-A) instrument, the scientific output in terms of EUV measurements has been significantly degraded. Building upon existing research, we propose to utilize deep learning for the virtualization of the defective instrument. Our architecture features a linear component and a convolutional neural network (CNN) -- with EfficientNet as a backbone. The architecture utilizes as input grayscale images of the Sun at multiple frequencies -- provided by the Atmospheric Imaging Assembly (AIA) -- as well as solar magnetograms produced by the Helioseismic and Magnetic Imager (HMI). Our findings highlight how AIA data are all that is needed for accurate predictions of solar irradiance. Additionally, our model constitutes an improvement with respect to the state-of-the-art in the field, further promoting the idea of deep learning as a viable option for the virtualization of scientific instruments. △ Less

Submitted 30 August, 2024; originally announced August 2024.

Comments: Submission to Machine Learning and the Physical Sciences Workshop, NeurIPS 2023

Two-neutrino double electron capture of $^{124}$Xe in the first LUX-ZEPLIN exposure

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, J. W. Bargemann, E. E. Barillier, K. Beattie, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. Bishop, G. M. Blockinger, B. Boxer, C. A. J. Brew , et al. (180 additional authors not shown)

Abstract: The broad physics reach of the LUX-ZEPLIN (LZ) experiment covers rare phenomena beyond the direct detection of dark matter. We report precise measurements of the extremely rare decay of $^{124}$Xe through the process of two-neutrino double electron capture (2$ν$2EC), utilizing a $1.39\,\mathrm{kg} \times \mathrm{yr}$ isotopic exposure from the first LZ science run. A half-life of… ▽ More The broad physics reach of the LUX-ZEPLIN (LZ) experiment covers rare phenomena beyond the direct detection of dark matter. We report precise measurements of the extremely rare decay of $^{124}$Xe through the process of two-neutrino double electron capture (2$ν$2EC), utilizing a $1.39\,\mathrm{kg} \times \mathrm{yr}$ isotopic exposure from the first LZ science run. A half-life of $T_{1/2}^{2\nu2\mathrm{EC}} = (1.09 \pm 0.14_{\text{stat}} \pm 0.05_{\text{sys}}) \times 10^{22}\,\mathrm{yr}$ is observed with a statistical significance of $8.3\,σ$, in agreement with literature. First empirical measurements of the KK capture fraction relative to other K-shell modes were conducted, and demonstrate consistency with respect to recent signal models at the $1.4\,σ$ level. △ Less

Submitted 30 August, 2024; originally announced August 2024.

Comments: 15 pages, 3 figures

Data-Driven Modeling of Telluric Features and Stellar Variability with StellarSpectraObservationFitting.jl

Authors: Christian Gilbertson, Eric B. Ford, Samuel Halverson, Evan Fitzmaurice, Cullen H. Blake, Guðmundur Stefánsson, Suvrath Mahadevan, Jason T. Wright, Jacob K. Luhn, Joe P. Ninan, Paul Robertson, Arpita Roy, Christian Schwab, Ryan C. Terrien

Abstract: A significant barrier to achieving the radial velocity (RV) measurement accuracy and precision required to characterize terrestrial mass exoplanets is the existence of time-variable features in the measured spectra, from both telluric absorption and stellar variability, which affect measured line shapes and can cause apparent RV shifts. Reaching the desired accuracy using traditional techniques of… ▽ More A significant barrier to achieving the radial velocity (RV) measurement accuracy and precision required to characterize terrestrial mass exoplanets is the existence of time-variable features in the measured spectra, from both telluric absorption and stellar variability, which affect measured line shapes and can cause apparent RV shifts. Reaching the desired accuracy using traditional techniques often requires avoiding lines contaminated by stellar variability and/or changing tellurics, and thus discarding a large fraction of the spectrum, lowering precision. New data-driven methods can help achieve extremely precise and accurate RVs by enabling the use of a larger fraction of the available data. While there exist methods for modeling telluric features or the stellar variability individually, there is a need for additional tools that are capable of modeling them simultaneously at the spectral level. Here we present StellarSpectraObservationFitting.jl (SSOF), a Julia package for measuring Doppler shifts and creating data-driven models (with fast, physically-motivated Gaussian Process regularization) for the time-variable spectral features for both the telluric transmission and stellar spectrum, while accounting for the wavelength-dependent instrumental line-spread function. We demonstrate SSOF's state-of-the-art performance on data from the NEID RV spectrograph on the WIYN 3.5m Telescope for multiple stars. We show SSOF's, ability to accurately identify and characterize spectral variability and provide $\sim$2-6x smaller photon-limited errors over the NEID CCF-based pipeline and match the performance of SERVAL, a leading template-based pipeline, using only observed EPRV spectra. △ Less

Submitted 30 August, 2024; originally announced August 2024.

Comments: 47 pages, 21 figures, 7 tables. In revision for AAS Journals. Code at https://github.com/christiangil/StellarSpectraObservationFitting.jl

Earths within Reach: Evaluation of Strategies for Mitigating Solar Variability using 3.5 years of NEID Sun-as-a-Star Observations

Authors: Eric B. Ford, Chad F. Bender, Cullen H. Blake, Arvind F. Gupta, Shubham Kanodia, Andrea S. J. Lin, Sarah E. Logsdon, Jacob K. Luhn, Suvrath Mahadevan, Michael L. Palumbo III, Ryan C. Terrien, Jason T. Wright, Jinglin Zhao, Samuel Halverson, Emily Hunting, Paul Robertson, Arpita Roy, Gudmundur Stefansson

Abstract: We present the results of Sun-as-a-star observations by the NEID Solar Telescope at WIYN Observatory, spanning January 1, 2021 through June 30, 2024. We identify 117,060 observations which are unlikely to be significantly affected by weather, hardware or major calibration issues. We describe several high-level data products being made available to the community to aid in the interpretation and int… ▽ More We present the results of Sun-as-a-star observations by the NEID Solar Telescope at WIYN Observatory, spanning January 1, 2021 through June 30, 2024. We identify 117,060 observations which are unlikely to be significantly affected by weather, hardware or major calibration issues. We describe several high-level data products being made available to the community to aid in the interpretation and inter comparisons of NEID solar observations. Solar observations demonstrate excellent performance of NEID, including radial velocity (RV) accuracy and long-term stability of better than $\simeq 0.37$ m s$^{-1}$ over $\simeq 3.5$ years, even though NEID was not originally designed or optimized for daytime observations of the Sun. Currently, intrinsic stellar variability is the primary barrier to detecting Earth-analog planets for most nearby, Sun-like stars. We present a comparison of the effectiveness of several methods proposed to mitigate the effects of solar variability on the Sun's estimated RV. We find that the Scalpels algorithm performs particularly well and substantially reduces the RMS RV of solar spectra from over 2 m s$^{-1}$ to 0.277 m s$^{-1}$. Even when training on a subset of days with NEID solar observations and testing on a held-out sample, the RMS of cleaned RV is 0.34-0.42 m s$^{-1}$. This is significantly better than previous attempts at removing solar variability and suggests that the current generation of EPRV instruments are technically capable of detecting Earth-mass planets orbiting a solar twin if provided with sufficient observing time allocations ($\sim 10^3$ nights of observations). △ Less

Submitted 23 August, 2024; originally announced August 2024.

Comments: 25 pages, 14 figures. Submitted to AAS Journals. Data release archived at https://zenodo.org/doi/10.5281/zenodo.13363761

EWOCS-II: X-ray properties of the Wolf-Rayet stars in the young Galactic super star cluster Westerlund 1

Authors: K. Anastasopoulou, M. G. Guarcello, E. Flaccomio, S. Sciortino, S. Benatti, M. De Becker, N. J. Wright, J. Drake, J. F. Albacete-Colombo, M. Andersen, C. Argiroffi, A. Bayo, R. Castellanos, M. Gennaro, E. K. Grebel, M. Miceli, F. Najarro, I. Negueruela, L. Prisinzano, B. Ritchie, M. Robberto, E. Sabbi, P. Zeidler

Abstract: We present the most comprehensive and deepest X-ray study to date of the properties of the richest Wolf-Rayet (WR) population observed in a single stellar cluster, Westerlund 1 (Wd1). This work is based on 36 Chandra observations obtained from the "Extended Westerlund 1 and 2 Open Clusters Survey" (EWOCS) project, plus 8 archival Chandra observations. The overall exposure depth (~1.1 Ms) and basel… ▽ More We present the most comprehensive and deepest X-ray study to date of the properties of the richest Wolf-Rayet (WR) population observed in a single stellar cluster, Westerlund 1 (Wd1). This work is based on 36 Chandra observations obtained from the "Extended Westerlund 1 and 2 Open Clusters Survey" (EWOCS) project, plus 8 archival Chandra observations. The overall exposure depth (~1.1 Ms) and baseline of the EWOCS observations extending over more than one year enable us to perform a detailed photometric, colour, and spectral analysis, as well as to search for short- and long-term periodicity. In X-rays, we detect 20 out of the 24 known Wolf-Rayet stars in Wd1 down to an observed luminosity of ~7$\times10^{29}$erg s$^{-1}$ (assuming a distance of 4.23 kpc to Wd1), with 8 WR stars being detected in X-rays for the first time. Nine stars show clear evidence of variability over the year-long baseline, with clear signs of periodicity. The X-ray colours and spectral analysis reveal that the vast majority of the WR stars are hard X-ray sources (kT$\geq$2.0keV). The Fe XXV emission line at ~6.7 keV, which commonly originates from the wind-wind collision zone in binary systems, is detected for the first time in the spectra of 17 WR stars in Wd1. In addition the ~6.4 keV fluorescent line is observed in the spectra of three stars, indicating that dense cold material coexists with the hot gas in these systems. Overall, our X-ray results alone suggest a very high binary fraction ($\geq$80%) for the WR star population in Wd1. When combining our results with properties of the WR population from other wavelengths, we estimate a binary fraction of $\geq$92%, which could even reach unity. This suggests that either all the most massive stars are found in binary systems within Wd1, or that binarity is essential for the formation of such a rich population of WR stars. △ Less

Submitted 20 August, 2024; originally announced August 2024.

Comments: 20 pages, 11 figures, 6 tables, accepted for publication in A&A

A reassessment of the "hard-steps" model for the evolution of intelligent life

Authors: Daniel B. Mills, Jennifer L. Macalady, Adam Frank, Jason T. Wright

Abstract: According to the "hard-steps" model, the origin of humanity required "successful passage through a number of intermediate steps" (so-called "hard" or "critical" steps) that were intrinsically improbable with respect to the total time available for biological evolution on Earth. This model similarly predicts that technological life analogous to human life on Earth is "exceedingly rare" in the unive… ▽ More According to the "hard-steps" model, the origin of humanity required "successful passage through a number of intermediate steps" (so-called "hard" or "critical" steps) that were intrinsically improbable with respect to the total time available for biological evolution on Earth. This model similarly predicts that technological life analogous to human life on Earth is "exceedingly rare" in the universe. Here, we critically reevaluate the core assumptions of the hard-steps model in light of recent advances in the Earth and life sciences. Specifically, we advance a potential alternative model where there are no hard steps, and evolutionary novelties (or singularities) required for human origins can be explained via mechanisms outside of intrinsic improbability. Furthermore, if Earth's surface environment was initially inhospitable not only to human life, but also to certain key intermediate steps in human evolution (e.g., the origin of eukaryotic cells, multicellular animals), then the "delay" in the appearance of humans can be best explained through the sequential opening of new global environmental windows of habitability over Earth history, with humanity arising relatively quickly once the right conditions were established. In this co-evolutionary (or geobiological) scenario, humans did not evolve "early" or "late" with respect to the total lifespan of the biosphere, but "on time." △ Less

Submitted 19 August, 2024; originally announced August 2024.

Apostle--Auriga: Effects of stellar feedback subgrid models on the evolution of angular momentum in disc galaxies

Authors: Hang Yang, Shihong Liao, Azadeh Fattahi, Carlos S. Frenk, Liang Gao, Qi Guo, Shi Shao, Lan Wang, Ruby J. Wright, Guangquan Zeng

Abstract: Utilizing the Apostle--Auriga simulations, which start from the same zoom-in initial conditions of Local Group-like systems but run with different galaxy formation subgrid models and hydrodynamic solvers, we study the impact of stellar feedback models on the evolution of angular momentum in disc galaxies. At $z = 0$, Auriga disc galaxies tend to exhibit higher specific angular momenta compared to… ▽ More Utilizing the Apostle--Auriga simulations, which start from the same zoom-in initial conditions of Local Group-like systems but run with different galaxy formation subgrid models and hydrodynamic solvers, we study the impact of stellar feedback models on the evolution of angular momentum in disc galaxies. At $z = 0$, Auriga disc galaxies tend to exhibit higher specific angular momenta compared to their cross-matched Apostle counterparts. By tracing the evolution history of the Lagrangian mass tracers of the in-situ star particles in the $z = 0$ galaxies, we find that the specific angular momentum distributions of the gas tracers from the two simulations at the halo accretion time are relatively similar. The present-day angular momentum difference is mainly driven by the physical processes occurring inside dark matter haloes, especially galactic fountains. Due to the different subgrid implementations of stellar feedback processes, Auriga galaxies contain a high fraction of recycled gas tracers (${\sim} 65$ per cent) which could acquire angular momentum through mixing with the high angular momentum circumgalactic medium (CGM). In Apostle, however, the fraction of recycled gas tracers is significantly lower (down to ${\sim} 20$ per cent for Milky Way-sized galaxies) and the angular momentum acquisition from the CGM is marginal. As a result, the present-day Auriga galaxies overall have higher specific angular momenta. △ Less

Submitted 19 August, 2024; originally announced August 2024.

Comments: 12 pages, 8 figures, submitted to MNRAS. comments are welcome

Me want cookie! Towards automated and transparent data governance on the Web

Authors: Jesse Wright, Beatriz Esteves, Rui Zhao

Abstract: This paper presents a sociotechnical vision for managing personal data, including cookies, within Web browsers. We first present our vision for a future of semi-automated data governance on the Web, using policy languages to describe data terms of use, and having browsers act on behalf of users to enact policy-based controls. Then, we present an overview of the technical research required to {prov… ▽ More This paper presents a sociotechnical vision for managing personal data, including cookies, within Web browsers. We first present our vision for a future of semi-automated data governance on the Web, using policy languages to describe data terms of use, and having browsers act on behalf of users to enact policy-based controls. Then, we present an overview of the technical research required to {prove} that existing policy languages express a sufficient range of concepts for describing cookie policies on the Web today. We view this work as a stepping stone towards a future of semi-automated data governance at Web-scale, which in the long term will also be used by next-generation Web technologies such as Web agents and Solid. △ Less

Submitted 16 August, 2024; originally announced August 2024.

Comments: Submitted to "NeXt-generation Data Governance workshop 2024"; available on OpenReview at https://openreview.net/forum?id=Bhia6mPaCF

Lightning declines over shipping lanes following regulation of fuel sulfur emissions

Authors: Chris J. Wright, Joel A. Thornton, Lyatt Jaeglé, Yang Cao, Yannian Zhu, Jihu Liu, Randall Jones II, Robert H Holzworth, Daniel Rosenfeld, Robert Wood, Peter Blossey, Daehyun Kim

Abstract: Aerosol interactions with clouds represent a significant uncertainty in our understanding of the Earth system. Deep convective clouds may respond to aerosol perturbations in several ways that have proven difficult to elucidate with observations. Here, we leverage the two busiest maritime shipping lanes in the world, which emit aerosol particles and their precursors into an otherwise relatively cle… ▽ More Aerosol interactions with clouds represent a significant uncertainty in our understanding of the Earth system. Deep convective clouds may respond to aerosol perturbations in several ways that have proven difficult to elucidate with observations. Here, we leverage the two busiest maritime shipping lanes in the world, which emit aerosol particles and their precursors into an otherwise relatively clean tropical marine boundary layer, to make headway on the influence of aerosol on deep convective clouds. The recent seven-fold change in allowable fuel sulfur by the International Maritime Organization allows us to test the sensitivity of the lightning to changes in ship plume aerosol size distributions. We find that, across a range of atmospheric thermodynamic conditions, the previously documented enhancement of lightning over the shipping lanes has fallen by over 40\%. The enhancement is therefore at least partially aerosol-mediated, a conclusion that is supported by observations of droplet number at cloud base, which show a similar decline over the shipping lane. These results have fundamental implications for our understanding of aerosol-cloud interactions, suggesting that deep convective clouds are impacted by the aerosol number distribution in the remote marine environment. △ Less

Submitted 13 August, 2024; originally announced August 2024.

Quiet Please: Detrending Radial Velocity Variations from Stellar Activity with a Physically Motivated Spot Model

Authors: Jared C. Siegel, Samuel Halverson, Jacob K. Luhn, Lily L. Zhao, Khaled Al Moulla, Paul Robertson, Chad F. Bender, Ryan C. Terrien, Arpita Roy, Suvrath Mahadevan, Fred Hearty, Joe P. Ninan, Jason T. Wright, Eric B. Ford, Christian Schwab, Guðmundur Stefánsson, Cullen H. Blake, Michael W. McElwain

Abstract: For solar-type stars, spots and their associated magnetic regions induce radial velocity perturbations through the Doppler rotation signal and the suppression of convective blueshift -- collectively known as rotation-modulation. We developed the Rotation-Convection (RC) model: a method of detrending and characterizing rotation-modulation, using only cross-correlation functions or 1-dimensional spe… ▽ More For solar-type stars, spots and their associated magnetic regions induce radial velocity perturbations through the Doppler rotation signal and the suppression of convective blueshift -- collectively known as rotation-modulation. We developed the Rotation-Convection (RC) model: a method of detrending and characterizing rotation-modulation, using only cross-correlation functions or 1-dimensional spectra, without the need for continuous high cadence measurements. The RC method uses a simple model for the anomalous radial velocity induced by an active region and has two inputs: stellar flux (or a flux proxy) and the relative radial velocity between strongly and weakly absorbed wavelengths (analogous to the bisector-inverse slope). On NEID solar data (three month baseline), the RC model lowers the amplitude of rotationally-modulated stellar activity to below the meter-per-second level. For the standard star HD 26965, the RC model detrends the activity signal to the meter-per-second level for HARPS, EXPRES, and NEID observations, even though the temporal density and timespan of the observations differs by an order of magnitude between the three datasets. In addition to detrending, the RC model also characterizes the rotation-modulation signal. From comparison with the Solar Dynamics Observatory, we confirmed that the model accurately recovers and separates the rotation and convection radial velocity components. We also mapped the amplitude of the rotation and convection perturbations as a function of height within the stellar atmosphere. Probing stellar atmospheres with our revised spot model will fuel future innovations in stellar activity mitigation, enabling robust exoplanet detection. △ Less

Submitted 13 August, 2024; originally announced August 2024.

Comments: Accepted to AJ, 22 pages, 9 figures

Utilizing Photometry from Multiple Sources to Mitigate Stellar Variability in Precise Radial Velocities: A Case Study of Kepler-21

Authors: Corey Beard, Paul Robertson, Mark R. Giovinazzi, Joseph M. Akana Murphy, Eric B. Ford, Samuel Halverson, Te Han, Rae Holcomb, Jack Lubin, Rafael Luque, Pranav Premnath, Chad F. Bender, Cullen H. Blake, Qian Gong, Howard Isaacson, Shubham Kanodia, Dan Li, Andrea S. J. Lin, 5 Sarah E. Logsdon, Emily Lubar, Michael W. McElwain, Andrew Monson, Joe P. Ninan, Jayadev Rajagopal, Arpita Roy , et al. (4 additional authors not shown)

Abstract: We present a new analysis of Kepler-21, the brightest (V = 8.5) Kepler system with a known transiting exoplanet, Kepler-21 b. Kepler-21 b is a radius valley planet ($R = 1.6\pm 0.2 R_{\oplus}$) with an Earth-like composition (8.38$\pm$1.62 g/cc), though its mass and radius fall in the regime of possible "water worlds." We utilize new Keck/HIRES and WIYN/NEID radial velocity (RV) data in conjunctio… ▽ More We present a new analysis of Kepler-21, the brightest (V = 8.5) Kepler system with a known transiting exoplanet, Kepler-21 b. Kepler-21 b is a radius valley planet ($R = 1.6\pm 0.2 R_{\oplus}$) with an Earth-like composition (8.38$\pm$1.62 g/cc), though its mass and radius fall in the regime of possible "water worlds." We utilize new Keck/HIRES and WIYN/NEID radial velocity (RV) data in conjunction with Kepler and TESS photometry to perform a detailed study of activity mitigation between photometry and RVs. We additionally refine the system parameters, and we utilize Gaia astrometry to place constraints on a long-term RV trend. Our activity analysis affirms the quality of Kepler photometry for removing correlated noise from RVs, despite its temporal distance, though we reveal some cases where TESS may be superior. Using refined orbital parameters and updated composition curves, we rule out a ``water world" scenario for Kepler-21 b, and we identify a long period super-Jupiter planetary candidate, Kepler-21 (c). △ Less

Submitted 5 August, 2024; originally announced August 2024.

Asteroseismology of the Nearby K-Dwarf $σ$ Draconis using the Keck Planet Finder and TESS

Authors: Marc Hon, Daniel Huber, Yaguang Li, Travis S. Metcalfe, Timothy R. Bedding, Joel Ong, Ashley Chontos, Ryan Rubenzahl, Samuel Halverson, Rafael A. García, Hans Kjeldsen, Dennis Stello, Daniel R. Hey, Tiago Campante, Andrew W. Howard, Steven R. Gibson, Kodi Rider, Arpita Roy, Ashley D. Baker, Jerry Edelstein, Chris Smith, Benjamin J. Fulton, Josh Walawender, Max Brodheim, Matt Brown , et al. (54 additional authors not shown)

Abstract: Asteroseismology of dwarf stars cooler than the Sun is very challenging due to the low amplitudes and rapid timescales of oscillations. Here, we present the asteroseismic detection of solar-like oscillations at 4-minute timescales ($ν_{\mathrm{max}}\sim4300μ$Hz) in the nearby K-dwarf $σ$ Draconis using extreme precision Doppler velocity observations from the Keck Planet Finder and 20-second cadenc… ▽ More Asteroseismology of dwarf stars cooler than the Sun is very challenging due to the low amplitudes and rapid timescales of oscillations. Here, we present the asteroseismic detection of solar-like oscillations at 4-minute timescales ($ν_{\mathrm{max}}\sim4300μ$Hz) in the nearby K-dwarf $σ$ Draconis using extreme precision Doppler velocity observations from the Keck Planet Finder and 20-second cadence photometry from NASA's Transiting Exoplanet Survey Satellite. The star is the coolest dwarf star to date with both velocity and luminosity observations of solar-like oscillations, having amplitudes of $5.9\pm0.8\,$cm$\,\text{s}^{-1}$ and $0.8\pm0.2$ ppm, respectively. These measured values are in excellent agreement with established luminosity-velocity amplitude relations for oscillations and provide further evidence that mode amplitudes for stars with $T_{\mathrm{eff}}<\,5500\,$K diminish in scale following a $(L/M)^{1.5}$ relation. By modeling the star's oscillation frequencies from photometric data, we measure an asteroseismic age of $4.5\pm0.9\,\rm{(ran)} \pm 1.2\,\rm{(sys)}$ Gyr. The observations demonstrate the capability of next-generation spectrographs and precise space-based photometry to extend observational asteroseismology to nearby cool dwarfs, which are benchmarks for stellar astrophysics and prime targets for directly imaging planets using future space-based telescopes. △ Less

Submitted 28 August, 2024; v1 submitted 30 July, 2024; originally announced July 2024.

Comments: Accepted for publication in The Astrophysical Journal

Obliquity Constraints for the Extremely Eccentric Sub-Saturn Kepler-1656 b

Authors: Ryan A. Rubenzahl, Andrew W. Howard, Samuel Halverson, Cristobal Petrovich, Isabel Angelo, Guðmundur Stefánsson, Fei Dai, Aaron Householder, Benjamin Fulton, Steven R. Gibson, Arpita Roy, Abby P. Shaum, Howard Isaacson, Max Brodheim, William Deich, Grant M. Hill, Bradford Holden, Daniel Huber, Russ R. Laher, Kyle Lanclos, Joel N. Payne, Erik A. Petigura, Christian Schwab, Josh Walawender, Sharon X. Wang , et al. (3 additional authors not shown)

Abstract: The orbits of close-in exoplanets provide clues to their formation and evolutionary history. Many close-in exoplanets likely formed far out in their protoplanetary disks and migrated to their current orbits, perhaps via high-eccentricity migration (HEM), a process that can also excite obliquities. A handful of known exoplanets are perhaps caught in the act of HEM, as they are observed on highly ec… ▽ More The orbits of close-in exoplanets provide clues to their formation and evolutionary history. Many close-in exoplanets likely formed far out in their protoplanetary disks and migrated to their current orbits, perhaps via high-eccentricity migration (HEM), a process that can also excite obliquities. A handful of known exoplanets are perhaps caught in the act of HEM, as they are observed on highly eccentric orbits with tidal circularization timescales shorter than their ages. One such exoplanet is Kepler-1656 b, which is also the only known non-giant exoplanet (<100 $M_\oplus$) with an extreme eccentricity (e=0.84). We measured the sky-projected obliquity of Kepler-1656 b by observing the Rossiter-McLaughlin effect during a transit with the Keck Planet Finder. Our data are consistent with an aligned orbit, but are also consistent with moderate misalignment with $λ< 50$ deg at 95% confidence, with the most likely solution of $35^{+14.9}_{-21.6}$ deg. A low obliquity would be an unlikely outcome of most eccentricity-exciting scenarios, but we show that the properties of the outer companion in the system are consistent with the coplanar HEM mechanism. Alternatively, if the system is not relatively coplanar (<20 deg mutual inclination), Kepler-1656 b may be presently at a rare snapshot of long-lived eccentricity oscillations that do not induce migration. Kepler-1656 b is only the fourth exoplanet with e>0.8 to have its obliquity constrained; expanding this population will help establish the degree to which orbital misalignment accompanies migration. Future work that constrains the mutual inclinations of outer perturbers will be key for distinguishing plausible mechanisms. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Comments: 11 pages, 5 figures, accepted to ApJL

An Earth-sized Planet on the Verge of Tidal Disruption

Authors: Fei Dai, Andrew W. Howard, Samuel Halverson, Jaume Orell-Miquel, Enric Palle, Howard Isaacson, Benjamin Fulton, Ellen M. Price, Mykhaylo Plotnykov, Leslie A. Rogers, Diana Valencia, Kimberly Paragas, Michael Greklek-McKeon, Jonathan Gomez Barrientos, Heather A. Knutson, Erik A. Petigura, Lauren M. Weiss, Rena Lee, Casey L. Brinkman, Daniel Huber, Gudmundur Steffansson, Kento Masuda, Steven Giacalone, Cicero X. Lu, Edwin S. Kite , et al. (73 additional authors not shown)

Abstract: TOI-6255~b (GJ 4256) is an Earth-sized planet (1.079$\pm0.065$ $R_\oplus$) with an orbital period of only 5.7 hours. With the newly commissioned Keck Planet Finder (KPF) and CARMENES spectrographs, we determined the planet's mass to be 1.44$\pm$0.14 $M_{\oplus}$. The planet is just outside the Roche limit, with $P_{\rm orb}/P_{\rm Roche}$ = 1.13 $\pm0.10$. The strong tidal force likely deforms the… ▽ More TOI-6255~b (GJ 4256) is an Earth-sized planet (1.079$\pm0.065$ $R_\oplus$) with an orbital period of only 5.7 hours. With the newly commissioned Keck Planet Finder (KPF) and CARMENES spectrographs, we determined the planet's mass to be 1.44$\pm$0.14 $M_{\oplus}$. The planet is just outside the Roche limit, with $P_{\rm orb}/P_{\rm Roche}$ = 1.13 $\pm0.10$. The strong tidal force likely deforms the planet into a triaxial ellipsoid with a long axis that is $\sim$10\% longer than the short axis. Assuming a reduced stellar tidal quality factor $Q_\star^\prime \approx10^7$, we predict that tidal orbital decay will cause TOI-6255 to reach the Roche limit in roughly 400 Myr. Such tidal disruptions may produce the possible signatures of planet engulfment that have been on stars with anomalously high refractory elemental abundances compared to its conatal binary companion. TOI-6255 b is also a favorable target for searching for star-planet magnetic interactions, which might cause interior melting and hasten orbital decay. TOI-6255 b is a top target (Emission Spectroscopy Metric of about 24) for phase curve observations with the James Webb Space Telescope. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Comments: 18 pages, 7 figures, 5 tables, accepted to AAS Journals. The first RV mass measurement from the Keck Planet Finder

Swift-BAT GUANO follow-up of gravitational-wave triggers in the third LIGO-Virgo-KAGRA observing run

Authors: Gayathri Raman, Samuele Ronchini, James Delaunay, Aaron Tohuvavohu, Jamie A. Kennea, Tyler Parsotan, Elena Ambrosi, Maria Grazia Bernardini, Sergio Campana, Giancarlo Cusumano, Antonino D'Ai, Paolo D'Avanzo, Valerio D'Elia, Massimiliano De Pasquale, Simone Dichiara, Phil Evans, Dieter Hartmann, Paul Kuin, Andrea Melandri, Paul O'Brien, Julian P. Osborne, Kim Page, David M. Palmer, Boris Sbarufatti, Gianpiero Tagliaferri , et al. (1797 additional authors not shown)

Abstract: We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA (LVK) network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received in low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wav… ▽ More We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA (LVK) network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received in low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wave Transient Catalogs (GWTC-3). Targeted searches were carried out on the entire GW sample using the maximum--likelihood NITRATES pipeline on the BAT data made available via the GUANO infrastructure. We do not detect any significant electromagnetic emission that is temporally and spatially coincident with any of the GW candidates. We report flux upper limits in the 15-350 keV band as a function of sky position for all the catalog candidates. For GW candidates where the Swift-BAT false alarm rate is less than 10$^{-3}$ Hz, we compute the GW--BAT joint false alarm rate. Finally, the derived Swift-BAT upper limits are used to infer constraints on the putative electromagnetic emission associated with binary black hole mergers. △ Less

Submitted 13 July, 2024; originally announced July 2024.

Comments: 50 pages, 10 figures, 4 tables

Private Heterogeneous Federated Learning Without a Trusted Server Revisited: Error-Optimal and Communication-Efficient Algorithms for Convex Losses

Authors: Changyu Gao, Andrew Lowy, Xingyu Zhou, Stephen J. Wright

Abstract: We revisit the problem of federated learning (FL) with private data from people who do not trust the server or other silos/clients. In this context, every silo (e.g. hospital) has data from several people (e.g. patients) and needs to protect the privacy of each person's data (e.g. health records), even if the server and/or other silos try to uncover this data. Inter-Silo Record-Level Differential… ▽ More We revisit the problem of federated learning (FL) with private data from people who do not trust the server or other silos/clients. In this context, every silo (e.g. hospital) has data from several people (e.g. patients) and needs to protect the privacy of each person's data (e.g. health records), even if the server and/or other silos try to uncover this data. Inter-Silo Record-Level Differential Privacy (ISRL-DP) prevents each silo's data from being leaked, by requiring that silo i's communications satisfy item-level differential privacy. Prior work arXiv:2106.09779 characterized the optimal excess risk bounds for ISRL-DP algorithms with homogeneous (i.i.d.) silo data and convex loss functions. However, two important questions were left open: (1) Can the same excess risk bounds be achieved with heterogeneous (non-i.i.d.) silo data? (2) Can the optimal risk bounds be achieved with fewer communication rounds? In this paper, we give positive answers to both questions. We provide novel ISRL-DP FL algorithms that achieve the optimal excess risk bounds in the presence of heterogeneous silo data. Moreover, our algorithms are more communication-efficient than the prior state-of-the-art. For smooth loss functions, our algorithm achieves the optimal excess risk bound and has communication complexity that matches the non-private lower bound. Additionally, our algorithms are more computationally efficient than the previous state-of-the-art. △ Less

Submitted 17 July, 2024; v1 submitted 12 July, 2024; originally announced July 2024.

Comments: The 41st International Conference on Machine Learning (ICML 2024)

Oscillatory integral operators and variable Schrödinger propagators: beyond the universal estimates

Authors: Mingfeng Chen, Shengwen Gan, Shaoming Guo, Jonathan Hickman, Marina Iliopoulou, James Wright

Abstract: We consider a class of Hörmander-type oscillatory integral operators in $\mathbb{R}^n$ for $n \geq 3$ odd with real analytic phase. We derive weak conditions on the phase which ensure $L^p$ bounds beyond the universal $p \geq 2 \cdot \frac{n+1}{n-1}$ range guaranteed by Stein's oscillatory integral theorem. This expands and elucidates pioneering work of Bourgain from the early 1990s. We also consi… ▽ More We consider a class of Hörmander-type oscillatory integral operators in $\mathbb{R}^n$ for $n \geq 3$ odd with real analytic phase. We derive weak conditions on the phase which ensure $L^p$ bounds beyond the universal $p \geq 2 \cdot \frac{n+1}{n-1}$ range guaranteed by Stein's oscillatory integral theorem. This expands and elucidates pioneering work of Bourgain from the early 1990s. We also consider a closely related class of variable coefficient Schrödinger propagator-type operators, and show that the corresponding theory differs significantly from that of the Hörmander-type operators. The main ingredient in the proof is a curved Kakeya/Nikodym maximal function estimate. This is established by combining the polynomial method with certain uniform sublevel set estimates for real analytic functions. The sublevel set estimates are the main novelty in the argument and can be interpreted as a form of quantification of linear independence in the $C^ω$ category. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 80 pages

MSC Class: 42B20

Does the Fundamental Metallicity Relation Evolve with Redshift? II: The Evolution in Normalisation of the Mass-Metallicity Relation

Authors: Alex M. Garcia, Paul Torrey, Sara L. Ellison, Kathryn Grasha, Qian-Hui Chen, Z. S. Hemler, Dhruv T. Zimmerman, Ruby J. Wright, Henry R. M. Zovaro, Erica J. Nelson, Ryan L. Sanders, Lisa J. Kewley, Lars Hernquist

Abstract: The metal content of galaxies is a direct probe of the baryon cycle. A hallmark example is the relationship between a galaxy's stellar mass, star formation rate (SFR), and gas-phase metallicity: the Fundamental Metallicity Relation (FMR). While low-redshift ($z\lesssim4$) observational studies suggest that the FMR is redshift-invariant, recent JWST data indicate deviations from this model. In this… ▽ More The metal content of galaxies is a direct probe of the baryon cycle. A hallmark example is the relationship between a galaxy's stellar mass, star formation rate (SFR), and gas-phase metallicity: the Fundamental Metallicity Relation (FMR). While low-redshift ($z\lesssim4$) observational studies suggest that the FMR is redshift-invariant, recent JWST data indicate deviations from this model. In this study, we utilize the FMR to predict the evolution of the normalisation of the mass-metallicity relation (MZR) using the cosmological simulations Illustris, IllustrisTNG, EAGLE, and SIMBA. Our findings demonstrate that a $z = 0$ calibrated FMR struggles to predict the evolution in the MZR of each simulation. To quantify the divergence of the predictions, we introduce the concepts of a ''static'' FMR, where the role of the SFR in setting the normalization of the MZR does not change with redshift, and a ''dynamic'' FMR, where the role of SFR evolves over time. We find static FMRs in Illustris and SIMBA and dynamic FMRs in IllustrisTNG and EAGLE. We suggest that the differences between these models likely points to the subtle differences in the implementation of the baryon cycle. Moreover, we echo recent JWST results at $z > 4$ by finding significant offsets from the FMR in IllustrisTNG and EAGLE, suggesting that the observed FMR may be dynamic as well. Overall, our findings imply that the current FMR framework neglects important variations in the baryon cycle through cosmic time. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted to MNRAS. 17 pages, 8 Figures. See also paper one of series (arXiv:2403.08856). Comments welcome

Surviving in the Hot Neptune Desert: The Discovery of the Ultra-Hot Neptune TOI-3261b

Authors: Emma Nabbie, Chelsea X. Huang, Jennifer A. Burt, David J. Armstrong, Eric E. Mamajek, Vardan Adibekyan, Sérgio G. Sousa, Eric D. Lopez, Daniel P. Thorngren, Jorge Fernández, Gongjie Li, James S. Jenkins, Jose I. Vines, João Gomes da Silva, Robert A. Wittenmyer, Daniel Bayliss, César Briceño, Karen A. Collins, Xavier Dumusque, Keith D. Horne, Marcelo F. Keniger, Nicholas Law, Jorge Lillo-Box, Shang-Fei Liu, Andrew W. Mann , et al. (23 additional authors not shown)

Abstract: The recent discoveries of Neptune-sized ultra-short period planets (USPs) challenge existing planet formation theories. It is unclear whether these residents of the Hot Neptune Desert have similar origins to smaller, rocky USPs, or if this discrete population is evidence of a different formation pathway altogether. We report the discovery of TOI-3261b, an ultra-hot Neptune with an orbital period… ▽ More The recent discoveries of Neptune-sized ultra-short period planets (USPs) challenge existing planet formation theories. It is unclear whether these residents of the Hot Neptune Desert have similar origins to smaller, rocky USPs, or if this discrete population is evidence of a different formation pathway altogether. We report the discovery of TOI-3261b, an ultra-hot Neptune with an orbital period $P$ = 0.88 days. The host star is a $V = 13.2$ magnitude, slightly super-solar metallicity ([Fe/H] $\simeq$ 0.15), inactive K1.5 main sequence star at $d = 300$ pc. Using data from the Transiting Exoplanet Survey Satellite and the Las Cumbres Observatory Global Telescope, we find that TOI-3261b has a radius of $3.82_{-0.35}^{+0.42}$ $R_{\oplus}$. Moreover, radial velocities from ESPRESSO and HARPS reveal a mass of $30.3_{-2.4}^{+2.2}$ $M_{\oplus}$, more than twice the median mass of Neptune-sized planets on longer orbits. We investigate multiple mechanisms of mass loss that can reproduce the current-day properties of TOI-3261b, simulating the evolution of the planet via tidal stripping and photoevaporation. Thermal evolution models suggest that TOI-3261b should retain an envelope potentially enriched with volatiles constituting $\sim$5% of its total mass. This is the second highest envelope mass fraction among ultra-hot Neptunes discovered to date, making TOI-3261b an ideal candidate for atmospheric follow-up observations. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: 20 pages, 11 figures, accepted to AJ

Gaia dr3 reveals the complex dynamical evolution within star clusters

Authors: C. Viscasillas Vázquez, L. Magrini, N. Miret-Roig, N. J. Wright, J. Alves, L. Spina, R. P. Church, G. Tautvaišienė, S. Randich

Abstract: Context: Star clusters, composed of stars born from the same molecular cloud, serve as invaluable natural laboratories for understanding the fundamental processes governing stellar formation and evolution. Aims: This study aims to investigate correlations between the Mean Interdistance ($\bar{D_{i}}$), Mean Closest Interdistance ($\bar{D_{c}}$) and Median Weighted Central Interdistance (… ▽ More Context: Star clusters, composed of stars born from the same molecular cloud, serve as invaluable natural laboratories for understanding the fundamental processes governing stellar formation and evolution. Aims: This study aims to investigate correlations between the Mean Interdistance ($\bar{D_{i}}$), Mean Closest Interdistance ($\bar{D_{c}}$) and Median Weighted Central Interdistance ($\bar{D_{cc}}$) with the age of star clusters, examining their evolutionary trends and assessing the robustness of these quantities as possible age indicators. Methods: We selected a sample of open clusters in the solar region and with a representative number of members (e.g. well populated and without outliers). The interdistances are derived from the spatial distribution of member stars within a cluster. Their evolution over time allows us to use them as an age indicators for star clusters. Results: Our investigation reveals a high-significant correlation between the interdistances and cluster age. Considering the full sample of clusters between 7 and 9 kpc, the relationship is very broad. This is due to uncertainties in parallax, which increase with increasing distance. In particular, we must limit the sample to a maximum distance from the Sun of about 200 pc to avoid artificial effects on cluster shape and on the spatial distribution of their stars along the line of sight. Conclusions: By conservatively restraining the distance to a maximum of $\sim$200 pc, we have established a relationship between the interdistances and the age of the clusters. In our sample, the relationship is mainly driven by the internal expansion of the clusters, and is marginally affected by external perturbative effects. Such relation might enhance our comprehension of cluster dynamics and might be used to derive cluster dynamical ages. △ Less

Submitted 30 June, 2024; originally announced July 2024.

Comments: This paper was accepted for publication in Astronomy & Astrophysics on June 28, 2024. The document includes 9 main figures and 4 additional figures in the appendix. Additionally, it contains one main table and an extract of another table in the appendix, available at the CDS

Local structure maturation in high entropy oxide (Mg,Co,Ni,Cu,Zn)1-x(Cr,Mn)xO thin films

Authors: Gabriela E. Niculescu, Gerald R. Bejger, John P. Barber, Joshua T. Wright, Saeed S. I. Almishal, Matthew Webb, Sai Venkata Gayathri Ayyagari, Jon-Paul Maria, Nasim Alem, John T. Heron, Christina M. Rost

Abstract: High entropy oxides (HEO)s have garnered much interest due to their available high degree of tunability. Here, we study the local structure of (MgNiCuCoZn)0.167(MnCr)0.083O, a composition based on the parent HEO (MgNiCuCoZn)0.2O.We synthesized a series of thin films via pulsed laser deposition at incremental oxygen partial pressures. X-ray diffraction shows lattice parameter to decrease with incre… ▽ More High entropy oxides (HEO)s have garnered much interest due to their available high degree of tunability. Here, we study the local structure of (MgNiCuCoZn)0.167(MnCr)0.083O, a composition based on the parent HEO (MgNiCuCoZn)0.2O.We synthesized a series of thin films via pulsed laser deposition at incremental oxygen partial pressures. X-ray diffraction shows lattice parameter to decrease with increased pO2 pressures until the onset of phase separation. X-ray absorption fine structure shows that specific atomic species in the composition dictate the global structure of the material as Cr, Co, and Mn shift to energetically favorable coordination with increasing pressure. Transmission electron microscopy analysis on a lower-pressure sample exhibits a rock salt structure, but the higher-pressure sample reveals reflections reminiscent of the spinel structure. In all, these findings give a more complete picture on how (MgNiCuCoZn)0.167(MnCr)0.083O forms with varying initial conditions and advances fundamental knowledge of cation behavior in high entropy oxides. △ Less

Submitted 19 June, 2024; originally announced June 2024.

The Design, Implementation, and Performance of the LZ Calibration Systems

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, E. E. Barillier, J. W. Bargemann, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. Bishop, G. M. Blockinger, B. Boxer , et al. (179 additional authors not shown)

Abstract: LUX-ZEPLIN (LZ) is a tonne-scale experiment searching for direct dark matter interactions and other rare events. It is located at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. The core of the LZ detector is a dual-phase xenon time projection chamber (TPC), designed with the primary goal of detecting Weakly Interacting Massive Particles (WIMPs) via their induced low e… ▽ More LUX-ZEPLIN (LZ) is a tonne-scale experiment searching for direct dark matter interactions and other rare events. It is located at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. The core of the LZ detector is a dual-phase xenon time projection chamber (TPC), designed with the primary goal of detecting Weakly Interacting Massive Particles (WIMPs) via their induced low energy nuclear recoils. Surrounding the TPC, two veto detectors immersed in an ultra-pure water tank enable reducing background events to enhance the discovery potential. Intricate calibration systems are purposely designed to precisely understand the responses of these three detector volumes to various types of particle interactions and to demonstrate LZ's ability to discriminate between signals and backgrounds. In this paper, we present a comprehensive discussion of the key features, requirements, and performance of the LZ calibration systems, which play a crucial role in enabling LZ's WIMP-search and its broad science program. The thorough description of these calibration systems, with an emphasis on their novel aspects, is valuable for future calibration efforts in direct dark matter and other rare-event search experiments. △ Less

Submitted 5 September, 2024; v1 submitted 2 May, 2024; originally announced June 2024.

Journal ref: JINST 19 P08027 (2024)

HD 21520 b: a warm sub-Neptune transiting a bright G dwarf

Authors: Molly Nies, Ismael Mireles, François Bouchy, Diana Dragomir, Belinda A. Nicholson, Nora L. Eisner, Sergio G. Sousa, Karen A. Collins, Steve B. Howell, Carl Ziegler, Coel Hellier, Brett Addison, Sarah Ballard, Brendan P. Bowler, César Briceño, Catherine A. Clark, Dennis M. Conti, Xavier Dumusque, Billy Edwards, Crystal L. Gnilka, Melissa Hobson, Jonathan Horner, Stephen R. Kane, John Kielkopf, Baptiste Lavie , et al. (27 additional authors not shown)

Abstract: We report the discovery and validation of HD 21520 b, a transiting planet found with TESS and orbiting a bright G dwarf (V=9.2, $T_{eff} = 5871 \pm 62$ K, $R_{\star} = 1.04\pm 0.02\, R_{\odot}$). HD 21520 b was originally alerted as a system (TOI-4320) consisting of two planet candidates with periods of 703.6 and 46.4 days. However, our analysis supports instead a single-planet system with an orbi… ▽ More We report the discovery and validation of HD 21520 b, a transiting planet found with TESS and orbiting a bright G dwarf (V=9.2, $T_{eff} = 5871 \pm 62$ K, $R_{\star} = 1.04\pm 0.02\, R_{\odot}$). HD 21520 b was originally alerted as a system (TOI-4320) consisting of two planet candidates with periods of 703.6 and 46.4 days. However, our analysis supports instead a single-planet system with an orbital period of $25.1292\pm0.0001$ days and radius of $2.70 \pm 0.09\, R_{\oplus}$. Three full transits in sectors 4, 30 and 31 match this period and have transit depths and durations in agreement with each other, as does a partial transit in sector 3. We also observe transits using CHEOPS and LCOGT. SOAR and Gemini high-resolution imaging do not indicate the presence of any nearby companions, and MINERVA-Australis and CORALIE radial velocities rule out an on-target spectroscopic binary. Additionally, we use ESPRESSO radial velocities to obtain a tentative mass measurement of $7.9^{+3.2}_{-3.0}\, M_{\oplus}$, with a 3-$σ$ upper limit of 17.7 $M_{\oplus}$. Due to the bright nature of its host and likely significant gas envelope of the planet, HD 21520 b is a promising candidate for further mass measurements and for atmospheric characterization. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: Submitted to MNRAS

Trials and Tribulations in the Reanalysis of KELT-24 b: a Case Study for the Importance of Stellar Modeling

Authors: Mark R. Giovinazzi, Bryson Cale, Jason D. Eastman, Joseph E. Rodriguez, Cullen H. Blake, Keivan G. Stassun, Thomas G. Beatty, Nate McCrady, Andrew Vanderburg, Michelle Kunimoto, Adam L. Kraus, Joseph Twicken, Cayla M. Dedrick, Jonathan Horner, John A. Johnson, Samson A. Johnson, Peter Plavchan, David H. Sliski, Maurice L. Wilson, Robert A. Wittenmyer, Jason T. Wright, Marshall C. Johnson, Mark E. Rose, Matthew Cornachione

Abstract: We present a new analysis of the KELT-24 system, comprising a well-aligned hot Jupiter, KELT-24~b, and a bright ($V=8.3$), nearby ($d=96.9~\mathrm{pc}$) F-type host star. KELT-24~b was independently discovered by two groups in 2019, with each reporting best-fit stellar parameters that were notably inconsistent. Here, we present three independent analyses of the KELT-24 system, each incorporating a… ▽ More We present a new analysis of the KELT-24 system, comprising a well-aligned hot Jupiter, KELT-24~b, and a bright ($V=8.3$), nearby ($d=96.9~\mathrm{pc}$) F-type host star. KELT-24~b was independently discovered by two groups in 2019, with each reporting best-fit stellar parameters that were notably inconsistent. Here, we present three independent analyses of the KELT-24 system, each incorporating a broad range of photometric and spectroscopic data, including eight sectors of TESS photometry and more than 200 new radial velocities (RVs) from MINERVA. Two of these analyses use KELT-24's observed spectral energy distribution (SED) through a direct comparison to stellar evolutionary models, while our third analysis assumes an unknown additional body contributing to the observed broadband photometry and excludes the SED. Ultimately, we find that the models that include the SED are a poor fit to the available data, so we adopt the system parameters derived without it. We also highlight a single transit-like event observed by TESS, deemed likely to be an eclipsing binary bound to KELT-24, that will require follow-up observations to confirm. We discuss the potential of these additional bodies in the KELT-24 system as a possible explanation for the discrepancies between the results of the different modeling approaches, and explore the system for longer-period planets that may be weakly evident in the RV observations. The comprehensive investigations that we present not only increase the fidelity of our understanding of the KELT-24 system, but also serve as a blueprint for future stellar modeling in global analyses of exoplanet systems. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 27 pages, 18 figures

Probing the Scalar WIMP-Pion Coupling with the first LUX-ZEPLIN data

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, E. E. Barillier, J. W. Bargemann, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. J. Bishop, G. M. Blockinger, B. Boxer , et al. (178 additional authors not shown)

Abstract: Weakly interacting massive particles (WIMPs) may interact with a virtual pion that is exchanged between nucleons. This interaction channel is important to consider in models where the spin-independent isoscalar channel is suppressed. Using data from the first science run of the LUX-ZEPLIN dark matter experiment, containing 60 live days of data in a 5.5~tonne fiducial mass of liquid xenon, we repor… ▽ More Weakly interacting massive particles (WIMPs) may interact with a virtual pion that is exchanged between nucleons. This interaction channel is important to consider in models where the spin-independent isoscalar channel is suppressed. Using data from the first science run of the LUX-ZEPLIN dark matter experiment, containing 60 live days of data in a 5.5~tonne fiducial mass of liquid xenon, we report the results on a search for WIMP-pion interactions. We observe no significant excess and set an upper limit of $1.5\times10^{-46}$~cm$^2$ at a 90\% confidence level for a WIMP mass of 33~GeV/c$^2$ for this interaction. △ Less

Submitted 4 June, 2024; originally announced June 2024.

The Widening Gap: The Benefits and Harms of Generative AI for Novice Programmers

Authors: James Prather, Brent Reeves, Juho Leinonen, Stephen MacNeil, Arisoa S. Randrianasolo, Brett Becker, Bailey Kimmel, Jared Wright, Ben Briggs

Abstract: Novice programmers often struggle through programming problem solving due to a lack of metacognitive awareness and strategies. Previous research has shown that novices can encounter multiple metacognitive difficulties while programming. Novices are typically unaware of how these difficulties are hindering their progress. Meanwhile, many novices are now programming with generative AI (GenAI), which… ▽ More Novice programmers often struggle through programming problem solving due to a lack of metacognitive awareness and strategies. Previous research has shown that novices can encounter multiple metacognitive difficulties while programming. Novices are typically unaware of how these difficulties are hindering their progress. Meanwhile, many novices are now programming with generative AI (GenAI), which can provide complete solutions to most introductory programming problems, code suggestions, hints for next steps when stuck, and explain cryptic error messages. Its impact on novice metacognition has only started to be explored. Here we replicate a previous study that examined novice programming problem solving behavior and extend it by incorporating GenAI tools. Through 21 lab sessions consisting of participant observation, interview, and eye tracking, we explore how novices are coding with GenAI tools. Although 20 of 21 students completed the assigned programming problem, our findings show an unfortunate divide in the use of GenAI tools between students who accelerated and students who struggled. Students who accelerated were able to use GenAI to create code they already intended to make and were able to ignore unhelpful or incorrect inline code suggestions. But for students who struggled, our findings indicate that previously known metacognitive difficulties persist, and that GenAI unfortunately can compound them and even introduce new metacognitive difficulties. Furthermore, struggling students often expressed cognitive dissonance about their problem solving ability, thought they performed better than they did, and finished with an illusion of competence. Based on our observations from both groups, we propose ways to scaffold the novice GenAI experience and make suggestions for future work. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: Accepted to ICER 2024

The Data Acquisition System of the LZ Dark Matter Detector: FADR

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, E. E. Barillier, J. W. Bargemann, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. Bishop, G. M. Blockinger, B. Boxer , et al. (191 additional authors not shown)

Abstract: The Data Acquisition System (DAQ) for the LUX-ZEPLIN (LZ) dark matter detector is described. The signals from 745 PMTs, distributed across three subsystems, are sampled with 100-MHz 32-channel digitizers (DDC-32s). A basic waveform analysis is carried out on the on-board Field Programmable Gate Arrays (FPGAs) to extract information about the observed scintillation and electroluminescence signals.… ▽ More The Data Acquisition System (DAQ) for the LUX-ZEPLIN (LZ) dark matter detector is described. The signals from 745 PMTs, distributed across three subsystems, are sampled with 100-MHz 32-channel digitizers (DDC-32s). A basic waveform analysis is carried out on the on-board Field Programmable Gate Arrays (FPGAs) to extract information about the observed scintillation and electroluminescence signals. This information is used to determine if the digitized waveforms should be preserved for offline analysis. The system is designed around the Kintex-7 FPGA. In addition to digitizing the PMT signals and providing basic event selection in real time, the flexibility provided by the use of FPGAs allows us to monitor the performance of the detector and the DAQ in parallel to normal data acquisition. The hardware and software/firmware of this FPGA-based Architecture for Data acquisition and Realtime monitoring (FADR) are discussed and performance measurements are described. △ Less

Submitted 16 August, 2024; v1 submitted 23 May, 2024; originally announced May 2024.

Comments: 18 pages, 24 figures

Microstructure and Stress Mapping in 3D at Industrially Relevant Degrees of Plastic Deformation

Authors: Axel Henningsson, Mustafacan Kutsal, Jonathan P. Wright, Wolfgang Ludwig, Henning Osholm Sørensen, Stephen A. Hall, Grethe Winther, Henning F. Poulsen

Abstract: Strength, ductility, and failure properties of metals are tailored by plastic deformation routes. Predicting these properties requires modeling of the structural dynamics and stress evolution taking place on several length scales. Progress has been hampered by a lack of representative 3D experimental data at industrially relevant degrees of deformation. We present an X-ray imaging based 3D mapping… ▽ More Strength, ductility, and failure properties of metals are tailored by plastic deformation routes. Predicting these properties requires modeling of the structural dynamics and stress evolution taking place on several length scales. Progress has been hampered by a lack of representative 3D experimental data at industrially relevant degrees of deformation. We present an X-ray imaging based 3D mapping of an aluminum polycrystal deformed to the ultimate tensile strength (32% elongation). The extensive dataset reveals significant intra-grain stress variations (36 MPa) up to at least half of the inter-grain variations (76 MPa), which are dominated by grain orientation effects. Local intra-grain stress concentrations are candidates for damage nucleation. Such data are important for models of structure-property relations and damage. △ Less

Submitted 21 May, 2024; v1 submitted 19 May, 2024; originally announced May 2024.

GRASS II: Simulations of Potential Granulation Noise Mitigation Methods

Authors: Michael L. Palumbo III, Eric B. Ford, Elizabeth B. Gonzalez, Jason T. Wright, Khaled Al Moulla, Rolf Schlichenmaier

Abstract: We present an updated version of GRASS (the GRanulation And Spectrum Simulator, Palumbo et al. 2022) which now uses an expanded library of 22 solar lines to empirically model time-resolved spectral variations arising from solar granulation. We show that our synthesis model accurately reproduces disk-integrated solar line profiles and bisectors, and we quantify the intrinsic granulation-driven radi… ▽ More We present an updated version of GRASS (the GRanulation And Spectrum Simulator, Palumbo et al. 2022) which now uses an expanded library of 22 solar lines to empirically model time-resolved spectral variations arising from solar granulation. We show that our synthesis model accurately reproduces disk-integrated solar line profiles and bisectors, and we quantify the intrinsic granulation-driven radial-velocity (RV) variability for each of the 22 lines studied. We show that summary statistics of bisector shape (e.g., bisector inverse slope) are strongly correlated with the measured anomalous, variability-driven RV at high pixel signal-to-noise ratio (SNR) and spectral resolution. Further, the strength of the correlations vary both line by line and with the summary statistic used. These correlations disappear for individual lines at the typical spectral resolutions and SNRs achieved by current EPRV spectrographs; so we use simulations from GRASS to demonstrate that they can, in principle, be recovered by selectively binning lines that are similarly affected by granulation. In the best-case scenario (high SNR and large number of binned lines), we find that a $\lesssim$30$\%$ reduction in the granulation-induced root mean square (RMS) RV can be achieved, but that the achievable reduction in variability is most strongly limited by the spectral resolution of the observing instrument. Based on our simulations, we predict that existing ultra-high-resolution spectrographs, namely ESPRESSO and PEPSI, should be able to resolve convective variability in other, bright stars. △ Less

Submitted 14 May, 2024; v1 submitted 13 May, 2024; originally announced May 2024.

Comments: 20 + 16 pages, 6 + 11 figures, 4 tables; in revision for AJ

Journal ref: AJ 168, 46 (2024)

Kinematic substructure in star clusters constrains star cluster formation

Authors: Becky Arnold, Nicholas J. Wright

Abstract: The spatial-kinematic structure of 48 young star clusters and associations is investigated. Moran's $I$ statistic is used to quantify the degree of kinematic substructure in each region, and the results are compared to those expected assuming the hierarchical or monolithic models of star cluster formation. Of the observed regions, 39 are found to have significant kinematic substructure, such that… ▽ More The spatial-kinematic structure of 48 young star clusters and associations is investigated. Moran's $I$ statistic is used to quantify the degree of kinematic substructure in each region, and the results are compared to those expected assuming the hierarchical or monolithic models of star cluster formation. Of the observed regions, 39 are found to have significant kinematic substructure, such that they are compatible with the hierarchical model and incompatible with the monolithic model. This includes multiple regions whose $Q$ parameter shows the region to be centrally concentrated and clustered. The remaining nine are compatible with both models. From this it is concluded that the kinematic substructure of the observed star clusters represents strong evidence in favour the hierarchical model of star cluster formation over the monolithic model. △ Less

Submitted 8 May, 2024; originally announced May 2024.

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

Logical Error Rates for a [[4,2,2]]-Encoded Variational Quantum Eigensolver Ansatz

Authors: Meenambika Gowrishankar, Daniel Claudino, Jerimiah Wright, Travis Humble

Abstract: Application benchmarks that run on noisy, intermediate-scale quantum (NISQ) computing devices require techniques for mitigating errors to improve accuracy and precision. Quantum error detection codes offer a framework by which to encode quantum computations and identify when errors occur. However, the subsequent logical error rate depends on the encoded application circuit as well as the underlyin… ▽ More Application benchmarks that run on noisy, intermediate-scale quantum (NISQ) computing devices require techniques for mitigating errors to improve accuracy and precision. Quantum error detection codes offer a framework by which to encode quantum computations and identify when errors occur. However, the subsequent logical error rate depends on the encoded application circuit as well as the underlying noise. Here, we quantify how the [[4,2,2]] quantum error detection code improves the logical error rate, accuracy, and precision of an encoded variational quantum eigensolver (VQE) application. We benchmark the performance of the encoded VQE for estimating the energy of the hydrogen molecule with a chemical accuracy of 1.6 mHa while managing the trade-off between probability of success of various post-selection methods. Using numerical simulation of the noisy mixed state preparation, we find that the most aggressive post-selection strategies improve the accuracy and precision of the encoded estimates even at the cost of increasing loss of samples. △ Less

Submitted 5 May, 2024; originally announced May 2024.

Project Hephaistos - II. Dyson sphere candidates from Gaia DR3, 2MASS, and WISE

Authors: Matías Suazo, Erik Zackrisson, Priyatam K. Mahto, Fabian Lundell, Carl Nettelblad, Andreas J. Korn, Jason T. Wright, Suman Majumdar

Abstract: The search for extraterrestrial intelligence is currently being pursued using multiple techniques and in different wavelength bands. Dyson spheres, megastructures that could be constructed by advanced civilizations to harness the radiation energy of their host stars, represent a potential technosignature, that in principle may be hiding in public data already collected as part of large astronomica… ▽ More The search for extraterrestrial intelligence is currently being pursued using multiple techniques and in different wavelength bands. Dyson spheres, megastructures that could be constructed by advanced civilizations to harness the radiation energy of their host stars, represent a potential technosignature, that in principle may be hiding in public data already collected as part of large astronomical surveys. In this study, we present a comprehensive search for partial Dyson spheres by analyzing optical and infrared observations from Gaia, 2MASS, and WISE. We develop a pipeline that employs multiple filters to identify potential candidates and reject interlopers in a sample of five million objects, which incorporates a convolutional neural network to help identify confusion in WISE data. Finally, the pipeline identifies 7 candidates deserving of further analysis. All of these objects are M-dwarfs, for which astrophysical phenomena cannot easily account for the observed infrared excess emission. △ Less

Submitted 5 May, 2024; originally announced May 2024.

Comments: Accepted to be published in MNRAS

Constraints On Covariant WIMP-Nucleon Effective Field Theory Interactions from the First Science Run of the LUX-ZEPLIN Experiment

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, E. E. Barillier, J. W. Bargemann, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. J. Bishop, G. M. Blockinger, B. Boxer , et al. (179 additional authors not shown)

Abstract: The first science run of the LUX-ZEPLIN (LZ) experiment, a dual-phase xenon time project chamber operating in the Sanford Underground Research Facility in South Dakota, USA, has reported leading limits on spin-independent WIMP-nucleon interactions and interactions described from a non-relativistic effective field theory (NREFT). Using the same 5.5~t fiducial mass and 60 live days of exposure we re… ▽ More The first science run of the LUX-ZEPLIN (LZ) experiment, a dual-phase xenon time project chamber operating in the Sanford Underground Research Facility in South Dakota, USA, has reported leading limits on spin-independent WIMP-nucleon interactions and interactions described from a non-relativistic effective field theory (NREFT). Using the same 5.5~t fiducial mass and 60 live days of exposure we report on the results of a relativistic extension to the NREFT. We present constraints on couplings from covariant interactions arising from the coupling of vector, axial currents, and electric dipole moments of the nucleon to the magnetic and electric dipole moments of the WIMP which cannot be described by recasting previous results described by an NREFT. Using a profile-likelihood ratio analysis, in an energy region between 0~keV$_\text{nr}$ to 270~keV$_\text{nr}$, we report 90% confidence level exclusion limits on the coupling strength of five interactions in both the isoscalar and isovector bases. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: 7 pages, 4 figures

Disentangling the co-evolution of galaxies and supermassive black holes with PRIMA

Authors: L. Bisigello, C. Gruppioni, A. Bolatto, L. Ciesla, A. Pope, L. Armus, L., J. D. Smith, R. Somerville, L. Y. A. Yung, R. J. Wright, C. M. Bradford, J. Glenn, A. Feltre

Abstract: The most active phases of star formation and black hole accretion are strongly affected by dust extinction, making far-infrared (far-IR) observations the best way to disentangle and study the co-evolution of galaxies and super massive black holes. The plethora of fine structure lines and emission features from dust, ionised and neutral atomic and warm molecular gas in the rest-frame mid- and far-I… ▽ More The most active phases of star formation and black hole accretion are strongly affected by dust extinction, making far-infrared (far-IR) observations the best way to disentangle and study the co-evolution of galaxies and super massive black holes. The plethora of fine structure lines and emission features from dust, ionised and neutral atomic and warm molecular gas in the rest-frame mid- and far-IR provide unmatched diagnostic power to determine the properties of gas and dust, measure gas-phase metallicities and map cold galactic outflows in even the most obscured galaxies. By combining multi-band photometric surveys with low and high-resolution far-IR spectroscopy, the PRobe far-Infrared Mission for Astrophysics (PRIMA), a concept for a far-IR, 1.8m-diameter, cryogenically cooled observatory, will revolutionise the field of galaxy evolution by taking advantage of this IR toolkit to find and study dusty galaxies across galactic time. In this work, we make use of the phenomenological simulation SPRITZ and the Santa Cruz semi-analytical model to describe how a moderately deep multi-band PRIMA photometric survey can easily reach beyond previous IR missions to detect and study galaxies down to $10^{11}\,L_{\odot}$ beyond cosmic noon and at least up to z=4, even in the absence of gravitational lensing. By decomposing the spectral energy distribution (SED) of these photometrically selected galaxies, we show that PRIMA can be used to accurately measure the relative AGN power, the mass fraction contributed by polycyclic aromatic hydrocarbon (PAH) and the total IR luminosity. At the same time, spectroscopic follow up with PRIMA will allow to trace both the star formation and black hole accretion rates (SFR, BHAR), the gas phase metallicities and the mass outflow rates of cold gas in hundreds to thousands of individual galaxies to z=2. △ Less

Submitted 3 June, 2024; v1 submitted 26 April, 2024; originally announced April 2024.

Comments: 16 pages, 17 figures, accepted for publication in A&A

The death of Vulcan: NEID reveals the planet candidate orbiting HD 26965 is stellar activity

Authors: Abigail Burrows, Samuel Halverson, Jared C. Siegel, Christian Gilbertson, Jacob Luhn, Jennifer Burt, Chad F. Bender, Arpita Roy, Ryan C. Terrien, Selma Vangstein, Suvrath Mahadevan, Jason T. Wright, Paul Robertson, Eric B. Ford, Guðmundur Stefánsson, Joe P. Ninan, Cullen H. Blake, Michael W. McElwain, Christian Schwab, Jinglin Zhao

Abstract: We revisit the long-studied radial velocity (RV) target HD26965 using recent observations from the NASA-NSF 'NEID' precision Doppler facility. Leveraging a suite of classical activity indicators, combined with line-by-line RV analyses, we demonstrate that the claimed 45-day signal previously identified as a planet candidate is most likely an activity-induced signal. Correlating the bulk (spectrall… ▽ More We revisit the long-studied radial velocity (RV) target HD26965 using recent observations from the NASA-NSF 'NEID' precision Doppler facility. Leveraging a suite of classical activity indicators, combined with line-by-line RV analyses, we demonstrate that the claimed 45-day signal previously identified as a planet candidate is most likely an activity-induced signal. Correlating the bulk (spectrally-averaged) RV with canonical line activity indicators confirms a multi-day 'lag' between the observed activity indicator time series and the measured RV. When accounting for this lag, we show that much of the observed RV signal can be removed by a linear detrending of the data. Investigating activity at the line-by-line level, we find a depth-dependent correlation between individual line RVs and the bulk RVs, further indicative of periodic suppression of convective blueshift causing the observed RV variability, rather than an orbiting planet. We conclude that the combined evidence of the activity correlations and depth dependence is consistent with a radial velocity signature dominated by a rotationally-modulated activity signal at a period of $\sim$42 days. We hypothesize that this activity signature is due to a combination of spots and convective blueshift suppression. The tools applied in our analysis are broadly applicable to other stars, and could help paint a more comprehensive picture of the manifestations of stellar activity in future Doppler RV surveys. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: 25 pages, 13 figures. Accepted in AJ

Journal ref: AJ 167 243 (2024)

Direct Zernike Coefficient Prediction from Point Spread Functions and Extended Images using Deep Learning

Authors: Yong En Kok, Alexander Bentley, Andrew Parkes, Amanda J. Wright, Michael G. Somekh, Michael Pound

Abstract: Optical imaging quality can be severely degraded by system and sample induced aberrations. Existing adaptive optics systems typically rely on iterative search algorithm to correct for aberrations and improve images. This study demonstrates the application of convolutional neural networks to characterise the optical aberration by directly predicting the Zernike coefficients from two to three phase-… ▽ More Optical imaging quality can be severely degraded by system and sample induced aberrations. Existing adaptive optics systems typically rely on iterative search algorithm to correct for aberrations and improve images. This study demonstrates the application of convolutional neural networks to characterise the optical aberration by directly predicting the Zernike coefficients from two to three phase-diverse optical images. We evaluated our network on 600,000 simulated Point Spread Function (PSF) datasets randomly generated within the range of -1 to 1 radians using the first 25 Zernike coefficients. The results show that using only three phase-diverse images captured above, below and at the focal plane with an amplitude of 1 achieves a low RMSE of 0.10 radians on the simulated PSF dataset. Furthermore, this approach directly predicts Zernike modes simulated extended 2D samples, while maintaining a comparable RMSE of 0.15 radians. We demonstrate that this approach is effective using only a single prediction step, or can be iterated a small number of times. This simple and straightforward technique provides rapid and accurate method for predicting the aberration correction using three or less phase-diverse images, paving the way for evaluation on real-world dataset. △ Less

Submitted 24 April, 2024; v1 submitted 23 April, 2024; originally announced April 2024.

Comments: 12 pages, 6 figures, 4 tables

Observation of Gravitational Waves from the Coalescence of a $2.5\text{-}4.5~M_\odot$ Compact Object and a Neutron Star

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, I. Abouelfettouh, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, D. Agarwal, M. Agathos, M. Aghaei Abchouyeh, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, S. Akçay, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Al-Jodah , et al. (1771 additional authors not shown)

Abstract: We report the observation of a coalescing compact binary with component masses $2.5\text{-}4.5~M_\odot$ and $1.2\text{-}2.0~M_\odot$ (all measurements quoted at the 90% credible level). The gravitational-wave signal GW230529_181500 was observed during the fourth observing run of the LIGO-Virgo-KAGRA detector network on 2023 May 29 by the LIGO Livingston Observatory. The primary component of the so… ▽ More We report the observation of a coalescing compact binary with component masses $2.5\text{-}4.5~M_\odot$ and $1.2\text{-}2.0~M_\odot$ (all measurements quoted at the 90% credible level). The gravitational-wave signal GW230529_181500 was observed during the fourth observing run of the LIGO-Virgo-KAGRA detector network on 2023 May 29 by the LIGO Livingston Observatory. The primary component of the source has a mass less than $5~M_\odot$ at 99% credibility. We cannot definitively determine from gravitational-wave data alone whether either component of the source is a neutron star or a black hole. However, given existing estimates of the maximum neutron star mass, we find the most probable interpretation of the source to be the coalescence of a neutron star with a black hole that has a mass between the most massive neutron stars and the least massive black holes observed in the Galaxy. We provisionally estimate a merger rate density of $55^{+127}_{-47}~\text{Gpc}^{-3}\,\text{yr}^{-1}$ for compact binary coalescences with properties similar to the source of GW230529_181500; assuming that the source is a neutron star-black hole merger, GW230529_181500-like sources constitute about 60% of the total merger rate inferred for neutron star-black hole coalescences. The discovery of this system implies an increase in the expected rate of neutron star-black hole mergers with electromagnetic counterparts and provides further evidence for compact objects existing within the purported lower mass gap. △ Less

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

Comments: 45 pages (10 pages author list, 13 pages main text, 1 page acknowledgements, 13 pages appendices, 8 pages bibliography), 17 figures, 16 tables. Update to match version published in The Astrophysical Journal Letters. Data products available from https://zenodo.org/records/10845779

Report number: LIGO-P2300352

Journal ref: ApJL 970, L34 (2024)

Robust Second-Order Nonconvex Optimization and Its Application to Low Rank Matrix Sensing

Authors: Shuyao Li, Yu Cheng, Ilias Diakonikolas, Jelena Diakonikolas, Rong Ge, Stephen J. Wright

Abstract: Finding an approximate second-order stationary point (SOSP) is a well-studied and fundamental problem in stochastic nonconvex optimization with many applications in machine learning. However, this problem is poorly understood in the presence of outliers, limiting the use of existing nonconvex algorithms in adversarial settings. In this paper, we study the problem of finding SOSPs in the strong c… ▽ More Finding an approximate second-order stationary point (SOSP) is a well-studied and fundamental problem in stochastic nonconvex optimization with many applications in machine learning. However, this problem is poorly understood in the presence of outliers, limiting the use of existing nonconvex algorithms in adversarial settings. In this paper, we study the problem of finding SOSPs in the strong contamination model, where a constant fraction of datapoints are arbitrarily corrupted. We introduce a general framework for efficiently finding an approximate SOSP with \emph{dimension-independent} accuracy guarantees, using $\widetilde{O}({D^2}/ε)$ samples where $D$ is the ambient dimension and $ε$ is the fraction of corrupted datapoints. As a concrete application of our framework, we apply it to the problem of low rank matrix sensing, developing efficient and provably robust algorithms that can tolerate corruptions in both the sensing matrices and the measurements. In addition, we establish a Statistical Query lower bound providing evidence that the quadratic dependence on $D$ in the sample complexity is necessary for computationally efficient algorithms. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Does the Fundamental Metallicity Relation Evolve with Redshift? I: The Correlation Between Offsets from the Mass-Metallicity Relation and Star Formation Rate

Authors: Alex M. Garcia, Paul Torrey, Sara Ellison, Kathryn Grasha, Lars Hernquist, Henry R. M. Zovaro, Qian-Hui Chen, Z. S. Hemler, Lisa J. Kewley, Erica J. Nelson, Ruby J. Wright

Abstract: The scatter about the mass-metallicity relation (MZR) has a correlation with the star formation rate (SFR) of galaxies. The lack of evidence of evolution in correlated scatter at $z\lesssim2.5$ leads many to refer to the relationship between mass, metallicity, and SFR as the Fundamental Metallicity Relation (FMR). Yet, recent high-redshift (z>3) JWST observations have challenged the fundamental (i… ▽ More The scatter about the mass-metallicity relation (MZR) has a correlation with the star formation rate (SFR) of galaxies. The lack of evidence of evolution in correlated scatter at $z\lesssim2.5$ leads many to refer to the relationship between mass, metallicity, and SFR as the Fundamental Metallicity Relation (FMR). Yet, recent high-redshift (z>3) JWST observations have challenged the fundamental (i.e., redshift-invariant) nature of the FMR. In this work, we show that the cosmological simulations Illustris, IllustrisTNG, and EAGLE all predict MZRs that exhibit scatter with a secondary dependence on SFR up to $z=8$. We introduce the concept of a "strong" FMR, where the strength of correlated scatter does not evolve with time, and a "weak" FMR, where there is some time evolution. We find that each simulation analysed has a weak FMR -- there is non-negligible evolution in the strength of the correlation with SFR. Furthermore, we show that the scatter is reduced an additional ~10-40% at $z\gtrsim3$ when using a weak FMR, compared to assuming a strong FMR. These results highlight the importance of avoiding coarse redshift binning when assessing the FMR. △ Less

Submitted 10 May, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Comments: 10 pages, 5 figures (Figure 5 is summary). Accepted to MNRAS

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

How to Make the Gradients Small Privately: Improved Rates for Differentially Private Non-Convex Optimization

Authors: Andrew Lowy, Jonathan Ullman, Stephen J. Wright

Abstract: We provide a simple and flexible framework for designing differentially private algorithms to find approximate stationary points of non-convex loss functions. Our framework is based on using a private approximate risk minimizer to "warm start" another private algorithm for finding stationary points. We use this framework to obtain improved, and sometimes optimal, rates for several classes of non-c… ▽ More We provide a simple and flexible framework for designing differentially private algorithms to find approximate stationary points of non-convex loss functions. Our framework is based on using a private approximate risk minimizer to "warm start" another private algorithm for finding stationary points. We use this framework to obtain improved, and sometimes optimal, rates for several classes of non-convex loss functions. First, we obtain improved rates for finding stationary points of smooth non-convex empirical loss functions. Second, we specialize to quasar-convex functions, which generalize star-convex functions and arise in learning dynamical systems and training some neural nets. We achieve the optimal rate for this class. Third, we give an optimal algorithm for finding stationary points of functions satisfying the Kurdyka-Lojasiewicz (KL) condition. For example, over-parameterized neural networks often satisfy this condition. Fourth, we provide new state-of-the-art rates for stationary points of non-convex population loss functions. Fifth, we obtain improved rates for non-convex generalized linear models. A modification of our algorithm achieves nearly the same rates for second-order stationary points of functions with Lipschitz Hessian, improving over the previous state-of-the-art for each of the above problems. △ Less

Submitted 19 August, 2024; v1 submitted 16 February, 2024; originally announced February 2024.

Comments: ICML 2024

New constraints on ultraheavy dark matter from the LZ experiment

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, A. Baker, S. Balashov, J. Bang, J. W. Bargemann, A. Baxter, K. Beattie, T. Benson, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, E. Bishop, G. M. Blockinger, B. Boxer, C. A. J. Brew , et al. (174 additional authors not shown)

Abstract: Searches for dark matter with liquid xenon time projection chamber experiments have traditionally focused on the region of the parameter space that is characteristic of weakly interacting massive particles, ranging from a few GeV/$c^2$ to a few TeV/$c^2$. Models of dark matter with a mass much heavier than this are well motivated by early production mechanisms different from the standard thermal f… ▽ More Searches for dark matter with liquid xenon time projection chamber experiments have traditionally focused on the region of the parameter space that is characteristic of weakly interacting massive particles, ranging from a few GeV/$c^2$ to a few TeV/$c^2$. Models of dark matter with a mass much heavier than this are well motivated by early production mechanisms different from the standard thermal freeze-out, but they have generally been less explored experimentally. In this work, we present a re-analysis of the first science run (SR1) of the LZ experiment, with an exposure of $0.9$ tonne$\times$year, to search for ultraheavy particle dark matter. The signal topology consists of multiple energy deposits in the active region of the detector forming a straight line, from which the velocity of the incoming particle can be reconstructed on an event-by-event basis. Zero events with this topology were observed after applying the data selection calibrated on a simulated sample of signal-like events. New experimental constraints are derived, which rule out previously unexplored regions of the dark matter parameter space of spin-independent interactions beyond a mass of 10$^{17}$ GeV/$c^2$. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 9 pages, 7 figures

Journal ref: Phys. Rev. D 109, 112010 (2024)

The baryon cycle in modern cosmological hydrodynamical simulations

Authors: Ruby J. Wright, Rachel S. Somerville, Claudia del P. Lagos, Matthieu Schaller, Romeel Davé, Daniel Anglés-Alcázar, Shy Genel

Abstract: In recent years, cosmological hydrodynamical simulations have proven their utility as key interpretative tools in the study of galaxy formation and evolution. In this work, we present a like-for-like comparison between the baryon cycle in three publicly available, leading cosmological simulation suites: EAGLE, IllustrisTNG, and SIMBA. While these simulations broadly agree in terms of their predict… ▽ More In recent years, cosmological hydrodynamical simulations have proven their utility as key interpretative tools in the study of galaxy formation and evolution. In this work, we present a like-for-like comparison between the baryon cycle in three publicly available, leading cosmological simulation suites: EAGLE, IllustrisTNG, and SIMBA. While these simulations broadly agree in terms of their predictions for the stellar mass content and star formation rates of galaxies at $z\approx0$, they achieve this result for markedly different reasons. In EAGLE and SIMBA, we demonstrate that at low halo masses ($M_{\rm 200c}\lesssim 10^{11.5}\, M_{\odot}$), stellar feedback (SF)-driven outflows can reach far beyond the scale of the halo, extending up to $2-3\times R_{\rm 200c}$. In contrast, in TNG, SF-driven outflows, while stronger at the scale of the ISM, recycle within the CGM (within $R_{\rm 200c}$). We find that AGN-driven outflows in SIMBA are notably potent, reaching several times $R_{\rm 200c}$ even at halo masses up to $M_{\rm 200c}\approx10^{13.5}\, M_{\odot}$. In both TNG and EAGLE, AGN feedback can eject gas beyond $R_{\rm 200c}$ at this mass scale, but seldom beyond $2-3\times R_{\rm 200c}$. We find that the scale of feedback-driven outflows can be directly linked with the prevention of cosmological inflow, as well as the total baryon fraction of haloes within $R_{\rm 200c}$. This work lays the foundation to develop targeted observational tests that can discriminate between feedback scenarios, and inform sub-grid feedback models in the next generation of simulations. △ Less

Submitted 9 July, 2024; v1 submitted 13 February, 2024; originally announced February 2024.

Comments: 20 pages, 8 figures. Accepted for publication in MNRAS 03/07/24

Grain boundary strain localization in CdTe solar cell revealed by Scanning 3D X-ray diffraction microscopy

Authors: A. Shukla, J. Wright, A. Henningsson, H. Stieglitz, E. Colegrove, L. Besley, C. Baur, S De Angelis, M. Stuckelberger, H. F. Poulsen, J. W. Andreasen

Abstract: Cadmium Telluride (CdTe) solar cell technology is a promising candidate to help boost green energy production. However, impurities and structural defects are major barriers to improving the solar power conversion efficiency. Grain boundaries often act as aggregation sites for impurities, resulting in strain localization in areas of high diffusion. In this study, we demonstrate the use of scanning… ▽ More Cadmium Telluride (CdTe) solar cell technology is a promising candidate to help boost green energy production. However, impurities and structural defects are major barriers to improving the solar power conversion efficiency. Grain boundaries often act as aggregation sites for impurities, resulting in strain localization in areas of high diffusion. In this study, we demonstrate the use of scanning 3D X-ray diffraction microscopy to non-destructively make 3D maps of the grains, their phase, orientation, and local strain within a CdTe solar cell absorber layer with a resolution of 100 nm. We quantify twin boundaries and suggest how they affect grain size and orientation distribution. Local strain analysis reveals that strain is primarily associated with high misorientation grain boundaries, whereas twin boundaries do not show high strain values. We also observe that high-strain grain boundaries form a continuous pathway connected to the CdS layer. Hence, this high-strain region is believed to be associated with the diffusion of sulfur from the CdS layer along grain boundaries. This hypothesis is supported by SEM EDS and X-ray fluorescence experiments. The method and analysis demonstrated in this work can be applied to different polycrystalline materials where the characterization of grain boundary properties is essential to understand the microstructural phenomena. △ Less

Submitted 11 February, 2024; originally announced February 2024.

Revisiting Inexact Fixed-Point Iterations for Min-Max Problems: Stochasticity and Structured Nonconvexity

Authors: Ahmet Alacaoglu, Donghwan Kim, Stephen J. Wright

Abstract: We focus on constrained, $L$-smooth, potentially stochastic and nonconvex-nonconcave min-max problems either satisfying $ρ$-cohypomonotonicity or admitting a solution to the $ρ$-weakly Minty Variational Inequality (MVI), where larger values of the parameter $ρ>0$ correspond to a greater degree of nonconvexity. These problem classes include examples in two player reinforcement learning, interaction… ▽ More We focus on constrained, $L$-smooth, potentially stochastic and nonconvex-nonconcave min-max problems either satisfying $ρ$-cohypomonotonicity or admitting a solution to the $ρ$-weakly Minty Variational Inequality (MVI), where larger values of the parameter $ρ>0$ correspond to a greater degree of nonconvexity. These problem classes include examples in two player reinforcement learning, interaction dominant min-max problems, and certain synthetic test problems on which classical min-max algorithms fail. It has been conjectured that first-order methods can tolerate a value of $ρ$ no larger than $\frac{1}{L}$, but existing results in the literature have stagnated at the tighter requirement $ρ< \frac{1}{2L}$. With a simple argument, we obtain optimal or best-known complexity guarantees with cohypomonotonicity or weak MVI conditions for $ρ< \frac{1}{L}$. First main insight for the improvements in the convergence analyses is to harness the recently proposed $\textit{conic nonexpansiveness}$ property of operators. Second, we provide a refined analysis for inexact Halpern iteration that relaxes the required inexactness level to improve some state-of-the-art complexity results even for constrained stochastic convex-concave min-max problems. Third, we analyze a stochastic inexact Krasnosel'skiĭ-Mann iteration with a multilevel Monte Carlo estimator when the assumptions only hold with respect to a solution. △ Less

Submitted 12 August, 2024; v1 submitted 7 February, 2024; originally announced February 2024.

Journal ref: Proceedings of the International Conference on Machine Learning (ICML) 2024

Small jet engine reservoir computing digital twin

Authors: C. J. Wright, N. Biederman, B. Gyovai, D. J. Gauthier, J. P. Wilhelm

Abstract: Machine learning was applied to create a digital twin of a numerical simulation of a single-scroll jet engine. A similar model based on the insights gained from this numerical study was used to create a digital twin of a JetCat P100-RX jet engine using only experimental data. Engine data was collected from a custom sensor system measuring parameters such as thrust, exhaust gas temperature, shaft s… ▽ More Machine learning was applied to create a digital twin of a numerical simulation of a single-scroll jet engine. A similar model based on the insights gained from this numerical study was used to create a digital twin of a JetCat P100-RX jet engine using only experimental data. Engine data was collected from a custom sensor system measuring parameters such as thrust, exhaust gas temperature, shaft speed, weather conditions, etc. Data was gathered while the engine was placed under different test conditions by controlling shaft speed. The machine learning model was generated (trained) using a next-generation reservoir computer, a best-in-class machine learning algorithm for dynamical systems. Once the model was trained, it was used to predict behavior it had never seen with an accuracy of better than 1.8% when compared to the testing data. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Approximation of Calogero-Moser lattices by Benjamin-Ono equations

Authors: J. Douglas Wright

Abstract: We provide a rigorous validation that the infinite Calogero-Moser lattice can be well-approximated by solutions of the Benjamin-Ono equation in a long-wave limit. We provide a rigorous validation that the infinite Calogero-Moser lattice can be well-approximated by solutions of the Benjamin-Ono equation in a long-wave limit. △ Less

Submitted 15 December, 2023; originally announced December 2023.

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 largely linear response of earth's ice volume to orbital forcing

Authors: Liam Wheen, Oscar Benjamin, Thomas Gernon, Cameron Hall, Jerry Wright

Abstract: Orbital forcing plays a key role in pacing the glacial-interglacial cycles. However, the mechanistic linkages between the orbital parameters - eccentricity, obliquity, and precession - and global ice volume remain unclear. Here, we investigate the effect of Earth's orbitally governed incoming solar radiation (that is, insolation) on global ice volume over the past 800,000 years. We consider a simp… ▽ More Orbital forcing plays a key role in pacing the glacial-interglacial cycles. However, the mechanistic linkages between the orbital parameters - eccentricity, obliquity, and precession - and global ice volume remain unclear. Here, we investigate the effect of Earth's orbitally governed incoming solar radiation (that is, insolation) on global ice volume over the past 800,000 years. We consider a simple linear model of ice volume that imposes minimal assumptions about its dynamics. We find that this model can adequately reproduce the observed ice volume variations for most of the past 800,000 years, with the notable exception of Marine Isotope Stage 11. This suggests that, aside from a few extrema, the ice volume dynamics primarily result from an approximately linear response to orbital forcing. We substantiate this finding by addressing some of the key criticisms of the orbitally forced hypothesis. In particular, we show that eccentricity can significantly vary the ocean temperature without the need for amplification on Earth. We also present a feasible mechanism to explain the absence of eccentricity's 400,000 year period in the ice volume data. This requires part of the forcing from eccentricity to be lagged via a slow-responding mechanism, resulting in a signal that closer approximates the change in eccentricity. A physical interpretation of our model is proposed, using bulk ocean and surface temperatures as intermediate mechanisms through which the orbital parameters affect ice volume. These show reasonable alignment with their relevant proxy data, though we acknowledge that these variables likely represent a combination of mechanisms. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: 25 Pages, 12 figures, reproducing code can be found at: https://github.com/liamwheen/Wheen_23