Showing 1–50 of 256 results for author: Harris, P

The Pinnacle Sets of a Graph

Authors: Chassidy Bozeman, Christine Cheng, Pamela E. Harris, Stephen Lasinis, Shanise Walker

Abstract: We introduce and study the pinnacle sets of a simple graph $G$ with $n$ vertices. Given a bijective vertex labeling $λ\,:\,V(G)\rightarrow [n]$, the label $λ(v)$ of vertex $v$ is a pinnacle of $(G, λ)$ if $λ(v)>λ(w)$ for all vertices $w$ in the neighborhood of $v$. The pinnacle set of $(G, λ)$ contains all the pinnacles of the labeled graph. A subset $S\subseteq[n]$ is a pinnacle set of $G$ if the… ▽ More We introduce and study the pinnacle sets of a simple graph $G$ with $n$ vertices. Given a bijective vertex labeling $λ\,:\,V(G)\rightarrow [n]$, the label $λ(v)$ of vertex $v$ is a pinnacle of $(G, λ)$ if $λ(v)>λ(w)$ for all vertices $w$ in the neighborhood of $v$. The pinnacle set of $(G, λ)$ contains all the pinnacles of the labeled graph. A subset $S\subseteq[n]$ is a pinnacle set of $G$ if there exists a labeling $λ$ such that $S$ is the pinnacle set of $(G,λ)$. Of interest to us is the question: Which subsets of $[n]$ are the pinnacle sets of $G$? Our main results are as follows. We show that when $G$ is connected, $G$ has a size-$k$ pinnacle set if and only if $G$ has an independent set of the same size. Consequently, determining if $G$ has a size-$k$ pinnacle set and determining if $G$ has a particular subset $S$ as a pinnacle set are NP-complete problems. Nonetheless, we completely identify all the pinnacle sets of complete graphs, complete bipartite graphs, cycles and paths. We also present two techniques for deriving new pinnacle sets from old ones that imply a typical graph has many pinnacle sets. Finally, we define a poset on all the size-$k$ pinnacle sets of $G$ and show that it is a join semilattice. If, additionally, the poset has a minimum element, then it is a distributive lattice. We conclude with some open problems for further study. △ Less

Submitted 27 June, 2024; originally announced June 2024.

MSC Class: 05C30; 05C78; 05C38; 06A06; 06A07

Reliable edge machine learning hardware for scientific applications

Authors: Tommaso Baldi, Javier Campos, Ben Hawks, Jennifer Ngadiuba, Nhan Tran, Daniel Diaz, Javier Duarte, Ryan Kastner, Andres Meza, Melissa Quinnan, Olivia Weng, Caleb Geniesse, Amir Gholami, Michael W. Mahoney, Vladimir Loncar, Philip Harris, Joshua Agar, Shuyu Qin

Abstract: Extreme data rate scientific experiments create massive amounts of data that require efficient ML edge processing. This leads to unique validation challenges for VLSI implementations of ML algorithms: enabling bit-accurate functional simulations for performance validation in experimental software frameworks, verifying those ML models are robust under extreme quantization and pruning, and enabling… ▽ More Extreme data rate scientific experiments create massive amounts of data that require efficient ML edge processing. This leads to unique validation challenges for VLSI implementations of ML algorithms: enabling bit-accurate functional simulations for performance validation in experimental software frameworks, verifying those ML models are robust under extreme quantization and pruning, and enabling ultra-fine-grained model inspection for efficient fault tolerance. We discuss approaches to developing and validating reliable algorithms at the scientific edge under such strict latency, resource, power, and area requirements in extreme experimental environments. We study metrics for developing robust algorithms, present preliminary results and mitigation strategies, and conclude with an outlook of these and future directions of research towards the longer-term goal of developing autonomous scientific experimentation methods for accelerated scientific discovery. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: IEEE VLSI Test Symposium 2024 (VTS)

Report number: FERMILAB-CONF-24-0116-CSAID

Metered Parking Functions

Authors: Spencer Daugherty, Pamela E. Harris, Ian Klein, Matt McClinton

Abstract: We introduce a generalization of parking functions called $t$-metered $(m,n)$-parking functions, in which one of $m$ cars parks among $n$ spots per hour then leaves after $t$ hours. We characterize and enumerate these sequences for $t=1$, $t=m-2$, and $t=n-1$, and provide data for other cases. We characterize the $1$-metered parking functions by decomposing them into sections based on which cars a… ▽ More We introduce a generalization of parking functions called $t$-metered $(m,n)$-parking functions, in which one of $m$ cars parks among $n$ spots per hour then leaves after $t$ hours. We characterize and enumerate these sequences for $t=1$, $t=m-2$, and $t=n-1$, and provide data for other cases. We characterize the $1$-metered parking functions by decomposing them into sections based on which cars are unlucky, and enumerate them using a Lucas sequence recursion. Additionally, we establish a new combinatorial interpretation of the numerator of the continued fraction $n-1/(n-1/\cdots)$ ($n$ times) as the number of $1$-metered $(n,n)$-parking functions. We introduce the $(m,n)$-parking function shuffle in order to count $(m-2)$-metered $(m,n)$-parking functions, which also yields an expression for the number of $(m,n)$-parking functions with any given first entry. As a special case, we find that the number of $(m-2)$-metered $(m, m-1)$-parking functions is equal to the sum of the first entries of classical parking function of length $m-1$. We enumerate the $(n-1)$-metered $(m,n)$-parking functions in terms of the number of classical parking functions of length $n$ with certain parking outcomes, which we show are periodic sequences with period $n$. We conclude with an array of open problems. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Arithmetical Structures on Coconut Trees

Authors: Alexander Diaz-Lopez, Brian Ha, Pamela E. Harris, Jonathan Rogers, Theo Koss, Dorian Smith

Abstract: If G is a finite connected graph, then an arithmetical structure on $G$ is a pair of vectors $(\mathbf{d}, \mathbf{r})$ with positive integer entries such that $(\diag(\mathbf{d}) - A)\cdot \mathbf{r} = \mathbf{0}$, where $A$ is the adjacency matrix of $G$ and the entries of $\mathbf{r}$ have no common factor other than $1$. In this paper, we generalize the result of Archer, Bishop, Diaz-Lopez, Ga… ▽ More If G is a finite connected graph, then an arithmetical structure on $G$ is a pair of vectors $(\mathbf{d}, \mathbf{r})$ with positive integer entries such that $(\diag(\mathbf{d}) - A)\cdot \mathbf{r} = \mathbf{0}$, where $A$ is the adjacency matrix of $G$ and the entries of $\mathbf{r}$ have no common factor other than $1$. In this paper, we generalize the result of Archer, Bishop, Diaz-Lopez, García Puente, Glass, and Louwsma on enumerating arithmetical structures on bidents (also called coconut tree graphs $\CT{p}{2}$) to all coconut tree graphs $\CT{p}{s}$ which consists of a path on $p>0$ vertices to which we append $s>0$ leaves to the right most vertex on the path. We also give a characterization of smooth arithmetical structures on coconut trees when given number assignments to the leaf nodes. △ Less

Submitted 16 June, 2024; originally announced June 2024.

Comments: 18 pages, 9 figures, comments are welcomed

MSC Class: 05C50; 05C30

Defective Parking Functions and Young Tableaux

Authors: Rebecca E. Garcia, Pamela E. Harris, Alex Moon, Aaron Ortiz, Lauren J. Quesada, Cynthia Marie Rivera SÁnchez, Dwight Anderson Williams II

Abstract: Recall that a defective $(m,n)$-parking function with defect $d$ is a parking function with $m$ cars attempting to park on a street with $n$ parking spots in which exactly $d$ cars fail to park. We establish a way to compute the defect of a defective $(m,n)$-parking function and show that the defect of a parking function is invariant under the action of $\mathfrak{S}_m$ the symmetric group on… ▽ More Recall that a defective $(m,n)$-parking function with defect $d$ is a parking function with $m$ cars attempting to park on a street with $n$ parking spots in which exactly $d$ cars fail to park. We establish a way to compute the defect of a defective $(m,n)$-parking function and show that the defect of a parking function is invariant under the action of $\mathfrak{S}_m$ the symmetric group on $[m]=\{1,2,\ldots,m\}$. We also show that the set of nondecreasing defective $(m,n)$-parking functions with defect $d$ are in bijection with the set of standard Young tableaux of shape $(n + d, m - d)$. This implies that the number of $\mathfrak{S}_m$-orbits of defective $(m,n)$-parking functions with defect $d$ is given by $\frac{n-m+2d+1}{n+d+1}\binom{m+n}{n+d}$. We also give a multinomial formula for the size of an $\mathfrak{S}_m$-orbit of a nondecreasing $(m,n)$-parking function with defect $d$. We conclude by using these results to give a new formula for the number of defective parking functions. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 17 pages, 3 figures, 1 table

MSC Class: 05A19 (primary) 05A05; 05A15 (secondary)

On the Correspondence Between Integer Sequences and Vacillating Tableaux

Authors: Zhanar Berikkyzy, Pamela E. Harris, Anna Pun, Catherine Yan, Chenchen Zhao

Abstract: A fundamental identity in the representation theory of the partition algebra is $n^k = \sum_λ f^λm_k^λ$ for $n \geq 2k$, where $λ$ ranges over integer partitions of $n$, $f^λ$ is the number of standard Young tableaux of shape $λ$, and $m_k^λ$ is the number of vacillating tableaux of shape $λ$ and length $2k$. Using a combination of RSK insertion and jeu de taquin, Halverson and Lewandowski constru… ▽ More A fundamental identity in the representation theory of the partition algebra is $n^k = \sum_λ f^λm_k^λ$ for $n \geq 2k$, where $λ$ ranges over integer partitions of $n$, $f^λ$ is the number of standard Young tableaux of shape $λ$, and $m_k^λ$ is the number of vacillating tableaux of shape $λ$ and length $2k$. Using a combination of RSK insertion and jeu de taquin, Halverson and Lewandowski constructed a bijection $DI_n^k$ that maps each integer sequence in $[n]^k$ to a pair of tableaux of the same shape, where one is a standard Young tableau and the other is a vacillating tableau. In this paper, we study the fine properties of Halverson and Lewandowski's bijection and explore the correspondence between integer sequences and the vacillating tableaux via the map $DI_n^k$ for general integers $n$ and $k$. In particular, we characterize the integer sequences $\boldsymbol{i}$ whose corresponding shape, $λ$, in the image $DI_n^k(\boldsymbol{i})$, satisfies $λ_1 = n$ or $λ_1 = n-k$. △ Less

Submitted 11 May, 2024; originally announced May 2024.

MSC Class: 05A05; 05E10

Flattened Catalan Words

Authors: Jean-Luc Baril, Pamela E. Harris, José L. Ramírez

Abstract: In this work, we define flattened Catalan words as Catalan words whose runs of weak ascents have leading terms that appear in weakly increasing order. We provide generating functions, formulas, and asymptotic expressions for the number of flattened Catalan words based on the number of runs of ascents (descents), runs of weak ascents (descents), $\ell$-valleys, valleys, symmetric valleys, $\ell$-pe… ▽ More In this work, we define flattened Catalan words as Catalan words whose runs of weak ascents have leading terms that appear in weakly increasing order. We provide generating functions, formulas, and asymptotic expressions for the number of flattened Catalan words based on the number of runs of ascents (descents), runs of weak ascents (descents), $\ell$-valleys, valleys, symmetric valleys, $\ell$-peaks, peaks, and symmetric peaks. △ Less

Submitted 8 May, 2024; originally announced May 2024.

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

MSC Class: 05A15; 05A19

Enumerating runs, valleys, and peaks in Catalan words

Authors: Jean-Luc Baril, Pamela E. Harris, Kimberly J. Harry, Matt McClinton, José L. Ramírez

Abstract: We provide generating functions, formulas, and asymptotic expressions for the number of Catalan words based on the number of runs of ascents (descents), runs of weak ascents (descents), $\ell$-valleys, valleys, symmetric valleys, $\ell$-peaks, peaks, and symmetric peaks. We also establish some bijections with restricted Dyck paths and ordered trees that transports some statistics. We provide generating functions, formulas, and asymptotic expressions for the number of Catalan words based on the number of runs of ascents (descents), runs of weak ascents (descents), $\ell$-valleys, valleys, symmetric valleys, $\ell$-peaks, peaks, and symmetric peaks. We also establish some bijections with restricted Dyck paths and ordered trees that transports some statistics. △ Less

Submitted 8 April, 2024; originally announced April 2024.

MSC Class: 05A15; 05A19

A machine-learning pipeline for real-time detection of gravitational waves from compact binary coalescences

Authors: Ethan Marx, William Benoit, Alec Gunny, Rafia Omer, Deep Chatterjee, Ricco C. Venterea, Lauren Wills, Muhammed Saleem, Eric Moreno, Ryan Raikman, Ekaterina Govorkova, Dylan Rankin, Michael W. Coughlin, Philip Harris, Erik Katsavounidis

Abstract: The promise of multi-messenger astronomy relies on the rapid detection of gravitational waves at very low latencies ($\mathcal{O}$(1\,s)) in order to maximize the amount of time available for follow-up observations. In recent years, neural-networks have demonstrated robust non-linear modeling capabilities and millisecond-scale inference at a comparatively small computational footprint, making them… ▽ More The promise of multi-messenger astronomy relies on the rapid detection of gravitational waves at very low latencies ($\mathcal{O}$(1\,s)) in order to maximize the amount of time available for follow-up observations. In recent years, neural-networks have demonstrated robust non-linear modeling capabilities and millisecond-scale inference at a comparatively small computational footprint, making them an attractive family of algorithms in this context. However, integration of these algorithms into the gravitational-wave astrophysics research ecosystem has proven non-trivial. Here, we present the first fully machine learning-based pipeline for the detection of gravitational waves from compact binary coalescences (CBCs) running in low-latency. We demonstrate this pipeline to have a fraction of the latency of traditional matched filtering search pipelines while achieving state-of-the-art sensitivity to higher-mass stellar binary black holes. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Boolean intervals in the weak Bruhat order of a finite Coxeter group

Authors: Ben Adenbaum, Jennifer Elder, Pamela E. Harris, J. Carlos Martínez Mori

Abstract: Given a Coxeter group $W$ with Coxeter system $(W,S)$, where $S$ is finite. We provide a complete characterization of Boolean intervals in the weak order of $W$ uniformly for all Coxeter groups in terms of independent sets of the Coxeter graph. Moreover, we establish that the number of Boolean intervals of rank $k$ in the weak order of $W$ is ${i_k(Γ_W)\cdot|W|}\,/\,2^{k}$, where $Γ_W$ is the Coxe… ▽ More Given a Coxeter group $W$ with Coxeter system $(W,S)$, where $S$ is finite. We provide a complete characterization of Boolean intervals in the weak order of $W$ uniformly for all Coxeter groups in terms of independent sets of the Coxeter graph. Moreover, we establish that the number of Boolean intervals of rank $k$ in the weak order of $W$ is ${i_k(Γ_W)\cdot|W|}\,/\,2^{k}$, where $Γ_W$ is the Coxeter graph of $W$ and $i_k(Γ_W)$ is the number of independent sets of size $k$ of $Γ_W$ when $W$ is finite. Specializing to $A_n$, we recover the characterizations and enumerations of Boolean intervals in the weak order of $A_n$ given in arXiv:2306.14734. We provide the analogous results for types $C_n$ and $D_n$, including the related generating functions and additional connections to well-known integer sequences. △ Less

Submitted 12 March, 2024; originally announced March 2024.

MSC Class: 05A05; 05A15; 05A19; 06A07; 20F55

Re-Simulation-based Self-Supervised Learning for Pre-Training Foundation Models

Authors: Philip Harris, Michael Kagan, Jeffrey Krupa, Benedikt Maier, Nathaniel Woodward

Abstract: Self-Supervised Learning (SSL) is at the core of training modern large machine learning models, providing a scheme for learning powerful representations that can be used in a variety of downstream tasks. However, SSL strategies must be adapted to the type of training data and downstream tasks required. We propose RS3L, a novel simulation-based SSL strategy that employs a method of re-simulation to… ▽ More Self-Supervised Learning (SSL) is at the core of training modern large machine learning models, providing a scheme for learning powerful representations that can be used in a variety of downstream tasks. However, SSL strategies must be adapted to the type of training data and downstream tasks required. We propose RS3L, a novel simulation-based SSL strategy that employs a method of re-simulation to drive data augmentation for contrastive learning. By intervening in the middle of the simulation process and re-running simulation components downstream of the intervention, we generate multiple realizations of an event, thus producing a set of augmentations covering all physics-driven variations available in the simulator. Using experiments from high-energy physics, we explore how this strategy may enable the development of a foundation model; we show how R3SL pre-training enables powerful performance in downstream tasks such as discrimination of a variety of objects and uncertainty mitigation. In addition to our results, we make the RS3L dataset publicly available for further studies on how to improve SSL strategies. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 24 pages, 9 figures

On the Lucky and Displacement Statistics of Stirling Permutations

Authors: Laura Colmenarejo, Aleyah Dawkins, Jennifer Elder, Pamela E. Harris, Kimberly J. Harry, Selvi Kara, Dorian Smith, Bridget Eileen Tenner

Abstract: Stirling permutations are parking functions, and we investigate two parking function statistics in the context of these objects: lucky cars and displacement. Among our results, we consider two extreme cases: extremely lucky Stirling permutations (those with maximally many lucky cars) and extremely unlucky Stirling permutations (those with exactly one lucky car). We show that the number of extremel… ▽ More Stirling permutations are parking functions, and we investigate two parking function statistics in the context of these objects: lucky cars and displacement. Among our results, we consider two extreme cases: extremely lucky Stirling permutations (those with maximally many lucky cars) and extremely unlucky Stirling permutations (those with exactly one lucky car). We show that the number of extremely lucky Stirling permutations of order $n$ is the Catalan number $C_n$, and the number of extremely unlucky Stirling permutations is $(n-1)!$. We also give some results for luck that lies between these two extremes. Further, we establish that the displacement of any Stirling permutation of order $n$ is $n^2$, and we prove several results about displacement composition vectors. We conclude with directions for further study. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 17 pages, 3 tables

MSC Class: 05A05; 05A15

Prime isogenous discriminant ideal twins

Authors: Alexander J. Barrios, Maila Brucal-Hallare, Alyson Deines, Piper Harris, Manami Roy

Abstract: Let $E_{1}$ and $E_{2}$ be elliptic curves defined over a number field $K$. We say that $E_{1}$ and $E_{2}$ are discriminant ideal twins if they are not $K$-isomorphic and have the same minimal discriminant ideal and conductor. Such curves are said to be discriminant twins if, for each prime $\mathfrak{p}$ of $K$, there are $\mathfrak{p}$-minimal models for $E_{1}$ and $E_{2}$ whose discriminants… ▽ More Let $E_{1}$ and $E_{2}$ be elliptic curves defined over a number field $K$. We say that $E_{1}$ and $E_{2}$ are discriminant ideal twins if they are not $K$-isomorphic and have the same minimal discriminant ideal and conductor. Such curves are said to be discriminant twins if, for each prime $\mathfrak{p}$ of $K$, there are $\mathfrak{p}$-minimal models for $E_{1}$ and $E_{2}$ whose discriminants are equal. This article explicitly classifies all prime-isogenous discriminant (ideal) twins over $\mathbb{Q}$. We obtain this classification as a consequence of our main results, which constructively gives all $p$-isogenous discriminant ideal twins over number fields where $p\in\left\{ 2,3,5,7,13\right\}$, i.e., where $X_0(p)$ has genus $0$. In particular, we find that up to twist, there are finitely many $p$-isogenous discriminant ideal twins if and only if $K$ is $\mathbb{Q}$ or an imaginary quadratic field. In the latter case, we provide instructions for finding the finitely many pairs of $j$-invariants that result in $p$-isogenous discriminant ideal twins. We prove our results by considering the local data of parameterized $p$-isogenous elliptic curves. △ Less

Submitted 29 February, 2024; originally announced February 2024.

Comments: 35 pages

MSC Class: 11G05; 11G07; 14K02; 14H10; 14H52

Graph Neural Network-based Tracking as a Service

Authors: Haoran Zhao, Andrew Naylor, Shih-Chieh Hsu, Paolo Calafiura, Steven Farrell, Yongbing Feng, Philip Coleman Harris, Elham E Khoda, William Patrick Mccormack, Dylan Sheldon Rankin, Xiangyang Ju

Abstract: Recent studies have shown promising results for track finding in dense environments using Graph Neural Network (GNN)-based algorithms. However, GNN-based track finding is computationally slow on CPUs, necessitating the use of coprocessors to accelerate the inference time. Additionally, the large input graph size demands a large device memory for efficient computation, a requirement not met by all… ▽ More Recent studies have shown promising results for track finding in dense environments using Graph Neural Network (GNN)-based algorithms. However, GNN-based track finding is computationally slow on CPUs, necessitating the use of coprocessors to accelerate the inference time. Additionally, the large input graph size demands a large device memory for efficient computation, a requirement not met by all computing facilities used for particle physics experiments, particularly those lacking advanced GPUs. Furthermore, deploying the GNN-based track-finding algorithm in a production environment requires the installation of all dependent software packages, exclusively utilized by this algorithm. These computing challenges must be addressed for the successful implementation of GNN-based track-finding algorithm into production settings. In response, we introduce a ``GNN-based tracking as a service'' approach, incorporating a custom backend within the NVIDIA Triton inference server to facilitate GNN-based tracking. This paper presents the performance of this approach using the Perlmutter supercomputer at NERSC. △ Less

Submitted 14 February, 2024; originally announced February 2024.

Comments: 7 pages, 4 figures, Proceeding of Connected the Dots Workshop (CTD 2023)

Report number: PROC-CTD2023-56

Polarization of circumstellar debris disk light echoes

Authors: Austin J. King, Benjamin C. Bromley, Preston W. Harris, Scott J. Kenyon

Abstract: Light echoes of debris disks around active stars can reveal disk structure and composition even when disks are not spatially resolved. Unfortunately, distinguishing reflected light from quiescent starlight and unexpected post-peak flare structure is challenging, especially for edge-on geometries where the time delay between observed flare photons and light scattered from the near side of the disk… ▽ More Light echoes of debris disks around active stars can reveal disk structure and composition even when disks are not spatially resolved. Unfortunately, distinguishing reflected light from quiescent starlight and unexpected post-peak flare structure is challenging, especially for edge-on geometries where the time delay between observed flare photons and light scattered from the near side of the disk is short. Here, we take advantage of the fact that scattered light from a dusty disk is polarized, depending on the location of the scattering site and the orientation of the disk relative to a distant observer. Filtering reflected light into its polarized components allows echoes to stand out in predictable ways. We test this idea with a simple model for a disk around an active M dwarf. Our results demonstrate that the use of polarimetric data of flaring stars can significantly enhance echo signals relative to starlight and yield more robust and accurate fits to disk parameters compared to analyses based on the total intensity alone. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 14 pages, 9 figures, 1 table, AJ accepted

Vacillating parking functions

Authors: Bruce Fang, Pamela E. Harris, Brian M. Kamau, David Wang

Abstract: For any integers $1\leq k\leq n$, we introduce a new family of parking functions called $k$-vacillating parking functions of length $n$. The parking rule for $k$-vacillating parking functions allows a car with preference $p$ to park in the first available spot in encounters among the parking spots numbered $p$, $p-k$, and $p+k$ (in that order and if those spots exists). In this way, $k$-vacillatin… ▽ More For any integers $1\leq k\leq n$, we introduce a new family of parking functions called $k$-vacillating parking functions of length $n$. The parking rule for $k$-vacillating parking functions allows a car with preference $p$ to park in the first available spot in encounters among the parking spots numbered $p$, $p-k$, and $p+k$ (in that order and if those spots exists). In this way, $k$-vacillating parking functions are a modification of Naples parking functions, which allow for backwards movement of a car, and of $\ell$-interval parking functions, which allow a car to park in its preference or up to $\ell$ spots in front of its preference. Among our results, we establish a combinatorial interpretation for the numerator of the $n$th convergent of the continued fraction of $\sqrt{2}$, as the number of non-decreasing $1$-vacillating parking functions of length~$n$. Our main result gives a product formula for the enumeration of $k$-vacillating parking functions of length $n$ based on the number of $1$-vacillating parking functions of smaller length. We also give closed formulas for the number of $k$-vacillating parking functions of length $n$ which are monotonically increasing or decreasing. We conclude with some directions for further research. △ Less

Submitted 4 February, 2024; originally announced February 2024.

Comments: 13 pages, 1 figure

MSC Class: 05A05; 05A15

Ultra Fast Transformers on FPGAs for Particle Physics Experiments

Authors: Zhixing Jiang, Dennis Yin, Elham E Khoda, Vladimir Loncar, Ekaterina Govorkova, Eric Moreno, Philip Harris, Scott Hauck, Shih-Chieh Hsu

Abstract: This work introduces a highly efficient implementation of the transformer architecture on a Field-Programmable Gate Array (FPGA) by using the \texttt{hls4ml} tool. Given the demonstrated effectiveness of transformer models in addressing a wide range of problems, their application in experimental triggers within particle physics becomes a subject of significant interest. In this work, we have imple… ▽ More This work introduces a highly efficient implementation of the transformer architecture on a Field-Programmable Gate Array (FPGA) by using the \texttt{hls4ml} tool. Given the demonstrated effectiveness of transformer models in addressing a wide range of problems, their application in experimental triggers within particle physics becomes a subject of significant interest. In this work, we have implemented critical components of a transformer model, such as multi-head attention and softmax layers. To evaluate the effectiveness of our implementation, we have focused on a particle physics jet flavor tagging problem, employing a public dataset. We recorded latency under 2 $μ$s on the Xilinx UltraScale+ FPGA, which is compatible with hardware trigger requirements at the CERN Large Hadron Collider experiments. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 6 pages, 2 figures

Journal ref: Machine Learning and the Physical Sciences Workshop, NeurIPS 2023

SymbolNet: Neural Symbolic Regression with Adaptive Dynamic Pruning

Authors: Ho Fung Tsoi, Vladimir Loncar, Sridhara Dasu, Philip Harris

Abstract: Contrary to the use of genetic programming, the neural network approach to symbolic regression can scale well with high input dimension and leverage gradient methods for faster equation searching. Common ways of constraining expression complexity have relied on multistage pruning methods with fine-tuning, but these often lead to significant performance loss. In this work, we propose SymbolNet, a n… ▽ More Contrary to the use of genetic programming, the neural network approach to symbolic regression can scale well with high input dimension and leverage gradient methods for faster equation searching. Common ways of constraining expression complexity have relied on multistage pruning methods with fine-tuning, but these often lead to significant performance loss. In this work, we propose SymbolNet, a neural network approach to symbolic regression in a novel framework that enables dynamic pruning of model weights, input features, and mathematical operators in a single training, where both training loss and expression complexity are optimized simultaneously. We introduce a sparsity regularization term per pruning type, which can adaptively adjust its own strength and lead to convergence to a target sparsity level. In contrast to most existing symbolic regression methods that cannot efficiently handle datasets with more than $O$(10) inputs, we demonstrate the effectiveness of our model on the LHC jet tagging task (16 inputs), MNIST (784 inputs), and SVHN (3072 inputs). △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: 11 pages. Submitted to IEEE TNNLS, under review

Unit interval parking functions and the $r$-Fubini numbers

Authors: S. Alex Bradt, Jennifer Elder, Pamela E. Harris, Gordon Rojas Kirby, Eva Reutercrona, Yuxuan, Wang, Juliet Whidden

Abstract: We recall that unit interval parking functions of length $n$ are a subset of parking functions in which every car parks in its preference or in the spot after its preference, and Fubini rankings of length $n$ are rankings of $n$ competitors allowing for ties. We present an independent proof of a result of Hadaway, which establishes that unit interval parking functions and Fubini rankings are in bi… ▽ More We recall that unit interval parking functions of length $n$ are a subset of parking functions in which every car parks in its preference or in the spot after its preference, and Fubini rankings of length $n$ are rankings of $n$ competitors allowing for ties. We present an independent proof of a result of Hadaway, which establishes that unit interval parking functions and Fubini rankings are in bijection. We also prove that the cardinality of these sets are given by Fubini numbers. In addition, we give a complete characterization of unit interval parking functions by determining when a rearrangement of a unit interval parking function is again a unit interval parking function. This yields an identity for the Fubini numbers as a sum of multinomials over compositions. Moreover, we introduce a generalization of Fubini rankings, which we call the $r$-Fubini rankings of length $n+r$. We show that this set is in bijection with unit interval parking functions of length $n+r$ where the first $r$ cars have distinct preferences. We conclude by establishing that these sets are enumerated by the $r$-Fubini numbers. △ Less

Submitted 12 January, 2024; originally announced January 2024.

Comments: 11 pages, 1 table

MSC Class: Primary: 05A05; Secondary 05A19

On some discrete statistics of parking functions

Authors: Ari Cruz, Pamela E. Harris, Kimberly J. Harry, Jan Kretschmann, Matt McClinton, Alex Moon, John O. Museus, Eric Redmon

Abstract: Recall that $α=(a_1,a_2,\ldots,a_n)\in[n]^n$ is a parking function if its nondecreasing rearrangement $β=(b_1,b_2,\ldots,b_n)$ satisfies $b_i\leq i$ for all $1\leq i\leq n$. In this article, we study parking functions based on their ascents (indices at which $a_i<a_{i+1}$), descents (indices at which $a_i>a_{i+1}$), and ties (indices at which $a_i=a_{i+1}$). By utilizing multiset Eulerian polynomi… ▽ More Recall that $α=(a_1,a_2,\ldots,a_n)\in[n]^n$ is a parking function if its nondecreasing rearrangement $β=(b_1,b_2,\ldots,b_n)$ satisfies $b_i\leq i$ for all $1\leq i\leq n$. In this article, we study parking functions based on their ascents (indices at which $a_i<a_{i+1}$), descents (indices at which $a_i>a_{i+1}$), and ties (indices at which $a_i=a_{i+1}$). By utilizing multiset Eulerian polynomials, we give a generating function for the number of parking functions of length $n$ with $i$ descents. We present a recursive formula for the number of parking functions of length $n$ with descents at a specified subset of $[n-1]$. We establish that the number of parking functions of length $n$ with descents at $I\subset[n-1]$ and descents at $J=\{n-i:i\in I\}$ are equinumerous. As a special case, we show that the number of parking functions of length $n$ with descents at the first $k$ indices is given by $f(n, n-k-1)=\frac{1}{n}\binom{n}{k}\binom{2n-k}{n-k-1}$. We prove this by bijecting to the set of standard Young tableaux of shape $((n-k)^2,1^k)$, which are enumerated by $f(n,n-k-1)$. We also study peaks of parking functions, which are indices at which $a_{i-1}<a_i>a_{i+1}$. We show that the set of parking functions with no peaks and no ties is enumerated by the Catalan numbers. We conclude our study by characterizing when a parking function is uniquely determined by their statistic encoding; a word indicating what indices in the parking function are ascents, descents, and ties. We provide open problems throughout. △ Less

Submitted 24 May, 2024; v1 submitted 27 December, 2023; originally announced December 2023.

Comments: 17 pages, 2 figures, 5 tables, version 2 provides a new proof of Theorem 3.8

MSC Class: 05A05; 05A15

A demonstrator for a real-time AI-FPGA-based triggering system for sPHENIX at RHIC

Authors: J. Kvapil, G. Borca-Tasciuc, H. Bossi, K. Chen, Y. Chen, Y. Corrales Morales, H. Da Costa, C. Da Silva, C. Dean, J. Durham, S. Fu, C. Hao, P. Harris, O. Hen, H. Jheng, Y. Lee, P. Li, X. Li, Y. Lin, M. X. Liu, A. Olvera, M. L. Purschke, M. Rigatti, G. Roland, J. Schambach , et al. (6 additional authors not shown)

Abstract: The RHIC interaction rate at sPHENIX will reach around 3 MHz in pp collisions and requires the detector readout to reject events by a factor of over 200 to fit the DAQ bandwidth of 15 kHz. Some critical measurements, such as heavy flavor production in pp collisions, often require the analysis of particles produced at low momentum. This prohibits adopting the traditional approach, where data rates… ▽ More The RHIC interaction rate at sPHENIX will reach around 3 MHz in pp collisions and requires the detector readout to reject events by a factor of over 200 to fit the DAQ bandwidth of 15 kHz. Some critical measurements, such as heavy flavor production in pp collisions, often require the analysis of particles produced at low momentum. This prohibits adopting the traditional approach, where data rates are reduced through triggering on rare high momentum probes. We explore a new approach based on real-time AI technology, adopt an FPGA-based implementation using a custom designed FELIX-712 board with the Xilinx Kintex Ultrascale FPGA, and deploy the system in the detector readout electronics loop for real-time trigger decision. △ Less

Submitted 27 December, 2023; v1 submitted 22 December, 2023; originally announced December 2023.

Comments: 7 pages, 5 figures, proceedings for TWEPP 2023 conference, v2: corrected Table 1 numbers

Report number: LA-UR-23-32546

New Constraints on Axion-Like Particles from IXPE Polarization Data for Magnetars

Authors: Ephraim Gau, Fazlollah Hajkarim, Steven P. Harris, P. S. Bhupal Dev, Jean-Francois Fortin, Henric Krawczynski, Kuver Sinha

Abstract: We derive new constraints on axion-like particles (ALPs) using precision $X$-ray polarization studies of magnetars. Specifically, we use the first detection of polarized $X$-rays from the magnetars 4U 0142+61 and 1RXS J170849.0-400910 by the Imaging $X$-ray Polarimetry Explorer (IXPE) to place bounds on the product of the ALP-photon and ALP-nucleon couplings, $g_{aγ}g_{aN}$, with $g_{aN}$ being re… ▽ More We derive new constraints on axion-like particles (ALPs) using precision $X$-ray polarization studies of magnetars. Specifically, we use the first detection of polarized $X$-rays from the magnetars 4U 0142+61 and 1RXS J170849.0-400910 by the Imaging $X$-ray Polarimetry Explorer (IXPE) to place bounds on the product of the ALP-photon and ALP-nucleon couplings, $g_{aγ}g_{aN}$, with $g_{aN}$ being responsible for ALP production in the core of the magnetar and $g_{aγ}$ controlling the ALP-photon conversion probability in the magnetosphere. These bounds are most sensitive to the magnetar core temperature, and we use two benchmark values of $1\times 10^8$ K and $5\times 10^8$ K to derive our constraints. For the latter choice, our bounds are competitive with the existing bounds on the coupling product coming from a combination of CAST (for $g_{aγ}$) and SN1987A (for $g_{aN}$). We advocate for more precise and extensive observational campaigns in the higher end of the $2~-~8~$keV spectral window, where ALP-induced polarization is the strongest. We further advocate for hard $X$-ray polarization studies of young, hot, near-Earth magnetars with strong magnetic fields. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 27 pages, 9 figures, 3 tables

Report number: INT-PUB-23-053

Modelling of Networked Measuring Systems -- From White-Box Models to Data Based Approaches

Authors: Klaus-Dieter Sommer, Peter Harris, Sascha Eichstädt, Roland Füssl, Tanja Dorst, Andreas Schütze, Michael Heizmann, Nadine Schiering, Andreas Maier, Yuhui Luo, Christos Tachtatzis, Ivan Andonovic, Gordon Gourlay

Abstract: Mathematical modelling is at the core of metrology as it transforms raw measured data into useful measurement results. A model captures the relationship between the measurand and all relevant quantities on which the measurand depends, and is used to design measuring systems, analyse measured data, make inferences and predictions, and is the basis for evaluating measurement uncertainties. Tradition… ▽ More Mathematical modelling is at the core of metrology as it transforms raw measured data into useful measurement results. A model captures the relationship between the measurand and all relevant quantities on which the measurand depends, and is used to design measuring systems, analyse measured data, make inferences and predictions, and is the basis for evaluating measurement uncertainties. Traditional modelling approaches are typically analytical, for example, based on principles of physics. But with the increasing use of digital technologies, large sensor networks and powerful computing hardware, these traditional approaches are being replaced more and more by data-driven methods. This paradigm shift holds true in particular for the digital future of measurement in all spheres of our lives and the environment, where data provided by large and complex interconnected systems of sensors are to be analysed. Additionally, there is a requirement for existing guidelines and standards in metrology to take the paradigm shift into account. In this paper we lay the foundation for the development from traditional to data-driven modelling approaches. We identify key aspects from traditional modelling approaches and discuss their transformation to data-driven modelling. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Optimizing Likelihood-free Inference using Self-supervised Neural Symmetry Embeddings

Authors: Deep Chatterjee, Philip C. Harris, Maanas Goel, Malina Desai, Michael W. Coughlin, Erik Katsavounidis

Abstract: Likelihood-free inference is quickly emerging as a powerful tool to perform fast/effective parameter estimation. We demonstrate a technique of optimizing likelihood-free inference to make it even faster by marginalizing symmetries in a physical problem. In this approach, physical symmetries, for example, time-translation are learned using joint-embedding via self-supervised learning with symmetry… ▽ More Likelihood-free inference is quickly emerging as a powerful tool to perform fast/effective parameter estimation. We demonstrate a technique of optimizing likelihood-free inference to make it even faster by marginalizing symmetries in a physical problem. In this approach, physical symmetries, for example, time-translation are learned using joint-embedding via self-supervised learning with symmetry data augmentations. Subsequently, parameter inference is performed using a normalizing flow where the embedding network is used to summarize the data before conditioning the parameters. We present this approach on two simple physical problems and we show faster convergence in a smaller number of parameters compared to a normalizing flow that does not use a pre-trained symmetry-informed representation. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: Accepted for Machine Learning and the Physical Sciences Workshop (submission 69) at NeurIPS 2023; for codes, see https://github.com/ML4GW/summer-projects-2023/blob/neurips-2023/symmetry-informed-flows/README.md

Interval and $\ell$-interval Rational Parking Functions

Authors: Tomás Aguilar-Fraga, Jennifer Elder, Rebecca E. Garcia, Kimberly P. Hadaway, Pamela E. Harris, Kimberly J. Harry, Imhotep B. Hogan, Jakeyl Johnson, Jan Kretschmann, Kobe Lawson-Chavanu, J. Carlos Martínez Mori, Casandra D. Monroe, Daniel Quiñonez, Dirk Tolson III, Dwight Anderson Williams II

Abstract: Interval parking functions are a generalization of parking functions in which cars have an interval preference for their parking. We generalize this definition to parking functions with $n$ cars and $m\geq n$ parking spots, which we call interval rational parking functions and provide a formula for their enumeration. By specifying an integer parameter $\ell\geq 0$, we then consider the subset of i… ▽ More Interval parking functions are a generalization of parking functions in which cars have an interval preference for their parking. We generalize this definition to parking functions with $n$ cars and $m\geq n$ parking spots, which we call interval rational parking functions and provide a formula for their enumeration. By specifying an integer parameter $\ell\geq 0$, we then consider the subset of interval rational parking functions in which each car parks at most $\ell$ spots away from their initial preference. We call these $\ell$-interval rational parking functions and provide recursive formulas to enumerate this set for all positive integers $m\geq n$ and $\ell$. We also establish formulas for the number of nondecreasing $\ell$-interval rational parking functions via the outcome map on rational parking functions. We also consider the intersection between $\ell$-interval parking functions and Fubini rankings and show the enumeration of these sets is given by generalized Fibonacci numbers. We conclude by specializing $\ell=1$, and establish that the set of $1$-interval rational parking functions with $n$ cars and $m$ spots are in bijection with the set of barred preferential arrangements of $[n]$ with $m-n$ bars. This readily implies enumerative formulas. Further, in the case where $\ell=1$, we recover the results of Hadaway and Harris that unit interval parking functions are in bijection with the set of Fubini rankings, which are enumerated by the Fubini numbers. △ Less

Submitted 24 May, 2024; v1 submitted 23 November, 2023; originally announced November 2023.

Comments: New updated version of Corollary 4.8

MSC Class: 05A05; 05A15; 05A18; 05A19

Knowledge Distillation for Anomaly Detection

Authors: Adrian Alan Pol, Ekaterina Govorkova, Sonja Gronroos, Nadezda Chernyavskaya, Philip Harris, Maurizio Pierini, Isobel Ojalvo, Peter Elmer

Abstract: Unsupervised deep learning techniques are widely used to identify anomalous behaviour. The performance of such methods is a product of the amount of training data and the model size. However, the size is often a limiting factor for the deployment on resource-constrained devices. We present a novel procedure based on knowledge distillation for compressing an unsupervised anomaly detection model int… ▽ More Unsupervised deep learning techniques are widely used to identify anomalous behaviour. The performance of such methods is a product of the amount of training data and the model size. However, the size is often a limiting factor for the deployment on resource-constrained devices. We present a novel procedure based on knowledge distillation for compressing an unsupervised anomaly detection model into a supervised deployable one and we suggest a set of techniques to improve the detection sensitivity. Compressed models perform comparably to their larger counterparts while significantly reducing the size and memory footprint. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Achieving ultra-low and -uniform residual magnetic fields in a very large magnetically shielded room for fundamental physics experiments

Authors: N. J. Ayres, G. Ban, G. Bison, K. Bodek, V. Bondar, T. Bouillaud, D. Bowles, E. Chanel, W. Chen, P. -J. Chiu, C. B. Crawford, O. Naviliat-Cuncic, C. B. Doorenbos, S. Emmenegger, M. Fertl, A. Fratangelo, W. C. Griffith, Z. D. Grujic, P. G. Harris, K. Kirch, V. Kletzl, J. Krempel, B. Lauss, T. Lefort, A. Lejuez , et al. (25 additional authors not shown)

Abstract: High-precision searches for an electric dipole moment of the neutron (nEDM) require stable and uniform magnetic field environments. We present the recent achievements of degaussing and equilibrating the magnetically shielded room (MSR) for the n2EDM experiment at the Paul Scherrer Institute. We present the final degaussing configuration that will be used for n2EDM after numerous studies. The optim… ▽ More High-precision searches for an electric dipole moment of the neutron (nEDM) require stable and uniform magnetic field environments. We present the recent achievements of degaussing and equilibrating the magnetically shielded room (MSR) for the n2EDM experiment at the Paul Scherrer Institute. We present the final degaussing configuration that will be used for n2EDM after numerous studies. The optimized procedure results in a residual magnetic field that has been reduced by a factor of two. The ultra-low field is achieved with the full magnetic-field-coil system, and a large vacuum vessel installed, both in the MSR. In the inner volume of ~1.4 m^3, the field is now more uniform and below 300 pT. In addition, the procedure is faster and dissipates less heat into the magnetic environment, which in turn, reduces its thermal relaxation time from 12 h down to ~1.5 h. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Chained Quantile Morphing with Normalizing Flows

Authors: Samuel Bright-Thonney, Philip Harris, Patrick McCormack, Simon Rothman

Abstract: Accounting for inaccuracies in Monte Carlo simulations is a crucial step in any high energy physics analysis. It becomes especially important when training machine learning models, which can amplify simulation inaccuracies and introduce large discrepancies and systematic uncertainties when the model is applied to data. In this paper, we introduce a method to transform simulated events to better ma… ▽ More Accounting for inaccuracies in Monte Carlo simulations is a crucial step in any high energy physics analysis. It becomes especially important when training machine learning models, which can amplify simulation inaccuracies and introduce large discrepancies and systematic uncertainties when the model is applied to data. In this paper, we introduce a method to transform simulated events to better match data using normalizing flows, a class of deep learning-based density estimation models. Our proposal uses a technique called chained quantile morphing, which corrects a set of observables by iteratively shifting each entry according to a conditonal cumulative density function. We demonstrate the technique on a realistic particle physics dataset, and compare it to a neural network-based reweighting method. We also introduce a new contrastive learning technique to correct high dimensional particle-level inputs, which naively cannot be efficiently corrected with morphing strategies. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: 13 pages, 7 figures, 3 tables

Cost-sharing in Parking Games

Authors: Jennifer Elder, Pamela E. Harris, Jan Kretschmann, J. Carlos Martínez Mori

Abstract: In this paper, we study the total displacement statistic of parking functions from the perspective of cooperative game theory. We introduce parking games, which are coalitional cost-sharing games in characteristic function form derived from the total displacement statistic. We show that parking games are supermodular cost-sharing games, indicating that cooperation is difficult (i.e., their core is… ▽ More In this paper, we study the total displacement statistic of parking functions from the perspective of cooperative game theory. We introduce parking games, which are coalitional cost-sharing games in characteristic function form derived from the total displacement statistic. We show that parking games are supermodular cost-sharing games, indicating that cooperation is difficult (i.e., their core is empty). Next, we study their Shapley value, which formalizes a notion of "fair" cost-sharing and amounts to charging each car for its expected marginal displacement under a random arrival order. Our main contribution is a polynomial-time algorithm to compute the Shapley value of parking games, in contrast with known hardness results on computing the Shapley value of arbitrary games. The algorithm leverages the permutation-invariance of total displacement, combinatorial enumeration, and dynamic programming. We conclude with open questions around alternative solution concepts for supermodular cost-sharing games and connections to other areas in combinatorics. △ Less

Submitted 14 November, 2023; v1 submitted 21 September, 2023; originally announced September 2023.

Comments: 12 pages

MSC Class: 05A05; 91A12; 91A46

GWAK: Gravitational-Wave Anomalous Knowledge with Recurrent Autoencoders

Authors: Ryan Raikman, Eric A. Moreno, Ekaterina Govorkova, Ethan J Marx, Alec Gunny, William Benoit, Deep Chatterjee, Rafia Omer, Muhammed Saleem, Dylan S Rankin, Michael W Coughlin, Philip C Harris, Erik Katsavounidis

Abstract: Matched-filtering detection techniques for gravitational-wave (GW) signals in ground-based interferometers rely on having well-modeled templates of the GW emission. Such techniques have been traditionally used in searches for compact binary coalescences (CBCs), and have been employed in all known GW detections so far. However, interesting science cases aside from compact mergers do not yet have ac… ▽ More Matched-filtering detection techniques for gravitational-wave (GW) signals in ground-based interferometers rely on having well-modeled templates of the GW emission. Such techniques have been traditionally used in searches for compact binary coalescences (CBCs), and have been employed in all known GW detections so far. However, interesting science cases aside from compact mergers do not yet have accurate enough modeling to make matched filtering possible, including core-collapse supernovae and sources where stochasticity may be involved. Therefore the development of techniques to identify sources of these types is of significant interest. In this paper, we present a method of anomaly detection based on deep recurrent autoencoders to enhance the search region to unmodeled transients. We use a semi-supervised strategy that we name Gravitational Wave Anomalous Knowledge (GWAK). While the semi-supervised nature of the problem comes with a cost in terms of accuracy as compared to supervised techniques, there is a qualitative advantage in generalizing experimental sensitivity beyond pre-computed signal templates. We construct a low-dimensional embedded space using the GWAK method, capturing the physical signatures of distinct signals on each axis of the space. By introducing signal priors that capture some of the salient features of GW signals, we allow for the recovery of sensitivity even when an unmodeled anomaly is encountered. We show that regions of the GWAK space can identify CBCs, detector glitches and also a variety of unmodeled astrophysical sources. △ Less

Submitted 20 September, 2023; originally announced September 2023.

Combinatorial Identities for Vacillating Tableaux

Authors: Zhanar Berikkyzy, Pamela E. Harris, Anna Pun, Catherine Yan, Chenchen Zhao

Abstract: Vacillating tableaux are sequences of integer partitions that satisfy specific conditions. The concept of vacillating tableaux stems from the representation theory of the partition algebra and the combinatorial theory of crossings and nestings of matchings and set partitions. In this paper, we further investigate the enumeration of vacillating tableaux and derive multiple combinatorial identities… ▽ More Vacillating tableaux are sequences of integer partitions that satisfy specific conditions. The concept of vacillating tableaux stems from the representation theory of the partition algebra and the combinatorial theory of crossings and nestings of matchings and set partitions. In this paper, we further investigate the enumeration of vacillating tableaux and derive multiple combinatorial identities and integer sequences relating to the number of vacillating tableaux, simplified vacillating tableaux, and limiting vacillating tableaux. △ Less

Submitted 27 August, 2023; originally announced August 2023.

MSC Class: 05A15; 05A18; 05E10

A large 'Active Magnetic Shield' for a high-precision experiment

Authors: C. Abel, N. J. Ayres, G. Ban, G. Bison, K. Bodek, V. Bondar, T. Bouillaud, E. Chanel, J. Chen, W. Chen, P. -J. Chiu, C. B. Crawford, M. Daum, C. B. Doorenbos, S. Emmenegger, L. Ferraris-Bouchez, M. Fertl, A. Fratangelo, W. C. Griffith, Z. D. Grujic, P. Harris, K. Kirch, V. Kletzl, P. A. Koss, J. Krempel , et al. (26 additional authors not shown)

Abstract: We present a novel Active Magnetic Shield (AMS), designed and implemented for the n2EDM experiment at the Paul Scherrer Institute. The experiment will perform a high-sensitivity search for the electric dipole moment of the neutron. Magnetic-field stability and control is of key importance for n2EDM. A large, cubic, 5m side length, magnetically shielded room (MSR) provides a passive, quasi-static s… ▽ More We present a novel Active Magnetic Shield (AMS), designed and implemented for the n2EDM experiment at the Paul Scherrer Institute. The experiment will perform a high-sensitivity search for the electric dipole moment of the neutron. Magnetic-field stability and control is of key importance for n2EDM. A large, cubic, 5m side length, magnetically shielded room (MSR) provides a passive, quasi-static shielding-factor of about 10^5 for its inner sensitive volume. The AMS consists of a system of eight complex, feedback-controlled compensation coils constructed on an irregular grid spanned on a volume of less than 1000m^3 around the MSR. The AMS is designed to provide a stable and uniform magnetic-field environment around the MSR, while being reasonably compact. The system can compensate static and variable magnetic fields up to +-50muT (homogeneous components) and +-5muT (first-order gradients), suppressing them to a few muT in the sub-Hertz frequency range. The presented design concept and implementation of the AMS fulfills the requirements of the n2EDM experiment and can be useful for other applications, where magnetically silent environments are important and spatial constraints inhibit simpler geometrical solutions. △ Less

Submitted 14 July, 2023; originally announced July 2023.

Parking functions, Fubini rankings, and Boolean intervals in the weak order of $\mathfrak{S}_n$

Authors: Jennifer Elder, Pamela E. Harris, Jan Kretschmann, J. Carlos Martínez Mori

Abstract: Let $\mathfrak{S}_n$ denote the symmetric group and let $W(\mathfrak{S}_n)$ denote the weak order of $\mathfrak{S}_n$. Through a surprising connection to a subset of parking functions, which we call unit Fubini rankings, we provide a complete characterization and enumeration for the total number of Boolean intervals in $W(\mathfrak{S}_n)$ and the total number of Boolean intervals of rank $k$ in… ▽ More Let $\mathfrak{S}_n$ denote the symmetric group and let $W(\mathfrak{S}_n)$ denote the weak order of $\mathfrak{S}_n$. Through a surprising connection to a subset of parking functions, which we call unit Fubini rankings, we provide a complete characterization and enumeration for the total number of Boolean intervals in $W(\mathfrak{S}_n)$ and the total number of Boolean intervals of rank $k$ in $W(\mathfrak{S}_n)$. Furthermore, for any $π\in\mathfrak{S}_n$, we establish that the number of Boolean intervals in $W(\mathfrak{S}_n)$ with minimal element $π$ is a product of Fibonacci numbers. We conclude with some directions for further study. △ Less

Submitted 27 May, 2024; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: 17 pages, 4 figures. To appear in Journal of Combinatorics

MSC Class: 05A05; 06A07; 05A05; 05A15; 05A19

Lucky Cars and the Quicksort Algorithm

Authors: Pamela E. Harris, Jan Kretschmann, J. Carlos Martínez Mori

Abstract: Quicksort is a classical divide-and-conquer sorting algorithm. It is a comparison sort that makes an average of $2(n+1)H_n - 4n$ comparisons on an array of size $n$ ordered uniformly at random, where $H_n = \sum_{i=1}^n\frac{1}{i}$ is the $n$th harmonic number. Therefore, it makes $n!\left[2(n+1)H_n - 4n\right]$ comparisons to sort all possible orderings of the array. In this article, we prove tha… ▽ More Quicksort is a classical divide-and-conquer sorting algorithm. It is a comparison sort that makes an average of $2(n+1)H_n - 4n$ comparisons on an array of size $n$ ordered uniformly at random, where $H_n = \sum_{i=1}^n\frac{1}{i}$ is the $n$th harmonic number. Therefore, it makes $n!\left[2(n+1)H_n - 4n\right]$ comparisons to sort all possible orderings of the array. In this article, we prove that this count also enumerates the parking preference lists of $n$ cars parking on a one-way street with $n$ parking spots resulting in exactly $n-1$ lucky cars (i.e., cars that park in their preferred spot). For $n\geq 2$, both counts satisfy the second order recurrence relation $ f_n=2nf_{n-1}-n(n-1)f_{n-2}+2(n-1)! $ with $f_0=f_1=0$. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 8 pages, and 2 figures, to appear in The American Mathematical Monthly

MSC Class: Primary 05A05; Secondary 68P10

Flattened Stirling Permutations

Authors: Adam Buck, Jennifer Elder, Azia A. Figueroa, Pamela E. Harris, Kimberly Harry, Anthony Simpson

Abstract: Recall that a Stirling permutation is a permutation on the multiset $\{1,1,2,2,\ldots,n,n\}$ such that any numbers appearing between repeated values of $i$ must be greater than $i$. We call a Stirling permutation ``flattened'' if the leading terms of maximal chains of ascents (called runs) are in weakly increasing order. Our main result establishes a bijection between flattened Stirling permutatio… ▽ More Recall that a Stirling permutation is a permutation on the multiset $\{1,1,2,2,\ldots,n,n\}$ such that any numbers appearing between repeated values of $i$ must be greater than $i$. We call a Stirling permutation ``flattened'' if the leading terms of maximal chains of ascents (called runs) are in weakly increasing order. Our main result establishes a bijection between flattened Stirling permutations and type $B$ set partitions of $\{0,\pm1,\pm2,\ldots,\pm (n-1)\}$, which are known to be enumerated by the Dowling numbers, and we give an independent proof of this fact. We also determine the maximal number of runs for any flattened Stirling permutation, and we enumerate flattened Stirling permutations with a small number of runs or with two runs of equal length. We conclude with some conjectures and generalizations worthy of future investigation. △ Less

Submitted 27 November, 2023; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: 17 pages, 1 figure, 2 table

MSC Class: Primary 05A05; Secondary 05A10; 05A15; 05A17; 05A18

Identifying and Extracting Rare Disease Phenotypes with Large Language Models

Authors: Cathy Shyr, Yan Hu, Paul A. Harris, Hua Xu

Abstract: Rare diseases (RDs) are collectively common and affect 300 million people worldwide. Accurate phenotyping is critical for informing diagnosis and treatment, but RD phenotypes are often embedded in unstructured text and time-consuming to extract manually. While natural language processing (NLP) models can perform named entity recognition (NER) to automate extraction, a major bottleneck is the devel… ▽ More Rare diseases (RDs) are collectively common and affect 300 million people worldwide. Accurate phenotyping is critical for informing diagnosis and treatment, but RD phenotypes are often embedded in unstructured text and time-consuming to extract manually. While natural language processing (NLP) models can perform named entity recognition (NER) to automate extraction, a major bottleneck is the development of a large, annotated corpus for model training. Recently, prompt learning emerged as an NLP paradigm that can lead to more generalizable results without any (zero-shot) or few labeled samples (few-shot). Despite growing interest in ChatGPT, a revolutionary large language model capable of following complex human prompts and generating high-quality responses, none have studied its NER performance for RDs in the zero- and few-shot settings. To this end, we engineered novel prompts aimed at extracting RD phenotypes and, to the best of our knowledge, are the first the establish a benchmark for evaluating ChatGPT's performance in these settings. We compared its performance to the traditional fine-tuning approach and conducted an in-depth error analysis. Overall, fine-tuning BioClinicalBERT resulted in higher performance (F1 of 0.689) than ChatGPT (F1 of 0.472 and 0.591 in the zero- and few-shot settings, respectively). Despite this, ChatGPT achieved similar or higher accuracy for certain entities (i.e., rare diseases and signs) in the one-shot setting (F1 of 0.776 and 0.725). This suggests that with appropriate prompt engineering, ChatGPT has the potential to match or outperform fine-tuned language models for certain entity types with just one labeled sample. While the proliferation of large language models may provide opportunities for supporting RD diagnosis and treatment, researchers and clinicians should critically evaluate model outputs and be well-informed of their limitations. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Journal ref: J Healthc Inform Res 8, 438-461 (2024)

Mesas of Stirling permutations

Authors: Nicolle González, Pamela E. Harris, Gordon Rojas Kirby, Mariana Smit Vega Garcia, Bridget Eileen Tenner

Abstract: Given a Stirling permutation w, we introduce the mesa set of w as the natural generalization of the pinnacle set of a permutation. Our main results characterize admissible mesa sets and give closed enumerative formulas in terms of rational Catalan numbers by providing an explicit bijection between mesa sets and rational Dyck paths. Given a Stirling permutation w, we introduce the mesa set of w as the natural generalization of the pinnacle set of a permutation. Our main results characterize admissible mesa sets and give closed enumerative formulas in terms of rational Catalan numbers by providing an explicit bijection between mesa sets and rational Dyck paths. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 11 pages, 3 figures

MSC Class: 05A05 (Primary); 05A15 (Secondary)

Demonstration of Machine Learning-assisted real-time noise regression in gravitational wave detectors

Authors: Muhammed Saleem, Alec Gunny, Chia-Jui Chou, Li-Cheng Yang, Shu-Wei Yeh, Andy H. Y. Chen, Ryan Magee, William Benoit, Tri Nguyen, Pinchen Fan, Deep Chatterjee, Ethan Marx, Eric Moreno, Rafia Omer, Ryan Raikman, Dylan Rankin, Ritwik Sharma, Michael Coughlin, Philip Harris, Erik Katsavounidis

Abstract: Real-time noise regression algorithms are crucial for maximizing the science outcomes of the LIGO, Virgo, and KAGRA gravitational-wave detectors. This includes improvements in the detectability, source localization and pre-merger detectability of signals thereby enabling rapid multi-messenger follow-up. In this paper, we demonstrate the effectiveness of \textit{DeepClean}, a convolutional neural n… ▽ More Real-time noise regression algorithms are crucial for maximizing the science outcomes of the LIGO, Virgo, and KAGRA gravitational-wave detectors. This includes improvements in the detectability, source localization and pre-merger detectability of signals thereby enabling rapid multi-messenger follow-up. In this paper, we demonstrate the effectiveness of \textit{DeepClean}, a convolutional neural network architecture that uses witness sensors to estimate and subtract non-linear and non-stationary noise from gravitational-wave strain data. Our study uses LIGO data from the third observing run with injected compact binary signals. As a demonstration, we use \textit{DeepClean} to subtract the noise at 60 Hz due to the power mains and their sidebands arising from non-linear coupling with other instrumental noise sources. Our parameter estimation study on the injected signals shows that \textit{DeepClean} does not do any harm to the underlying astrophysical signals in the data while it can enhances the signal-to-noise ratio of potential signals. We show that \textit{DeepClean} can be used for low-latency noise regression to produce cleaned output data at latencies $\sim 1-2$\, s. We also discuss various considerations that may be made while training \textit{DeepClean} for low latency applications. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Unit-Interval Parking Functions and the Permutohedron

Authors: Lucas Chaves Meyles, Pamela E. Harris, Richter Jordaan, Gordon Rojas Kirby, Sam Sehayek, Ethan Spingarn

Abstract: Unit-interval parking functions are subset of parking functions in which cars park at most one spot away from their preferred parking spot. In this paper, we characterize unit-interval parking functions by understanding how they decompose into prime parking functions and count unit-interval parking functions when exactly $k<n$ cars do not park in their preference. This count yields an alternate pr… ▽ More Unit-interval parking functions are subset of parking functions in which cars park at most one spot away from their preferred parking spot. In this paper, we characterize unit-interval parking functions by understanding how they decompose into prime parking functions and count unit-interval parking functions when exactly $k<n$ cars do not park in their preference. This count yields an alternate proof of a result of Hadaway and Harris establishing that unit-interval parking functions are enumerated by the Fubini numbers. Then, our main result, establishes that for all integers $0\leq k<n$, the unit-interval parking functions of length $n$ with displacement $k$ are in bijection with the $k$-dimensional faces of the permutohedron of order $n$. We conclude with some consequences of this result. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 15 pages, 5 figures

A linearization for stable and fast geographically weighted Poisson regression

Authors: Daisuke Murakami, Narumasa Tsutsumida, Takahiro Yoshida, Tomoki Nakaya, Binbin Lu, Paul Harris

Abstract: Although geographically weighted Poisson regression (GWPR) is a popular regression for spatially indexed count data, its development is relatively limited compared to that found for linear geographically weighted regression (GWR), where many extensions (e.g., multiscale GWR, scalable GWR) have been proposed. The weak development of GWPR can be attributed to the computational cost and identificatio… ▽ More Although geographically weighted Poisson regression (GWPR) is a popular regression for spatially indexed count data, its development is relatively limited compared to that found for linear geographically weighted regression (GWR), where many extensions (e.g., multiscale GWR, scalable GWR) have been proposed. The weak development of GWPR can be attributed to the computational cost and identification problem in the underpinning Poisson regression model. This study proposes linearized GWPR (L-GWPR) by introducing a log-linear approximation into the GWPR model to overcome these bottlenecks. Because the L-GWPR model is identical to the Gaussian GWR model, it is free from the identification problem, easily implemented, computationally efficient, and offers similar potential for extension. Specifically, L-GWPR does not require a double-loop algorithm, which makes GWPR slow for large samples. Furthermore, we extended L-GWPR by introducing ridge regularization to enhance its stability (regularized L-GWPR). The results of the Monte Carlo experiments confirmed that regularized L-GWPR estimates local coefficients accurately and computationally efficiently. Finally, we compared GWPR and regularized L-GWPR through a crime analysis in Tokyo. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Symbolic Regression on FPGAs for Fast Machine Learning Inference

Authors: Ho Fung Tsoi, Adrian Alan Pol, Vladimir Loncar, Ekaterina Govorkova, Miles Cranmer, Sridhara Dasu, Peter Elmer, Philip Harris, Isobel Ojalvo, Maurizio Pierini

Abstract: The high-energy physics community is investigating the potential of deploying machine-learning-based solutions on Field-Programmable Gate Arrays (FPGAs) to enhance physics sensitivity while still meeting data processing time constraints. In this contribution, we introduce a novel end-to-end procedure that utilizes a machine learning technique called symbolic regression (SR). It searches the equati… ▽ More The high-energy physics community is investigating the potential of deploying machine-learning-based solutions on Field-Programmable Gate Arrays (FPGAs) to enhance physics sensitivity while still meeting data processing time constraints. In this contribution, we introduce a novel end-to-end procedure that utilizes a machine learning technique called symbolic regression (SR). It searches the equation space to discover algebraic relations approximating a dataset. We use PySR (a software to uncover these expressions based on an evolutionary algorithm) and extend the functionality of hls4ml (a package for machine learning inference in FPGAs) to support PySR-generated expressions for resource-constrained production environments. Deep learning models often optimize the top metric by pinning the network size because the vast hyperparameter space prevents an extensive search for neural architecture. Conversely, SR selects a set of models on the Pareto front, which allows for optimizing the performance-resource trade-off directly. By embedding symbolic forms, our implementation can dramatically reduce the computational resources needed to perform critical tasks. We validate our method on a physics benchmark: the multiclass classification of jets produced in simulated proton-proton collisions at the CERN Large Hadron Collider. We show that our approach can approximate a 3-layer neural network using an inference model that achieves up to a 13-fold decrease in execution time, down to 5 ns, while still preserving more than 90% approximation accuracy. △ Less

Submitted 17 January, 2024; v1 submitted 6 May, 2023; originally announced May 2023.

Comments: 9 pages. Accepted to 26th International Conference on Computing in High Energy & Nuclear Physics (CHEP 2023)

Journal ref: EPJ Web of Conferences 295, 09036 (2024)

Feebly Interacting Particles: FIPs 2022 workshop report

Authors: C. Antel, M. Battaglieri, J. Beacham, C. Boehm, O. Buchmüller, F. Calore, P. Carenza, B. Chauhan, P. Cladè, P. Coloma, P. Crivelli, V. Dandoy, L. Darmé, B. Dey, F. F. Deppisch, A. De Roeck, M. Drewes, B. Echenard, V. V. Flambaum, P. Foldenauer, C. Gatti, M. Giannotti, A. Golutvin, M. C. Gonzalez-Garcia, S. Gori , et al. (53 additional authors not shown)

Abstract: Particle physics today faces the challenge of explaining the mystery of dark matter, the origin of matter over anti-matter in the Universe, the origin of the neutrino masses, the apparent fine-tuning of the electro-weak scale, and many other aspects of fundamental physics. Perhaps the most striking frontier to emerge in the search for answers involves new physics at mass scales comparable to famil… ▽ More Particle physics today faces the challenge of explaining the mystery of dark matter, the origin of matter over anti-matter in the Universe, the origin of the neutrino masses, the apparent fine-tuning of the electro-weak scale, and many other aspects of fundamental physics. Perhaps the most striking frontier to emerge in the search for answers involves new physics at mass scales comparable to familiar matter, below the GeV-scale, or even radically below, down to sub-eV scales, and with very feeble interaction strength. New theoretical ideas to address dark matter and other fundamental questions predict such feebly interacting particles (FIPs) at these scales, and indeed, existing data provide numerous hints for such possibility. A vibrant experimental program to discover such physics is under way, guided by a systematic theoretical approach firmly grounded on the underlying principles of the Standard Model. This document represents the report of the FIPs 2022 workshop, held at CERN between the 17 and 21 October 2022 and aims to give an overview of these efforts, their motivations, and the decadal goals that animate the community involved in the search for FIPs. △ Less

Submitted 2 May, 2023; originally announced May 2023.

Comments: 383 pp, 176 figures

Report number: CERN-TH-2023-061 DESY-23-050 FERMILAB-PUB-23-149-PPD INFN-23-14-LNF JLAB-PHY-23-3789 LA-UR-23-21432 MITP-23-015

First Constraints on the Photon Coupling of Axion-like Particles from Multimessenger Studies of the Neutron Star Merger GW170817

Authors: P. S. Bhupal Dev, Jean-François Fortin, Steven P. Harris, Kuver Sinha, Yongchao Zhang

Abstract: We use multimessenger observations of the neutron star merger event GW170817 to derive new constraints on axion-like particles (ALPs) coupling to photons. ALPs are produced via Primakoff and photon coalescence processes in the merger, escape the remnant and decay back into two photons, giving rise to a photon signal approximately along the line-of-sight to the merger. We analyze the spectral and t… ▽ More We use multimessenger observations of the neutron star merger event GW170817 to derive new constraints on axion-like particles (ALPs) coupling to photons. ALPs are produced via Primakoff and photon coalescence processes in the merger, escape the remnant and decay back into two photons, giving rise to a photon signal approximately along the line-of-sight to the merger. We analyze the spectral and temporal information of the ALP-induced photon signal, and use the Fermi-LAT observations of GW170817 to derive our new ALP constraints. We also show the improved prospects with future MeV gamma-ray missions, taking the spectral and temporal coverage of Fermi-LAT as an example. △ Less

Submitted 27 February, 2024; v1 submitted 1 May, 2023; originally announced May 2023.

Comments: 8+8 pages, 4+5 figures, some updates and minor changes, version to appear in PRL

Report number: INT-PUB-23-014

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

ECG Feature Importance Rankings: Cardiologists vs. Algorithms

Authors: Temesgen Mehari, Ashish Sundar, Alen Bosnjakovic, Peter Harris, Steven E. Williams, Axel Loewe, Olaf Doessel, Claudia Nagel, Nils Strodthoff, Philip J. Aston

Abstract: Feature importance methods promise to provide a ranking of features according to importance for a given classification task. A wide range of methods exist but their rankings often disagree and they are inherently difficult to evaluate due to a lack of ground truth beyond synthetic datasets. In this work, we put feature importance methods to the test on real-world data in the domain of cardiology,… ▽ More Feature importance methods promise to provide a ranking of features according to importance for a given classification task. A wide range of methods exist but their rankings often disagree and they are inherently difficult to evaluate due to a lack of ground truth beyond synthetic datasets. In this work, we put feature importance methods to the test on real-world data in the domain of cardiology, where we try to distinguish three specific pathologies from healthy subjects based on ECG features comparing to features used in cardiologists' decision rules as ground truth. Some methods generally performed well and others performed poorly, while some methods did well on some but not all of the problems considered. △ Less

Submitted 5 April, 2023; originally announced April 2023.

Search for an interaction mediated by axion-like particles with ultracold neutrons at the PSI

Authors: N. J. Ayres, G. Bison, K. Bodek, V. Bondar, T. Bouillaud, E. Chanel, P. -J. Chiu, B. Clement, C. B. Crawford, M. Daum, C. B. Doorenbos, S. Emmenegger, M. Fertl, P. Flaux, W. C. Griffith, P. G. Harris, N. Hild, M. Kasprzak, K. Kirch, V. Kletzl, P. A. Koss, J. Krempel, B. Lauss, T. Lefort, P. Mohanmurthy , et al. (22 additional authors not shown)

Abstract: We report on a search for a new, short-range, spin-dependent interaction using a modified version of the experimental apparatus used to measure the permanent neutron electric dipole moment at the Paul Scherrer Institute. This interaction, which could be mediated by axion-like particles, concerned the unpolarized nucleons (protons and neutrons) near the material surfaces of the apparatus and polari… ▽ More We report on a search for a new, short-range, spin-dependent interaction using a modified version of the experimental apparatus used to measure the permanent neutron electric dipole moment at the Paul Scherrer Institute. This interaction, which could be mediated by axion-like particles, concerned the unpolarized nucleons (protons and neutrons) near the material surfaces of the apparatus and polarized ultracold neutrons stored in vacuum. The dominant systematic uncertainty resulting from magnetic-field gradients was controlled to an unprecedented level of approximately 4 pT/cm using an array of optically-pumped cesium vapor magnetometers and magnetic-field maps independently recorded using a dedicated measurement device. No signature of a theoretically predicted new interaction was found, and we set a new limit on the product of the scalar and the pseudoscalar couplings $g_sg_pλ^2 < 8.3 \times 10^{-28}\,\text{m}^2$ (95% C.L.) in a range of $5\,μ\text{m} < λ< 25\,\text{mm}$ for the monopole-dipole interaction. This new result confirms and improves our previous limit by a factor of 2.7 and provides the current tightest limit obtained with free neutrons. △ Less

Submitted 31 March, 2023; originally announced March 2023.

Dense Nuclear Matter Equation of State from Heavy-Ion Collisions

Authors: Agnieszka Sorensen, Kshitij Agarwal, Kyle W. Brown, Zbigniew Chajęcki, Paweł Danielewicz, Christian Drischler, Stefano Gandolfi, Jeremy W. Holt, Matthias Kaminski, Che-Ming Ko, Rohit Kumar, Bao-An Li, William G. Lynch, Alan B. McIntosh, William G. Newton, Scott Pratt, Oleh Savchuk, Maria Stefaniak, Ingo Tews, ManYee Betty Tsang, Ramona Vogt, Hermann Wolter, Hanna Zbroszczyk, Navid Abbasi, Jörg Aichelin , et al. (111 additional authors not shown)

Abstract: The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of mu… ▽ More The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. This white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS. △ Less

Submitted 25 January, 2024; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: White paper prepared for the 2023 Long Range Plan. v3: Updated version as published in Progress in Particle and Nuclear Physics. Note: the published version does not include the executive summary; in the updated arXiv version, the executive summary is included as an appendix. v4: Corrected list of authors

Report number: INT-PUB-23-001, LA-UR-23-20514, LLNL-TR-844629

Journal ref: Prog.Part.Nucl.Phys. 134 (2024) 104080

Counting Parking Sequences and Parking Assortments Through Permutations

Authors: Spencer J. Franks, Pamela E. Harris, Kimberly Harry, Jan Kretschmann, Megan Vance

Abstract: Parking sequences (a generalization of parking functions) are defined by specifying car lengths and requiring that a car attempts to park in the first available spot after its preference. If it does not fit there, then a collision occurs and the car fails to park. In contrast, parking assortments generalize parking sequences (and parking functions) by allowing cars (also of assorted lengths) to se… ▽ More Parking sequences (a generalization of parking functions) are defined by specifying car lengths and requiring that a car attempts to park in the first available spot after its preference. If it does not fit there, then a collision occurs and the car fails to park. In contrast, parking assortments generalize parking sequences (and parking functions) by allowing cars (also of assorted lengths) to seek forward from their preference to identify a set of contiguous unoccupied spots in which they fit. We consider both parking sequences and parking assortments and establish that the number of preferences resulting in a fixed parking order $σ$ is related to the lengths of cars indexed by certain subsequences in $σ$. The sum of these numbers over all parking orders (i.e. permutations of $[n]$) yields new formulas for the total number of parking sequences and of parking assortments. △ Less

Submitted 25 January, 2023; originally announced January 2023.

Comments: 9 pages, 2 figures, 1 table

MSC Class: 05A15; 05A19

Accelerating Machine Learning Inference with GPUs in ProtoDUNE Data Processing

Authors: Tejin Cai, Kenneth Herner, Tingjun Yang, Michael Wang, Maria Acosta Flechas, Philip Harris, Burt Holzman, Kevin Pedro, Nhan Tran

Abstract: We study the performance of a cloud-based GPU-accelerated inference server to speed up event reconstruction in neutrino data batch jobs. Using detector data from the ProtoDUNE experiment and employing the standard DUNE grid job submission tools, we attempt to reprocess the data by running several thousand concurrent grid jobs, a rate we expect to be typical of current and future neutrino physics e… ▽ More We study the performance of a cloud-based GPU-accelerated inference server to speed up event reconstruction in neutrino data batch jobs. Using detector data from the ProtoDUNE experiment and employing the standard DUNE grid job submission tools, we attempt to reprocess the data by running several thousand concurrent grid jobs, a rate we expect to be typical of current and future neutrino physics experiments. We process most of the dataset with the GPU version of our processing algorithm and the remainder with the CPU version for timing comparisons. We find that a 100-GPU cloud-based server is able to easily meet the processing demand, and that using the GPU version of the event processing algorithm is two times faster than processing these data with the CPU version when comparing to the newest CPUs in our sample. The amount of data transferred to the inference server during the GPU runs can overwhelm even the highest-bandwidth network switches, however, unless care is taken to observe network facility limits or otherwise distribute the jobs to multiple sites. We discuss the lessons learned from this processing campaign and several avenues for future improvements. △ Less

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

Comments: 13 pages, 9 figures, matches accepted version

Report number: FERMILAB-PUB-22-944-ND-PPD-SCD

Journal ref: Comput Softw Big Sci 7, 11 (2023)

Pinnacle sets of signed permutations

Authors: Nicolle González, Pamela E. Harris, Gordon Rojas Kirby, Mariana Smit Vega Garcia, Bridget Eileen Tenner

Abstract: Pinnacle sets record the values of the local maxima for a given family of permutations. They were introduced by Davis-Nelson-Petersen-Tenner as a dual concept to that of peaks, previously defined by Billey-Burdzy-Sagan. In recent years pinnacles and admissible pinnacles sets for the type $A$ symmetric group have been widely studied. In this article we define the pinnacle set of signed permutations… ▽ More Pinnacle sets record the values of the local maxima for a given family of permutations. They were introduced by Davis-Nelson-Petersen-Tenner as a dual concept to that of peaks, previously defined by Billey-Burdzy-Sagan. In recent years pinnacles and admissible pinnacles sets for the type $A$ symmetric group have been widely studied. In this article we define the pinnacle set of signed permutations of types $B$ and $D$. We give a closed formula for the number of type $B$/$D$ admissible pinnacle sets and answer several other related enumerative questions. △ Less

Submitted 24 March, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: 15 pages, 3 figures, to appear in Discrete Mathematics

MSC Class: 05A05 (Primary); 05A15 (Secondary)