Showing 1–49 of 49 results for author: Bonilla, J

Measurement of Electron Neutrino and Antineutrino Cross Sections at Low Momentum Transfer

Authors: S. Henry, H. Su, S. Akhter, Z. Ahmad Dar, V. Ansari, M. V. Ascencio, M. Sajjad Athar, A. Bashyal, M. Betancourt, J. L. Bonilla, A. Bravar, G. Caceres, G. A. Díaz, J. Felix, L. Fields, R. Fine, P. K. Gaur, S. M. Gilligan, R. Gran, E. Granados, D. A. Harris, A. L. Hart, J. Kleykamp, A. Klustová, M. Kordosky , et al. (31 additional authors not shown)

Abstract: Accelerator based neutrino oscillation experiments seek to measure the relative number of electron and muon neutrinos and antineutrinos at different $L/E$ values. However high statistics studies of neutrino interactions are almost exclusively measured using muon neutrinos and antineutrinos since the dominant flavor of neutrinos produced by accelerator based beams are of the muon type. This work re… ▽ More Accelerator based neutrino oscillation experiments seek to measure the relative number of electron and muon neutrinos and antineutrinos at different $L/E$ values. However high statistics studies of neutrino interactions are almost exclusively measured using muon neutrinos and antineutrinos since the dominant flavor of neutrinos produced by accelerator based beams are of the muon type. This work reports new measurements of electron neutrino and antineutrino interactions in hydrocarbon, obtained by strongly suppressing backgrounds initiated by muon flavor neutrinos and antineutrinos. Double differential cross sections as a function of visible energy transfer, $E_\text{avail}$, and transverse momentum transfer, $p_T$, or three momentum transfer, $q_3$ are presented. △ Less

Submitted 16 April, 2024; v1 submitted 27 December, 2023; originally announced December 2023.

Comments: 25 pages, 32 figures and 7 tables, accepted for publication in Physical Review D. Revised to add content updated in review process

Report number: FERMILAB-PUB-23-0830-PPD

Journal ref: Physical Review D109, 092008 (2024)

Measuring the CMB primordial B-modes with Bolometric Interferometry

Authors: A. Mennella, P. Ade, A. Almela, G. Amico, L. H. Arnaldi, J. Aumont, S. Banfi, E. S. Battistelli, B. Bélier, L. Bergé, J. -Ph. Bernard, P. de Bernardis, M. Bersanelli, J. Bonaparte, J. D. Bonilla, E. Bunn, D. Buzi, F. Cacciotti, D. Camilieri, F. Cavaliere, P. Chanial, C. Chapron, L. Colombo, F. Columbro, A. Coppolecchia , et al. (89 additional authors not shown)

Abstract: The Q&U Bolometric Interferometer for Cosmology (QUBIC) is the first bolometric interferometer designed to measure the primordial B-mode polarization of the Cosmic Microwave Background (CMB). Bolometric interferometry is a novel technique that combines the sensitivity of bolometric detectors with the control of systematic effects that is typical of interferometry, both key features in the quest fo… ▽ More The Q&U Bolometric Interferometer for Cosmology (QUBIC) is the first bolometric interferometer designed to measure the primordial B-mode polarization of the Cosmic Microwave Background (CMB). Bolometric interferometry is a novel technique that combines the sensitivity of bolometric detectors with the control of systematic effects that is typical of interferometry, both key features in the quest for the faint signal of the primordial B-modes. A unique feature is the so-called "spectral imaging", i.e., the ability to recover the sky signal in several sub-bands within the physical band during data analysis. This feature provides an in-band spectral resolution of Δν/ν \sim 0.04 that is unattainable by a traditional imager. This is a key tool for controlling the Galactic foregrounds contamination. In this paper, we describe the principles of bolometric interferometry, the current status of the QUBIC experiment and future prospects. △ Less

Submitted 5 November, 2023; originally announced November 2023.

Comments: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Measurement of the Multi-Neutron $\barν_μ$ Charged Current Differential Cross Section at Low Available Energy on Hydrocarbon

Authors: A. Olivier, T. Cai, S. Akhter, Z. Ahmad Dar, V. Ansari, M. V. Ascencio, M. Sajjad Athar, A. Bashyal, A. Bercellie, M. Betancourt, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, G. A. Díaz, J. Felix, L. Fields, A. Filkins, R. Fine, A. M. Gago, P. K. Gaur, S. M. Gilligan, R. Gran, E. Granados, D. A. Harris , et al. (36 additional authors not shown)

Abstract: Neutron production in antineutrino interactions can lead to bias in energy reconstruction in neutrino oscillation experiments, but these interactions have rarely been studied. MINERvA previously studied neutron production at an average antineutrino energy of ~3 GeV in 2016 and found deficiencies in leading models. In this paper, the MINERvA 6 GeV average antineutrino energy data set is shown to ha… ▽ More Neutron production in antineutrino interactions can lead to bias in energy reconstruction in neutrino oscillation experiments, but these interactions have rarely been studied. MINERvA previously studied neutron production at an average antineutrino energy of ~3 GeV in 2016 and found deficiencies in leading models. In this paper, the MINERvA 6 GeV average antineutrino energy data set is shown to have similar disagreements. A measurement of the cross section for an antineutrino to produce two or more neutrons and have low visible energy is presented as an experiment-independent way to explore neutron production modeling. This cross section disagrees with several leading models' predictions. Neutron modeling techniques from nuclear physics are used to quantify neutron detection uncertainties on this result. △ Less

Submitted 21 November, 2023; v1 submitted 25 October, 2023; originally announced October 2023.

Comments: 25 pages, 11 figures; Added ancillary files with cross section values as .csv Matches preprint accepted by publisher

Report number: PUB-23-610-ND

Limits on ALP-neutrino couplings from loop-level processes

Authors: J. Bonilla, B. Gavela, J. Machado-Rodríguez

Abstract: We obtain new bounds on effective ALP-neutrino interactions from their loop-level impact on the ALP couplings to electroweak gauge bosons $γZ$, $ZZ$ and $W^+W^-$. Both types of interactions are connected via the chiral anomaly, which allows to overcome the fermion mass suppression of tree-level analyses. ALP couplings to charged leptons are considered simultaneously with ALP-neutrino ones, as requ… ▽ More We obtain new bounds on effective ALP-neutrino interactions from their loop-level impact on the ALP couplings to electroweak gauge bosons $γZ$, $ZZ$ and $W^+W^-$. Both types of interactions are connected via the chiral anomaly, which allows to overcome the fermion mass suppression of tree-level analyses. ALP couplings to charged leptons are considered simultaneously with ALP-neutrino ones, as required by gauge invariance. Complete one-loop computations are presented together with relevant kinematic limits. We leverage data obtained from rare meson decay measurements alongside observations from collider experiments, among others. The novel constraints are particularly competitive for ALP masses $\gtrsim 100$ MeV, evidencing the powerful potential of future ALP searches at colliders which target its couplings to electroweak gauge bosons. △ Less

Submitted 19 October, 2023; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: 24 pages, 6 figures. References and bounds have been adjusted. Acknowledgments have been extended

Report number: IFT-UAM/CSIC-23-121

Journal ref: Phys. Rev. D 109, 055023, 2024

High Energy Physics in Africa, Latin America and other developing regions

Authors: Kétévi A. Assamagan, Johan Sebastian Bonilla, Claudio Dib, Azwinndini Muronga, Heath B. O'Connell, Rogerio Rosenfeld, Suyog Shrestha

Abstract: We summarize the current status of high energy physics (HEP) in Africa, Latin America, and other developing regions We summarize the current status of high energy physics (HEP) in Africa, Latin America, and other developing regions △ Less

Submitted 17 October, 2023; v1 submitted 29 August, 2023; originally announced August 2023.

Comments: 6 pages, 2 figures

Exploring numerical blow-up phenomena for the Keller-Segel-Navier-Stokes equations

Authors: Jesús Bonilla, Juan Vicente Gutiérrez-Santacreu

Abstract: The Keller-Segel-Navier-Stokes system governs chemotaxis in liquid environments. This system is to be solved for the organism and chemoattractant densities and for the fluid velocity and pressure. It is known that if the total initial cell density mass is below $2π$ there exist globally defined generalised solutions, but what is less understood is whether there are blow-up solutions beyond such a… ▽ More The Keller-Segel-Navier-Stokes system governs chemotaxis in liquid environments. This system is to be solved for the organism and chemoattractant densities and for the fluid velocity and pressure. It is known that if the total initial cell density mass is below $2π$ there exist globally defined generalised solutions, but what is less understood is whether there are blow-up solutions beyond such a threshold and its optimality. Motivated by this issue, a numerical blow-up scenario is investigated. Approximate solutions computed via a stabilised finite element method founded on a shock capturing technique are such that they satisfy \emph{a priori} bounds as well as lower and $L^1(Ω)$ bounds for the cell and chemoattractant densities. In particular, this latter properties are essential in detecting numerical blow-up configurations, since the non-satisfaction of these two requirements might trigger numerical oscillations leading to non-realistic finite-time collapses into persistent Dirac-type measures. Our findings show that the existence threshold value $2π$ encountered for the cell density mass may not be optimal and hence it is conjectured that the critical threshold value $4π$ may be inherited from the fluid-free Keller-Segel equations. Additionally it is observed that the formation of singular points can be neglected if the fluid flow is intensified. △ Less

Submitted 31 January, 2023; originally announced February 2023.

Comments: 30 pages, 23 figures

MSC Class: 35Q35; 65N30; 92C17

Simultaneous measurement of muon neutrino quasielastic-like cross sections on CH, C, water, Fe, and Pb as a function of muon kinematics at MINERvA

Authors: J. Kleykamp, S. Akhter, Z. Ahmad Dar, V. Ansari, M. V. Ascencio, M. Sajjad Athar, A. Bashyal, A. Bercellie, M. Betancourt, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, T. Cai, M. F. Carneiro, G. A. Díaz, H. da Motta, S. A. Dytman, J. Felix, L. Fields, A. Filkins, R. Fine, A. M. Gago, H. Gallagher , et al. (43 additional authors not shown)

Abstract: This paper presents the first simultaneous measurement of the quasielastic-like neutrino-nucleus cross sections on C, water, Fe, Pb and scintillator (hydrocarbon or CH) as a function of longitudinal and transverse muon momentum. The ratio of cross sections per nucleon between Pb and CH is always above unity and has a characteristic shape as a function of transverse muon momentum that evolves slowl… ▽ More This paper presents the first simultaneous measurement of the quasielastic-like neutrino-nucleus cross sections on C, water, Fe, Pb and scintillator (hydrocarbon or CH) as a function of longitudinal and transverse muon momentum. The ratio of cross sections per nucleon between Pb and CH is always above unity and has a characteristic shape as a function of transverse muon momentum that evolves slowly as a function of longitudinal muon momentum. The ratio is constant versus longitudinal momentum within uncertainties above a longitudinal momentum of 4.5GeV/c. The cross section ratios to CH for C, water, and Fe remain roughly constant with increasing longitudinal momentum, and the ratios between water or C to CH do not have any significant deviation from unity. Both the overall cross section level and the shape for Pb and Fe as a function of transverse muon momentum are not reproduced by current neutrino event generators. These measurements provide a direct test of nuclear effects in quasielastic-like interactions, which are major contributors to long-baseline neutrino oscillation data samples. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 9 pages, 8 flgures, including supplemental material

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

Muon capture on the deuteron in chiral effective field theory

Authors: Jose Bonilla, Bijaya Acharya, Lucas Platter

Abstract: We consider the capture of a muon on a deuteron. An uncertainty analysis of the dominant channels is important for a careful analysis of forthcoming experimental data. We quantify the theoretical uncertainties of chiral effective-field-theory predictions of the muon-deuteron capture rate from the relevant neutron-neutron partial wave channels in the final state. We study the dependence on the cuto… ▽ More We consider the capture of a muon on a deuteron. An uncertainty analysis of the dominant channels is important for a careful analysis of forthcoming experimental data. We quantify the theoretical uncertainties of chiral effective-field-theory predictions of the muon-deuteron capture rate from the relevant neutron-neutron partial wave channels in the final state. We study the dependence on the cutoff used to regularize the interactions, low-energy constants calibrated using different fitting data and strategies, and truncation of the effective-field-theory expansion of the currents. Combining these approaches gives as an estimate of $Γ^{1/2}_{μd} = 399.1 \pm 7.6 \pm 4.4$ s$^{-1}$ for capture from the atomic doublet state, and $Γ^{3/2}_{μd} = 12.31 \pm 0.47 \pm 0.04$ s$^{-1}$ for capture from the quartet state. △ Less

Submitted 15 December, 2022; originally announced December 2022.

Comments: 17 pages, 3 figures

Community Engagement Frontier

Authors: Ketevi A. Assamagan, Breese Quinn, Kenneth Bloom, Veronique Boisvert, Carla Bonifazi, Johan S. Bonilla, Mu-Chun Chen, Sarah M. Demers, Farah Fahim, Rob Fine, Mike Headley, Julie Hogan, Kathryn Jepsen, Sijbrand de Jong, Aneliya Karadzhinova-Ferrer, Yi-Hsuan Lin, Don Lincoln, Sudhir Malik, Alex Murokh, Azwinndini Muronga, Randal Ruchti, Louise Suter, Koji Yoshimura

Abstract: This is the summary report of the Community Engagement Frontier for the Snowmass 2021 study of the future of particle physics. The report discusses a number of general issues of importance to the particle physics community, including (1) the relation of universities, national laboratories, and industry, (2) career paths for scientists engaged in particle physics, (3) diversity, equity, and inclusi… ▽ More This is the summary report of the Community Engagement Frontier for the Snowmass 2021 study of the future of particle physics. The report discusses a number of general issues of importance to the particle physics community, including (1) the relation of universities, national laboratories, and industry, (2) career paths for scientists engaged in particle physics, (3) diversity, equity, and inclusion, (4) physics education, (5) public education and outreach, (6) engagement with the government and public policy, and (7) the environmental and social impacts of particle physics. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 51 pages, 2 figures; to appear in the proceedings of the 2021 Community Study on the Future of U.S. Particle Physics (Snowmass 2021)

High-Statistics Measurement of Antineutrino Quasielastic-like scattering at $E_ν\sim$ 6~GeV on a Hydrocarbon Target

Authors: A. Bashyal, S. Akhter, Z. Ahmad Dar, F. Akbar, V. Ansari, M. V. Ascencio, M. Sajjad Athar, A. Bercellie, M. Betancourt, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, M. F. Carneiro, G. A. Díaz, J. Felix, L. Fields, A. Filkins, R. Fine, A. M. Gago, H. Gallagher, P. K. Gaur, S. M. Gilligan, R. Gran , et al. (44 additional authors not shown)

Abstract: We present measurements of the cross section for anti-neutrino charged-current quasielastic-like scattering on hydrocarbon using the medium energy (ME) NuMI wide-band neutrino beam peaking at $<E_ν>\sim 6$ GeV. The cross section measurements are presented as a function of the longitudinal momentum ($p_{||}$) and transverse momentum ($p_{T}$) of the final state muon. This work complements our previ… ▽ More We present measurements of the cross section for anti-neutrino charged-current quasielastic-like scattering on hydrocarbon using the medium energy (ME) NuMI wide-band neutrino beam peaking at $<E_ν>\sim 6$ GeV. The cross section measurements are presented as a function of the longitudinal momentum ($p_{||}$) and transverse momentum ($p_{T}$) of the final state muon. This work complements our previously reported high statistics measurement in the neutrino channel and extends the previous anti-neutrino measurement made in the low energy (LE) beam at neutrino energy($<E_ν>$) $\sim$ 3.5 GeV to $p_{T}$ of 2.5 GeV/c. Current theoretical models do not completely describe the data in this previously unexplored high $p_{T}$ region. The single differential cross section as a function of four momentum transfer ($Q^{2}_{QE}$) now extends to 4 GeV$^2$ with high statistics. The cross section as a function of $Q^{2}_{QE}$ shows that the tuned simulations developed by the MINERvA collaboration that agreed well with the low energy beam measurements do not agree as well with the medium energy beam measurements. Newer neutrino interaction models such as the GENIE 3 tunes are better able to simulate the high $Q^{2}_{QE}$. △ Less

Submitted 25 June, 2023; v1 submitted 18 November, 2022; originally announced November 2022.

Comments: 17 pages, 14 figures

Report number: FERMILAB-PUB-22-835

Neutral $B$-anomalies from an $\mathit{on\text{-}shell}$ scalar exchange

Authors: Jesus Bonilla, Arturo de Giorgi, Maria Ramos

Abstract: The neutral $B$-anomalies are analysed in terms of the tree-level exchange of a (pseudo)scalar gauge singlet $a$. Solutions to both $R_{K^{(*)}}$ central bin anomalies are found within $1σ$ for ${m_{a}^2 \in [{1.1},\,{6}]\,\text{GeV}^2}$, while the low $q^2$-bin anomaly can also be accounted for with masses close to the bin threshold. The impact of these solutions on other ${b\to s e^+ e^-}$ obser… ▽ More The neutral $B$-anomalies are analysed in terms of the tree-level exchange of a (pseudo)scalar gauge singlet $a$. Solutions to both $R_{K^{(*)}}$ central bin anomalies are found within $1σ$ for ${m_{a}^2 \in [{1.1},\,{6}]\,\text{GeV}^2}$, while the low $q^2$-bin anomaly can also be accounted for with masses close to the bin threshold. The impact of these solutions on other ${b\to s e^+ e^-}$ observables is discussed in detail. Due to the $\mathit{on\text{-}shell}$ enhancement, sizable effects are expected in null tests of the SM, such as the flat term, $F_H$, of the $B\to K e^+ e^-$ angular distribution. At the same time, the observable sensitive to the $K^\ast$ polarisation, $F_L\, (B\to K^* e^+ e^-)$, and the lepton forward-backward asymmetry, ${A_{FB}\, (B\to K^* e^+ e^-)}$, can be suppressed with respect to their SM values. Corrections from the new physics to $\mathcal{B}(B_s \to e^+ e^-)$ are, on the other hand, negligible. Along with the previous observables, improved measurements of the cross section ${σ(e^+ e^- \to a (e^+ e^-) γ)}$ could potentially probe the relevant parameter space of the model. A comparison between our results and those stemming from an axion-like particle exchange is also discussed, showing that the exchange of a general scalar singlet offers a noticeably wider parameter space. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 8 pages, 5 figures

Status of QUBIC, the Q&U Bolometer for Cosmology

Authors: L. Mousset, P. Ade, A. Almela, G. Amico, L. H. Arnaldi, J. Aumont, S. Banfi, E. S. Battistelli, B. Bélier, L. Bergé, J. -Ph. Bernard, P. de Bernardis, M. Bersanelli, J. Bonaparte, J. D. Bonilla, E. Bunn, D. Buzi, D. Camilieri, F. Cavaliere, P. Chanial, C. Chapron, S. Colombo, F. Columbro, A. Coppolecchia, B. Costanza , et al. (86 additional authors not shown)

Abstract: The Q&U Bolometric Interferometer for Cosmology (QUBIC) is a novel kind of polarimeter optimized for the measurement of the B-mode polarization of the Cosmic Microwave Back-ground (CMB), which is one of the major challenges of observational cosmology. The signal is expected to be of the order of a few tens of nK, prone to instrumental systematic effects and polluted by various astrophysical foregr… ▽ More The Q&U Bolometric Interferometer for Cosmology (QUBIC) is a novel kind of polarimeter optimized for the measurement of the B-mode polarization of the Cosmic Microwave Back-ground (CMB), which is one of the major challenges of observational cosmology. The signal is expected to be of the order of a few tens of nK, prone to instrumental systematic effects and polluted by various astrophysical foregrounds which can only be controlled through multichroic observations. QUBIC is designed to address these observational issues with a novel approach that combines the advantages of interferometry in terms of control of instrumental systematics with those of bolometric detectors in terms of wide-band, background-limited sensitivity. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: Contribution to the 2022 Cosmology session of the 33rd Rencontres de Blois. arXiv admin note: substantial text overlap with arXiv:2203.08947

Neutrino-induced coherent $π^{+}$ production in C, CH, Fe and Pb at $\langle E_ν\rangle \sim 6$ GeV

Authors: M. A. Ramírez, S. Akhter, Z. Ahmad Dar, F. Akbar, V. Ansari, M. V. Ascencio, M. Sajjad Athar, A. Bashyal, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, T. Cai, G. A. Díaz, H. da Motta, S. A. Dytman, J. Felix, L. Fields, A. Filkins, R. Fine, H. Gallagher , et al. (41 additional authors not shown)

Abstract: MINERvA has measured the $ν_μ$-induced coherent $π^{+}$ cross section simultaneously in hydrocarbon (CH), graphite (C), iron (Fe) and lead (Pb) targets using neutrinos from 2 to 20 GeV. The measurements exceed the predictions of the Rein-Sehgal and Berger-Sehgal PCAC based models at multi-GeV $ν_μ$ energies and at produced $π^{+}$ energies and angles, $E_π>1$ GeV and $θ_π<10^{\circ}$. Measurements… ▽ More MINERvA has measured the $ν_μ$-induced coherent $π^{+}$ cross section simultaneously in hydrocarbon (CH), graphite (C), iron (Fe) and lead (Pb) targets using neutrinos from 2 to 20 GeV. The measurements exceed the predictions of the Rein-Sehgal and Berger-Sehgal PCAC based models at multi-GeV $ν_μ$ energies and at produced $π^{+}$ energies and angles, $E_π>1$ GeV and $θ_π<10^{\circ}$. Measurements of the cross-section ratios of Fe and Pb relative to CH reveal the effective $A$-scaling to increase from an approximate $A^{1/3}$ scaling at few GeV to an $A^{2/3}$ scaling for $E_ν>10$ GeV. △ Less

Submitted 26 June, 2023; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: 8 pages, 6 figures and 60 pages of supplemental material. Updated to accepted version, with ancillary files of the cross section data

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

Diversity, Equity, and Inclusion in Particle Physics

Authors: C. Bonifazi, J. S. Bonilla, M. -C. Chen, Y. H. Lin, K. A. Assamagan, E. V. Hansen, S. Meehan, E. Smith

Abstract: To achieve the highest level of intellectual excellence calls for the greatest extent of diversity. However, due to the unjust institutional and societal barriers, the field of particle physics remains as one of the least diverse fields, severely limiting the potential of our scientific achievements. In order for the US Particle Physics Community, including the accelerator science and engineering… ▽ More To achieve the highest level of intellectual excellence calls for the greatest extent of diversity. However, due to the unjust institutional and societal barriers, the field of particle physics remains as one of the least diverse fields, severely limiting the potential of our scientific achievements. In order for the US Particle Physics Community, including the accelerator science and engineering fields, to remain at the forefront of global scientific leadership, it is imperative for our community to act urgently and diligently to improve the status quo of diversity, equity, inclusion, and accessibility (DEIA). This report summarizes the community contributions and recommendations within the Diversity and Inclusion Topical Group of the Community Engagement Frontier (CEF3) at the Snowmass 2021. △ Less

Submitted 25 September, 2022; originally announced September 2022.

Comments: Snowmass 2021 CEF3 TG Report

The cost of an ALP solution to the neutral $B$-anomalies

Authors: J. Bonilla, A. de Giorgi, B. Gavela, L. Merlo, M. Ramos

Abstract: The neutral anomalies in $B$ decays are analysed in terms of the tree-level exchange of an axion-like-particle (ALP), within the effective field theory framework. The complete two-dimensional parameter space for ALP couplings to electrons and muons is explored. The solutions to $R_K$ and to the two energy bins of $R_{K^\ast}$ are confronted with the impact of ALP exchange on other observables (mes… ▽ More The neutral anomalies in $B$ decays are analysed in terms of the tree-level exchange of an axion-like-particle (ALP), within the effective field theory framework. The complete two-dimensional parameter space for ALP couplings to electrons and muons is explored. The solutions to $R_K$ and to the two energy bins of $R_{K^\ast}$ are confronted with the impact of ALP exchange on other observables (meson oscillations, leptonic and semileptonic decays of $B$ mesons including searches for new resonances, astrophysical constraints), as well as with the theoretical domain of validity of the effective theory. Solutions based on ALPs heavier than $B$ mesons, or lighter than twice the muon mass, are shown to be excluded. In contrast, the exchange of on-shell ALPs provides solutions to $R_K$ and/or $R_{K^\ast}$ within $2σ$ sensitivity which are technically compatible with those constraints. Furthermore, a ''golden ALP mass'' is identified at the frontier between the two energy bin windows of $R_{K^\ast}$, which could simultaneously explain these two $R_{K^\ast}$ anomalies together with $R_K$; this calls for the convenience of different energy binning which would easily clear up this (unlikely) possibility. The impact of smearing on data analysis is also discussed. When loop effects are taken into account, the solutions found can be in addition compatible with the data on the $g-2$ of the electron but not simultaneously with those on the $g-2$ of the muon. Furthermore, loop effects may require fine-tunings of some coupling values. △ Less

Submitted 17 February, 2023; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: 50 pages, 15 figures

Simultaneous measurement of muon neutrino $ν_μ$ charged-current single $π^+$ production in CH, C, H$_2$O, Fe, and Pb targets in MINERvA

Authors: A. Bercellie, K. A. Kroma-Wiley, S. Akhter, Z. Ahmad Dar, F. Akbar, V. Ansari, M. V. Ascencio, M. Sajjad Athar, L. Bellantoni, M. Betancourt, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, T. Cai, G. A. Díaz, H. da Motta, S. A. Dytman, J. Felix, L. Fields, A. Filkins, R. Fine, A. M. Gago, H. Gallagher , et al. (47 additional authors not shown)

Abstract: Neutrino-induced charged-current single $π^+$ production in the $Δ(1232)$ resonance region is of considerable interest to accelerator-based neutrino oscillation experiments. In this work, high statistics differential cross sections are reported for the semi-exclusive reaction $ν_μA \to μ^- π^+ +$ nucleon(s) on scintillator, carbon, water, iron, and lead targets recorded by MINERvA using a wide-ban… ▽ More Neutrino-induced charged-current single $π^+$ production in the $Δ(1232)$ resonance region is of considerable interest to accelerator-based neutrino oscillation experiments. In this work, high statistics differential cross sections are reported for the semi-exclusive reaction $ν_μA \to μ^- π^+ +$ nucleon(s) on scintillator, carbon, water, iron, and lead targets recorded by MINERvA using a wide-band $ν_μ$ beam with $\left< E_ν\right> \approx 6$~GeV. Suppression of the cross section at low $Q^2$ and enhancement of low $T_π$ are observed in both light and heavy nuclear targets compared to phenomenological models used in current neutrino interaction generators. The cross-section ratios for iron and lead compared to CH across the kinematic variables probed are 0.8 and 0.5 respectively, a scaling which is also not predicted by current generators. △ Less

Submitted 12 July, 2023; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 6 pages, 6 figures, 117 pages of supplementary material; updated to accepted version and added ancillary files with cross section data; updated again to add ancillary files with neutrino flux used for result

Report number: FERMILAB-PUB-22-716-ND-QIS-SCD

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

Improved constraint on the MINERvA medium energy neutrino flux using $\barνe^{-} \!\rightarrow \barνe^{-}$ data

Authors: L. Zazueta, S. Akhter, Z. Ahmad Dar, F. Akbar, V. Ansari, M. V. Ascencio, M. Sajjad Athar, A. Bashyal, A. Bercellie, M. Betancourt, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, T. Cai, G. A. Díaz, H. da Motta, J. Felix, L. Fields, A. Filkins, R. Fine, A. M. Gago, H. Gallagher, A. Ghosh, S. M. Gilligan , et al. (36 additional authors not shown)

Abstract: Processes with precisely known cross sections, like neutrino electron elastic scattering ($νe^{-} \!\rightarrow νe^{-}$) and inverse muon decay ($ν_μe^{-} \!\rightarrow μ^{-} ν_e$) have been used by MINERvA to constrain the uncertainty on the NuMI neutrino beam flux. This work presents a new measurement of neutrino elastic scattering with electrons using the medium energy \numubar enhanced NuMI be… ▽ More Processes with precisely known cross sections, like neutrino electron elastic scattering ($νe^{-} \!\rightarrow νe^{-}$) and inverse muon decay ($ν_μe^{-} \!\rightarrow μ^{-} ν_e$) have been used by MINERvA to constrain the uncertainty on the NuMI neutrino beam flux. This work presents a new measurement of neutrino elastic scattering with electrons using the medium energy \numubar enhanced NuMI beam. A sample of 578 events after background subtraction is used in combination with the previous measurement on the \numu beam and the inverse muon decay measurement to reduce the uncertainty on the \numu flux in the \numu-enhanced beam from 7.6\% to 3.3\% and the \numubar flux in the \numubar-enhanced beam from 7.8\% to 4.7\%. △ Less

Submitted 12 September, 2022; originally announced September 2022.

Comments: 12 pages, 16 figures

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

Bound-preserving finite element approximations of the Keller-Segel equations

Authors: Santiago Badia, Jesús Bonilla, Juan Vicente Gutiérrez-Santacreu

Abstract: This paper aims to develop numerical approximations of the Keller--Segel equations that mimic at the discrete level the lower bounds and the energy law of the continuous problem. We solve these equations for two unknowns: the organism (or cell) density, which is a positive variable, and the chemoattractant density, which is a nonnegative variable. We propose two algorithms, which combine a stabili… ▽ More This paper aims to develop numerical approximations of the Keller--Segel equations that mimic at the discrete level the lower bounds and the energy law of the continuous problem. We solve these equations for two unknowns: the organism (or cell) density, which is a positive variable, and the chemoattractant density, which is a nonnegative variable. We propose two algorithms, which combine a stabilized finite element method and a semi-implicit time integration. The stabilization consists of a nonlinear artificial diffusion that employs a graph-Laplacian operator and a shock detector that localizes local extrema. As a result, both algorithms turn out to be nonlinear.Both algorithms can generate cell and chemoattractant numerical densities fulfilling lower bounds. However, the first algorithm requires a suitable constraint between the space and time discrete parameters, whereas the second one does not. We design the latter to attain a discrete energy law on acute meshes. We report some numerical experiments to validate the theoretical results on blowup and non-blowup phenomena. In the blowup setting, we identify a \textit{locking} phenomenon that relates the $L^\infty(Ω)$-norm to the $L^1(Ω)$-norm limiting the growth of the singularity when supported on a macroelement. △ Less

Submitted 22 July, 2022; originally announced July 2022.

Comments: 27 pages, 22 figures

MSC Class: 35K20; 35K55; 65N30

ALP one-loop corrections and nonresonant searches

Authors: Jesús Bonilla

Abstract: We consider the complete leading-order -- dimension five -- effective linear Lagrangian for ALPs. The full set of one-loop contributions to ALP-SM effective couplings are derived, including all finite corrections. A few phenomenological consequences of these computations are also explored as illustration, with flavour diagonal channels in the case of fermions: in particular, we explore one-loop co… ▽ More We consider the complete leading-order -- dimension five -- effective linear Lagrangian for ALPs. The full set of one-loop contributions to ALP-SM effective couplings are derived, including all finite corrections. A few phenomenological consequences of these computations are also explored as illustration, with flavour diagonal channels in the case of fermions: in particular, we explore one-loop constraints on the coupling of the ALP to top quarks. Additionally, we propose a new search for ALPs, targeting VBS processes at the LHC. For this, we consider the tree-level diboson production in VBS, where the ALP participates as an off-shell mediator. Upper limits on ALP EW couplings are obtained from a reinterpretation of Run 2 public CMS VBS analyses. Simple projections for LHC Run 3 and HL-LHC are also calculated, demonstrating the power of future dedicated analyses at ATLAS and CMS. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: 4 pages, 4 figures. Contribution to the 2022 EW session of the 56th Rencontres de Moriond. arXiv admin note: substantial text overlap with arXiv:2111.14750

Accessibility in High Energy Physics: Lessons from the Snowmass Process

Authors: K. A. Assamagan, C. Bonifazi, J. S. Bonilla, P. A. Breur, M. -C. Chen, T. Y. Chen, A. Roepe-Gier, Y. H. Lin, S. Meehan, M. E. Monzani, E. Novitski, G. Stark

Abstract: Accessibility to participation in the high energy physics community can be impeded by many barriers. These barriers must be acknowledged and addressed to make access more equitable in the future. An accessibility survey, the Snowmass Summer Study attendance survey, and an improved accessibility survey were sent to the Snowmass2021 community. This paper will summarize and present the barriers that… ▽ More Accessibility to participation in the high energy physics community can be impeded by many barriers. These barriers must be acknowledged and addressed to make access more equitable in the future. An accessibility survey, the Snowmass Summer Study attendance survey, and an improved accessibility survey were sent to the Snowmass2021 community. This paper will summarize and present the barriers that prevent people from participating in the Snowmass2021 process,recommendations for the various barriers, and discussions of resources and funding needed to enact these recommendations, based on the results of all three surveys, along with personal experiences of the community members. △ Less

Submitted 1 March, 2023; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021

Simultaneous measurement of proton and lepton kinematics in quasielastic-like $ν_μ$-hydrocarbon interactions from 2 to 20 GeV

Authors: The MINERvA Collaboration, D. Ruterbories, S. Akhter, Z. Ahmad Dar, F. Akbar, V. Ansari, M. V. Ascencio, M. Sajjad Athar, A. Bashyal, A. Bercellie, M. Betancourt, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, T. Cai, M. F. Carneiro, G. A. Díaz, H. da Motta, J. Felix, L. Fields, A. Filkins, R. Fine, A. M. Gago , et al. (49 additional authors not shown)

Abstract: Neutrino charged-current quasielastic-like scattering, a reaction category extensively used in neutrino oscillation measurements, probes nuclear effects that govern neutrino-nucleus interactions. This Letter reports the first measurement of the triple-differential cross section for $ν_μ$ quasielastic-like reactions using the hydrocarbon medium of the MINERvA detector exposed to a wide-band beam sp… ▽ More Neutrino charged-current quasielastic-like scattering, a reaction category extensively used in neutrino oscillation measurements, probes nuclear effects that govern neutrino-nucleus interactions. This Letter reports the first measurement of the triple-differential cross section for $ν_μ$ quasielastic-like reactions using the hydrocarbon medium of the MINERvA detector exposed to a wide-band beam spanning 2 $\leq$ E$_ν\leq$ 20 GeV. The measurement maps the correlations among transverse and longitudinal muon momenta and summed proton kinetic energies, and compares them to predictions from a state-of-art simulation. Discrepancies are observed that likely reflect shortfalls with modeling of pion and nucleon intranuclear scattering and/or spectator nucleon ejection from struck nuclei. The separate determination of leptonic and hadronic variables can inform experimental approaches to neutrino-energy estimation. △ Less

Submitted 25 May, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: 6 pages, 4 figures. Supplementary 18 figures

Report number: FERMILAB-PUB-22-130-PPD-QIS-SCD

Jets and Jet Substructure at Future Colliders

Authors: Ben Nachman, Salvatore Rappoccio, Nhan Tran, Johan Bonilla, Grigorios Chachamis, Barry M. Dillon, Sergei V. Chekanov, Robin Erbacher, Loukas Gouskos, Andreas Hinzmann, Stefan Höche, B. Todd Huffman, Ashutosh. V. Kotwal, Deepak Kar, Roman Kogler, Clemens Lange, Matt LeBlanc, Roy Lemmon, Christine McLean, Mark S. Neubauer, Tilman Plehn, Debarati Roy, Giordan Stark, Jennifer Roloff, Marcel Vos , et al. (2 additional authors not shown)

Abstract: Even though jet substructure was not an original design consideration for the Large Hadron Collider (LHC) experiments, it has emerged as an essential tool for the current physics program. We examine the role of jet substructure on the motivation for and design of future energy frontier colliders. In particular, we discuss the need for a vibrant theory and experimental research and development prog… ▽ More Even though jet substructure was not an original design consideration for the Large Hadron Collider (LHC) experiments, it has emerged as an essential tool for the current physics program. We examine the role of jet substructure on the motivation for and design of future energy frontier colliders. In particular, we discuss the need for a vibrant theory and experimental research and development program to extend jet substructure physics into the new regimes probed by future colliders. Jet substructure has organically evolved with a close connection between theorists and experimentalists and has catalyzed exciting innovations in both communities. We expect such developments will play an important role in the future energy frontier physics program. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Nonresonant Searches for Axion-Like Particles in Vector Boson Scattering Processes at the LHC

Authors: J. Bonilla, I. Brivio, J. Machado-Rodríguez, J. F. de Trocóniz

Abstract: We propose a new search for Axion-Like Particles (ALPs), targeting Vector Boson Scattering (VBS) processes at the LHC. We consider nonresonant ALP-mediated VBS, where the ALP participates as an off-shell mediator. This process occurs whenever the ALP is too light to be produced resonantly, and it takes advantage of the derivative nature of ALP interactions with the electroweak Standard Model boson… ▽ More We propose a new search for Axion-Like Particles (ALPs), targeting Vector Boson Scattering (VBS) processes at the LHC. We consider nonresonant ALP-mediated VBS, where the ALP participates as an off-shell mediator. This process occurs whenever the ALP is too light to be produced resonantly, and it takes advantage of the derivative nature of ALP interactions with the electroweak Standard Model bosons. We study the production of $ZZ$, $Zγ$, $W^\pm γ$, $W^\pm Z$ and $W^\pm W^\pm$ pairs with large diboson invariant masses in association with two jets. Working in a gauge-invariant framework, upper limits on ALP couplings to electroweak bosons are obtained from a reinterpretation of Run 2 public CMS VBS analyses. The constraints inferred on ALP couplings to $ZZ$, $Zγ$ and $W^\pm W^\pm$ pairs are very competitive for ALP masses up to 100 GeV. They have the advantage of being independent of the ALP coupling to gluons and of the ALP decay width. Simple projections for LHC Run 3 and HL-LHC are also calculated, demonstrating the power of future dedicated analyses at ATLAS and CMS. △ Less

Submitted 18 May, 2022; v1 submitted 7 February, 2022; originally announced February 2022.

Comments: 27 pages, 11 figures, 6 tables. Version accepted for publication in JHEP

Report number: IFT-UAM/CSIC-22-7, VBSCAN-PUB-01-22

Journal ref: JHEP 06 (2022) 113

Vertex finding in neutrino-nucleus interaction: A Model Architecture Comparison

Authors: F. Akbar, A. Ghosh, S. Young, S. Akhter, Z. Ahmad Dar, V. Ansari, M. V. Ascencio, M. Sajjad Athar, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, T. Cai, M. F. Carneiro, G. A. Díaz, J. Felix, L. Fields, A. Filkins, R. Fine, P. K. Gaura, R. Gran, D. A. Harris, D. Jena, S. Jena , et al. (26 additional authors not shown)

Abstract: We compare different neural network architectures for Machine Learning (ML) algorithms designed to identify the neutrino interaction vertex position in the MINERvA detector. The architectures developed and optimized by hand are compared with the architectures developed in an automated way using the package "Multi-node Evolutionary Neural Networks for Deep Learning" (MENNDL), developed at Oak Ridge… ▽ More We compare different neural network architectures for Machine Learning (ML) algorithms designed to identify the neutrino interaction vertex position in the MINERvA detector. The architectures developed and optimized by hand are compared with the architectures developed in an automated way using the package "Multi-node Evolutionary Neural Networks for Deep Learning" (MENNDL), developed at Oak Ridge National Laboratory (ORNL). The two architectures resulted in a similar performance which suggests that the systematics associated with the optimized network architecture are small. Furthermore, we find that while the domain expert hand-tuned network was the best performer, the differences were negligible and the auto-generated networks performed well. There is always a trade-off between human, and computer resources for network optimization and this work suggests that automated optimization, assuming resources are available, provides a compelling way to save significant expert time. △ Less

Submitted 7 January, 2022; originally announced January 2022.

One-loop corrections to ALP couplings

Authors: J. Bonilla, I. Brivio, M. B. Gavela, V. Sanz

Abstract: We derive the one-loop contributions to ALP-SM effective couplings, including all finite corrections. The complete leading-order -- dimension five -- effective linear Lagrangian is considered. The ALP is left off-shell, which is of particular impact on LHC and accelerator searches of ALP couplings to $γγ$, $ZZ$, $Zγ$, $WW$, gluons and fermions. All results are obtained in the covariant $R_ξ$ gauge… ▽ More We derive the one-loop contributions to ALP-SM effective couplings, including all finite corrections. The complete leading-order -- dimension five -- effective linear Lagrangian is considered. The ALP is left off-shell, which is of particular impact on LHC and accelerator searches of ALP couplings to $γγ$, $ZZ$, $Zγ$, $WW$, gluons and fermions. All results are obtained in the covariant $R_ξ$ gauge. A few phenomenological consequences are also explored as illustration, with flavour diagonal channels in the case of fermions: in particular, we explore constraints on the coupling of the ALP to top quarks, that can be extracted from LHC data, from astrophysical sources and from Dark Matter direct detection experiments such as PandaX, LUX and XENON1T. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 4 pages, 4 figures. Contribution to the proceedings of EPS-HEP 2021, European Physical Society conference on High Energy Physics, 26-30 July 2021

Measurement of inclusive charged-current $ν_{\numu}$ scattering on hydrocarbon at {<Enu>} 6 GeV with low three-momentum transfer

Authors: M. V. Ascencio, D. A. Andrade, I. Mahbub, Z. Ahmad Dar, F. Akbar, A. Bashyal, S. Bender, A. Bercellie, M. Betancourt, A. Bodek, J. L. Bonilla, K. Bonin, H. Budd, T. Cai, M. F. Carneiro, G. A. Diaz, H. da Motta, J. Felix, L. Fields, A. Filkins, R. Fine, N. Fuad, A. M. Gago, H. Gallagher, A. Ghosh , et al. (41 additional authors not shown)

Abstract: The \minerva experiment reports double-differential cross-section measurements for $ν_μ$-carbon interactions with three-momentum transfer $|\vec{q}| < 1.2$ GeV obtained with medium energy exposures in the NuMI beam. These measurements are performed as a function of the three-momentum transfer and an energy transfer estimator called the available energy defined as the energy that would be visible i… ▽ More The \minerva experiment reports double-differential cross-section measurements for $ν_μ$-carbon interactions with three-momentum transfer $|\vec{q}| < 1.2$ GeV obtained with medium energy exposures in the NuMI beam. These measurements are performed as a function of the three-momentum transfer and an energy transfer estimator called the available energy defined as the energy that would be visible in the detector. The double differential cross sections are compared to the GENIE and NuWro predictions along with the modified version of GENIE which incorporates new models for better agreement with earlier measurements from MINERvA. In these measurements, the quasi-elastic, resonance, and multi-nucleon knockout processes appear at different kinematics in this two-dimensional space. The results can be used to improve models for neutrino interactions needed by neutrino oscillation experiments. △ Less

Submitted 25 July, 2022; v1 submitted 25 October, 2021; originally announced October 2021.

One-loop corrections to ALP couplings

Authors: J. Bonilla, I. Brivio, M. B. Gavela, V. Sanz

Abstract: The plethora of increasingly precise experiments which hunt for axion-like particles (ALPs), as well as their widely different energy reach, call for the theoretical understanding of ALP couplings at loop-level. We derive the one-loop contributions to ALP-SM effective couplings, including finite corrections. The complete leading-order -- dimension five -- effective linear Lagrangian is considered.… ▽ More The plethora of increasingly precise experiments which hunt for axion-like particles (ALPs), as well as their widely different energy reach, call for the theoretical understanding of ALP couplings at loop-level. We derive the one-loop contributions to ALP-SM effective couplings, including finite corrections. The complete leading-order -- dimension five -- effective linear Lagrangian is considered. The ALP is left off-shell, which is of particular impact on LHC and accelerator searches of ALP couplings to $γγ$, $ZZ$, $Zγ$, $WW$, gluons and fermions. All results are obtained in the covariant $R_ξ$ gauge. A few phenomenological consequences are also explored as illustration, with flavour diagonal channels in the case of fermions: in particular, we explore constraints on the coupling of the ALP to top quarks, that can be extracted from LHC data, from astrophysical sources and from Dark Matter direct detection experiments such as PandaX, LUX and XENON1T. Furthermore, we clarify the relation between alternative ALP bases, the role of gauge anomalous couplings and their interface with chirality-conserving and chirality-flip fermion interactions, and we briefly discuss renormalization group aspects. △ Less

Submitted 29 October, 2021; v1 submitted 23 July, 2021; originally announced July 2021.

Comments: 58 pages, 11 figures. Full one-loop results available at https://notebookarchive.org/2021-07-9otlr9o . Algebraic error corrected, which does not affect the bulk of the computations

Report number: IFT-UAM/CSIC-21-82

Journal ref: JHEP11(2021)168

Exploring Neutrino-Nucleus Interactions in the GeV Regime using MINERvA

Authors: X. -G. Lu, Z. Ahmad Dar, F. Akbar, D. A. Andrade, M. V. Ascencio, G. D. Barr, A. Bashyal, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, J. L. Bonilla, H. Budd, G. Caceres, T. Cai, M. F. Carneiro, H. da Motta, G. A. Diaz, J. Felix, L. Fields, A. Filkins, R. Fine, A. M. Gago, H. Gallagher, S. M. Gilligan , et al. (42 additional authors not shown)

Abstract: With the advance of particle accelerator and detector technologies, the neutrino physics landscape is rapidly expanding. As neutrino oscillation experiments enter the intensity and precision frontiers, neutrino-nucleus interaction measurements are providing crucial input. MINERvA is an experiment at Fermilab dedicated to the study of neutrino-nucleus interactions in the regime of incident neutrino… ▽ More With the advance of particle accelerator and detector technologies, the neutrino physics landscape is rapidly expanding. As neutrino oscillation experiments enter the intensity and precision frontiers, neutrino-nucleus interaction measurements are providing crucial input. MINERvA is an experiment at Fermilab dedicated to the study of neutrino-nucleus interactions in the regime of incident neutrino energies from one to few GeV. The experiment recorded neutrino and antineutrino scattering data with the NuMI beamline from 2009 to 2019 using the Low-Energy and Medium-Energy beams that peak at 3 GeV and 6 GeV, respectively. This article reviews the broad spectrum of interesting nuclear and particle physics that MINERvA investigations have illuminated. The newfound, detailed knowledge of neutrino interactions with nuclear targets thereby obtained is proving essential to continued progress in the neutrino physics sector. △ Less

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

Comments: 22 pages, 17 figures. Accepted as invited contribution to EPJ special issue "Neutrino Interactions in the Intermediate and High Energy Region. Guest editors: Mohammad Sajjad Athar and Shri Krishna Singh"

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

Constraining the NuMI neutrino flux using inverse muon decay reactions in MINERvA

Authors: D. Ruterbories, Z. Ahmad Dar, F. Akbar, M. V. Ascencio, A. Bashyal, A. Bercellie, M. Betancourt, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, T. Cai, M. F. Carneiro, G. A. DÍaz, H. da Motta, J. Felix, L. Fields, A. Filkins, R. Fine, A. M. Gago, H. Gallagher, A. Ghosh, R. Gran, D. A. Harris , et al. (39 additional authors not shown)

Abstract: Inverse muon decay, $ν_μe^-\toμ^-ν_e$, is a reaction whose cross-section can be predicted with very small uncertainties. It has a neutrino energy threshold of $\approx 11$ GeV and can be used to constrain the high-energy part of the flux in the NuMI neutrino beam. This reaction is the dominant source of events which only contain high-energy muons nearly parallel to the direction of the neutrino be… ▽ More Inverse muon decay, $ν_μe^-\toμ^-ν_e$, is a reaction whose cross-section can be predicted with very small uncertainties. It has a neutrino energy threshold of $\approx 11$ GeV and can be used to constrain the high-energy part of the flux in the NuMI neutrino beam. This reaction is the dominant source of events which only contain high-energy muons nearly parallel to the direction of the neutrino beam. We have isolated a sample of hundreds of such events in neutrino and anti-neutrino enhanced beams, and have constrained the predicted high-energy flux. △ Less

Submitted 23 November, 2021; v1 submitted 2 July, 2021; originally announced July 2021.

Comments: 9 pages, 8 figures, submitted to Physical Review D. Minor corrections and clarifications

Report number: FERMILAB-PUB-21-312-E-SCD

Journal ref: Phys. Rev. D 104, 092010 (2021)

Measurement of inclusive charged-current $ν_μ$ cross sections as a function of muon kinematics at $<E_ν>\sim6~GeV$ on hydrocarbon

Authors: D. Ruterbories, A. Filkins, Z. Ahmad Dar, F. Akbar, D. A. Andrade, M. V. Ascencio, A. Bashyal, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, T. Cai, M. F. Carneiro, G. A. Díaz, H. da Motta, S. A. Dytman, J. Felix, L. Fields, A. M. Gago, H. Gallagher, R. Gran , et al. (38 additional authors not shown)

Abstract: MINERvA presents a new analysis of inclusive charged-current neutrino interactions on a hydrocarbon target. We report single and double-differential cross sections in muon transverse and longitudinal momentum. These measurements are compared to neutrino interaction generator predictions from GENIE, NuWro, GiBUU, and NEUT. In addition, comparisons against models with different treatments of multi-n… ▽ More MINERvA presents a new analysis of inclusive charged-current neutrino interactions on a hydrocarbon target. We report single and double-differential cross sections in muon transverse and longitudinal momentum. These measurements are compared to neutrino interaction generator predictions from GENIE, NuWro, GiBUU, and NEUT. In addition, comparisons against models with different treatments of multi-nucleon correlations, nuclear effects, resonant pion production, and deep inelastic scattering are presented. The data recorded corresponds to $10.61\times10^{20}$ protons on target with a peak neutrino energy of approximately 6 GeV. The higher energy and larger statistics of these data extend the kinematic range for model testing beyond previous MINERvA inclusive charged-current measurements. The results are not well modeled by several generator predictions using a variety of input models. △ Less

Submitted 2 November, 2022; v1 submitted 30 June, 2021; originally announced June 2021.

Comments: 20 pages, 19 figures

Report number: FERMILAB-PUB-21-300-E-QIS-SCD

Use of Neutrino Scattering Events with Low Hadronic Recoil to Inform Neutrino Flux and Detector Energy Scale

Authors: A. Bashyal, D. Rimal, B. Messerly, Z. Ahmad Dar, F. Akbar, M. V. Ascencio, A. Bercellie, M. Betancourt, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, T. Cai, M. F. Carneiro, H. da Motta, S. A. Dytman, G. A. Díaz, J. Felix, L. Fields, A. Filkins, R. Fine, A. M. Gago, H. Gallagher, A. Ghosh , et al. (38 additional authors not shown)

Abstract: Charged-current neutrino interactions with low hadronic recoil ("low-nu") have a cross-section that is approximately constant versus neutrino energy. These interactions have been used to measure the shape of neutrino fluxes as a function of neutrino energy at accelerator-based neutrino experiments such as CCFR, NuTeV, MINOS and MINERvA. In this paper, we demonstrate that low-nu events can be used… ▽ More Charged-current neutrino interactions with low hadronic recoil ("low-nu") have a cross-section that is approximately constant versus neutrino energy. These interactions have been used to measure the shape of neutrino fluxes as a function of neutrino energy at accelerator-based neutrino experiments such as CCFR, NuTeV, MINOS and MINERvA. In this paper, we demonstrate that low-nu events can be used to measure parameters of neutrino flux and detector models and that utilization of event distributions over the upstream detector face can discriminate among parameters that affect the neutrino flux model. From fitting a large sample of low-nu events obtained by exposing MINERvA to the NuMI medium-energy beam, we find that the best-fit flux parameters are within their a priori uncertainties, but the energy scale of muons reconstructed in the MINOS detector is shifted by 3.6% (or 1.8 times the a priori uncertainty on that parameter). These fit results are now used in all MINERvA cross-section measurements, and this technique can be applied by other experiments operating at MINERvA energies, such as DUNE. △ Less

Submitted 17 May, 2022; v1 submitted 12 April, 2021; originally announced April 2021.

Comments: 18 pages, 10 figures, Replaced with accepted version

Report number: FNAL PUB-21-181-ND

An Error Analysis Toolkit for Binned Counting Experiments

Authors: B. Messerly, R. Fine, A. Olivier, Z. Ahmad Dar, F. Akbar, M. V. Ascencio, A. Bashyal, L. Bellantoni, A. Bercellie, J. L. Bonilla, G. Caceres, T. Cai, M. F. Carneiro, G. A. Díaz, J. Felix, L. Fields, A. Filkins, A. Ghosh, S. Gilligan, R. Gran, H. Haider, D. A. Harris, S. Henry, S. Jena, D. Jena , et al. (20 additional authors not shown)

Abstract: We introduce the MINERvA Analysis Toolkit (MAT), a utility for centralizing the handling of systematic uncertainties in HEP analyses. The fundamental utilities of the toolkit are the MnvHnD, a powerful histogram container class, and the systematic Universe classes, which provide a modular implementation of the many universe error analysis approach. These products can be used stand-alone or as part… ▽ More We introduce the MINERvA Analysis Toolkit (MAT), a utility for centralizing the handling of systematic uncertainties in HEP analyses. The fundamental utilities of the toolkit are the MnvHnD, a powerful histogram container class, and the systematic Universe classes, which provide a modular implementation of the many universe error analysis approach. These products can be used stand-alone or as part of a complete error analysis prescription. They support the propagation of systematic uncertainty through all stages of analysis, and provide flexibility for an arbitrary level of user customization. This extensible solution to error analysis enables the standardization of systematic uncertainty definitions across an experiment and a transparent user interface to lower the barrier to entry for new analyzers. △ Less

Submitted 15 March, 2021; originally announced March 2021.

Neutral pion reconstruction using machine learning in the MINERvA experiment at $\langle E_ν\rangle \sim 6$ GeV

Authors: A. Ghosh, B. Yaeggy, R. Galindo, Z. Ahmad Dar, F. Akbar, M. V. Ascencio, A. Bashyal, A. Bercellie, J. L. Bonilla, G. Caceres, T. Cai, M. F. Carneiro, H. da Motta, G. A. Díaz, J. Felix, A. Filkins, R. Fine, A. M. Gago, T. Golan, R. Gran, D. A. Harris, S. Henry, S. Jena, D. Jena, J. Kleykamp , et al. (31 additional authors not shown)

Abstract: This paper presents a novel neutral-pion reconstruction that takes advantage of the machine learning technique of semantic segmentation using MINERvA data collected between 2013-2017, with an average neutrino energy of $6$ GeV. Semantic segmentation improves the purity of neutral pion reconstruction from two gammas from 71\% to 89\% and improves the efficiency of the reconstruction by approximatel… ▽ More This paper presents a novel neutral-pion reconstruction that takes advantage of the machine learning technique of semantic segmentation using MINERvA data collected between 2013-2017, with an average neutrino energy of $6$ GeV. Semantic segmentation improves the purity of neutral pion reconstruction from two gammas from 71\% to 89\% and improves the efficiency of the reconstruction by approximately 40\%. We demonstrate our method in a charged current neutral pion production analysis where a single neutral pion is reconstructed. This technique is applicable to modern tracking calorimeters, such as the new generation of liquid-argon time projection chambers, exposed to neutrino beams with $\langle E_ν\rangle$ between 1-10 GeV. In such experiments it can facilitate the identification of ionization hits which are associated with electromagnetic showers, thereby enabling improved reconstruction of charged-current $ν_e$ events arising from $ν_μ \rightarrow ν_{e}$ appearance. △ Less

Submitted 10 April, 2022; v1 submitted 11 March, 2021; originally announced March 2021.

Comments: 26 pages, v2 matches published version

Journal ref: JINST 16 P07060 2021

Software Training in HEP

Authors: Sudhir Malik, Samuel Meehan, Kilian Lieret, Meirin Oan Evans, Michel H. Villanueva, Daniel S. Katz, Graeme A. Stewart, Peter Elmer, Sizar Aziz, Matthew Bellis, Riccardo Maria Bianchi, Gianluca Bianco, Johan Sebastian Bonilla, Angela Burger, Jackson Burzynski, David Chamont, Matthew Feickert, Philipp Gadow, Bernhard Manfred Gruber, Daniel Guest, Stephan Hageboeck, Lukas Heinrich, Maximilian M. Horzela, Marc Huwiler, Clemens Lange , et al. (22 additional authors not shown)

Abstract: Long term sustainability of the high energy physics (HEP) research software ecosystem is essential for the field. With upgrades and new facilities coming online throughout the 2020s this will only become increasingly relevant throughout this decade. Meeting this sustainability challenge requires a workforce with a combination of HEP domain knowledge and advanced software skills. The required softw… ▽ More Long term sustainability of the high energy physics (HEP) research software ecosystem is essential for the field. With upgrades and new facilities coming online throughout the 2020s this will only become increasingly relevant throughout this decade. Meeting this sustainability challenge requires a workforce with a combination of HEP domain knowledge and advanced software skills. The required software skills fall into three broad groups. The first is fundamental and generic software engineering (e.g. Unix, version control,C++, continuous integration). The second is knowledge of domain specific HEP packages and practices (e.g., the ROOT data format and analysis framework). The third is more advanced knowledge involving more specialized techniques. These include parallel programming, machine learning and data science tools, and techniques to preserve software projects at all scales. This paper dis-cusses the collective software training program in HEP and its activities led by the HEP Software Foundation (HSF) and the Institute for Research and Innovation in Software in HEP (IRIS-HEP). The program equips participants with an array of software skills that serve as ingredients from which solutions to the computing challenges of HEP can be formed. Beyond serving the community by ensuring that members are able to pursue research goals, this program serves individuals by providing intellectual capital and transferable skills that are becoming increasingly important to careers in the realm of software and computing, whether inside or outside HEP △ Less

Submitted 6 August, 2021; v1 submitted 28 February, 2021; originally announced March 2021.

Comments: For CHEP 2021 conference,sent for publication to CSBS Springer

MSC Class: HEP; software; training

Data Preservation at MINERvA

Authors: R. Fine, B. Messerly, K. S. McFarland, S. Akhter, V. Ansari, M. V. Ascencio, J. L. Bonilla, H. da Motta, L. Fields, R. Gran, E. Granados, D. A. Harris, D. Jena, S. Jena, J. Kleykamp, A. Klustová, X. -G. Lu, W. A. Mann, A. Olivier, K. -J. Plows, D. Ruterbories, H. Schellman, H. Su, E. Valencia, B. Yaeggy

Abstract: Between 2013 and 2019 MINERvA collected an accelerator neutrino interaction dataset that is uniquely relevant to the energy range of DUNE. These are the only currently available data at intermediate and high momentum transfers for multiple nuclear targets in the same beam. MINERvA is undertaking a campaign to preserve these data and make them publicly available so that they may be analyzed beyond… ▽ More Between 2013 and 2019 MINERvA collected an accelerator neutrino interaction dataset that is uniquely relevant to the energy range of DUNE. These are the only currently available data at intermediate and high momentum transfers for multiple nuclear targets in the same beam. MINERvA is undertaking a campaign to preserve these data and make them publicly available so that they may be analyzed beyond the end of the MINERvA collaboration. We encourage the community to consider the development of centralized resources to enable long-term access to these data and analysis tools for the entire HEP community. △ Less

Submitted 27 September, 2022; v1 submitted 9 September, 2020; originally announced September 2020.

Comments: Snowmass 2021 Contributed Paper

Double-Differential Inclusive Charged-Current $ν_μ$ Cross Sections on Hydrocarbon in MINERvA at $\langle E_ν \rangle \sim$ 3.5 GeV

Authors: A. Filkins, D. Ruterbories, Y. Liu, Z. Ahmad Dar, F. Akbar, O. Altinok, D. A. Andrade, M. V. Ascencio, A. Bashyal, A. Bercellie, M. Betancourt, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, T. Cai, M. F. Carneiro, H. da Motta, S. A. Dytman, G. A. Díaz, J. Felix, L. Fields, R. Fine, A. M. Gago , et al. (42 additional authors not shown)

Abstract: MINERvA reports inclusive charged-current cross sections for muon neutrinos on hydrocarbon in the NuMI beamline. We measured the double-differential cross section in terms of the longitudinal and transverse muon momenta, as well as the single-differential cross sections in those variables. The data used in this analysis correspond to an exposure of $3.34 \times 10^{20}$ protons on target with a pe… ▽ More MINERvA reports inclusive charged-current cross sections for muon neutrinos on hydrocarbon in the NuMI beamline. We measured the double-differential cross section in terms of the longitudinal and transverse muon momenta, as well as the single-differential cross sections in those variables. The data used in this analysis correspond to an exposure of $3.34 \times 10^{20}$ protons on target with a peak neutrino energy of approximately 3.5 GeV. Measurements are compared to the GENIE, NuWro and GiBUU neutrino cross-section predictions, as well as a version of GENIE modified to produce better agreement with prior exclusive MINERvA measurements. None of the models or variants were able to successfully reproduce the data across the entire phase space, which includes areas dominated by each interaction channel. △ Less

Submitted 23 June, 2020; v1 submitted 27 February, 2020; originally announced February 2020.

Comments: 20 pages, 16 figures

Report number: FERMILAB-PUB-20-098-ND-SCD

Journal ref: Phys. Rev. D 101, 112007 (2020)

Monotonicity-preserving finite element schemes with adaptive mesh refinement for hyperbolic problems

Authors: Jesus Bonilla, Santiago Badia

Abstract: This work is focused on the extension and assessment of the monotonicity-preserving scheme in [3] and the local bounds preserving scheme in [5] to hierarchical octree adaptive mesh refinement (AMR). Whereas the former can readily be used on this kind of meshes, the latter requires some modifications. A key question that we want to answer in this work is whether to move from a linear to a nonlinear… ▽ More This work is focused on the extension and assessment of the monotonicity-preserving scheme in [3] and the local bounds preserving scheme in [5] to hierarchical octree adaptive mesh refinement (AMR). Whereas the former can readily be used on this kind of meshes, the latter requires some modifications. A key question that we want to answer in this work is whether to move from a linear to a nonlinear stabilization mechanism pays the price when combined with shock-adapted meshes. Whereas nonlinear (or shock-capturing) stabilization leads to improved accuracy compared to linear schemes, it also negatively hinders nonlinear convergence, increasing computational cost. We compare linear and nonlinear schemes in terms of the required computational time versus accuracy for several steady benchmark problems. Numerical results indicate that, in general, nonlinear schemes can be cost-effective for sufficiently refined meshes. Besides, it is also observed that it is better to refine further around shocks rather than use sharper shock capturing terms, which usually yield stiffer nonlinear problems. In addition, a new refinement criterion has been proposed. The proposed criterion is based on the graph Laplacian used in the definition of the stabilization method. Numerical results show that this shock detector performs better than the well-known Kelly estimator for problems with shocks or discontinuities. △ Less

Submitted 27 March, 2020; v1 submitted 24 December, 2019; originally announced December 2019.

Comments: arXiv admin note: text overlap with arXiv:1911.06792

On differentiable local bounds preserving stabilization for Euler equations

Authors: Santiago Badia, Jesús Bonilla, Sibusiso Mabuza, John N. Shadid

Abstract: This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A differentiable local bounds preserving method has been developed, which combines a Rusanov artificial diffusion operator and a differentiable shock detector. Nonlinear s… ▽ More This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A differentiable local bounds preserving method has been developed, which combines a Rusanov artificial diffusion operator and a differentiable shock detector. Nonlinear stabilization schemes are usually stiff and highly nonlinear. This issue is mitigated by the differentiability properties of the proposed method. Moreover, in order to further improve the nonlinear convergence, we also propose a continuation method for a subset of the stabilization parameters. The resulting method has been successfully applied to steady and transient problems with complex shock patterns. Numerical experiments show that it is able to provide sharp and well resolved shocks. The importance of the differentiability is assessed by comparing the new scheme with its non-differentiable counterpart. Numerical experiments suggest that, for up to moderate nonlinear tolerances, the method exhibits improved robustness and nonlinear convergence behavior for steady problems. In the case of transient problem, we also observe a reduction in the computational cost. △ Less

Submitted 7 July, 2020; v1 submitted 15 November, 2019; originally announced November 2019.

Comments: arXiv admin note: text overlap with arXiv:1912.11487

Nuclear binding energy and transverse momentum imbalance in neutrino-nucleus reactions

Authors: T. Cai, X. -G. Lu, L. A. Harewood, C. Wret, F. Akbar, D. A. Andrade, M. V. Ascencio, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, M. F. Carneiro, D. Coplowe, H. da Motta, Zubair Ahmad Dar, G. A. Díaz, J. Felix, L. Fields, A. Filkins, R. Fine, A. M. Gago , et al. (42 additional authors not shown)

Abstract: We have measured new observables based on the final state kinematic imbalances in the mesonless production of $ν_μ+A\rightarrowμ^-+p+X$ in the $\text{MINER}ν\text{A}$ tracker. Components of the muon-proton momentum imbalances parallel ($δp_\mathrm{Ty}$) and perpendicular($δp_\mathrm{Tx}$) to the momentum transfer in the transverse plane are found to be sensitive to the nuclear effects such as Ferm… ▽ More We have measured new observables based on the final state kinematic imbalances in the mesonless production of $ν_μ+A\rightarrowμ^-+p+X$ in the $\text{MINER}ν\text{A}$ tracker. Components of the muon-proton momentum imbalances parallel ($δp_\mathrm{Ty}$) and perpendicular($δp_\mathrm{Tx}$) to the momentum transfer in the transverse plane are found to be sensitive to the nuclear effects such as Fermi motion, binding energy and non-QE contributions. The QE peak location in $δp_\mathrm{Ty}$ is particularly sensitive to the binding energy. Differential cross sections are compared to predictions from different neutrino interaction models. The Fermi gas models presented in this study cannot simultaneously describe features such as QE peak location, width and the non-QE events contributing to the signal process. Correcting the GENIE's binding energy implementation according to theory causes better agreement with data. Hints of proton left-right asymmetry are observed in $δp_\mathrm{Tx}$. Better modeling of the binding energy can reduce bias in neutrino energy reconstruction and these observables can be applied in current and future experiments to better constrain nuclear effects. △ Less

Submitted 3 May, 2020; v1 submitted 18 October, 2019; originally announced October 2019.

Comments: 14 pages, 12 figures; Added data release file in .root format

Report number: FERMILAB-PUB-19-539-ND

Journal ref: Phys. Rev. D 101, 092001 (2020)

Measurement of $\barν_μ$ charged-current single $π^{-}$ production on hydrocarbon in the few-GeV region using MINERvA

Authors: T. Le, F. Akbar, L. Aliaga, D. A. Andrade, M. V. Ascencio, A. Bashyal, A. Bercellie, M. Betancourt, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, T. Cai, M. F. Carneiro, D. Coplowe, S. A. Dytman, G. A. Díaz, 5 J. Felix, L. Fields, A. Filkins, R. Fine, N. Fiza, A. M. Gago, H. Gallagher , et al. (41 additional authors not shown)

Abstract: The antineutrino scattering channel $\barν_μ \,\text{CH} \rightarrow μ^{+} \,π^{-} \,X$(nucleon(s)) is analyzed in the incident energy range 1.5 to 10 GeV using the MINERvA detector at Fermilab. Differential cross sections are reported as functions of $μ^{+}$ momentum and production angle, $π^{-}$ kinetic energy and production angle, and antineutrino energy and squared four-momentum transfer. Dist… ▽ More The antineutrino scattering channel $\barν_μ \,\text{CH} \rightarrow μ^{+} \,π^{-} \,X$(nucleon(s)) is analyzed in the incident energy range 1.5 to 10 GeV using the MINERvA detector at Fermilab. Differential cross sections are reported as functions of $μ^{+}$ momentum and production angle, $π^{-}$ kinetic energy and production angle, and antineutrino energy and squared four-momentum transfer. Distribution shapes are generally reproduced by simulations based on the GENIE, NuWro, and GiBUU event generators, however GENIE (GiBUU) overestimates (underestimates) the cross-section normalizations by 8% (10%). Comparisons of data with the GENIE-based reference simulation probe conventional treatments of cross sections and pion intranuclear rescattering. The distribution of non-track vertex energy is used to decompose the signal sample into reaction categories, and cross sections are determined for the exclusive reactions $μ^{+} π^{-} n$ and $ μ^+ π^{-} p$. A similar treatment applied to the published MINERvA sample $\barν_μ \,\text{CH} \rightarrow μ^{+} \,π^{0} \,X$(nucleon(s)) has determined the $μ^{+} π^{0} n$ cross section, and the latter is used with $σ(π^{-} n)$ and $σ(π^{-} p)$ to carry out an isospin decomposition of $\barν_μ$-induced CC($π$). The ratio of magnitudes and relative phase for isospin amplitudes $A_{3}$ and $A_{1}$ thereby obtained are: $R^{\barν} = 0.99 \pm 0.19$ and $φ^{\barν} = 93^{\circ} \pm 7^{\circ}$. Our results are in agreement with bubble chamber measurements made four decades ago. △ Less

Submitted 27 August, 2019; v1 submitted 19 June, 2019; originally announced June 2019.

Comments: 22 pages, 17 figures

Report number: FERMILAB-PUB-18-578-ND

Journal ref: Phys. Rev. D 100, 052008 (2019)

Constraint of the MINERvA Medium Energy Neutrino Flux using Neutrino-Electron Elastic Scattering

Authors: E. Valencia, D. Jena, Nuruzzaman, F. Akbar, L. Aliaga, D. A. Andrade, M. V. Ascencio, A. Bashyal, L. Bellantoni, A. Bercellie, A. Bodek, J. L. Bonilla, A. Bravar, H. Budd, G. Caceres, T. Cai, M. F. Carneiro, J. Chaves, D. Coplowe, H. da Motta, S. A. Dytman, G. A. Diaz, J. Felix, L. Fields, A. Filkins , et al. (46 additional authors not shown)

Abstract: Elastic neutrino scattering on electrons is a precisely-known purely leptonic process that provides a standard candle for measuring neutrino flux in conventional neutrino beams. Using a total sample of 810 neutino-electron scatters after background subtraction, the measurement reduces the normalization uncertainty on the muon neutrino NuMI flux between 2 and 20 GeV from 7.5% to 3.9%. This is the m… ▽ More Elastic neutrino scattering on electrons is a precisely-known purely leptonic process that provides a standard candle for measuring neutrino flux in conventional neutrino beams. Using a total sample of 810 neutino-electron scatters after background subtraction, the measurement reduces the normalization uncertainty on the muon neutrino NuMI flux between 2 and 20 GeV from 7.5% to 3.9%. This is the most precise measurement of neutrino-electron scattering to date, will reduce uncertainties on MINERvA's absolute cross section measurements, and demonstrates a technique that can be used in future neutrino beams such as LBNF. △ Less

Submitted 24 February, 2022; v1 submitted 31 May, 2019; originally announced June 2019.

Comments: 12 pages, 12 figures. Final accepted version, including ancillary files with flux data tables. Ancillary data files updated October 2021 to include the neutrino-electron measured event rate and covariance

Report number: FERMILAB-PUB-19-256-CD-ND

Journal ref: Phys. Rev. D 100, 092001 (2019)

Renormalization of a Finite Range Inverse Cube Potential

Authors: Daniel Odell, Arnoldas Deltuva, Jose Bonilla, Lucas Platter

Abstract: We study the regularization and renormalization of a finite range inverse cube potential in the two- and three-body sectors. Specifically, we compare and contrast three different regulation schemes frequently used to study few-body systems as well as the associated renormalization group flows. We also calculate bound state and scattering observables over a wide range of cutoffs, demonstrating the… ▽ More We study the regularization and renormalization of a finite range inverse cube potential in the two- and three-body sectors. Specifically, we compare and contrast three different regulation schemes frequently used to study few-body systems as well as the associated renormalization group flows. We also calculate bound state and scattering observables over a wide range of cutoffs, demonstrating the sufficiency of a two-body contact interaction to renormalize two- and three-body observables. We supplement these plots with quantified analyses of the observables' residual cutoff dependence. △ Less

Submitted 20 September, 2019; v1 submitted 28 February, 2019; originally announced March 2019.

Comments: 26 pages, 14 figures

Journal ref: Phys. Rev. C 100, 054001 (2019)

Maximum-principle preserving space-time isogeometric analysis

Authors: Jesús Bonilla, Santiago Badia

Abstract: In this work we propose a nonlinear stabilization technique for convection-diffusion-reaction and pure transport problems discretized with space-time isogeometric analysis. The stabilization is based on a graph-theoretic artificial diffusion operator and a novel shock detector for isogeometric analysis. Stabilization in time and space directions are performed similarly, which allow us to use high-… ▽ More In this work we propose a nonlinear stabilization technique for convection-diffusion-reaction and pure transport problems discretized with space-time isogeometric analysis. The stabilization is based on a graph-theoretic artificial diffusion operator and a novel shock detector for isogeometric analysis. Stabilization in time and space directions are performed similarly, which allow us to use high-order discretizations in time without any CFL-like condition. The method is proven to yield solutions that satisfy the discrete maximum principle (DMP) unconditionally for arbitrary order. In addition, the stabilization is linearity preserving in a space-time sense. Moreover, the scheme is proven to be Lipschitz continuous ensuring that the nonlinear problem is well-posed. Solving large problems using a space-time discretization can become highly costly. Therefore, we also propose a partitioned space-time scheme that allows us to select the length of every time slab, and solve sequentially for every subdomain. As a result, the computational cost is reduced while the stability and convergence properties of the scheme remain unaltered. In addition, we propose a twice differentiable version of the stabilization scheme, which enjoys the same stability properties while the nonlinear convergence is significantly improved. Finally, the proposed schemes are assessed with numerical experiments. In particular, we considered steady and transient pure convection and convection-diffusion problems in one and two dimensions. △ Less

Submitted 15 November, 2019; v1 submitted 13 December, 2018; originally announced December 2018.

A Novel Heat Exchanger Design Method Using a Delayed Rejection Adaptive Metropolis Hasting Algorithm

Authors: Ahad Mohammadi, Javier Bonilla, Reza Zarghami, Shahab Golshan

Abstract: In this study, a shell-and-tube heat exchanger (STHX) design based on seven continuous independent design variables is proposed. Delayed Rejection Adaptive Metropolis hasting (DRAM) was utilized as a powerful tool in the Markov chain Monte Carlo (MCMC) sampling method. This Reverse Sampling (RS) method was used to find the probability distribution of design variables of the shell and tube heat exc… ▽ More In this study, a shell-and-tube heat exchanger (STHX) design based on seven continuous independent design variables is proposed. Delayed Rejection Adaptive Metropolis hasting (DRAM) was utilized as a powerful tool in the Markov chain Monte Carlo (MCMC) sampling method. This Reverse Sampling (RS) method was used to find the probability distribution of design variables of the shell and tube heat exchanger. Thanks to this probability distribution, an uncertainty analysis was also performed to find the quality of these variables. In addition, a decision-making strategy based on confidence intervals of design variables and on the Total Annual Cost (TAC) provides the final selection of design variables. Results indicated high accuracies for the estimation of design variables which leads to marginally improved performance compared to commonly used optimization methods. In order to verify the capability of the proposed method, a case of study is also presented, it shows that a significant cost reduction is feasible with respect to multi-objective and single-objective optimization methods. Furthermore, the selected variables have good quality (in terms of probability distribution) and a lower TAC was also achieved. Results show that the costs of the proposed design are lower than those obtained from optimization method reported in previous studies. The algorithm was also used to determine the impact of using probability values for the design variables rather than single values to obtain the best heat transfer area and pumping power. In particular, a reduction of the TAC up to 3.5% was achieved in the case considered. △ Less

Submitted 19 September, 2018; originally announced October 2018.

Journal ref: Applied Thermal Engineering, Volume 137, 2018, Pages 808-821

New approach for solar tracking systems based on computer vision, low cost hardware and deep learning

Authors: Jose A. Carballo, Javier Bonilla, Manuel Berenguel, Jesús Fernández-Reche, Ginés García

Abstract: In this work, a new approach for Sun tracking systems is presented. Due to the current system limitations regarding costs and operational problems, a new approach based on low cost, computer vision open hardware and deep learning has been developed. The preliminary tests carried out successfully in Plataforma solar de Almeria (PSA), reveal the great potential and show the new approach as a good al… ▽ More In this work, a new approach for Sun tracking systems is presented. Due to the current system limitations regarding costs and operational problems, a new approach based on low cost, computer vision open hardware and deep learning has been developed. The preliminary tests carried out successfully in Plataforma solar de Almeria (PSA), reveal the great potential and show the new approach as a good alternative to traditional systems. The proposed approach can provide key variables for the Sun tracking system control like cloud movements prediction, block and shadow detection, atmospheric attenuation or measures of concentrated solar radiation, which can improve the control strategies of the system and therefore the system performance. △ Less

Submitted 19 September, 2018; originally announced September 2018.

Comments: 12 pages, 10 figures,

MSC Class: 68T45

Journal ref: Carballo, J. A., Bonilla, J., Berenguel, M., Fernandez-Reche, J., & Garcia, G. (2018). New approach for solar tracking systems based on computer vision, low cost hardware and deep learning. Renewable Energy

Simulation on various topocentric points of the zodiac constellations

Authors: Harold Yesid Laserna Díaz, Sergio Libardo Silva Sandobal, Yesid Javier Cruz Bonilla

Abstract: On the sighting of the constellations and their stars, they seem to be near each other at first sight but some of thema are actually quite further from the others. Froma the geocentric point view can the Zodiac Costellations be perceived as referential marks divided into terrestrial ecliptic, which are organized according to their shapes. In this paper are these star groups simulated in order to c… ▽ More On the sighting of the constellations and their stars, they seem to be near each other at first sight but some of thema are actually quite further from the others. Froma the geocentric point view can the Zodiac Costellations be perceived as referential marks divided into terrestrial ecliptic, which are organized according to their shapes. In this paper are these star groups simulated in order to compare the shapes of these groups and the diverse coordinate system of the astronomy. Within the simulation is the distance measure light-year remplaced with a so-called program's default "unity" that uses a compilation in Python, where such distances are escalated to 1:10 △ Less

Submitted 27 December, 2017; originally announced December 2017.

Differentiable monotonicity-preserving schemes for discontinuous Galerkin methods on arbitrary meshes

Authors: Santiago Badia, Jesús Bonilla, Alba Hierro

Abstract: This work is devoted to the design of interior penalty discontinuous Galerkin (dG) schemes that preserve maximum principles at the discrete level for the steady transport and convection-diffusion problems and the respective transient problems with implicit time integration. Monotonic schemes that combine explicit time stepping with dG space discretization are very common, but the design of such sc… ▽ More This work is devoted to the design of interior penalty discontinuous Galerkin (dG) schemes that preserve maximum principles at the discrete level for the steady transport and convection-diffusion problems and the respective transient problems with implicit time integration. Monotonic schemes that combine explicit time stepping with dG space discretization are very common, but the design of such schemes for implicit time stepping is rare, and it had only been attained so far for 1D problems. The proposed scheme is based on an artificial diffusion that linearly depends on a shock detector that identifies the troublesome areas. In order to define the new shock detector, we have introduced the concept of discrete local extrema. The diffusion operator is a graph-Laplacian, instead of the more common finite element discretization of the Laplacian operator, which is essential to keep monotonicity on general meshes and in multi-dimension. The resulting nonlinear stabilization is non-smooth and nonlinear solvers can fail to converge. As a result, we propose a smoothed (twice differentiable) version of the nonlinear stabilization, which allows us to use Newton with line search nonlinear solvers and dramatically improve nonlinear convergence. A theoretical numerical analysis of the proposed schemes show that they satisfy the desired monotonicity properties. Further, the resulting operator is Lipschitz continuous and there exists at least one solution of the discrete problem, even in the non-smooth version. We provide a set of numerical results to support our findings. △ Less

Submitted 13 December, 2018; v1 submitted 27 December, 2016; originally announced December 2016.

Automatic Identification of Scenedesmus Polymorphic Microalgae from Microscopic Images

Authors: Jhony-Heriberto Giraldo-Zuluaga, Geman Diez, Alexander Gomez, Tatiana Martinez, Mariana Peñuela Vasquez, Jesus Francisco Vargas Bonilla, Augusto Salazar

Abstract: Microalgae counting is used to measure biomass quantity. Usually, it is performed in a manual way using a Neubauer chamber and expert criterion, with the risk of a high error rate. This paper addresses the methodology for automatic identification of Scenedesmus microalgae (used in the methane production and food industry) and applies it to images captured by a digital microscope. The use of contra… ▽ More Microalgae counting is used to measure biomass quantity. Usually, it is performed in a manual way using a Neubauer chamber and expert criterion, with the risk of a high error rate. This paper addresses the methodology for automatic identification of Scenedesmus microalgae (used in the methane production and food industry) and applies it to images captured by a digital microscope. The use of contrast adaptive histogram equalization for pre-processing, and active contours for segmentation are presented. The calculation of statistical features (Histogram of Oriented Gradients, Hu and Zernike moments) with texture features (Haralick and Local Binary Patterns descriptors) are proposed for algae characterization. Scenedesmus algae can build coenobia consisting of 1, 2, 4 and 8 cells. The amount of algae of each coenobium helps to determine the amount of lipids, proteins, and other substances in a given sample of a algae crop. The knowledge of the quantity of those elements improves the quality of bioprocess applications. Classification of coenobia achieves accuracies of 98.63% and 97.32% with Support Vector Machine (SVM) and Artificial Neural Network (ANN), respectively. According to the results it is possible to consider the proposed methodology as an alternative to the traditional technique for algae counting. The database used in this paper is publicly available for download. △ Less

Submitted 23 October, 2017; v1 submitted 21 December, 2016; originally announced December 2016.

Comments: This is a pre-print of an article published in Pattern Analysis and Applications. The final authenticated version is available online at: https://doi.org/10.1007/s10044-017-0662-3, Pattern Anal Applic (2017)

Monotonicity-preserving finite element schemes based on differentiable nonlinear stabilization

Authors: Santiago Badia, Jesús Bonilla

Abstract: In this work, we propose a nonlinear stabilization technique for scalar conservation laws with implicit time stepping. The method relies on an artificial diffusion method, based on a graph-Laplacian operator. It is nonlinear, since it depends on a shock detector. The same shock detector is used to gradually lump the mass matrix. The resulting method is LED, positivity preserving, linearity preserv… ▽ More In this work, we propose a nonlinear stabilization technique for scalar conservation laws with implicit time stepping. The method relies on an artificial diffusion method, based on a graph-Laplacian operator. It is nonlinear, since it depends on a shock detector. The same shock detector is used to gradually lump the mass matrix. The resulting method is LED, positivity preserving, linearity preserving, and also satisfies a global DMP. Lipschitz continuity has also been proved. However, the resulting scheme is highly nonlinear, leading to very poor nonlinear convergence rates. We propose a smooth version of the scheme, which leads to twice differentiable nonlinear stabilization schemes. It allows one to straightforwardly use Newton's method and obtain quadratic convergence. In the numerical experiments, steady and transient linear transport, and transient Burgers' equation have been considered in 2D. Using the Newton method with a smooth version of the scheme we can reduce 10 to 20 times the number of iterations of Anderson acceleration with the original non-smooth scheme. In any case, these properties are only true for the converged solution, but not for iterates. In this sense, we have also proposed the concept of projected nonlinear solvers, where a projection step is performed at the end of every nonlinear iteration onto a FE space of admissible solutions. The space of admissible solutions is the one that satisfies the desired monotonic properties (maximum principle or positivity). △ Less

Submitted 22 December, 2016; v1 submitted 28 June, 2016; originally announced June 2016.