Showing 1–50 of 150 results for author: Bergman, D

Network Relaxations for Discrete Bilevel Optimization under Linear Interactions

Authors: Leonardo Lozano, David Bergman, Andre Augusto Cire

Abstract: We investigate relaxations for a class of discrete bilevel programs where the interaction constraints linking the leader and the follower are linear. Our approach reformulates the upper-level optimality constraints by projecting the leader's decisions onto vectors that map to distinct follower solution values, each referred to as a state. Based on such a state representation, we develop a network-… ▽ More We investigate relaxations for a class of discrete bilevel programs where the interaction constraints linking the leader and the follower are linear. Our approach reformulates the upper-level optimality constraints by projecting the leader's decisions onto vectors that map to distinct follower solution values, each referred to as a state. Based on such a state representation, we develop a network-flow linear program via a decision diagram that captures the convex hull of the follower's value function graph, leading to a new single-level reformulation of the bilevel problem. We also present a reduction procedure that exploits symmetry to identify the reformulation of minimal size. For large networks, we introduce parameterized relaxations that aggregate states by considering tractable hyperrectangles based on lower and upper bounds associated with the interaction constraints, and can be integrated into existing mixed-integer bilevel linear programming (MIBLP) solvers. Numerical experiments suggest that the new relaxations, whether used within a simple cutting-plane procedure or integrated into state-of-the-art MIBLP solvers, significantly reduce runtimes or solve additional benchmark instances. Our findings also highlight the correlation between the quality of relaxations and the properties of the interaction matrix, underscoring the potential of our approach in enhancing solution methods for structured bilevel optimization instances. △ Less

Submitted 25 July, 2024; originally announced July 2024.

The Madness of Multiple Entries in March Madness

Authors: Jeff Decary, David Bergman, Carlos Cardonha, Jason Imbrogno, Andrea Lodi

Abstract: This paper explores multi-entry strategies for betting pools related to single-elimination tournaments. In such betting pools, participants select winners of games, and their respective score is a weighted sum of the number of correct selections. Most betting pools have a top-heavy payoff structure, so the paper focuses on strategies that maximize the expected score of the best-performing entry. T… ▽ More This paper explores multi-entry strategies for betting pools related to single-elimination tournaments. In such betting pools, participants select winners of games, and their respective score is a weighted sum of the number of correct selections. Most betting pools have a top-heavy payoff structure, so the paper focuses on strategies that maximize the expected score of the best-performing entry. There is no known closed-formula expression for the estimation of this metric, so the paper investigates the challenges associated with the estimation and the optimization of multi-entry solutions. We present an exact dynamic programming approach for calculating the maximum expected score of any given fixed solution, which is exponential in the number of entries. We explore the structural properties of the problem to develop several solution techniques. In particular, by extracting insights from the solutions produced by one of our algorithms, we design a simple yet effective problem-specific heuristic that was the best-performing technique in our experiments, which were based on real-world data extracted from recent March Madness tournaments. In particular, our results show that the best 100-entry solution identified by our heuristic had a 2.2% likelihood of winning a $1 million prize in a real-world betting pool. △ Less

Submitted 18 July, 2024; originally announced July 2024.

MSC Class: 90-10 (Primary) 91A60 (Secondary)

Isotropy of cosmic rays beyond $10^{20}$ eV favors their heavy mass composition

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the resul… ▽ More We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the results are consistent with a relatively heavy injected composition at E ~ 10 EeV that becomes lighter up to E ~ 100 EeV, while the composition at E > 100 EeV is very heavy. The latter is true even in the presence of highest experimentally allowed extra-galactic magnetic fields, while the composition at lower energies can be light if a strong EGMF is present. The effect of the uncertainty in the galactic magnetic field on these results is subdominant. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: 8 pages, 3 figures, accepted for publication in PRL

Mass composition of ultra-high energy cosmic rays from distribution of their arrival directions with the Telescope Array

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale struc… ▽ More We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale structure (LSS) of the Universe. As we report in the companion letter, the TA data show large deflections with respect to the LSS which can be explained, assuming small extra-galactic magnetic fields (EGMF), by an intermediate composition changing to a heavy one (iron) in the highest energy bin. Here we show that these results are robust to uncertainties in UHECR injection spectra, the energy scale of the experiment and galactic magnetic fields (GMF). The assumption of weak EGMF, however, strongly affects this interpretation at all but the highest energies E > 100 EeV, where the remarkable isotropy of the data implies a heavy injected composition even in the case of strong EGMF. This result also holds if UHECR sources are as rare as $2 \times 10^{-5}$ Mpc$^{-3}$, that is the conservative lower limit for the source number density. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: 18 pages, 11 figures, accepted for publication in PRD

Observation of Declination Dependence in the Cosmic Ray Energy Spectrum

Authors: The Telescope Array Collaboration, R. U. Abbasi, T. Abu-Zayyad, M. Allen, J. W. Belz, D. R. Bergman, I. Buckland, W. Campbell, B. G. Cheon, K. Endo, A. Fedynitch, T. Fujii, K. Fujisue, K. Fujita, M. Fukushima, G. Furlich, Z. Gerber, N. Globus, W. Hanlon, N. Hayashida, H. He, K. Hibino, R. Higuchi, D. Ikeda, T. Ishii , et al. (101 additional authors not shown)

Abstract: We report on an observation of the difference between northern and southern skies of the ultrahigh energy cosmic ray energy spectrum with a significance of ${\sim}8σ$. We use measurements from the two largest experiments$\unicode{x2014}$the Telescope Array observing the northern hemisphere and the Pierre Auger Observatory viewing the southern hemisphere. Since the comparison of two measurements fr… ▽ More We report on an observation of the difference between northern and southern skies of the ultrahigh energy cosmic ray energy spectrum with a significance of ${\sim}8σ$. We use measurements from the two largest experiments$\unicode{x2014}$the Telescope Array observing the northern hemisphere and the Pierre Auger Observatory viewing the southern hemisphere. Since the comparison of two measurements from different observatories introduces the issue of possible systematic differences between detectors and analyses, we validate the methodology of the comparison by examining the region of the sky where the apertures of the two observatories overlap. Although the spectra differ in this region, we find that there is only a $1.8σ$ difference between the spectrum measurements when anisotropic regions are removed and a fiducial cut in the aperture is applied. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 8 pages, 6 figures

MORBDD: Multiobjective Restricted Binary Decision Diagrams by Learning to Sparsify

Authors: Rahul Patel, Elias B. Khalil, David Bergman

Abstract: In multicriteria decision-making, a user seeks a set of non-dominated solutions to a (constrained) multiobjective optimization problem, the so-called Pareto frontier. In this work, we seek to bring a state-of-the-art method for exact multiobjective integer linear programming into the heuristic realm. We focus on binary decision diagrams (BDDs) which first construct a graph that represents all feas… ▽ More In multicriteria decision-making, a user seeks a set of non-dominated solutions to a (constrained) multiobjective optimization problem, the so-called Pareto frontier. In this work, we seek to bring a state-of-the-art method for exact multiobjective integer linear programming into the heuristic realm. We focus on binary decision diagrams (BDDs) which first construct a graph that represents all feasible solutions to the problem and then traverse the graph to extract the Pareto frontier. Because the Pareto frontier may be exponentially large, enumerating it over the BDD can be time-consuming. We explore how restricted BDDs, which have already been shown to be effective as heuristics for single-objective problems, can be adapted to multiobjective optimization through the use of machine learning (ML). MORBDD, our ML-based BDD sparsifier, first trains a binary classifier to eliminate BDD nodes that are unlikely to contribute to Pareto solutions, then post-processes the sparse BDD to ensure its connectivity via optimization. Experimental results on multiobjective knapsack problems show that MORBDD is highly effective at producing very small restricted BDDs with excellent approximation quality, outperforming width-limited restricted BDDs and the well-known evolutionary algorithm NSGA-II. △ Less

Submitted 4 March, 2024; originally announced March 2024.

BDD for Complete Characterization of a Safety Violation in Linear Systems with Inputs

Authors: Manish Goyal, David Bergman, Parasara Sridhar Duggirala

Abstract: The control design tools for linear systems typically involves pole placement and computing Lyapunov functions which are useful for ensuring stability. But given higher requirements on control design, a designer is expected to satisfy other specification such as safety or temporal logic specification as well, and a naive control design might not satisfy such specification. A control designer can e… ▽ More The control design tools for linear systems typically involves pole placement and computing Lyapunov functions which are useful for ensuring stability. But given higher requirements on control design, a designer is expected to satisfy other specification such as safety or temporal logic specification as well, and a naive control design might not satisfy such specification. A control designer can employ model checking as a tool for checking safety and obtain a counterexample in case of a safety violation. While several scalable techniques for verification have been developed for safety verification of linear dynamical systems, such tools merely act as decision procedures to evaluate system safety and, consequently, yield a counterexample as an evidence to safety violation. However these model checking methods are not geared towards discovering corner cases or re-using verification artifacts for another sub-optimal safety specification. In this paper, we describe a technique for obtaining complete characterization of counterexamples for a safety violation in linear systems. The proposed technique uses the reachable set computed during safety verification for a given temporal logic formula, performs constraint propagation, and represents all modalities of counterexamples using a binary decision diagram (BDD). We introduce an approach to dynamically determine isomorphic nodes for obtaining a considerably reduced (in size) decision diagram. A thorough experimental evaluation on various benchmarks exhibits that the reduction technique achieves up to $67\%$ reduction in the number of nodes and $75\%$ reduction in the width of the decision diagram. △ Less

Submitted 26 November, 2023; originally announced November 2023.

Comments: 16 pages, 5 figures, 2 tables

Neutrino propagation through Earth: modeling uncertainties using nuPyProp

Authors: Diksha Garg, Mary Hall Reno, Sameer Patel, Alexander Ruestle, Yosui Akaike, Luis A. Anchordoqui, Douglas R. Bergman, Isaac Buckland, Austin L. Cummings, Johannes Eser, Fred Garcia, Claire Guépin, Tobias Heibges, Andrew Ludwig, John F. Krizmanic, Simon Mackovjak, Eric Mayotte, Sonja Mayotte, Angela V. Olinto, Thomas C. Paul, Andrés Romero-Wolf, Frédéric Sarazin, Tonia M. Venters, Lawrence Wiencke, Stephanie Wissel

Abstract: Using the Earth as a neutrino converter, tau neutrino fluxes from astrophysical point sources can be detected by tau-lepton-induced extensive air showers (EASs). Both muon neutrino and tau neutrino induced upward-going EAS signals can be detected by terrestrial, sub-orbital and satellite-based instruments. The sensitivity of these neutrino telescopes can be evaluated with the nuSpaceSim package, w… ▽ More Using the Earth as a neutrino converter, tau neutrino fluxes from astrophysical point sources can be detected by tau-lepton-induced extensive air showers (EASs). Both muon neutrino and tau neutrino induced upward-going EAS signals can be detected by terrestrial, sub-orbital and satellite-based instruments. The sensitivity of these neutrino telescopes can be evaluated with the nuSpaceSim package, which includes the nuPyProp simulation package. The nuPyProp package propagates neutrinos ($ν_μ$, $ν_τ$) through the Earth to produce the corresponding charged leptons (muons and tau-leptons). We use nuPyProp to quantify the uncertainties from Earth density models, tau depolarization effects and photo-nuclear electromagnetic energy loss models in the charged lepton exit probabilities and their spectra. The largest uncertainties come from electromagnetic energy loss modeling, with as much as a 20-50% difference between the models. We compare nuPyProp results with other simulation package results. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: 9 pages, 6 figures, Proceedings for International Cosmic Ray Conference (ICRC) 2023

Universality of Cherenkov Light in EAS

Authors: Isaac Buckland, Douglas Bergman

Abstract: The reconstruction of cosmic-ray-induced extensive air showers with a non-imaging Cherenkov detector array requires knowledge of the Cherenkov yield of any given air shower for a given set of shower parameters. Although air showers develop in a stochastic cascade, certain characteristics of the particles in the shower have been shown to come from universal probability distributions, a property kno… ▽ More The reconstruction of cosmic-ray-induced extensive air showers with a non-imaging Cherenkov detector array requires knowledge of the Cherenkov yield of any given air shower for a given set of shower parameters. Although air showers develop in a stochastic cascade, certain characteristics of the particles in the shower have been shown to come from universal probability distributions, a property known as shower universality. Both the energy and the angular distributions of charged particles within a shower have been parameterized. One can use these distributions to calculate the Cherenkov photon yield as an angular distribution from the Cherenkov cones of charged particles at various stages of shower development. This Cherenkov photon yield can then be tabulated for use in the reconstruction of air showers. In this work, we develop the calculation of both the Cherenkov angular distribution and Cherenkov yield per shower particle, and show how a look-up table was constructed to capture the relevant features of these distributions for general use. We compare the results of our calculations with the results of full, particle-stack, Monte Carlo simulation of the Cherenkov light produced in extensive air showers using CORSIKA-IACT. We make comparisons of both the lateral distribution of the Cherenkov photon flux amongst several detectors and of the arrival-time distribution of the Cherenkov photons in a single detector. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Neutrino propagation in the Earth and emerging charged leptons with $\texttt{nuPyProp}$

Authors: Diksha Garg, Sameer Patel, Mary Hall Reno, Alexander Reustle, Yosui Akaike, Luis A. Anchordoqui, Douglas R. Bergman, Isaac Buckland, Austin L. Cummings, Johannes Eser, Fred Garcia, Claire Guépin, Tobias Heibges, Andrew Ludwig, John F. Krizmanic, Simon Mackovjak, Eric Mayotte, Sonja Mayotte, Angela V. Olinto, Thomas C. Paul, Andrés Romero-Wolf, Frédéric Sarazin, Tonia M. Venters, Lawrence Wiencke, Stephanie Wissel

Abstract: Ultra-high-energy neutrinos serve as messengers of some of the highest energy astrophysical environments. Given that neutrinos are neutral and only interact via weak interactions, neutrinos can emerge from sources, traverse astronomical distances, and point back to their origins. Their weak interactions require large target volumes for neutrino detection. Using the Earth as a neutrino converter, t… ▽ More Ultra-high-energy neutrinos serve as messengers of some of the highest energy astrophysical environments. Given that neutrinos are neutral and only interact via weak interactions, neutrinos can emerge from sources, traverse astronomical distances, and point back to their origins. Their weak interactions require large target volumes for neutrino detection. Using the Earth as a neutrino converter, terrestrial, sub-orbital, and satellite-based instruments are able to detect signals of neutrino-induced extensive air showers. In this paper, we describe the software code $\texttt{nuPyProp}$ that simulates tau neutrino and muon neutrino interactions in the Earth and predicts the spectrum of the $τ$-lepton and muons that emerge. The $\texttt{nuPyProp}$ outputs are lookup tables of charged lepton exit probabilities and energies that can be used directly or as inputs to the $\texttt{nuSpaceSim}$ code designed to simulate optical and radio signals from extensive air showers induced by the emerging charged leptons. We describe the inputs to the code, demonstrate its flexibility and show selected results for $τ$-lepton and muon exit probabilities and energy distributions. The $\texttt{nuPyProp}$ code is open source, available on Github. △ Less

Submitted 13 February, 2023; v1 submitted 30 September, 2022; originally announced September 2022.

Comments: 42 pages, 21 figures, fixed typo in one of the author's name. Code available at https://github.com/NuSpaceSim/nupyprop

Journal ref: Journal of Cosmology and Astroparticle physics, Volume 2023, January 2023, 041

Recursive McCormick Linearization of Multilinear Programs

Authors: Arvind U Raghunathan, Carlos Cardonha, David Bergman, Carlos J Nohra

Abstract: Linear programming (LP) relaxations are widely employed in exact solution methods for multilinear programs (MLP). One example is the family of Recursive McCormick Linearization (RML) strategies, where bilinear products are substituted for artificial variables, which deliver a relaxation of the original problem when introduced together with concave and convex envelopes. In this article, we introduc… ▽ More Linear programming (LP) relaxations are widely employed in exact solution methods for multilinear programs (MLP). One example is the family of Recursive McCormick Linearization (RML) strategies, where bilinear products are substituted for artificial variables, which deliver a relaxation of the original problem when introduced together with concave and convex envelopes. In this article, we introduce the first systematic approach for identifying RMLs, in which we focus on the identification of linear relaxation with a small number of artificial variables and with strong LP bounds. We present a novel mechanism for representing all the possible RMLs, which we use to design an exact mixed-integer programming (MIP) formulation for the identification of minimum-size RMLs; we show that this problem is NP-hard in general, whereas a special case is fixed-parameter tractable. Moreover, we explore structural properties of our formulation to derive an exact MIP model that identifies RMLs of a given size with the best possible relaxation bound is optimal. Our numerical results on a collection of benchmarks indicate that our algorithms outperform the RML strategy implemented in state-of-the-art global optimization solvers. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Comments: 22 pages, 11 figures, Under Review

Constrained Shortest-Path Reformulations via Decision Diagrams for Structured Two-stage Optimization Problems

Authors: Leonardo Lozano, David Bergman, Andre A. Cire

Abstract: Many discrete optimization problems are amenable to constrained shortest-path reformulations in an extended network space, a technique that has been key in convexification, bound strengthening, and search. In this paper, we propose a constrained variant of these models for two challenging classes of discrete two-stage optimization problems, where traditional methods (e.g., dualize-and-combine) are… ▽ More Many discrete optimization problems are amenable to constrained shortest-path reformulations in an extended network space, a technique that has been key in convexification, bound strengthening, and search. In this paper, we propose a constrained variant of these models for two challenging classes of discrete two-stage optimization problems, where traditional methods (e.g., dualize-and-combine) are not applicable compared to their continuous counterparts. Specifically, we propose a framework that models problems as decision diagrams and introduces side constraints either as linear inequalities in the underlying polyhedral representation, or as state variables in shortest-path dynamic programming models. For our first structured class, we investigate two-stage problems with interdiction constraints. We show that such constraints can be formulated as indicator functions in the arcs of the diagram, providing an alternative single-level reformulation of the problem via a network-flow representation. Our second structured class is classical robust optimization, where we leverage the decision diagram network to iteratively identify label variables, akin to an L-shaped method. We evaluate these strategies on a competitive project selection problem and the robust traveling salesperson with time windows, observing considerable improvements in computational efficiency as compared to general methods in the respective areas. △ Less

Submitted 7 July, 2024; v1 submitted 26 June, 2022; originally announced June 2022.

Ultra-High-Energy Cosmic Rays: The Intersection of the Cosmic and Energy Frontiers

Authors: A. Coleman, J. Eser, E. Mayotte, F. Sarazin, F. G. Schröder, D. Soldin, T. M. Venters, R. Aloisio, J. Alvarez-Muñiz, R. Alves Batista, D. Bergman, M. Bertaina, L. Caccianiga, O. Deligny, H. P. Dembinski, P. B. Denton, A. di Matteo, N. Globus, J. Glombitza, G. Golup, A. Haungs, J. R. Hörandel, T. R. Jaffe, J. L. Kelley, J. F. Krizmanic , et al. (73 additional authors not shown)

Abstract: The present white paper is submitted as part of the "Snowmass" process to help inform the long-term plans of the United States Department of Energy and the National Science Foundation for high-energy physics. It summarizes the science questions driving the Ultra-High-Energy Cosmic-Ray (UHECR) community and provides recommendations on the strategy to answer them in the next two decades. The present white paper is submitted as part of the "Snowmass" process to help inform the long-term plans of the United States Department of Energy and the National Science Foundation for high-energy physics. It summarizes the science questions driving the Ultra-High-Energy Cosmic-Ray (UHECR) community and provides recommendations on the strategy to answer them in the next two decades. △ Less

Submitted 15 April, 2023; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: Prepared as a solicited white paper for the 2021 Snowmass process. To be published in the Journal of High Energy Astrophysics. v2: fixed typos in author list. v3: included all community feedback received by July 1st 2022 and added the list of endorsers. v4 is the post-reviewer preprint accepted to Astroparticle Physics vol. 149

Journal ref: Astroparticle Physics Volume 149, July 2023, 102819

First High-speed Video Camera Observations of a Lightning Flash Associated with a Downward Terrestrial Gamma-ray Flash

Authors: R. U. Abbasi, M. M. F. Saba, J. W. Belz, P. R. Krehbiel, W. Rison, N. Kieu, D. R. da Silva, Dan Rodeheffer, M. A. Stanley, J. Remington, J. Mazich, R. LeVon, K. Smout, A. Petrizze, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii , et al. (127 additional authors not shown)

Abstract: In this paper, we present the first high-speed video observation of a cloud-to-ground lightning flash and its associated downward-directed Terrestrial Gamma-ray Flash (TGF). The optical emission of the event was observed by a high-speed video camera running at 40,000 frames per second in conjunction with the Telescope Array Surface Detector, Lightning Mapping Array, interferometer, electric-field… ▽ More In this paper, we present the first high-speed video observation of a cloud-to-ground lightning flash and its associated downward-directed Terrestrial Gamma-ray Flash (TGF). The optical emission of the event was observed by a high-speed video camera running at 40,000 frames per second in conjunction with the Telescope Array Surface Detector, Lightning Mapping Array, interferometer, electric-field fast antenna, and the National Lightning Detection Network. The cloud-to-ground flash associated with the observed TGF was formed by a fast downward leader followed by a very intense return stroke peak current of -154 kA. The TGF occurred while the downward leader was below cloud base, and even when it was halfway in its propagation to ground. The suite of gamma-ray and lightning instruments, timing resolution, and source proximity offer us detailed information and therefore a unique look at the TGF phenomena. △ Less

Submitted 9 August, 2023; v1 submitted 10 May, 2022; originally announced May 2022.

Journal ref: Geophysical Research Letters, 50, e2023GL102958 (2023)

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

Authors: The ANTARES collaboration, A. Albert, S. Alves, M. André, M. Anghinolfi, M. Ardid, S. Ardid, J. -J. Aubert, J. Aublin, B. Baret, S. Basa, B. Belhorma, M. Bendahman, V. Bertin, S. Biagi, M. Bissinger, J. Boumaaza, M. Bouta, M. C. Bouwhuis, H. Brânzaş, R. Bruijn, J. Brunner, J. Busto, B. Caiffi, D. Calvo , et al. (1025 additional authors not shown)

Abstract: For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for corre… ▽ More For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs. △ Less

Submitted 23 August, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 39 pages, 7 figures, 4 tables; updated source files including xml authorlist

Report number: FERMILAB-PUB-22-033-AD-PPD-SCD-TD

Journal ref: ApJ 934 164 (2022)

Constraint Learning to Define Trust Regions in Predictive-Model Embedded Optimization

Authors: Chenbo Shi, Mohsen Emadikhiav, Leonardo Lozano, David Bergman

Abstract: There is a recent proliferation of research on the integration of machine learning and optimization. One expansive area within this research stream is predictive-model embedded optimization, which proposes the use of pre-trained predictive models as surrogates for uncertain or highly complex objective functions. In this setting, features of the predictive models become decision variables in the op… ▽ More There is a recent proliferation of research on the integration of machine learning and optimization. One expansive area within this research stream is predictive-model embedded optimization, which proposes the use of pre-trained predictive models as surrogates for uncertain or highly complex objective functions. In this setting, features of the predictive models become decision variables in the optimization problem. Despite a recent surge in publications in this area, only a few papers note the importance of incorporating trust region considerations in this decision-making pipeline, i.e., enforcing solutions to be similar to the data used to train the predictive models. Without such constraints, the evaluation of the predictive model at solutions obtained from optimization cannot be trusted and the practicality of the solutions may be unreasonable. In this paper, we provide an overview of the approaches appearing in the literature to construct a trust region, and propose three alternative approaches. Our numerical evaluation highlights that trust-region constraints learned through isolation forests, one of the newly proposed approaches, outperform all previously suggested approaches, both in terms of solution quality and computational time. △ Less

Submitted 19 October, 2022; v1 submitted 12 January, 2022; originally announced January 2022.

Optimizing over an ensemble of neural networks

Authors: Keliang Wang, Leonardo Lozano, Carlos Cardonha, David Bergman

Abstract: We study optimization problems where the objective function is modeled through feedforward neural networks with rectified linear unit (ReLU) activation. Recent literature has explored the use of a single neural network to model either uncertain or complex elements within an objective function. However, it is well known that ensembles of neural networks produce more stable predictions and have bett… ▽ More We study optimization problems where the objective function is modeled through feedforward neural networks with rectified linear unit (ReLU) activation. Recent literature has explored the use of a single neural network to model either uncertain or complex elements within an objective function. However, it is well known that ensembles of neural networks produce more stable predictions and have better generalizability than models with single neural networks, which motivates the investigation of ensembles of neural networks rather than single neural networks in decision-making pipelines. We study how to incorporate a neural network ensemble as the objective function of an optimization model and explore computational approaches for the ensuing problem. We present a mixed-integer linear program based on existing popular big-M formulations for optimizing over a single neural network. We develop a two-phase approach for our model that combines preprocessing procedures to tighten bounds for critical neurons in the neural networks with a Lagrangian relaxation-based branch-and-bound approach. Experimental evaluations of our solution methods suggest that using ensembles of neural networks yields more stable and higher quality solutions, compared to single neural networks, and that our optimization algorithm outperforms (the adaption of) a state-of-the-art approach in terms of computational time and optimality gaps. △ Less

Submitted 10 May, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: 29 pages, 9 tables, 4 figures

Optimizing the expected maximum of two linear functions defined on a multivariate Gaussian distribution

Authors: David Bergman, Carlos Cardonha, Jason Imbrogno, Leonardo Lozano

Abstract: We study stochastic optimization problems with objective function given by the expectation of the maximum of two linear functions defined on the component random variables of a multivariate Gaussian distribution. We consider random variables that are arbitrarily correlated, and we show that the problem is NP-hard even if the space of feasible solutions is unconstrained. We exploit a closed-form ex… ▽ More We study stochastic optimization problems with objective function given by the expectation of the maximum of two linear functions defined on the component random variables of a multivariate Gaussian distribution. We consider random variables that are arbitrarily correlated, and we show that the problem is NP-hard even if the space of feasible solutions is unconstrained. We exploit a closed-form expression for the objective function from the literature to construct a cutting-plane algorithm that can be seen as an extension of the integer L-shaped method for a highly nonlinear function, which includes the evaluation of the c.d.f and p.d.f of a standard normal random variable with decision variables as part of the arguments. To exhibit the model's applicability, we consider two featured applications. The first is daily fantasy sports, where the algorithm identifies entries with positive returns during the 2018-2019 National Football League season. The second is a special case of makespan minimization for two parallel machines and jobs with uncertain processing times; for the special case where the jobs are uncorrelated, we prove the equivalence between its deterministic and stochastic versions and show that our algorithm can deliver a constant-factor approximation guarantee for the problem. The results of our computational evaluation involving synthetic and real-world data suggest that our discretization and upper bounding techniques lead to significant computational improvements and that the proposed algorithm outperforms sub-optimal solutions approaches. △ Less

Submitted 13 December, 2021; originally announced December 2021.

Comments: 29 pages, 4 tables, 1 figure

Observation of Variations in Cosmic Ray Single Count Rates During Thunderstorms and Implications for Large-Scale Electric Field Changes

Authors: R. U. Abbasi, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, R. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, M. Hayashi , et al. (140 additional authors not shown)

Abstract: We present the first observation by the Telescope Array Surface Detector (TASD) of the effect of thunderstorms on the development of cosmic ray single count rate intensity over a 700 km$^{2}$ area. Observations of variations in the secondary low-energy cosmic ray counting rate, using the TASD, allow us to study the electric field inside thunderstorms, on a large scale, as it progresses on top of t… ▽ More We present the first observation by the Telescope Array Surface Detector (TASD) of the effect of thunderstorms on the development of cosmic ray single count rate intensity over a 700 km$^{2}$ area. Observations of variations in the secondary low-energy cosmic ray counting rate, using the TASD, allow us to study the electric field inside thunderstorms, on a large scale, as it progresses on top of the 700 km$^{2}$ detector, without dealing with the limitation of narrow exposure in time and space using balloons and aircraft detectors. In this work, variations in the cosmic ray intensity (single count rate) using the TASD, were studied and found to be on average at the $\sim(0.5-1)\%$ and up to 2\% level. These observations were found to be both in excess and in deficit. They were also found to be correlated with lightning in addition to thunderstorms. These variations lasted for tens of minutes; their footprint on the ground ranged from 6 to 24 km in diameter and moved in the same direction as the thunderstorm. With the use of simple electric field models inside the cloud and between cloud to ground, the observed variations in the cosmic ray single count rate were recreated using CORSIKA simulations. Depending on the electric field model used and the direction of the electric field in that model, the electric field magnitude that reproduces the observed low-energy cosmic ray single count rate variations was found to be approximately between 0.2-0.4 GV. This in turn allows us to get a reasonable insight on the electric field and its effect on cosmic ray air showers inside thunderstorms. △ Less

Submitted 18 November, 2021; originally announced November 2021.

Indications of a Cosmic Ray Source in the Perseus-Pisces Supercluster

Authors: Telescope Array Collaboration, R. U. Abbasi, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, R. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon , et al. (135 additional authors not shown)

Abstract: The Telescope Array Collaboration has observed an excess of events with $E \ge 10^{19.4} ~{\rm eV}$ in the data which is centered at (RA, dec) = ($19^\circ$, $35^\circ$). This is near the center of the Perseus-Pisces supercluster (PPSC). The PPSC is about $70 ~{\rm Mpc}$ distant and is the closest supercluster in the Northern Hemisphere (other than the Virgo supercluster of which we are a part). A… ▽ More The Telescope Array Collaboration has observed an excess of events with $E \ge 10^{19.4} ~{\rm eV}$ in the data which is centered at (RA, dec) = ($19^\circ$, $35^\circ$). This is near the center of the Perseus-Pisces supercluster (PPSC). The PPSC is about $70 ~{\rm Mpc}$ distant and is the closest supercluster in the Northern Hemisphere (other than the Virgo supercluster of which we are a part). A Li-Ma oversampling analysis with $20^\circ$-radius circles indicates an excess in the arrival direction of events with a local significance of about 4 standard deviations. The probability of having such excess close to the PPSC by chance is estimated to be 3.5 standard deviations. This result indicates that a cosmic ray source likely exists in that supercluster. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 8 pages, 4 figures, 1 table

Monte Carlo simulations of neutrino and charged lepton propagation in the Earth with nuPyProp

Authors: Sameer Patel, Mary Hall Reno, Yosui Akaike, Luis Anchordoqui, Douglas Bergman, Isaac Buckland, Austin Cummings, Johannes Eser, Claire Guépin, John F. Krizmanic, Simon Mackovjak, Angela Olinto, Thomas Paul, Alex Reustle, Andrew Romero-Wolf, Fred Sarazin, Tonia Venters, Lawrence Wiencke, Stephanie Wissel

Abstract: An accurate modeling of neutrino flux attenuation and the distribution of leptons they produce in transit through the Earth is an essential component to determine neutrino flux sensitivities of underground, sub-orbital and space-based detectors. Through neutrino oscillations over cosmic distances, astrophysical neutrino sources are expected to produce nearly equal fluxes of electron, muon and tau… ▽ More An accurate modeling of neutrino flux attenuation and the distribution of leptons they produce in transit through the Earth is an essential component to determine neutrino flux sensitivities of underground, sub-orbital and space-based detectors. Through neutrino oscillations over cosmic distances, astrophysical neutrino sources are expected to produce nearly equal fluxes of electron, muon and tau neutrinos. Of particular interest are tau neutrinos that interact in the Earth at modest slant depths to produce $τ$-leptons. Some $τ$-leptons emerge from the Earth and decay in the atmosphere to produce extensive air showers. Future balloon-borne and satellite-based optical Cherenkov neutrino telescopes will be sensitive to upward air showers from tau neutrino induced $τ$-lepton decays. We present nuPyProp, a python code that is part of the nuSpaceSim package. nuPyProp generates look-up tables for exit probabilities and energy distributions for $ν_τ\to τ$ and $ν_μ\to μ$ propagation in the Earth. This flexible code runs with either stochastic or continuous electromagnetic energy losses for the lepton transit through the Earth. Current neutrino cross section models and energy loss models are included along with templates for user input of other models. Results from nuPyProp are compared with other recent simulation packages for neutrino and charged lepton propagation. Sources of modeling uncertainties are described and quantified. △ Less

Submitted 16 September, 2021; originally announced September 2021.

Comments: 9 pages, 6 figures, PoS(ICRC2021)1203; for associated code, see https://github.com/NuSpaceSim/nupyprop

Surface detectors of the TAx4 experiment

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, R. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, N. Hayashida, K. Hibino , et al. (124 additional authors not shown)

Abstract: Telescope Array (TA) is the largest ultrahigh energy cosmic-ray (UHECR) observatory in the Northern Hemisphere. It explores the origin of UHECRs by measuring their energy spectrum, arrival-direction distribution, and mass composition using a surface detector (SD) array covering approximately 700 km$^2$ and fluorescence detector (FD) stations. TA has found evidence for a cluster of cosmic rays with… ▽ More Telescope Array (TA) is the largest ultrahigh energy cosmic-ray (UHECR) observatory in the Northern Hemisphere. It explores the origin of UHECRs by measuring their energy spectrum, arrival-direction distribution, and mass composition using a surface detector (SD) array covering approximately 700 km$^2$ and fluorescence detector (FD) stations. TA has found evidence for a cluster of cosmic rays with energies greater than 57 EeV. In order to confirm this evidence with more data, it is necessary to increase the data collection rate.We have begun building an expansion of TA that we call TAx4. In this paper, we explain the motivation, design, technical features, and expected performance of the TAx4 SD. We also present TAx4's current status and examples of the data that have already been collected. △ Less

Submitted 1 March, 2021; originally announced March 2021.

Comments: 26 pages, 17 figures, submitted to Nuclear Inst. and Methods in Physics Research, A

The POEMMA (Probe of Extreme Multi-Messenger Astrophysics) Observatory

Authors: A. V. Olinto, J. Krizmanic, J. H. Adams, R. Aloisio, L. A. Anchordoqui, A. Anzalone, M. Bagheri, D. Barghini, M. Battisti, D. R. Bergman, M. E. Bertaina, P. F. Bertone, F. Bisconti, M. Bustamante, F. Cafagna, R. Caruso, M. Casolino, K. Černý, M. J. Christl, A. L. Cummings, I. De Mitri, R. Diesing, R. Engel, J. Eser, K. Fang , et al. (51 additional authors not shown)

Abstract: The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to accurately observe ultra-high-energy cosmic rays (UHECRs) and cosmic neutrinos from space with sensitivity over the full celestial sky. POEMMA will observe the extensive air showers (EASs) from UHECRs and UHE neutrinos above 20 EeV via air fluorescence. Additionally, POEMMA will observe the Cherenkov signal from upward-movin… ▽ More The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to accurately observe ultra-high-energy cosmic rays (UHECRs) and cosmic neutrinos from space with sensitivity over the full celestial sky. POEMMA will observe the extensive air showers (EASs) from UHECRs and UHE neutrinos above 20 EeV via air fluorescence. Additionally, POEMMA will observe the Cherenkov signal from upward-moving EASs induced by Earth-interacting tau neutrinos above 20 PeV. The POEMMA spacecraft are designed to quickly re-orientate to follow up transient neutrino sources and obtain unparalleled neutrino flux sensitivity. Developed as a NASA Astrophysics Probe-class mission, POEMMA consists of two identical satellites flying in loose formation in 525 km altitude orbits. Each POEMMA instrument incorporates a wide field-of-view (45$^\circ$) Schmidt telescope with over 6 m$^2$ of collecting area. The hybrid focal surface of each telescope includes a fast (1~$μ$s) near-ultraviolet camera for EAS fluorescence observations and an ultrafast (10~ns) optical camera for Cherenkov EAS observations. In a 5-year mission, POEMMA will provide measurements that open new multi-messenger windows onto the most energetic events in the universe, enabling the study of new astrophysics and particle physics at these otherwise inaccessible energies. △ Less

Submitted 24 May, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Comments: 66 pages, 53 figures; JCAP accepted version

Journal ref: JCAP 06 (2021) 007

Observations of the Origin of Downward Terrestrial Gamma-Ray Flashes

Authors: J. W. Belz, P. R. Krehbiel, J. Remington, M. A. Stanley, R. U. Abbasi, R. LeVon, W. Rison, D. Rodeheffer, the Telescope Array Scientific Collaboration, :, T. Abu-Zayyad, M. Allen, E. Barcikowski, D. R. Bergman, S. A. Blake, M. Byrne, R. Cady, B. G. Cheon, M. Chikawa, A. di Matteo, T. Fujii, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich , et al. (116 additional authors not shown)

Abstract: In this paper we report the first close, high-resolution observations of downward-directed terrestrial gamma-ray flashes (TGFs) detected by the large-area Telescope Array cosmic ray observatory, obtained in conjunction with broadband VHF interferometer and fast electric field change measurements of the parent discharge. The results show that the TGFs occur during strong initial breakdown pulses (I… ▽ More In this paper we report the first close, high-resolution observations of downward-directed terrestrial gamma-ray flashes (TGFs) detected by the large-area Telescope Array cosmic ray observatory, obtained in conjunction with broadband VHF interferometer and fast electric field change measurements of the parent discharge. The results show that the TGFs occur during strong initial breakdown pulses (IBPs) in the first few milliseconds of negative cloud-to-ground and low-altitude intracloud flashes, and that the IBPs are produced by a newly-identified streamer-based discharge process called fast negative breakdown. The observations indicate the relativistic runaway electron avalanches (RREAs) responsible for producing the TGFs are initiated by embedded spark-like transient conducting events (TCEs) within the fast streamer system, and potentially also by individual fast streamers themselves. The TCEs are inferred to be the cause of impulsive sub-pulses that are characteristic features of classic IBP sferics. Additional development of the avalanches would be facilitated by the enhanced electric field ahead of the advancing front of the fast negative breakdown. In addition to showing the nature of IBPs and their enigmatic sub-pulses, the observations also provide a possible explanation for the unsolved question of how the streamer to leader transition occurs during the initial negative breakdown, namely as a result of strong currents flowing in the final stage of successive IBPs, extending backward through both the IBP itself and the negative streamer breakdown preceding the IBP. △ Less

Submitted 12 October, 2020; v1 submitted 29 September, 2020; originally announced September 2020.

Comments: Typo fixed and reference added. Manuscript is 36 pages. Supplemental Information is 42 pages. This paper is to be published in the Journal of Geophysical Research: Atmospheres. Online data repository: Open Science Framework DOI: 10.17605/OSF.IO/Z3XDA

Snowmass 2021 Letter of Interest: The Probe Of Multi-Messenger Astrophysics (POEMMA)

Authors: A. V. Olinto, F. Sarazin, J. H. Adams, R. Aloisio, L. A. Anchordoqui, M. Bagheri, D. Barghini, M. Battisti, D. R. Bergman, M. E. Bertaina, P. F. Bertone, F. Bisconti, M. Bustamante, M. Casolino, M. J. Christl, A. L. Cummings, I. De Mitri, R. Diesing, R. Engel, J. Eser, K. Fang, G. Fillipatos, F. Fenu, E. Gazda, C. Guepin , et al. (39 additional authors not shown)

Abstract: The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to identify the sources of Ultra-High-Energy Cosmic Rays (UHECRs) and to observe cosmic neutrinos, both with full-sky coverage. Developed as a NASA Astrophysics Probe-class mission, POEMMA consists of two spacecraft flying in a loose formation at 525 km altitude, 28.5 deg inclination orbits. Each spacecraft hosts a Schmidt tele… ▽ More The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to identify the sources of Ultra-High-Energy Cosmic Rays (UHECRs) and to observe cosmic neutrinos, both with full-sky coverage. Developed as a NASA Astrophysics Probe-class mission, POEMMA consists of two spacecraft flying in a loose formation at 525 km altitude, 28.5 deg inclination orbits. Each spacecraft hosts a Schmidt telescope with a large collecting area and wide field of view. A novel focal plane is optimized to observe both the UV fluorescence signal from extensive air showers (EASs) and the beamed optical Cherenkov signals from EASs. In POEMMA-stereo fluorescence mode, POEMMA will measure the spectrum, composition, and full-sky distribution of the UHECRs above 20 EeV with high statistics along with remarkable sensitivity to UHE neutrinos. The spacecraft are designed to quickly re-orient to a POEMMA-limb mode to observe neutrino emission from Target-of-Opportunity (ToO) transient astrophysical sources viewed just below the Earth's limb. In this mode, POEMMA will have unique sensitivity to cosmic neutrino tau events above 20 PeV by measuring the upward-moving EASs induced by the decay of the emerging tau leptons following the interactions of neutrino tau inside the Earth. △ Less

Submitted 1 September, 2020; v1 submitted 29 August, 2020; originally announced August 2020.

Comments: 3 pages, no figures, LoI for Snowmass 2021

Radio Detection of Ultra-high Energy Cosmic Rays with Low Lunar Orbiting SmallSats

Authors: Andrés Romero-Wolf, Jaime Alvarez-Muñiz, Luis A. Anchordoqui, Douglas Bergman, Washington Carvalho Jr., Austin L. Cummings, Peter Gorham, Casey J. Handmer, Nate Harvey, John Krizmanic, Kurtis Nishimura, Remy Prechelt, Mary Hall Reno, Harm Schoorlemmer, Gary Varner, Tonia Venters, Stephanie Wissel, Enrique Zas

Abstract: Ultra-high energy cosmic rays (UHECRs) are the most energetic particles observed and serve as a probe of the extreme universe. A key question to understanding the violent processes responsible for their acceleration is identifying which classes of astrophysical objects (active galactic nuclei or starburst galaxies, for example) correlate to their arrival directions. While source clustering is limi… ▽ More Ultra-high energy cosmic rays (UHECRs) are the most energetic particles observed and serve as a probe of the extreme universe. A key question to understanding the violent processes responsible for their acceleration is identifying which classes of astrophysical objects (active galactic nuclei or starburst galaxies, for example) correlate to their arrival directions. While source clustering is limited by deflections in the Galactic magnetic field, at the highest energies the scattering angles are sufficiently low to retain correlation with source catalogues. While there have been several studies attempting to identify source catalogue correlations with data from the Pierre Auger Observatory and the Telescope Array, the significance above an isotropic background has not yet reached the threshold for discovery. It has been known for several decades that a full-sky UHECR observatory would provide a substantial increase in sensitivity to the anisotropic component of UHECRs. There have been several concepts developed in that time targeting the identification of UHECR sources such as OWL, JEM-EUSO, and POEMMA, using fluorescence detection in the Earth's atmosphere from orbit. In this white paper, we present a concept called the Zettavolt Askaryan Polarimeter (ZAP), designed to identify the source of UHECRs using radio detection of the Askaryan radio emissions produced by UHECRs interacting in the Moon's regolith from low lunar orbit. △ Less

Submitted 25 August, 2020; originally announced August 2020.

Comments: 6 pages, 2 figures, letter of interest submitted to Snowmass 2021

Search for Large-scale Anisotropy on Arrival Directions of Ultra-high-energy Cosmic Rays Observed with the Telescope Array Experiment

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, N. Hayashida, K. Hibino , et al. (121 additional authors not shown)

Abstract: Motivated by the detection of a significant dipole structure in the arrival directions of ultrahigh-energy cosmic rays above 8 EeV reported by the Pierre Auger Observatory (Auger), we search for a large-scale anisotropy using data collected with the surface detector array of the Telescope Array Experiment (TA). With 11 years of TA data, a dipole structure in a projection of the right ascension is… ▽ More Motivated by the detection of a significant dipole structure in the arrival directions of ultrahigh-energy cosmic rays above 8 EeV reported by the Pierre Auger Observatory (Auger), we search for a large-scale anisotropy using data collected with the surface detector array of the Telescope Array Experiment (TA). With 11 years of TA data, a dipole structure in a projection of the right ascension is fitted with an amplitude of 3.3+- 1.9% and a phase of 131 +- 33 degrees. The corresponding 99% confidence-level upper limit on the amplitude is 7.3%. At the current level of statistics, the fitted result is compatible with both an isotropic distribution and the dipole structure reported by Auger. △ Less

Submitted 27 July, 2020; v1 submitted 30 June, 2020; originally announced July 2020.

Comments: 6 pages, 3 figures, 1 table, Proofed title. Added journal reference and DOI

Journal ref: The Astrophysical Journal Letters 898, L28 (2020)

Measurement of the Proton-Air Cross Section with Telescope Array's Black Rock Mesa and Long Ridge Fluorescence Detectors, and Surface Array in Hybrid Mode

Authors: R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, N. Hayashida, K. Hibino, R. Higuchi , et al. (120 additional authors not shown)

Abstract: Ultra high energy cosmic rays provide the highest known energy source in the universe to measure proton cross sections. Though conditions for collecting such data are less controlled than an accelerator environment, current generation cosmic ray observatories have large enough exposures to collect significant statistics for a reliable measurement for energies above what can be attained in the lab.… ▽ More Ultra high energy cosmic rays provide the highest known energy source in the universe to measure proton cross sections. Though conditions for collecting such data are less controlled than an accelerator environment, current generation cosmic ray observatories have large enough exposures to collect significant statistics for a reliable measurement for energies above what can be attained in the lab. Cosmic ray measurements of cross section use atmospheric calorimetry to measure depth of air shower maximum ($X_{\mathrm{max}}$), which is related to the primary particle's energy and mass. The tail of the $X_{\mathrm{max}}$ distribution is assumed to be dominated by showers generated by protons, allowing measurement of the inelastic proton-air cross section. In this work the proton-air inelastic cross section measurement, $σ^{\mathrm{inel}}_{\mathrm{p-air}}$, using data observed by Telescope Array's Black Rock Mesa and Long Ridge fluorescence detectors and surface detector array in hybrid mode is presented. $σ^{\mathrm{inel}}_{\mathrm{p-air}}$ is observed to be $520.1 \pm 35.8$[Stat.] $^{+25.0}_{-40}$[Sys.]~mb at $\sqrt{s} = 73$ TeV. The total proton-proton cross section is subsequently inferred from Glauber formalism and is found to be $σ^{\mathrm{tot}}_{\mathrm{pp}} = 139.4 ^{+23.4}_{-21.3}$ [Stat.]$ ^{+15.0}_{-24.0}$[Sys.]~mb. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Journal ref: Phys. Rev. D 102, 062004 (2020)

Evidence for a Supergalactic Structure of Magnetic Deflection Multiplets of Ultra-High Energy Cosmic Rays

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, N. Hayashida, K. Hibino , et al. (119 additional authors not shown)

Abstract: Evidence for a large-scale supergalactic cosmic ray multiplet (arrival directions correlated with energy) structure is reported for ultra-high energy cosmic ray (UHECR) energies above 10$^{19}$ eV using seven years of data from the Telescope Array (TA) surface detector and updated to 10 years. Previous energy-position correlation studies have made assumptions regarding magnetic field shapes and st… ▽ More Evidence for a large-scale supergalactic cosmic ray multiplet (arrival directions correlated with energy) structure is reported for ultra-high energy cosmic ray (UHECR) energies above 10$^{19}$ eV using seven years of data from the Telescope Array (TA) surface detector and updated to 10 years. Previous energy-position correlation studies have made assumptions regarding magnetic field shapes and strength, and UHECR composition. Here the assumption tested is that, since the supergalactic plane is a fit to the average matter density of the local Large Scale Structure (LSS), UHECR sources and intervening extragalactic magnetic fields are correlated with this plane. This supergalactic deflection hypothesis is tested by the entire field-of-view (FOV) behavior of the strength of intermediate-scale energy-angle correlations. These multiplets are measured in spherical cap section bins (wedges) of the FOV to account for coherent and random magnetic fields. The structure found is consistent with supergalactic deflection, the previously published energy spectrum anisotropy results of TA (the hotspot and coldspot), and toy-model simulations of a supergalactic magnetic sheet. The seven year data post-trial significance of this supergalactic structure of multiplets appearing by chance, on an isotropic sky, is found by Monte Carlo simulation to be 4.2$σ$. The ten years of data post-trial significance is 4.1$σ$. Furthermore, the starburst galaxy M82 is shown to be a possible source of the TA Hotspot, and an estimate of the supergalactic magnetic field using UHECR measurements is presented. △ Less

Submitted 2 July, 2020; v1 submitted 14 May, 2020; originally announced May 2020.

White Paper: ARIANNA-200 high energy neutrino telescope

Authors: A. Anker, P. Baldi, S. W. Barwick, D. Bergman, H. Bernhoff, D. Z. Besson, N. Bingefors, O. Botner, P. Chen, Y. Chen, D. García-Fernández, G. Gaswint, C. Glaser, A. Hallgren, J. C. Hanson, J. J. Huang, S. R. Klein, S. A. Kleinfelder, C. -Y. Kuo, R. Lahmann, U. Latif, T. Liu, Y. Lyu, S. McAleer, J. Nam , et al. (11 additional authors not shown)

Abstract: The proposed ARIANNA-200 neutrino detector, located at sea-level on the Ross Ice Shelf, Antarctica, consists of 200 autonomous and independent detector stations separated by 1 kilometer in a uniform triangular mesh, and serves as a pathfinder mission for the future IceCube-Gen2 project. The primary science mission of ARIANNA-200 is to search for sources of neutrinos with energies greater than 10^1… ▽ More The proposed ARIANNA-200 neutrino detector, located at sea-level on the Ross Ice Shelf, Antarctica, consists of 200 autonomous and independent detector stations separated by 1 kilometer in a uniform triangular mesh, and serves as a pathfinder mission for the future IceCube-Gen2 project. The primary science mission of ARIANNA-200 is to search for sources of neutrinos with energies greater than 10^17 eV, complementing the reach of IceCube. An ARIANNA observation of a neutrino source would provide strong insight into the enigmatic sources of cosmic rays. ARIANNA observes the radio emission from high energy neutrino interactions in the Antarctic ice. Among radio based concepts under current investigation, ARIANNA-200 would uniquely survey the vast majority of the southern sky at any instant in time, and an important region of the northern sky, by virtue of its location on the surface of the Ross Ice Shelf in Antarctica. The broad sky coverage is specific to the Moore's Bay site, and makes ARIANNA-200 ideally suited to contribute to the multi-messenger thrust by the US National Science Foundation, Windows on the Universe - Multi-Messenger Astrophysics, providing capabilities to observe explosive sources from unknown directions. The ARIANNA architecture is designed to measure the angular direction to within 3 degrees for every neutrino candidate, which too plays an important role in the pursuit of multi-messenger observations of astrophysical sources. △ Less

Submitted 21 April, 2020; originally announced April 2020.

JANOS: An Integrated Predictive and Prescriptive Modeling Framework

Authors: David Bergman, Teng Huang, Philip Brooks, Andrea Lodi, Arvind U. Raghunathan

Abstract: Business research practice is witnessing a surge in the integration of predictive modeling and prescriptive analysis. We describe a modeling framework JANOS that seamlessly integrates the two streams of analytics, for the first time allowing researchers and practitioners to embed machine learning models in an optimization framework. JANOS allows for specifying a prescriptive model using standard o… ▽ More Business research practice is witnessing a surge in the integration of predictive modeling and prescriptive analysis. We describe a modeling framework JANOS that seamlessly integrates the two streams of analytics, for the first time allowing researchers and practitioners to embed machine learning models in an optimization framework. JANOS allows for specifying a prescriptive model using standard optimization modeling elements such as constraints and variables. The key novelty lies in providing modeling constructs that allow for the specification of commonly used predictive models and their features as constraints and variables in the optimization model. The framework considers two sets of decision variables; regular and predicted. The relationship between the regular and the predicted variables are specified by the user as pre-trained predictive models. JANOS currently supports linear regression, logistic regression, and neural network with rectified linear activation functions, but we plan to expand on this set in the future. In this paper, we demonstrate the flexibility of the framework through an example on scholarship allocation in a student enrollment problem and provide a numeric performance evaluation. △ Less

Submitted 21 November, 2019; originally announced November 2019.

Comments: 12 pages, 3 figures

Template-based Minor Embedding for Adiabatic Quantum Optimization

Authors: Thiago Serra, Teng Huang, Arvind Raghunathan, David Bergman

Abstract: Quantum Annealing (QA) can be used to quickly obtain near-optimal solutions for Quadratic Unconstrained Binary Optimization (QUBO) problems. In QA hardware, each decision variable of a QUBO should be mapped to one or more adjacent qubits in such a way that pairs of variables defining a quadratic term in the objective function are mapped to some pair of adjacent qubits. However, qubits have limited… ▽ More Quantum Annealing (QA) can be used to quickly obtain near-optimal solutions for Quadratic Unconstrained Binary Optimization (QUBO) problems. In QA hardware, each decision variable of a QUBO should be mapped to one or more adjacent qubits in such a way that pairs of variables defining a quadratic term in the objective function are mapped to some pair of adjacent qubits. However, qubits have limited connectivity in existing QA hardware. This has spurred work on preprocessing algorithms for embedding the graph representing problem variables with quadratic terms into the hardware graph representing qubits adjacencies, such as the Chimera graph in hardware produced by D-Wave Systems. In this paper, we use integer linear programming to search for an embedding of the problem graph into certain classes of minors of the Chimera graph, which we call template embeddings. One of these classes corresponds to complete bipartite graphs, for which we show the limitation of the existing approach based on minimum Odd Cycle Transversals (OCTs). One of the formulations presented is exact, and thus can be used to certify the absence of a minor embedding using that template. On an extensive test set consisting of random graphs from five different classes of varying size and sparsity, we can embed more graphs than a state-of-the-art OCT-based approach, our approach scales better with the hardware size, and the runtime is generally orders of magnitude smaller. △ Less

Submitted 19 January, 2021; v1 submitted 4 October, 2019; originally announced October 2019.

Comments: INFORMS Journal on Computing (to appear)

The POEMMA (Probe of Extreme Multi-Messenger Astrophysics) mission

Authors: A. V. Olinto, J. H. Adams, R. Aloisio, L. A. Anchordoqui, D. R. Bergman, M. E. Bertaina, P. Bertone, F. Bisconti, M. Bustamante, M. Casolino, M. J. Christl, A. L. Cummings, I. De Mitri, R. Diesing, J. B. Eser, F. Fenu, C. Guépin, E. A. Hays, E. Judd, J. F. Krizmanic, E. Kuznetsov, A. Liberatore, S. Mackovjak, J. McEnery, J. W. Mitchell , et al. (20 additional authors not shown)

Abstract: The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to observe cosmic neutrinos (CNs) above 20 PeV and ultra-high energy cosmic rays (UHECRs) above 20 EeV over the full sky. The POEMMA mission calls for two identical satellites flying in loose formation, each comprised of a 4-meter wide field-of-view (45 degrees) Schmidt photometer. The hybrid focal surface includes a fast (1… ▽ More The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to observe cosmic neutrinos (CNs) above 20 PeV and ultra-high energy cosmic rays (UHECRs) above 20 EeV over the full sky. The POEMMA mission calls for two identical satellites flying in loose formation, each comprised of a 4-meter wide field-of-view (45 degrees) Schmidt photometer. The hybrid focal surface includes a fast (1 $μ$s) ultraviolet camera for fluorescence observations and an ultrafast (10 ns) optical camera for Cherenkov observations. POEMMA will provide new multi-messenger windows onto the most energetic events in the universe, enabling the study of new astrophysics and particle physics at these otherwise inaccessible energies. △ Less

Submitted 18 September, 2019; originally announced September 2019.

Comments: 8 pages, 6 figures, presented in 36th International Cosmic Ray Conference (Madison). arXiv admin note: substantial text overlap with arXiv:1907.06217

Journal ref: Proceedings of the 36th International Cosmic Ray Conference (Madison; 2019) Pos(ICRC2019) 378

POEMMA (Probe of Extreme Multi-Messenger Astrophysics) design

Authors: A. V. Olinto, J. H. Adams, R. Aloisio, L. A. Anchordoqui, D. R. Bergman, M. E. Bertaina, P. Bertone, F. Bisconti, M. Bustamante, M. Casolino, M. J. Christl, A. L. Cummings, I. De Mitri, R. Diesing, J. Eser, F. Fenu, C. Guepin, E. A. Hays, E. G. Judd, J. F. Krizmanic, E. Kuznetsov, A. Liberatore, S. Mackovjak, J. McEnery, J. W. Mitchell , et al. (20 additional authors not shown)

Abstract: The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is a NASA Astrophysics probe-class mission designed to observe ultra-high energy cosmic rays (UHECRs) and cosmic neutrinos from space. Astro2020 APC white paper: Medium-class Space Particle Astrophysics Project. The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is a NASA Astrophysics probe-class mission designed to observe ultra-high energy cosmic rays (UHECRs) and cosmic neutrinos from space. Astro2020 APC white paper: Medium-class Space Particle Astrophysics Project. △ Less

Submitted 14 July, 2019; originally announced July 2019.

Comments: Astro2020 APC white paper: Medium-class Space Particle Astrophysics Project

Performance and science reach of POEMMA for ultrahigh-energy particles

Authors: Luis A. Anchordoqui, Douglas R. Bergman, Mario E. Bertaina, Francesco Fenu, John F. Krizmanic, Alessandro Liberatore, Angela V. Olinto, Mary Hall Reno, Fred Sarazin, Kenji Shinozaki, Jorge F. Soriano, Ralf Ulrich, Michael Unger, Tonia M. Venters, Lawrence Wiencke

Abstract: The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is a potential NASA Astrophysics Probe-class mission designed to observe ultra-high energy cosmic rays (UHECRs) and cosmic neutrinos from space. POEMMA will monitor colossal volumes of the Earth's atmosphere to detect extensive air showers (EASs) produced by extremely energetic cosmic messengers: UHECRs above 20 EeV over the full sky and c… ▽ More The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is a potential NASA Astrophysics Probe-class mission designed to observe ultra-high energy cosmic rays (UHECRs) and cosmic neutrinos from space. POEMMA will monitor colossal volumes of the Earth's atmosphere to detect extensive air showers (EASs) produced by extremely energetic cosmic messengers: UHECRs above 20 EeV over the full sky and cosmic neutrinos above 20 PeV. We focus most of this study on the impact of POEMMA for UHECR science by simulating the detector response and mission performance for EAS from UHECRs. We show that POEMMA will provide a significant increase in the statistics of observed UHECRs at the highest energies over the entire sky. POEMMA will be the first UHECR fluorescence detector deployed in space that will provide high-quality stereoscopic observations of the longitudinal development of air showers. Therefore, it will be able to provide event-by-event estimates of the calorimetric energy and nuclear mass of UHECRs. The particle physics in the interactions limits the interpretation of the shower maximum on an event by event basis. In contrast, the calorimetric energy measurement is significantly less sensitive to the different possible final states in the early interactions. We study the prospects to discover the origin and nature of UHECRs using expectations for measurements of the energy spectrum, the distribution of arrival direction, and the atmospheric column depth at which the EAS longitudinal development reaches maximum. We also explore supplementary science capabilities of POEMMA through its sensitivity to particle interactions at extreme energies and its ability to detect ultra-high energy neutrinos and photons produced by top-down models including cosmic strings and super-heavy dark matter particle decay in the halo of the Milky Way. △ Less

Submitted 25 October, 2019; v1 submitted 8 July, 2019; originally announced July 2019.

Comments: 40 pages revtex, with 42 figures

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

Search for Ultra-High-Energy Neutrinos with the Telescope Array Surface Detector

Authors: R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino, K. Honda , et al. (112 additional authors not shown)

Abstract: We present an upper limit on the flux of ultra-high-energy down-going neutrinos for $E > 10^{18}\ \mbox{eV}$ derived with the nine years of data collected by the Telescope Array surface detector (05-11-2008 -- 05-10-2017). The method is based on the multivariate analysis technique, so-called Boosted Decision Trees (BDT). Proton-neutrino classifier is built upon 16 observables related to both the p… ▽ More We present an upper limit on the flux of ultra-high-energy down-going neutrinos for $E > 10^{18}\ \mbox{eV}$ derived with the nine years of data collected by the Telescope Array surface detector (05-11-2008 -- 05-10-2017). The method is based on the multivariate analysis technique, so-called Boosted Decision Trees (BDT). Proton-neutrino classifier is built upon 16 observables related to both the properties of the shower front and the lateral distribution function. △ Less

Submitted 12 May, 2020; v1 submitted 9 May, 2019; originally announced May 2019.

Comments: 10 pages, 4 figures, accepted to JETP

Journal ref: JETP Vol. 158 (8(2)) (2020)

Search for point sources of ultra-high energy photons with the Telescope Array surface detector

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. diMatteo, T. Fujii, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino , et al. (114 additional authors not shown)

Abstract: The surface detector (SD) of the Telescope Array (TA) experiment allows one to indirectly detect photons with energies of order $10^{18}$ eV and higher and to separate photons from the cosmic-ray background. In this paper we present the results of a blind search for point sources of ultra-high energy (UHE) photons in the Northern sky using the TA SD data. The photon-induced extensive air showers (… ▽ More The surface detector (SD) of the Telescope Array (TA) experiment allows one to indirectly detect photons with energies of order $10^{18}$ eV and higher and to separate photons from the cosmic-ray background. In this paper we present the results of a blind search for point sources of ultra-high energy (UHE) photons in the Northern sky using the TA SD data. The photon-induced extensive air showers (EAS) are separated from the hadron-induced EAS background by means of a multivariate classifier based upon 16 parameters that characterize the air shower events. No significant evidence for the photon point sources is found. The upper limits are set on the flux of photons from each particular direction in the sky within the TA field of view, according to the experiment's angular resolution for photons. Average 95% C.L. upper limits for the point-source flux of photons with energies greater than $10^{18}$, $10^{18.5}$, $10^{19}$, $10^{19.5}$ and $10^{20}$ eV are $0.094$, $0.029$, $0.010$, $0.0073$ and $0.0058$ km$^{-2}$yr$^{-1}$ respectively. For the energies higher than $10^{18.5}$ eV, the photon point-source limits are set for the first time. Numerical results for each given direction in each energy range are provided as a supplement to this paper. △ Less

Submitted 9 March, 2020; v1 submitted 30 March, 2019; originally announced April 2019.

Comments: accepted to MNRAS, 11 pages, 4 figures, 2 tables; results in text-file format are supplemented to paper source

Report number: INR-TH-2019-005

Journal ref: MNRAS 492 (2020), 3984

High-Energy Galactic Cosmic Rays (Astro2020 Science White Paper)

Authors: Frank G. Schröder, Tareq AbuZayyad, Luis Anchordoqui, Karen Andeen, Xinhua Bai, Segev BenZvi, Doug Bergman, Alan Coleman, Hans Dembinski, Michael DuVernois, Tom Gaisser, Francis Halzen, Andreas Haungs, John Kelley, Hermann Kolanoski, Frank McNally, Markus Roth, Frederic Sarazin, Dave Seckel, Radomir Smida, Dennis Soldin, Delia Tosi

Abstract: The origin of the highest energy Galactic cosmic rays is still not understood, nor is the transition to EeV extragalactic particles. Scientific progress requires enhancements of existing air-shower arrays, such as: IceCube with its surface detector IceTop, and the low-energy extensions of both the Telescope Array and the Pierre Auger Observatory. The origin of the highest energy Galactic cosmic rays is still not understood, nor is the transition to EeV extragalactic particles. Scientific progress requires enhancements of existing air-shower arrays, such as: IceCube with its surface detector IceTop, and the low-energy extensions of both the Telescope Array and the Pierre Auger Observatory. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: Science White Paper for the Astronomy and Astrophysics Decadal Survey (Astro2020)

Journal ref: 2019 BAAS 51(3) 131

What is the nature and origin of the highest-energy particles in the universe?

Authors: Fred Sarazin, Luis Anchordoqui, James Beatty, Douglas Bergman, Corbin Covault, Glennys Farrar, John Krizmanic, David Nitz, Angela Olinto, Michael Unger, Peter Tinyakov, Lawrence Wiencke

Abstract: This white paper was submitted to the US Astronomy and Astrophysics Decadal Survey (Astro2020) and defines the science questions to be answered in the next decade in the field of Ultra-High Energy Cosmic-Rays. Following a review of the recent experimental and theoretical advances in the field, the paper outlines strategies and requirements desirable for the design of future experiments. This white paper was submitted to the US Astronomy and Astrophysics Decadal Survey (Astro2020) and defines the science questions to be answered in the next decade in the field of Ultra-High Energy Cosmic-Rays. Following a review of the recent experimental and theoretical advances in the field, the paper outlines strategies and requirements desirable for the design of future experiments. △ Less

Submitted 23 April, 2019; v1 submitted 10 March, 2019; originally announced March 2019.

Comments: 5 pages, 3 figures

Binary Decision Diagrams for Bin Packing with Minimum Color Fragmentation

Authors: David Bergman, Carlos Cardonha, Saharnaz Mehrani

Abstract: Bin Packing with Minimum Color Fragmentation (BPMCF) is an extension of the Bin Packing Problem in which each item has a size and a color and the goal is to minimize the sum of the number of bins containing items of each color. In this work, we introduce BPMCF and present a decomposition strategy to solve the problem, where the assignment of items to bins is formulated as a binary decision diagram… ▽ More Bin Packing with Minimum Color Fragmentation (BPMCF) is an extension of the Bin Packing Problem in which each item has a size and a color and the goal is to minimize the sum of the number of bins containing items of each color. In this work, we introduce BPMCF and present a decomposition strategy to solve the problem, where the assignment of items to bins is formulated as a binary decision diagram and an optimal integrated solutions is identified through a mixed-integer linear programming model. Our computational experiments show that the proposed approach greatly outperforms a direct formulation of BPMCF and that its performance is suitable for large instances of the problem. △ Less

Submitted 30 November, 2018; originally announced December 2018.

Comments: 10 pages, 2 figures, 1 table, currently under review

MSC Class: 90B80

Symmetry constrained machine learning

Authors: Doron L. Bergman

Abstract: Symmetry, a central concept in understanding the laws of nature, has been used for centuries in physics, mathematics, and chemistry, to help make mathematical models tractable. Yet, despite its power, symmetry has not been used extensively in machine learning, until rather recently. In this article we show a general way to incorporate symmetries into machine learning models. We demonstrate this wi… ▽ More Symmetry, a central concept in understanding the laws of nature, has been used for centuries in physics, mathematics, and chemistry, to help make mathematical models tractable. Yet, despite its power, symmetry has not been used extensively in machine learning, until rather recently. In this article we show a general way to incorporate symmetries into machine learning models. We demonstrate this with a detailed analysis on a rather simple real world machine learning system - a neural network for classifying handwritten digits, lacking bias terms for every neuron. We demonstrate that ignoring symmetries can have dire over-fitting consequences, and that incorporating symmetry into the model reduces over-fitting, while at the same time reducing complexity, ultimately requiring less training data, and taking less time and resources to train. △ Less

Submitted 10 September, 2019; v1 submitted 16 November, 2018; originally announced November 2018.

Journal ref: In: Bi Y., Bhatia R., Kapoor S. (eds) Intelligent Systems and Applications. IntelliSys 2019. Advances in Intelligent Systems and Computing, vol 1038. Springer, Cham

Constraints on the diffuse photon flux with energies above $10^{18}$ eV using the surface detector of the Telescope Array experiment

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujita, M. Fukushima, G. Furlich, T. Goto, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino , et al. (118 additional authors not shown)

Abstract: We present the results of the search for ultra-high-energy photons with nine years of data from the Telescope Array surface detector. A multivariate classifier is built upon 16 reconstructed parameters of the extensive air shower. These parameters are related to the curvature and the width of the shower front, the steepness of the lateral distribution function, and the timing parameters of the wav… ▽ More We present the results of the search for ultra-high-energy photons with nine years of data from the Telescope Array surface detector. A multivariate classifier is built upon 16 reconstructed parameters of the extensive air shower. These parameters are related to the curvature and the width of the shower front, the steepness of the lateral distribution function, and the timing parameters of the waveforms sensitive to the shower muon content. A total number of two photon candidates found in the search is fully compatible with the expected background. The $95\%\,$CL limits on the diffuse flux of the photons with energies greater than $10^{18.0}$, $10^{18.5}$, $10^{19.0}$, $10^{19.5}$ and $10^{20.0}$ eV are set at the level of $0.067$, $0.012$, $0.0036$, $0.0013$, $0.0013~\mbox{km}^{-2}\mbox{yr}^{-1}\mbox{sr}^{-1}$ correspondingly. △ Less

Submitted 19 March, 2019; v1 submitted 9 November, 2018; originally announced November 2018.

Comments: 7 pages, 2 figures

Report number: INR-TH-2018-026

Journal ref: Astroparticle Physics, 110 (2019) 8-14

Improving Optimization Bounds using Machine Learning: Decision Diagrams meet Deep Reinforcement Learning

Authors: Quentin Cappart, Emmanuel Goutierre, David Bergman, Louis-Martin Rousseau

Abstract: Finding tight bounds on the optimal solution is a critical element of practical solution methods for discrete optimization problems. In the last decade, decision diagrams (DDs) have brought a new perspective on obtaining upper and lower bounds that can be significantly better than classical bounding mechanisms, such as linear relaxations. It is well known that the quality of the bounds achieved th… ▽ More Finding tight bounds on the optimal solution is a critical element of practical solution methods for discrete optimization problems. In the last decade, decision diagrams (DDs) have brought a new perspective on obtaining upper and lower bounds that can be significantly better than classical bounding mechanisms, such as linear relaxations. It is well known that the quality of the bounds achieved through this flexible bounding method is highly reliant on the ordering of variables chosen for building the diagram, and finding an ordering that optimizes standard metrics is an NP-hard problem. In this paper, we propose an innovative and generic approach based on deep reinforcement learning for obtaining an ordering for tightening the bounds obtained with relaxed and restricted DDs. We apply the approach to both the Maximum Independent Set Problem and the Maximum Cut Problem. Experimental results on synthetic instances show that the deep reinforcement learning approach, by achieving tighter objective function bounds, generally outperforms ordering methods commonly used in the literature when the distribution of instances is known. To the best knowledge of the authors, this is the first paper to apply machine learning to directly improve relaxation bounds obtained by general-purpose bounding mechanisms for combinatorial optimization problems. △ Less

Submitted 27 February, 2019; v1 submitted 10 September, 2018; originally announced September 2018.

Comments: Accepted and presented at AAAI'19

Testing a reported correlation between arrival directions of ultrahigh-energy cosmic rays and a flux pattern from nearby starburst galaxies using Telescope Array data

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujita, M. Fukushima, G. Furlich, T. Goto, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino , et al. (117 additional authors not shown)

Abstract: The Pierre Auger Collaboration (Auger) recently reported a correlation between the arrival directions of cosmic rays with energies above 39 EeV and the flux pattern of 23 nearby starburst galaxies (SBGs). In this Letter, we tested the same hypothesis using cosmic rays detected by the Telescope Array experiment (TA) in the 9-year period from May 2008 to May 2017. Unlike the Auger analysis, we did n… ▽ More The Pierre Auger Collaboration (Auger) recently reported a correlation between the arrival directions of cosmic rays with energies above 39 EeV and the flux pattern of 23 nearby starburst galaxies (SBGs). In this Letter, we tested the same hypothesis using cosmic rays detected by the Telescope Array experiment (TA) in the 9-year period from May 2008 to May 2017. Unlike the Auger analysis, we did not optimize the parameter values but kept them fixed to the best-fit values found by Auger, namely 9.7% for the anisotropic fraction of cosmic rays assumed to originate from the SBGs in the list and 12.9° for the angular scale of the correlations. The energy threshold we adopted is 43 EeV, corresponding to 39 EeV in Auger when taking into account the energy-scale difference between two experiments. We find that the TA data is compatible with isotropy to within 1.1σ and with the Auger result to within 1.4σ, meaning that it is not capable to discriminate between these two hypotheses. △ Less

Submitted 22 October, 2018; v1 submitted 5 September, 2018; originally announced September 2018.

Comments: More accurate title, fixed an inconsistency in the declination range used, with a few figures slightly changing as a result, other minor revisions

Journal ref: R. U. Abbasi et al 2018 ApJL 867 L27

Mass composition of ultra-high-energy cosmic rays with the Telescope Array Surface Detector Data

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. di Matteo, T. Fujii, K. Fujita, M. Fukushima, G. Furlich, T. Goto, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino , et al. (118 additional authors not shown)

Abstract: The results on ultra-high-energy cosmic rays (UHECR) mass composition obtained with the Telescope Array surface detector are presented. The analysis employs the boosted decision tree (BDT) multivariate analysis built upon 14 observables related to both the properties of the shower front and the lateral distribution function. The multivariate classifier is trained with Monte-Carlo sets of events in… ▽ More The results on ultra-high-energy cosmic rays (UHECR) mass composition obtained with the Telescope Array surface detector are presented. The analysis employs the boosted decision tree (BDT) multivariate analysis built upon 14 observables related to both the properties of the shower front and the lateral distribution function. The multivariate classifier is trained with Monte-Carlo sets of events induced by the primary protons and iron. An average atomic mass of UHECR is presented for energies $10^{18.0}-10^{20.0}\ \mbox{eV}$. The average atomic mass of primary particles shows no significant energy dependence and corresponds to $\langle \ln A \rangle = 2.0 \pm 0.1 (stat.) \pm 0.44 (syst.)$. The result is compared to the mass composition obtained by the Telescope Array with $\mbox{X}_{\mbox{max}}$ technique along with the results of other experiments. Possible systematic errors of the method are discussed. △ Less

Submitted 24 January, 2019; v1 submitted 10 August, 2018; originally announced August 2018.

Comments: 11 pages, 9 figures

Journal ref: Phys. Rev. D 99, 022002 (2019)

Seamless Multimodal Transportation Scheduling

Authors: Arvind U Raghunathan, David Bergman, John Hooker, Thiago Serra, Shingo Kobori

Abstract: Ride-hailing services have expanded the role of shared mobility in passenger transportation systems, creating new markets and creative planning solutions for major urban centers. In this paper, we consider their use for the first-mile or last-mile passenger transportation in coordination with a mass transit service to provide a seamless multimodal transportation experience for the user. A system t… ▽ More Ride-hailing services have expanded the role of shared mobility in passenger transportation systems, creating new markets and creative planning solutions for major urban centers. In this paper, we consider their use for the first-mile or last-mile passenger transportation in coordination with a mass transit service to provide a seamless multimodal transportation experience for the user. A system that provides passengers with predictable information on travel and waiting times in their commutes is immensely valuable. We envision that the passengers will inform the system of their desired travel and arrival windows so that the system can jointly optimize the schedules of passengers. The problem we study balances minimizing travel time and the number of trips taken by the last-mile vehicles, so that long-term planning, maintenance, and environmental impact are all taken into account. We focus on the case where the last-mile service aggregates passengers by destination. We show that this problem is NP-hard, and propose a decision diagram-based branch-and-price decomposition model that can solve instances of real-world size (10,000 passengers spread over an hour, 50 last-mile destinations, 600 last-mile vehicles) in computational time (~1 minute) that is orders-of-magnitude faster than other methods appearing in the literature. Our experiments also indicate that aggregating passengers by destination on the last-mile service provides high-quality solutions to more general settings. △ Less

Submitted 27 March, 2022; v1 submitted 25 July, 2018; originally announced July 2018.

Comments: 41 pages

MSC Class: 90B06; 90B10; 90B50; 90C06; 90C90

Study of muons from ultra-high energy cosmic ray air showers measured with the Telescope Array experiment

Authors: Telescope Array Collaboration, R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. Di Matteo, T. Fujii, K. Fujita, M. Fukushima, G. Furlich, T. Goto, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino , et al. (117 additional authors not shown)

Abstract: One of the uncertainties in interpretation of ultra-high energy cosmic ray (UHECR) data comes from the hadronic interaction models used for air shower Monte Carlo (MC) simulations. The number of muons observed at the ground from UHECR-induced air showers is expected to depend upon the hadronic interaction model. One may therefore test the hadronic interaction models by comparing the measured numbe… ▽ More One of the uncertainties in interpretation of ultra-high energy cosmic ray (UHECR) data comes from the hadronic interaction models used for air shower Monte Carlo (MC) simulations. The number of muons observed at the ground from UHECR-induced air showers is expected to depend upon the hadronic interaction model. One may therefore test the hadronic interaction models by comparing the measured number of muons with the MC prediction. In this paper, we present the results of studies of muon densities in UHE extensive air showers obtained by analyzing the signal of surface detector stations which should have high $\it{muon \, purity}$. The muon purity of a station will depend on both the inclination of the shower and the relative position of the station. In 7 years' data from the Telescope Array experiment, we find that the number of particles observed for signals with an expected muon purity of $\sim$65% at a lateral distance of 2000 m from the shower core is $1.72 \pm 0.10{\rm (stat.)} \pm 0.37 {\rm (syst.)}$ times larger than the MC prediction value using the QGSJET II-03 model for proton-induced showers. A similar effect is also seen in comparisons with other hadronic models such as QGSJET II-04, which shows a $1.67 \pm 0.10 \pm 0.36$ excess. We also studied the dependence of these excesses on lateral distances and found a slower decrease of the lateral distribution of muons in the data as compared to the MC, causing larger discrepancy at larger lateral distances. △ Less

Submitted 11 April, 2018; originally announced April 2018.

Journal ref: Phys. Rev. D 98, 022002 (2018)

Predictive and Prescriptive Analytics for Location Selection of Add-on Retail Products

Authors: Teng Huang, David Bergman, Ram Gopal

Abstract: In this paper, we study an analytical approach to selecting expansion locations for retailers selling add-on products whose demand is derived from the demand of another base product. Demand for the add-on product is realized only as a supplement to the demand of the base product. In our context, either of the two products could be subject to spatial autocorrelation where demand at a given location… ▽ More In this paper, we study an analytical approach to selecting expansion locations for retailers selling add-on products whose demand is derived from the demand of another base product. Demand for the add-on product is realized only as a supplement to the demand of the base product. In our context, either of the two products could be subject to spatial autocorrelation where demand at a given location is impacted by demand at other locations. Using data from an industrial partner selling add-on products, we build predictive models for understanding the derived demand of the add-on product and establish an optimization framework for automating expansion decisions to maximize expected sales. Interestingly, spatial autocorrelation and the complexity of the predictive model impact the complexity and the structure of the prescriptive optimization model. Our results indicate that the models formulated are highly effective in predicting add-on product sales, and that using the optimization framework built on the predictive model can result in substantial increases in expected sales over baseline policies. △ Less

Submitted 3 April, 2018; originally announced April 2018.

Comments: 29 pages, 3 figures

On Finding Stable and Efficient Solutions for the Team Formation Problem

Authors: Hoda Atef Yekta, David Bergman, Robert Day

Abstract: The assignment of personnel to teams is a fundamental and ubiquitous managerial function, typically involving several objectives and a variety of idiosyncratic practical constraints. Despite the prevalence of this task in practice, the process is seldom approached as a precise optimization problem over the reported preferences of all agents. This is due in part to the underlying computational comp… ▽ More The assignment of personnel to teams is a fundamental and ubiquitous managerial function, typically involving several objectives and a variety of idiosyncratic practical constraints. Despite the prevalence of this task in practice, the process is seldom approached as a precise optimization problem over the reported preferences of all agents. This is due in part to the underlying computational complexity that occurs when quadratic (i.e., intra-team interpersonal) interactions are taken into consideration, and also due to game-theoretic considerations, when those taking part in the process are self-interested agents. Variants of this fundamental decision problem arise in a number of settings, including, for example, human resources and project management, military platooning, sports-league management, ride sharing, data clustering, and in assigning students to group projects. In this paper, we study a mathematical-programming approach to "team formation" focused on the interplay between two of the most common objectives considered in the related literature: economic efficiency (i.e., the maximization of social welfare) and game-theoretic stability (e.g., finding a core solution when one exists). With a weighted objective across these two goals, the problem is modeled as a bi-level binary optimization problem, and transformed into a single-level, exponentially sized binary integer program. We then devise a branch-cut-and-price algorithms and demonstrate its efficacy through an extensive set of simulations, with favorable comparisons to other algorithms from the literature. △ Less

Submitted 1 April, 2018; originally announced April 2018.

Comments: 41 pages (32 pages, plus references and appendix); 10 figures

The Cosmic-Ray Energy Spectrum between 2 PeV and 2 EeV Observed with the TALE detector in monocular mode

Authors: R. U. Abbasi, M. Abe, T. Abu-Zayyad, M. Allen, R. Azuma, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, R. Cady, B. G. Cheon, J. Chiba, M. Chikawa, A. Di Matteo, T. Fujii, K. Fujita, M. Fukushima, G. Furlich, T. Goto, W. Hanlon, M. Hayashi, Y. Hayashi, N. Hayashida, K. Hibino, K. Honda , et al. (116 additional authors not shown)

Abstract: We report on a measurement of the cosmic ray energy spectrum by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector. The TALE air fluorescence detector is also sensitive to the Cherenkov light produced by shower particles. Low energy cosmic rays, in the PeV energy range, are detectable by TALE as "Cherenkov Events". Using these events, we measure the energy spectrum from a lo… ▽ More We report on a measurement of the cosmic ray energy spectrum by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector. The TALE air fluorescence detector is also sensitive to the Cherenkov light produced by shower particles. Low energy cosmic rays, in the PeV energy range, are detectable by TALE as "Cherenkov Events". Using these events, we measure the energy spectrum from a low energy of $\sim 2$ PeV to an energy greater than 100 PeV. Above 100 PeV TALE can detect cosmic rays using air fluorescence. This allows for the extension of the measurement to energies greater than a few EeV. In this paper, we will describe the detector, explain the technique, and present results from a measurement of the spectrum using $\sim 1000$ hours of observation. The observed spectrum shows a clear steepening near $10^{17.1}$ eV, along with an ankle-like structure at $10^{16.2}$ eV. These features present important constraints on galactic cosmic rays origin and propagation models. The feature at $10^{17.1}$ eV may also mark the end of the galactic cosmic rays flux and the start of the transition to extra-galactic sources. △ Less

Submitted 3 March, 2018; originally announced March 2018.

Comments: 26 pages, 24 figures, submitted to ApJ

Journal ref: ApJ 865 74 (2018)