Stephen J. Parke
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Research Interests
Stephen J. Parke has a total of 1000 co-authors. Here are the top co-authors by number of shared publications:
- Peter B. Denton
Brookhaven National Laboratory
29 shared publications - Renata Zukanovich Funchal
Universidade de São Paulo
18 shared publications - Olga Mena
Instituto de Fisica Corpuscular
13 shared publications - Pedro Accioly Nogueira Machado
Particle Theory Department, Fermilab
13 shared publications Gregory Mahlon
10 shared publications- Hisakazu Minakata
Department of Physics Tokyo Metropolitan University 1-1 Minami-Osawa 192-0397 Hachioji Tokyo Japan
10 shared publications
Latest Publications
- Open Access
- 1 August 2026
- Published by IOP Publishing in Journal of Instrumentation
- Vol. 21 (08) , P08018
Charge readout electronics for the DUNE horizontal drift far detector: design and performance in ProtoDUNE-HD
DUNE (Deep Underground Neutrino Experiment) is a long-baseline neutrino oscillation experiment currently under construction, whose far detectors will be the largest liquid argon time projection chambers ever built. This detector design calls for custom-built cryogenic front-end electronics to meet its performance requirements. This paper describes the charge readout electronics that will be used in the DUNE horizontal drift (HD) far detector and presents performance results using data from the ProtoDUNE-HD detector, a 770 ton liquid argon time projection chamber operated at the CERN Neutrino Platform in 2024 that served as the final prototype of the DUNE HD design.
- 1 August 2026
- Published by IOP Publishing in Journal of Instrumentation
- Vol. 21 (08) , P08018
Open AccessCharge readout electronics for the DUNE horizontal drift far detector: design and performance in ProtoDUNE-HD
DUNE (Deep Underground Neutrino Experiment) is a long-baseline neutrino oscillation experiment currently under construction, whose far detectors will be the largest liquid argon time projection chambers ever built. This detector design calls for custom-built cryogenic front-end electronics to meet its performance requirements. This paper describes the charge readout electronics that will be used in the DUNE horizontal drift (HD) far detector and presents performance results using data from the ProtoDUNE-HD detector, a 770 ton liquid argon time projection chamber operated at the CERN Neutrino Platform in 2024 that served as the final prototype of the DUNE HD design.
- Open Access
- 3 June 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 113 (11)
Measurement of exclusive π + -argon interactions using ProtoDUNE-SP
We present the measurement of π + -argon inelastic cross sections using the ProtoDUNE single-phase liquid argon time projection chamber in the incident π + kinetic energy range of 500–800 MeV in multiple exclusive channels (absorption, charge exchange, and the remaining inelastic interactions). The results of this analysis are important inputs to simulations of liquid argon neutrino experiments such as the Deep Underground Neutrino Experiment and the Short Baseline Neutrino program at Fermi National Accelerator Laboratory. They will be employed to improve the modeling of final state interactions within neutrino event generators used by these experiments, as well as the modeling of π + -argon secondary interactions within the liquid argon. This is the first measurement of π + -argon absorption at this kinetic energy range as well as the first ever measurement of π + -argon charge exchange.
- 3 June 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 113 (11)
Open AccessMeasurement of exclusive π + -argon interactions using ProtoDUNE-SP
We present the measurement of π + -argon inelastic cross sections using the ProtoDUNE single-phase liquid argon time projection chamber in the incident π + kinetic energy range of 500–800 MeV in multiple exclusive channels (absorption, charge exchange, and the remaining inelastic interactions). The results of this analysis are important inputs to simulations of liquid argon neutrino experiments such as the Deep Underground Neutrino Experiment and the Short Baseline Neutrino program at Fermi National Accelerator Laboratory. They will be employed to improve the modeling of final state interactions within neutrino event generators used by these experiments, as well as the modeling of π + -argon secondary interactions within the liquid argon. This is the first measurement of π + -argon absorption at this kinetic energy range as well as the first ever measurement of π + -argon charge exchange.
- Open Access
- 10 March 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 113 (5)
Efficient atmospheric, solar, and supernova neutrino propagation through the Earth
Algorithms for computing neutrino oscillation probabilities in sharply varying matter potentials such as the Earth are becoming increasingly important. As the next generation of experiments, DUNE and HyperK as well as the IceCube upgrade and KM3NeT, come online, the computational cost for atmospheric and solar neutrinos will continue to increase. To address these issues, we expand upon our previous algorithm for long-baseline calculations to efficiently handle probabilities through the Earth for atmospheric, nighttime solar, and supernova neutrinos. The algorithm is fast, flexible, and accurate. It can handle arbitrary Earth models with two different schemes for varying density profiles. We also provide a ++ implementation of the code called uast-arth along with a detailed user manual. The code intelligently keeps track of repeated calculations and only recalculates what is needed on each successive call which can also help provide significant speed-ups.
- 10 March 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 113 (5)
Open AccessEfficient atmospheric, solar, and supernova neutrino propagation through the Earth
Algorithms for computing neutrino oscillation probabilities in sharply varying matter potentials such as the Earth are becoming increasingly important. As the next generation of experiments, DUNE and HyperK as well as the IceCube upgrade and KM3NeT, come online, the computational cost for atmospheric and solar neutrinos will continue to increase. To address these issues, we expand upon our previous algorithm for long-baseline calculations to efficiently handle probabilities through the Earth for atmospheric, nighttime solar, and supernova neutrinos. The algorithm is fast, flexible, and accurate. It can handle arbitrary Earth models with two different schemes for varying density profiles. We also provide a ++ implementation of the code called uast-arth along with a detailed user manual. The code intelligently keeps track of repeated calculations and only recalculates what is needed on each successive call which can also help provide significant speed-ups.
- Open Access
- 30 September 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 112 (5)
Exploring quantum statistics for massive Dirac and Majorana neutrinos using spinor-helicity techniques
Recently, there has been interest in the applicability of quantum statistics to distinguish Dirac from Majorana neutrinos in multineutrino final states. In particular, debate has arisen over the validity of the Dirac-Majorana confusion theorem in these processes, i.e., that any distinction between the Dirac and Majorana processes goes to zero as the neutrino mass goes to zero. Here we approach this problem equipped with spinor-helicity methods generalized for massive Dirac and Majorana fermions. We explicitly calculate all helicity amplitudes, and their squares, for the decay of a light scalar particle to two neutrinos and two oppositely charged leptons. This allows us to pinpoint the crucial steps which could lead to claims of a violation of the confusion theorem. We show that, if the correct antisymmetrization of Dirac to Majorana amplitudes is used, identification of which is clear in this framework, and all relevant contributions are appropriately summed, a scalar decay into two charged leptons and two neutrinos satisfies the Dirac-Majorana confusion theorem.
- 30 September 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 112 (5)
Open AccessExploring quantum statistics for massive Dirac and Majorana neutrinos using spinor-helicity techniques
Recently, there has been interest in the applicability of quantum statistics to distinguish Dirac from Majorana neutrinos in multineutrino final states. In particular, debate has arisen over the validity of the Dirac-Majorana confusion theorem in these processes, i.e., that any distinction between the Dirac and Majorana processes goes to zero as the neutrino mass goes to zero. Here we approach this problem equipped with spinor-helicity methods generalized for massive Dirac and Majorana fermions. We explicitly calculate all helicity amplitudes, and their squares, for the decay of a light scalar particle to two neutrinos and two oppositely charged leptons. This allows us to pinpoint the crucial steps which could lead to claims of a violation of the confusion theorem. We show that, if the correct antisymmetrization of Dirac to Majorana amplitudes is used, identification of which is clear in this framework, and all relevant contributions are appropriately summed, a scalar decay into two charged leptons and two neutrinos satisfies the Dirac-Majorana confusion theorem.
- Open Access
Exploring Quantum Statistics for Dirac and Majorana Neutrinos using Spinor-Helicity techniques
Recently, there has been interest in the applicability of quantum statistics to distinguish Dirac from Majorana neutrinos in multi-neutrino final states. In particular, debate has arisen over the validity of the Dirac-Majorana confusion theorem in these processes, i.e. that any distinction between the Dirac and Majorana processes goes to zero as the neutrino mass goes to zero. Here we approach this problem equipped with spinor-helicity methods generalized for massive Dirac and Majorana fermions. We explicitly calculate all helicity amplitudes for the decay of a light scalar particle to two neutrinos and two oppositely charged leptons. This allows us to pinpoint the crucial steps which could lead to claims of a violation of the confusion theorem. We show that if the correct anti-symmetrization of Dirac to Majorana amplitudes is used, identification of which is clear in this framework, and all relevant contributions are appropriately summed, a scalar decay into two charged leptons and two neutrinos satisfies the Dirac-Majorana confusion theorem.
Open AccessExploring Quantum Statistics for Dirac and Majorana Neutrinos using Spinor-Helicity techniques
Recently, there has been interest in the applicability of quantum statistics to distinguish Dirac from Majorana neutrinos in multi-neutrino final states. In particular, debate has arisen over the validity of the Dirac-Majorana confusion theorem in these processes, i.e. that any distinction between the Dirac and Majorana processes goes to zero as the neutrino mass goes to zero. Here we approach this problem equipped with spinor-helicity methods generalized for massive Dirac and Majorana fermions. We explicitly calculate all helicity amplitudes for the decay of a light scalar particle to two neutrinos and two oppositely charged leptons. This allows us to pinpoint the crucial steps which could lead to claims of a violation of the confusion theorem. We show that if the correct anti-symmetrization of Dirac to Majorana amplitudes is used, identification of which is clear in this framework, and all relevant contributions are appropriately summed, a scalar decay into two charged leptons and two neutrinos satisfies the Dirac-Majorana confusion theorem.
- Open Access
- 12 May 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (9)
Supernova pointing capabilities of DUNE
The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electron-neutrino charged-current absorption on Ar 40 and elastic scattering of neutrinos on electrons. Procedures to reconstruct individual interactions, including a newly developed technique called “brems flipping,” as well as the burst direction from an ensemble of interactions are described. Performance of the burst direction reconstruction is evaluated for supernovae happening at a distance of 10 kpc for a specific supernova burst flux model. The pointing resolution is found to be 3.4 degrees at 68% coverage for a perfect interaction-channel classification and a fiducial mass of 40 kton, and 6.6 degrees for a 10 kton fiducial mass respectively. Assuming a 4% rate of charged-current interactions being misidentified as elastic scattering, DUNE’s burst pointing resolution is found to be 4.3 degrees (8.7 degrees) at 68% coverage. Published by the American Physical Society 2025
- 12 May 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (9)
Open AccessSupernova pointing capabilities of DUNE
The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electron-neutrino charged-current absorption on Ar 40 and elastic scattering of neutrinos on electrons. Procedures to reconstruct individual interactions, including a newly developed technique called “brems flipping,” as well as the burst direction from an ensemble of interactions are described. Performance of the burst direction reconstruction is evaluated for supernovae happening at a distance of 10 kpc for a specific supernova burst flux model. The pointing resolution is found to be 3.4 degrees at 68% coverage for a perfect interaction-channel classification and a fiducial mass of 40 kton, and 6.6 degrees for a 10 kton fiducial mass respectively. Assuming a 4% rate of charged-current interactions being misidentified as elastic scattering, DUNE’s burst pointing resolution is found to be 4.3 degrees (8.7 degrees) at 68% coverage. Published by the American Physical Society 2025
- Open Access
- 31 January 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (1)
Mass ordering sum rule for the neutrino disappearance channels in T2K, NOvA, and JUNO
We revisit, reformulate, and extend a method for determining the neutrino mass ordering by using precision measurements of the atmospheric Δ m 2 s in both electron and muon neutrino disappearance channels, first proposed by the authors in 2005 [H. Nunokawa , ]. The mass ordering is a very important outstanding question for our understanding of the elusive neutrino and determination of the mass ordering has consequences to particle physics, nuclear physics, and cosmology. The JUNO reactor experiment will start data taking this year, and the precision of the atmospheric Δ m 2 s from electron antineutrino measurements will improve by a factor of 3 from Daya Bay’s 2.4% to 0.8% within a year. This measurement, when combined with the atmospheric Δ m 2 ’s measurements from T2K and NOvA for muon neutrino disappearance, will contribute substantially to the Δ χ 2 between the two remaining neutrino mass orderings. In this paper we derive for the first time a mass ordering sum rule that can be used to address the possibility that JUNO’s atmospheric Δ m 2 ’s measurement, when combined with other experiments in particular T2K and NOvA, can determine the neutrino mass ordering at the 3 σ confidence level within one year of operation. For a confidence level of 5 σ in a single experiment, we will have to wait until the middle of the next decade when the DUNE experiment is operating.
- 31 January 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (1)
Open AccessMass ordering sum rule for the neutrino disappearance channels in T2K, NOvA, and JUNO
We revisit, reformulate, and extend a method for determining the neutrino mass ordering by using precision measurements of the atmospheric Δ m 2 s in both electron and muon neutrino disappearance channels, first proposed by the authors in 2005 [H. Nunokawa , ]. The mass ordering is a very important outstanding question for our understanding of the elusive neutrino and determination of the mass ordering has consequences to particle physics, nuclear physics, and cosmology. The JUNO reactor experiment will start data taking this year, and the precision of the atmospheric Δ m 2 s from electron antineutrino measurements will improve by a factor of 3 from Daya Bay’s 2.4% to 0.8% within a year. This measurement, when combined with the atmospheric Δ m 2 ’s measurements from T2K and NOvA for muon neutrino disappearance, will contribute substantially to the Δ χ 2 between the two remaining neutrino mass orderings. In this paper we derive for the first time a mass ordering sum rule that can be used to address the possibility that JUNO’s atmospheric Δ m 2 ’s measurement, when combined with other experiments in particular T2K and NOvA, can determine the neutrino mass ordering at the 3 σ confidence level within one year of operation. For a confidence level of 5 σ in a single experiment, we will have to wait until the middle of the next decade when the DUNE experiment is operating.
- Open Access
- 29 January 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (1)
Exploring the interference between the atmospheric and solar neutrino oscillation subamplitudes
The interference between the atmospheric and solar neutrino oscillation subamplitudes is said to be responsible for C P violation () in neutrino appearance channels. More precisely, is generated by the interference between the parts of the neutrino oscillation amplitude that are C P even and C P odd: even or odd when the neutrino mixing matrix is replaced with its complex conjugate. This is the interference term, as it gives a contribution to the oscillation probability, the square of the amplitude, which is opposite in sign for neutrinos and antineutrinos and is unique. For this interference to be nonzero, at least two subamplitudes are required. There are, however, other interference terms, which are even under the above exchange, and these are the C P conserving () interference terms. In this paper, we explore in detail these interference terms and show that they cannot be uniquely defined, as one can move pieces of the amplitude from the atmospheric subamplitude to the solar subamplitude and vice versa. This freedom allows one to move the interference terms around, but does not let you eliminate them completely. We also show that there is a reasonable definition of the atmospheric and solar subamplitudes for the appearance channels such that in neutrino disappearance probability there is no atmospheric-solar interference term. However, with this choice, there is a interference term within the atmospheric sector.
- 29 January 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (1)
Open AccessExploring the interference between the atmospheric and solar neutrino oscillation subamplitudes
The interference between the atmospheric and solar neutrino oscillation subamplitudes is said to be responsible for C P violation () in neutrino appearance channels. More precisely, is generated by the interference between the parts of the neutrino oscillation amplitude that are C P even and C P odd: even or odd when the neutrino mixing matrix is replaced with its complex conjugate. This is the interference term, as it gives a contribution to the oscillation probability, the square of the amplitude, which is opposite in sign for neutrinos and antineutrinos and is unique. For this interference to be nonzero, at least two subamplitudes are required. There are, however, other interference terms, which are even under the above exchange, and these are the C P conserving () interference terms. In this paper, we explore in detail these interference terms and show that they cannot be uniquely defined, as one can move pieces of the amplitude from the atmospheric subamplitude to the solar subamplitude and vice versa. This freedom allows one to move the interference terms around, but does not let you eliminate them completely. We also show that there is a reasonable definition of the atmospheric and solar subamplitudes for the appearance channels such that in neutrino disappearance probability there is no atmospheric-solar interference term. However, with this choice, there is a interference term within the atmospheric sector.