Showing 1–50 of 277 results for author: Miller, K

Galaxy cluster profiles: A Gaussian mixture model approach to halo miscentering

Authors: Matthew Currie, Kyle Miller, Tae-Hyeon Shin, Eric Baxter, Bhuvnesh Jain

Abstract: Measurements of the galaxy density and weak-lensing profiles of galaxy clusters typically rely on an assumed cluster center, which is taken to be the brightest cluster galaxy or other proxies for the true halo center. Departure of the assumed cluster center from the true halo center bias the resultant profile measurements, an effect known as miscentering bias. Currently, miscentering is typically… ▽ More Measurements of the galaxy density and weak-lensing profiles of galaxy clusters typically rely on an assumed cluster center, which is taken to be the brightest cluster galaxy or other proxies for the true halo center. Departure of the assumed cluster center from the true halo center bias the resultant profile measurements, an effect known as miscentering bias. Currently, miscentering is typically modeled in stacked profiles of clusters with a two parameter model. We use an alternate approach in which the profiles of individual clusters are used with the corresponding likelihood computed using a Gaussian mixture model. We test the approach using halos and the corresponding subhalo profiles from the IllustrisTNG hydrodynamic simulations. We obtain significantly improved estimates of the miscentering parameters for both 3D and projected 2D profiles relevant for imaging surveys. We discuss applications to upcoming cosmological surveys. △ Less

Submitted 19 August, 2024; originally announced August 2024.

Comments: 7 pages, 6 figures, to be submitted to A&A The Python package used in this paper is publicly available at https://github.com/KyleMiller1/Halo-Miscentering-Mixture-Model

Demonstration of hybrid foreground removal on CHIME data

Authors: Haochen Wang, Kiyoshi Masui, Kevin Bandura, Arnab Chakraborty, Matt Dobbs, Simon Foreman, Liam Gray, Mark Halpern, Albin Joseph, Joshua MacEachern, Juan Mena-Parra, Kyle Miller, Laura Newburgh, Sourabh Paul, Alex Reda, Pranav Sanghavi, Seth Siegel, Dallas Wulf

Abstract: The main challenge of 21 cm cosmology experiments is astrophysical foregrounds which are difficult to separate from the signal due to telescope systematics. An earlier study has shown that foreground residuals induced by antenna gain errors can be estimated and subtracted using the hybrid foreground residual subtraction (HyFoReS) technique which relies on cross-correlating linearly filtered data.… ▽ More The main challenge of 21 cm cosmology experiments is astrophysical foregrounds which are difficult to separate from the signal due to telescope systematics. An earlier study has shown that foreground residuals induced by antenna gain errors can be estimated and subtracted using the hybrid foreground residual subtraction (HyFoReS) technique which relies on cross-correlating linearly filtered data. In this paper, we apply a similar technique to the CHIME stacking analysis to subtract beam-induced foreground contamination. Using a linear high-pass delay filter for foreground suppression, the CHIME collaboration reported a $11.1σ$ detection in the 21 cm signal stacked on eBOSS quasar locations, despite foreground residual contamination mostly due to the instrument chromatic transfer function. We cross-correlate the foreground-dominated data at low delay with the contaminated signal at high delay to estimate residual foregrounds and subtract them from the signal. We find foreground residual subtraction can improve the signal-to-noise ratio of the stacked 21 cm signal by $ 10 - 20\%$ after the delay foreground filter, although some of the improvement can also be achieved with an alternative flagging technique. We have shown that it is possible to use HyFoReS to reduce beam-induced foreground contamination, benefiting the analysis of the HI auto power spectrum with CHIME and enabling the recovery of large scale modes. △ Less

Submitted 16 August, 2024; originally announced August 2024.

On endosplit $p$-permutation resolutions and Broué's conjecture for $p$-solvable groups

Authors: Sam K. Miller

Abstract: Endosplit $p$-permutation resolutions play an instrumental role in verifying Broué's abelian defect group conjecture in numerous cases. In this article, we give a complete classification of endosplit $p$-permutation resolutions and reduce the question of Galois descent of an endosplit $p$-permutation resolution to the Galois descent of the module it resolves. This is shown using techniques from th… ▽ More Endosplit $p$-permutation resolutions play an instrumental role in verifying Broué's abelian defect group conjecture in numerous cases. In this article, we give a complete classification of endosplit $p$-permutation resolutions and reduce the question of Galois descent of an endosplit $p$-permutation resolution to the Galois descent of the module it resolves. This is shown using techniques from the study of endotrivial complexes, the invertible objects of the bounded homotopy category of $p$-permutation modules. As an application, we show that a refinement of Broué's conjecture proposed by Kessar and Linckelmann holds for all blocks of $p$-solvable groups. △ Less

Submitted 7 August, 2024; originally announced August 2024.

Comments: 18 pages

MSC Class: 20C20; 20J05

Differentiable Optimization of Similarity Scores Between Models and Brains

Authors: Nathan Cloos, Moufan Li, Markus Siegel, Scott L. Brincat, Earl K. Miller, Guangyu Robert Yang, Christopher J. Cueva

Abstract: What metrics should guide the development of more realistic models of the brain? One proposal is to quantify the similarity between models and brains using methods such as linear regression, Centered Kernel Alignment (CKA), and angular Procrustes distance. To better understand the limitations of these similarity measures we analyze neural activity recorded in five experiments on nonhuman primates,… ▽ More What metrics should guide the development of more realistic models of the brain? One proposal is to quantify the similarity between models and brains using methods such as linear regression, Centered Kernel Alignment (CKA), and angular Procrustes distance. To better understand the limitations of these similarity measures we analyze neural activity recorded in five experiments on nonhuman primates, and optimize synthetic datasets to become more similar to these neural recordings. How similar can these synthetic datasets be to neural activity while failing to encode task relevant variables? We find that some measures like linear regression and CKA, differ from angular Procrustes, and yield high similarity scores even when task relevant variables cannot be linearly decoded from the synthetic datasets. Synthetic datasets optimized to maximize similarity scores initially learn the first principal component of the target dataset, but angular Procrustes captures higher variance dimensions much earlier than methods like linear regression and CKA. We show in both theory and simulations how these scores change when different principal components are perturbed. And finally, we jointly optimize multiple similarity scores to find their allowed ranges, and show that a high angular Procrustes similarity, for example, implies a high CKA score, but not the converse. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 16 pages, 6 figures

An efficient pipeline to compute patient-specific cerebral aneurysm wall tension

Authors: Mostafa Jamshidian, Benjamin Zwick, Arosha S Dissanayake, Adam Wittek, Timothy J Phillips, Stephen Honeybul, Graeme J Hankey, Karol Miller

Abstract: Cerebral aneurysm rupture, leading to subarachnoid hemorrhage with a high mortality rate, disproportionately affects younger populations, resulting in a significant loss of productive life years. A significant proportion of these deaths is due to aneurysmal re-bleeding within the first three days following the initial bleed, prior to treatment. While early aneurysm treatment is recommended, there… ▽ More Cerebral aneurysm rupture, leading to subarachnoid hemorrhage with a high mortality rate, disproportionately affects younger populations, resulting in a significant loss of productive life years. A significant proportion of these deaths is due to aneurysmal re-bleeding within the first three days following the initial bleed, prior to treatment. While early aneurysm treatment is recommended, there is no consensus on the ideal timing, and emergency treatment offers only an incremental benefit at a significant cost. Although various multivariable prediction models have been proposed to provide personalized risk assessments, no validated patient-specific predictor is available to rationalize emergency treatment. Furthermore, no model has yet incorporated emerging computational biomechanics-based biomarkers such as wall tension. In this paper, we proposed and validated an efficient semi-automatic pipeline to compute patient-specific cerebral aneurysm wall tension as a potential biomarker for the likelihood of re-bleeding. Our pipeline uses the patient's computed tomography angiography (CTA) image obtained at the time of subarachnoid hemorrhage diagnosis to create a patient-specific biomechanical model of the cerebral aneurysm using the finite element method. A distinctive feature of our approach is the straightforward model creation and wall tension computation using shell finite elements, without requiring patient-specific material properties or aneurysm wall thickness. Our non-invasive, patient-specific method for cerebral aneurysm wall tension can potentially provide individualized risk prediction and enhance clinical decision-making. △ Less

Submitted 6 July, 2024; originally announced July 2024.

Abdominal aortic aneurysm wall stress: A 7-line code in MATLAB and a one-click software application

Authors: Mostafa Jamshidian, Saeideh Sekhavat, Adam Wittek, Karol Miller

Abstract: An abdominal aortic aneurysm (AAA) is a life-threatening condition characterized by the irreversible dilation of the lower aorta, usually detected incidentally during imaging for other health issues. Current clinical practice for managing AAA relies on a one-size-fits-all approach, based on the aneurysm's maximum diameter and growth rate, which can lead to underestimation or overestimation of AAA… ▽ More An abdominal aortic aneurysm (AAA) is a life-threatening condition characterized by the irreversible dilation of the lower aorta, usually detected incidentally during imaging for other health issues. Current clinical practice for managing AAA relies on a one-size-fits-all approach, based on the aneurysm's maximum diameter and growth rate, which can lead to underestimation or overestimation of AAA rupture risk. Patient-specific AAA wall stress, computed using biomechanical models derived from medical images without needing patient-specific material properties, has been widely investigated for developing individualized AAA rupture risk predictors. Therefore, AAA wall stress, determined reliably and quickly, has the potential to enhance patient-specific treatment plans. This paper presents a 7-line code, written in MATLAB using the Partial Differential Equation Toolbox, for AAA wall stress computations via finite element analysis. The code takes AAA wall geometry as input and outputs stress components over the AAA wall domain. Additionally, we present a one-click standalone software application for AAA wall stress computation, developed based on our 7-line code using MATLAB Compiler. After verification, we used our code to compute AAA wall stress in ten patients. Our analysis indicated that the 99th percentile of maximum principal stress across all patients ranged from 0.320 MPa to 0.522 MPa, with an average of 0.401 MPa and a standard deviation of 0.056 MPa. Moreover, for every case, the MATLAB simulation time was less than a minute on a laptop workstation. △ Less

Submitted 6 July, 2024; originally announced July 2024.

Enhance the Image: Super Resolution using Artificial Intelligence in MRI

Authors: Ziyu Li, Zihan Li, Haoxiang Li, Qiuyun Fan, Karla L. Miller, Wenchuan Wu, Akshay S. Chaudhari, Qiyuan Tian

Abstract: This chapter provides an overview of deep learning techniques for improving the spatial resolution of MRI, ranging from convolutional neural networks, generative adversarial networks, to more advanced models including transformers, diffusion models, and implicit neural representations. Our exploration extends beyond the methodologies to scrutinize the impact of super-resolved images on clinical an… ▽ More This chapter provides an overview of deep learning techniques for improving the spatial resolution of MRI, ranging from convolutional neural networks, generative adversarial networks, to more advanced models including transformers, diffusion models, and implicit neural representations. Our exploration extends beyond the methodologies to scrutinize the impact of super-resolved images on clinical and neuroscientific assessments. We also cover various practical topics such as network architectures, image evaluation metrics, network loss functions, and training data specifics, including downsampling methods for simulating low-resolution images and dataset selection. Finally, we discuss existing challenges and potential future directions regarding the feasibility and reliability of deep learning-based MRI super-resolution, with the aim to facilitate its wider adoption to benefit various clinical and neuroscientific applications. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: A book chapter in Machine Learning in MRI: From methods to clinical translation. Copyright may be transferred without notice, after which this version may no longer be accessible

Demonstration of neutron identification in neutrino interactions in the MicroBooNE liquid argon time projection chamber

Authors: MicroBooNE collaboration, P. Abratenko, O. Alterkait, D. Andrade Aldana, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, A. Barnard, G. Barr, D. Barrow, J. Barrow, V. Basque, J. Bateman, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bhattacharya, M. Bishai, A. Blake, B. Bogart, T. Bolton, J. Y. Book , et al. (165 additional authors not shown)

Abstract: A significant challenge in measurements of neutrino oscillations is reconstructing the incoming neutrino energies. While modern fully-active tracking calorimeters such as liquid argon time projection chambers in principle allow the measurement of all final state particles above some detection threshold, undetected neutrons remain a considerable source of missing energy with little to no data const… ▽ More A significant challenge in measurements of neutrino oscillations is reconstructing the incoming neutrino energies. While modern fully-active tracking calorimeters such as liquid argon time projection chambers in principle allow the measurement of all final state particles above some detection threshold, undetected neutrons remain a considerable source of missing energy with little to no data constraining their production rates and kinematics. We present the first demonstration of tagging neutrino-induced neutrons in liquid argon time projection chambers using secondary protons emitted from neutron-argon interactions in the MicroBooNE detector. We describe the method developed to identify neutrino-induced neutrons and demonstrate its performance using neutrons produced in muon-neutrino charged current interactions. The method is validated using a small subset of MicroBooNE's total dataset. The selection yields a sample with $60\%$ of selected tracks corresponding to neutron-induced secondary protons. △ Less

Submitted 15 June, 2024; originally announced June 2024.

Report number: FERMILAB-PUB-24-0301

Improving neutrino energy estimation of charged-current interaction events with recurrent neural networks in MicroBooNE

Authors: MicroBooNE collaboration, P. Abratenko, O. Alterkait, D. Andrade Aldana, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, A. Barnard, G. Barr, D. Barrow, J. Barrow, V. Basque, J. Bateman, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bhattacharya, M. Bishai, A. Blake, B. Bogart, T. Bolton, J. Y. Book , et al. (164 additional authors not shown)

Abstract: We present a deep learning-based method for estimating the neutrino energy of charged-current neutrino-argon interactions. We employ a recurrent neural network (RNN) architecture for neutrino energy estimation in the MicroBooNE experiment, utilizing liquid argon time projection chamber (LArTPC) detector technology. Traditional energy estimation approaches in LArTPCs, which largely rely on reconstr… ▽ More We present a deep learning-based method for estimating the neutrino energy of charged-current neutrino-argon interactions. We employ a recurrent neural network (RNN) architecture for neutrino energy estimation in the MicroBooNE experiment, utilizing liquid argon time projection chamber (LArTPC) detector technology. Traditional energy estimation approaches in LArTPCs, which largely rely on reconstructing and summing visible energies, often experience sizable biases and resolution smearing because of the complex nature of neutrino interactions and the detector response. The estimation of neutrino energy can be improved after considering the kinematics information of reconstructed final-state particles. Utilizing kinematic information of reconstructed particles, the deep learning-based approach shows improved resolution and reduced bias for the muon neutrino Monte Carlo simulation sample compared to the traditional approach. In order to address the common concern about the effectiveness of this method on experimental data, the RNN-based energy estimator is further examined and validated with dedicated data-simulation consistency tests using MicroBooNE data. We also assess its potential impact on a neutrino oscillation study after accounting for all statistical and systematic uncertainties and show that it enhances physics sensitivity. This method has good potential to improve the performance of other physics analyses. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Report number: FERMILAB-PUB-24-0287

FlowMM: Generating Materials with Riemannian Flow Matching

Authors: Benjamin Kurt Miller, Ricky T. Q. Chen, Anuroop Sriram, Brandon M Wood

Abstract: Crystalline materials are a fundamental component in next-generation technologies, yet modeling their distribution presents unique computational challenges. Of the plausible arrangements of atoms in a periodic lattice only a vanishingly small percentage are thermodynamically stable, which is a key indicator of the materials that can be experimentally realized. Two fundamental tasks in this area ar… ▽ More Crystalline materials are a fundamental component in next-generation technologies, yet modeling their distribution presents unique computational challenges. Of the plausible arrangements of atoms in a periodic lattice only a vanishingly small percentage are thermodynamically stable, which is a key indicator of the materials that can be experimentally realized. Two fundamental tasks in this area are to (a) predict the stable crystal structure of a known composition of elements and (b) propose novel compositions along with their stable structures. We present FlowMM, a pair of generative models that achieve state-of-the-art performance on both tasks while being more efficient and more flexible than competing methods. We generalize Riemannian Flow Matching to suit the symmetries inherent to crystals: translation, rotation, permutation, and periodic boundary conditions. Our framework enables the freedom to choose the flow base distributions, drastically simplifying the problem of learning crystal structures compared with diffusion models. In addition to standard benchmarks, we validate FlowMM's generated structures with quantum chemistry calculations, demonstrating that it is about 3x more efficient, in terms of integration steps, at finding stable materials compared to previous open methods. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: https://github.com/facebookresearch/flowmm

Journal ref: ICML 2024

Dephasingless two-color terahertz generation

Authors: Tanner T. Simpson, Jeremy J. Pigeon, Kyle G. Miller, Dillon Ramsey, Dustin H. Froula, John P. Palastro

Abstract: A laser pulse composed of a fundamental and an appropriately phased second harmonic can drive a time-dependent current of photoionized electrons that generates broadband THz radiation. Over the propagation distances relevant to many experiments, dispersion causes the relative phase between the harmonics to evolve. This "dephasing" slows the accumulation of THz energy and results in a multi-cycle T… ▽ More A laser pulse composed of a fundamental and an appropriately phased second harmonic can drive a time-dependent current of photoionized electrons that generates broadband THz radiation. Over the propagation distances relevant to many experiments, dispersion causes the relative phase between the harmonics to evolve. This "dephasing" slows the accumulation of THz energy and results in a multi-cycle THz pulse with significant angular dispersion. Here, we introduce a novel optical configuration that compensates the relative phase evolution, allowing for the formation of a half-cycle THz pulse with almost no angular dispersion. The configuration uses the spherical aberration of an axilens to map a prescribed radial phase variation in the near field to a desired longitudinal phase variation in the far field. Simulations that combine this configuration with an ultrashort flying focus demonstrate the formation of a half-cycle THz pulse with a controlled emission angle and 1/4 the angular divergence of the multi-cycle pulse created by a conventional optical configuration. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Galois descent of splendid Rickard equivalences for blocks of $p$-nilpotent groups

Authors: Sam K. Miller

Abstract: We strengthen the results of Boltje and Yilmaz regarding the Galois descent of equivalences of blocks of $p$-nilpotent groups and a result of Kessar and Linckelmann regarding Galois descent of splendid Rickard equivalences for blocks with compatible Galois stabilizers. A more general descent criteria for chain complexes is proven along the way, which requires the adaptation of a theorem of Reiner… ▽ More We strengthen the results of Boltje and Yilmaz regarding the Galois descent of equivalences of blocks of $p$-nilpotent groups and a result of Kessar and Linckelmann regarding Galois descent of splendid Rickard equivalences for blocks with compatible Galois stabilizers. A more general descent criteria for chain complexes is proven along the way, which requires the adaptation of a theorem of Reiner for chain complexes. This verifies Kessar and Linckelmann's refinement of Broué's abelian defect group conjecture for blocks of $p$-nilpotent groups with abelian Sylow $p$-subgroup. △ Less

Submitted 25 May, 2024; originally announced May 2024.

Comments: 7 pages, 1 figure, comments welcome!

MSC Class: 20C20; 20J05; 20C05; 20F19; 18G35

Kinematics of Abdominal Aortic Aneurysms

Authors: Mostafa Jamshidian, Adam Wittek, Saeideh Sekhavat, Karol Miller

Abstract: A search in Scopus within "Article title, Abstract, Keywords" unveils 2,444 documents focused on the biomechanics of Abdominal Aortic Aneurysm (AAA), mostly on AAA wall stress. Only 24 documents investigated AAA kinematics, an important topic that could potentially offer insights into the biomechanics of AAA. In this paper, we present an image-based approach for patient-specific, in vivo, and non-… ▽ More A search in Scopus within "Article title, Abstract, Keywords" unveils 2,444 documents focused on the biomechanics of Abdominal Aortic Aneurysm (AAA), mostly on AAA wall stress. Only 24 documents investigated AAA kinematics, an important topic that could potentially offer insights into the biomechanics of AAA. In this paper, we present an image-based approach for patient-specific, in vivo, and non-invasive AAA kinematic analysis using patient's time-resolved 3D computed tomography angiography (4D CTA) images. Our approach relies on regularized deformable image registration for estimating wall displacement, estimation of the local wall strain as the ratio of its normal displacement to its local radius of curvature, and local surface fitting with non-deterministic outlier detection for estimating the wall radius of curvature. We verified our approach against synthetic ground truth image data created by warping a 3D CTA image of AAA using a realistic displacement field obtained from a finite element biomechanical model. We applied our approach to assess AAA wall displacements and strains in ten patients. Our kinematic analysis results indicated that the 99th percentile of circumferential wall strain, among all patients, ranged from 3.16% to 7.31%, with an average of 5.36% and a standard deviation of 1.28%. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Human-interpretable clustering of short-text using large language models

Authors: Justin K. Miller, Tristram J. Alexander

Abstract: Large language models have seen extraordinary growth in popularity due to their human-like content generation capabilities. We show that these models can also be used to successfully cluster human-generated content, with success defined through the measures of distinctiveness and interpretability. This success is validated by both human reviewers and ChatGPT, providing an automated means to close… ▽ More Large language models have seen extraordinary growth in popularity due to their human-like content generation capabilities. We show that these models can also be used to successfully cluster human-generated content, with success defined through the measures of distinctiveness and interpretability. This success is validated by both human reviewers and ChatGPT, providing an automated means to close the 'validation gap' that has challenged short-text clustering. Comparing the machine and human approaches we identify the biases inherent in each, and question the reliance on human-coding as the 'gold standard'. We apply our methodology to Twitter bios and find characteristic ways humans describe themselves, agreeing well with prior specialist work, but with interesting differences characteristic of the medium used to express identity. △ Less

Submitted 12 May, 2024; originally announced May 2024.

Comments: Main text: 18 pages, 8 figures. Supplementary: 21 pages, 15 figures, 3 tables

ACM Class: I.2.7

Inhomogeneous wave kinetic equation and its hierarchy in polynomially weighted $L^\infty$ spaces

Authors: Ioakeim Ampatzoglou, Joseph K. Miller, Nataša Pavlović, Maja Tasković

Abstract: Inspired by ideas stemming from the analysis of the Boltzmann equation, in this paper we expand well-posedness theory of the spatially inhomogeneous 4-wave kinetic equation, and also analyze an infinite hierarchy of PDE associated with this nonlinear equation. More precisely, we show global in time well-posedness of the spatially inhomogeneous 4-wave kinetic equation for polynomially decaying init… ▽ More Inspired by ideas stemming from the analysis of the Boltzmann equation, in this paper we expand well-posedness theory of the spatially inhomogeneous 4-wave kinetic equation, and also analyze an infinite hierarchy of PDE associated with this nonlinear equation. More precisely, we show global in time well-posedness of the spatially inhomogeneous 4-wave kinetic equation for polynomially decaying initial data. For the associated infinite hierarchy, we construct global in time solutions using the solutions of the wave kinetic equation and the Hewitt-Savage theorem. Uniqueness of these solutions is proved by using a combinatorial board game argument tailored to this context, which allows us to control the factorial growth of the Dyson series. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 40 pages

First double-differential cross section measurement of neutral-current $π^0$ production in neutrino-argon scattering in the MicroBooNE detector

Authors: MicroBooNE collaboration, P. Abratenko, O. Alterkait, D. Andrade Aldana, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, A. Barnard, G. Barr, D. Barrow, J. Barrow, V. Basque, J. Bateman, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bhattacharya, M. Bishai, A. Blake, B. Bogart, T. Bolton, J. Y. Book , et al. (166 additional authors not shown)

Abstract: We report the first double-differential cross section measurement of neutral-current neutral pion (NC$π^0$) production in neutrino-argon scattering, as well as single-differential measurements of the same channel in terms of final states with and without protons. The kinematic variables of interest for these measurements are the $π^0$ momentum and the $π^0$ scattering angle with respect to the neu… ▽ More We report the first double-differential cross section measurement of neutral-current neutral pion (NC$π^0$) production in neutrino-argon scattering, as well as single-differential measurements of the same channel in terms of final states with and without protons. The kinematic variables of interest for these measurements are the $π^0$ momentum and the $π^0$ scattering angle with respect to the neutrino beam. A total of 4971 candidate NC$π^0$ events fully-contained within the MicroBooNE detector are selected using data collected at a mean neutrino energy of $\sim 0.8$ GeV from $6.4\times10^{20}$ protons on target from the Booster Neutrino Beam at the Fermi National Accelerator Laboratory. After extensive data-driven model validation to ensure unbiased unfolding, the Wiener-SVD method is used to extract nominal flux-averaged cross sections. The results are compared to predictions from commonly used neutrino event generators, which tend to overpredict the measured NC$π^0$ cross section, especially in the 0.2-0.5 GeV/c $π^0$ momentum range, at forward scattering angles, and when at least one proton is present in the final state. These measurements show sensitivity to a variety of features that complicate the description of NC$π^0$ production including the form factors describing the elementary neutrino interaction and the final state interactions of the outgoing particles in the residual argon nucleus. This data will help improve the modeling of NC$π^0$ production, which represents a major background in measurements of charge-parity violation in the neutrino sector and in searches for new physics beyond the Standard Model. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Report number: FERMILAB-PUB-24-0125

Measurement of the differential cross section for neutral pion production in charged-current muon neutrino interactions on argon with the MicroBooNE detector

Authors: MicroBooNE collaboration, P. Abratenko, O. Alterkait, D. Andrade Aldana, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, G. Barr, D. Barrow, J. Barrow, V. Basque, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bhattacharya, M. Bishai, A. Blake, B. Bogart, T. Bolton, J. Y. Book, M. B. Brunetti, L. Camilleri , et al. (163 additional authors not shown)

Abstract: We present a measurement of neutral pion production in charged-current interactions using data recorded with the MicroBooNE detector exposed to Fermilab's booster neutrino beam. The signal comprises one muon, one neutral pion, any number of nucleons, and no charged pions. Studying neutral pion production in the MicroBooNE detector provides an opportunity to better understand neutrino-argon interac… ▽ More We present a measurement of neutral pion production in charged-current interactions using data recorded with the MicroBooNE detector exposed to Fermilab's booster neutrino beam. The signal comprises one muon, one neutral pion, any number of nucleons, and no charged pions. Studying neutral pion production in the MicroBooNE detector provides an opportunity to better understand neutrino-argon interactions, and is crucial for future accelerator-based neutrino oscillation experiments. Using a dataset corresponding to $6.86 \times 10^{20}$ protons on target, we present single-differential cross sections in muon and neutral pion momenta, scattering angles with respect to the beam for the outgoing muon and neutral pion, as well as the opening angle between the muon and neutral pion. Data extracted cross sections are compared to generator predictions. We report good agreement between the data and the models for scattering angles, except for an over-prediction by generators at muon forward angles. Similarly, the agreement between data and the models as a function of momentum is good, except for an underprediction by generators in the medium momentum ranges, $200-400$ MeV for muons and $100-200$ MeV for pions. △ Less

Submitted 6 May, 2024; v1 submitted 15 April, 2024; originally announced April 2024.

Report number: FERMILAB-PUB-24-0142-CSAID-PPD

Measurement of double-differential cross sections for mesonless charged-current muon neutrino interactions on argon with final-state protons using the MicroBooNE detector

Authors: MicroBooNE collaboration, P. Abratenko, O. Alterkait, D. Andrade Aldana, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, G. Barr, D. Barrow, J. Barrow, V. Basque, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bhattacharya, M. Bishai, A. Blake, B. Bogart, T. Bolton, J. Y. Book, M. B. Brunetti, L. Camilleri , et al. (163 additional authors not shown)

Abstract: Charged-current neutrino interactions with final states containing zero mesons and at least one proton are of high interest for current and future accelerator-based neutrino oscillation experiments. Using the Booster Neutrino Beam and the MicroBooNE detector at Fermi National Accelerator Laboratory, we have obtained the first double-differential cross section measurements of this channel for muon… ▽ More Charged-current neutrino interactions with final states containing zero mesons and at least one proton are of high interest for current and future accelerator-based neutrino oscillation experiments. Using the Booster Neutrino Beam and the MicroBooNE detector at Fermi National Accelerator Laboratory, we have obtained the first double-differential cross section measurements of this channel for muon neutrino scattering on an argon target with a proton momentum threshold of 0.25 GeV/c. We also report a flux-averaged total cross section of $σ= (11.8 \pm 1.2) \times 10^{-38}$ cm$^2$ / Ar and several single-differential measurements which extend and improve upon previous results. Statistical and systematic uncertainties are quantified with a full treatment of correlations across 359 kinematic bins, including correlations between distributions describing different observables. The resulting data set provides the most detailed information obtained to date for testing models of mesonless neutrino-argon scattering. △ Less

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

Comments: 83 pages, 67 figures (including supplemental material). For v2, added oversized files in extended data release

Report number: FERMILAB-PUB-24-0120-AD-CSAID-LBNF-PPD-TD

Active Learning of Dynamics Using Prior Domain Knowledge in the Sampling Process

Authors: Kevin S. Miller, Adam J. Thorpe, Ufuk Topcu

Abstract: We present an active learning algorithm for learning dynamics that leverages side information by explicitly incorporating prior domain knowledge into the sampling process. Our proposed algorithm guides the exploration toward regions that demonstrate high empirical discrepancy between the observed data and an imperfect prior model of the dynamics derived from side information. Through numerical exp… ▽ More We present an active learning algorithm for learning dynamics that leverages side information by explicitly incorporating prior domain knowledge into the sampling process. Our proposed algorithm guides the exploration toward regions that demonstrate high empirical discrepancy between the observed data and an imperfect prior model of the dynamics derived from side information. Through numerical experiments, we demonstrate that this strategy explores regions of high discrepancy and accelerates learning while simultaneously reducing model uncertainty. We rigorously prove that our active learning algorithm yields a consistent estimate of the underlying dynamics by providing an explicit rate of convergence for the maximum predictive variance. We demonstrate the efficacy of our approach on an under-actuated pendulum system and on the half-cheetah MuJoCo environment. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Towards Full Automation of Geometry Extraction for Biomechanical Analysis of Abdominal Aortic Aneurysm; Neural Network-Based versus Classical Methodologies

Authors: Farah Alkhatib, Mostafa Jamshidian, Donatien Le Liepvre, Florian Bernard, Ludovic Minvielle, Adam Wittek, Karol Miller

Abstract: In this study we investigated the impact of image segmentation methods on the results of stress computation in the wall of abdominal aortic aneurysms (AAAs). We compared wall stress distributions and magnitudes calculated from geometry models obtained from classical semi-automated segmentation versus automated neural network-based segmentation. Ten different AAA contrast-enhanced computed tomograp… ▽ More In this study we investigated the impact of image segmentation methods on the results of stress computation in the wall of abdominal aortic aneurysms (AAAs). We compared wall stress distributions and magnitudes calculated from geometry models obtained from classical semi-automated segmentation versus automated neural network-based segmentation. Ten different AAA contrast-enhanced computed tomography (CT) images were semi-automatically segmented by an analyst, taking, depending on the quality of an image, between 15 and 40 minutes of human effort per patient. The same images were automatically segmented using PRAEVAorta 2, commercial software by NUREA (https://www.nurea-soft.com/), developed based on artificial intelligence (AI) algorithms, requiring only 1-2 minutes of computer time per patient. Aneurysm wall stress calculations performed using the BioPARR software (https://bioparr.mech.uwa.edu.au/) revealed that, compared to the classical semi-automated segmentation, the automatic neural network-based segmentation leads to equivalent stress distributions, and slightly higher peak and 99th percentile maximum principal stress values. This difference is due to consistently larger lumen surface areas in automatically segmented models as compared to classical semi-automated segmentations, resulting in greater total pressure load on the wall. Our findings are a steppingstone toward a fully automated pipeline for biomechanical analysis of AAAs, starting with CT scans and concluding with wall stress assessment, while at the same time highlighting the critical importance of the repeatable and accurate segmentation of the lumen, the difficult problem often underestimated by the literature. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 32 pages, 9 figures

On endotrivial complexes and the generalized Dade group

Authors: Sam K. Miller

Abstract: Let $p$ be a prime, $G$ a $p$-group, and $k$ a field of characteristic $p$. We give a complete description of the group $\mathcal{E}_k(G)$ of endotrivial complexes, identifying $\mathcal{E}_k(G)$ with the additive subgroup of integer-valued class functions satisfying the Borel-Smith conditions. This is done by constructing a short exact sequence isomorphic to a short exact sequence of rational… ▽ More Let $p$ be a prime, $G$ a $p$-group, and $k$ a field of characteristic $p$. We give a complete description of the group $\mathcal{E}_k(G)$ of endotrivial complexes, identifying $\mathcal{E}_k(G)$ with the additive subgroup of integer-valued class functions satisfying the Borel-Smith conditions. This is done by constructing a short exact sequence isomorphic to a short exact sequence of rational $p$-biset functors constructed by Bouc and Yalçin, endowing $\mathcal{E}_k$ with rational $p$-biset structure. We also determine all the possible h-marks which arise in $\mathcal{E}_k(G)$ when $G$ is a finite group via the induced restriction map to a Sylow $p$-subgroup. As a consequence, we prove that every $p$-permutation autoequivalence of a $p$-group extends to a splendid Rickard autoequivalence. Additionally, we give a positive answer to a question of Gelvin and Yalçin, showing the kernel of the Bouc homomorphism for an arbitrary finite group $G$ is described by integer-valued class functions on $p$-subgroups satisfying the oriented Artin-Borel-Smith conditions. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 25 pages, comments welcome!

MSC Class: 20J05; 19A22; 20C05; 20C20

Consumers' Perceived Privacy Violations in Online Advertising

Authors: Kinshuk Jerath, Klaus M. Miller

Abstract: In response to privacy concerns about collecting and using personal data, the online advertising industry has been developing privacy-enhancing technologies (PETs), e.g., under Google's Privacy Sandbox initiative. In this research, we use the dual-privacy framework, which postulates that consumers have intrinsic and instrumental preferences for privacy, to understand consumers' perceived privacy v… ▽ More In response to privacy concerns about collecting and using personal data, the online advertising industry has been developing privacy-enhancing technologies (PETs), e.g., under Google's Privacy Sandbox initiative. In this research, we use the dual-privacy framework, which postulates that consumers have intrinsic and instrumental preferences for privacy, to understand consumers' perceived privacy violations (PPVs) for current and proposed online advertising practices. The key idea is that different practices differ in whether individual data leaves the consumer's machine or not and in how they track and target consumers; these affect, respectively, the intrinsic and instrumental components of privacy preferences differently, leading to different PPVs for different practices. We conducted online studies focused on consumers in the United States to elicit PPVs for various advertising practices. Our findings confirm the intuition that tracking and targeting consumers under the industry status quo of behavioral targeting leads to high PPV. New technologies or proposals that ensure that data are kept on the consumer's machine lower PPV relative to behavioral targeting but, importantly, this decrease is small. Furthermore, group-level targeting does not differ significantly from individual-level targeting in reducing PPV. Under contextual targeting, where there is no tracking, PPV is significantly reduced. Interestingly, with respect to PPV, consumers are indifferent between seeing untargeted ads and no ads when they are not being tracked. We find that consumer perceptions of privacy violations under different tracking and targeting practices may differ from what technical definitions suggest. Therefore, rather than relying solely on technical perspectives, a consumer-centric approach to privacy is needed, based on, for instance, the dual-privacy framework. △ Less

Submitted 28 May, 2024; v1 submitted 6 March, 2024; originally announced March 2024.

Paying for Privacy: Pay-or-Tracking Walls

Authors: Timo Mueller-Tribbensee, Klaus M. Miller, Bernd Skiera

Abstract: Prestigious news publishers, and more recently, Meta, have begun to request that users pay for privacy. Specifically, users receive a notification banner, referred to as a pay-or-tracking wall, that requires them to (i) pay money to avoid being tracked or (ii) consent to being tracked. These walls have invited concerns that privacy might become a luxury. However, little is known about pay-or-track… ▽ More Prestigious news publishers, and more recently, Meta, have begun to request that users pay for privacy. Specifically, users receive a notification banner, referred to as a pay-or-tracking wall, that requires them to (i) pay money to avoid being tracked or (ii) consent to being tracked. These walls have invited concerns that privacy might become a luxury. However, little is known about pay-or-tracking walls, which prevents a meaningful discussion about their appropriateness. This paper conducts several empirical studies and finds that top EU publishers use pay-or-tracking walls. Their implementations involve various approaches, including bundling the pay option with advertising-free access or additional content. The price for not being tracked exceeds the advertising revenue that publishers generate from a user who consents to being tracked. Notably, publishers' traffic does not decline when implementing a pay-or-tracking wall and most users consent to being tracked; only a few users pay. In short, pay-or-tracking walls seem to provide the means for expanding the practice of tracking. Publishers profit from pay-or-tracking walls and may observe a revenue increase of 16.4% due to tracking more users than under a cookie consent banner. △ Less

Submitted 6 March, 2024; originally announced March 2024.

First simultaneous measurement of differential muon-neutrino charged-current cross sections on argon for final states with and without protons using MicroBooNE data

Authors: MicroBooNE collaboration, P. Abratenko, O. Alterkait, D. Andrade Aldana, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, G. Barr, D. Barrow, J. Barrow, V. Basque, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bhattacharya, M. Bishai, A. Blake, B. Bogart, T. Bolton, J. Y. Book, M. B. Brunetti, L. Camilleri , et al. (163 additional authors not shown)

Abstract: We report the first double-differential neutrino-argon cross section measurement made simultaneously for final states with and without protons for the inclusive muon neutrino charged-current interaction channel. The proton kinematics of this channel are further explored with a differential cross section measurement as a function of the leading proton's kinetic energy that extends across the detect… ▽ More We report the first double-differential neutrino-argon cross section measurement made simultaneously for final states with and without protons for the inclusive muon neutrino charged-current interaction channel. The proton kinematics of this channel are further explored with a differential cross section measurement as a function of the leading proton's kinetic energy that extends across the detection threshold. These measurements utilize data collected using the MicroBooNE detector from 6.4$\times10^{20}$ protons on target from the Fermilab Booster Neutrino Beam with a mean neutrino energy of $\sim$0.8 GeV. Extensive data-driven model validation utilizing the conditional constraint formalism is employed. This motivates enlarging the uncertainties with an empirical reweighting approach to minimize the possibility of extracting biased cross section results. The extracted nominal flux-averaged cross sections are compared to widely used event generator predictions revealing severe mismodeling of final states without protons for muon neutrino charged-current interactions, possibly from insufficient treatment of final state interactions. These measurements provide a wealth of new information useful for improving event generators which will enhance the sensitivity of precision measurements in neutrino experiments. △ Less

Submitted 27 July, 2024; v1 submitted 29 February, 2024; originally announced February 2024.

Report number: FERMILAB-PUB-24-0045

Journal ref: Phys. Rev. Lett. 133, 041801 (2024)

Inclusive cross section measurements in final states with and without protons for charged-current $ν_μ$-Ar scattering in MicroBooNE

Authors: MicroBooNE collaboration, P. Abratenko, O. Alterkait, D. Andrade Aldana, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, G. Barr, D. Barrow, J. Barrow, V. Basque, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bhattacharya, M. Bishai, A. Blake, B. Bogart, T. Bolton, J. Y. Book, M. B. Brunetti, L. Camilleri , et al. (164 additional authors not shown)

Abstract: A detailed understanding of inclusive muon neutrino charged-current interactions on argon is crucial to the study of neutrino oscillations in current and future experiments using liquid argon time projection chambers. To that end, we report a comprehensive set of differential cross section measurements for this channel that simultaneously probe the leptonic and hadronic systems by dividing the cha… ▽ More A detailed understanding of inclusive muon neutrino charged-current interactions on argon is crucial to the study of neutrino oscillations in current and future experiments using liquid argon time projection chambers. To that end, we report a comprehensive set of differential cross section measurements for this channel that simultaneously probe the leptonic and hadronic systems by dividing the channel into final states with and without protons. Measurements of the proton kinematics and proton multiplicity of the final state are also presented. For these measurements, we utilize data collected with the MicroBooNE detector from 6.4$\times10^{20}$ protons on target from the Fermilab Booster Neutrino Beam at a mean neutrino energy of approximately 0.8 GeV. We present in detail the cross section extraction procedure, including the unfolding, and model validation that uses data to model comparisons and the conditional constraint formalism to detect mismodeling that may introduce biases to extracted cross sections that are larger than their uncertainties. The validation exposes insufficiencies in the overall model, motivating the inclusion of an additional data-driven reweighting systematic to ensure the accuracy of the unfolding. The extracted results are compared to a number of event generators and their performance is discussed with a focus on the regions of phase-space that indicate the greatest need for modeling improvements. △ Less

Submitted 27 July, 2024; v1 submitted 29 February, 2024; originally announced February 2024.

Report number: FERMILAB-PUB-24-0044

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

Open Meshed Anatomy: Towards a comprehensive finite element hexahedral mesh derived from open atlases

Authors: Andy Trung Huynh, Benjamin Zwick, Michael Halle, Adam Wittek, Karol Miller

Abstract: Computational simulations using methods such as the finite element (FE) method rely on high-quality meshes for achieving accurate results. This study introduces a method for creating a high-quality hexahedral mesh using the Open Anatomy Project's brain atlas. Our atlas-based FE hexahedral mesh of the brain mitigates potential inaccuracies and uncertainties due to segmentation - a process that ofte… ▽ More Computational simulations using methods such as the finite element (FE) method rely on high-quality meshes for achieving accurate results. This study introduces a method for creating a high-quality hexahedral mesh using the Open Anatomy Project's brain atlas. Our atlas-based FE hexahedral mesh of the brain mitigates potential inaccuracies and uncertainties due to segmentation - a process that often requires input of an inexperienced analyst. It accomplishes this by leveraging existing segmentation from the atlas. We further extend the mesh's usability by forming a two-way correspondence between the atlas and mesh. This feature facilitates property assignment for computational simulations and enhances result analysis within an anatomical context. We demonstrate the application of the mesh by solving the electroencephalography (EEG) forward problem. Our method simplifies the mesh creation process, reducing time and effort, and provides a more comprehensive and contextually enriched visualisation of simulation outcomes. △ Less

Submitted 22 February, 2024; originally announced February 2024.

Relatively endotrivial complexes

Authors: Sam K. Miller

Abstract: Let $G$ be a finite group and $k$ be a field of characteristic $p > 0$. In prior work, we studied endotrivial complexes, the invertible objects of the bounded homotopy category $K^b({}_{kG}\mathbf{triv})$ of $p$-permutation $kG$-modules. Using the notion of projectivity relative to a $kG$-module, we expand on this study by defining notions of "relatively" endotrivial chain complexes, analogous to… ▽ More Let $G$ be a finite group and $k$ be a field of characteristic $p > 0$. In prior work, we studied endotrivial complexes, the invertible objects of the bounded homotopy category $K^b({}_{kG}\mathbf{triv})$ of $p$-permutation $kG$-modules. Using the notion of projectivity relative to a $kG$-module, we expand on this study by defining notions of "relatively" endotrivial chain complexes, analogous to Lassueur's construction of relatively endotrivial $kG$-modules. We obtain equivalent characterizations of relative endotriviality and find corresponding local homological data which almost completely determine the isomorphism class of a relatively endotrivial complex. We show this local data must partially satisfy the Borel-Smith conditions, and consider the behavior of restriction to subgroups containing Sylow $p$-subgroups $S$ of $G$. △ Less

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

Comments: v2: Removed numerous minor results and added a new section on Borel-Smith functions. 31 pages

MSC Class: 20J05; 20C05; 20C20

On the global in time existence and uniqueness of solutions to the Boltzmann hierarchy

Authors: Ioakeim Ampatzoglou, Joseph K. Miller, Nataša Pavlović, Maja Tasković

Abstract: In this paper we establish the global in time existence and uniqueness of solutions to the Boltzmann hierarchy, a hierarchy of equations instrumental for the rigorous derivation of the Boltzmann equation from many particles. Inspired by available $L^{\infty}$-based a-priori estimate for solutions to the Boltzmann equation, we develop the polynomially weighted $L^\infty$ a-priori bounds for solutio… ▽ More In this paper we establish the global in time existence and uniqueness of solutions to the Boltzmann hierarchy, a hierarchy of equations instrumental for the rigorous derivation of the Boltzmann equation from many particles. Inspired by available $L^{\infty}$-based a-priori estimate for solutions to the Boltzmann equation, we develop the polynomially weighted $L^\infty$ a-priori bounds for solutions to the Boltzmann hierarchy and handle the factorial growth of the number of terms in the Dyson's series by reorganizing the sum through a combinatorial technique known as the Klainerman-Machedon board game argument. This paper is the first work that exploits such a combinatorial technique in conjunction with an $L^{\infty}$-based estimate to prove uniqueness of the mild solutions to the Boltzmann hierarchy. Our proof of existence of global in time mild solutions to the Boltzmann hierarchy for admissible initial data is constructive and it employs known global in time solutions to the Boltzmann equation via a Hewitt-Savage type theorem. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 41 pages

Self-navigated 3D diffusion MRI using an optimized CAIPI sampling and structured low-rank reconstruction

Authors: Ziyu Li, Karla L. Miller, Xi Chen, Mark Chiew, Wenchuan Wu

Abstract: 3D multi-slab acquisitions are an appealing approach for diffusion MRI because they are compatible with the imaging regime delivering optimal SNR efficiency. In conventional 3D multi-slab imaging, shot-to-shot phase variations caused by motion pose challenges due to the use of multi-shot k-space acquisition. Navigator acquisition after each imaging echo is typically employed to correct phase varia… ▽ More 3D multi-slab acquisitions are an appealing approach for diffusion MRI because they are compatible with the imaging regime delivering optimal SNR efficiency. In conventional 3D multi-slab imaging, shot-to-shot phase variations caused by motion pose challenges due to the use of multi-shot k-space acquisition. Navigator acquisition after each imaging echo is typically employed to correct phase variations, which prolongs scan time and increases the specific absorption rate (SAR). The aim of this study is to develop a highly efficient, self-navigated method to correct for phase variations in 3D multi-slab diffusion MRI without explicitly acquiring navigators. The sampling of each shot is carefully designed to intersect with the central kz plane of each slab, and the multi-shot sampling is optimized for self-navigation performance while retaining decent reconstruction quality. The central kz intersections from all shots are jointly used to reconstruct a 2D phase map for each shot using a structured low-rank constrained reconstruction that leverages the redundancy in shot and coil dimensions. The phase maps are used to eliminate the shot-to-shot phase inconsistency in the final 3D multi-shot reconstruction. We demonstrate the method's efficacy using retrospective simulations and prospectively acquired in-vivo experiments at 1.22 mm and 1.09 mm isotropic resolutions. Compared to conventional navigated 3D multi-slab imaging, the proposed self-navigated method achieves comparable image quality while shortening the scan time by 31.7% and improving the SNR efficiency by 15.5%. The proposed method produces comparable quality of DTI and white matter tractography to conventional navigated 3D multi-slab acquisition with a much shorter scan time. △ Less

Submitted 11 January, 2024; originally announced January 2024.

Comments: 10 pages, 11 figures, 2 tables. This work has been submitted to the IEEE for possible publication. Copyright may be transferred without notice, after which this version may no longer be accessible

First search for dark-trident processes using the MicroBooNE detector

Authors: MicroBooNE collaboration, P. Abratenko, O. Alterkait, D. Andrade Aldana, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, G. Barr, D. Barrow, J. Barrow, V. Basque, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bhattacharya, M. Bishai, A. Blake, B. Bogart, T. Bolton, J. Y. Book, M. B. Brunetti, L. Camilleri , et al. (163 additional authors not shown)

Abstract: We present a first search for dark-trident scattering in a neutrino beam using a data set corresponding to $7.2 \times 10^{20}$ protons on target taken with the MicroBooNE detector at Fermilab. Proton interactions in the neutrino target at the Main Injector produce $π^0$ and $η$ mesons, which could decay into dark-matter (DM) particles mediated via a dark photon $A^\prime$. A convolutional neural… ▽ More We present a first search for dark-trident scattering in a neutrino beam using a data set corresponding to $7.2 \times 10^{20}$ protons on target taken with the MicroBooNE detector at Fermilab. Proton interactions in the neutrino target at the Main Injector produce $π^0$ and $η$ mesons, which could decay into dark-matter (DM) particles mediated via a dark photon $A^\prime$. A convolutional neural network is trained to identify interactions of the DM particles in the liquid-argon time projection chamber (LArTPC) exploiting its image-like reconstruction capability. In the absence of a DM signal, we provide limits at the $90\%$ confidence level on the squared kinematic mixing parameter $\varepsilon^2$ as a function of the dark-photon mass in the range $10\le M_{A^\prime}\le 400$ MeV. The limits cover previously unconstrained parameter space for the production of fermion or scalar DM particles $χ$ for two benchmark models with mass ratios $M_χ/M_{A^\prime}=0.6$ and $2$ and for dark fine-structure constants $0.1\leα_D\le 1$. △ Less

Submitted 16 May, 2024; v1 submitted 21 December, 2023; originally announced December 2023.

Spatiotemporal control of high-intensity laser pulses with a plasma lens

Authors: D. Li, K. G. Miller, J. R. Pierce, W. B. Mori, A. G. R. Thomas, J. P. Palastro

Abstract: Spatiotemporal control encompasses a variety of techniques for producing laser pulses with dynamic intensity peaks that move independently of the group velocity. This controlled motion of the intensity peak offers a new approach to optimizing laser-based applications and enhancing signatures of fundamental phenomena. Here, we demonstrate spatiotemporal control with a plasma optic. A chirped laser… ▽ More Spatiotemporal control encompasses a variety of techniques for producing laser pulses with dynamic intensity peaks that move independently of the group velocity. This controlled motion of the intensity peak offers a new approach to optimizing laser-based applications and enhancing signatures of fundamental phenomena. Here, we demonstrate spatiotemporal control with a plasma optic. A chirped laser pulse focused by a plasma lens exhibits a moving focal point, or "flying focus," that can travel at an arbitrary, predetermined velocity. Unlike currently used conventional or adaptive optics, a plasma lens can be located close to the interaction region and can operate at an orders of magnitude higher, near-relativistic intensity. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Brauer pairs for splendid Rickard complexes

Authors: Jadyn V. Breland, Sam K. Miller

Abstract: We define the notion of a Brauer pair of a chain complex, extending the notion of a Brauer pair of a $p$-permutation module introduced by Boltje and Perepelitsky. In fact, the Brauer pairs of a splendid Rickard complex $C$ coincide with the set of Brauer pairs of the corresponding $p$-permutation equivalence $Λ(C)$ induced by $C$. As a result, we derive structural results for splendid Rickard comp… ▽ More We define the notion of a Brauer pair of a chain complex, extending the notion of a Brauer pair of a $p$-permutation module introduced by Boltje and Perepelitsky. In fact, the Brauer pairs of a splendid Rickard complex $C$ coincide with the set of Brauer pairs of the corresponding $p$-permutation equivalence $Λ(C)$ induced by $C$. As a result, we derive structural results for splendid Rickard complexes which are analogous to known structural properties for $p$-permutation equivalences. In particular, we show splendid Rickard equivalences induce local splendid Rickard equivalences between normalizer block algebras as well as centralizer block algebras. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 22 pages, comments welcome! arXiv admin note: text overlap with arXiv:2007.09253 by other authors

On the generation of attosecond gigawatt soft X-ray pulses through coherent Thomson backscattering

Authors: Qianyi Ma, Jiaxin Liu, Zhuo Pan, Xuezhi Wu, Huangang Lu, Zhenan Wang, Yuhui Xia, Yuekai Chen, Kyle Miller, Xinlu Xu, Xueqing Yan

Abstract: Collision between relativistic electron sheets and counter-propagating laser pulses is recognized as a promising way to produce intense attosecond X-rays through coherent Thomson backscattering (TBS). In a double-layer scheme, the electrons in an ultrathin solid foil are first pushed out by an intense laser driver and then interact with the laser reflected off a second foil to form a high-density… ▽ More Collision between relativistic electron sheets and counter-propagating laser pulses is recognized as a promising way to produce intense attosecond X-rays through coherent Thomson backscattering (TBS). In a double-layer scheme, the electrons in an ultrathin solid foil are first pushed out by an intense laser driver and then interact with the laser reflected off a second foil to form a high-density relativistic electron sheet with vanishing transverse momentum. However, the repulsion between these concentrated electrons can increase the thickness of the layer, reducing both its density and subsequently the coherent TBS. Here, we present a systematic study on the evolution of the flying electron layer and find that its resulting thickness is determined by the interplay between the intrinsic space-charge expansion and the velocity compression induced by the drive laser. How the laser driver, the target areal density, the reflector and the collision laser intensity affect the properties of the produced X-rays is explored. Multi-dimensional particle-in-cell simulations indicate that employing this scheme in the nonlinear regime has the potential to stably produce soft X-rays with several GW peak power in hundreds of TW ultrafast laser facilities. The pulse duration can be tuned to tens of attoseconds. This compact and intense attosecond X-ray source may have broad applications in attosecond science. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 10 pages, 6 figures

The composition of Saturn's rings

Authors: Kelly E. Miller, Gianrico Filacchione, Jeffrey Cuzzi, Philip D. Nicholson, Matthew M. Hedman, Kevin Baillie, Robert E. Johnson, Wei-Ling Tseng, Paul R. Estrada, J. Hunter Waite, Mauro Ciarniello, Cécile Ferrari, Zhimeng Zhang, Amanda Hendrix, Julianne I. Moses

Abstract: The origin and evolution of Saturn's rings is critical to understanding the Saturnian system as a whole. Here, we discuss the physical and chemical composition of the rings, as a foundation for evolutionary models described in subsequent chapters. We review the physical characteristics of the main rings, and summarize current constraints on their chemical composition. Radial trends are observed in… ▽ More The origin and evolution of Saturn's rings is critical to understanding the Saturnian system as a whole. Here, we discuss the physical and chemical composition of the rings, as a foundation for evolutionary models described in subsequent chapters. We review the physical characteristics of the main rings, and summarize current constraints on their chemical composition. Radial trends are observed in temperature and to a limited extent in particle size distribution, with the C ring exhibiting higher temperatures and a larger population of small particles. The C ring also shows evidence for the greatest abundance of silicate material, perhaps indicative of formation from a rocky body. The C ring and Cassini Division have lower optical depths than the A and B rings, which contributes to the higher abundance of the exogenous neutral absorber in these regions. Overall, the main ring composition is strongly dominated by water ice, with minor silicate, UV absorber, and neutral absorber components. Sampling of the innermost D ring during Cassini's Grand Finale provides a new set of in situ constraints on the ring composition, and we explore ongoing work to understand the linkages between the main rings and the D ring. The D ring material is organic- and silicate-rich and water-poor relative to the main rings, with a large population of small grains. This composition may be explained in part by volatile losses in the D ring, and current constraints suggest some degree of fractionation rather than sampling of the bulk D ring material. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: Submitted to SSR for publication in the collection "New Vision of the Saturnian System in the Context of a Highly Dissipative Saturn"

Dirichlet Active Learning

Authors: Kevin Miller, Ryan Murray

Abstract: This work introduces Dirichlet Active Learning (DiAL), a Bayesian-inspired approach to the design of active learning algorithms. Our framework models feature-conditional class probabilities as a Dirichlet random field and lends observational strength between similar features in order to calibrate the random field. This random field can then be utilized in learning tasks: in particular, we can use… ▽ More This work introduces Dirichlet Active Learning (DiAL), a Bayesian-inspired approach to the design of active learning algorithms. Our framework models feature-conditional class probabilities as a Dirichlet random field and lends observational strength between similar features in order to calibrate the random field. This random field can then be utilized in learning tasks: in particular, we can use current estimates of mean and variance to conduct classification and active learning in the context where labeled data is scarce. We demonstrate the applicability of this model to low-label rate graph learning by constructing ``propagation operators'' based upon the graph Laplacian, and offer computational studies demonstrating the method's competitiveness with the state of the art. Finally, we provide rigorous guarantees regarding the ability of this approach to ensure both exploration and exploitation, expressed respectively in terms of cluster exploration and increased attention to decision boundaries. △ Less

Submitted 9 November, 2023; originally announced November 2023.

Comments: 66 pages, 16 figures

Conditional Randomization Tests for Behavioral and Neural Time Series

Authors: Kenneth D. Harris, Kevin J. Miller

Abstract: Randomization tests allow simple and unambiguous tests of null hypotheses, by comparing observed data to a null ensemble in which experimentally-controlled variables are randomly resampled. In behavioral and neuroscience experiments, however, the stimuli presented often depend on the subject's previous actions, so simple randomization tests are not possible. We describe how conditional randomizati… ▽ More Randomization tests allow simple and unambiguous tests of null hypotheses, by comparing observed data to a null ensemble in which experimentally-controlled variables are randomly resampled. In behavioral and neuroscience experiments, however, the stimuli presented often depend on the subject's previous actions, so simple randomization tests are not possible. We describe how conditional randomization can be used to perform exact hypothesis tests in this situation, and illustrate it with two examples. We contrast conditional randomization with a related approach of tangent randomization, in which stimuli are resampled based only on events occurring in the past, which is not valid for all choices of test statistic. We discuss how to design experiments that allow conditional randomization tests to be used. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 9 pages, 5 figures

Search for heavy neutral leptons in electron-positron and neutral-pion final states with the MicroBooNE detector

Authors: MicroBooNE collaboration, P. Abratenko, O. Alterkait, D. Andrade Aldana, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, G. Barr, D. Barrow, J. Barrow, V. Basque, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bhattacharya, M. Bishai, A. Blake, B. Bogart, T. Bolton, J. Y. Book, M. B. Brunetti, L. Camilleri , et al. (163 additional authors not shown)

Abstract: We present the first search for heavy neutral leptons (HNL) decaying into $νe^+e^-$ or $νπ^0$ final states in a liquid-argon time projection chamber using data collected with the MicroBooNE detector. The data were recorded synchronously with the NuMI neutrino beam from Fermilab's Main Injector corresponding to a total exposure of $7.01 \times 10^{20}$ protons on target. We set upper limits at the… ▽ More We present the first search for heavy neutral leptons (HNL) decaying into $νe^+e^-$ or $νπ^0$ final states in a liquid-argon time projection chamber using data collected with the MicroBooNE detector. The data were recorded synchronously with the NuMI neutrino beam from Fermilab's Main Injector corresponding to a total exposure of $7.01 \times 10^{20}$ protons on target. We set upper limits at the $90\%$ confidence level on the mixing parameter $\lvert U_{μ4}\rvert^2$ in the mass ranges $10\le m_{\rm HNL}\le 150$ MeV for the $νe^+e^-$ channel and $150\le m_{\rm HNL}\le 245$ MeV for the $νπ^0$ channel, assuming $\lvert U_{e 4}\rvert^2 = \lvert U_{τ4}\rvert^2 = 0$. These limits represent the most stringent constraints in the mass range $35<m_{\rm HNL}<175$ MeV and the first constraints from a direct search for $νπ^0$ decays. △ Less

Submitted 12 January, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: Version as accepted by Physical Review Letters, some presentational changes and updated references, no changes to results

Report number: FERMILAB-PUB-23-574-ND

Measurement of nuclear effects in neutrino-argon interactions using generalized kinematic imbalance variables with the MicroBooNE detector

Authors: MicroBooNE collaboration, P. Abratenko, O. Alterkait, D. Andrade Aldana, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, G. Barr, D. Barrow, J. Barrow, V. Basque, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bhattacharya, M. Bishai, A. Blake, B. Bogart, T. Bolton, J. Y. Book, M. B. Brunetti, L. Camilleri , et al. (163 additional authors not shown)

Abstract: We present a set of new generalized kinematic imbalance variables that can be measured in neutrino scattering. These variables extend previous measurements of kinematic imbalance on the transverse plane, and are more sensitive to modeling of nuclear effects. We demonstrate the enhanced power of these variables using simulation, and then use the MicroBooNE detector to measure them for the first tim… ▽ More We present a set of new generalized kinematic imbalance variables that can be measured in neutrino scattering. These variables extend previous measurements of kinematic imbalance on the transverse plane, and are more sensitive to modeling of nuclear effects. We demonstrate the enhanced power of these variables using simulation, and then use the MicroBooNE detector to measure them for the first time. We report flux-integrated single- and double-differential measurements of charged-current muon neutrino scattering on argon using a topolgy with one muon and one proton in the final state as a function of these novel kinematic imbalance variables. These measurements allow us to demonstrate that the treatment of charged current quasielastic interactions in GENIE version 2 is inadequate to describe data. Further, they reveal tensions with more modern generator predictions particularly in regions of phase space where final state interactions are important. △ Less

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

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

Searching for Optimal Runtime Assurance via Reachability and Reinforcement Learning

Authors: Kristina Miller, Christopher K. Zeitler, William Shen, Kerianne Hobbs, Sayan Mitra, John Schierman, Mahesh Viswanathan

Abstract: A runtime assurance system (RTA) for a given plant enables the exercise of an untrusted or experimental controller while assuring safety with a backup (or safety) controller. The relevant computational design problem is to create a logic that assures safety by switching to the safety controller as needed, while maximizing some performance criteria, such as the utilization of the untrusted controll… ▽ More A runtime assurance system (RTA) for a given plant enables the exercise of an untrusted or experimental controller while assuring safety with a backup (or safety) controller. The relevant computational design problem is to create a logic that assures safety by switching to the safety controller as needed, while maximizing some performance criteria, such as the utilization of the untrusted controller. Existing RTA design strategies are well-known to be overly conservative and, in principle, can lead to safety violations. In this paper, we formulate the optimal RTA design problem and present a new approach for solving it. Our approach relies on reward shaping and reinforcement learning. It can guarantee safety and leverage machine learning technologies for scalability. We have implemented this algorithm and present experimental results comparing our approach with state-of-the-art reachability and simulation-based RTA approaches in a number of scenarios using aircraft models in 3D space with complex safety requirements. Our approach can guarantee safety while increasing utilization of the experimental controller over existing approaches. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Simulation-based Inference with the Generalized Kullback-Leibler Divergence

Authors: Benjamin Kurt Miller, Marco Federici, Christoph Weniger, Patrick Forré

Abstract: In Simulation-based Inference, the goal is to solve the inverse problem when the likelihood is only known implicitly. Neural Posterior Estimation commonly fits a normalized density estimator as a surrogate model for the posterior. This formulation cannot easily fit unnormalized surrogates because it optimizes the Kullback-Leibler divergence. We propose to optimize a generalized Kullback-Leibler di… ▽ More In Simulation-based Inference, the goal is to solve the inverse problem when the likelihood is only known implicitly. Neural Posterior Estimation commonly fits a normalized density estimator as a surrogate model for the posterior. This formulation cannot easily fit unnormalized surrogates because it optimizes the Kullback-Leibler divergence. We propose to optimize a generalized Kullback-Leibler divergence that accounts for the normalization constant in unnormalized distributions. The objective recovers Neural Posterior Estimation when the model class is normalized and unifies it with Neural Ratio Estimation, combining both into a single objective. We investigate a hybrid model that offers the best of both worlds by learning a normalized base distribution and a learned ratio. We also present benchmark results. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: Accepted at Synergy of Scientific and Machine Learning Modeling ICML 2023 Workshop https://syns-ml.github.io/2023/contributions/

CARLA: Adjusted common average referencing for cortico-cortical evoked potential data

Authors: Harvey Huang, Gabriela Ojeda Valencia, Nicholas M. Gregg, Gamaleldin M. Osman, Morgan N. Montoya, Gregory A. Worrell, Kai J. Miller, Dora Hermes

Abstract: Human brain connectivity can be mapped by single pulse electrical stimulation during intracranial EEG measurements. The raw cortico-cortical evoked potentials (CCEP) are often contaminated by noise. Common average referencing (CAR) removes common noise and preserves response shapes but can introduce bias from responsive channels. We address this issue with an adjusted, adaptive CAR algorithm terme… ▽ More Human brain connectivity can be mapped by single pulse electrical stimulation during intracranial EEG measurements. The raw cortico-cortical evoked potentials (CCEP) are often contaminated by noise. Common average referencing (CAR) removes common noise and preserves response shapes but can introduce bias from responsive channels. We address this issue with an adjusted, adaptive CAR algorithm termed "CAR by Least Anticorrelation (CARLA)". CARLA was tested on simulated CCEP data and real CCEP data collected from four human participants. In CARLA, the channels are ordered by increasing mean cross-trial covariance, and iteratively added to the common average until anticorrelation between any single channel and all re-referenced channels reaches a minimum, as a measure of shared noise. We simulated CCEP data with true responses in 0 to 45 of 50 total channels. We quantified CARLA's error and found that it erroneously included 0 (median) truly responsive channels in the common average with less than or equal to 42 responsive channels, and erroneously excluded less than or equal to 2.5 (median) unresponsive channels at all responsiveness levels. On real CCEP data, signal quality was quantified with the mean R-squared between all pairs of channels, which represents inter-channel dependency and is low for well-referenced data. CARLA re-referencing produced significantly lower mean R-squared than standard CAR, CAR using a fixed bottom quartile of channels by covariance, and no re-referencing. CARLA minimizes bias in re-referenced CCEP data by adaptively selecting the optimal subset of non-responsive channels. It showed high specificity and sensitivity on simulated CCEP data and lowered inter-channel dependency compared to CAR on real CCEP data. △ Less

Submitted 29 September, 2023; originally announced October 2023.

Comments: 29 pages, 8 main figures, 3 supplemental figures. For associated code, see https://github.com/hharveygit/CARLA_JNM

Endotrivial complexes

Authors: Sam K. Miller

Abstract: Let $G$ be a finite group, $p$ a prime, and $k$ a field of characteristic $p$. We introduce the notion of an endotrivial chain complex of $p$-permutation $kG$-modules, which are the invertible objects in the bounded homotopy category of $p$-permutation $kG$-modules, and study the corresponding Picard group $\mathcal{E}_k(G)$ of endotrivial complexes. Such complexes are shown to induce splendid Ric… ▽ More Let $G$ be a finite group, $p$ a prime, and $k$ a field of characteristic $p$. We introduce the notion of an endotrivial chain complex of $p$-permutation $kG$-modules, which are the invertible objects in the bounded homotopy category of $p$-permutation $kG$-modules, and study the corresponding Picard group $\mathcal{E}_k(G)$ of endotrivial complexes. Such complexes are shown to induce splendid Rickard autoequivalences of $kG$. The elements of $\mathcal{E}_k(G)$ are determined uniquely by integral invariants arising from the Brauer construction and a degree one character $G \to k^\times$. Using ideas from Bouc's theory of biset functors, we provide a canonical decomposition of $\mathcal{E}_k(G)$, and as an application, give complete descriptions of $\mathcal{E}_k(G)$ for abelian groups and $p$-groups of normal $p$-rank 1. Taking Lefschetz invariants of endotrivial complexes induces a group homomorphism $Λ: \mathcal{E}_k(G) \to O(T(kG))$, where $O(T(kG))$ is the orthogonal unit group of the trivial source ring. Using recent results of Boltje and Carman, we give a Frobenius stability condition elements in the image of $Λ$ must satisfy. △ Less

Submitted 26 May, 2024; v1 submitted 21 September, 2023; originally announced September 2023.

Comments: Accepted to Journal of Algebra. 29 pages

Report number: Journal of Algebra 650 (2024) 173--218 MSC Class: 20J05; 19A22; 20C05; 20C20

Space-time structured plasma waves

Authors: J. P. Palastro, K. G. Miller, R. K. Follett, D. Ramsey, K. Weichman, A. V. Arefiev, D. H. Froula

Abstract: Electrostatic waves play a critical role in nearly every branch of plasma physics from fusion to advanced accelerators, to astro, solar, and ionospheric physics. The properties of planar electrostatic waves are fully determined by the plasma conditions, such as density, temperature, ionization state, or details of the distribution functions. Here we demonstrate that electrostatic wavepackets struc… ▽ More Electrostatic waves play a critical role in nearly every branch of plasma physics from fusion to advanced accelerators, to astro, solar, and ionospheric physics. The properties of planar electrostatic waves are fully determined by the plasma conditions, such as density, temperature, ionization state, or details of the distribution functions. Here we demonstrate that electrostatic wavepackets structured with space-time correlations can have properties that are independent of the plasma conditions. For instance, an appropriately structured electrostatic wavepacket can travel at any group velocity, even backward with respect to its phase fronts, while maintaining a localized energy density. These linear, propagation-invariant wavepackets can be constructed with or without orbital angular momentum by superposing natural modes of the plasma and can be ponderomotively excited by space-time structured laser pulses like the flying focus. △ Less

Submitted 15 September, 2023; v1 submitted 12 September, 2023; originally announced September 2023.

On the effective dynamics of Bose-Fermi mixtures

Authors: Esteban Cárdenas, Joseph K. Miller, Nataša Pavlović

Abstract: In this work, we describe the dynamics of a Bose-Einstein condensate interacting with a degenerate Fermi gas, at zero temperature. First, we analyze the mean-field approximation of the many-body Schrödinger dynamics and prove emergence of a coupled Hartree-type system of equations. We obtain rigorous error control that yields a non-trivial scaling window in which the approximation is meaningful. S… ▽ More In this work, we describe the dynamics of a Bose-Einstein condensate interacting with a degenerate Fermi gas, at zero temperature. First, we analyze the mean-field approximation of the many-body Schrödinger dynamics and prove emergence of a coupled Hartree-type system of equations. We obtain rigorous error control that yields a non-trivial scaling window in which the approximation is meaningful. Second, starting from this Hartree system, we identify a novel scaling regime in which the fermion distribution behaves semi-clasically, but the boson field remains quantum-mechanical; this is one of the main contributions of the present article. In this regime, the bosons are much lighter and more numerous than the fermions. We then prove convergence to a coupled Vlasov-Hartee system of equations with an explicit convergence rate. △ Less

Submitted 21 March, 2024; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 55 pages, 0 figures

MSC Class: 35Q40 35Q41 35Q49 35Q82 81V73 81V74 82C10 82D05

Dephasingless laser wakefield acceleration in the bubble regime

Authors: Kyle G. Miller, Jacob R. Pierce, Manfred V. Ambat, Jessica L. Shaw, Kale Weichman, Warren B. Mori, Dustin H. Froula, John P. Palastro

Abstract: Laser wakefield accelerators (LWFAs) have electric fields that are orders of magnitude larger than those of conventional accelerators, promising an attractive, small-scale alternative for next-generation light sources and lepton colliders. The maximum energy gain in a single-stage LWFA is limited by dephasing, which occurs when the trapped particles outrun the accelerating phase of the wakefield.… ▽ More Laser wakefield accelerators (LWFAs) have electric fields that are orders of magnitude larger than those of conventional accelerators, promising an attractive, small-scale alternative for next-generation light sources and lepton colliders. The maximum energy gain in a single-stage LWFA is limited by dephasing, which occurs when the trapped particles outrun the accelerating phase of the wakefield. Here, we demonstrate that a single space-time structured laser pulse can be used for ionization injection and electron acceleration over many dephasing lengths in the bubble regime. Simulations of a dephasingless laser wakefield accelerator driven by a 6.2-J laser pulse show 25 pC of injected charge accelerated over 20 dephasing lengths (1.3 cm) to a maximum energy of 2.1 GeV. The space-time structured laser pulse features an ultrashort, programmable-trajectory focus. Accelerating the focus, reducing the focused spot-size variation, and mitigating unwanted self-focusing stabilize the electron acceleration, which improves beam quality and leads to projected energy gains of 125 GeV in a single, sub-meter stage driven by a 500-J pulse. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: 18 pages, 4 figures

First application of a liquid argon time projection chamber for the search for intranuclear neutron-antineutron transitions and annihilation in $^{40}$Ar using the MicroBooNE detector

Authors: MicroBooNE collaboration, P. Abratenko, O. Alterkait, D. Andrade Aldana, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, G. Barr, D. Barrow, J. Barrow, V. Basque, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bhattacharya, M. Bishai, A. Blake, B. Bogart, T. Bolton, J. Y. Book, L. Camilleri, Y. Cao , et al. (164 additional authors not shown)

Abstract: We present a novel methodology to search for intranuclear neutron-antineutron transition ($n\rightarrow\bar{n}$) followed by $\bar{n}$-nucleon annihilation within an $^{40}$Ar nucleus, using the MicroBooNE liquid argon time projection chamber (LArTPC) detector. A discovery of $n\rightarrow\bar{n}$ transition or a new best limit on the lifetime of this process would either constitute physics beyond… ▽ More We present a novel methodology to search for intranuclear neutron-antineutron transition ($n\rightarrow\bar{n}$) followed by $\bar{n}$-nucleon annihilation within an $^{40}$Ar nucleus, using the MicroBooNE liquid argon time projection chamber (LArTPC) detector. A discovery of $n\rightarrow\bar{n}$ transition or a new best limit on the lifetime of this process would either constitute physics beyond the Standard Model or greatly constrain theories of baryogenesis, respectively. The approach presented in this paper makes use of deep learning methods to select $n\rightarrow\bar{n}$ events based on their unique features and differentiate them from cosmogenic backgrounds. The achieved signal and background efficiencies are (70.22$\pm$6.04)\% and (0.0020$\pm$0.0003)\%, respectively. A demonstration of a search is performed with a data set corresponding to an exposure of $3.32 \times10^{26}\,$neutron-years, and where the background rate is constrained through direct measurement, assuming the presence of a negligible signal. With this approach, no excess of events over the background prediction is observed, setting a demonstrative lower bound on the $n\rightarrow\bar{n}$ lifetime in $^{40}$Ar of $τ_{\textrm{m}} \gtrsim 1.1\times10^{26}\,$years, and on the free $n\rightarrow\bar{n}$ transition time of $τ_{\textrm{\nnbar}} \gtrsim 2.6\times10^{5}\,$s, each at the $90\%$ confidence level. This analysis represents a first-ever proof-of-principle demonstration of the ability to search for this rare process in LArTPCs with high efficiency and low background. △ Less

Submitted 27 June, 2024; v1 submitted 7 August, 2023; originally announced August 2023.

Comments: 22 pages, 7 figures, 6 tables

Report number: FERMILAB-PUB-23-410-ND

Cluster-aware Semi-supervised Learning: Relational Knowledge Distillation Provably Learns Clustering

Authors: Yijun Dong, Kevin Miller, Qi Lei, Rachel Ward

Abstract: Despite the empirical success and practical significance of (relational) knowledge distillation that matches (the relations of) features between teacher and student models, the corresponding theoretical interpretations remain limited for various knowledge distillation paradigms. In this work, we take an initial step toward a theoretical understanding of relational knowledge distillation (RKD), wit… ▽ More Despite the empirical success and practical significance of (relational) knowledge distillation that matches (the relations of) features between teacher and student models, the corresponding theoretical interpretations remain limited for various knowledge distillation paradigms. In this work, we take an initial step toward a theoretical understanding of relational knowledge distillation (RKD), with a focus on semi-supervised classification problems. We start by casting RKD as spectral clustering on a population-induced graph unveiled by a teacher model. Via a notion of clustering error that quantifies the discrepancy between the predicted and ground truth clusterings, we illustrate that RKD over the population provably leads to low clustering error. Moreover, we provide a sample complexity bound for RKD with limited unlabeled samples. For semi-supervised learning, we further demonstrate the label efficiency of RKD through a general framework of cluster-aware semi-supervised learning that assumes low clustering errors. Finally, by unifying data augmentation consistency regularization into this cluster-aware framework, we show that despite the common effect of learning accurate clusterings, RKD facilitates a "global" perspective through spectral clustering, whereas consistency regularization focuses on a "local" perspective via expansion. △ Less

Submitted 23 October, 2023; v1 submitted 20 July, 2023; originally announced July 2023.

Comments: NeurIPS 2023

Novel Batch Active Learning Approach and Its Application to Synthetic Aperture Radar Datasets

Authors: James Chapman, Bohan Chen, Zheng Tan, Jeff Calder, Kevin Miller, Andrea L. Bertozzi

Abstract: Active learning improves the performance of machine learning methods by judiciously selecting a limited number of unlabeled data points to query for labels, with the aim of maximally improving the underlying classifier's performance. Recent gains have been made using sequential active learning for synthetic aperture radar (SAR) data arXiv:2204.00005. In each iteration, sequential active learning s… ▽ More Active learning improves the performance of machine learning methods by judiciously selecting a limited number of unlabeled data points to query for labels, with the aim of maximally improving the underlying classifier's performance. Recent gains have been made using sequential active learning for synthetic aperture radar (SAR) data arXiv:2204.00005. In each iteration, sequential active learning selects a query set of size one while batch active learning selects a query set of multiple datapoints. While batch active learning methods exhibit greater efficiency, the challenge lies in maintaining model accuracy relative to sequential active learning methods. We developed a novel, two-part approach for batch active learning: Dijkstra's Annulus Core-Set (DAC) for core-set generation and LocalMax for batch sampling. The batch active learning process that combines DAC and LocalMax achieves nearly identical accuracy as sequential active learning but is more efficient, proportional to the batch size. As an application, a pipeline is built based on transfer learning feature embedding, graph learning, DAC, and LocalMax to classify the FUSAR-Ship and OpenSARShip datasets. Our pipeline outperforms the state-of-the-art CNN-based methods. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 16 pages, 7 figures, Preprint

ACM Class: I.2.6; I.2.10; I.4.0; I.4.9

Journal ref: Proc. SPIE. Algorithms for Synthetic Aperture Radar Imagery XXX (Vol. 12520, pp. 96-111). 13 June 2023

A voxelized immersed boundary (VIB) finite element method for accurate and efficient blood flow simulation

Authors: G. C. Bourantas, B. F. Zwick, D. S. Lampropoulos, V. C. Loukopoulos, K. Katsanos, A. A. Dimas, V. N. Burganos, A. Wittek, K. Miller

Abstract: We present an efficient and accurate immersed boundary (IB) finite element (FE) method for internal flow problems with complex geometries (e.g., blood flow in the vascular system). In this study, we use a voxelized flow domain (discretized with hexahedral and tetrahedral elements) instead of a box domain, which is frequently used in IB methods. The proposed method utilizes the well-established inc… ▽ More We present an efficient and accurate immersed boundary (IB) finite element (FE) method for internal flow problems with complex geometries (e.g., blood flow in the vascular system). In this study, we use a voxelized flow domain (discretized with hexahedral and tetrahedral elements) instead of a box domain, which is frequently used in IB methods. The proposed method utilizes the well-established incremental pressure correction scheme (IPCS) FE solver, and the boundary condition-enforced IB (BCE-IB) method to numerically solve the transient, incompressible Navier--Stokes flow equations. We verify the accuracy of our numerical method using the analytical solution for the Poiseuille flow in a cylinder, and the available experimental data (laser Doppler velocimetry) for the flow in a three-dimensional 90° angle tube bend. We further examine the accuracy and applicability of the proposed method by considering flow within complex geometries, such as blood flow in aneurysmal vessels and the aorta, flow configurations that would otherwise be difficult to solve by most IB methods. Our method offers high accuracy, as demonstrated by the verification examples, and high applicability, as demonstrated through the solution of blood flow within complex geometry. The proposed method is efficient, since it is as fast as the traditional finite element method used to solve the Navier--Stokes flow equations, with a small overhead (not more than 5$\%$) due to the numerical solution of a linear system formulated for the IB method. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: arXiv admin note: substantial text overlap with arXiv:2007.02082

Measurement of triple-differential inclusive muon-neutrino charged-current cross section on argon with the MicroBooNE detector

Authors: MicroBooNE Collaboration, P. Abratenko, O. Alterkait, D. Andrade Aldana, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, G. Barr, D. Barrow, J. Barrow, V. Basque, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bhattacharya, M. Bishai, A. Blake, B. Bogart, T. Bolton, J. Y. Book, L. Camilleri, Y. Cao , et al. (165 additional authors not shown)

Abstract: We report the first measurement of the differential cross section $d^{2}σ(E_ν)/ d\cos(θ_μ) dP_μ$ for inclusive muon-neutrino charged-current scattering on argon. This measurement utilizes data from 6.4$\times10^{20}$ protons on target of exposure collected using the MicroBooNE liquid argon time projection chamber located along the Fermilab Booster Neutrino Beam with a mean neutrino energy of appro… ▽ More We report the first measurement of the differential cross section $d^{2}σ(E_ν)/ d\cos(θ_μ) dP_μ$ for inclusive muon-neutrino charged-current scattering on argon. This measurement utilizes data from 6.4$\times10^{20}$ protons on target of exposure collected using the MicroBooNE liquid argon time projection chamber located along the Fermilab Booster Neutrino Beam with a mean neutrino energy of approximately 0.8~GeV. The mapping from reconstructed kinematics to truth quantities, particularly from reconstructed to true neutrino energy, is validated by comparing the distribution of reconstructed hadronic energy in data to that of the model prediction in different muon scattering angle bins after conditional constraint from the muon momentum distribution in data. The success of this validation gives confidence that the missing energy in the MicroBooNE detector is well-modeled in simulation, enabling the unfolding to a triple-differential measurement over muon momentum, muon scattering angle, and neutrino energy. The unfolded measurement covers an extensive phase space, providing a wealth of information useful for future liquid argon time projection chamber experiments measuring neutrino oscillations. Comparisons against a number of commonly used model predictions are included and their performance in different parts of the available phase-space is discussed. △ Less

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

Report number: FERMILAB-PUB-23-368-ND