We demonstrate an all-optical, minimally invasive method for electron beam (e-beam) characterization using Rydberg electrometry. The e-beam passes through a dilute Rb vapor prepared in a quantum superposition of ground and Rydberg states that reduces resonant absorption in a narrow spectral region. Imaging the modifications of Rb fluorescence due to shifts in the Rydberg state from the e-beam electric field allows us to reconstruct e-beam width, centroid position, and current. We experimentally demonstrate this technique using a 20 keV e-beam in the range of currents down to 20 $\mu$A, and discuss technical challenges produced by environmental electric potentials in the detection chamber. Overall, we demonstrate the promising potential of such an approach as a minimally invasive diagnostic for charged particle beams.

We show that numerical computations based on tensor renormalization group (TRG) methods can be significantly accelerated with PyTorch on graphics processing units (GPUs) by leveraging NVIDIA's Compute Unified Device Architecture (CUDA). We find improvement in the runtime and its scaling with bond dimension for two-dimensional systems. Our results establish that the utilization of GPU resources is essential for future precision computations with TRG.

We review lattice results related to pion, kaon, $D$-meson, $B$-meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor $f_+(0)$ arising in the semileptonic $K \to \pi$ transition at zero momentum transfer, as well as the decay constant ratio $f_K/f_\pi$ and its consequences for the CKM matrix elements $V_{us}$ and $V_{ud}$. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of $SU(2)_L\times SU(2)_R$ and $SU(3)_L\times SU(3)_R$ Chiral Perturbation Theory. We review the determination of the $B_K$ parameter of neutral kaon mixing as well as the additional four $B$ parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for $m_c$ and $m_b$ as well as those for the decay constants, form factors, and mixing parameters of charmed and bottom mesons and baryons. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. We review the status of lattice determinations of the strong coupling constant $\alpha_s$. We consider nucleon matrix elements, and review the determinations of the axial, scalar and tensor bilinears, both isovector and flavor diagonal. Finally, in this review we have added a new section reviewing determinations of scale-setting quantities.

Accelerating cavities are an integral part of the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Laboratory. When any of the over 400 cavities in CEBAF experiences a fault, it disrupts beam delivery to experimental user halls. In this study, we propose the use of a deep learning model to predict slowly developing cavity faults. By utilizing pre-fault signals, we train a LSTM-CNN binary classifier to distinguish between radio-frequency (RF) signals during normal operation and RF signals indicative of impending faults. We optimize the model by adjusting the fault confidence threshold and implementing a multiple consecutive window criterion to identify fault events, ensuring a low false positive rate. Results obtained from analysis of a real dataset collected from the accelerating cavities simulating a deployed scenario demonstrate the model's ability to identify normal signals with 99.99% accuracy and correctly predict 80% of slowly developing faults. Notably, these achievements were achieved in the context of a highly imbalanced dataset, and fault predictions were made several hundred milliseconds before the onset of the fault. Anticipating faults enables preemptive measures to improve operational efficiency by preventing or mitigating their occurrence.

The circuit complexity for Hamiltonian simulation of the sparsified SYK model with $N$ Majorana fermions and $q=4$ (quartic interactions) which retains holographic features (referred to as `minimal holographic sparsified SYK') with $k\ll N^{3}/24$ (where $k$ is the total number of interaction terms times 1/$N$) using second-order Trotter method and Jordan-Wigner encoding is found to be $\widetilde{\mathcal{O}}(k^{p}N^{3/2} \log N (\mathcal{J}t)^{3/2}\varepsilon^{-1/2})$ where $t$ is the simulation time, $\varepsilon$ is the desired error in the implementation of the unitary $U = \exp(-iHt)$, $\mathcal{J}$ is the disorder strength, and p < 1. This complexity implies that with less than a hundred logical qubits and about $10^{6}$ gates, it will be possible to achieve an advantage in this model and simulate real-time dynamics up to scrambling time.

One of the portals to new physics is a light scalar coupled to the Standard Model (SM) Higgs. In this paper we focus on hadronic decays of such a scalar in the regime where QCD dynamics is nonperturbative, resulting, e.g., in decays to pairs of pions or kaons, while also allowing for scalar couplings to the SM fermions to deviate from the Higgs-mixed light scalar limit. Representations of the corresponding form factors can be obtained using dispersive techniques, however, several sources of uncertainty affect the final results. We reexamine these decays, paying special attention to the quantification of uncertainties. For the light Higgs-mixed scalar scenario, we compare our results with previous works. For a general set of couplings of the light scalar to Standard Model fields, we provide a public code, {\tt hipsofcobra}, to compute the decay widths.

We present a quantum computational framework for SU(2) lattice gauge theory, leveraging continuous variables instead of discrete qubits to represent the infinite-dimensional Hilbert space of the gauge fields. We consider a ladder as well as a two-dimensional grid of plaquettes, detailing the use of gauge fixing to reduce the degrees of freedom and simplify the Hamiltonian. We demonstrate how the system dynamics, ground states, and energy gaps can be computed using the continuous-variable approach to quantum computing. Our results indicate that it is feasible to study non-Abelian gauge theories with continuous variables, providing new avenues for understanding the real-time dynamics of quantum field theories.

Integrals involving derivatives of Legendre polynomials frequently arise in applications ranging from multipole expansions for processes involving electromagnetic probes to spectral methods in numerical physics. Despite their practical relevance, closed-form expressions for such integrals - particularly involving arbitrary derivative orders - are not readily accessible in standard references or symbolic tools. In this note, we derive and present general analytic expressions for integrals of the form $\int_{-1}^{+1} dx P^{(q)}_{n} (x) P^{(k)}_{m} (x)$, where $P_{n} (x)$ and $P_{m} (x)$ are Legendre polynomials and $q$, $k$ denote their order of differentiation. Using repeated integration by parts, parity arguments, and closed-form boundary evaluations, we obtain explicit binomial and Gamma-function representations valid for all non-negative integers $n$, $m$, $q$, $k$. These results unify and extend known orthogonality relations and provide ready-to-use tools for analytic and computational contexts.

We use scanning superconducting quantum interference device (SQUID) microscopy to image vortices in superconducting strips fabricated from NbTiN thin films. We repeatedly cool superconducting strips with different width in an applied magnetic field and image the individual vortices. From these images we determine the threshold field at which the first vortex enters a strip, as well as the number and spatial configuration of vortices beyond this threshold field. We model vortex behavior with and without considering the effect of pinning by numercially minimizing the Gibbs free energy of vortices in the strips. Our measurements provide a first benchmark to understand the flux trapping properties of NbTiN thin films directly relevant to NbTiN-based superconducting circuits and devices.

We investigate the relativistic scattering of three identical scalar bosons interacting via pair-wise interactions. Extending techniques from the non-relativistic three-body scattering theory, we provide a detailed and general prescription for solving and analytically continuing integral equations describing the three-body reactions. We use these techniques to study a system with zero angular momenta described by a single scattering length leading to a bound state in a two-body sub-channel. We obtain bound-state--particle and three-particle amplitudes in the previously unexplored kinematical regime; in particular, for real energies below elastic thresholds and complex energies in the physical and unphysical Riemann sheets. We extract positions of three-particle bound-states that agree with previous finite-volume studies, providing further evidence for the consistency of the relativistic finite-volume three-body quantization conditions. We also determine previously unobserved virtual bound states in this theory. Finally, we find numerical evidence of the breakdown of the two-body finite-volume formalism in the vicinity of the left-hand cuts and argue for the generalization of the existing formalism.

We present our results on the electromagnetic form factor of pion over a wide range of $Q^2$ using lattice QCD simulations with Wilson-clover valence quarks and HISQ sea quarks. We study the form factor at the physical point with a lattice spacing $a=0.076$ fm. To study the lattice spacing and quark mass effects, we also present results for 300 MeV pion at two different lattice spacings $a=0.04$ and 0.06 fm. The lattice calculations at the physical quark mass appear to agree with the experimental results. Through fits to the form factor, we estimate the charge radius of pion for physical pion mass to be $\langle r_{\pi}^2 \rangle=0.42(2)~{\rm fm}^2$.

Dark sectors, consisting of new, light, weakly-coupled particles that do not interact with the known strong, weak, or electromagnetic forces, are a particularly compelling possibility for new physics. Nature may contain numerous dark sectors, each with their own beautiful structure, distinct particles, and forces. This review summarizes the physics motivation for dark sectors and the exciting opportunities for experimental exploration. It is the summary of the Intensity Frontier subgroup "New, Light, Weakly-coupled Particles" of the Community Summer Study 2013 (Snowmass). We discuss axions, which solve the strong CP problem and are an excellent dark matter candidate, and their generalization to axion-like particles. We also review dark photons and other dark-sector particles, including sub-GeV dark matter, which are theoretically natural, provide for dark matter candidates or new dark matter interactions, and could resolve outstanding puzzles in particle and astro-particle physics. In many cases, the exploration of dark sectors can proceed with existing facilities and comparatively modest experiments. A rich, diverse, and low-cost experimental program has been identified that has the potential for one or more game-changing discoveries. These physics opportunities should be vigorously pursued in the US and elsewhere.

In this paper we perform an amplitude analysis of essentially all published pion and kaon pair production data from two photon collisions below 1.5 GeV. This includes all the high statistics results from Belle, as well as older data from Mark II at SLAC, CELLO at DESY, Crystal Ball at SLAC. The purpose of this analysis is to provide as close to a model-independent determination of the $\gamma\gamma$ to meson pair amplitudes as possible. Having data with limited angular coverage, typically |\cos \theta| < 0.6-0.8, and no polarization information for reactions in which spin is an essential complication, the determination of the underlying amplitudes might appear an intractable problem. However, imposing the basic constraints required by analyticity, unitarity, and crossing-symmetry makes up for the experimentally missing information. Final state interactions among the meson pairs are critical to this analysis. To fix these, we include the latest $\pi\pi\to\pi\pi$, ${\overline K}K$ scattering amplitudes given by dispersive analyses, supplemented in the ${\overline K}K$ threshold region by the recent precision Dalitz plot analysis from BaBar. With these hadronic amplitudes built into unitarity, we can constrain the overall description of $\gamma\gamma\to\pi\pi$ and $\overline{K}K$ datasets, both integrated and differential cross-sections, including the high statistics charged and neutral pion, as well as $K_sK_s$ data from Belle. Since this analysis invokes coupled hadronic channels, having data on both $\gamma\gamma\to\pi\pi$ and $\overline{K}K$ reduces the solution space to essentially a single form. We present the partial wave amplitudes, show how well they fit all the available data, and give the two photon couplings of scalar and tensor resonances that appear. These partial waves are important inputs into forthcoming dispersive calculations of hadronic light-by-light scattering.

The performance of superconducting radio-frequency (SRF) cavities made of Niobium is tied to the quality of their inner surfaces exposed to the radio frequency (RF) waves. Future superconducting particle accelerators, because of their dimensions or the unprecedented stringent technical requirements, require the development of innovative surface processing techniques to improve processing reliability and if possible ecological footprint and cost, compared to conventional chemical processes. Metallographic polishing (MP) has emerged as a promising polishing technology to address these challenges. Previous studies focused on the characterization of the processed material surface at room temperature in the absence of RF waves. However, the evaluation of material properties, such as surface resistance under RF, at cryogenic temperature has failed, primarily due to the unavailability of devices capable of achieving the necessary resolution in the nanohm range. To overcome this limitation, a quadrupole resonator (QPR) has been utilized. The RF results demonstrate that the MP polishing, developed to preserve a high-quality niobium surface with very low surface resistance, is highly effective compared to conventional polishing. This conclusion is further supported by topography and microstructural analysis of the QPR top-hat samples, which revealed the clear superiority of the metallographic approach.

Electropolishing is the premier surface preparation method for high-Q, high-gradient superconducting RF cavities made of Nb. This leaves behind an apparently smooth surface, yet the achievable peak magnetic fields fall well below the superheating field of Nb, in most cases. In this work, the ultimate surface finish of electropolishing was investigated by studying its effect on highly polished Nb samples. Electropolishing introduces high slope angle sloped-steps at grain boundaries. The magnetic field enhancement and superheating field suppression factors associated with such a geometry are calculated in the London theory. Despite the by-eye smoothness of electropolished Nb, such defects compromise the stability of the low-loss Meissner state, likely limiting the achievable peak accelerating fields in superconducting RF cavities. Finally, the impact of surface roughness on impurity diffusion is investigated which can link surface roughness to the effectiveness of heat treatments like low-temperature baking or nitrogen infusion in the vortex nucleation or hydride hypotheses. Surface roughness tends to decrease the effective dose of impurities as a result of the expansion of impurities into regions with greater internal angle. The effective dose of impurities can be protected by minimizing slope angles and step heights, ensuring uniformity.

Particle accelerator operation requires simultaneous optimization of multiple objectives. Multi-Objective Optimization (MOO) is particularly challenging due to trade-offs between the objectives. Evolutionary algorithms, such as genetic algorithm (GA), have been leveraged for many optimization problems, however, they do not apply to complex control problems by design. This paper demonstrates the power of differentiability for solving MOO problems using a Deep Differentiable Reinforcement Learning (DDRL) algorithm in particle accelerators. We compare DDRL algorithm with Model Free Reinforcement Learning (MFRL), GA and Bayesian Optimization (BO) for simultaneous optimization of heat load and trip rates in the Continuous Electron Beam Accelerator Facility (CEBAF). The underlying problem enforces strict constraints on both individual states and actions as well as cumulative (global) constraint for energy requirements of the beam. A physics-based surrogate model based on real data is developed. This surrogate model is differentiable and allows back-propagation of gradients. The results are evaluated in the form of a Pareto-front for two objectives. We show that the DDRL outperforms MFRL, BO, and GA on high dimensional problems.

We explore the feasibility of gate-based hybrid quantum computing using both discrete (qubit) and continuous (qumode) variables on trapped-ion platforms. Trapped-ion systems have demonstrated record one- and two-qubit gate fidelities and long qubit coherence times, while qumodes, which can be represented by the collective vibrational modes of the ion chain, have remained relatively unexplored for their use in computing. Using numerical simulations, we show that high-fidelity hybrid gates and measurement operations can be achieved for existing trapped-ion quantum platforms. As an exemplary application, we consider quantum simulations of the Jaynes-Cummings-Hubbard model, which is given by a one-dimensional chain of interacting spin and boson degrees of freedom. Using classical simulations, we study its real-time evolution and develop a suitable variational quantum algorithm for ground state preparation. Our results motivate further studies of hybrid quantum computing in this context, which may lead to direct applications in condensed matter and fundamental particle and nuclear physics.

Spatially non-local matrix elements are useful lattice-QCD observables in a variety of contexts, for example in determining hadron structure. To quote credible estimates of the systematic uncertainties in these calculations, one must understand, among other things, the size of the finite-volume effects when such matrix elements are extracted from numerical lattice calculations. In this work, we estimate finite-volume effects for matrix elements of non-local operators, composed of two currents displaced in a spatial direction by a distance $\xi$. We find that the finite-volume corrections depend on the details of the matrix element. If the external state is the lightest degree of freedom in the theory, e.g.~the pion in QCD, then the volume corrections scale as $e^{-m_\pi (L- \xi)}$, where $m_\pi$ is the mass of the light state. For heavier external states the usual $e^{- m_\pi L}$ form is recovered, but with a polynomial prefactor of the form $L^m/|L - \xi|^n$ that can lead to enhanced volume effects. These observations are potentially relevant to a wide variety of observables being studied using lattice QCD, including parton distribution functions, double-beta-decay and Compton-scattering matrix elements, and long-range weak matrix elements.

Measurements of SIDIS multiplicities for $\pi^+$ and $\pi^-$ from proton and deuteron targets are reported on a grid of hadron kinematic variables $z$, $P_{T}$, and $\phi^{*}$ for leptonic kinematic variables in the range $0.3