Neutrino Physics at Fermilab

The Liquid Scintillator Neutrino Detector (LSND), located at the LANSCE (formerly LAMPF) linear accelerator at Los Alamos National Laboratory, has seen evidence for the oscillation of neutrinos, and hence neutrino mass. That discovery was the impetus for this LDRD project, begun in 1996. The goal of this project was to define the appropriate technologies to use in a follow up experiment and to set in place the requirements for such an experiment

First Observation of Coherent π0\pi^0 Production in Neutrino Nucleus Interactions with Eν<E_{\nu}< 2 GeV

The MiniBooNE experiment at Fermilab has amassed the largest sample to date of $\pi^0$s produced in neutral current (NC) neutrino-nucleus interactions at low energy. This paper reports a measurement of the momentum distribution of $\pi^0$s produced in mineral oil (CH$_2$) and the first observation of coherent $\pi^0$ production below 2 GeV. In the forward direction, the yield of events observed above the expectation for resonant production is attributed primarily to coherent production off carbon, but may also include a small contribution from diffractive production on hydrogen. Integrated over the MiniBooNE neutrino flux, the sum of the NC coherent and diffractive modes is found to be (19.5 $\pm$1.1 (stat) $\pm$2.5 (sys))% of all exclusive NC $\pi^0$ production at MiniBooNE. These measurements are of immediate utility because they quantify an important background to MiniBooNE's search for $\nu_{\mu} \to \nu_e$ oscillations.Comment: Submitted to Phys. Lett.

Test of Lorentz and CPT violation with Short Baseline Neutrino Oscillation Excesses

The sidereal time dependence of MiniBooNE electron neutrino and anti-electron neutrino appearance data are analyzed to search for evidence of Lorentz and CPT violation. An unbinned Kolmogorov-Smirnov test shows both the electron neutrino and anti-electron neutrino appearance data are compatible with the null sidereal variation hypothesis to more than 5%. Using an unbinned likelihood fit with a Lorentz-violating oscillation model derived from the Standard Model Extension (SME) to describe any excess events over background, we find that the electron neutrino appearance data prefer a sidereal time-independent solution, and the anti-electron neutrino appearance data slightly prefer a sidereal time-dependent solution. Limits of order 10E-20 GeV are placed on combinations of SME coefficients. These limits give the best limits on certain SME coefficients for muon neutrino to electron neutrino and anti-muon neutrino to anti-electron neutrino oscillations. The fit values and limits of combinations of SME coefficients are provided.Comment: 14 pages, 3 figures, and 2 tables, submitted to Physics Letters

All-sky search for long-duration gravitational wave transients with initial LIGO

We present the results of a search for long-duration gravitational wave transients in two sets of data collected by the LIGO Hanford and LIGO Livingston detectors between November 5, 2005 and September 30, 2007, and July 7, 2009 and October 20, 2010, with a total observational time of 283.0 days and 132.9 days, respectively. The search targets gravitational wave transients of duration 10-500 s in a frequency band of 40-1000 Hz, with minimal assumptions about the signal waveform, polarization, source direction, or time of occurrence. All candidate triggers were consistent with the expected background; as a result we set 90% confidence upper limits on the rate of long-duration gravitational wave transients for different types of gravitational wave signals. For signals from black hole accretion disk instabilities, we set upper limits on the source rate density between 3.4×10-5 and 9.4×10-4 Mpc-3 yr-1 at 90% confidence. These are the first results from an all-sky search for unmodeled long-duration transient gravitational waves. © 2016 American Physical Society

All-sky search for long-duration gravitational wave transients with initial LIGO

We present the results of a search for long-duration gravitational wave transients in two sets of data collected by the LIGO Hanford and LIGO Livingston detectors between November 5, 2005 and September 30, 2007, and July 7, 2009 and October 20, 2010, with a total observational time of 283.0 days and 132.9 days, respectively. The search targets gravitational wave transients of duration 10-500 s in a frequency band of 40-1000 Hz, with minimal assumptions about the signal waveform, polarization, source direction, or time of occurrence. All candidate triggers were consistent with the expected background; as a result we set 90% confidence upper limits on the rate of long-duration gravitational wave transients for different types of gravitational wave signals. For signals from black hole accretion disk instabilities, we set upper limits on the source rate density between 3.4×10-5 and 9.4×10-4 Mpc-3 yr-1 at 90% confidence. These are the first results from an all-sky search for unmodeled long-duration transient gravitational waves. © 2016 American Physical Society

Search for gravitational waves from Scorpius X-1 in the second Advanced LIGO observing run with an improved hidden Markov model

We present results from a semicoherent search for continuous gravitational waves from the low-mass x-ray binary Scorpius X-1, using a hidden Markov model (HMM) to track spin wandering. This search improves on previous HMM-based searches of LIGO data by using an improved frequency domain matched filter, the J-statistic, and by analyzing data from Advanced LIGO's second observing run. In the frequency range searched, from 60 to 650 Hz, we find no evidence of gravitational radiation. At 194.6 Hz, the most sensitive search frequency, we report an upper limit on gravitational wave strain (at 95% confidence) of h095%=3.47×10-25 when marginalizing over source inclination angle. This is the most sensitive search for Scorpius X-1, to date, that is specifically designed to be robust in the presence of spin wandering. © 2019 American Physical Society

Search for Tensor, Vector, and Scalar Polarizations in the Stochastic Gravitational-Wave Background

The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational waves is encoded in the spectral shape of the stochastic gravitational-wave background, formed by the superposition of cosmological and individually unresolved astrophysical sources. Using data recorded by Advanced LIGO during its first observing run, we search for a stochastic background of generically polarized gravitational waves. We find no evidence for a background of any polarization, and place the first direct bounds on the contributions of vector and scalar polarizations to the stochastic background. Under log-uniform priors for the energy in each polarization, we limit the energy densities of tensor, vector, and scalar modes at 95% credibility to Ω0T<5.58×10-8, Ω0V<6.35×10-8, and Ω0S<1.08×10-7 at a reference frequency f0=25 Hz. © 2018 American Physical Society

On the progenitor of binary neutron star merger GW170817

On 2017 August 17 the merger of two compact objects with masses consistent with two neutron stars was discovered through gravitational-wave (GW170817), gamma-ray (GRB 170817A), and optical (SSS17a/AT 2017gfo) observations. The optical source was associated with the early-type galaxy NGC 4993 at a distance of just ∼40 Mpc, consistent with the gravitational-wave measurement, and the merger was localized to be at a projected distance of ∼2 kpc away from the galaxy's center. We use this minimal set of facts and the mass posteriors of the two neutron stars to derive the first constraints on the progenitor of GW170817 at the time of the second supernova (SN). We generate simulated progenitor populations and follow the three-dimensional kinematic evolution from binary neutron star (BNS) birth to the merger time, accounting for pre-SN galactic motion, for considerably different input distributions of the progenitor mass, pre-SN semimajor axis, and SN-kick velocity. Though not considerably tight, we find these constraints to be comparable to those for Galactic BNS progenitors. The derived constraints are very strongly influenced by the requirement of keeping the binary bound after the second SN and having the merger occur relatively close to the center of the galaxy. These constraints are insensitive to the galaxy's star formation history, provided the stellar populations are older than 1 Gyr

THE RATE OF BINARY BLACK HOLE MERGERS INFERRED FROM ADVANCED LIGO OBSERVATIONS SURROUNDING GW150914

A transient gravitational-wave signal, GW150914, was identi fi ed in the twin Advanced LIGO detectors on 2015 September 2015 at 09:50:45 UTC. To asse ss the implications of this discovery, the detectors remained in operation with unchanged con fi gurations over a period of 39 days around the time of t he signal. At the detection statistic threshold corresponding to that observed for GW150914, our search of the 16 days of simultaneous two-detector observational data is estimated to have a false-alarm rate ( FAR ) of < ́ -- 4.9 10 yr 61 , yielding a p -value for GW150914 of < ́ - 210 7 . Parameter estimation follo w-up on this trigger identi fi es its source as a binary black hole ( BBH ) merger with component masses ( )( ) = - + - + mm M ,36,29 12 4 5 4 4 at redshift = - + z 0.09 0.04 0.03 ( median and 90% credible range ) . Here, we report on the constraints these observations place on the rate of BBH coalescences. Considering only GW150914, assuming that all BBHs in the universe have the same masses and spins as this event, imposing a search FAR threshold of 1 per 100 years, and assuming that the BBH merger rate is constant in the comoving frame, we infer a 90% credible range of merger rates between – -- 2 53 Gpc yr 31 ( comoving frame ) . Incorporating all search triggers that pass a much lower threshold while accounting for the uncerta inty in the astrophysical origin of each trigger, we estimate a higher rate, ranging from – -- 13 600 Gpc yr 31 depending on assumptions about the BBH mass distribution. All together, our various rate estimat es fall in the conservative range – -- 2 600 Gpc yr 31

Highly-parallelized simulation of a pixelated LArTPC on a GPU

The rapid development of general-purpose computing on graphics processing units (GPGPU) is allowing the implementation of highly-parallelized Monte Carlo simulation chains for particle physics experiments. This technique is particularly suitable for the simulation of a pixelated charge readout for time projection chambers, given the large number of channels that this technology employs. Here we present the first implementation of a full microphysical simulator of a liquid argon time projection chamber (LArTPC) equipped with light readout and pixelated charge readout, developed for the DUNE Near Detector. The software is implemented with an end-to-end set of GPU-optimized algorithms. The algorithms have been written in Python and translated into CUDA kernels using Numba, a just-in-time compiler for a subset of Python and NumPy instructions. The GPU implementation achieves a speed up of four orders of magnitude compared with the equivalent CPU version. The simulation of the current induced on 10^3 pixels takes around 1 ms on the GPU, compared with approximately 10 s on the CPU. The results of the simulation are compared against data from a pixel-readout LArTPC prototype