Showing 1–50 of 280 results for author: Johnson, R

Point containment algorithms for constructive solid geometry with unbounded primitives

Authors: Paul K. Romano, Patrick A. Myers, Seth R. Johnson, Aljaž Kolšek, Patrick C. Shriwise

Abstract: We present several algorithms for evaluating point containment in constructive solid geometry (CSG) trees with unbounded primitives. Three algorithms are presented based on postfix, prefix, and infix notations of the CSG binary expression tree. We show that prefix and infix notations enable short-circuiting logic, which reduces the number of primitives that must be checked during point containment… ▽ More We present several algorithms for evaluating point containment in constructive solid geometry (CSG) trees with unbounded primitives. Three algorithms are presented based on postfix, prefix, and infix notations of the CSG binary expression tree. We show that prefix and infix notations enable short-circuiting logic, which reduces the number of primitives that must be checked during point containment. To evaluate the performance of the algorithms, each algorithm was implemented in the OpenMC Monte Carlo particle transport code, which relies on CSG to represent solid bodies through which subatomic particles travel. Two sets of tests were carried out. First, the execution time to generate a high-resolution rasterized image of a 2D slice of a detailed CSG model of the ITER tokamak was measured. Use of both prefix and infix notations offered significant speedup over the postfix notation that has traditionally been used in particle transport codes, with infix resulting in a 6$\times$ reduction in execution time relative to postfix. We then measured the execution time of neutron transport simulations of the same ITER model using each of the algorithms. The results and performance improvements reveal the same trends as for the rasterization test, with a 4.59$\times$ overall speedup using the infix notation relative to the original postfix notation in OpenMC. △ Less

Submitted 18 June, 2024; originally announced June 2024.

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

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

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

Submitted 14 June, 2024; originally announced June 2024.

Report number: FERMILAB-PUB-24-0287

Targeting spectroscopic accuracy for dispersion bound systems from ab initio techniques: translational eigenstates of Ne@C$_{70}$

Authors: K. Panchagnula, D. Graf, E. R. Johnson, A. J. W. Thom

Abstract: We investigate the endofullerene system Ne@C$_{70}$, by constructing a three-dimensional Potential Energy Surface (PES) describing the translational motion of the Ne atom. We compare a plethora of electronic structure methods including: MP2, SCS-MP2, SOS-MP2, RPA@PBE, C(HF)-RPA, which were previously used for He@C$_{60}$ in J. Chem. Phys. 160, 104303 (2024), alongside B86bPBE-25X-XDM and B86bPBE-5… ▽ More We investigate the endofullerene system Ne@C$_{70}$, by constructing a three-dimensional Potential Energy Surface (PES) describing the translational motion of the Ne atom. We compare a plethora of electronic structure methods including: MP2, SCS-MP2, SOS-MP2, RPA@PBE, C(HF)-RPA, which were previously used for He@C$_{60}$ in J. Chem. Phys. 160, 104303 (2024), alongside B86bPBE-25X-XDM and B86bPBE-50X-XDM. The reduction in symmetry moving from C$_{60}$ to C$_{70}$ introduces a double well potential along the anisotropic direction, which forms a test of the sensitivity and effectiveness of the methods. Due to the large cost of these calculations, the PES is interpolated using Gaussian Process Regression due to its effectiveness with sparse training data. The nuclear Hamiltonian is diagonalised using a symmetrised double minimum basis set outlined in J. Chem. Phys. 159, 164308 (2023), with translational energies having error bars $\pm 1$ and $\pm 2$ cm$^{-1}$. We quantify the shape of the ground state wavefunction by considering its prolateness and kurtosis, and compare the eigenfunctions between electronic structure methods from their Hellinger distance. We find no consistency between electronic structure methods as they find a range of barrier heights and minima positions of the double well, and different translational eigenspectra which also differ from the Lennard-Jones (LJ) PES given in J. Chem. Phys. 101, 2126,2140 (1994). We find that generating effective LJ parameters for each electronic structure method cannot reproduce the full PES, nor recreate the eigenstates and this suggests that the LJ form of the PES, while simple, may not be best suited to describe these systems. Even though MP2 and RPA@PBE performed best for He@C$_{60}$, due to the lack of concordance between all electronic structure methods we require more experimental data in order to properly validate the choice. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 12 pages, 4 figures; SI 6 pages, 3 figures

Spectral Methods for Coastal-Trapped Waves and Instabilities in a Background Flow

Authors: Matthew N. Crowe, Edward R. Johnson

Abstract: Here we present a numerical method for finding non-hydrostatic coastal-trapped wave and instability solutions to the non-hydrostatic Boussinesq equations in the presence of a background flow and complicated coastal topography. We use spectral methods to discretise the two-dimensional eigenvalue problem and solve the resulting discrete problem by standard methods. Our approach is applied to three e… ▽ More Here we present a numerical method for finding non-hydrostatic coastal-trapped wave and instability solutions to the non-hydrostatic Boussinesq equations in the presence of a background flow and complicated coastal topography. We use spectral methods to discretise the two-dimensional eigenvalue problem and solve the resulting discrete problem by standard methods. Our approach is applied to three examples and shown to be consistent with previous numerical and analytical results. In particular, we show that our method is able to reliably identify coastal-trapped wave solutions that correspond to waves seen in realistic simulations of the Southeast Greenland shelf. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 12 pages, 7 figures

On forced RF generation of CW magnetrons for SRF accelerators

Authors: G. Kazakevich, R. P. Johnson, T. Khabiboulline, G. Romanov, V. Yakovlev, Ya Derbenev, Yu. Eidelman

Abstract: CW magnetrons, initially developed for industrial RF heaters, were suggested to power RF cavities of superconducting accelerators due to their higher efficiency and lower cost than traditionally used klystrons, IOTs or solid-state amplifiers. RF amplifiers driven by a master oscillator serve as coherent RF sources. CW magnetrons are regenerative RF generators with a huge regenerative gain. This ca… ▽ More CW magnetrons, initially developed for industrial RF heaters, were suggested to power RF cavities of superconducting accelerators due to their higher efficiency and lower cost than traditionally used klystrons, IOTs or solid-state amplifiers. RF amplifiers driven by a master oscillator serve as coherent RF sources. CW magnetrons are regenerative RF generators with a huge regenerative gain. This causes regenerative instability with a large noise when a magnetron operates with the anode voltage above the threshold of self-excitation. Traditionally for stabilization of magnetrons is used injection locking by a quite small signal. Then the magnetron except the injection locked oscillations may generate noise. This may preclude use of standard CW magnetrons in some SRF accelerators. Recently we developed briefly described below a mode for forced RF generation of CW magnetrons when the magnetron startup is provided by the injected forcing signal and the regenerative noise is suppressed. The mode is most suitable for powering high Q-factor SRF cavities. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 7 pages, 12 figures, one table

Report number: FERMILAB-PUB-24-0129-TD

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

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

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

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

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

Modon solutions in an N-layer quasi-geostrophic model

Authors: Matthew N. Crowe, Edward R. Johnson

Abstract: Modons, or dipolar vortices, are common and long-lived features of the upper ocean, consisting of a pair of monopolar vortices moving through self-advection. Such structures remain stable over long times and may be important for fluid transport over large distances. Here we present a semi-analytical method for finding fully nonlinear modon solutions in a multi-layer quasi-geostrophic model with ar… ▽ More Modons, or dipolar vortices, are common and long-lived features of the upper ocean, consisting of a pair of monopolar vortices moving through self-advection. Such structures remain stable over long times and may be important for fluid transport over large distances. Here we present a semi-analytical method for finding fully nonlinear modon solutions in a multi-layer quasi-geostrophic model with arbitrarily many layers. Our approach works by reducing the problem to a multi-parameter linear eigenvalue problem which can be solved using numerical techniques from linear algebra. The method is shown to replicate previous results for one and two-layer models and is applied to a three-layer model to find a solution describing a mid-depth propagating, topographic vortex. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 11 pages, 3 figures

Attached and separated rotating flow over a finite height ridge

Authors: Stefan Frei, Erik Burman, Edward R Johnson

Abstract: This paper discusses the effect of rotation on the boundary layer in high Reynolds number flow over a ridge using a numerical method based on stabilised finite elements that captures steady solutions up to Reynolds number of order $10^6$. The results are validated against boundary layer computations in shallow flows and for deep flows against experimental observations reported in Machicoane et al.… ▽ More This paper discusses the effect of rotation on the boundary layer in high Reynolds number flow over a ridge using a numerical method based on stabilised finite elements that captures steady solutions up to Reynolds number of order $10^6$. The results are validated against boundary layer computations in shallow flows and for deep flows against experimental observations reported in Machicoane et al. (Phys. Rev. Fluids, 2018). In all cases considered the boundary layer remains attached, even at large Reynolds numbers, provided the Rossby number of the flow is sufficiently small. At any fixed Rossby number the flow detaches at sufficiently high Reynolds number to form a steady recirculating region in the lee of the ridge. At even higher Reynolds numbers no steady flow is found. This disappearance of steady solutions closely reproduces the transition to unsteadiness seen in the laboratory. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Prediction of storm surge on evolving landscapes under climate change

Authors: David R Johnson, Mohammad Ahmadi Gharehtoragh

Abstract: Planners who wish to manage coastal flood risk with long-lived infrastructure (e.g., levees, floodwalls) under a constrained computational budget face a tradeoff. Simulating a large number of future time periods or scenarios with different assumptions about land subsidence, sea level rise, land accretion, imposes a limit on how many storm simulations can be run in each scenario and time period. Ma… ▽ More Planners who wish to manage coastal flood risk with long-lived infrastructure (e.g., levees, floodwalls) under a constrained computational budget face a tradeoff. Simulating a large number of future time periods or scenarios with different assumptions about land subsidence, sea level rise, land accretion, imposes a limit on how many storm simulations can be run in each scenario and time period. Machine learning approaches have been developed to reduce the computational burden of predicting storm surge from many tropical cyclone events, but prior efforts focus on predicting surge as a function of storm parameters on a single landscape. In this analysis, we present a deep learning model that also incorporates landscape characteristics and boundary conditions (e.g., mean sea level). The model is informed by a dataset of peak surge elevations from Advanced Circulation (ADCIRC) hydrodynamic simulations of coastal Louisiana in eleven scenarios: a 2020 baseline and decadal time slices from 2030 to 2070 under two scenarios varying land subsidence and sea level rise rates. Training on ten scenarios to make predictions on the eleventh yields a grand RMSE of 0.086 m and grand MAE of 0.050 m over 90 storms per scenario and over 94,000 geospatial locations. We also aggregated the 90 storms in each scenario to generate an annual exceedance probability distribution; a two-sided Kolmogorov-Smirnov test comparing AEP estimates from the model predictions to the original ADCIRC simulations rejected the null hypothesis that the predictions and ADCIRC AEP values were drawn from the same distribution only 1.1% of the time. △ Less

Submitted 23 December, 2023; originally announced January 2024.

A note on reducing spurious pressure oscillations in fully conservative discontinuous Galerkin simulations of multicomponent flows

Authors: Eric J. Ching, Ryan F. Johnson, Andrew D. Kercher

Abstract: A well-known issue associated with the use of fully conservative schemes in multicomponent-flow simulations is the generation of spurious pressure oscillations at contact interfaces. These oscillations can rapidly lead to solver divergence even in the presence of smooth interfaces that are not fully resolved. In this note, we compare various strategies for reducing such oscillations that do not (a… ▽ More A well-known issue associated with the use of fully conservative schemes in multicomponent-flow simulations is the generation of spurious pressure oscillations at contact interfaces. These oscillations can rapidly lead to solver divergence even in the presence of smooth interfaces that are not fully resolved. In this note, we compare various strategies for reducing such oscillations that do not (a) introduce conservation error, (b) rely on artificial viscosity or limiting, or (c) degrade order of accuracy in smooth regions of the flow. The considered test case is one-dimensional advection of a high-pressure nitrogen/n-dodecane thermal bubble using the thermally perfect gas model. Several results are presented that contradict those corresponding to the more conventional hydrogen/oxygen thermal-bubble case. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Navier-Stokes equations

Authors: Eric J. Ching, Ryan F. Johnson, Sarah Burrows, Jacklyn Higgs, Andrew D. Kercher

Abstract: This article concerns the development of a fully conservative, positivity-preserving, and entropy-bounded discontinuous Galerkin scheme for the multicomponent, chemically reacting, compressible Navier-Stokes equations with complex thermodynamics. In particular, we extend to viscous flows the fully conservative, positivity-preserving, and entropy-bounded discontinuous Galerkin method for the chemic… ▽ More This article concerns the development of a fully conservative, positivity-preserving, and entropy-bounded discontinuous Galerkin scheme for the multicomponent, chemically reacting, compressible Navier-Stokes equations with complex thermodynamics. In particular, we extend to viscous flows the fully conservative, positivity-preserving, and entropy-bounded discontinuous Galerkin method for the chemically reacting Euler equations that we previously introduced. An important component of the formulation is the positivity-preserving Lax-Friedrichs-type viscous flux function devised by Zhang [J. Comput. Phys., 328 (2017), pp. 301-343], which was adapted to multicomponent flows by Du and Yang [J. Comput. Phys., 469 (2022), pp. 111548] in a manner that treats the inviscid and viscous fluxes as a single flux. Here, we similarly extend the aforementioned flux function to multicomponent flows but separate the inviscid and viscous fluxes, resulting in a different dissipation coefficient. This separation of the fluxes allows for use of other inviscid flux functions, as well as enforcement of entropy boundedness on only the convective contribution to the evolved state, as motivated by physical and mathematical principles. We also detail how to account for boundary conditions and incorporate previously developed techniques to reduce spurious pressure oscillations into the positivity-preserving framework. Furthermore, potential issues associated with the Lax-Friedrichs-type viscous flux function in the case of zero species concentrations are discussed and addressed. The resulting formulation is compatible with curved, multidimensional elements and general quadrature rules with positive weights. A variety of multicomponent, viscous flows is computed, ranging from a one-dimensional shock tube problem to multidimensional detonation waves and shock/mixing-layer interaction. △ Less

Submitted 28 June, 2024; v1 submitted 26 October, 2023; originally announced October 2023.

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

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

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

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

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

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

Radial Variation of Microstructure in a Direct-Chill Cast AA7050 Billet on Homogenization

Authors: Pikee Priya, Kyle Fezi, D. R. Johnson, M. J. M. Krane

Abstract: A "through-process" solidification homogenization numerical model to study the microstructural evolution has been developed for AA7050 alloy. A continuum scale Direct Chill Casting (DCC) solidification model has been coupled with a meso-scale homogenization precipitation model to evaluate the radial variation of microstructure in a cylindrical DC cast billet after homogenization. The Local Solidif… ▽ More A "through-process" solidification homogenization numerical model to study the microstructural evolution has been developed for AA7050 alloy. A continuum scale Direct Chill Casting (DCC) solidification model has been coupled with a meso-scale homogenization precipitation model to evaluate the radial variation of microstructure in a cylindrical DC cast billet after homogenization. The Local Solidification Times predicted by the DCC numerical model is used to predict the Secondary Dendritic Arm Spacing (SDAS) across the radius from an empirical relationship. It also predicts the macrosegregation used to estimate the radial as-cast microstructures (phases and phase fractions) using Thermo-Calc. Macrosegregation affects the initial and homogenized microstructures across the radius, making the surface prone to recrystallization, due to lesser dispersoids and larger precipitated particles during cooling, caused by lesser extent of Zenner pinning and Particle Stimulated Nucleation (PSN) of recrystallized grains respectively. △ Less

Submitted 26 August, 2023; originally announced August 2023.

Comments: 20 Pages, 8 Figures, part of the PhD Dissertation work by Pikee Priya at Purdue University (2016). Cited for reference

MSC Class: 74N25 ACM Class: G.0

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

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

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

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

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

Measurement of ambient radon progeny decay rates and energy spectra in liquid argon using the MicroBooNE detector

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

Abstract: We report measurements of radon progeny in liquid argon within the MicroBooNE time projection chamber (LArTPC). The presence of specific radon daughters in MicroBooNE's 85 metric tons of active liquid argon bulk is probed with newly developed charge-based low-energy reconstruction tools and analysis techniques to detect correlated $^{214}$Bi-$^{214}$Po radioactive decays. Special datasets taken du… ▽ More We report measurements of radon progeny in liquid argon within the MicroBooNE time projection chamber (LArTPC). The presence of specific radon daughters in MicroBooNE's 85 metric tons of active liquid argon bulk is probed with newly developed charge-based low-energy reconstruction tools and analysis techniques to detect correlated $^{214}$Bi-$^{214}$Po radioactive decays. Special datasets taken during periods of active radon doping enable new demonstrations of the calorimetric capabilities of single-phase neutrino LArTPCs for $β$ and $α$ particles with electron-equivalent energies ranging from 0.1 to 3.0 MeV. By applying $^{214}$Bi-$^{214}$Po detection algorithms to data recorded over a 46-day period, no statistically significant presence of radioactive $^{214}$Bi is detected, and a limit on the activity is placed at $<0.35$ mBq/kg at the 95% confidence level. This bulk $^{214}$Bi radiopurity limit -- the first ever reported for a liquid argon detector incorporating liquid-phase purification -- is then further discussed in relation to the targeted upper limit of 1 mBq/kg on bulk $^{222}$Rn activity for the DUNE neutrino detector. △ Less

Submitted 22 March, 2024; v1 submitted 6 July, 2023; originally announced July 2023.

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

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

First measurement of $η$ production in neutrino interactions on argon with MicroBooNE

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

Abstract: We present a measurement of $η$ production from neutrino interactions on argon with the MicroBooNE detector. The modeling of resonant neutrino interactions on argon is a critical aspect of the neutrino oscillation physics program being carried out by the DUNE and Short Baseline Neutrino programs. $η$ production in neutrino interactions provides a powerful new probe of resonant interactions, comple… ▽ More We present a measurement of $η$ production from neutrino interactions on argon with the MicroBooNE detector. The modeling of resonant neutrino interactions on argon is a critical aspect of the neutrino oscillation physics program being carried out by the DUNE and Short Baseline Neutrino programs. $η$ production in neutrino interactions provides a powerful new probe of resonant interactions, complementary to pion channels, and is particularly suited to the study of higher-order resonances beyond the $Δ(1232)$. We measure a flux-integrated cross section for neutrino-induced $η$ production on argon of $3.22 \pm 0.84 \; \textrm{(stat.)} \pm 0.86 \; \textrm{(syst.)}$ $10^{-41}{\textrm{cm}}^{2}$/nucleon. By demonstrating the successful reconstruction of the two photons resulting from $η$ production, this analysis enables a novel calibration technique for electromagnetic showers in GeV accelerator neutrino experiments. △ Less

Submitted 4 May, 2024; v1 submitted 25 May, 2023; originally announced May 2023.

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

Journal ref: Phys. Rev. Lett. 132, 151801 Published 10 April 2024

First demonstration of $\mathcal{O}(1\,\text{ns})$ timing resolution in the MicroBooNE liquid argon time projection chamber

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

Abstract: MicroBooNE is a neutrino experiment located in the Booster Neutrino Beamline (BNB) at Fermilab, which collected data from 2015 to 2021. MicroBooNE's liquid argon time projection chamber (LArTPC) is accompanied by a photon detection system consisting of 32 photomultiplier tubes used to measure the argon scintillation light and determine the timing of neutrino interactions. Analysis techniques combi… ▽ More MicroBooNE is a neutrino experiment located in the Booster Neutrino Beamline (BNB) at Fermilab, which collected data from 2015 to 2021. MicroBooNE's liquid argon time projection chamber (LArTPC) is accompanied by a photon detection system consisting of 32 photomultiplier tubes used to measure the argon scintillation light and determine the timing of neutrino interactions. Analysis techniques combining light signals and reconstructed tracks are applied to achieve a neutrino interaction time resolution of $\mathcal{O}(1\,\text{ns})$. The result obtained allows MicroBooNE to access the ns neutrino pulse structure of the BNB for the first time. The timing resolution achieved will enable significant enhancement of cosmic background rejection for all neutrino analyses. Furthermore, the ns timing resolution opens new avenues to search for long-lived-particles such as heavy neutral leptons in MicroBooNE, as well as in future large LArTPC experiments, namely the SBN program and DUNE. △ Less

Submitted 29 August, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Deflected Beam Method for Absolute Current Density Determination

Authors: Richard H. Mattish, Timothy J. Burke, Patrick R. Johnson, Chad E. Sosolik, Joan P. Marler

Abstract: We present a broadly applicable in situ method for profiling ion beams using electrostatic deflectors and a Faraday cup. By deconvolving the detector geometry from the resulting current profiles, spatially resolved absolute current density profiles are obtained. We demonstrate this method's efficacy with low-density highly charged ion beams (specifically, Ne$^{8+}$). Details on experimental design… ▽ More We present a broadly applicable in situ method for profiling ion beams using electrostatic deflectors and a Faraday cup. By deconvolving the detector geometry from the resulting current profiles, spatially resolved absolute current density profiles are obtained. We demonstrate this method's efficacy with low-density highly charged ion beams (specifically, Ne$^{8+}$). Details on experimental design are provided as well as the link to the deconvolution routine on Github. △ Less

Submitted 1 March, 2023; originally announced March 2023.

Comments: The following article has been submitted to Review of Scientific Instruments. After it is published, it will be found at https://aip.scitation.org/journal/rsi

The evolution of surface quasi-geostrophic modons on sloping topography

Authors: Matthew N. Crowe, Edward R. Johnson

Abstract: This work discusses modons, or dipolar vortices, propagating along sloping topography. Two different regimes exist which are studied separately using the surface quasi-geostrophic equations. First, when the modon propagates in the opposite direction to topographic Rossby waves, steady solutions exist and a semi-analytical method is presented for calculating these solutions. Second, when the modon… ▽ More This work discusses modons, or dipolar vortices, propagating along sloping topography. Two different regimes exist which are studied separately using the surface quasi-geostrophic equations. First, when the modon propagates in the opposite direction to topographic Rossby waves, steady solutions exist and a semi-analytical method is presented for calculating these solutions. Second, when the modon propagates in the same direction of the Rossby waves, a wave wake is generated. This wake removes energy from the modon causing it to decay slowly. Asymptotic predictions are presented for this decay and found to agree closely with numerical simulations. Over long times, decaying vortices are found to break down due to an asymmetry resulting from the generation of waves inside the vortex. A monopolar vortex moving along a wall is shown to behave in a similar way to a dipole, though the presence of the wall is found to stabilise the vortex and prevent the long-time breakdown. The problem is mathematically equivalent to a dipolar vortex moving along a density front hence our results apply directly to this case. △ Less

Submitted 10 February, 2023; originally announced February 2023.

Comments: 19 pages, 7 figures, 2 supplementary movies

Electron Impact Ionization in the Icy Galilean Satellites' Atmospheres

Authors: Shane R. Carberry Mogan, Robert E. Johnson, Audrey Vorburger, Lorenz Roth

Abstract: Electron impact ionization is critical in producing the ionospheres on many planetary bodies and, as discussed here, is critical for interpreting spacecraft and telescopic observations of the tenuous atmospheres of the icy Galilean satellites of Jupiter (Europa, Ganymede, and Callisto), which form an interesting planetary system. Fortunately, laboratory measurements, extrapolated by theoretical mo… ▽ More Electron impact ionization is critical in producing the ionospheres on many planetary bodies and, as discussed here, is critical for interpreting spacecraft and telescopic observations of the tenuous atmospheres of the icy Galilean satellites of Jupiter (Europa, Ganymede, and Callisto), which form an interesting planetary system. Fortunately, laboratory measurements, extrapolated by theoretical models, were developed and published over a number of years by K. Becker and colleagues (see Deutsch et al. 2009) to provide accurate electron impact ionization cross sections for atoms and molecules, which are crucial to correctly interpret these measurements. Because of their relevance for the Jovian icy satellites we provide useful fits to the complex, semi-empirical Deutsch-Mark formula for energy-dependent electron impact ionization cross-sections of gas-phase water products (i.e., H2O, H2, O2, H, O). These are then used with measurements of the thermal plasma in the Jovian magnetosphere to produce ionization rates for comparison with solar photo-ionization rates at the icy Galilean satellites. △ Less

Submitted 26 January, 2023; originally announced January 2023.

Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part II: The multidimensional case

Authors: Eric J. Ching, Ryan F. Johnson, Andrew D. Kercher

Abstract: In this second part of our two-part paper, we extend to multiple spatial dimensions the one-dimensional, fully conservative, positivity-preserving, and entropy-bounded discontinuous Galerkin scheme developed in the first part for the chemically reacting Euler equations. Our primary objective is to enable robust and accurate solutions to complex reacting-flow problems using the high-order discontin… ▽ More In this second part of our two-part paper, we extend to multiple spatial dimensions the one-dimensional, fully conservative, positivity-preserving, and entropy-bounded discontinuous Galerkin scheme developed in the first part for the chemically reacting Euler equations. Our primary objective is to enable robust and accurate solutions to complex reacting-flow problems using the high-order discontinuous Galerkin method without requiring extremely high resolution. Variable thermodynamics and detailed chemistry are considered. Our multidimensional framework can be regarded as a further generalization of similar positivity-preserving and/or entropy-bounded discontinuous Galerkin schemes in the literature. In particular, the proposed formulation is compatible with curved elements of arbitrary shape, a variety of numerical flux functions, general quadrature rules with positive weights, and mixtures of thermally perfect gases. Preservation of pressure equilibrium between adjacent elements, especially crucial in simulations of multicomponent flows, is discussed. Complex detonation waves in two and three dimensions are accurately computed using high-order polynomials. Enforcement of an entropy bound, as opposed to solely the positivity property, is found to significantly improve stability. Mass, total energy, and atomic elements are shown to be discretely conserved. △ Less

Submitted 15 January, 2024; v1 submitted 29 November, 2022; originally announced November 2022.

Journal ref: J. Comput. Phys. 505 (2024) 112878

Positivity-preserving and entropy-bounded discontinuous Galerkin method for the chemically reacting, compressible Euler equations. Part I: The one-dimensional case

Authors: Eric J. Ching, Ryan F. Johnson, Andrew D. Kercher

Abstract: In this paper, we develop a fully conservative, positivity-preserving, and entropy-bounded discontinuous Galerkin scheme for simulating the chemically reacting, compressible Euler equations with complex thermodynamics. The proposed formulation is an extension of the conservative, high-order numerical method previously developed by Johnson and Kercher [J. Comput. Phys., 423 (2020), 109826] that mai… ▽ More In this paper, we develop a fully conservative, positivity-preserving, and entropy-bounded discontinuous Galerkin scheme for simulating the chemically reacting, compressible Euler equations with complex thermodynamics. The proposed formulation is an extension of the conservative, high-order numerical method previously developed by Johnson and Kercher [J. Comput. Phys., 423 (2020), 109826] that maintains pressure equilibrium between adjacent elements. In this first part of our two-part paper, we focus on the one-dimensional case. Our methodology is rooted in the minimum entropy principle satisfied by entropy solutions to the multicomponent, compressible Euler equations, which was proved by Gouasmi et al. [ESAIM: Math. Model. Numer. Anal., 54 (2020), 373--389] for nonreacting flows. We first show that the minimum entropy principle holds in the reacting case as well. Next, we introduce the ingredients required for the solution to have nonnegative species concentrations, positive density, positive pressure, and bounded entropy. We also discuss how to retain the aforementioned ability to preserve pressure equilibrium between elements. Operator splitting is employed to handle stiff chemical reactions. To guarantee satisfaction of the minimum entropy principle in the reaction step, we develop an entropy-stable discontinuous Galerkin method based on diagonal-norm summation-by-parts operators for solving ordinary differential equations. The developed formulation is used to compute canonical one-dimensional test cases, namely thermal-bubble advection, multicomponent shock-tube flow, and a moving hydrogen-oxygen detonation wave with detailed chemistry. We find that the enforcement of an entropy bound can considerably reduce the large-scale nonlinear instabilities that emerge when only the positivity property is enforced, to an even greater extent than in the monocomponent, calorically perfect case. △ Less

Submitted 5 December, 2022; v1 submitted 29 November, 2022; originally announced November 2022.

Journal ref: J. Comput. Phys. 505 (2024) 112881

Oceanic dipoles in a surface quasigeostrophic model

Authors: E. R. Johnson, M. N. Crowe

Abstract: Analysis of satellite altimetry and Argo float data leads Ni et al. (2020, JGR Oceans) to argue that mesoscale dipoles are widespread features of the global ocean having a relatively uniform three-dimensional structure that can lead to strong vertical exchanges. Almost all the features of the composite dipole they construct can be derived from a model for multipoles in the surface quasigeostrophic… ▽ More Analysis of satellite altimetry and Argo float data leads Ni et al. (2020, JGR Oceans) to argue that mesoscale dipoles are widespread features of the global ocean having a relatively uniform three-dimensional structure that can lead to strong vertical exchanges. Almost all the features of the composite dipole they construct can be derived from a model for multipoles in the surface quasigeostrophic equations for which we present a straightforward novel solution in terms of an explicit linear eigenvalue problem, allowing simple evaluation of the higher radial modes that appear to be present in the observations and suggesting that mass conservation may explain the observed frontogenetic velocities. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Extended analysis of the effects of the Sumatra topography on downstream low-level vortex development over the Indian Ocean

Authors: Paul E. Ciesielski, Richard H. Johnson

Abstract: Fine et al. (2016, hereafter F16) investigated the potential role of Sumatra Island, as well as the Malay Peninsula and Java, in creating terrain-induced circulations over the Indian Ocean (IO) that subsequently develop into tropical cyclones (TCs). Applying sophisticated vortex tracking software to 2.5 yrs of model analyses, F16 found four regions downstream of topographic features in the Maritim… ▽ More Fine et al. (2016, hereafter F16) investigated the potential role of Sumatra Island, as well as the Malay Peninsula and Java, in creating terrain-induced circulations over the Indian Ocean (IO) that subsequently develop into tropical cyclones (TCs). Applying sophisticated vortex tracking software to 2.5 yrs of model analyses, F16 found four regions downstream of topographic features in the Maritime Continent to be prolific generators of low-level cyclonic vortices (123/yr). This present study extends the limited analyses of F16 by applying a similar approach to 10 yrs (2008-2017) of ERA5 analyses. While the 2.5-yr period which F16 studied was slightly (8%) more active in terms of vortex production than the 10-yr period, in general the findings of F16 are representative of the longer-term record with 80% of all shed vortices occurring with easterly low-level (925 hPa) flow over the hotspot regions. Additional analysis of the 10-yr record found that vortex counts are highest near MJO phase 1 when low-level easterlies are strongest over the maritime continent region. A secondary peak in vortex counts occurs during MJO phase 4 when low-level westerlies were present near the equator west of Sumatra. This suggests that low-level westerly surges on the equator impinging on Sumatra associated with the MJO contribute to an increase in wake vortex development. While the frequency of vortex genesis over the four hotspot regions is strongly tied to the annual cycle of winds, periods with anomalous zonal flow are shown to impact vortex counts. This is most apparent in the region off the southern tip of Sumatra where vortex counts were 4 times higher during periods of anomalous easterlies compared to periods with anomalous westerlies. In this region positive (negative) ENSO and Indian-Ocean Dipole conditions drive anomalous easterlies (westerlies) which impact vortex formation rates. △ Less

Submitted 25 October, 2022; originally announced October 2022.

Three-sided pyramid wavefront sensor. II. Preliminary demonstration on the new CACTI testbed

Authors: Lauren Schatz, Johanan Codona, Joseph D. Long, Jared R. Males, Weslin Pullen, Jennifer Lumbres, Kyle Van Gorkom, Vincent Chambouleyron, Laird M. Close, Carlos Correia, Olivier Fauvarque, Thierry Fusco, Olivier Guyon, Michael Hart, Pierre Janin-Potiron, Robert Johnson, Nemanja Jovanovic, Mala Mateen, Jean-Francois Sauvage, Benoit Neichel

Abstract: The next generation of giant ground and space telescopes will have the light-collecting power to detect and characterize potentially habitable terrestrial exoplanets using high-contrast imaging for the first time. This will only be achievable if the performance of Giant Segmented Mirror Telescopes (GSMTs) extreme adaptive optics (ExAO) systems are optimized to their full potential. A key component… ▽ More The next generation of giant ground and space telescopes will have the light-collecting power to detect and characterize potentially habitable terrestrial exoplanets using high-contrast imaging for the first time. This will only be achievable if the performance of Giant Segmented Mirror Telescopes (GSMTs) extreme adaptive optics (ExAO) systems are optimized to their full potential. A key component of an ExAO system is the wavefront sensor (WFS), which measures aberrations from atmospheric turbulence. A common choice in current and next-generation instruments is the pyramid wavefront sensor (PWFS). ExAO systems require high spatial and temporal sampling of wavefronts to optimize performance, and as a result, require large detectors for the WFS. We present a closed-loop testbed demonstration of a three-sided pyramid wavefront sensor (3PWFS) as an alternative to the conventional four-sided pyramid wavefront (4PWFS) sensor for GSMT-ExAO applications on the new Comprehensive Adaptive Optics and Coronagraph Test Instrument (CACTI). The 3PWFS is less sensitive to read noise than the 4PWFS because it uses fewer detector pixels. The 3PWFS has further benefits: a high-quality three-sided pyramid optic is easier to manufacture than a four-sided pyramid. We detail the design of the two components of the CACTI system, the adaptive optics simulator and the PWFS testbed that includes both a 3PWFS and 4PWFS. A preliminary experiment was performed on CACTI to study the performance of the 3PWFS to the 4PWFS in varying strengths of turbulence using both the Raw Intensity and Slopes Map signal processing methods. This experiment was repeated for a modulation radius of 1.6 lambda/D and 3.25 lambda/D. We found that the performance of the two wavefront sensors is comparable if modal loop gains are tuned. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: 28 Pages, 15 Figures, and 4 Tables

A Scalable Cryogenic LED Module for Selectively Illuminating Kinetic Inductance Detector Arrays

Authors: Jordan E. Shroyer, Matt Nelson, Liam Walters, Bradley R. Johnson

Abstract: We present the design and measured performance of a light emitting diode (LED) module for spatially mapping kinetic inductance detector (KID) arrays in the laboratory. Our novel approach uses a multiplexing scheme that only requires seven wires to control 480 red LEDs, and the number of LEDs can be scaled up without adding any additional wires. This multiplexing approach relies on active surface m… ▽ More We present the design and measured performance of a light emitting diode (LED) module for spatially mapping kinetic inductance detector (KID) arrays in the laboratory. Our novel approach uses a multiplexing scheme that only requires seven wires to control 480 red LEDs, and the number of LEDs can be scaled up without adding any additional wires. This multiplexing approach relies on active surface mount components that can operate at cryogenic temperatures down to 10 K. Cryogenic tests in liquid nitrogen and inside our cryostat demonstrate that the multiplexer circuit works at 77 and 10 K, respectively. The LED module presented here is tailored for our millimeter-wave detector modules, but the approach could be adapted for use with other KID-based detector systems. △ Less

Submitted 11 January, 2023; v1 submitted 20 June, 2022; originally announced June 2022.

Comments: Published in Review of Scientific Instruments as of 15 November 2022

Journal ref: Review of Scientific Instruments 93, 113107 (2022)

Separation of track- and shower-like energy deposits in ProtoDUNE-SP using a convolutional neural network

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Aimard, B. Ali-Mohammadzadeh, T. Alion, K. Allison, S. Alonso Monsalve, M. AlRashed, C. Alt, A. Alton, R. Alvarez, P. Amedo, J. Anderson , et al. (1204 additional authors not shown)

Abstract: Liquid argon time projection chamber detector technology provides high spatial and calorimetric resolutions on the charged particles traversing liquid argon. As a result, the technology has been used in a number of recent neutrino experiments, and is the technology of choice for the Deep Underground Neutrino Experiment (DUNE). In order to perform high precision measurements of neutrinos in the det… ▽ More Liquid argon time projection chamber detector technology provides high spatial and calorimetric resolutions on the charged particles traversing liquid argon. As a result, the technology has been used in a number of recent neutrino experiments, and is the technology of choice for the Deep Underground Neutrino Experiment (DUNE). In order to perform high precision measurements of neutrinos in the detector, final state particles need to be effectively identified, and their energy accurately reconstructed. This article proposes an algorithm based on a convolutional neural network to perform the classification of energy deposits and reconstructed particles as track-like or arising from electromagnetic cascades. Results from testing the algorithm on data from ProtoDUNE-SP, a prototype of the DUNE far detector, are presented. The network identifies track- and shower-like particles, as well as Michel electrons, with high efficiency. The performance of the algorithm is consistent between data and simulation. △ Less

Submitted 30 June, 2022; v1 submitted 31 March, 2022; originally announced March 2022.

Comments: 31 pages, 15 figures

Report number: FERMILAB-PUB-22-240-AD-ESH-LBNF-ND-SCD, CERN-EP-2022-077

Journal ref: Eur.Phys.J.C 82 (2022) 10, 903

Scintillation light detection in the 6-m drift-length ProtoDUNE Dual Phase liquid argon TPC

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Aimard, B. Ali-Mohammadzadeh, T. Alion, K. Allison, S. Alonso Monsalve, M. AlRashed, C. Alt, A. Alton, R. Alvarez, P. Amedo, J. Anderson , et al. (1202 additional authors not shown)

Abstract: DUNE is a dual-site experiment for long-baseline neutrino oscillation studies, neutrino astrophysics and nucleon decay searches. ProtoDUNE Dual Phase (DP) is a 6x6x6m3 liquid argon time-projection-chamber (LArTPC) that recorded cosmic-muon data at the CERN Neutrino Platform in 2019-2020 as a prototype of the DUNE Far Detector. Charged particles propagating through the LArTPC produce ionization and… ▽ More DUNE is a dual-site experiment for long-baseline neutrino oscillation studies, neutrino astrophysics and nucleon decay searches. ProtoDUNE Dual Phase (DP) is a 6x6x6m3 liquid argon time-projection-chamber (LArTPC) that recorded cosmic-muon data at the CERN Neutrino Platform in 2019-2020 as a prototype of the DUNE Far Detector. Charged particles propagating through the LArTPC produce ionization and scintillation light. The scintillation light signal in these detectors can provide the trigger for non-beam events. In addition, it adds precise timing capabilities and improves the calorimetry measurements. In ProtoDUNE-DP, scintillation and electroluminescence light produced by cosmic muons in the LArTPC is collected by photomultiplier tubes placed up to 7 m away from the ionizing track. In this paper, the ProtoDUNE-DP photon detection system performance is evaluated with a particular focus on the different wavelength shifters, such as PEN and TPB, and the use of Xe-doped LAr, considering its future use in giant LArTPCs. The scintillation light production and propagation processes are analyzed and a comparison of simulation to data is performed, improving understanding of the liquid argon properties △ Less

Submitted 3 June, 2022; v1 submitted 30 March, 2022; originally announced March 2022.

Comments: 31 pages, 29 figures

Report number: CERN-EP-DRAFT-MISC-2022-003; FERMILAB-PUB-22-242-LBNF

Observation of Radon Mitigation in MicroBooNE by a Liquid Argon Filtration System

Authors: MicroBooNE collaboration, P. Abratenko, J. Anthony, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, C. Barnes, G. Barr, J. Barrow, V. Basque, L. Bathe-Peters, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bhattacharya, M. Bishai, A. Blake, T. Bolton, J. Y. Book, L. Camilleri, D. Caratelli, I. Caro Terrazas , et al. (168 additional authors not shown)

Abstract: The MicroBooNE liquid argon time projection chamber (LArTPC) maintains a high level of liquid argon purity through the use of a filtration system that removes electronegative contaminants in continuously-circulated liquid, recondensed boil off, and externally supplied argon gas. We use the MicroBooNE LArTPC to reconstruct MeV-scale radiological decays. Using this technique we measure the liquid ar… ▽ More The MicroBooNE liquid argon time projection chamber (LArTPC) maintains a high level of liquid argon purity through the use of a filtration system that removes electronegative contaminants in continuously-circulated liquid, recondensed boil off, and externally supplied argon gas. We use the MicroBooNE LArTPC to reconstruct MeV-scale radiological decays. Using this technique we measure the liquid argon filtration system's efficacy at removing radon. This is studied by placing a 500 kBq $^{222}$Rn source upstream of the filters and searching for a time-dependent increase in the number of radiological decays in the LArTPC. In the context of two models for radon mitigation via a liquid argon filtration system, a slowing mechanism and a trapping mechanism, MicroBooNE data supports a radon reduction factor of greater than 99.999% or 97%, respectively. Furthermore, a radiological survey of the filters found that the copper-based filter material was the primary medium that removed the $^{222}$Rn. This is the first observation of radon mitigation in liquid argon with a large-scale copper-based filter and could offer a radon mitigation solution for future large LArTPCs. △ Less

Submitted 26 October, 2022; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: 14 pages, 7 figures, accepted by the Journal of Instrumentation

Report number: FERMILAB-PUB-22-203-ND

Journal ref: JINST 17 P11022 (2022)

Celeritas: GPU-accelerated particle transport for detector simulation in High Energy Physics experiments

Authors: S. C. Tognini, P. Canal, T. M. Evans, G. Lima, A. L. Lund, S. R. Johnson, S. Y. Jun, V. R. Pascuzzi, P. K. Romano

Abstract: Within the next decade, experimental High Energy Physics (HEP) will enter a new era of scientific discovery through a set of targeted programs recommended by the Particle Physics Project Prioritization Panel (P5), including the upcoming High Luminosity Large Hadron Collider (LHC) HL-LHC upgrade and the Deep Underground Neutrino Experiment (DUNE). These efforts in the Energy and Intensity Frontiers… ▽ More Within the next decade, experimental High Energy Physics (HEP) will enter a new era of scientific discovery through a set of targeted programs recommended by the Particle Physics Project Prioritization Panel (P5), including the upcoming High Luminosity Large Hadron Collider (LHC) HL-LHC upgrade and the Deep Underground Neutrino Experiment (DUNE). These efforts in the Energy and Intensity Frontiers will require an unprecedented amount of computational capacity on many fronts including Monte Carlo (MC) detector simulation. In order to alleviate this impending computational bottleneck, the Celeritas MC particle transport code is designed to leverage the new generation of heterogeneous computer architectures, including the exascale computing power of U.S. Department of Energy (DOE) Leadership Computing Facilities (LCFs), to model targeted HEP detector problems at the full fidelity of Geant4. This paper presents the planned roadmap for Celeritas, including its proposed code architecture, physics capabilities, and strategies for integrating it with existing and future experimental HEP computing workflows. △ Less

Submitted 22 March, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021

Report number: FERMILAB-FN-1159-SCD

Detector and Beamline Simulation for Next-Generation High Energy Physics Experiments

Authors: Sunanda Banerjee, D. N. Brown, David N. Brown, Paolo Calafiura, Jacob Calcutt, Philippe Canal, Miriam Diamond, Daniel Elvira, Thomas Evans, Renee Fatemi, Krzysztof Genser, Robert Hatcher, Alexander Himmel, Seth R. Johnson, Soon Yung Jun, Michael Kelsey, Evangelos Kourlitis, Robert K. Kutschke, Guilherme Lima, Kevin Lynch, Kendall Mahn, Zachary Marshall, Michael Mooney, Adam Para, Vincent R. Pascuzzi , et al. (9 additional authors not shown)

Abstract: The success of high energy physics programs relies heavily on accurate detector simulations and beam interaction modeling. The increasingly complex detector geometries and beam dynamics require sophisticated techniques in order to meet the demands of current and future experiments. Common software tools used today are unable to fully utilize modern computational resources, while data-recording rat… ▽ More The success of high energy physics programs relies heavily on accurate detector simulations and beam interaction modeling. The increasingly complex detector geometries and beam dynamics require sophisticated techniques in order to meet the demands of current and future experiments. Common software tools used today are unable to fully utilize modern computational resources, while data-recording rates are often orders of magnitude larger than what can be produced via simulation. In this paper, we describe the state, current and future needs of high energy physics detector and beamline simulations and related challenges, and we propose a number of possible ways to address them. △ Less

Submitted 20 April, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021

Report number: FERMILAB-FN-1151-ND-PPD-SCD

The Simons Observatory: Design and Measured Performance of a Carbon Fiber Strut for a Cryogenic Truss

Authors: Kevin D. Crowley, Peter Dow, Jordan E. Shroyer, John C. Groh, Bradley Dober, Jacob Spisak, Nicholas Galitzki, Tanay Bhandarkar, Mark J. Devlin, Simon Dicker, Patricio A. Gallardo, Kathleen Harrington, Bradley R. Johnson, Delwin Johnson, Anna M. Kofman, Akito Kusaka, Adrian Lee, Michele Limon, Jeffrey Iuliano, Federico Nati, John Orlowski-Scherer, Lyman Page, Michael Randall, Grant Teply, Tran Tsan , et al. (3 additional authors not shown)

Abstract: We present the design and measured performance of a new carbon fiber strut design that is used in a cryogenically cooled truss for the Simons Observatory Small Aperture Telescope (SAT). The truss consists of two aluminum 6061 rings separated by 24 struts. Each strut consists of a central carbon fiber tube fitted with two aluminum end caps. We tested the performance of the strut and truss by (i) cr… ▽ More We present the design and measured performance of a new carbon fiber strut design that is used in a cryogenically cooled truss for the Simons Observatory Small Aperture Telescope (SAT). The truss consists of two aluminum 6061 rings separated by 24 struts. Each strut consists of a central carbon fiber tube fitted with two aluminum end caps. We tested the performance of the strut and truss by (i) cryogenically cycling and destructively pull-testing strut samples, (ii) non-destructively pull-testing the final truss, and (iii) measuring the thermal conductivity of the carbon fiber tubes. We found that the strut strength is limited by the mounting fasteners and the strut end caps, not the epoxy adhesive or the carbon fiber tube. This result is consistent with our numerical predictions. Our thermal measurements suggest that the conductive heat load through the struts (from 4 K to 1 K) will be less than 1 mW. This strut design may be a promising candidate for use in other cryogenic support structures. △ Less

Submitted 18 January, 2022; v1 submitted 16 January, 2022; originally announced January 2022.

Journal ref: Review of Scientific Instruments 93, 055106 (2022)

Laser guide star return-flux gain from frequency chirping

Authors: J. Hellemeier, M. Enderlein, M. Hager, D. Bonaccini Calia, R. L. Johnson, F. Lison, M. O. Byrd, L. A. Kann, M. Centrone, P. Hickson

Abstract: Spectral hole burning reduces sodium laser guide star efficiency. Due to photon recoil, atoms that are initially resonant with the single-frequency laser get Doppler shifted out of resonance, which reduces the return flux. Frequency-chirped (also known as frequency-swept) continuous-wave lasers have the potential to mitigate the effect of spectral hole burning and even increase the laser guide sta… ▽ More Spectral hole burning reduces sodium laser guide star efficiency. Due to photon recoil, atoms that are initially resonant with the single-frequency laser get Doppler shifted out of resonance, which reduces the return flux. Frequency-chirped (also known as frequency-swept) continuous-wave lasers have the potential to mitigate the effect of spectral hole burning and even increase the laser guide star efficiency beyond the theoretical limit of a single-frequency laser. On-sky measurements of a frequency-chirped, single-frequency laser guide star are performed at the Roque de los Muchachos Observatory on La Palma. In the experiment, a 35-cm telescope and a fast photon counting receiver system are employed to resolve the return flux response during laser frequency sweeps gaining insights into the population dynamics of the sodium layer. At a launched laser power of 16.5 W, we find a maximum gain in return flux of 22\% compared to a fixed-frequency laser. Our results suggest a strong dependence of chirping gain on power density at the mesosphere, i.e. laser power and seeing. Maximum gains are recorded at a chirping amplitude on the order of 150 MHz and a chirping rate of 0.8 MHz $μ$s$^{-1}$, as predicted by theory. Time-resolved measurements during the chirping period confirm our understanding of the population dynamics in the sodium layer. To our knowledge these are the first measurements of return flux enhancement for laser guide stars excited by a single frequency-chirped continuous-wave laser. For higher laser powers, the effectiveness of chirping is expected to increase, which could be highly beneficial for telescopes equipped with high-power laser guide star adaptive optics systems, also for emerging space awareness applications using laser guide stars such as satellite imaging and ground-to-space optical communications. △ Less

Submitted 15 December, 2021; v1 submitted 30 November, 2021; originally announced December 2021.

The Simons Observatory: Magnetic Shielding Measurements for the Universal Multiplexing Module

Authors: Zachary B. Huber, Yaqiong Li, Eve M. Vavagiakis, Steve K. Choi, Jake Connors, Nicholas F. Cothard, Cody J. Duell, Nicholas Galitzki, Erin Healy, Johannes Hubmayr, Bradley R. Johnson, Benjamin Keller, Heather McCarrick, Michael D. Niemack, Yuhan Wang, Zhilei Xu, Kaiwen Zheng

Abstract: The Simons Observatory (SO) includes four telescopes that will measure the temperature and polarization of the cosmic microwave background using over 60,000 highly sensitive transition-edge bolometers (TES). These multichroic TES bolometers are read out by a microwave RF SQUID multiplexing system with a multiplexing factor of 910. Given that both TESes and SQUIDs are susceptible to magnetic field… ▽ More The Simons Observatory (SO) includes four telescopes that will measure the temperature and polarization of the cosmic microwave background using over 60,000 highly sensitive transition-edge bolometers (TES). These multichroic TES bolometers are read out by a microwave RF SQUID multiplexing system with a multiplexing factor of 910. Given that both TESes and SQUIDs are susceptible to magnetic field pickup and that it is hard to predict how they will respond to such fields, it is important to characterize the magnetic response of these systems empirically. This information can then be used to limit spurious signals by informing magnetic shielding designs for the detectors and readout. This paper focuses on measurements of magnetic pickup with different magnetic shielding configurations for the SO universal multiplexing module (UMM), which contains the SQUIDs, associated resonators, and TES bias circuit. The magnetic pickup of a prototype UMM was tested under three shielding configurations: no shielding (copper packaging), aluminum packaging for the UMM, and a tin/lead-plated shield surrounding the entire dilution refrigerator 100 mK cold stage. The measurements show that the aluminum packaging outperforms the copper packaging by a shielding factor of 8-10, and adding the tin/lead-plated 1K shield further increases the relative shielding factor in the aluminum configuration by 1-2 orders of magnitude. △ Less

Submitted 1 March, 2023; v1 submitted 22 November, 2021; originally announced November 2021.

Comments: 7 pages, 4 figure, conference proceedings submitted to the Journal of Low Temperature Physics; includes updates in response to reviewer comments

Journal ref: J Low Temp Phys 209, 667-676 (2022)

Simons Observatory Focal-Plane Module: In-lab Testing and Characterization Program

Authors: Yuhan Wang, Kaiwen Zheng, Zachary Atkins, Jason Austermann, Tanay Bhandarkar, Steve K. Choi, Shannon M. Duff, Daniel Dutcher, Nicholas Galitzki, Erin Healy, Zachary B. Huber, Johannes Hubmayr, Bradley R. Johnson, Jack Lashner, Yaqiong Li, Heather McCarrick, Michael D. Niemack, Joseph Seibert, Maximiliano Silva-Feaver, Rita Sonka, Suzanne T. Staggs, Eve Vavagiakis, Zhilei Xu

Abstract: The Simons Observatory (SO) is a ground-based cosmic microwave background instrument to be sited in the Atacama Desert in Chile. SO will deploy 60,000 transition-edge sensor bolometers in 49 separate focal-plane modules across a suite of four telescopes covering three dichroic bands termed low frequency (LF), mid frequency (MF) and ultra-high frequency (UHF). Each MF and UHF focal-plane module pac… ▽ More The Simons Observatory (SO) is a ground-based cosmic microwave background instrument to be sited in the Atacama Desert in Chile. SO will deploy 60,000 transition-edge sensor bolometers in 49 separate focal-plane modules across a suite of four telescopes covering three dichroic bands termed low frequency (LF), mid frequency (MF) and ultra-high frequency (UHF). Each MF and UHF focal-plane module packages 1720 optical detectors and corresponding 100 mK microwave SQUID multiplexing readout components. In this paper we describe the testing program we have developed for high-throughput validation of the modules after they are assembled. The validation requires measurements of the yield, saturation powers, time constants, noise properties and optical efficiencies. Additional measurements will be performed for further characterizations as needed. We describe the methods developed and demonstrate preliminary results from initial testing of prototype modules. △ Less

Submitted 5 July, 2022; v1 submitted 22 November, 2021; originally announced November 2021.

Journal ref: J Low Temp Phys 209 (2022), 944-952

Novel Approach for Evaluating Detector-Related Uncertainties in a LArTPC Using MicroBooNE Data

Authors: MicroBooNE collaboration, P. Abratenko, R. An, J. Anthony, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, C. Barnes, G. Barr, V. Basque, L. Bathe-Peters, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bishai, A. Blake, T. Bolton, J. Y. Book, L. Camilleri, D. Caratelli, I. Caro Terrazas, F. Cavanna , et al. (161 additional authors not shown)

Abstract: Primary challenges for current and future precision neutrino experiments using liquid argon time projection chambers (LArTPCs) include understanding detector effects and quantifying the associated systematic uncertainties. This paper presents a novel technique for assessing and propagating LArTPC detector-related systematic uncertainties. The technique makes modifications to simulation waveforms b… ▽ More Primary challenges for current and future precision neutrino experiments using liquid argon time projection chambers (LArTPCs) include understanding detector effects and quantifying the associated systematic uncertainties. This paper presents a novel technique for assessing and propagating LArTPC detector-related systematic uncertainties. The technique makes modifications to simulation waveforms based on a parameterization of observed differences in ionization signals from the TPC between data and simulation, while remaining insensitive to the details of the detector model. The modifications are then used to quantify the systematic differences in low- and high-level reconstructed quantities. This approach could be applied to future LArTPC detectors, such as those used in SBN and DUNE. △ Less

Submitted 16 June, 2022; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: 14 pages, 12 figures, version accepted by EPJ C

Report number: FERMILAB-PUB-21-501-ND-SCD

Journal ref: Eur. Phys. J. C 82, 454 (2022)

Wire-Cell 3D Pattern Recognition Techniques for Neutrino Event Reconstruction in Large LArTPCs: Algorithm Description and Quantitative Evaluation with MicroBooNE Simulation

Authors: MicroBooNE collaboration, P. Abratenko, R. An, J. Anthony, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, C. Barnes, G. Barr, V. Basque, L. Bathe-Peters, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bishai, A. Blake, T. Bolton, J. Y. Book, L. Camilleri, D. Caratelli, I. Caro Terrazas, R. Castillo Fernandez , et al. (163 additional authors not shown)

Abstract: Wire-Cell is a 3D event reconstruction package for liquid argon time projection chambers. Through geometry, time, and drifted charge from multiple readout wire planes, 3D space points with associated charge are reconstructed prior to the pattern recognition stage. Pattern recognition techniques, including track trajectory and $dQ/dx$ (ionization charge per unit length) fitting, 3D neutrino vertex… ▽ More Wire-Cell is a 3D event reconstruction package for liquid argon time projection chambers. Through geometry, time, and drifted charge from multiple readout wire planes, 3D space points with associated charge are reconstructed prior to the pattern recognition stage. Pattern recognition techniques, including track trajectory and $dQ/dx$ (ionization charge per unit length) fitting, 3D neutrino vertex fitting, track and shower separation, particle-level clustering, and particle identification are then applied on these 3D space points as well as the original 2D projection measurements. A deep neural network is developed to enhance the reconstruction of the neutrino interaction vertex. Compared to traditional algorithms, the deep neural network boosts the vertex efficiency by a relative 30\% for charged-current $ν_e$ interactions. This pattern recognition achieves 80-90\% reconstruction efficiencies for primary leptons, after a 65.8\% (72.9\%) vertex efficiency for charged-current $ν_e$ ($ν_μ$) interactions. Based on the resulting reconstructed particles and their kinematics, we also achieve 15-20\% energy reconstruction resolutions for charged-current neutrino interactions. △ Less

Submitted 26 December, 2021; v1 submitted 26 October, 2021; originally announced October 2021.

Report number: FERMILAB-PUB-21-509-ND

Exploratory Lagrangian-Based Particle Tracing Using Deep Learning

Authors: Mengjiao Han, Sudhanshu Sane, Chris R. Johnson

Abstract: Time-varying vector fields produced by computational fluid dynamics simulations are often prohibitively large and pose challenges for accurate interactive analysis and exploration. To address these challenges, reduced Lagrangian representations have been increasingly researched as a means to improve scientific time-varying vector field exploration capabilities. This paper presents a novel deep neu… ▽ More Time-varying vector fields produced by computational fluid dynamics simulations are often prohibitively large and pose challenges for accurate interactive analysis and exploration. To address these challenges, reduced Lagrangian representations have been increasingly researched as a means to improve scientific time-varying vector field exploration capabilities. This paper presents a novel deep neural network-based particle tracing method to explore time-varying vector fields represented by Lagrangian flow maps. In our workflow, in situ processing is first utilized to extract Lagrangian flow maps, and deep neural networks then use the extracted data to learn flow field behavior. Using a trained model to predict new particle trajectories offers a fixed small memory footprint and fast inference. To demonstrate and evaluate the proposed method, we perform an in-depth study of performance using a well-known analytical data set, the Double Gyre. Our study considers two flow map extraction strategies as well as the impact of the number of training samples and integration durations on efficacy, evaluates multiple sampling options for training and testing and informs hyperparameter settings. Overall, we find our method requires a fixed memory footprint of 10.5 MB to encode a Lagrangian representation of a time-varying vector field while maintaining accuracy. For post hoc analysis, loading the trained model costs only two seconds, significantly reducing the burden of I/O when reading data for visualization. Moreover, our parallel implementation can infer one hundred locations for each of two thousand new pathlines across the entire temporal resolution in 1.3 seconds using one NVIDIA Titan RTX GPU. △ Less

Submitted 6 January, 2022; v1 submitted 15 October, 2021; originally announced October 2021.

Comments: The paper has been accepted to publish by Journal of Flow Visualization and Image Processing

On the Phase Grouping Mechanism for a Magnetron Coherent Oscillation

Authors: G. Kazakevich, R. P. Johnson, V. Yakovlev

Abstract: CW magnetrons, developed for industrial heaters, but driven by an injection-locking signal were suggested to power Superconducting RF (SRF) cavities due to higher efficiency and lower cost of generated RF power per Watt than traditionally used RF sources (klystrons, IOTs, solid-state amplifiers). When the magnetrons are intended to feed Room Temperature (RT) cavities, the injected phase or frequen… ▽ More CW magnetrons, developed for industrial heaters, but driven by an injection-locking signal were suggested to power Superconducting RF (SRF) cavities due to higher efficiency and lower cost of generated RF power per Watt than traditionally used RF sources (klystrons, IOTs, solid-state amplifiers). When the magnetrons are intended to feed Room Temperature (RT) cavities, the injected phase or frequency locking signal may provide required phase or frequency stability of the accelerating field. However, when the magnetron RF sources are intended to feed high Q-factor SRF cavities, the sources must be controlled in phase and power in a wide bandwidth to compensate parasitic phase and amplitude modulations caused by microphonics. In dependence on parameters of magnetron and the injection-locking signal one can choose regime most suitable for feeding SRF cavities, enabling magnetron almost coherent oscillation at the wide bandwidth of control. A novel approach considering magnetrons as quasi-coherent or coherent RF generators enables choosing the tube parameters and operation most suitable for various SRF accelerators. △ Less

Submitted 15 September, 2021; originally announced September 2021.

Comments: 6 pages, 7 figures

First Measurement of Inclusive Electron-Neutrino and Antineutrino Charged Current Differential Cross Sections in Charged Lepton Energy on Argon in MicroBooNE

Authors: MicroBooNE collaboration, P. Abratenko, R. An, J. Anthony, L. Arellano, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, C. Barnes, G. Barr, V. Basque, L. Bathe-Peters, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bishai, A. Blake, T. Bolton, J. Y. Book, L. Camilleri, D. Caratelli, I. Caro Terrazas, R. Castillo Fernandez , et al. (163 additional authors not shown)

Abstract: We present the first measurement of the single-differential $ν_e + \barν_e$ charged-current inclusive cross sections on argon in electron or positron energy and in electron or positron scattering cosine over the full angular range. Data were collected using the MicroBooNE liquid argon time projection chamber located off-axis from the Fermilab Neutrinos at the Main Injector beam over an exposure of… ▽ More We present the first measurement of the single-differential $ν_e + \barν_e$ charged-current inclusive cross sections on argon in electron or positron energy and in electron or positron scattering cosine over the full angular range. Data were collected using the MicroBooNE liquid argon time projection chamber located off-axis from the Fermilab Neutrinos at the Main Injector beam over an exposure of $2.0\times10^{20}$ protons on target. The signal definition includes a 60 MeV threshold on the $ν_e$ or $\barν_e$ energy and a 120 MeV threshold on the electron or positron energy. The measured total and differential cross sections are found to be in agreement with the GENIE, NuWro, and GiBUU neutrino generators. △ Less

Submitted 3 February, 2022; v1 submitted 14 September, 2021; originally announced September 2021.

Comments: 7 Pages, 2 figures. Supplemental materials include cross section values, efficiencies, purity, flux, fractional uncertainties, additional smearing matrices, and unfolded covariance matrices

Report number: FERMILAB-PUB-21-443-ND

Journal ref: Phys. Rev. D 105, L051102 (2022)

Calorimetric classification of track-like signatures in liquid argon TPCs using MicroBooNE data

Authors: MicroBooNE collaboration, P. Abratenko, R. An, J. Anthony, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, C. Barnes, G. Barr, V. Basque, L. Bathe-Peters, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bishai, A. Blake, T. Bolton, L. Camilleri, D. Caratelli, I. Caro Terrazas, R. Castillo Fernandez, F. Cavanna, G. Cerati , et al. (157 additional authors not shown)

Abstract: The MicroBooNE liquid argon time projection chamber located at Fermilab is a neutrino experiment dedicated to the study of short-baseline oscillations, the measurements of neutrino cross sections in liquid argon, and to the research and development of this novel detector technology. Accurate and precise measurements of calorimetry are essential to the event reconstruction and are achieved by lever… ▽ More The MicroBooNE liquid argon time projection chamber located at Fermilab is a neutrino experiment dedicated to the study of short-baseline oscillations, the measurements of neutrino cross sections in liquid argon, and to the research and development of this novel detector technology. Accurate and precise measurements of calorimetry are essential to the event reconstruction and are achieved by leveraging the TPC to measure deposited energy per unit length along the particle trajectory, with mm resolution. We describe the non-uniform calorimetric reconstruction performance in the detector, showing dependence on the angle of the particle trajectory. Such non-uniform reconstruction directly affects the performance of the particle identification algorithms which infer particle type from calorimetric measurements. This work presents a new particle identification method which accounts for and effectively addresses such non-uniformity. The newly developed method shows improved performance compared to previous algorithms, illustrated by a 94% proton selection efficiency and a 10% muon mis-identification rate, with a fairly loose selection of tracks performed on beam data. The performance is further demonstrated by identifying exclusive final states in $ν_μ CC$ interactions. While developed using MicroBooNE data and simulation, this method is easily applicable to future LArTPC experiments, such as SBND, ICARUS, and DUNE. △ Less

Submitted 4 January, 2022; v1 submitted 31 August, 2021; originally announced September 2021.

Comments: 17 pages, 9 figures The updated version contains a clearer fig 1, some better quantification of physics reach in section 6.3, while several typos have been fixed

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

Journal ref: Journal of High Energy Physics volume 2021, Article number: 153 (2021)

Low exposure long-baseline neutrino oscillation sensitivity of the DUNE experiment

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, D. Adams, M. Adinolfi, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Aimard, B. Ali-Mohammadzadeh, T. Alion, K. Allison, S. Alonso Monsalve, M. AlRashed, C. Alt, A. Alton, P. Amedo, J. Anderson, C. Andreopoulos, M. Andreotti , et al. (1132 additional authors not shown)

Abstract: The Deep Underground Neutrino Experiment (DUNE) will produce world-leading neutrino oscillation measurements over the lifetime of the experiment. In this work, we explore DUNE's sensitivity to observe charge-parity violation (CPV) in the neutrino sector, and to resolve the mass ordering, for exposures of up to 100 kiloton-megawatt-years (kt-MW-yr). The analysis includes detailed uncertainties on t… ▽ More The Deep Underground Neutrino Experiment (DUNE) will produce world-leading neutrino oscillation measurements over the lifetime of the experiment. In this work, we explore DUNE's sensitivity to observe charge-parity violation (CPV) in the neutrino sector, and to resolve the mass ordering, for exposures of up to 100 kiloton-megawatt-years (kt-MW-yr). The analysis includes detailed uncertainties on the flux prediction, the neutrino interaction model, and detector effects. We demonstrate that DUNE will be able to unambiguously resolve the neutrino mass ordering at a 3$σ$ (5$σ$) level, with a 66 (100) kt-MW-yr far detector exposure, and has the ability to make strong statements at significantly shorter exposures depending on the true value of other oscillation parameters. We also show that DUNE has the potential to make a robust measurement of CPV at a 3$σ$ level with a 100 kt-MW-yr exposure for the maximally CP-violating values $δ_{\rm CP}} = \pmπ/2$. Additionally, the dependence of DUNE's sensitivity on the exposure taken in neutrino-enhanced and antineutrino-enhanced running is discussed. An equal fraction of exposure taken in each beam mode is found to be close to optimal when considered over the entire space of interest. △ Less

Submitted 3 September, 2021; originally announced September 2021.

Report number: FERMILAB-PUB-21-391-ND

Design, construction and operation of the ProtoDUNE-SP Liquid Argon TPC

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, D. Adams, M. Adinolfi, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Ali-Mohammadzadeh, T. Alion, K. Allison, S. Alonso Monsalve, M. Alrashed, C. Alt, A. Alton, P. Amedo, J. Anderson, C. Andreopoulos, M. Andreotti, M. P. Andrews , et al. (1158 additional authors not shown)

Abstract: The ProtoDUNE-SP detector is a single-phase liquid argon time projection chamber (LArTPC) that was constructed and operated in the CERN North Area at the end of the H4 beamline. This detector is a prototype for the first far detector module of the Deep Underground Neutrino Experiment (DUNE), which will be constructed at the Sandford Underground Research Facility (SURF) in Lead, South Dakota, USA.… ▽ More The ProtoDUNE-SP detector is a single-phase liquid argon time projection chamber (LArTPC) that was constructed and operated in the CERN North Area at the end of the H4 beamline. This detector is a prototype for the first far detector module of the Deep Underground Neutrino Experiment (DUNE), which will be constructed at the Sandford Underground Research Facility (SURF) in Lead, South Dakota, USA. The ProtoDUNE-SP detector incorporates full-size components as designed for DUNE and has an active volume of $7\times 6\times 7.2$~m$^3$. The H4 beam delivers incident particles with well-measured momenta and high-purity particle identification. ProtoDUNE-SP's successful operation between 2018 and 2020 demonstrates the effectiveness of the single-phase far detector design. This paper describes the design, construction, assembly and operation of the detector components. △ Less

Submitted 23 September, 2021; v1 submitted 4 August, 2021; originally announced August 2021.

Building a Quantum Engineering Undergraduate Program

Authors: Abraham Asfaw, Alexandre Blais, Kenneth R. Brown, Jonathan Candelaria, Christopher Cantwell, Lincoln D. Carr, Joshua Combes, Dripto M. Debroy, John M. Donohue, Sophia E. Economou, Emily Edwards, Michael F. J. Fox, Steven M. Girvin, Alan Ho, Hilary M. Hurst, Zubin Jacob, Blake R. Johnson, Ezekiel Johnston-Halperin, Robert Joynt, Eliot Kapit, Judith Klein-Seetharaman, Martin Laforest, H. J. Lewandowski, Theresa W. Lynn, Corey Rae H. McRae , et al. (12 additional authors not shown)

Abstract: The rapidly growing quantum information science and engineering (QISE) industry will require both quantum-aware and quantum-proficient engineers at the bachelor's level. We provide a roadmap for building a quantum engineering education program to satisfy this need. For quantum-aware engineers, we describe how to design a first quantum engineering course accessible to all STEM students. For the edu… ▽ More The rapidly growing quantum information science and engineering (QISE) industry will require both quantum-aware and quantum-proficient engineers at the bachelor's level. We provide a roadmap for building a quantum engineering education program to satisfy this need. For quantum-aware engineers, we describe how to design a first quantum engineering course accessible to all STEM students. For the education and training of quantum-proficient engineers, we detail both a quantum engineering minor accessible to all STEM majors, and a quantum track directly integrated into individual engineering majors. We propose that such programs typically require only three or four newly developed courses that complement existing engineering and science classes available on most larger campuses. We describe a conceptual quantum information science course for implementation at any post-secondary institution, including community colleges and military schools. QISE presents extraordinary opportunities to work towards rectifying issues of inclusivity and equity that continue to be pervasive within engineering. We present a plan to do so and describe how quantum engineering education presents an excellent set of education research opportunities. Finally, we outline a hands-on training plan on quantum hardware, a key component of any quantum engineering program, with a variety of technologies including optics, atoms and ions, cryogenic and solid-state technologies, nanofabrication, and control and readout electronics. Our recommendations provide a flexible framework that can be tailored for academic institutions ranging from teaching and undergraduate-focused two- and four-year colleges to research-intensive universities. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Comments: 25 pages, 2 figures

Journal ref: IEEE Transactions on Education 65, 220 (2022)

The influence of upper boundary conditions on molecular kinetic atmospheric escape simulations

Authors: S. R. Carberry Mogan, O. J. Tucker, R. E. Johnson

Abstract: Molecular kinetic simulations are typically used to accurately describe the tenuous regions of the upper atmospheres on planetary bodies. These simulations track the motion of particles representing real atmospheric atoms and/or molecules subject to collisions, the object's gravity, and external influences. Because particles can end up in very large ballistic orbits, upper boundary conditions (UBC… ▽ More Molecular kinetic simulations are typically used to accurately describe the tenuous regions of the upper atmospheres on planetary bodies. These simulations track the motion of particles representing real atmospheric atoms and/or molecules subject to collisions, the object's gravity, and external influences. Because particles can end up in very large ballistic orbits, upper boundary conditions (UBC) are typically used to limit the domain size thereby reducing the time for the atmosphere to reach steady-state. In the absence of a clear altitude at which all molecules are removed, such as a Hill sphere, an often used condition is to choose an altitude at which collisions become infrequent so that particles on escape trajectories are removed. The remainder are then either specularly reflected back into the simulation domain or their ballistic trajectories are calculated analytically or explicitly tracked so they eventually re-enter the domain. Here we examine the effect of the choice of the UBC on the escape rate and the structure of the atmosphere near the nominal exobase in the convenient and frequently used 1D spherically symmetric approximation. Using Callisto as the example body, we show that the commonly used specular reflection UBC can lead to significant uncertainties when simulating a species with a lifetime comparable to or longer than a dynamical time scale, such as an overestimation of escape rates and an inflated exosphere. Therefore, although specular reflection is convenient, the molecular lifetimes and body's dynamical time scales need to be considered even when implementing the convenient 1D spherically symmetric simulations in order to accurately estimate the escape rate and the density and temperature structure in the transition regime. △ Less

Submitted 26 July, 2021; originally announced July 2021.

A tenuous, collisional atmosphere on Callisto

Authors: Shane Carberry Mogan, Orenthal Tucker, Robert Johnson, Audrey Vorburger, Andre Galli, Benoit Marchand, Angelo Tafuni, Sunil Kumar, Iskender Sahin, Katepalli Sreenivasan

Abstract: A simulation tool which utilizes parallel processing is developed to describe molecular kinetics in 2D, single-and multi-component atmospheres on Callisto. This expands on our previous study on the role of collisions in 1D atmospheres on Callisto composed of radiolytic products (Carberry Mogan et al., 2020) by implementing a temperature gradient from noon to midnight across Callisto's surface and… ▽ More A simulation tool which utilizes parallel processing is developed to describe molecular kinetics in 2D, single-and multi-component atmospheres on Callisto. This expands on our previous study on the role of collisions in 1D atmospheres on Callisto composed of radiolytic products (Carberry Mogan et al., 2020) by implementing a temperature gradient from noon to midnight across Callisto's surface and introducing sublimated water vapor. We compare single-species, ballistic and collisional O2, H2 and H2O atmospheres, as well as an O2+H2O atmosphere to 3-species atmospheres which contain H2 in varying amounts. Because the H2O vapor pressure is extremely sensitive to the surface temperatures, the density drops several order of magnitude with increasing distance from the subsolar point, and the flow transitions from collisional to ballistic accordingly. In an O2+H2O atmosphere the local temperatures are determined by H2O near the subsolar point and transition with increasing distance from the subsolar point to being determined by O2 When radiolytically produced H2 is not negligible in O2+H2O+H2 atmospheres, this much lighter molecule, with a scale height roughly an order of magnitude larger than that for the heavier species, can cool the local temperatures via collisions. In addition, if the H2 component is dense enough, particles originating on the day-side and precipitating into the night-side atmosphere deposit energy via collisions, which in turn heats the local atmosphere relative to the surface temperature. Finally, we discuss the potential implications of this study on the presence of H2 in Callisto's atmosphere and how the simulated densities correlate with expected detection thresholds at flyby altitudes of the proposed JUpiter ICy moons Explorer (JUICE) spacecraft. △ Less

Submitted 26 July, 2021; originally announced July 2021.

The Simons Observatory: HoloSim-ML: machine learning applied to the efficient analysis of radio holography measurements of complex optical systems

Authors: Grace E. Chesmore, Alexandre E. Adler, Nicholas F. Cothard, Nadia Dachlythra, Patricio A. Gallardo, Jon Gudmundsson, Bradley R. Johnson, Michele Limon, Jeff McMahon, Federico Nati, Michael D. Niemack, Giuseppe Puglisi, Sara M. Simon, Edward J. Wollack, Kevin Wolz, Zhilei Xu, Ningfeng Zhu

Abstract: Near-field radio holography is a common method for measuring and aligning mirror surfaces for millimeter and sub-millimeter telescopes. In instruments with more than a single mirror, degeneracies arise in the holography measurement, requiring multiple measurements and new fitting methods. We present HoloSim-ML, a Python code for beam simulation and analysis of radio holography data from complex op… ▽ More Near-field radio holography is a common method for measuring and aligning mirror surfaces for millimeter and sub-millimeter telescopes. In instruments with more than a single mirror, degeneracies arise in the holography measurement, requiring multiple measurements and new fitting methods. We present HoloSim-ML, a Python code for beam simulation and analysis of radio holography data from complex optical systems. This code uses machine learning to efficiently determine the position of hundreds of mirror adjusters on multiple mirrors with few micron accuracy. We apply this approach to the example of the Simons Observatory 6m telescope. △ Less

Submitted 5 October, 2021; v1 submitted 8 July, 2021; originally announced July 2021.

Comments: Software is publicly available at: https://github.com/McMahonCosmologyGroup/holosim-ml

Statistical Rendering for Visualization of Red Sea Eddy Simulation Data

Authors: Tushar M. Athawale, Alireza Entezari, Bei Wang, Chris R. Johnson

Abstract: Analyzing the effects of ocean eddies is important in oceanology for gaining insights into transport of energy and biogeochemical particles. We present an application of statistical visualization algorithms for the analysis of the Red Sea eddy simulation ensemble. Specifically, we demonstrate the applications of statistical volume rendering and statistical Morse complex summary maps to a velocity… ▽ More Analyzing the effects of ocean eddies is important in oceanology for gaining insights into transport of energy and biogeochemical particles. We present an application of statistical visualization algorithms for the analysis of the Red Sea eddy simulation ensemble. Specifically, we demonstrate the applications of statistical volume rendering and statistical Morse complex summary maps to a velocity magnitude field for studying the eddy positions in the flow dataset. In statistical volume rendering, we model per-voxel data uncertainty using noise models, such as parametric and nonparametric, and study the propagation of uncertainty into the volume rendering pipeline. In the statistical Morse complex summary maps, we derive histograms charactering uncertainty of gradient flow destinations to understand Morse complex topological variations across the ensemble. We demonstrate the utility of our statistical visualizations for an effective analysis of the potential eddy positions and their spatial uncertainty. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Role of the Solar Minimum in the Waiting Time Distribution Throughout the Heliosphere

Authors: Yosia I. Nurhan, Jay R. Johnson, Jonathan R. Homan, Simon Wing

Abstract: We explore the tail of various waiting time datasets of processes that follow a nonstationary Poisson distribution with a sinusoidal driver. Analytically, we find that the distribution of large waiting times of such processes can be described using a power law slope of -2.5. We show that this result applies more broadly to any nonstationary Poisson process driven periodically. Examples of such pro… ▽ More We explore the tail of various waiting time datasets of processes that follow a nonstationary Poisson distribution with a sinusoidal driver. Analytically, we find that the distribution of large waiting times of such processes can be described using a power law slope of -2.5. We show that this result applies more broadly to any nonstationary Poisson process driven periodically. Examples of such processes include solar flares, coronal mass ejections, geomagnetic storms, and substorms. We also discuss how the power law specifically relates to the behavior of driver near its minima. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Measurement of the Longitudinal Diffusion of Ionization Electrons in the MicroBooNE Detector

Authors: P. Abratenko, R. An, J. Anthony, J. Asaadi, A. Ashkenazi, S. Balasubramanian, B. Baller, C. Barnes, G. Barr, V. Basque, L. Bathe-Peters, O. Benevides Rodrigues, S. Berkman, A. Bhanderi, A. Bhat, M. Bishai, A. Blake, T. Bolton, L. Camilleri, D. Caratelli, I. Caro Terrazas, R. Castillo Fernandez, F. Cavanna, G. Cerati, Y. Chen , et al. (157 additional authors not shown)

Abstract: Accurate knowledge of electron transport properties is vital to understanding the information provided by liquid argon time projection chambers (LArTPCs). Ionization electron drift-lifetime, local electric field distortions caused by positive ion accumulation, and electron diffusion can all significantly impact the measured signal waveforms. This paper presents a measurement of the effective longi… ▽ More Accurate knowledge of electron transport properties is vital to understanding the information provided by liquid argon time projection chambers (LArTPCs). Ionization electron drift-lifetime, local electric field distortions caused by positive ion accumulation, and electron diffusion can all significantly impact the measured signal waveforms. This paper presents a measurement of the effective longitudinal electron diffusion coefficient, $D_L$, in MicroBooNE at the nominal electric field strength of 273.9 V/cm. Historically, this measurement has been made in LArTPC prototype detectors. This represents the first measurement in a large-scale (85 tonne active volume) LArTPC operating in a neutrino beam. This is the largest dataset ever used for this measurement. Using a sample of $\sim$70,000 through-going cosmic ray muon tracks tagged with MicroBooNE's cosmic ray tagger system, we measure $D_L = 3.74^{+0.28}_{-0.29}$ cm$^2$/s. △ Less

Submitted 25 June, 2021; v1 submitted 13 April, 2021; originally announced April 2021.

Comments: Changes following feedback from journal referee. Some wording changes, clarifications. Merged several figures to save on space