Showing 1–50 of 73 results for author: Turner, W

New measurements of the Lyman-$α$ forest continuum and effective optical depth with LyCAN and DESI Y1 data

Authors: Wynne Turner, Paul Martini, Naim Göksel Karaçaylı, J. Aguilar, S. Ahlen, D. Brooks, T. Claybaugh, A. de la Macorra, A. Dey, P. Doel, K. Fanning, J. E. Forero-Romero, S. Gontcho A Gontcho, A. X. Gonzalez-Morales, G. Gutierrez, J. Guy, H. K. Herrera-Alcantar, K. Honscheid, S. Juneau, T. Kisner, A. Kremin, A. Lambert, M. Landriau, L. Le Guillou, A. Meisner , et al. (20 additional authors not shown)

Abstract: We present the Lyman-$α$ Continuum Analysis Network (LyCAN), a Convolutional Neural Network that predicts the unabsorbed quasar continuum within the rest-frame wavelength range of $1040-1600$ Angstroms based on the red side of the Lyman-$α$ emission line ($1216-1600$ Angstroms). We developed synthetic spectra based on a Gaussian Mixture Model representation of Nonnegative Matrix Factorization (NMF… ▽ More We present the Lyman-$α$ Continuum Analysis Network (LyCAN), a Convolutional Neural Network that predicts the unabsorbed quasar continuum within the rest-frame wavelength range of $1040-1600$ Angstroms based on the red side of the Lyman-$α$ emission line ($1216-1600$ Angstroms). We developed synthetic spectra based on a Gaussian Mixture Model representation of Nonnegative Matrix Factorization (NMF) coefficients. These coefficients were derived from high-resolution, low-redshift ($z<0.2$) Hubble Space Telescope/Cosmic Origins Spectrograph quasar spectra. We supplemented this COS-based synthetic sample with an equal number of DESI Year 5 mock spectra. LyCAN performs extremely well on testing sets, achieving a median error in the forest region of 1.5% on the DESI mock sample, 2.0% on the COS-based synthetic sample, and 4.1% on the original COS spectra. LyCAN outperforms Principal Component Analysis (PCA)- and NMF-based prediction methods using the same training set by a factor of two or more. We predict the intrinsic continua of 83,635 DESI Year 1 spectra in the redshift range of $2.1 \leq z \leq 4.2$ and perform an absolute measurement of the evolution of the effective optical depth. This is the largest sample employed to measure the optical depth evolution to date. We fit a power-law of the form $τ(z) = τ_0 (1+z)^γ$ to our measurements and find $τ_0 = (2.46 \pm 0.14)\times10^{-3}$ and $γ= 3.62 \pm 0.04$. Our results show particular agreement with high-resolution, ground-based observations around $z = 2$, indicating that LyCAN is able to predict the quasar continuum in the forest region with only spectral information outside the forest. △ Less

Submitted 10 May, 2024; originally announced May 2024.

Comments: 23 pages, 15 figures, 3 tables; submitted to ApJ

Towards a Stallings-type theorem for finite groups

Authors: Johannes Carmesin, George Kontogeorgiou, Jan Kurkofka, Will J. Turner

Abstract: A recent development in graph theory is to study local separators, vertex sets that need not separate graphs globally but just locally. We use this idea to conjecture an extension of Stallings' theorem to finite nilpotent groups. We provide an example demonstrating that the assumption of nilpotency is necessary, and prove a stronger form of this conjecture for low connectivities, as follows. A f… ▽ More A recent development in graph theory is to study local separators, vertex sets that need not separate graphs globally but just locally. We use this idea to conjecture an extension of Stallings' theorem to finite nilpotent groups. We provide an example demonstrating that the assumption of nilpotency is necessary, and prove a stronger form of this conjecture for low connectivities, as follows. A finite nilpotent group that has a Cayley graph with a local separator of size at most two is cyclic, dihedral or a direct product of a cyclic group with the group of order two. △ Less

Submitted 12 March, 2024; originally announced March 2024.

Comments: 33 pages, 18 figures

MSC Class: 05C40 (Primary); 05C83; 20F65; 05E18; 20E34; 20F18 (Secondary)

Detailed Report on the Measurement of the Positive Muon Anomalous Magnetic Moment to 0.20 ppm

Authors: D. P. Aguillard, T. Albahri, D. Allspach, A. Anisenkov, K. Badgley, S. Baeßler, I. Bailey, L. Bailey, V. A. Baranov, E. Barlas-Yucel, T. Barrett, E. Barzi, F. Bedeschi, M. Berz, M. Bhattacharya, H. P. Binney, P. Bloom, J. Bono, E. Bottalico, T. Bowcock, S. Braun, M. Bressler, G. Cantatore, R. M. Carey, B. C. K. Casey , et al. (168 additional authors not shown)

Abstract: We present details on a new measurement of the muon magnetic anomaly, $a_μ= (g_μ-2)/2$. The result is based on positive muon data taken at Fermilab's Muon Campus during the 2019 and 2020 accelerator runs. The measurement uses $3.1$ GeV$/c$ polarized muons stored in a $7.1$-m-radius storage ring with a $1.45$ T uniform magnetic field. The value of $ a_μ$ is determined from the measured difference b… ▽ More We present details on a new measurement of the muon magnetic anomaly, $a_μ= (g_μ-2)/2$. The result is based on positive muon data taken at Fermilab's Muon Campus during the 2019 and 2020 accelerator runs. The measurement uses $3.1$ GeV$/c$ polarized muons stored in a $7.1$-m-radius storage ring with a $1.45$ T uniform magnetic field. The value of $ a_μ$ is determined from the measured difference between the muon spin precession frequency and its cyclotron frequency. This difference is normalized to the strength of the magnetic field, measured using Nuclear Magnetic Resonance (NMR). The ratio is then corrected for small contributions from beam motion, beam dispersion, and transient magnetic fields. We measure $a_μ= 116 592 057 (25) \times 10^{-11}$ (0.21 ppm). This is the world's most precise measurement of this quantity and represents a factor of $2.2$ improvement over our previous result based on the 2018 dataset. In combination, the two datasets yield $a_μ(\text{FNAL}) = 116 592 055 (24) \times 10^{-11}$ (0.20 ppm). Combining this with the measurements from Brookhaven National Laboratory for both positive and negative muons, the new world average is $a_μ$(exp) $ = 116 592 059 (22) \times 10^{-11}$ (0.19 ppm). △ Less

Submitted 22 May, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: 48 pages, 29 figures; 4 pages of Supplement Material; version accepted for publication in Physical Review D

Report number: FERMILAB-PUB-24-0084-AD-CSAID-PPD

Free probability, path developments and signature kernels as universal scaling limits

Authors: Thomas Cass, William F. Turner

Abstract: Random developments of a path into a matrix Lie group $G_N$ have recently been used to construct signature-based kernels on path space. Two examples include developments into GL$(N;\mathbb{R})$ and $U(N;\mathbb{C})$, the general linear and unitary groups of dimension $N$. For the former, [MLS23] showed that the signature kernel is obtained via a scaling limit of developments with Gaussian vector f… ▽ More Random developments of a path into a matrix Lie group $G_N$ have recently been used to construct signature-based kernels on path space. Two examples include developments into GL$(N;\mathbb{R})$ and $U(N;\mathbb{C})$, the general linear and unitary groups of dimension $N$. For the former, [MLS23] showed that the signature kernel is obtained via a scaling limit of developments with Gaussian vector fields. The second instance was used in [LLN23] to construct a metric between probability measures on path space. We present a unified treatment to obtaining large $N$ limits by leveraging the tools of free probability theory. An important conclusion is that the limiting kernels, while dependent on the choice of Lie group, are nonetheless universal limits with respect to how the development map is randomised. For unitary developments, the limiting kernel is given by the contraction of a signature against the monomials of freely independent semicircular random variables. Using the Schwinger-Dyson equations, we show that this kernel can be obtained by solving a novel quadratic functional equation. We provide a convergent numerical scheme for this equation, together with rates, which does not require computation of signatures themselves. △ Less

Submitted 19 February, 2024; originally announced February 2024.

MSC Class: 60L10 (Primary)

New constraints on ultraheavy dark matter from the LZ experiment

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

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

Submitted 13 February, 2024; originally announced February 2024.

Comments: 9 pages, 7 figures

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

Integration of 4D BIM and Robot Task Planning: Creation and Flow of Construction-Related Information for Action-Level Simulation of Indoor Wall Frame Installation

Authors: Hafiz Oyediran, William Turner, Kyungki Kim, Matthew Barrows

Abstract: An obstacle toward construction robotization is the lack of methods to plan robot operations within the entire construction planning process. Despite the strength in modeling construction site conditions, 4D BIM technologies cannot perform construction robot task planning considering the contexts of given work environments. To address this limitation, this study presents a framework that integrate… ▽ More An obstacle toward construction robotization is the lack of methods to plan robot operations within the entire construction planning process. Despite the strength in modeling construction site conditions, 4D BIM technologies cannot perform construction robot task planning considering the contexts of given work environments. To address this limitation, this study presents a framework that integrates 4D BIM and robot task planning, presents an information flow for the integration, and performs high-level robot task planning and detailed simulation. The framework uniquely incorporates a construction robot knowledge base that derives robot-related modeling requirements to augment a 4D BIM model. Then, the 4D BIM model is converted into a robot simulation world where a robot performs a sequence of actions retrieving construction-related information. A case study focusing on the interior wall frame installation demonstrates the potential of systematic integration in achieving context-aware robot task planning and simulation in construction environments. △ Less

Submitted 5 February, 2024; originally announced February 2024.

First Constraints on WIMP-Nucleon Effective Field Theory Couplings in an Extended Energy Region From LUX-ZEPLIN

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

Abstract: Following the first science results of the LUX-ZEPLIN (LZ) experiment, a dual-phase xenon time projection chamber operating from the Sanford Underground Research Facility in Lead, South Dakota, USA, we report the initial limits on a model-independent non-relativistic effective field theory describing the complete set of possible interactions of a weakly interacting massive particle (WIMP) with a n… ▽ More Following the first science results of the LUX-ZEPLIN (LZ) experiment, a dual-phase xenon time projection chamber operating from the Sanford Underground Research Facility in Lead, South Dakota, USA, we report the initial limits on a model-independent non-relativistic effective field theory describing the complete set of possible interactions of a weakly interacting massive particle (WIMP) with a nucleon. These results utilize the same 5.5 t fiducial mass and 60 live days of exposure collected for the LZ spin-independent and spin-dependent analyses while extending the upper limit of the energy region of interest by a factor of 7.5 to 270 keVnr. No significant excess in this high energy region is observed. Using a profile-likelihood ratio analysis, we report 90% confidence level exclusion limits on the coupling of each individual non-relativistic WIMP-nucleon operator for both elastic and inelastic interactions in the isoscalar and isovector bases. △ Less

Submitted 26 February, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 17 pages 11 figures

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

ChipNeMo: Domain-Adapted LLMs for Chip Design

Authors: Mingjie Liu, Teodor-Dumitru Ene, Robert Kirby, Chris Cheng, Nathaniel Pinckney, Rongjian Liang, Jonah Alben, Himyanshu Anand, Sanmitra Banerjee, Ismet Bayraktaroglu, Bonita Bhaskaran, Bryan Catanzaro, Arjun Chaudhuri, Sharon Clay, Bill Dally, Laura Dang, Parikshit Deshpande, Siddhanth Dhodhi, Sameer Halepete, Eric Hill, Jiashang Hu, Sumit Jain, Ankit Jindal, Brucek Khailany, George Kokai , et al. (17 additional authors not shown)

Abstract: ChipNeMo aims to explore the applications of large language models (LLMs) for industrial chip design. Instead of directly deploying off-the-shelf commercial or open-source LLMs, we instead adopt the following domain adaptation techniques: domain-adaptive tokenization, domain-adaptive continued pretraining, model alignment with domain-specific instructions, and domain-adapted retrieval models. We e… ▽ More ChipNeMo aims to explore the applications of large language models (LLMs) for industrial chip design. Instead of directly deploying off-the-shelf commercial or open-source LLMs, we instead adopt the following domain adaptation techniques: domain-adaptive tokenization, domain-adaptive continued pretraining, model alignment with domain-specific instructions, and domain-adapted retrieval models. We evaluate these methods on three selected LLM applications for chip design: an engineering assistant chatbot, EDA script generation, and bug summarization and analysis. Our evaluations demonstrate that domain-adaptive pretraining of language models, can lead to superior performance in domain related downstream tasks compared to their base LLaMA2 counterparts, without degradations in generic capabilities. In particular, our largest model, ChipNeMo-70B, outperforms the highly capable GPT-4 on two of our use cases, namely engineering assistant chatbot and EDA scripts generation, while exhibiting competitive performance on bug summarization and analysis. These results underscore the potential of domain-specific customization for enhancing the effectiveness of large language models in specialized applications. △ Less

Submitted 4 April, 2024; v1 submitted 31 October, 2023; originally announced November 2023.

Comments: Updated results for ChipNeMo-70B model

Data fusion for a multi-scale model of a wheat leaf surface: a unifying approach using a radial basis function partition of unity method

Authors: Riley M. Whebell, Timothy J. Moroney, Ian W. Turner, Ravindra Pethiyagoda, Marie-Luise Wille, Justin J. Cooper-White, Arvind Kumar, Philip Taylor, Scott W. McCue

Abstract: Realistic digital models of plant leaves are crucial to fluid dynamics simulations of droplets for optimising agrochemical spray technologies. The presence and nature of small features (on the order of 100$\mathrm{μm}$) such as ridges and hairs on the surface have been shown to significantly affect the droplet evaporation, and thus the leaf's potential uptake of active ingredients. We show that th… ▽ More Realistic digital models of plant leaves are crucial to fluid dynamics simulations of droplets for optimising agrochemical spray technologies. The presence and nature of small features (on the order of 100$\mathrm{μm}$) such as ridges and hairs on the surface have been shown to significantly affect the droplet evaporation, and thus the leaf's potential uptake of active ingredients. We show that these microstructures can be captured by implicit radial basis function partition of unity (RBFPU) surface reconstructions from micro-CT scan datasets. However, scanning a whole leaf ($20\mathrm{cm^2}$) at micron resolutions is infeasible due to both extremely large data storage requirements and scanner time constraints. Instead, we micro-CT scan only a small segment of a wheat leaf ($4\mathrm{mm^2}$). We fit a RBFPU implicit surface to this segment, and an explicit RBFPU surface to a lower resolution laser scan of the whole leaf. Parameterising the leaf using a locally orthogonal coordinate system, we then replicate the now resolved microstructure many times across a larger, coarser, representation of the leaf surface that captures important macroscale features, such as its size, shape, and orientation. The edge of one segment of the microstructure model is blended into its neighbour naturally by the partition of unity method. The result is one implicit surface reconstruction that captures the wheat leaf's features at both the micro- and macro-scales. △ Less

Submitted 17 September, 2023; originally announced September 2023.

Comments: 23 pages, 11 figures

MSC Class: 65D05 (Primary); 65D18; 65Z05; 68U05 (Secondary)

Measurement of the Positive Muon Anomalous Magnetic Moment to 0.20 ppm

Authors: D. P. Aguillard, T. Albahri, D. Allspach, A. Anisenkov, K. Badgley, S. Baeßler, I. Bailey, L. Bailey, V. A. Baranov, E. Barlas-Yucel, T. Barrett, E. Barzi, F. Bedeschi, M. Berz, M. Bhattacharya, H. P. Binney, P. Bloom, J. Bono, E. Bottalico, T. Bowcock, S. Braun, M. Bressler, G. Cantatore, R. M. Carey, B. C. K. Casey , et al. (166 additional authors not shown)

Abstract: We present a new measurement of the positive muon magnetic anomaly, $a_μ\equiv (g_μ- 2)/2$, from the Fermilab Muon $g\!-\!2$ Experiment using data collected in 2019 and 2020. We have analyzed more than 4 times the number of positrons from muon decay than in our previous result from 2018 data. The systematic error is reduced by more than a factor of 2 due to better running conditions, a more stable… ▽ More We present a new measurement of the positive muon magnetic anomaly, $a_μ\equiv (g_μ- 2)/2$, from the Fermilab Muon $g\!-\!2$ Experiment using data collected in 2019 and 2020. We have analyzed more than 4 times the number of positrons from muon decay than in our previous result from 2018 data. The systematic error is reduced by more than a factor of 2 due to better running conditions, a more stable beam, and improved knowledge of the magnetic field weighted by the muon distribution, $\tildeω'^{}_p$, and of the anomalous precession frequency corrected for beam dynamics effects, $ω_a$. From the ratio $ω_a / \tildeω'^{}_p$, together with precisely determined external parameters, we determine $a_μ= 116\,592\,057(25) \times 10^{-11}$ (0.21 ppm). Combining this result with our previous result from the 2018 data, we obtain $a_μ\text{(FNAL)} = 116\,592\,055(24) \times 10^{-11}$ (0.20 ppm). The new experimental world average is $a_μ(\text{Exp}) = 116\,592\,059(22)\times 10^{-11}$ (0.19 ppm), which represents a factor of 2 improvement in precision. △ Less

Submitted 4 October, 2023; v1 submitted 11 August, 2023; originally announced August 2023.

Comments: 8 pages, 3 figures

Report number: FERMILAB-PUB-23-385-AD-CSAID-PPD

Journal ref: Phys. Rev. Lett. 131, 161802 (2023)

Generative Modelling of Lévy Area for High Order SDE Simulation

Authors: Andraž Jelinčič, Jiajie Tao, William F. Turner, Thomas Cass, James Foster, Hao Ni

Abstract: It is well known that, when numerically simulating solutions to SDEs, achieving a strong convergence rate better than O(\sqrt{h}) (where h is the step size) requires the use of certain iterated integrals of Brownian motion, commonly referred to as its "Lévy areas". However, these stochastic integrals are difficult to simulate due to their non-Gaussian nature and for a d-dimensional Brownian motion… ▽ More It is well known that, when numerically simulating solutions to SDEs, achieving a strong convergence rate better than O(\sqrt{h}) (where h is the step size) requires the use of certain iterated integrals of Brownian motion, commonly referred to as its "Lévy areas". However, these stochastic integrals are difficult to simulate due to their non-Gaussian nature and for a d-dimensional Brownian motion with d > 2, no fast almost-exact sampling algorithm is known. In this paper, we propose LévyGAN, a deep-learning-based model for generating approximate samples of Lévy area conditional on a Brownian increment. Due to our "Bridge-flipping" operation, the output samples match all joint and conditional odd moments exactly. Our generator employs a tailored GNN-inspired architecture, which enforces the correct dependency structure between the output distribution and the conditioning variable. Furthermore, we incorporate a mathematically principled characteristic-function based discriminator. Lastly, we introduce a novel training mechanism termed "Chen-training", which circumvents the need for expensive-to-generate training data-sets. This new training procedure is underpinned by our two main theoretical results. For 4-dimensional Brownian motion, we show that LévyGAN exhibits state-of-the-art performance across several metrics which measure both the joint and marginal distributions. We conclude with a numerical experiment on the log-Heston model, a popular SDE in mathematical finance, demonstrating that high-quality synthetic Lévy area can lead to high order weak convergence and variance reduction when using multilevel Monte Carlo (MLMC). △ Less

Submitted 4 August, 2023; originally announced August 2023.

MSC Class: 65C30

A search for new physics in low-energy electron recoils from the first LZ exposure

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

Abstract: The LUX-ZEPLIN (LZ) experiment is a dark matter detector centered on a dual-phase xenon time projection chamber. We report searches for new physics appearing through few-keV-scale electron recoils, using the experiment's first exposure of 60 live days and a fiducial mass of 5.5t. The data are found to be consistent with a background-only hypothesis, and limits are set on models for new physics inc… ▽ More The LUX-ZEPLIN (LZ) experiment is a dark matter detector centered on a dual-phase xenon time projection chamber. We report searches for new physics appearing through few-keV-scale electron recoils, using the experiment's first exposure of 60 live days and a fiducial mass of 5.5t. The data are found to be consistent with a background-only hypothesis, and limits are set on models for new physics including solar axion electron coupling, solar neutrino magnetic moment and millicharge, and electron couplings to galactic axion-like particles and hidden photons. Similar limits are set on weakly interacting massive particle (WIMP) dark matter producing signals through ionized atomic states from the Migdal effect. △ Less

Submitted 9 September, 2023; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: 13 pages, 10 figures. See https://tinyurl.com/LZDataReleaseRun1ER for a data release related to this paper

Journal ref: Phys. Rev. D 108, 072006 (2023)

Measurement of the forward $η$ meson production rate in p-p collisions at $\sqrt{s}$=13 TeV with the LHCf-Arm2 detector

Authors: Giuseppe Piparo, Oscar Adriani, Eugenio Berti, Pietro Betti, Lorenzo Bonechi, Massimo Bongi, Raffaello D'Alessandro, Sebastiano Detti, Maurice Haguenauer, Yoshitaka Itow, Katsuaki Kasahara, Yuga Kitagami, Moe Kondo, Yutaka Matsubara, Hiroaki Menjo, Yasushi Muraki, Ken Ohashi, Paolo Papini, Sergio Ricciarini, Takashi Sako, Nobuyuki Sakurai, Monica Scaringella, Yuki Shimizu, Tadashi Tamura, Alessio Tiberio , et al. (4 additional authors not shown)

Abstract: The forward $η$ mesons production has been observed by the Large Hadron Collider forward (LHCf) experiment in proton-proton collision at $\sqrt{s}$=13 TeV. This paper presents the measurement of the inclusive production rate of $η$ in $p_T<$ 1.1 GeV/c, expressed as a function of the Feynman-x variable. These results are compared with the predictions of several hadronic interaction models commonly… ▽ More The forward $η$ mesons production has been observed by the Large Hadron Collider forward (LHCf) experiment in proton-proton collision at $\sqrt{s}$=13 TeV. This paper presents the measurement of the inclusive production rate of $η$ in $p_T<$ 1.1 GeV/c, expressed as a function of the Feynman-x variable. These results are compared with the predictions of several hadronic interaction models commonly used for the modelling of the air showers produced by ultra-high energy cosmic rays. This is both the first measurement of $η$ mesons from LHCf and the first time a particle containing strange quarks has been observed in the forward region for high-energy collisions. These results will provide a powerful constraint on hadronic interaction models for the purpose of improving the understanding of the processes underlying the air showers produced in the Earth's atmosphere by ultra-energetic cosmic rays. △ Less

Submitted 11 May, 2023; originally announced May 2023.

Comments: 18 pages, 7 figures, 4 tables

Report number: CERN-EP-2023-076

Journal ref: J. High Energ. Phys. 2023, 169 (2023)

Background Determination for the LUX-ZEPLIN (LZ) Dark Matter Experiment

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

Abstract: The LUX-ZEPLIN experiment recently reported limits on WIMP-nucleus interactions from its initial science run, down to $9.2\times10^{-48}$ cm$^2$ for the spin-independent interaction of a 36 GeV/c$^2$ WIMP at 90% confidence level. In this paper, we present a comprehensive analysis of the backgrounds important for this result and for other upcoming physics analyses, including neutrinoless double-bet… ▽ More The LUX-ZEPLIN experiment recently reported limits on WIMP-nucleus interactions from its initial science run, down to $9.2\times10^{-48}$ cm$^2$ for the spin-independent interaction of a 36 GeV/c$^2$ WIMP at 90% confidence level. In this paper, we present a comprehensive analysis of the backgrounds important for this result and for other upcoming physics analyses, including neutrinoless double-beta decay searches and effective field theory interpretations of LUX-ZEPLIN data. We confirm that the in-situ determinations of bulk and fixed radioactive backgrounds are consistent with expectations from the ex-situ assays. The observed background rate after WIMP search criteria were applied was $(6.3\pm0.5)\times10^{-5}$ events/keV$_{ee}$/kg/day in the low-energy region, approximately 60 times lower than the equivalent rate reported by the LUX experiment. △ Less

Submitted 17 July, 2023; v1 submitted 30 November, 2022; originally announced November 2022.

Comments: 25 pages, 15 figures

Journal ref: Phys. Rev. D 108, 012010 (2023)

Energetic electron precipitation driven by electromagnetic ion cyclotron waves from ELFIN's low altitude perspective

Authors: V. Angelopoulos, X. -J. Zhang, A. V. Artemyev, D. Mourenas, E. Tsai, C. Wilkins, A. Runov, J. Liu, D. L. Turner, W. Li, K. Khurana, R. E. Wirz, V. A. Sergeev, X. Meng, J. Wu, M. D. Hartinger, T. Raita, Y. Shen, X. An, X. Shi, M. F. Bashir, X. Shen, L. Gan, M. Qin, L. Capannolo , et al. (61 additional authors not shown)

Abstract: We review comprehensive observations of electromagnetic ion cyclotron (EMIC) wave-driven energetic electron precipitation using data from the energetic electron detector on the Electron Losses and Fields InvestigatioN (ELFIN) mission, two polar-orbiting low-altitude spinning CubeSats, measuring 50-5000 keV electrons with good pitch-angle and energy resolution. EMIC wave-driven precipitation exhibi… ▽ More We review comprehensive observations of electromagnetic ion cyclotron (EMIC) wave-driven energetic electron precipitation using data from the energetic electron detector on the Electron Losses and Fields InvestigatioN (ELFIN) mission, two polar-orbiting low-altitude spinning CubeSats, measuring 50-5000 keV electrons with good pitch-angle and energy resolution. EMIC wave-driven precipitation exhibits a distinct signature in energy-spectrograms of the precipitating-to-trapped flux ratio: peaks at 0.5 MeV which are abrupt (bursty) with significant substructure (occasionally down to sub-second timescale). Multiple ELFIN passes over the same MLT sector allow us to study the spatial and temporal evolution of the EMIC wave - electron interaction region. Using two years of ELFIN data, we assemble a statistical database of 50 events of strong EMIC wave-driven precipitation. Most reside at L=5-7 at dusk, while a smaller subset exists at L=8-12 at post-midnight. The energies of the peak-precipitation ratio and of the half-peak precipitation ratio (our proxy for the minimum resonance energy) exhibit an L-shell dependence in good agreement with theoretical estimates based on prior statistical observations of EMIC wave power spectra. The precipitation ratio's spectral shape for the most intense events has an exponential falloff away from the peak (i.e., on either side of 1.45 MeV). It too agrees well with quasi-linear diffusion theory based on prior statistics of wave spectra. Sub-MeV electron precipitation observed concurrently with strong EMIC wave-driven 1MeV precipitation has a spectral shape that is consistent with efficient pitch-angle scattering down to 200-300 keV by much less intense higher frequency EMIC waves. These results confirm the critical role of EMIC waves in driving relativistic electron losses. Nonlinear effects may abound and require further investigation. △ Less

Submitted 28 November, 2022; originally announced November 2022.

TAG: Learning Circuit Spatial Embedding From Layouts

Authors: Keren Zhu, Hao Chen, Walker J. Turner, George F. Kokai, Po-Hsuan Wei, David Z. Pan, Haoxing Ren

Abstract: Analog and mixed-signal (AMS) circuit designs still rely on human design expertise. Machine learning has been assisting circuit design automation by replacing human experience with artificial intelligence. This paper presents TAG, a new paradigm of learning the circuit representation from layouts leveraging text, self-attention and graph. The embedding network model learns spatial information with… ▽ More Analog and mixed-signal (AMS) circuit designs still rely on human design expertise. Machine learning has been assisting circuit design automation by replacing human experience with artificial intelligence. This paper presents TAG, a new paradigm of learning the circuit representation from layouts leveraging text, self-attention and graph. The embedding network model learns spatial information without manual labeling. We introduce text embedding and a self-attention mechanism to AMS circuit learning. Experimental results demonstrate the ability to predict layout distances between instances with industrial FinFET technology benchmarks. The effectiveness of the circuit representation is verified by showing the transferability to three other learning tasks with limited data in the case studies: layout matching prediction, wirelength estimation, and net parasitic capacitance prediction. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Comments: Accepted by ICCAD 2022

Burgers' equation in the complex plane

Authors: Daniel J. VandenHeuvel, Christopher J. Lustri, John R. King, Ian W. Turner, Scott W. McCue

Abstract: Burgers' equation is a well-studied model in applied mathematics with connections to the Navier-Stokes equations in one spatial direction and traffic flow, for example. Following on from previous work, we analyse solutions to Burgers' equation in the complex plane, concentrating on the dynamics of the complex singularities and their relationship to the solution on the real line. For an initial con… ▽ More Burgers' equation is a well-studied model in applied mathematics with connections to the Navier-Stokes equations in one spatial direction and traffic flow, for example. Following on from previous work, we analyse solutions to Burgers' equation in the complex plane, concentrating on the dynamics of the complex singularities and their relationship to the solution on the real line. For an initial condition with a simple pole in each of the upper- and lower-half planes, we apply formal asymptotics in the small- and large-time limits in order to characterise the initial and later motion of the singularities. The small-time limit highlights how infinitely many singularities are born at $t=0$ and how they orientate themselves to lie increasingly close to anti-Stokes lines in the far-field of the inner problem. This inner problem also reveals whether or not the closest singularity to the real axis moves toward the axis or away. For intermediate times, we use the exact solution, apply method of steepest descents, and implement the AAA approximation to track the complex singularities. Connections are made between the motion of the closest singularity to the real axis and the steepness of the solution on the real line. While Burgers' equation has an exact solution, we deliberately apply a mix of techniques in our analysis in an attempt to develop methodology that can be applied to other nonlinear partial differential equations that do not. △ Less

Submitted 26 February, 2023; v1 submitted 11 August, 2022; originally announced August 2022.

Efficient inference and identifiability analysis for differential equation models with random parameters

Authors: Alexander P. Browning, Christopher Drovandi, Ian W. Turner, Adrianne L. Jenner, Matthew J. Simpson

Abstract: Heterogeneity is a dominant factor in the behaviour of many biological processes. Despite this, it is common for mathematical and statistical analyses to ignore biological heterogeneity as a source of variability in experimental data. Therefore, methods for exploring the identifiability of models that explicitly incorporate heterogeneity through variability in model parameters are relatively under… ▽ More Heterogeneity is a dominant factor in the behaviour of many biological processes. Despite this, it is common for mathematical and statistical analyses to ignore biological heterogeneity as a source of variability in experimental data. Therefore, methods for exploring the identifiability of models that explicitly incorporate heterogeneity through variability in model parameters are relatively underdeveloped. We develop a new likelihood-based framework, based on moment matching, for inference and identifiability analysis of differential equation models that capture biological heterogeneity through parameters that vary according to probability distributions. As our novel method is based on an approximate likelihood function, it is highly flexible; we demonstrate identifiability analysis using both a frequentist approach based on profile likelihood, and a Bayesian approach based on Markov-chain Monte Carlo. Through three case studies, we demonstrate our method by providing a didactic guide to inference and identifiability analysis of hyperparameters that relate to the statistical moments of model parameters from independent observed data. Our approach has a computational cost comparable to analysis of models that neglect heterogeneity, a significant improvement over many existing alternatives. We demonstrate how analysis of random parameter models can aid better understanding of the sources of heterogeneity from biological data. △ Less

Submitted 27 October, 2022; v1 submitted 20 July, 2022; originally announced July 2022.

Comments: Minor changes to text. Additional results in supplementary material. Additional statistics regarding results given in main and supplementary material

MSC Class: 92B99; 62F99

First Dark Matter Search Results from the LUX-ZEPLIN (LZ) Experiment

Authors: J. Aalbers, D. S. Akerib, C. W. Akerlof, A. K. Al Musalhi, F. Alder, A. Alqahtani, S. K. Alsum, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, S. Azadi, A. J. Bailey, A. Baker, J. Balajthy, S. Balashov, J. Bang, J. W. Bargemann, M. J. Barry, J. Barthel, D. Bauer, A. Baxter , et al. (322 additional authors not shown)

Abstract: The LUX-ZEPLIN experiment is a dark matter detector centered on a dual-phase xenon time projection chamber operating at the Sanford Underground Research Facility in Lead, South Dakota, USA. This Letter reports results from LUX-ZEPLIN's first search for weakly interacting massive particles (WIMPs) with an exposure of 60~live days using a fiducial mass of 5.5 t. A profile-likelihood ratio analysis s… ▽ More The LUX-ZEPLIN experiment is a dark matter detector centered on a dual-phase xenon time projection chamber operating at the Sanford Underground Research Facility in Lead, South Dakota, USA. This Letter reports results from LUX-ZEPLIN's first search for weakly interacting massive particles (WIMPs) with an exposure of 60~live days using a fiducial mass of 5.5 t. A profile-likelihood ratio analysis shows the data to be consistent with a background-only hypothesis, setting new limits on spin-independent WIMP-nucleon, spin-dependent WIMP-neutron, and spin-dependent WIMP-proton cross sections for WIMP masses above 9 GeV/c$^2$. The most stringent limit is set for spin-independent scattering at 36 GeV/c$^2$, rejecting cross sections above 9.2$\times 10^{-48}$ cm$^2$ at the 90% confidence level. △ Less

Submitted 2 August, 2023; v1 submitted 8 July, 2022; originally announced July 2022.

Comments: 9 pages, 8 figures. See https://doi.org/10.1103/PhysRevLett.131.041002 for a data release related to this paper

Journal ref: Phys. Rev. Lett. 131, 041002 (2023)

Topologies on unparameterised path space

Authors: Thomas Cass, William F. Turner

Abstract: The signature of a path, introduced by K.T. Chen [5] in $1954$, has been extensively studied in recent years. The $2010$ paper [12] of Hambly and Lyons showed that the signature is injective on the space of continuous finite-variation paths up to a general notion of reparameterisation called tree-like equivalence. The signature has been widely used in applications, underpinned by the result [15] t… ▽ More The signature of a path, introduced by K.T. Chen [5] in $1954$, has been extensively studied in recent years. The $2010$ paper [12] of Hambly and Lyons showed that the signature is injective on the space of continuous finite-variation paths up to a general notion of reparameterisation called tree-like equivalence. The signature has been widely used in applications, underpinned by the result [15] that guarantees uniform approximation of a continuous function on a compact set by a linear functional of the signature. We study in detail, and for the first time, the properties of three candidate topologies on the set of unparameterised paths (the tree-like equivalence classes). These are obtained through properties of the signature and are: (1) the product topology, obtained by equipping the tensor algebra with the product topology and requiring $S$ to be an embedding, (2) the quotient topology derived from the 1-variation topology on the underlying path space, and (3) the metric topology associated to $d( [ γ] ,[ σ] ) := \vert\vert γ^*-σ^*\vert\vert_{1}$ using suitable representatives $γ^*$ and $σ^*$ of the equivalence classes. The topologies are ordered by strict inclusion, (1) being the weakest and (3) the strongest. Each is separable and Hausdorff, (1) being both metrisable and $σ$-compact, but not a Baire space and so neither Polish nor locally compact. The quotient topology (2) is not metrisable and the metric $d$ is not complete. An important function on (unparameterised) path space is the (fixed-time) solution map of a controlled differential equation. For a broad class of such equations, we prove measurability of this map for each topology. Under stronger regularity assumptions, we show continuity on explicit compact subsets of the product topology (1). We relate these results to the expected signature model of [15]. △ Less

Submitted 29 June, 2022; v1 submitted 22 June, 2022; originally announced June 2022.

A Next-Generation Liquid Xenon Observatory for Dark Matter and Neutrino Physics

Authors: J. Aalbers, K. Abe, V. Aerne, F. Agostini, S. Ahmed Maouloud, D. S. Akerib, D. Yu. Akimov, J. Akshat, A. K. Al Musalhi, F. Alder, S. K. Alsum, L. Althueser, C. S. Amarasinghe, F. D. Amaro, A. Ames, T. J. Anderson, B. Andrieu, N. Angelides, E. Angelino, J. Angevaare, V. C. Antochi, D. Antón Martin, B. Antunovic, E. Aprile, H. M. Araújo , et al. (572 additional authors not shown)

Abstract: The nature of dark matter and properties of neutrinos are among the most pressing issues in contemporary particle physics. The dual-phase xenon time-projection chamber is the leading technology to cover the available parameter space for Weakly Interacting Massive Particles (WIMPs), while featuring extensive sensitivity to many alternative dark matter candidates. These detectors can also study neut… ▽ More The nature of dark matter and properties of neutrinos are among the most pressing issues in contemporary particle physics. The dual-phase xenon time-projection chamber is the leading technology to cover the available parameter space for Weakly Interacting Massive Particles (WIMPs), while featuring extensive sensitivity to many alternative dark matter candidates. These detectors can also study neutrinos through neutrinoless double-beta decay and through a variety of astrophysical sources. A next-generation xenon-based detector will therefore be a true multi-purpose observatory to significantly advance particle physics, nuclear physics, astrophysics, solar physics, and cosmology. This review article presents the science cases for such a detector. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 77 pages, 40 figures, 1262 references

Report number: INT-PUB-22-003

Journal ref: J. Phys. G: Nucl. Part. Phys. 50 (2023) 013001

Cosmogenic production of $^{37}$Ar in the context of the LUX-ZEPLIN experiment

Authors: J. Aalbers, D. S. Akerib, A. K. Al Musalhi, F. Alder, S. K. Alsum, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, X. Bai, A. Baker, J. Balajthy, S. Balashov, J. Bang, J. W. Bargemann, D. Bauer, A. Baxter, K. Beattie, E. P. Bernard, A. Bhatti, A. Biekert, T. P. Biesiadzinski , et al. (183 additional authors not shown)

Abstract: We estimate the amount of $^{37}$Ar produced in natural xenon via cosmic ray-induced spallation, an inevitable consequence of the transportation and storage of xenon on the Earth's surface. We then calculate the resulting $^{37}$Ar concentration in a 10-tonne payload~(similar to that of the LUX-ZEPLIN experiment) assuming a representative schedule of xenon purification, storage and delivery to the… ▽ More We estimate the amount of $^{37}$Ar produced in natural xenon via cosmic ray-induced spallation, an inevitable consequence of the transportation and storage of xenon on the Earth's surface. We then calculate the resulting $^{37}$Ar concentration in a 10-tonne payload~(similar to that of the LUX-ZEPLIN experiment) assuming a representative schedule of xenon purification, storage and delivery to the underground facility. Using the spallation model by Silberberg and Tsao, the sea level production rate of $^{37}$Ar in natural xenon is estimated to be 0.024~atoms/kg/day. Assuming the xenon is successively purified to remove radioactive contaminants in 1-tonne batches at a rate of 1~tonne/month, the average $^{37}$Ar activity after 10~tonnes are purified and transported underground is 0.058--0.090~$μ$Bq/kg, depending on the degree of argon removal during above-ground purification. Such cosmogenic $^{37}$Ar will appear as a noticeable background in the early science data, while decaying with a 35~day half-life. This newly-noticed production mechanism of $^{37}$Ar should be considered when planning for future liquid xenon-based experiments. △ Less

Submitted 22 March, 2022; v1 submitted 8 January, 2022; originally announced January 2022.

The Straw Tracking Detector for the Fermilab Muon $g-2$ Experiment

Authors: B. T. King, T. Albahri, S. Al-Kilani, D. Allspach, D. Beckner, A. Behnke, T. J. V. Bowcock, D. Boyden, R. M. Carey, J. Carroll, B. C. K. Casey, S. Charity, R. Chislett, M. Eads, A. Epps, S. B. Foster, D. Gastler, S. Grant, T. Halewood-Leagas, K. Hardin, E. Hazen, G. Hesketh, D. J. Hollywood, T. Jones, C. Kenziora , et al. (32 additional authors not shown)

Abstract: The Muon $g-2$ Experiment at Fermilab uses a gaseous straw tracking detector to make detailed measurements of the stored muon beam profile, which are essential for the experiment to achieve its uncertainty goals. Positrons from muon decays spiral inward and pass through the tracking detector before striking an electromagnetic calorimeter. The tracking detector is therefore located inside the vacuu… ▽ More The Muon $g-2$ Experiment at Fermilab uses a gaseous straw tracking detector to make detailed measurements of the stored muon beam profile, which are essential for the experiment to achieve its uncertainty goals. Positrons from muon decays spiral inward and pass through the tracking detector before striking an electromagnetic calorimeter. The tracking detector is therefore located inside the vacuum chamber in a region where the magnetic field is large and non-uniform. As such, the tracking detector must have a low leak rate to maintain a high-quality vacuum, must be non-magnetic so as not to perturb the magnetic field and, to minimize energy loss, must have a low radiation length. The performance of the tracking detector has met or surpassed the design requirements, with adequate electronic noise levels, an average straw hit resolution of $(110 \pm 20) \,μ$m, a detection efficiency of 97% or higher, and no performance degradation or signs of aging. The tracking detector's measurements result in an otherwise unachievable understanding of the muon's beam motion, particularly at early times in the experiment's measurement period when there are a significantly greater number of muons decaying. This is vital to the statistical power of the experiment, as well as facilitating the precise extraction of several systematic corrections and uncertainties. This paper describes the design, construction, testing, commissioning, and performance of the tracking detector. △ Less

Submitted 24 February, 2022; v1 submitted 3 November, 2021; originally announced November 2021.

Comments: 37 pages, 27 figures

Journal ref: 2022 JINST 17 P02035

Signature asymptotics, empirical processes, and optimal transport

Authors: Thomas Cass, William F. Turner, Remy Messadene

Abstract: Rough path theory provides one with the notion of signature, a graded family of tensors which characterise, up to a negligible equivalence class, and ordered stream of vector-valued data. In the last few years, use of the signature has gained traction in time-series analysis, machine learning , deep learning and more recently in kernel methods. In this article, we lay down the theoretical foundati… ▽ More Rough path theory provides one with the notion of signature, a graded family of tensors which characterise, up to a negligible equivalence class, and ordered stream of vector-valued data. In the last few years, use of the signature has gained traction in time-series analysis, machine learning , deep learning and more recently in kernel methods. In this article, we lay down the theoretical foundations for a connection between signature asymptotics, the theory of empirical processes, and Wasserstein distances, opening up the landscape and toolkit of the second and third in the study of the first. Our main contribution is to show that the Hambly-Lyons limit can be reinterpreted as a statement about the asymptotic behaviour of Wasserstein distances between two independent empirical measures of samples from the same underlying distribution. In the setting studied here, these measures are derived from samples from a probability distribution which is determined by geometrical properties of the underlying path. The general question of rates of convergence for these objects has been studied in depth in the recent monograph of Bobkov and Ledoux. By using these results, we generalise the original result of Hambly and Lyons from $C^3$ curves to a broad class of $C^2$ ones. We conclude by providing an explicit way to compute the limit in terms of a second-order differential equation. △ Less

Submitted 5 May, 2023; v1 submitted 23 July, 2021; originally announced July 2021.

Comments: 22 pages, 1 figure

MSC Class: 60L10

Projected sensitivity of the LUX-ZEPLIN (LZ) experiment to the two-neutrino and neutrinoless double beta decays of $^{134}$Xe

Authors: The LUX-ZEPLIN, Collaboration, :, D. S. Akerib, A. K. Al Musalhi, S. K. Alsum, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araujo, J. E. Armstrong, M. Arthurs, X. Bai, J. Balajthy, S. Balashov, J. Bang, J. W. Bargemann, D. Bauer, A. Baxter, P. Beltrame, E. P. Bernard, A. Bernstein, A. Bhatti, A. Biekert , et al. (172 additional authors not shown)

Abstract: The projected sensitivity of the LUX-ZEPLIN (LZ) experiment to two-neutrino and neutrinoless double beta decay of $^{134}$Xe is presented. LZ is a 10-tonne xenon time projection chamber optimized for the detection of dark matter particles, that is expected to start operating in 2021 at Sanford Underground Research Facility, USA. Its large mass of natural xenon provides an exceptional opportunity t… ▽ More The projected sensitivity of the LUX-ZEPLIN (LZ) experiment to two-neutrino and neutrinoless double beta decay of $^{134}$Xe is presented. LZ is a 10-tonne xenon time projection chamber optimized for the detection of dark matter particles, that is expected to start operating in 2021 at Sanford Underground Research Facility, USA. Its large mass of natural xenon provides an exceptional opportunity to search for the double beta decay of $^{134}$Xe, for which xenon detectors enriched in $^{136}$Xe are less effective. For the two-neutrino decay mode, LZ is predicted to exclude values of the half-life up to 1.7$\times$10$^{24}$ years at 90% confidence level (CL), and has a three-sigma observation potential of 8.7$\times$10$^{23}$ years, approaching the predictions of nuclear models. For the neutrinoless decay mode LZ, is projected to exclude values of the half-life up to 7.3$\times$10$^{24}$ years at 90% CL. △ Less

Submitted 22 November, 2021; v1 submitted 26 April, 2021; originally announced April 2021.

Comments: Version accepted for publication in Phys. Rev. C

Measurement of the Positive Muon Anomalous Magnetic Moment to 0.46 ppm

Authors: B. Abi, T. Albahri, S. Al-Kilani, D. Allspach, L. P. Alonzi, A. Anastasi, A. Anisenkov, F. Azfar, K. Badgley, S. Baeßler, I. Bailey, V. A. Baranov, E. Barlas-Yucel, T. Barrett, E. Barzi, A. Basti, F. Bedeschi, A. Behnke, M. Berz, M. Bhattacharya, H. P. Binney, R. Bjorkquist, P. Bloom, J. Bono, E. Bottalico , et al. (212 additional authors not shown)

Abstract: We present the first results of the Fermilab Muon g-2 Experiment for the positive muon magnetic anomaly $a_μ\equiv (g_μ-2)/2$. The anomaly is determined from the precision measurements of two angular frequencies. Intensity variation of high-energy positrons from muon decays directly encodes the difference frequency $ω_a$ between the spin-precession and cyclotron frequencies for polarized muons in… ▽ More We present the first results of the Fermilab Muon g-2 Experiment for the positive muon magnetic anomaly $a_μ\equiv (g_μ-2)/2$. The anomaly is determined from the precision measurements of two angular frequencies. Intensity variation of high-energy positrons from muon decays directly encodes the difference frequency $ω_a$ between the spin-precession and cyclotron frequencies for polarized muons in a magnetic storage ring. The storage ring magnetic field is measured using nuclear magnetic resonance probes calibrated in terms of the equivalent proton spin precession frequency ${\tildeω'^{}_p}$ in a spherical water sample at 34.7$^{\circ}$C. The ratio $ω_a / {\tildeω'^{}_p}$, together with known fundamental constants, determines $a_μ({\rm FNAL}) = 116\,592\,040(54)\times 10^{-11}$ (0.46\,ppm). The result is 3.3 standard deviations greater than the standard model prediction and is in excellent agreement with the previous Brookhaven National Laboratory (BNL) E821 measurement. After combination with previous measurements of both $μ^+$ and $μ^-$, the new experimental average of $a_μ({\rm Exp}) = 116\,592\,061(41)\times 10^{-11}$ (0.35\,ppm) increases the tension between experiment and theory to 4.2 standard deviations △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 10 pages; 4 figures

Report number: FERMILAB-PUB-21-132-E

Journal ref: Phys. Rev. Lett. 126, 141801 (2021)

Measurement of the anomalous precession frequency of the muon in the Fermilab Muon g-2 experiment

Authors: T. Albahri, A. Anastasi, A. Anisenkov, K. Badgley, S. Baeßler, I. Bailey, V. A. Baranov, E. Barlas-Yucel, T. Barrett, A. Basti, F. Bedeschi, M. Berz, M. Bhattacharya, H. P. Binney, P. Bloom, J. Bono, E. Bottalico, T. Bowcock, G. Cantatore, R. M. Carey, B. C. K. Casey, D. Cauz, R. Chakraborty, S. P. Chang, A. Chapelain , et al. (153 additional authors not shown)

Abstract: The Muon g-2 Experiment at Fermi National Accelerator Laboratory (FNAL) has measured the muon anomalous precession frequency $ω_a$ to an uncertainty of 434 parts per billion (ppb), statistical, and 56 ppb, systematic, with data collected in four storage ring configurations during its first physics run in 2018. When combined with a precision measurement of the magnetic field of the experiment's muo… ▽ More The Muon g-2 Experiment at Fermi National Accelerator Laboratory (FNAL) has measured the muon anomalous precession frequency $ω_a$ to an uncertainty of 434 parts per billion (ppb), statistical, and 56 ppb, systematic, with data collected in four storage ring configurations during its first physics run in 2018. When combined with a precision measurement of the magnetic field of the experiment's muon storage ring, the precession frequency measurement determines a muon magnetic anomaly of $a_μ({\rm FNAL}) = 116\,592\,040(54) \times 10^{-11}$ (0.46 ppm). This article describes the multiple techniques employed in the reconstruction, analysis and fitting of the data to measure the precession frequency. It also presents the averaging of the results from the eleven separate determinations of ω_a, and the systematic uncertainties on the result. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 29 pages, 19 figures. Published in Physical Review D

Report number: FERMILAB-PUB-21-183-E

Journal ref: Phys. Rev. D 103, 072002 (2021)

Beam dynamics corrections to the Run-1 measurement of the muon anomalous magnetic moment at Fermilab

Authors: T. Albahri, A. Anastasi, K. Badgley, S. Baeßler, I. Bailey, V. A. Baranov, E. Barlas-Yucel, T. Barrett, F. Bedeschi, M. Berz, M. Bhattacharya, H. P. Binney, P. Bloom, J. Bono, E. Bottalico, T. Bowcock, G. Cantatore, R. M. Carey, B. C. K. Casey, D. Cauz, R. Chakraborty, S. P. Chang, A. Chapelain, S. Charity, R. Chislett , et al. (152 additional authors not shown)

Abstract: This paper presents the beam dynamics systematic corrections and their uncertainties for the Run-1 data set of the Fermilab Muon g-2 Experiment. Two corrections to the measured muon precession frequency $ω_a^m$ are associated with well-known effects owing to the use of electrostatic quadrupole (ESQ) vertical focusing in the storage ring. An average vertically oriented motional magnetic field is fe… ▽ More This paper presents the beam dynamics systematic corrections and their uncertainties for the Run-1 data set of the Fermilab Muon g-2 Experiment. Two corrections to the measured muon precession frequency $ω_a^m$ are associated with well-known effects owing to the use of electrostatic quadrupole (ESQ) vertical focusing in the storage ring. An average vertically oriented motional magnetic field is felt by relativistic muons passing transversely through the radial electric field components created by the ESQ system. The correction depends on the stored momentum distribution and the tunes of the ring, which has relatively weak vertical focusing. Vertical betatron motions imply that the muons do not orbit the ring in a plane exactly orthogonal to the vertical magnetic field direction. A correction is necessary to account for an average pitch angle associated with their trajectories. A third small correction is necessary because muons that escape the ring during the storage time are slightly biased in initial spin phase compared to the parent distribution. Finally, because two high-voltage resistors in the ESQ network had longer than designed RC time constants, the vertical and horizontal centroids and envelopes of the stored muon beam drifted slightly, but coherently, during each storage ring fill. This led to the discovery of an important phase-acceptance relationship that requires a correction. The sum of the corrections to $ω_a^m$ is 0.50 $\pm$ 0.09 ppm; the uncertainty is small compared to the 0.43 ppm statistical precision of $ω_a^m$. △ Less

Submitted 23 April, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: 35 pages, 29 figures. Accepted by Phys. Rev. Accel. Beams

Report number: FERMILAB-PUB-21-133-E

Journal ref: Phys. Rev. Accel. Beams 24, 044002 (2021)

Magnetic Field Measurement and Analysis for the Muon g-2 Experiment at Fermilab

Authors: T. Albahri, A. Anastasi, K. Badgley, S. Baeßler, I. Bailey, V. A. Baranov, E. Barlas-Yucel, T. Barrett, F. Bedeschi, M. Berz, M. Bhattacharya, H. P. Binney, P. Bloom, J. Bono, E. Bottalico, T. Bowcock, G. Cantatore, R. M. Carey, B. C. K. Casey, D. Cauz, R. Chakraborty, S. P. Chang, A. Chapelain, S. Charity, R. Chislett , et al. (148 additional authors not shown)

Abstract: The Fermi National Accelerator Laboratory has measured the anomalous precession frequency $a^{}_μ= (g^{}_μ-2)/2$ of the muon to a combined precision of 0.46 parts per million with data collected during its first physics run in 2018. This paper documents the measurement of the magnetic field in the muon storage ring. The magnetic field is monitored by nuclear magnetic resonance systems and calibrat… ▽ More The Fermi National Accelerator Laboratory has measured the anomalous precession frequency $a^{}_μ= (g^{}_μ-2)/2$ of the muon to a combined precision of 0.46 parts per million with data collected during its first physics run in 2018. This paper documents the measurement of the magnetic field in the muon storage ring. The magnetic field is monitored by nuclear magnetic resonance systems and calibrated in terms of the equivalent proton spin precession frequency in a spherical water sample at 34.7$^\circ$C. The measured field is weighted by the muon distribution resulting in $\tildeω'^{}_p$, the denominator in the ratio $ω^{}_a$/$\tildeω'^{}_p$ that together with known fundamental constants yields $a^{}_μ$. The reported uncertainty on $\tildeω'^{}_p$ for the Run-1 data set is 114 ppb consisting of uncertainty contributions from frequency extraction, calibration, mapping, tracking, and averaging of 56 ppb, and contributions from fast transient fields of 99 ppb. △ Less

Submitted 17 June, 2022; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: Added one citation and corrected missing normalization in Eqs (35) and (36)

Report number: FERMILAB-PUB-21-109-E

Journal ref: Phys. Rev. A 103, 042208 (2021)

Projected sensitivities of the LUX-ZEPLIN (LZ) experiment to new physics via low-energy electron recoils

Authors: The LZ Collaboration, D. S. Akerib, A. K. Al Musalhi, S. K. Alsum, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, X. Bai, J. Balajthy, S. Balashov, J. Bang, J. W. Bargemann, D. Bauer, A. Baxter, P. Beltrame, E. P. Bernard, A. Bernstein, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch , et al. (172 additional authors not shown)

Abstract: LUX-ZEPLIN (LZ) is a dark matter detector expected to obtain world-leading sensitivity to weakly interacting massive particles (WIMPs) interacting via nuclear recoils with a ~7-tonne xenon target mass. This manuscript presents sensitivity projections to several low-energy signals of the complementary electron recoil signal type: 1) an effective neutrino magnetic moment and 2) an effective neutrino… ▽ More LUX-ZEPLIN (LZ) is a dark matter detector expected to obtain world-leading sensitivity to weakly interacting massive particles (WIMPs) interacting via nuclear recoils with a ~7-tonne xenon target mass. This manuscript presents sensitivity projections to several low-energy signals of the complementary electron recoil signal type: 1) an effective neutrino magnetic moment and 2) an effective neutrino millicharge, both for pp-chain solar neutrinos, 3) an axion flux generated by the Sun, 4) axion-like particles forming the galactic dark matter, 5) hidden photons, 6) mirror dark matter, and 7) leptophilic dark matter. World-leading sensitivities are expected in each case, a result of the large 5.6t 1000d exposure and low expected rate of electron recoil backgrounds in the $<$100keV energy regime. A consistent signal generation, background model and profile-likelihood analysis framework is used throughout. △ Less

Submitted 18 May, 2021; v1 submitted 23 February, 2021; originally announced February 2021.

Comments: v2 updates exclusion sensitivities from single-sided to two-sided

Optical Calibration System for the LUX-ZEPLIN (LZ) Outer Detector

Authors: W. Turner, A. Baxter, H. J. Birch, B. Boxer, S. Burdin, E. Fraser, A. Greenall, S. Powel, P. Sutcliffe

Abstract: The LUX-ZEPLIN experiment will search for dark matter particle interactions with a detector containing a total of 10 tonnes of liquid xenon. Surrounding the liquid xenon cryostat is an outer detector veto system with the primary aim of vetoing neutron single-scatter events in the liquid xenon that could mimic a weakly interacting massive particle (WIMP) dark matter signal. The outer detector consi… ▽ More The LUX-ZEPLIN experiment will search for dark matter particle interactions with a detector containing a total of 10 tonnes of liquid xenon. Surrounding the liquid xenon cryostat is an outer detector veto system with the primary aim of vetoing neutron single-scatter events in the liquid xenon that could mimic a weakly interacting massive particle (WIMP) dark matter signal. The outer detector consists of approximately 17 tonnes of gadolinium-loaded liquid scintillator confined to 10 acrylic tanks surrounding the cryostat and 228,000 litres of water as the outermost layer. It will be monitored by 120 inward-facing 8-inch photomultiplier tubes. An optical calibration system has been designed and built to calibrate and monitor these photomultiplier tubes allowing the veto system to reach its required efficiency and thus ensuring that LUX-ZEPLIN meets its target sensitivity. △ Less

Submitted 16 February, 2021; v1 submitted 11 February, 2021; originally announced February 2021.

Enhancing the sensitivity of the LUX-ZEPLIN (LZ) dark matter experiment to low energy signals

Authors: D. S. Akerib, A. K. Al Musalhi, S. K. Alsum, C. S. Amarasinghe, A. Ames, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, X. Bai, J. Balajthy, S. Balashov, J. Bang, J. W. Bargemann, D. Bauer, A. Baxter, P. Beltrame, E. P. Bernard, A. Bernstein, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, G. M. Blockinger , et al. (162 additional authors not shown)

Abstract: Two-phase xenon detectors, such as that at the core of the forthcoming LZ dark matter experiment, use photomultiplier tubes to sense the primary (S1) and secondary (S2) scintillation signals resulting from particle interactions in their liquid xenon target. This paper describes a simulation study exploring two techniques to lower the energy threshold of LZ to gain sensitivity to low-mass dark matt… ▽ More Two-phase xenon detectors, such as that at the core of the forthcoming LZ dark matter experiment, use photomultiplier tubes to sense the primary (S1) and secondary (S2) scintillation signals resulting from particle interactions in their liquid xenon target. This paper describes a simulation study exploring two techniques to lower the energy threshold of LZ to gain sensitivity to low-mass dark matter and astrophysical neutrinos, which will be applicable to other liquid xenon detectors. The energy threshold is determined by the number of detected S1 photons; typically, these must be recorded in three or more photomultiplier channels to avoid dark count coincidences that mimic real signals. To lower this threshold: a) we take advantage of the double photoelectron emission effect, whereby a single vacuum ultraviolet photon has a $\sim20\%$ probability of ejecting two photoelectrons from a photomultiplier tube photocathode; and b) we drop the requirement of an S1 signal altogether, and use only the ionization signal, which can be detected more efficiently. For both techniques we develop signal and background models for the nominal exposure, and explore accompanying systematic effects, including the dependence on the free electron lifetime in the liquid xenon. When incorporating double photoelectron signals, we predict a factor of $\sim 4$ sensitivity improvement to the dark matter-nucleon scattering cross-section at $2.5$ GeV/c$^2$, and a factor of $\sim1.6$ increase in the solar $^8$B neutrino detection rate. Dropping the S1 requirement may allow sensitivity gains of two orders of magnitude in both cases. Finally, we apply these techniques to even lower masses by taking into account the atomic Migdal effect; this could lower the dark matter particle mass threshold to $80$ MeV/c$^2$. △ Less

Submitted 21 January, 2021; originally announced January 2021.

Comments: 14 pages, 6 figures

Homogenisation for the monodomain model in the presence of microscopic fibrotic structures

Authors: Brodie A. J. Lawson, Rodrigo Weber dos Santos, Ian W. Turner, Alfonso Bueno-Orovio, Pamela Burrage, Kevin Burrage

Abstract: Computational models in cardiac electrophysiology are notorious for long runtimes, restricting the numbers of nodes and mesh elements in the numerical discretisations used for their solution. This makes it particularly challenging to incorporate structural heterogeneities on small spatial scales, preventing a full understanding of the critical arrhythmogenic effects of conditions such as cardiac f… ▽ More Computational models in cardiac electrophysiology are notorious for long runtimes, restricting the numbers of nodes and mesh elements in the numerical discretisations used for their solution. This makes it particularly challenging to incorporate structural heterogeneities on small spatial scales, preventing a full understanding of the critical arrhythmogenic effects of conditions such as cardiac fibrosis. In this work, we explore the technique of homogenisation by volume averaging for the inclusion of non-conductive micro-structures into larger-scale cardiac meshes with minor computational overhead. Importantly, our approach is not restricted to periodic patterns, enabling homogenised models to represent, for example, the intricate patterns of collagen deposition present in different types of fibrosis. We first highlight the importance of appropriate boundary condition choice for the closure problems that define the parameters of homogenised models. Then, we demonstrate the technique's ability to correctly upscale the effects of fibrotic patterns with a spatial resolution of 10 $μ$m into much larger numerical mesh sizes of 100-250 $μ$m. The homogenised models using these coarser meshes correctly predict critical pro-arrhythmic effects of fibrosis, including slowed conduction, source/sink mismatch, and stabilisation of re-entrant activation patterns. As such, this approach to homogenisation represents a significant step towards whole organ simulations that unravel the effects of microscopic cardiac tissue heterogeneities. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: 17 pages, 7 figures in main

MSC Class: 35B27

Implicit reconstructions of thin leaf surfaces from large, noisy point clouds

Authors: Riley M. Whebell, Timothy J. Moroney, Ian W. Turner, Ravindra Pethiyagoda, Scott W. McCue

Abstract: Thin surfaces, such as the leaves of a plant, pose a significant challenge for implicit surface reconstruction techniques, which typically assume a closed, orientable surface. We show that by approximately interpolating a point cloud of the surface (augmented with off-surface points) and restricting the evaluation of the interpolant to a tight domain around the point cloud, we need only require an… ▽ More Thin surfaces, such as the leaves of a plant, pose a significant challenge for implicit surface reconstruction techniques, which typically assume a closed, orientable surface. We show that by approximately interpolating a point cloud of the surface (augmented with off-surface points) and restricting the evaluation of the interpolant to a tight domain around the point cloud, we need only require an orientable surface for the reconstruction. We use polyharmonic smoothing splines to fit approximate interpolants to noisy data, and a partition of unity method with an octree-like strategy for choosing subdomains. This method enables us to interpolate an N-point dataset in O(N) operations. We present results for point clouds of capsicum and tomato plants, scanned with a handheld device. An important outcome of the work is that sufficiently smooth leaf surfaces are generated that are amenable for droplet spreading simulations. △ Less

Submitted 12 April, 2021; v1 submitted 21 September, 2020; originally announced September 2020.

Numerical investigation into coarse-scale models of diffusion in complex heterogeneous media

Authors: Nathan G. March, Elliot J. Carr, Ian W. Turner

Abstract: Computational modelling of diffusion in heterogeneous media is prohibitively expensive for problems with fine-scale heterogeneities. A common strategy for resolving this issue is to decompose the domain into a number of non-overlapping sub-domains and homogenize the spatially-dependent diffusivity within each sub-domain (homogenization cell). This process yields a coarse-scale model for approximat… ▽ More Computational modelling of diffusion in heterogeneous media is prohibitively expensive for problems with fine-scale heterogeneities. A common strategy for resolving this issue is to decompose the domain into a number of non-overlapping sub-domains and homogenize the spatially-dependent diffusivity within each sub-domain (homogenization cell). This process yields a coarse-scale model for approximating the solution behaviour of the original fine-scale model at a reduced computational cost. In this paper, we study coarse-scale diffusion models in block heterogeneous media and investigate, for the first time, the effect that various factors have on the accuracy of resulting coarse-scale solutions. We present new findings on the error associated with homogenization as well as confirm via numerical experimentation that periodic boundary conditions are the best choice for the homogenization cell and demonstrate that the smallest homogenization cell that is computationally feasible should be used in numerical simulations. △ Less

Submitted 24 August, 2021; v1 submitted 12 September, 2020; originally announced September 2020.

Comments: 22 pages, 10 figures, accepted version

Journal ref: Transport in Porous Media (2021)

The ELFIN Mission

Authors: V. Angelopoulos, E. Tsai, L. Bingley, C. Shaffer, D. L. Turner, A. Runov, W. Li, J. Liu, A. V. Artemyev, X. -J. Zhang, R. J. Strangeway, R. E. Wirz, Y. Y. Shprits, V. A. Sergeev, R. P. Caron, M. Chung, P. Cruce, W. Greer, E. Grimes, K. Hector, M. J. Lawson, D. Leneman, E. V. Masongsong, C. L. Russell, C. Wilkins , et al. (57 additional authors not shown)

Abstract: The Electron Loss and Fields Investigation with a Spatio-Temporal Ambiguity-Resolving option (ELFIN-STAR, or simply: ELFIN) mission comprises two identical 3-Unit (3U) CubeSats on a polar (~93deg inclination), nearly circular, low-Earth (~450 km altitude) orbit. Launched on September 15, 2018, ELFIN is expected to have a >2.5 year lifetime. Its primary science objective is to resolve the mechanism… ▽ More The Electron Loss and Fields Investigation with a Spatio-Temporal Ambiguity-Resolving option (ELFIN-STAR, or simply: ELFIN) mission comprises two identical 3-Unit (3U) CubeSats on a polar (~93deg inclination), nearly circular, low-Earth (~450 km altitude) orbit. Launched on September 15, 2018, ELFIN is expected to have a >2.5 year lifetime. Its primary science objective is to resolve the mechanism of storm-time relativistic electron precipitation, for which electromagnetic ion cyclotron (EMIC) waves are a prime candidate. From its ionospheric vantage point, ELFIN uses its unique pitch-angle-resolving capability to determine whether measured relativistic electron pitch-angle and energy spectra within the loss cone bear the characteristic signatures of scattering by EMIC waves or whether such scattering may be due to other processes. Pairing identical ELFIN satellites with slowly-variable along-track separation allows disambiguation of spatial and temporal evolution of the precipitation over minutes-to-tens-of-minutes timescales, faster than the orbit period of a single low-altitude satellite (~90min). Each satellite carries an energetic particle detector for electrons (EPDE) that measures 50keV to 5MeV electrons with deltaE/E<40% and a fluxgate magnetometer (FGM) on a ~72cm boom that measures magnetic field waves (e.g., EMIC waves) in the range from DC to 5Hz Nyquist (nominally) with <0.3nT/sqrt(Hz) noise at 1Hz. The spinning satellites (T_spin~3s) are equipped with magnetorquers that permit spin-up/down and reorientation maneuvers. The spin axis is placed normal to the orbit plane, allowing full pitch-angle resolution twice per spin. An energetic particle detector for ions (EPDI) measures 250keV-5MeV ions, addressing secondary science. Funded initially by CalSpace and the University Nanosat Program, ELFIN was selected for flight with joint support from NSF and NASA between 2014 and 2018. △ Less

Submitted 16 June, 2020; v1 submitted 13 June, 2020; originally announced June 2020.

Comments: Submitted to Space Science Reviews April 2020. 51 pages, 7 tables, 21 figures

The LUX-ZEPLIN (LZ) radioactivity and cleanliness control programs

Authors: D. S. Akerib, C. W. Akerlof, D. Yu. Akimov, A. Alquahtani, S. K. Alsum, T. J. Anderson, N. Angelides, H. M. Araújo, A. Arbuckle, J. E. Armstrong, M. Arthurs, H. Auyeung, S. Aviles, X. Bai, A. J. Bailey, J. Balajthy, S. Balashov, J. Bang, M. J. Barry, D. Bauer, P. Bauer, A. Baxter, J. Belle, P. Beltrame, J. Bensinger , et al. (365 additional authors not shown)

Abstract: LUX-ZEPLIN (LZ) is a second-generation direct dark matter experiment with spin-independent WIMP-nucleon scattering sensitivity above $1.4 \times 10^{-48}$ cm$^{2}$ for a WIMP mass of 40 GeV/c$^{2}$ and a 1000 d exposure. LZ achieves this sensitivity through a combination of a large 5.6 t fiducial volume, active inner and outer veto systems, and radio-pure construction using materials with inherent… ▽ More LUX-ZEPLIN (LZ) is a second-generation direct dark matter experiment with spin-independent WIMP-nucleon scattering sensitivity above $1.4 \times 10^{-48}$ cm$^{2}$ for a WIMP mass of 40 GeV/c$^{2}$ and a 1000 d exposure. LZ achieves this sensitivity through a combination of a large 5.6 t fiducial volume, active inner and outer veto systems, and radio-pure construction using materials with inherently low radioactivity content. The LZ collaboration performed an extensive radioassay campaign over a period of six years to inform material selection for construction and provide an input to the experimental background model against which any possible signal excess may be evaluated. The campaign and its results are described in this paper. We present assays of dust and radon daughters depositing on the surface of components as well as cleanliness controls necessary to maintain background expectations through detector construction and assembly. Finally, examples from the campaign to highlight fixed contaminant radioassays for the LZ photomultiplier tubes, quality control and quality assurance procedures through fabrication, radon emanation measurements of major sub-systems, and bespoke detector systems to assay scintillator are presented. △ Less

Submitted 28 February, 2022; v1 submitted 3 June, 2020; originally announced June 2020.

Comments: 45 pages (79 inc. tables), 7 figures, 9 tables

Journal ref: The European Physical Journal C, Volume 80, Article number: 1044 (2020)

Measurement of energy flow, cross section and average inelasticity of forward neutrons produced in $\mathrm{\sqrt{s} = 13 TeV}$ proton-proton collisions with the LHCf Arm2 detector

Authors: O. Adriani, E. Berti, L. Bonechi, M. Bongi, R. D'Alessandro, S. Detti, M. Haguenauer, Y. Itow, K. Kasahara, H. Menjo, Y. Muraki, K. Ohashi, P. Papini, S. Ricciarini, T. Sako, N. Sakurai, K. Sato, T. Tamura, A. Tiberio, S. Torii, A. Tricomi, W. C. Turner, M. Ueno

Abstract: In this paper, we report the measurement of the energy flow, the cross section and the average inelasticity of forward neutrons (+ antineutrons) produced in $\sqrt{s} = 13$ TeV proton-proton collisions. These quantities are obtained from the inclusive differential production cross section, measured using the LHCf Arm2 detector at the CERN Large Hadron Collider. The measurements are performed in si… ▽ More In this paper, we report the measurement of the energy flow, the cross section and the average inelasticity of forward neutrons (+ antineutrons) produced in $\sqrt{s} = 13$ TeV proton-proton collisions. These quantities are obtained from the inclusive differential production cross section, measured using the LHCf Arm2 detector at the CERN Large Hadron Collider. The measurements are performed in six pseudorapidity regions: three of them ($η> 10.75$, $8.99 < η< 9.21$ and $8.80 < η< 8.99$), albeit with smaller acceptance and larger uncertainties, were already published in a previous work, whereas the remaining three ($10.06 < η< 10.75$, $9.65 < η< 10.06$ and $8.65 < η< 8.80$) are presented here for the first time. The analysis was carried out using a data set acquired in June 2015 with a corresponding integrated luminosity of $\mathrm{0.194~nb^{-1}}$. Comparing the experimental measurements with the expectations of several hadronic interaction models used to simulate cosmic ray air showers, none of these generators resulted to have a satisfactory agreement in all the phase space selected for the analysis. The inclusive differential production cross section for $η> 10.75$ is not reproduced by any model, whereas the results still indicate a significant but less serious deviation at lower pseudorapidities. Depending on the pseudorapidity region, the generators showing the best overall agreement with data are either SIBYLL 2.3 or EPOS-LHC. Furthermore, apart from the most forward region, the derived energy flow and cross section distributions are best reproduced by EPOS-LHC. Finally, even if none of the models describe the elasticity distribution in a satisfactory way, the extracted average inelasticity is consistent with the QGSJET II-04 value, while most of the other generators give values that lie just outside the experimental uncertainties. △ Less

Submitted 10 July, 2020; v1 submitted 4 March, 2020; originally announced March 2020.

Comments: 22 pages, 11 figures

Report number: CERN-EP-2020-029

Simulations of Events for the LUX-ZEPLIN (LZ) Dark Matter Experiment

Authors: The LUX-ZEPLIN Collaboration, :, D. S. Akerib, C. W. Akerlof, A. Alqahtani, S. K. Alsum, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, X. Bai, J. Balajthy, S. Balashov, J. Bang, D. Bauer, A. Baxter, J. Bensinger, E. P. Bernard, A. Bernstein, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, K. E. Boast , et al. (173 additional authors not shown)

Abstract: The LUX-ZEPLIN dark matter search aims to achieve a sensitivity to the WIMP-nucleon spin-independent cross-section down to (1--2)$\times10^{-12}$\,pb at a WIMP mass of 40 GeV/$c^2$. This paper describes the simulations framework that, along with radioactivity measurements, was used to support this projection, and also to provide mock data for validating reconstruction and analysis software. Of par… ▽ More The LUX-ZEPLIN dark matter search aims to achieve a sensitivity to the WIMP-nucleon spin-independent cross-section down to (1--2)$\times10^{-12}$\,pb at a WIMP mass of 40 GeV/$c^2$. This paper describes the simulations framework that, along with radioactivity measurements, was used to support this projection, and also to provide mock data for validating reconstruction and analysis software. Of particular note are the event generators, which allow us to model the background radiation, and the detector response physics used in the production of raw signals, which can be converted into digitized waveforms similar to data from the operational detector. Inclusion of the detector response allows us to process simulated data using the same analysis routines as developed to process the experimental data. △ Less

Submitted 23 June, 2020; v1 submitted 25 January, 2020; originally announced January 2020.

Comments: 24 pages, 19 figures; Corresponding Authors: A. Cottle, V. Kudryavtsev, D. Woodward

Projected sensitivity of the LUX-ZEPLIN experiment to the $0νββ$ decay of $^{136}$Xe

Authors: D. S. Akerib, C. W. Akerlof, A. Alqahtani, S. K. Alsum, T. J. Anderson, N. Angelides, H. M. Araújo, J. E. Armstrong, M. Arthurs, X. Bai, J. Balajthy, S. Balashov, J. Bang, A. Baxter, J. Bensinger, E. P. Bernard, A. Bernstein, A. Bhatti, A. Biekert, T. P. Biesiadzinski, H. J. Birch, K. E. Boast, B. Boxer, P. Brás, J. H. Buckley , et al. (167 additional authors not shown)

Abstract: The LUX-ZEPLIN (LZ) experiment will enable a neutrinoless double beta decay search in parallel to the main science goal of discovering dark matter particle interactions. We report the expected LZ sensitivity to $^{136}$Xe neutrinoless double beta decay, taking advantage of the significant ($>$600 kg) $^{136}$Xe mass contained within the active volume of LZ without isotopic enrichment. After 1000 l… ▽ More The LUX-ZEPLIN (LZ) experiment will enable a neutrinoless double beta decay search in parallel to the main science goal of discovering dark matter particle interactions. We report the expected LZ sensitivity to $^{136}$Xe neutrinoless double beta decay, taking advantage of the significant ($>$600 kg) $^{136}$Xe mass contained within the active volume of LZ without isotopic enrichment. After 1000 live-days, the median exclusion sensitivity to the half-life of $^{136}$Xe is projected to be 1.06$\times$10$^{26}$ years (90% confidence level), similar to existing constraints. We also report the expected sensitivity of a possible subsequent dedicated exposure using 90% enrichment with $^{136}$Xe at 1.06$\times$10$^{27}$ years. △ Less

Submitted 24 April, 2020; v1 submitted 9 December, 2019; originally announced December 2019.

Comments: 13 pages, 7 figures, 2 tables, version 2 changes: additional clarifications requested by referee on Sections II.A, III.C, III.E, III.F and IV.B

Journal ref: Phys. Rev. C 102, 014602 (2020)

The LUX-ZEPLIN (LZ) Experiment

Authors: The LZ Collaboration, D. S. Akerib, C. W. Akerlof, D. Yu. Akimov, A. Alquahtani, S. K. Alsum, T. J. Anderson, N. Angelides, H. M. Araújo, A. Arbuckle, J. E. Armstrong, M. Arthurs, H. Auyeung, X. Bai, A. J. Bailey, J. Balajthy, S. Balashov, J. Bang, M. J. Barry, J. Barthel, D. Bauer, P. Bauer, A. Baxter, J. Belle, P. Beltrame , et al. (357 additional authors not shown)

Abstract: We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient n… ▽ More We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements. △ Less

Submitted 3 November, 2019; v1 submitted 20 October, 2019; originally announced October 2019.

Non-universal Casimir forces at Bose-Einstein condensation of an ideal gas: effect of Dirichlet boundary conditions

Authors: M. Napiorkowski, J. Piasecki, J. W. Turner

Abstract: We analyze the Casimir forces for an ideal Bose gas enclosed between two infinite parallel walls separated by the distance D. The walls are characterized by the Dirichlet boundary conditions. We show that if the thermodynamic state with Bose-Einstein condensate present is correctly approached along the path pertinent to the Dirichlet b.c. then the leading term describing the large-distance decay o… ▽ More We analyze the Casimir forces for an ideal Bose gas enclosed between two infinite parallel walls separated by the distance D. The walls are characterized by the Dirichlet boundary conditions. We show that if the thermodynamic state with Bose-Einstein condensate present is correctly approached along the path pertinent to the Dirichlet b.c. then the leading term describing the large-distance decay of thermal Casimir force between the walls is 1/D^2 with a non-universal amplitude. The next order correction is lnD/D^3. These observations remain in contrast with the decay law for both the periodic and Neumann boundary conditions for which the leading term is 1/D^3 with a universal amplitude. We associate this discrepancy with the non-zero D-dependent positive value of the one-particle ground state energy in the case of Dirichlet boundary conditions. △ Less

Submitted 10 October, 2019; originally announced October 2019.

A fast algorithm for semi-analytically solving the homogenization boundary value problem for block locally-isotropic heterogeneous media

Authors: Nathan G. March, Elliot J. Carr, Ian W. Turner

Abstract: Direct numerical simulation of diffusion through heterogeneous media can be difficult due to the computational cost of resolving fine-scale heterogeneities. One method to overcome this difficulty is to homogenize the model by replacing the spatially-varying fine-scale diffusivity with an effective diffusivity calculated from the solution of an appropriate boundary value problem. In this paper, we… ▽ More Direct numerical simulation of diffusion through heterogeneous media can be difficult due to the computational cost of resolving fine-scale heterogeneities. One method to overcome this difficulty is to homogenize the model by replacing the spatially-varying fine-scale diffusivity with an effective diffusivity calculated from the solution of an appropriate boundary value problem. In this paper, we present a new semi-analytical method for solving this boundary value problem and computing the effective diffusivity for pixellated, locally-isotropic, heterogeneous media. We compare our new solution method to a standard finite volume method and show that equivalent accuracy can be achieved in less computational time for several standard test cases. We also demonstrate how the new solution method can be applied to complex heterogeneous geometries represented by a two-dimensional grid of rectangular blocks. These results indicate that our new semi-analytical method has the potential to significantly speed up simulations of diffusion in heterogeneous media. △ Less

Submitted 22 November, 2020; v1 submitted 17 December, 2018; originally announced December 2018.

Comments: 28 pages, 4 figures, 4 tables, accepted version of paper published in Applied Mathematical Modelling

Journal ref: Appl. Math. Model. 92 (2021) 23-43

Computational modelling of cardiac ischaemia using a variable-order fractional Laplacian

Authors: Megan E. Farquhar, Timothy J. Moroney, Qianqian Yang, Ian W. Turner, Kevin Burrage

Abstract: Heart failure is one of the most common causes of death in the western world. Many heart problems are linked to disturbances in cardiac electrical activity, such as wave re-entry caused by ischaemia. In terms of mathematical modelling, the monodomain equation is widely used to model electrical activity in the heart. Recently, Bueno-Orovio et al. [J. R. Soc. Interface 11: 20140352, 2014] pioneered… ▽ More Heart failure is one of the most common causes of death in the western world. Many heart problems are linked to disturbances in cardiac electrical activity, such as wave re-entry caused by ischaemia. In terms of mathematical modelling, the monodomain equation is widely used to model electrical activity in the heart. Recently, Bueno-Orovio et al. [J. R. Soc. Interface 11: 20140352, 2014] pioneered the use of a fractional Laplacian operator in the monodomain equation to account for the complex heterogeneous structures in heart tissue. In this work we consider how to extend this approach to apply to hearts with regions of damaged tissue. This requires the use of a fractional Laplacian operator whose fractional order varies spatially. We develop efficient numerical methods capable of solving this challenging problem on domains ranging from simple one-dimensional intervals with uniform meshes, through to full three-dimensional geometries on unstructured meshes. Results are presented for several test problems in one dimension, demonstrating the effects of different fractional orders in regions of healthy and damaged tissue. Then we showcase some new results for a three-dimensional fractional monodomain equation with a Beeler-Reuter ionic current model on a rabbit heart mesh. These simulation results are found to exhibit wave re-entry behaviour, brought about only by varying the value of the fractional order in a region representing damaged tissue. △ Less

Submitted 21 September, 2018; originally announced September 2018.

Measurement of inclusive forward neutron production cross section in proton-proton collisions at $\mathrm{\sqrt{s} = 13~TeV}$ with the LHCf Arm2 detector

Authors: O. Adriani, E. Berti, L. Bonechi, M. Bongi, R. D'Alessandro, S. Detti, M. Haguenauer, Y. Itow, K. Kasahara, Y. Makino, K. Masuda, H. Menjo, Y. Muraki, K. Ohashi, P. Papini, S. Ricciarini, T. Sako, N. Sakurai, K. Sato, M. Shinoda, T. Suzuki, T. Tamura, A. Tiberio, S. Torii, A. Tricomi , et al. (3 additional authors not shown)

Abstract: In this paper, we report the measurement relative to the production of forward neutrons in proton-proton collisions at $\mathrm{\sqrt{s} = 13~TeV}$ obtained using the LHCf Arm2 detector at the Large Hadron Collider. The results for the inclusive differential production cross section are presented as a function of energy in three different pseudorapidity regions: $η> 10.76$, $8.99 < η< 9.22$ and… ▽ More In this paper, we report the measurement relative to the production of forward neutrons in proton-proton collisions at $\mathrm{\sqrt{s} = 13~TeV}$ obtained using the LHCf Arm2 detector at the Large Hadron Collider. The results for the inclusive differential production cross section are presented as a function of energy in three different pseudorapidity regions: $η> 10.76$, $8.99 < η< 9.22$ and $8.81 < η< 8.99$. The analysis was performed using a data set acquired in June 2015 that corresponds to an integrated luminosity of $\mathrm{0.194~nb^{-1}}$. The measurements were compared with the predictions of several hadronic interaction models used to simulate air showers generated by Ultra High Energy Cosmic Rays. None of these generators showed good agreement with the data for all pseudorapidity intervals. For $η> 10.76$, no model is able to reproduce the observed peak structure at around $\mathrm{5~TeV}$ and all models underestimate the total production cross section: among them, QGSJET II-04 shows the smallest deficit with respect to data for the whole energy range. For $8.99 < η< 9.22$ and $8.81 < η< 8.99$, the models having the best overall agreement with data are SIBYLL 2.3 and EPOS-LHC, respectively: in particular, in both regions SIBYLL 2.3 is able to reproduce the observed peak structure at around $\mathrm{1.5-2.5~TeV}$. △ Less

Submitted 29 November, 2018; v1 submitted 29 August, 2018; originally announced August 2018.

Report number: CERN-EP-2018-239

Journal ref: The LHCf collaboration, Adriani, O., Berti, E. et al. J. High Energ. Phys. (2018) 2018: 73

On the analysis of mixed-index time fractional differential equation systems

Authors: Kevin Burrage, Pamela M. Burrage, Ian W. Turner, Fanhai Zeng

Abstract: In this paper we study the class of mixed-index time fractional differential equations in which different components of the problem have different time fractional derivatives on the left hand side. We prove a theorem on the solution of the linear system of equations, which collapses to the well-known Mittag-Leffler solution in the case the indices are the same, and also generalises the solution of… ▽ More In this paper we study the class of mixed-index time fractional differential equations in which different components of the problem have different time fractional derivatives on the left hand side. We prove a theorem on the solution of the linear system of equations, which collapses to the well-known Mittag-Leffler solution in the case the indices are the same, and also generalises the solution of the so-called linear sequential class of time fractional problems. We also investigate the asymptotic stability properties of this class of problems using Laplace transforms and show how Laplace transforms can be used to write solutions as linear combinations of generalised Mittag-Leffler functions in some cases. Finally we illustrate our results with some numerical simulations. △ Less

Submitted 4 September, 2017; originally announced September 2017.

Comments: 21 pages, 6 figures (some are made up of sub-figures - there are 15 figures or sub-figures)

MSC Class: 34A08 (Primary) 65L05 (Secondary)

Measurement of forward photon production cross-section in proton-proton collisions at $\sqrt{s}$ = 13 TeV with the LHCf detector

Authors: O. Adriani, E. Berti, L. Bonechi, M. Bongi, R. D'Alessandro, M. Haguenauer, Y. Itow, T. Iwata, K. Kasahara, Y. Makino, K. Masuda, E. Matsubayashi, H. Menjo, Y. Muraki, P. Papini, S. Ricciarini, T. Sako, N. Sakurai, M. Shinoda, T. Suzuki, T. Tamura, A. Tiberio, S. Torii, A. Tricomi, W. C. Turner , et al. (2 additional authors not shown)

Abstract: In this paper, we report the production cross-section of forward photons in the pseudorapidity regions of $η\,>\,10.94$ and $8.99\,>\,η\,>\,8.81$, measured by the LHCf experiment with proton--proton collisions at $\sqrt{s}$ = 13 TeV. The results from the analysis of 0.191 $\mathrm{nb^{-1}}$ of data obtained in June 2015 are compared to the predictions of several hadronic interaction models that ar… ▽ More In this paper, we report the production cross-section of forward photons in the pseudorapidity regions of $η\,>\,10.94$ and $8.99\,>\,η\,>\,8.81$, measured by the LHCf experiment with proton--proton collisions at $\sqrt{s}$ = 13 TeV. The results from the analysis of 0.191 $\mathrm{nb^{-1}}$ of data obtained in June 2015 are compared to the predictions of several hadronic interaction models that are used in air-shower simulations for ultra-high-energy cosmic rays. Although none of the models agree perfectly with the data, EPOS-LHC shows the best agreement with the experimental data among the models. △ Less

Submitted 24 November, 2017; v1 submitted 22 March, 2017; originally announced March 2017.

Comments: 21 pages, 4 figures

Report number: CERN-EP-DRAFT-LHCf-2017-001

The BRAN luminosity detectors for the LHC

Authors: Howard S. Matis, Massimo Placidi, Alessandro Ratti, William C. Turner, Enrico Bravin, Ryoichi Miyamoto

Abstract: This paper describes the several phases which led, from the conceptual design, prototyping, construction and tests with beam, to the installation and operation of the BRAN (Beam RAte of Neutrals) relative luminosity monitors for the LHC. The detectors have been operating since 2009 to contribute, optimize and maintain the accelerator performance in the two high luminosity interaction regions (IR),… ▽ More This paper describes the several phases which led, from the conceptual design, prototyping, construction and tests with beam, to the installation and operation of the BRAN (Beam RAte of Neutrals) relative luminosity monitors for the LHC. The detectors have been operating since 2009 to contribute, optimize and maintain the accelerator performance in the two high luminosity interaction regions (IR), the IR1 (ATLAS) and the IR5 (CMS). The devices are gas ionization chambers installed inside a neutral particle absorber 140 m away from the Interaction Points in IR1 and IR5 and monitor the energy deposited by electromagnetic showers produced by high-energy neutral particles from the collisions. The detectors have the capability to resolve the bunch-by-bunch luminosity at the 40 MHz bunch rate, as well as to survive the extreme level of radiation during the nominal LHC operation. The devices have operated since the early commissioning phase of the accelerator over a broad range of luminosities reaching 1.4*10^34 cm-2 s-1 with a peak pileup of 45 events per bunch crossing. Even though the nominal design luminosity of the LHC has been exceeded, the BRAN is operating well. After describing the multiple applications that the BRAN can be used to monitor the luminosity of the accelerator, we discuss the technical choices that led to its construction and the different tests performed prior to the installation in two IRs of the LHC. Performance simulations are presented together with operational results obtained during p-p operations, including runs at 40 MHz bunch rate, Pb-Pb operations and p-Pb operations. △ Less

Submitted 4 December, 2016; originally announced December 2016.

Comments: 22 pages, 32 Figures

Scanning Quantum Cryogenic Atom Microscope

Authors: Fan Yang, Alicia J. Kollár, Stephen F. Taylor, Richard W. Turner, Benjamin L. Lev

Abstract: Microscopic imaging of local magnetic fields provides a window into the organizing principles of complex and technologically relevant condensed matter materials. However, a wide variety of intriguing strongly correlated and topologically nontrivial materials exhibit poorly understood phenomena outside the detection capability of state-of-the-art high-sensitivity, high-resolution scanning probe mag… ▽ More Microscopic imaging of local magnetic fields provides a window into the organizing principles of complex and technologically relevant condensed matter materials. However, a wide variety of intriguing strongly correlated and topologically nontrivial materials exhibit poorly understood phenomena outside the detection capability of state-of-the-art high-sensitivity, high-resolution scanning probe magnetometers. We introduce a quantum-noise-limited scanning probe magnetometer that can operate from room to cryogenic temperatures with unprecedented DC-field sensitivity and micron-scale resolution. The Scanning Quantum Cryogenic Atom Microscope (SQCRAMscope) employs a magnetically levitated atomic Bose-Einstein condensate (BEC), thereby providing immunity to conductive and blackbody radiative heating. It has a field sensitivity of 1.4 nT per resolution-limited point ($\sim$2 $μ$m), or 6 nT/$\sqrt{\text{Hz}}$ per point at its duty cycle. Compared to point-by-point sensors, the long length of the BEC provides a naturally parallel measurement, allowing one to measure nearly one-hundred points with an effective field sensitivity of 600 pT$/\sqrt{\text{Hz}}$ for each point during the same time as a point-by-point scanner would measure these points sequentially. Moreover, it has a noise floor of 300 pT and provides nearly two orders of magnitude improvement in magnetic flux sensitivity (down to $10^{-6}$ $Φ_0/\sqrt{\text{Hz}}$) over previous atomic probe magnetometers capable of scanning near samples. These capabilities are, for the first time, carefully benchmarked by imaging magnetic fields arising from microfabricated wire patterns, in a system where samples may be scanned, cryogenically cooled, and easily exchanged. The SQCRAMscope will provide charge transport images at temperatures from room to 4 K in unconventional superconductors and topologically nontrivial materials. △ Less

Submitted 28 March, 2017; v1 submitted 24 August, 2016; originally announced August 2016.

Comments: 17 pages, 18 figures

Journal ref: Phys. Rev. Applied 7, 034026 (2017)

Measurements of longitudinal and transverse momentum distributions for neutral pions in the forward-rapidity region with the LHCf detector

Authors: O. Adriani, E. Berti, L. Bonechi, M. Bongi, R. D'Alessandro, M. Del Prete, M. Haguenauer, Y. Itow, K. Kasahara, K. Kawade, Y. Makino, K. Masuda, E. Matsubayashi, H. Menjo, G. Mitsuka, Y. Muraki, P. Papini, A. -L. Perrot, S. Ricciarini, T. Sako, N. Sakurai, T. Suzuki, T. Tamura, A. Tiberio, S. Torii , et al. (3 additional authors not shown)

Abstract: The differential cross sections for inclusive neutral pions as a function of transverse and longitudinal momentum in the very forward rapidity region have been measured at the Large Hadron Collider (LHC) with the Large Hadron Collider forward detector (LHCf) in proton-proton collisions at $\sqrt{s}=$ 2.76 and 7 TeV and in proton-lead collisions at nucleon-nucleon center-of-mass energies of… ▽ More The differential cross sections for inclusive neutral pions as a function of transverse and longitudinal momentum in the very forward rapidity region have been measured at the Large Hadron Collider (LHC) with the Large Hadron Collider forward detector (LHCf) in proton-proton collisions at $\sqrt{s}=$ 2.76 and 7 TeV and in proton-lead collisions at nucleon-nucleon center-of-mass energies of $\sqrt{s_\text{NN}}=$ 5.02 TeV. Such differential cross sections in proton-proton collisions are compatible with the hypotheses of limiting fragmentation and Feynman scaling. Comparing proton-proton with proton-lead collisions, we find a sizable suppression of the production of neutral pions in the differential cross sections after subtraction of ultra-peripheral proton-lead collisions. This suppression corresponds to the nuclear modification factor value of about 0.1-0.3. The experimental measurements presented in this paper provide a benchmark for the hadronic interaction Monte Carlo simulation codes that are used for the simulation of cosmic ray air showers. △ Less

Submitted 2 August, 2016; v1 submitted 31 July, 2015; originally announced July 2015.

Comments: 47 pages, 23 figures, 59 tables, accepted in PRD

Report number: CERN-PH-EP-2015-201

Journal ref: Phys. Rev. D 94, 032007 (2016)