Showing 1–50 of 779 results for author: Kaplan, D

Volume-optimal persistence homological scaffolds of hemodynamic networks covary with MEG theta-alpha aperiodic dynamics

Authors: Nghi Nguyen, Tao Hou, Enrico Amico, Jingyi Zheng, Huajun Huang, Alan D. Kaplan, Giovanni Petri, Joaquín Goñi, Yize Zhao, Duy Duong-Tran, Li Shen

Abstract: Higher-order properties of functional magnetic resonance imaging (fMRI) induced connectivity have been shown to unravel many exclusive topological and dynamical insights beyond pairwise interactions. Nonetheless, whether these fMRI-induced higher-order properties play a role in disentangling other neuroimaging modalities' insights remains largely unexplored and poorly understood. In this work, by… ▽ More Higher-order properties of functional magnetic resonance imaging (fMRI) induced connectivity have been shown to unravel many exclusive topological and dynamical insights beyond pairwise interactions. Nonetheless, whether these fMRI-induced higher-order properties play a role in disentangling other neuroimaging modalities' insights remains largely unexplored and poorly understood. In this work, by analyzing fMRI data from the Human Connectome Project Young Adult dataset using persistent homology, we discovered that the volume-optimal persistence homological scaffolds of fMRI-based functional connectomes exhibited conservative topological reconfigurations from the resting state to attentional task-positive state. Specifically, while reflecting the extent to which each cortical region contributed to functional cycles following different cognitive demands, these reconfigurations were constrained such that the spatial distribution of cavities in the connectome is relatively conserved. Most importantly, such level of contributions covaried with powers of aperiodic activities mostly within the theta-alpha (4-12 Hz) band measured by magnetoencephalography (MEG). This comprehensive result suggests that fMRI-induced hemodynamics and MEG theta-alpha aperiodic activities are governed by the same functional constraints specific to each cortical morpho-structure. Methodologically, our work paves the way toward an innovative computing paradigm in multimodal neuroimaging topological learning. △ Less

Submitted 6 July, 2024; originally announced July 2024.

Comments: The code for our analyses is provided in https://github.com/ngcaonghi/scaffold_noise

An Image-Based Search for Pulsar Candidates in the MeerKAT Bulge Survey

Authors: Dale A. Frail, Emil Polisensky, Scott D. Hyman, W. M. Cotton, Namir E. Kassim, Michele L. Silverstein, Rahul Sengar, David L. Kaplan, Francesca Calore, Joanna Berteaud, Maica Clavel, Marisa Geyer, Samuel Legodi, Vasaant Krishnan, Sarah Buchner, Fernando Camilo

Abstract: We report on the results of an image-based search for pulsar candidates toward the Galactic bulge. We used mosaic images from the MeerKAT radio telescope, that were taken as part of a 173 deg**2 survey of the bulge and Galactic center of our Galaxy at L band (856-1712 MHz) in all four Stokes I, Q, U and V. The image root-mean-square noise levels of 12-17 uJy/ba represent a significant increase in… ▽ More We report on the results of an image-based search for pulsar candidates toward the Galactic bulge. We used mosaic images from the MeerKAT radio telescope, that were taken as part of a 173 deg**2 survey of the bulge and Galactic center of our Galaxy at L band (856-1712 MHz) in all four Stokes I, Q, U and V. The image root-mean-square noise levels of 12-17 uJy/ba represent a significant increase in sensitivity over past image-based pulsar searches. Our primary search criterion was circular polarization, but we used other criteria including linear polarization, in-band spectral index, compactness, variability and multi-wavelength counterparts to select pulsar candidates. We first demonstrate the efficacy of this technique by searching for polarized emission from known pulsars, and comparing our results with measurements from the literature. Our search resulted in a sample of 75 polarized pulsar candidates. Bright stars or young stellar objects were associated with 28 of these sources, including a small sample of highly polarized dwarf stars with pulsar-like steep spectra. Comparing the properties of this sample with the known pulsars, we identified 30 compelling candidates for pulsation follow-up, including two sources with both strong circular and linear polarization. The remaining 17 sources are either pulsars or stars, but we cannot rule out an extragalactic origin or image artifacts among the brighter, flat spectrum objects. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 28 pages, 10 figures, 6 tables. ApJ, submitted

A principled framework to assess information theoretical fitness of brain functional sub-circuits

Authors: Duy Duong-Tran, Nghi Nguyen, Shizhuo Mu, Jiong Chen, Jingxuan Bao, Frederick Xu, Sumita Garai, Jose Cadena-Pico, Alan David Kaplan, Tianlong Chen, Yize Zhao, Li Shen, Joaquín Goñi

Abstract: In systems and network neuroscience, many common practices in brain connectomic analysis are often not properly scrutinized. One such practice is mapping a predetermined set of sub-circuits, like functional networks (FNs), onto subjects' functional connectomes (FCs) without adequately assessing the information-theoretic appropriateness of the partition. Another practice that goes unchallenged is t… ▽ More In systems and network neuroscience, many common practices in brain connectomic analysis are often not properly scrutinized. One such practice is mapping a predetermined set of sub-circuits, like functional networks (FNs), onto subjects' functional connectomes (FCs) without adequately assessing the information-theoretic appropriateness of the partition. Another practice that goes unchallenged is thresholding weighted FCs to remove spurious connections without justifying the chosen threshold. This paper leverages recent theoretical advances in Stochastic Block Models (SBMs) to formally define and quantify the information-theoretic fitness (e.g., prominence) of a predetermined set of FNs when mapped to individual FCs under different fMRI task conditions. Our framework allows for evaluating any combination of FC granularity, FN partition, and thresholding strategy, thereby optimizing these choices to preserve important topological features of the human brain connectomes. Our results pave the way for the proper use of predetermined FNs and thresholding methods and provide insights for future research in individualized parcellations. △ Less

Submitted 26 June, 2024; originally announced June 2024.

A Precision Gyroscope from the Spin of Light

Authors: Michael A. Fedderke, Roni Harnik, David E. Kaplan, Sam Posen, Surjeet Rajendran, Francesco Serra, Vyacheslav P. Yakovlev

Abstract: We describe a gyroscope that measures rotation based on the effects of the rotation on the polarization of light. Rotation induces a differential phase shift in the propagation of left and right circularly polarized light and this phase shift can be measured in suitably designed interferometric setups. The signal in this setup is independent of the frequency of light, unlike various sources of noi… ▽ More We describe a gyroscope that measures rotation based on the effects of the rotation on the polarization of light. Rotation induces a differential phase shift in the propagation of left and right circularly polarized light and this phase shift can be measured in suitably designed interferometric setups. The signal in this setup is independent of the frequency of light, unlike various sources of noise such as vibrations, which cause phase shifts that depend on the frequency. Such vibrations are the practical limit on the sensitivity of conventional Sagnac-style optical interferometers that are typically used as gyroscopes. In the proposed setup, one can potentially mitigate this source of noise by simultaneously using two (or more) sources of light that have different frequencies. The signal in this setup scales with the total storage time of the light. Due to its frequency independence, it is thus most optimal to measure the signal using superconducting RF systems where the high finesse of the available cavities enables considerably longer storage times than is possible in an optical setup. △ Less

Submitted 23 June, 2024; originally announced June 2024.

Comments: 15 pages, 1 figure

A two-minute burst of highly polarised radio emission originating from low Galactic latitude

Authors: Dougal Dobie, Andrew Zic, Lucy S. Oswald, Joshua Pritchard, Marcus E. Lower, Ziteng Wang, Hao Qiu, Natasha Hurley-Walker, Yuanming Wang, Emil Lenc, David L. Kaplan, Akash Anumarlapudi, Katie Auchettl, Matthew Bailes, Andrew D. Cameron, Jeffrey Cooke, Adam Deller, Laura N. Driessen, James Freeburn, Tara Murphy, Ryan M. Shannon, Adam J. Stewart

Abstract: Several sources of repeating coherent bursts of radio emission with periods of many minutes have now been reported in the literature. These ``ultra-long period'' (ULP) sources have no clear multi-wavelength counterparts and challenge canonical pulsar emission models, leading to debate regarding their nature. In this work we report the discovery of a bright, highly-polarised burst of radio emission… ▽ More Several sources of repeating coherent bursts of radio emission with periods of many minutes have now been reported in the literature. These ``ultra-long period'' (ULP) sources have no clear multi-wavelength counterparts and challenge canonical pulsar emission models, leading to debate regarding their nature. In this work we report the discovery of a bright, highly-polarised burst of radio emission at low Galactic latitude as part of a wide-field survey for transient and variable radio sources. ASKAP\,J175534.9$-$252749.1 does not appear to repeat, with only a single intense two-minute $\sim 200$\,mJy burst detected from 60~hours of observations. The burst morphology and polarisation properties are comparable to those of classical pulsars but the duration is more than one hundred times longer, analogous to ULPs. No comparable bursts are detected in the rest of our widefield survey to date. Combined with the existing ULP population, this suggests that these sources have a strong Galactic latitude dependence and hints at an unexplored population of transient and variable radio sources in the thin disk of the Milky Way. The resemblance of this burst with both ULPs and pulsars calls for a unified coherent emission model for objects with spin periods from milliseconds to tens of minutes. However, whether or not these are all neutron stars or have the same underlying power source remains open for debate. △ Less

Submitted 18 June, 2024; originally announced June 2024.

An Untargeted Search for Radio-Emitting Tidal Disruption Events in the VAST Pilot Survey

Authors: Hannah Dykaar, Maria R. Drout, B. M. Gaensler, David L. Kaplan, Tara Murphy, Assaf Horesh, Akash Anumarlapudi, Dougal Dobie, Laura N. Driessen, Emil Lenc, Adam Stewart

Abstract: We present a systematic search for tidal disruption events (TDEs) using radio data from the Variables and Slow Transients (VAST) Pilot Survey conducted using the Australian Square Kilometre Array Pathfinder (ASKAP). Historically, TDEs have been identified using observations at X-ray, optical, and ultraviolet wavelengths. After discovery, a few dozen TDEs have been shown to have radio counterparts… ▽ More We present a systematic search for tidal disruption events (TDEs) using radio data from the Variables and Slow Transients (VAST) Pilot Survey conducted using the Australian Square Kilometre Array Pathfinder (ASKAP). Historically, TDEs have been identified using observations at X-ray, optical, and ultraviolet wavelengths. After discovery, a few dozen TDEs have been shown to have radio counterparts through follow-up observations. With systematic time-domain radio surveys becoming available, we can now identify new TDEs in the radio regime. A population of radio-discovered TDEs has the potential to provide several key insights including an independent constraint on their volumetric rate. We conducted a search to select variable radio sources with a single prominent radio flare and a position consistent within 2$σ$ of the nucleus of a known galaxy. While TDEs were the primary target of our search, sources identified in this search may also be consistent with active galactic nuclei exhibiting unusual flux density changes at the timescales probed, uncharacteristically bright supernovae, or a population of gamma-ray bursts. We identify a sample of 12 radio-bright candidate TDEs. The timescales and luminosities range from ~6 to 230 days and ~10$^{38}$ to 10$^{41}$ erg s$^{-1}$, consistent with models of radio emission from TDEs that launch relativistic jets. After calculating the detection efficiency of our search using a Monte Carlo simulation of TDEs, and assuming all 12 sources are jetted TDEs, we derive a volumetric rate for jetted TDEs of 0.80$^{+0.31}_{-0.23}$ Gpc$^{-3}$ yr$^{-1}$, consistent with previous empirically estimated rates. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 37 pages, 37 figures, accepted to ApJ

Wrong Signs are Alright

Authors: David E. Kaplan, Surjeet Rajendran, Francesco Serra

Abstract: It has been shown that some Lorentz-invariant quantum field theories, such as those with higher-dimensional operators with negative coefficients, lead to superluminality on some classical backgrounds. While superluminality by itself is not logically inconsistent, these theories also predict the formation of closed time-like curves at the classical level, starting from initial conditions without su… ▽ More It has been shown that some Lorentz-invariant quantum field theories, such as those with higher-dimensional operators with negative coefficients, lead to superluminality on some classical backgrounds. While superluminality by itself is not logically inconsistent, these theories also predict the formation of closed time-like curves at the classical level, starting from initial conditions without such curves. This leads to the formation of a Cauchy Horizon which prevents a complete description of the time evolution of such systems. Inspired by the chronology protection arguments of General Relativity, we show that quantum mechanical effects from low energy quanta strongly backreact on such configurations, exciting unknown short-distance degrees of freedom and invalidating the classical predictions. Thus, there is no obvious low-energy obstruction to the existence of these operators. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 5 pages, 2 figures

Even Integer Quantum Hall Effect in Materials with Hidden Spin Texture

Authors: Daniel Kaplan, Ady Stern, Binghai Yan

Abstract: Because spin-orbit coupling (SOC) is invisible in the band structure when inversion symmetry exists, whether spins are trivially degenerate or strongly coupled to momentum due to SOC is presumed to make little difference in transport measurements, such as magnetoresistance and quantum oscillations. In this work, however, we show that hidden Rashba SOC in a centrosymmetric two-dimensional material… ▽ More Because spin-orbit coupling (SOC) is invisible in the band structure when inversion symmetry exists, whether spins are trivially degenerate or strongly coupled to momentum due to SOC is presumed to make little difference in transport measurements, such as magnetoresistance and quantum oscillations. In this work, however, we show that hidden Rashba SOC in a centrosymmetric two-dimensional material can lead to the quantum Hall effect with only even-integer plateaus, unlike a spinless electron gas. Here, two Rashba layers that are degenerate but with opposite SOC due to inversion symmetry, hybridize with each other and create two doubly-degenerate bands with hidden spin texture. Correspondingly, two branches of Landau levels interact, resulting in significant suppression of spin splitting due to the balancing of intralayer SOC and interlayer hybridization. Furthermore, we show that breaking inversion symmetry restores the ordinary quantum Hall fluid by introducing spin-split Fermi surfaces. Our theory can apply to centrosymmetric materials with strong SOC, as demonstrated in a recent experiment on the two-dimensional semiconductor Bi$_2$O$_2$Se. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: 11 pages + appendices

Electrical control of intrinsic nonlinear Hall effect in antiferromagnetic topological insulator sandwiches

Authors: Ruobing Mei, Daniel Kaplan, Binghai Yan, Cui-Zu Chang, Chao-Xing Liu

Abstract: Nonlinear Hall effect (NHE) can originate from the quantum metric mechanism in antiferromagnetic topological materials with PT symmetry, which has been experimentally observed in MnBi2Te4. In this work, we propose that breaking PT symmetry via external electric fields can lead to a dramatic enhancement of NHE, thus allowing for an electric control of NHE. Microscopically, this is because breaking… ▽ More Nonlinear Hall effect (NHE) can originate from the quantum metric mechanism in antiferromagnetic topological materials with PT symmetry, which has been experimentally observed in MnBi2Te4. In this work, we propose that breaking PT symmetry via external electric fields can lead to a dramatic enhancement of NHE, thus allowing for an electric control of NHE. Microscopically, this is because breaking PT symmetry can lift spin degeneracy of a Kramers' pair, giving rise to additional contributions within one Kramers' pair of bands. We demonstrate this enhancement through a model Hamiltonian that describes an antiferromagnetic topological insulator sandwich structure. △ Less

Submitted 4 June, 2024; originally announced June 2024.

The NANOGrav 15 yr Data Set: Chromatic Gaussian Process Noise Models for Six Pulsars

Authors: Bjorn Larsen, Chiara M. F. Mingarelli, Jeffrey S. Hazboun, Aurelien Chalumeau, Deborah C. Good, Joseph Simon, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Paul R. Brook, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels, Peter A. Gentile, Joseph Glaser, Ross J. Jennings , et al. (39 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio-frequency. However, pulse arrival times are also affected by radio-frequency dependent "chromatic" noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characteriz… ▽ More Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio-frequency. However, pulse arrival times are also affected by radio-frequency dependent "chromatic" noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characterization of GW signals may be influenced by the choice of chromatic noise model for each pulsar. To better understand this effect, we assess if and how different chromatic noise models affect achromatic noise properties in each pulsar. The models we compare include existing DM models used by NANOGrav and noise models used for the European PTA Data Release 2 (EPTA DR2). We perform this comparison using a subsample of six pulsars from the NANOGrav 15 yr data set, selecting the same six pulsars as from the EPTA DR2 six-pulsar dataset. We find that the choice of chromatic noise model noticeably affects the achromatic noise properties of several pulsars. This is most dramatic for PSR J1713+0747, where the amplitude of its achromatic red noise lowers from $\log_{10}A_{\text{RN}} = -14.1^{+0.1}_{-0.1}$ to $-14.7^{+0.3}_{-0.5}$, and the spectral index broadens from $γ_{\text{RN}} = 2.6^{+0.5}_{-0.4}$ to $γ_{\text{RN}} = 3.5^{+1.2}_{-0.9}$. We also compare each pulsar's noise properties with those inferred from the EPTA DR2, using the same models. From the discrepancies, we identify potential areas where the noise models could be improved. These results highlight the potential for custom chromatic noise models to improve PTA sensitivity to GWs. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Searching for gravitational wave optical counterparts with the Zwicky Transient Facility: summary of O4a

Authors: Tomás Ahumada, Shreya Anand, Michael W. Coughlin, Vaidehi Gupta, Mansi M. Kasliwal, Viraj R. Karambelkar, Robert D. Stein, Gaurav Waratkar, Vishwajeet Swain, Theophile Jegou du Laz, Akash Anumarlapudi, Igor Andreoni, Mattia Bulla, Gokul P. Srinivasaragavan, Andrew Toivonen, Avery Wold, Eric C. Bellm, S. Bradley Cenko, David L. Kaplan, Jesper Sollerman, Varun Bhalerao, Daniel Perley, Anirudh Salgundi, Aswin Suresh, K-Ryan Hinds , et al. (27 additional authors not shown)

Abstract: During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network (IGWN), the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star-black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (FAR < 1 per year) BNS and NSBH candida… ▽ More During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network (IGWN), the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star-black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (FAR < 1 per year) BNS and NSBH candidates in O4a. Our follow-up campaigns relied on both target-of-opportunity observations (ToO) and re-weighting of the nominal survey schedule to maximize coverage. We describe the toolkit we have been developing, Fritz, an instance of SkyPortal, instrumental in coordinating and managing our telescope scheduling, candidate vetting, and follow-up observations through a user-friendly interface. ZTF covered a total of 2841 deg$^2$ within the skymaps of the high-significance GW events, reaching a median depth of g~20.2 mag. We circulated 15 candidates, but found no viable KN counterpart to any of the GW events. Based on the ZTF non-detections of the high-significance events in O4a, we used a Bayesian approach, nimbus, to quantify the posterior probability of KN model parameters that are consistent with our non-detections. Our analysis favors KNe with initial absolute magnitude fainter than -16 mag. The joint posterior probability of a GW170817-like KN associated with all our O4a follow-ups was 64%. Additionally, we use a survey simulation software, simsurvey, to determine that our combined filtered efficiency to detect a GW170817-like KN is 36%, when considering the 5 confirmed astrophysical events in O3 (1 BNS and 4 NSBH), along with our O4a follow-ups. Following Kasliwal et al. (2020), we derived joint constraints on the underlying KN luminosity function based on our O3 and O4a follow-ups, determining that no more than 76% of KNe fading at 1 mag/day can peak at a magnitude brighter than -17.5 mag. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: submitted

Tunable Spatiotemporal Orders in Driven Insulators

Authors: Daniel Kaplan, Pavel A. Volkov, Ahana Chakraborty, Zekun Zhuang, Premala Chandra

Abstract: We show that driving optical phonons above a threshold fluence induces spatiotemporal orders, where material properties oscillate at an incommensurate wavevector $q_0$ in space and at half the drive frequency in time. The order is robust against temperature on timescales much larger than the lifetime of the excited modes and can be accompanied by a static $2q_0$ modulation. We make predictions for… ▽ More We show that driving optical phonons above a threshold fluence induces spatiotemporal orders, where material properties oscillate at an incommensurate wavevector $q_0$ in space and at half the drive frequency in time. The order is robust against temperature on timescales much larger than the lifetime of the excited modes and can be accompanied by a static $2q_0$ modulation. We make predictions for time-resolved diffraction and provide estimates for candidate materials. Our results show the possibility of using THz waves in solids to realize tunable incommensurate order on the nanoscale. △ Less

Submitted 24 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 5 pages, 4 figures + supplementary information

Cosmological Consequences of Unconstrained Gravity and Electromagnetism

Authors: Loris Del Grosso, David E Kaplan, Tom Melia, Vivian Poulin, Surjeet Rajendran, Tristan L Smith

Abstract: Motivated by the quantum description of gauge theories, we study the cosmological effects of relaxing the Hamiltonian and momentum constraints in general relativity and Gauss' law in electromagnetism. We show that the unconstrained theories have new source terms that mimic a pressureless dust and a charge density that only follows geodesics. The source terms may be the simplest explanation for dar… ▽ More Motivated by the quantum description of gauge theories, we study the cosmological effects of relaxing the Hamiltonian and momentum constraints in general relativity and Gauss' law in electromagnetism. We show that the unconstrained theories have new source terms that mimic a pressureless dust and a charge density that only follows geodesics. The source terms may be the simplest explanation for dark matter and generically predict a charged component. We comment that discovery of such terms would rule out inflation and be a direct probe of the initial conditions of the universe. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 6 pages

PINT: Maximum-likelihood estimation of pulsar timing noise parameters

Authors: Abhimanyu Susobhanan, David Kaplan, Anne Archibald, Jing Luo, Paul Ray, Timothy Pennucci, Scott Ransom, Gabriella Agazie, William Fiore, Bjorn Larsen, Patrick O'Neill, Rutger van Haasteren, Akash Anumarlapudi, Matteo Bachetti, Deven Bhakta, Chloe Champagne, H. Thankful Cromartie, Paul Demorest, Ross Jennings, Matthew Kerr, Sasha Levina, Alexander McEwen, Brent Shapiro-Albert, Joseph Swiggum

Abstract: PINT is a pure-Python framework for high-precision pulsar timing developed on top of widely used and well-tested Python libraries, supporting both interactive and programmatic data analysis workflows. We present a new frequentist framework within PINT to characterize the single-pulsar noise processes present in pulsar timing datasets. This framework enables the parameter estimation for both uncorr… ▽ More PINT is a pure-Python framework for high-precision pulsar timing developed on top of widely used and well-tested Python libraries, supporting both interactive and programmatic data analysis workflows. We present a new frequentist framework within PINT to characterize the single-pulsar noise processes present in pulsar timing datasets. This framework enables the parameter estimation for both uncorrelated and correlated noise processes as well as the model comparison between different timing and noise models in a computationally inexpensive way. We demonstrate the efficacy of the new framework by applying it to simulated datasets as well as a real dataset of PSR B1855+09. We also describe the new features implemented in PINT since it was first described in the literature. △ Less

Submitted 19 June, 2024; v1 submitted 3 May, 2024; originally announced May 2024.

Comments: Accepted for publication in ApJ

The Sydney Radio Star Catalogue: properties of radio stars at megahertz to gigahertz frequencies

Authors: Laura N. Driessen, Joshua Pritchard, Tara Murphy, George Heald, Jan Robrade, Barnali Das, Stefan Duchesne, David L. Kaplan, Emil Lenc, Christene R. Lynch, Benjamin J. S. Pope, Kovi Rose, Beate Stelzer, Yuanming Wang, Andrew Zic

Abstract: We present the Sydney Radio Star Catalogue, a new catalogue of stars detected at megahertz to gigahertz radio frequencies. It consists of 839 unique stars with 3,405 radio detections, more than doubling the previously known number of radio stars. We have included stars from large area searches for radio stars found using circular polarisation searches, cross-matching, variability searches, and pro… ▽ More We present the Sydney Radio Star Catalogue, a new catalogue of stars detected at megahertz to gigahertz radio frequencies. It consists of 839 unique stars with 3,405 radio detections, more than doubling the previously known number of radio stars. We have included stars from large area searches for radio stars found using circular polarisation searches, cross-matching, variability searches, and proper motion searches as well as presenting hundreds of newly detected stars from our search of Australian SKA Pathfinder observations. The focus of this first version of the catalogue is on objects detected in surveys using SKA precursor instruments; however we will expand this scope in future versions. The 839 objects in the Sydney Radio Star Catalogue are distributed across the whole sky and range from ultracool dwarfs to Wolf-Rayet stars. We find that the radio luminosities of cool dwarfs are lower than the radio luminosities of more evolved sub-giant and giant stars. We use X-ray detections of 530 radio stars by the eROSITA soft X-ray instrument onboard the SRG spacecraft to show that almost all of the radio stars in the catalogue are over-luminous in the radio, indicating that the majority of stars at these radio frequencies are coherent radio emitters. The Sydney Radio Star Catalogue can be found in Vizier or at https://radiostars.org. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 21 pages, 8 Figures, 7 Tables, Submitted to PASA

The NANOGrav 15 yr Data Set: Looking for Signs of Discreteness in the Gravitational-wave Background

Authors: Gabriella Agazie, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Lucas Brown, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Heling Deng, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels , et al. (58 additional authors not shown)

Abstract: The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analyt… ▽ More The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 realizations we study the populations' characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power-law. The first, at $2 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 σ- 1.9 σ$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 σ- 2.1 σ$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the $26$~nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 10 pages, 8 figures, 1 appendix, submitted to ApJ

Even-integer Quantum Hall Effect in an Oxide Caused by Hidden Rashba Effect

Authors: Jingyue Wang, Junwei Huang, Daniel Kaplan, Xuehan Zhou, Congwei Tan, Jing Zhang, Gangjian Jin, Xuzhong Cong, Yongchao Zhu, Xiaoyin Gao, Yan Liang, Huakun Zuo, Zengwei Zhu, Ruixue Zhu, Ady Stern, Hongtao Liu, Peng Gao, Binghai Yan, Hongtao Yuan, Hailin Peng

Abstract: In the presence of high magnetic field, quantum Hall systems usually host both even- and odd-integer quantized states because of lifted band degeneracies. Selective control of these quantized states is challenging but essential to understand the exotic ground states and manipulate the spin textures. Here, we study the quantum Hall effect in Bi2O2Se thin films. In magnetic fields as high as 50 T, w… ▽ More In the presence of high magnetic field, quantum Hall systems usually host both even- and odd-integer quantized states because of lifted band degeneracies. Selective control of these quantized states is challenging but essential to understand the exotic ground states and manipulate the spin textures. Here, we study the quantum Hall effect in Bi2O2Se thin films. In magnetic fields as high as 50 T, we observe only even-integer quantum Hall states, but no sign of odd-integer states. However, when reducing the thickness of the epitaxial Bi2O2Se film to one unit cell, we observe both odd- and even-integer states in this Janus (asymmetric) film grown on SrTiO3. By means of a Rashba bilayer model based on ab initio band structures of Bi2O2Se thin films, we can ascribe the absence of odd-integer states in thicker films to the hidden Rasbha effect, where the local inversion symmetry breaking in two sectors of the [Bi2O2]2+ layer yields opposite Rashba spin polarizations, which compensate with each other. In the one unit cell Bi2O2Se film grown on SrTiO3, the asymmetry introduced by top surface and bottom interface induces a net polar field. The resulting global Rashba effect lifts the band degeneracies present in the symmetric case of thicker films. △ Less

Submitted 28 June, 2024; v1 submitted 31 March, 2024; originally announced April 2024.

Comments: 6 Figures, 23 pages

Sequential Inference of Hospitalization Electronic Health Records Using Probabilistic Models

Authors: Alan D. Kaplan, Priyadip Ray, John D. Greene, Vincent X. Liu

Abstract: In the dynamic hospital setting, decision support can be a valuable tool for improving patient outcomes. Data-driven inference of future outcomes is challenging in this dynamic setting, where long sequences such as laboratory tests and medications are updated frequently. This is due in part to heterogeneity of data types and mixed-sequence types contained in variable length sequences. In this work… ▽ More In the dynamic hospital setting, decision support can be a valuable tool for improving patient outcomes. Data-driven inference of future outcomes is challenging in this dynamic setting, where long sequences such as laboratory tests and medications are updated frequently. This is due in part to heterogeneity of data types and mixed-sequence types contained in variable length sequences. In this work we design a probabilistic unsupervised model for multiple arbitrary-length sequences contained in hospitalization Electronic Health Record (EHR) data. The model uses a latent variable structure and captures complex relationships between medications, diagnoses, laboratory tests, neurological assessments, and medications. It can be trained on original data, without requiring any lossy transformations or time binning. Inference algorithms are derived that use partial data to infer properties of the complete sequences, including their length and presence of specific values. We train this model on data from subjects receiving medical care in the Kaiser Permanente Northern California integrated healthcare delivery system. The results are evaluated against held-out data for predicting the length of sequences and presence of Intensive Care Unit (ICU) in hospitalization bed sequences. Our method outperforms a baseline approach, showing that in these experiments the trained model captures information in the sequences that is informative of their future values. △ Less

Submitted 24 April, 2024; v1 submitted 27 March, 2024; originally announced March 2024.

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

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

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

Submitted 11 March, 2024; originally announced March 2024.

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

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

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

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

Submitted 5 March, 2024; originally announced March 2024.

Comments: 47 pages, 41 figures

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

The GMRT High-Resolution Southern Sky Survey for pulsars and transients -- VII: Timing of Spider MSP J1242-4712, A Bridge Between Redback and Black Widow Pulsars

Authors: Ankita Ghosh, Bhaswati Bhattacharyya, Andrew Lyne, David L. Kaplan, Jayanta Roy, Paul S. Ray, Ben Stappers, Sangita Kumari, Shubham Singh, Rahul Sharan

Abstract: We present the timing solution for the 5.31-ms spider millisecond pulsar (MSP) J1242-4712, discovered with the GMRT. PSR J1242-4712 orbits a companion of minimum mass 0.08 M$_{\odot}$ with an orbital period of 7.7 hrs and occupies a relatively unexplored region in the orbital period versus companion mass space. We did not detect gamma-ray pulsations for this MSP, and also could not identify the op… ▽ More We present the timing solution for the 5.31-ms spider millisecond pulsar (MSP) J1242-4712, discovered with the GMRT. PSR J1242-4712 orbits a companion of minimum mass 0.08 M$_{\odot}$ with an orbital period of 7.7 hrs and occupies a relatively unexplored region in the orbital period versus companion mass space. We did not detect gamma-ray pulsations for this MSP, and also could not identify the optical counterpart for PSR J1242-4712 in available optical/near-infrared data. The profile of J1242-4712 evolves with frequency showing a clear single component at lower frequencies and a three-component profile at 650 MHz. PSR J1242-4712 eclipses for a very short duration near superior conjunction (orbital phase ~ 0.23-0.25) below 360 MHz. Moreover, significant DM delays and errors in pulse times of arrivals are observed near inferior conjunction (orbital phase ~ 0.7), along with an observed eclipse in one epoch at 650 MHz. Observed eclipses and significant orbital period variability suggest that PSR J1242-4712 is possibly not a He-WD binary, but has a semi or non-degenerate companion, indicating that this is a ``spider" MSP lying in a region between typical black widows and redbacks. This system may represent a distinct category of spider MSPs, displaying characteristics that bridge the gap between known black widow and redback MSPs. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 13 pages, 8 figures, 2 tables. Accepted for publication in The Astrophysical Journal

On the Improvement of Predictive Modeling Using Bayesian Stacking and Posterior Predictive Checking

Authors: Mariana Nold, Florian Meinfelder, David Kaplan

Abstract: Model uncertainty is pervasive in real world analysis situations and is an often-neglected issue in applied statistics. However, standard approaches to the research process do not address the inherent uncertainty in model building and, thus, can lead to overconfident and misleading analysis interpretations. One strategy to incorporate more flexible models is to base inferences on predictive modeli… ▽ More Model uncertainty is pervasive in real world analysis situations and is an often-neglected issue in applied statistics. However, standard approaches to the research process do not address the inherent uncertainty in model building and, thus, can lead to overconfident and misleading analysis interpretations. One strategy to incorporate more flexible models is to base inferences on predictive modeling. This approach provides an alternative to existing explanatory models, as inference is focused on the posterior predictive distribution of the response variable. Predictive modeling can advance explanatory ambitions in the social sciences and in addition enrich the understanding of social phenomena under investigation. Bayesian stacking is a methodological approach rooted in Bayesian predictive modeling. In this paper, we outline the method of Bayesian stacking but add to it the approach of posterior predictive checking (PPC) as a means of assessing the predictive quality of those elements of the stacking ensemble that are important to the research question. Thus, we introduce a viable workflow for incorporating PPC into predictive modeling using Bayesian stacking without presuming the existence of a true model. We apply these tools to the PISA 2018 data to investigate potential inequalities in reading competency with respect to gender and socio-economic background. Our empirical example serves as rough guideline for practitioners who want to implement the concepts of predictive modeling and model uncertainty in their work to similar research questions. △ Less

Submitted 29 February, 2024; originally announced February 2024.

Comments: 40 pages including abstract and references (23 pages without), 3 figures

MSC Class: 62P25

Results of the follow-up of ANTARES neutrino alerts

Authors: A. Albert, S. Alves, M. André, M. Ardid, S. Ardid, J. -J. Aubert, J. Aublin, B. Baret, S. Basa, Y. Becherini, B. Belhorma, M. Bendahman, F. Benfenati, V. Bertin, S. Biagi, M. Bissinger, J. Boumaaza, M. Bouta, M. C. Bouwhuis, H. Brânzas, R. Bruijn, J. Brunner, J. Busto, B. Caiffi, D. Calvo , et al. (166 additional authors not shown)

Abstract: High-energy neutrinos could be produced in the interaction of charged cosmic rays with matter or radiation surrounding astrophysical sources. To look for transient sources associated with neutrino emission, a follow-up program of neutrino alerts has been operating within the ANTARES Collaboration since 2009. This program, named TAToO, has triggered robotic optical telescopes (MASTER, TAROT, ROTSE… ▽ More High-energy neutrinos could be produced in the interaction of charged cosmic rays with matter or radiation surrounding astrophysical sources. To look for transient sources associated with neutrino emission, a follow-up program of neutrino alerts has been operating within the ANTARES Collaboration since 2009. This program, named TAToO, has triggered robotic optical telescopes (MASTER, TAROT, ROTSE and the SVOM ground based telescopes) immediately after the detection of any relevant neutrino candidate and scheduled several observations in the weeks following the detection. A subset of ANTARES events with highest probabilities of being of cosmic origin has also been followed by the Swift and the INTEGRAL satellites, the Murchison Widefield Array radio telescope and the H.E.S.S. high-energy gamma-ray telescope. The results of twelve years of observations are reported. No optical counterpart has been significantly associated with an ANTARES candidate neutrino signal during image analysis. Constraints on transient neutrino emission have been set. In September 2015, ANTARES issued a neutrino alert and during the follow-up, a potential transient counterpart was identified by Swift and MASTER. A multi-wavelength follow-up campaign has allowed to identify the nature of this source and has proven its fortuitous association with the neutrino. The return of experience is particularly important for the design of the alert system of KM3NeT, the next generation neutrino telescope in the Mediterranean Sea. △ Less

Submitted 26 February, 2024; originally announced February 2024.

Comments: 27 pages, 14 figures, submitted to JCAP

Fireball anti-nucleosynthesis

Authors: Michael A. Fedderke, David E. Kaplan, Anubhav Mathur, Surjeet Rajendran, Erwin H. Tanin

Abstract: The tentative identification of approximately ten relativistic anti-helium cosmic-ray events at AMS-02 would, if confirmed, challenge our understanding of the astrophysical synthesis of heavy anti-nuclei. We propose a novel scenario for the enhanced production of such anti-nuclei that is triggered by isolated, catastrophic injections of large quantities of energetic Standard Model (SM) anti-quarks… ▽ More The tentative identification of approximately ten relativistic anti-helium cosmic-ray events at AMS-02 would, if confirmed, challenge our understanding of the astrophysical synthesis of heavy anti-nuclei. We propose a novel scenario for the enhanced production of such anti-nuclei that is triggered by isolated, catastrophic injections of large quantities of energetic Standard Model (SM) anti-quarks in our galaxy by physics beyond the Standard Model (BSM). We demonstrate that SM anti-nucleosynthetic processes that occur in the resulting rapidly expanding, thermalized fireballs of SM plasma can, for a reasonable range of parameters, produce the reported tentative $\sim 2:1$ ratio of anti-helium-3 to anti-helium-4 events at AMS-02, as well as their relativistic boosts. Moreover, we show that this can be achieved without violating anti-deuterium or anti-proton flux constraints for the appropriate anti-helium fluxes. A plausible BSM paradigm for the catastrophic injections is the collision of macroscopic composite dark-matter objects carrying large net anti-baryon number. Such a scenario would require these objects to be cosmologically stable, but to destabilize upon collision, promptly releasing a fraction of their mass energy into SM anti-particles within a tiny volume. We show that, in principle, the injection rate needed to attain the necessary anti-helium fluxes and the energetic conditions required to seed the fireballs appear possible to obtain in such a paradigm. We leave open the question of constructing a BSM particle physics model to realize this, but we suggest two concrete scenarios as promising targets for further investigation. △ Less

Submitted 18 June, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: 42 pages, 10 figures, journal version

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

A Case for a Binary Black Hole System Revealed via Quasi-Periodic Outflows

Authors: Dheeraj R. Pasham, Francesco Tombesi, Petra Sukova, Michal Zajacek, Suvendu Rakshit, Eric Coughlin, Peter Kosec, Vladimir Karas, Megan Masterson, Andrew Mummery, Thomas W. -S. Holoien, Muryel Guolo, Jason Hinkle, Bart Ripperda, Vojtech Witzany, Ben Shappee, Erin Kara, Assaf Horesh, Sjoert van Velzen, Itai Sfaradi, David L. Kaplan, Noam Burger, Tara Murphy, Ronald Remillard, James F. Steiner , et al. (11 additional authors not shown)

Abstract: Binaries containing a compact object orbiting a supermassive black hole are thought to be precursors of gravitational wave events, but their identification has been extremely challenging. Here, we report quasi-periodic variability in X-ray absorption which we interpret as quasi-periodic outflows (QPOuts) from a previously low-luminosity active galactic nucleus after an outburst, likely caused by a… ▽ More Binaries containing a compact object orbiting a supermassive black hole are thought to be precursors of gravitational wave events, but their identification has been extremely challenging. Here, we report quasi-periodic variability in X-ray absorption which we interpret as quasi-periodic outflows (QPOuts) from a previously low-luminosity active galactic nucleus after an outburst, likely caused by a stellar tidal disruption. We rule out several models based on observed properties and instead show using general relativistic magnetohydrodynamic simulations that QPOuts, separated by roughly 8.3 days, can be explained with an intermediate-mass black hole secondary on a mildly eccentric orbit at a mean distance of about 100 gravitational radii from the primary. Our work suggests that QPOuts could be a new way to identify intermediate/extreme-mass ratio binary candidates. △ Less

Submitted 15 February, 2024; originally announced February 2024.

Comments: Accepted for publication in Science Advances. We report a new supermassive black hole phenomenon that we call quasi-periodic outflows (QPOuts)

Low-Rank Extragradient Methods for Scalable Semidefinite Optimization

Authors: Dan Garber. Atara Kaplan

Abstract: We consider several classes of highly important semidefinite optimization problems that involve both a convex objective function (smooth or nonsmooth) and additional linear or nonlinear smooth and convex constraints, which are ubiquitous in statistics, machine learning, combinatorial optimization, and other domains. We focus on high-dimensional and plausible settings in which the problem admits a… ▽ More We consider several classes of highly important semidefinite optimization problems that involve both a convex objective function (smooth or nonsmooth) and additional linear or nonlinear smooth and convex constraints, which are ubiquitous in statistics, machine learning, combinatorial optimization, and other domains. We focus on high-dimensional and plausible settings in which the problem admits a low-rank solution which also satisfies a low-rank complementarity condition. We provide several theoretical results proving that, under these circumstances, the well-known Extragradient method, when initialized in the proximity of an optimal primal-dual solution, converges to a solution of the constrained optimization problem with its standard convergence rates guarantees, using only low-rank singular value decompositions (SVD) to project onto the positive semidefinite cone, as opposed to computationally-prohibitive full-rank SVDs required in worst-case. Our approach is supported by numerical experiments conducted with a dataset of Max-Cut instances. △ Less

Submitted 14 February, 2024; originally announced February 2024.

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

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

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

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

Comments: 35 pages, 20 figures

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

Intensive Swift and LCO monitoring of PG 1302$-$102: AGN disk reverberation mapping of a supermassive black hole binary candidate

Authors: Tingting Liu, Rick Edelson, Juan V. Hernández Santisteban, Erin Kara, John Montano, Jonathan Gelbord, Keith Horne, Aaron J. Barth, Edward M. Cackett, David L. Kaplan

Abstract: We present an intensive multiwavelength monitoring campaign of the quasar PG 1302$-$102 with Swift and the Las Cumbres Observatory network telescopes. At $z\sim0.3$, it tests the limits of the reverberation mapping (RM) technique in probing the accretion disk around a supermassive black hole (SMBH) and extends the parameter space to high masses and high accretion rates. This is also the first time… ▽ More We present an intensive multiwavelength monitoring campaign of the quasar PG 1302$-$102 with Swift and the Las Cumbres Observatory network telescopes. At $z\sim0.3$, it tests the limits of the reverberation mapping (RM) technique in probing the accretion disk around a supermassive black hole (SMBH) and extends the parameter space to high masses and high accretion rates. This is also the first time the RM technique has been applied to test disk structures predicted in the SMBH binary model that has been suggested for this source. PG 1302$-$102 was observed at a $\sim$daily cadence for $\sim 9$ months in 14 bands spanning from X-ray to UV and optical wavelengths, and it shows moderate to significant levels of variability correlated between wavelengths. We measure the inter-band time lags which are consistent with a $τ\propto λ^{4/3}$ relation as expected from standard disk reprocessing, albeit with large errors. The disk size implied by the lag spectrum is consistent with the expected disk size for its black hole mass within uncertainties. While the source resembles other reverberation-mapped AGN in many respects, and we do not find evidence supporting the prevalent hypothesis that it hosts an SMBH binary, we demonstrate the feasibility of studying SMBH binaries from this novel angle and suggest possibilities for the LSST Deep Drilling Fields. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: 11 pages, 5 figures, 4 tables. Accepted for publication in ApJ. The full dataset will be available with the ApJ article

Scalable High-Resolution Pixel-Space Image Synthesis with Hourglass Diffusion Transformers

Authors: Katherine Crowson, Stefan Andreas Baumann, Alex Birch, Tanishq Mathew Abraham, Daniel Z. Kaplan, Enrico Shippole

Abstract: We present the Hourglass Diffusion Transformer (HDiT), an image generative model that exhibits linear scaling with pixel count, supporting training at high-resolution (e.g. $1024 \times 1024$) directly in pixel-space. Building on the Transformer architecture, which is known to scale to billions of parameters, it bridges the gap between the efficiency of convolutional U-Nets and the scalability of… ▽ More We present the Hourglass Diffusion Transformer (HDiT), an image generative model that exhibits linear scaling with pixel count, supporting training at high-resolution (e.g. $1024 \times 1024$) directly in pixel-space. Building on the Transformer architecture, which is known to scale to billions of parameters, it bridges the gap between the efficiency of convolutional U-Nets and the scalability of Transformers. HDiT trains successfully without typical high-resolution training techniques such as multiscale architectures, latent autoencoders or self-conditioning. We demonstrate that HDiT performs competitively with existing models on ImageNet $256^2$, and sets a new state-of-the-art for diffusion models on FFHQ-$1024^2$. △ Less

Submitted 21 January, 2024; originally announced January 2024.

Comments: 20 pages, 13 figures, project page and code available at https://crowsonkb.github.io/hourglass-diffusion-transformers/

A foundation model for atomistic materials chemistry

Authors: Ilyes Batatia, Philipp Benner, Yuan Chiang, Alin M. Elena, Dávid P. Kovács, Janosh Riebesell, Xavier R. Advincula, Mark Asta, Matthew Avaylon, William J. Baldwin, Fabian Berger, Noam Bernstein, Arghya Bhowmik, Samuel M. Blau, Vlad Cărare, James P. Darby, Sandip De, Flaviano Della Pia, Volker L. Deringer, Rokas Elijošius, Zakariya El-Machachi, Fabio Falcioni, Edvin Fako, Andrea C. Ferrari, Annalena Genreith-Schriever , et al. (51 additional authors not shown)

Abstract: Machine-learned force fields have transformed the atomistic modelling of materials by enabling simulations of ab initio quality on unprecedented time and length scales. However, they are currently limited by: (i) the significant computational and human effort that must go into development and validation of potentials for each particular system of interest; and (ii) a general lack of transferabilit… ▽ More Machine-learned force fields have transformed the atomistic modelling of materials by enabling simulations of ab initio quality on unprecedented time and length scales. However, they are currently limited by: (i) the significant computational and human effort that must go into development and validation of potentials for each particular system of interest; and (ii) a general lack of transferability from one chemical system to the next. Here, using the state-of-the-art MACE architecture we introduce a single general-purpose ML model, trained on a public database of 150k inorganic crystals, that is capable of running stable molecular dynamics on molecules and materials. We demonstrate the power of the MACE-MP-0 model - and its qualitative and at times quantitative accuracy - on a diverse set problems in the physical sciences, including the properties of solids, liquids, gases, chemical reactions, interfaces and even the dynamics of a small protein. The model can be applied out of the box and as a starting or "foundation model" for any atomistic system of interest and is thus a step towards democratising the revolution of ML force fields by lowering the barriers to entry. △ Less

Submitted 1 March, 2024; v1 submitted 29 December, 2023; originally announced January 2024.

Comments: 119 pages, 63 figures, 37MB PDF

Multi-epoch sampling of the radio star population with the Australian SKA Pathfinder

Authors: Joshua Pritchard, Tara Murphy, George Heald, Michael S. Wheatland, David L. Kaplan, Emil Lenc, Andrew O'Brien, Ziteng Wang

Abstract: The population of radio-loud stars has to date been studied primarily through either targeted observations of a small number of highly active stars or widefield, single-epoch surveys that cannot easily distinguish stellar emission from background extra-Galactic sources. As a result it has been difficult to constrain population statistics such as the surface density and fraction of the population p… ▽ More The population of radio-loud stars has to date been studied primarily through either targeted observations of a small number of highly active stars or widefield, single-epoch surveys that cannot easily distinguish stellar emission from background extra-Galactic sources. As a result it has been difficult to constrain population statistics such as the surface density and fraction of the population producing radio emission in a particular variable or spectral class. In this paper we present a sample of 36 radio stars detected in a circular polarisation search of the multi-epoch Variables and Slow Transients (VAST) pilot survey with ASKAP at 887.5~MHz. Through repeat sampling of the VAST pilot survey footprint we find an upper limit to the duty cycle of M-dwarf radio bursts of 8.5 per cent, and that at least $10 \pm 3$ per cent of the population should produce radio bursts more luminous than $10^{15}$ erg s$^{-1}$ Hz$^{-1}$. We infer a lower limit on the long-term surface density of such bursts in a shallow 1.25 mJy PSF$^{-1}$ sensitivity survey of $9^{+11}_{-7} \times 10^{-3}$ deg$^{-2}$ and an instantaneous radio star surface density of $1.7 \pm 0.2 \times 10^{-3}$ deg$^{-2}$ on 12 min timescales. Based on these rates we anticipate ${\sim}200 \pm 50$ new radio star detections per year over the full VAST survey and ${\sim}41\,000^{+10\,000}_{-9\,000}$ in next-generation all-sky surveys with the Square Kilometre Array. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: 14 pages, 9 figures, submitted to MNRAS

The Green Bank North Celestial Cap Survey IX: Timing Follow-up for 128 Pulsars

Authors: A. E. McEwen, J. K. Swiggum, D. L. Kaplan, C. M. Tan, B. W. Meyers, E. Fonseca, G. Y. Agazie, P. Chawla, K. Crowter, M. E. DeCesar, T. Dolch, F. A. Dong, W. Fiore, E. Fonseca, D. C. Good, A. G. Istrate, V. M. Kaspi, V. I. Kondratiev, J. van Leeuwen, L. Levin, E. F. Lewis, R. S. Lynch, K. W. Masui, J. W. McKee, M. A. McLaughlin , et al. (6 additional authors not shown)

Abstract: The Green Bank North Celestial Cap survey is one of the largest and most sensitive searches for pulsars and transient radio objects. Observations for the survey have finished; priorities have shifted toward long-term monitoring of its discoveries. In this study, we have developed a pipeline to handle large datasets of archival observations and connect them to recent, high-cadence observations take… ▽ More The Green Bank North Celestial Cap survey is one of the largest and most sensitive searches for pulsars and transient radio objects. Observations for the survey have finished; priorities have shifted toward long-term monitoring of its discoveries. In this study, we have developed a pipeline to handle large datasets of archival observations and connect them to recent, high-cadence observations taken using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. This pipeline handles data for 128 pulsars and has produced measurements of spin, positional, and orbital parameters that connect data over observation gaps as large as 2000 days. We have also measured glitches in the timing residuals for five of the pulsars included and proper motion for 19 sources (13 new). We include updates to orbital parameters for 19 pulsars, including 9 previously unpublished binaries. For two of these binaries, we provide updated measurements of post-Keplerian binary parameters, which result in much more precise estimates of the total masses of both systems. For PSR J0509+3801, the much improved measurement of the Einstein delay yields much improved mass measurements for the pulsar and its companion, 1.399(6)\Msun and 1.412(6)\Msun, respectively. For this system, we have also obtained a measurement of the orbital decay due to the emission of gravitational waves: $\dot{P}_{\rm B} = -1.37(7)\times10^{-12}$, which is in agreement with the rate predicted by general relativity for these masses. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: accepted for publication in The Astrophysical Journal

Weyl fermions on a finite lattice

Authors: David B. Kaplan, Srimoyee Sen

Abstract: The phenomenon of unpaired Weyl fermions appearing on the sole 2n-dimensional boundary of a (2n+1)-dimensional manifold with massive Dirac fermions was recently analyzed in a companion paper by one of the authors. In this Letter we show that similar unpaired Weyl edge states can be seen on a finite lattice. In particular, we consider the discretized Hamiltonian for a Wilson fermion in (2+1) dimens… ▽ More The phenomenon of unpaired Weyl fermions appearing on the sole 2n-dimensional boundary of a (2n+1)-dimensional manifold with massive Dirac fermions was recently analyzed in a companion paper by one of the authors. In this Letter we show that similar unpaired Weyl edge states can be seen on a finite lattice. In particular, we consider the discretized Hamiltonian for a Wilson fermion in (2+1) dimensions with a 1+1 dimensional boundary and continuous time. We demonstrate that the low lying boundary spectrum is indeed Weyl-like: it has a linear dispersion relation and definite chirality and circulates in only one direction around the boundary. We comment on how our results are consistent with Nielsen-Ninomiya theorem. This work removes one potential obstacle facing the program for regulating chiral gauge theories recently proposed by one of the authors. △ Less

Submitted 13 March, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: 5 pages, 4 figures. This version coincides with the version accepted for publication in Phys. Rev. Lett. on Jan. 26, 2024. Major changes include performing computations on a disk-shaped lattice instead of square, with corresponding changes in the figures. Supplemental material of the publication version has been included here as an appendix. Scientific conclusions have not changed

Report number: INT-PUB-23-046

The DUNE Far Detector Vertical Drift Technology, Technical Design Report

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

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

Submitted 5 December, 2023; originally announced December 2023.

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

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

Chiral gauge theory at the boundary between topological phases

Authors: David B. Kaplan

Abstract: I demonstrate how chiral fermions with an exact gauge symmetry can appear on the d-dimensional boundary of a finite volume (d+1)-dimensional manifold, without any light mirror partners. The condition for the d-dimensional boundary theory to be local is that gauge anomalies cancel and that the volume be large. This can likely be achieved on a lattice and provides a new paradigm for the lattice regu… ▽ More I demonstrate how chiral fermions with an exact gauge symmetry can appear on the d-dimensional boundary of a finite volume (d+1)-dimensional manifold, without any light mirror partners. The condition for the d-dimensional boundary theory to be local is that gauge anomalies cancel and that the volume be large. This can likely be achieved on a lattice and provides a new paradigm for the lattice regularization of chiral gauge theories. △ Less

Submitted 13 March, 2024; v1 submitted 3 December, 2023; originally announced December 2023.

Comments: V3: updated to version accepted for publication by Phys. Rev. Lett. on Jan. 21, 2024. Rendered more concise, no changes in the physics

Report number: INT-PUB-23-045

Discovery of radio eclipses from 4FGL J1646.5$-$4406: a new candidate redback pulsar binary

Authors: Andrew Zic, Ziteng Wang, Emil Lenc, David L. Kaplan, Tara Murphy, Alessandro Ridolfi, Rahul Sengar, Natasha Hurley-Walker, Dougal Dobie, James K. Leung, Joshua Pritchard, Yuanming Wang

Abstract: Large widefield surveys make possible the serendipitous discovery of rare sub-classes of pulsars. One such class are "spider"-type pulsar binaries, comprised of a pulsar in a compact orbit with a low-mass (sub)stellar companion. In a search for circularly-polarized radio sources in ASKAP Pilot Survey observations, we discovered highly variable and circularly polarized emission from a radio source… ▽ More Large widefield surveys make possible the serendipitous discovery of rare sub-classes of pulsars. One such class are "spider"-type pulsar binaries, comprised of a pulsar in a compact orbit with a low-mass (sub)stellar companion. In a search for circularly-polarized radio sources in ASKAP Pilot Survey observations, we discovered highly variable and circularly polarized emission from a radio source within the error region of the $γ$-ray source {4FGL}~J1646.5$-$4406. The variability is consistent with the eclipse of a compact, steep-spectrum source behind ablated material from a companion in a $\sim 5.3\,$h binary orbit. Based on the eclipse properties and spatial coincidence with {4FGL} J1646.5$-$4406, we argue that the source is likely a recycled pulsar in a "redback" binary system. Using properties of the eclipses from ASKAP and Murchison Widefield Array observations, we provide broad constraints on the properties of the eclipse medium. We identified a potential optical/infra-red counterpart in archival data consistent with a variable low-mass star. Using the Parkes Radio Telescope "Murriyang" and MeerKAT, we searched extensively for radio pulsations but yielded no viable detections of pulsed emission. We suggest that the non-detection of pulses is due to scattering in the intra-binary material, but scattering from the ISM can also plausibly explain the pulse non-detections if the interstellar dispersion measure exceeds $\sim$600$\,$pc$\,$cm$^{-3}$. Orbital constraints derived from optical observations of the counterpart would be highly valuable for future $γ$-ray pulsation searches, which may confirm the source nature as a pulsar. △ Less

Submitted 30 November, 2023; originally announced December 2023.

Comments: Submitted to MNRAS. 13 Pages, 10 figures, 3 tables

Discovery of a young, highly scattered pulsar PSR J1032-5804 with the Australian SKA Pathfinder

Authors: Ziteng Wang, David L. Kaplan, Rahul Sengar, Emil Lenc, Andrew Zic, Akash Anumarlapudi, B. M. Gaensler, Natasha Hurley-Walker, Tara Murphy, Yuanming Wang

Abstract: We report the discovery of a young, highly scattered pulsar in a search for highly circularly polarized radio sources as part of the Australian Square Kilometre Array Pathfinder (ASKAP) Variables and Slow Transients (VAST) survey. In follow-up observations with Murriyang/Parkes, we identified PSR J1032-5804 and measured a period of 78.7 ms, dispersion measure (DM) of 819$\pm$4 pc cm$^{-3}$, rotati… ▽ More We report the discovery of a young, highly scattered pulsar in a search for highly circularly polarized radio sources as part of the Australian Square Kilometre Array Pathfinder (ASKAP) Variables and Slow Transients (VAST) survey. In follow-up observations with Murriyang/Parkes, we identified PSR J1032-5804 and measured a period of 78.7 ms, dispersion measure (DM) of 819$\pm$4 pc cm$^{-3}$, rotation measure of -2000$\pm$1 rad m$^{-2}$, and a characteristic age of 34.6 kyr. We found a pulse scattering timescale at 3 GHz of ~22 ms, implying a timescale at 1 GHz of ~3845 ms, which is the third most scattered pulsar known and explains its non-detection in previous pulsar surveys. We discuss the identification of a possible pulsar wind nebula and supernova remnant in the pulsar's local environment by analyzing the pulsar spectral energy distribution and the surrounding extended emission from multiwavelength images. Our result highlights the possibility of identifying extremely scattered pulsars from radio continuum images. Ongoing and future large-scale radio continuum surveys will offer us an unprecedented opportunity to find more extreme pulsars (e.g., highly scattered, highly intermittent, highly accelerated), which will enhance our understanding of the characteristics of pulsars and the interstellar medium. △ Less

Submitted 24 November, 2023; originally announced November 2023.

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

Expanding neutrino oscillation parameter measurements in NOvA using a Bayesian approach

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

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

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

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

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

Crepant resolutions of stratified varieties via gluing

Authors: Daniel Kaplan, Travis Schedler

Abstract: Let $X$ be a variety with a stratification $\mathcal{S}$ into smooth locally closed subvarieties such that $X$ is locally a product along each stratum (e.g., a symplectic singularity). We prove that assigning to each open subset $U \subset X$ the set of isomorphism classes of locally projective crepant resolutions of $U$ defines an $\mathcal{S}$-constructible sheaf of sets. Thus, for each stratum… ▽ More Let $X$ be a variety with a stratification $\mathcal{S}$ into smooth locally closed subvarieties such that $X$ is locally a product along each stratum (e.g., a symplectic singularity). We prove that assigning to each open subset $U \subset X$ the set of isomorphism classes of locally projective crepant resolutions of $U$ defines an $\mathcal{S}$-constructible sheaf of sets. Thus, for each stratum $S$ and basepoint $s \in S$, the fundamental group acts on the set of germs of projective crepant resolutions at $s$, leaving invariant the germs extending to the entire stratum. Global locally projective crepant resolutions correspond to compatible such choices for all strata. For example, if the local projective crepant resolutions are unique, they automatically glue uniquely. We give criteria for a locally projective crepant resolution $ρ: \tilde X \to X$ to be globally projective. We show that the sheafification of the presheaf $U \mapsto \operatorname{Pic}(ρ^{-1}(U)/U)$ of relative Picard classes is also constructible. The resolution is globally projective only if there exist local relatively ample bundles whose classes glue to a global section of this sheaf. The obstruction to lifting this section to a global ample line bundle is encoded by a gerbe on the singularity $X$. We show the gerbes are automatically trivial if $X$ is a symplectic quotient singularity. Our main results hold in the more general setting of partial crepant resolutions, that need not have smooth source. We apply the theory to symmetric powers and Hilbert schemes of surfaces with du Val singularities, finite quotients of tori, multiplicative and Nakajima quiver varieties, as well as to canonical threefold singularities. △ Less

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

Comments: 32 pages, comments welcome, fixed issue with proof that locally projective resolutions glue to a globally projective resolution, added subsection 3.3 to further explore the global projectivity issue

The NANOGrav 15-year data set: Search for Transverse Polarization Modes in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy Baier, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest , et al. (74 additional authors not shown)

Abstract: Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correl… ▽ More Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correlations. In this work we search the NANOGrav 15-year data set for evidence of a gravitational wave background with quadrupolar Hellings and Downs (HD) and Scalar Transverse (ST) correlations. We find that HD correlations are the best fit to the data, and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors $\sim 2$ when comparing HD to ST correlations, and $\sim 1$ for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise-ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise-ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise-ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 11 pages, 5 figures

ForceGen: End-to-end de novo protein generation based on nonlinear mechanical unfolding responses using a protein language diffusion model

Authors: Bo Ni, David L. Kaplan, Markus J. Buehler

Abstract: Through evolution, nature has presented a set of remarkable protein materials, including elastins, silks, keratins and collagens with superior mechanical performances that play crucial roles in mechanobiology. However, going beyond natural designs to discover proteins that meet specified mechanical properties remains challenging. Here we report a generative model that predicts protein designs to m… ▽ More Through evolution, nature has presented a set of remarkable protein materials, including elastins, silks, keratins and collagens with superior mechanical performances that play crucial roles in mechanobiology. However, going beyond natural designs to discover proteins that meet specified mechanical properties remains challenging. Here we report a generative model that predicts protein designs to meet complex nonlinear mechanical property-design objectives. Our model leverages deep knowledge on protein sequences from a pre-trained protein language model and maps mechanical unfolding responses to create novel proteins. Via full-atom molecular simulations for direct validation, we demonstrate that the designed proteins are novel, and fulfill the targeted mechanical properties, including unfolding energy and mechanical strength, as well as the detailed unfolding force-separation curves. Our model offers rapid pathways to explore the enormous mechanobiological protein sequence space unconstrained by biological synthesis, using mechanical features as target to enable the discovery of protein materials with superior mechanical properties. △ Less

Submitted 15 December, 2023; v1 submitted 16 October, 2023; originally announced October 2023.

A Computational Approach to Style in American Poetry

Authors: David M. Kaplan, David M. Blei

Abstract: We develop a quantitative method to assess the style of American poems and to visualize a collection of poems in relation to one another. Qualitative poetry criticism helped guide our development of metrics that analyze various orthographic, syntactic, and phonemic features. These features are used to discover comprehensive stylistic information from a poem's multi-layered latent structure, and to… ▽ More We develop a quantitative method to assess the style of American poems and to visualize a collection of poems in relation to one another. Qualitative poetry criticism helped guide our development of metrics that analyze various orthographic, syntactic, and phonemic features. These features are used to discover comprehensive stylistic information from a poem's multi-layered latent structure, and to compute distances between poems in this space. Visualizations provide ready access to the analytical components. We demonstrate our method on several collections of poetry, showing that it better delineates poetry style than the traditional word-occurrence features that are used in typical text analysis algorithms. Our method has potential applications to academic research of texts, to research of the intuitive personal response to poetry, and to making recommendations to readers based on their favorite poems. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Comments: accepted manuscript; see doi for version of record

ACM Class: J.5; I.2.7

Journal ref: Seventh IEEE International Conference on Data Mining (ICDM 2007), Omaha, NE, USA, 2007, pp. 553-558

Smoothed instrumental variables quantile regression

Authors: David M. Kaplan

Abstract: In this article, I introduce the sivqr command, which estimates the coefficients of the instrumental variables (IV) quantile regression model introduced by Chernozhukov and Hansen (2005). The sivqr command offers several advantages over the existing ivqreg and ivqreg2 commands for estimating this IV quantile regression model, which complements the alternative "triangular model" behind cqiv and the… ▽ More In this article, I introduce the sivqr command, which estimates the coefficients of the instrumental variables (IV) quantile regression model introduced by Chernozhukov and Hansen (2005). The sivqr command offers several advantages over the existing ivqreg and ivqreg2 commands for estimating this IV quantile regression model, which complements the alternative "triangular model" behind cqiv and the "local quantile treatment effect" model of ivqte. Computationally, sivqr implements the smoothed estimator of Kaplan and Sun (2017), who show that smoothing improves both computation time and statistical accuracy. Standard errors are computed analytically or by Bayesian bootstrap; for non-iid sampling, sivqr is compatible with bootstrap. I discuss syntax and the underlying methodology, and I compare sivqr with other commands in an example. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Comments: accepted manuscript

Journal ref: Stata Journal 22 (2022) 379-403

Transverse Emittance Reduction in Muon Beams by Ionization Cooling

Authors: The MICE Collaboration, M. Bogomilov, R. Tsenov, G. Vankova-Kirilova, Y. P. Song, J. Y. Tang, Z. H. Li, R. Bertoni, M. Bonesini, F. Chignoli, R. Mazza, A. de Bari, D. Orestano, L. Tortora, Y. Kuno, H. Sakamoto, A. Sato, S. Ishimoto, M. Chung, C. K. Sung, F. Filthaut, M. Fedorov, D. Jokovic, D. Maletic, M. Savic , et al. (112 additional authors not shown)

Abstract: Accelerated muon beams have been considered for next-generation studies of high-energy lepton-antilepton collisions and neutrino oscillations. However, high-brightness muon beams have not yet been produced. The main challenge for muon acceleration and storage stems from the large phase-space volume occupied by the beam, derived from the muon production mechanism through the decay of pions from pro… ▽ More Accelerated muon beams have been considered for next-generation studies of high-energy lepton-antilepton collisions and neutrino oscillations. However, high-brightness muon beams have not yet been produced. The main challenge for muon acceleration and storage stems from the large phase-space volume occupied by the beam, derived from the muon production mechanism through the decay of pions from proton collisions. Ionization cooling is the technique proposed to decrease the muon beam phase-space volume. Here we demonstrate a clear signal of ionization cooling through the observation of transverse emittance reduction in beams that traverse lithium hydride or liquid hydrogen absorbers in the Muon Ionization Cooling Experiment (MICE). The measurement is well reproduced by the simulation of the experiment and the theoretical model. The results shown here represent a substantial advance towards the realization of muon-based facilities that could operate at the energy and intensity frontiers. △ Less

Submitted 13 October, 2023; v1 submitted 9 October, 2023; originally announced October 2023.

Comments: 23 pages and 5 figures

Report number: STFC-P-2023-004

Realistic non-collinear ground states of solids with source-free exchange correlation functional

Authors: Guy C. Moore, Matthew K. Horton, Aaron D. Kaplan, Sinéad M. Griffin, Kristin A. Persson

Abstract: In this work, we extend the source-free (SF) exchange correlation (XC) functional developed by Sangeeta Sharma and co-workers to plane-wave density functional theory (DFT) based on the projector augmented wave (PAW) method. This constraint is implemented by the current authors within the VASP source code, using a fast Poisson solver that capitalizes on the parallel three-dimensional fast Fourier t… ▽ More In this work, we extend the source-free (SF) exchange correlation (XC) functional developed by Sangeeta Sharma and co-workers to plane-wave density functional theory (DFT) based on the projector augmented wave (PAW) method. This constraint is implemented by the current authors within the VASP source code, using a fast Poisson solver that capitalizes on the parallel three-dimensional fast Fourier transforms (FFTs) implemented in VASP. Using this modified XC functional, we explore the improved convergence behavior that results from applying this constraint to the GGA-PBE+$U$+$J$ functional. In the process, we compare the non-collinear magnetic ground state computed by each functional and their SF counterpart for a select number of magnetic materials in order to provide a metric for comparing with experimentally determined magnetic orderings. We observe significantly improved agreement with experimentally measured magnetic ground state structures after applying the source-free constraint. Furthermore, we explore the importance of considering probability current densities in spin polarized systems, even under no applied field. We analyze the XC torque as well, in order to provide theoretical and computational analyses of the net XC magnetic torque induced by the source-free constraint. Along these lines, we highlight the importance of properly considering the real-space integral of the source-free local magnetic XC field. Our analysis on probability currents, net torque, and constant terms draws additional links to the rich body of previous research on spin-current density functional theory (SCDFT), and paves the way for future extensions and corrections to the SF corrected XC functional. △ Less

Submitted 24 May, 2024; v1 submitted 29 September, 2023; originally announced October 2023.

The NANOGrav 12.5-year data set: A computationally efficient eccentric binary search pipeline and constraints on an eccentric supermassive binary candidate in 3C 66B

Authors: Gabriella Agazie, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Belinda D. Cheeseboro, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Lankeswar Dey, Timothy Dolch, Justin A. Ellis, Robert D. Ferdman, Elizabeth C. Ferrara , et al. (63 additional authors not shown)

Abstract: The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for contin… ▽ More The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it on simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5-year data set employing PTA signal models containing Earth term-only as well as Earth+Pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricity $e_0<0.5$ and the symmetric mass ratio $η>0.1$. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of $88.1\pm3.7$ ns for the Earth term-only and $81.74\pm0.86$ ns for the Earth+Pulsar term searches for $e_0<0.5$ and $η>0.1$. Similar 95% upper limits on the chirp mass are $(1.98 \pm 0.05) \times 10^9\,M_{\odot}$ and $(1.81 \pm 0.01) \times 10^9\,M_{\odot}$. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5-year data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations. △ Less

Submitted 15 January, 2024; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 27 Pages, 10 Figures, 1 Table, Accepted for publication in ApJ

MWA rapid follow-up of gravitational wave transients: prospects for detecting prompt radio counterparts

Authors: J. Tian, G. E. Anderson, A. J. Cooper, K. Gourdji, M. Sokolowski, A. Rowlinson, A. Williams, G. Sleap, D. Dobie, D. L. Kaplan, Tara Murphy, S. J. Tingay, F. H. Panther, P. D. Lasky, A. Bahramian, J. C. A. Miller-Jones, C. W. James, B. W. Meyers, S. J. McSweeney, P. J. Hancock

Abstract: We present and evaluate the prospects for detecting coherent radio counterparts to gravitational wave (GW) events using Murchison Widefield Array (MWA) triggered observations. The MWA rapid-response system, combined with its buffering mode ($\sim4$ minutes negative latency), enables us to catch any radio signals produced from seconds prior to hours after a binary neutron star (BNS) merger. The lar… ▽ More We present and evaluate the prospects for detecting coherent radio counterparts to gravitational wave (GW) events using Murchison Widefield Array (MWA) triggered observations. The MWA rapid-response system, combined with its buffering mode ($\sim4$ minutes negative latency), enables us to catch any radio signals produced from seconds prior to hours after a binary neutron star (BNS) merger. The large field of view of the MWA ($\sim1000\,\text{deg}^2$ at 120\,MHz) and its location under the high sensitivity sky region of the LIGO-Virgo-KAGRA (LVK) detector network, forecast a high chance of being on-target for a GW event. We consider three observing configurations for the MWA to follow up GW BNS merger events, including a single dipole per tile, the full array, and four sub-arrays. We then perform a population synthesis of BNS systems to predict the radio detectable fraction of GW events using these configurations. We find that the configuration with four sub-arrays is the best compromise between sky coverage and sensitivity as it is capable of placing meaningful constraints on the radio emission from 12.6\% of GW BNS detections. Based on the timescales of four BNS merger coherent radio emission models, we propose an observing strategy that involves triggering the buffering mode to target coherent signals emitted prior to, during or shortly following the merger, which is then followed by continued recording for up to three hours to target later time post-merger emission. We expect MWA to trigger on $\sim5\text{--}22$ BNS merger events during the LVK O4 observing run, which could potentially result in two detections of predicted coherent emission. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Comments: Accepted for publication in PASA

Generative modeling, design and analysis of spider silk protein sequences for enhanced mechanical properties

Authors: Wei Lu, David L. Kaplan, Markus J. Buehler

Abstract: Spider silks are remarkable materials characterized by superb mechanical properties such as strength, extensibility and lightweightedness. Yet, to date, limited models are available to fully explore sequence-property relationships for analysis and design. Here we propose a custom generative large-language model to enable design of novel spider silk protein sequences to meet complex combinations of… ▽ More Spider silks are remarkable materials characterized by superb mechanical properties such as strength, extensibility and lightweightedness. Yet, to date, limited models are available to fully explore sequence-property relationships for analysis and design. Here we propose a custom generative large-language model to enable design of novel spider silk protein sequences to meet complex combinations of target mechanical properties. The model, pretrained on a large set of protein sequences, is fine-tuned on ~1,000 major ampullate spidroin (MaSp) sequences for which associated fiber-level mechanical properties exist, to yield an end-to-end forward and inverse generative strategy. Performance is assessed through: (1), a novelty analysis and protein type classification for generated spidroin sequences through BLAST searches, (2) property evaluation and comparison with similar sequences, (3) comparison of molecular structures, as well as, and (4) a detailed sequence motif analyses. We generate silk sequences with property combinations that do not exist in nature, and develop a deep understanding the mechanistic roles of sequence patterns in achieving overarching key mechanical properties (elastic modulus, strength, toughness, failure strain). The model provides an efficient approach to expand the silkome dataset, facilitating further sequence-structure analyses of silks, and establishes a foundation for synthetic silk design and optimization. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Searching for axion forces with spin precession in atoms and molecules

Authors: Prateek Agrawal, Nicholas R. Hutzler, David E. Kaplan, Surjeet Rajendran, Mario Reig

Abstract: We propose to use atoms and molecules as quantum sensors of axion-mediated monopole-dipole forces. We show that electron spin precession experiments using atomic and molecular beams are well-suited for axion searches thanks to the presence of co-magnetometer states and single-shot temporal resolution. Experimental strategies to detect axion gradients from localised sources and the earth are presen… ▽ More We propose to use atoms and molecules as quantum sensors of axion-mediated monopole-dipole forces. We show that electron spin precession experiments using atomic and molecular beams are well-suited for axion searches thanks to the presence of co-magnetometer states and single-shot temporal resolution. Experimental strategies to detect axion gradients from localised sources and the earth are presented, taking ACME III as a prototype example. Other possibilities including atomic beams, and laser-cooled atoms and molecules are discussed. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 9 pages, 2 figures. Comments welcome

Generalized Ginsparg-Wilson relations

Authors: Michael Clancy, David B. Kaplan, Hersh Singh

Abstract: We give a general derivation of Ginsparg-Wilson relations for both Dirac and Majorana fermions in any dimension. These relations encode continuous and discrete chiral, parity and time reversal anomalies and will apply to the various classes of free fermion topological insulators and superconductors (in the framework of a relativistic quantum field theory in Euclidean spacetime). We show how to for… ▽ More We give a general derivation of Ginsparg-Wilson relations for both Dirac and Majorana fermions in any dimension. These relations encode continuous and discrete chiral, parity and time reversal anomalies and will apply to the various classes of free fermion topological insulators and superconductors (in the framework of a relativistic quantum field theory in Euclidean spacetime). We show how to formulate the exact symmetries of the lattice action and the relevant index theorems for the anomalies. △ Less

Submitted 6 October, 2023; v1 submitted 15 September, 2023; originally announced September 2023.

Comments: 14 pages. Submitted Version. Section III.D revised, and other minor improvements

Report number: INT-PUB-23-024,IQuS@UW-21-064,FERMILAB-PUB-23-541-T