skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Unfolding neutron spectrum with Markov Chain Monte Carlo at MIT research Reactor with He-3 Neutral Current Detectors [Measuring neutron spectrum at MIT research reactor utilizing He-3 Bonner Cylinder Approach with an unfolding analysis]

Abstract

Here, the Ricochet experiment seeks to measure Coherent (neutral-current) Elastic Neutrino-Nucleus Scattering (CEνNS) using dark-matter-style detectors with sub-keV thresholds placed near a neutrino source, such as the MIT (research) Reactor (MITR), which operates at 5.5 MW generating approximately 2.2 × 10 18 ν/second in its core. Currently, Ricochet is characterizing the backgrounds at MITR, the main component of which comes in the form of neutrons emitted from the core simultaneous with the neutrino signal. To characterize this background, we wrapped Bonner cylinders around a 3 2He thermal neutron detector, whose data was then unfolded via a Markov Chain Monte Carlo (MCMC) to produce a neutron energy spectrum across several orders of magnitude. We discuss the resulting spectrum and its implications for deploying Ricochet at the MITR site as well as the feasibility of reducing this background level via the addition of polyethylene shielding around the detector setup.

Authors:
ORCiD logo [1]; ORCiD logo [2];  [3];  [4];  [1];  [5];  [1];  [5];  [1];  [1];  [1]
  1. Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States)
  2. Fermi National Accelerator Lab. (FNAL), Batavia, IL (United States); Univ. of Chicago, Chicago, IL (United States)
  3. Univ. Lyon (France)
  4. Northwestern Univ., Evanston, IL (United States)
  5. Univ. of Wisconsin, Madison, WI (United States)
Publication Date:
Research Org.:
Fermi National Accelerator Lab. (FNAL), Batavia, IL (United States)
Sponsoring Org.:
USDOE Office of Science (SC), High Energy Physics (HEP) (SC-25)
OSTI Identifier:
1407124
Report Number(s):
FERMILAB-PUB-17-441-AE; arXiv:1710.00802
Journal ID: ISSN 1748-0221; 1628106
Grant/Contract Number:
AC02-07CH11359
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Journal of Instrumentation
Additional Journal Information:
Journal Volume: 13; Journal Issue: 02; Journal ID: ISSN 1748-0221
Publisher:
Institute of Physics (IOP)
Country of Publication:
United States
Language:
English
Subject:
73 NUCLEAR PHYSICS AND RADIATION PHYSICS; 46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND TECHNOLOGY; 72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; Neutrino detectors; Neutron detectors (cold, thermal, fast neutrons); Detector modelling and simulations I (interaction of radiation with matter, interaction of photons with matter, interaction of hadrons with matter, etc); Gaseous detectors

Citation Formats

Leder, A., Anderson, A. J., Billard, J., Figueroa-Feliciano, E., Formaggio, J. A., Hasselkus, C., Newman, E., Palladino, K., Phuthi, M., Winslow, L., and Zhang, L. Unfolding neutron spectrum with Markov Chain Monte Carlo at MIT research Reactor with He-3 Neutral Current Detectors [Measuring neutron spectrum at MIT research reactor utilizing He-3 Bonner Cylinder Approach with an unfolding analysis]. United States: N. p., 2018. Web. doi:10.1088/1748-0221/13/02/P02004.
Leder, A., Anderson, A. J., Billard, J., Figueroa-Feliciano, E., Formaggio, J. A., Hasselkus, C., Newman, E., Palladino, K., Phuthi, M., Winslow, L., & Zhang, L. Unfolding neutron spectrum with Markov Chain Monte Carlo at MIT research Reactor with He-3 Neutral Current Detectors [Measuring neutron spectrum at MIT research reactor utilizing He-3 Bonner Cylinder Approach with an unfolding analysis]. United States. doi:10.1088/1748-0221/13/02/P02004.
Leder, A., Anderson, A. J., Billard, J., Figueroa-Feliciano, E., Formaggio, J. A., Hasselkus, C., Newman, E., Palladino, K., Phuthi, M., Winslow, L., and Zhang, L. Fri . "Unfolding neutron spectrum with Markov Chain Monte Carlo at MIT research Reactor with He-3 Neutral Current Detectors [Measuring neutron spectrum at MIT research reactor utilizing He-3 Bonner Cylinder Approach with an unfolding analysis]". United States. doi:10.1088/1748-0221/13/02/P02004.
@article{osti_1407124,
title = {Unfolding neutron spectrum with Markov Chain Monte Carlo at MIT research Reactor with He-3 Neutral Current Detectors [Measuring neutron spectrum at MIT research reactor utilizing He-3 Bonner Cylinder Approach with an unfolding analysis]},
author = {Leder, A. and Anderson, A. J. and Billard, J. and Figueroa-Feliciano, E. and Formaggio, J. A. and Hasselkus, C. and Newman, E. and Palladino, K. and Phuthi, M. and Winslow, L. and Zhang, L.},
abstractNote = {Here, the Ricochet experiment seeks to measure Coherent (neutral-current) Elastic Neutrino-Nucleus Scattering (CEνNS) using dark-matter-style detectors with sub-keV thresholds placed near a neutrino source, such as the MIT (research) Reactor (MITR), which operates at 5.5 MW generating approximately 2.2 × 1018 ν/second in its core. Currently, Ricochet is characterizing the backgrounds at MITR, the main component of which comes in the form of neutrons emitted from the core simultaneous with the neutrino signal. To characterize this background, we wrapped Bonner cylinders around a 32He thermal neutron detector, whose data was then unfolded via a Markov Chain Monte Carlo (MCMC) to produce a neutron energy spectrum across several orders of magnitude. We discuss the resulting spectrum and its implications for deploying Ricochet at the MITR site as well as the feasibility of reducing this background level via the addition of polyethylene shielding around the detector setup.},
doi = {10.1088/1748-0221/13/02/P02004},
journal = {Journal of Instrumentation},
number = 02,
volume = 13,
place = {United States},
year = {Fri Feb 02 00:00:00 EST 2018},
month = {Fri Feb 02 00:00:00 EST 2018}
}

Journal Article:
Free Publicly Available Full Text
This content will become publicly available on February 2, 2019
Publisher's Version of Record

Save / Share:
  • Monte Carlo methods of coupled neutron/photon transport are being used in the design of filtered beams for Neutron Capture Therapy (NCT). This method of beam analysis provides segregation of each individual dose component, and thereby facilitates beam optimization. The Monte Carlo method is discussed in some detail in relation to NCT epithermal beam design. Ideal neutron beams (i.e., plane-wave monoenergetic neutron beams with no primary gamma-ray contamination) have been modeled both for comparison and to establish target conditions for a practical NCT epithermal beam design. Detailed models of the 5 MWt Massachusetts Institute of Technology Research Reactor (MITR-II) together withmore » a polyethylene head phantom have been used to characterize approximately 100 beam filter and moderator configurations. Using the Monte Carlo methodology of beam design and benchmarking/calibrating our computations with measurements, has resulted in an epithermal beam design which is useful for therapy of deep-seated brain tumors. This beam is predicted to be capable of delivering a dose of 2000 RBE-cGy (cJ/kg) to a therapeutic advantage depth of 5.7 cm in polyethylene assuming 30 micrograms/g 10B in tumor with a ten-to-one tumor-to-blood ratio, and a beam diameter of 18.4 cm. The advantage ratio (AR) is predicted to be 2.2 with a total irradiation time of approximately 80 minutes. Further optimization work on the MITR-II epithermal beams is expected to improve the available beams. 20 references.« less
  • We describe a stochastic inversion method for mapping subsurface regions where the electrical resistivity is changing. The technique combines prior information, electrical resistance data and forward models to produce subsurface resistivity models that are most consistent with all available data. Bayesian inference and a Metropolis simulation algorithm form the basis for this approach. Attractive features include its ability to: (1) provide quantitative measures of the uncertainty of a generated estimate and, (2) allow alternative model estimates to be identified, compared and ranked. Methods that monitor convergence and summarize important trends of the posterior distribution are introduced. Results from a physicalmore » model test and a field experiment were used to assess performance. The stochastic inversions presented provide useful estimates of the most probable location, shape, and volume of the changing region, and the most likely resistivity change. The proposed method is computationally expensive, requiring the use of extensive computational resources to make its application practical.« less
  • The Laser Interferometer Space Antenna (LISA) defines new demands on data analysis efforts in its all-sky gravitational wave survey, recording simultaneously thousands of galactic compact object binary foreground sources and tens to hundreds of background sources like binary black hole mergers and extreme-mass ratio inspirals. We approach this problem with an adaptive and fully automatic Reversible Jump Markov Chain Monte Carlo sampler, able to sample from the joint posterior density function (as established by Bayes theorem) for a given mixture of signals ''out of the box'', handling the total number of signals as an additional unknown parameter beside the unknownmore » parameters of each individual source and the noise floor. We show in examples from the LISA Mock Data Challenge implementing the full response of LISA in its TDI description that this sampler is able to extract monochromatic Double White Dwarf signals out of colored instrumental noise and additional foreground and background noise successfully in a global fitting approach. We introduce 2 examples with fixed number of signals (MCMC sampling), and 1 example with unknown number of signals (RJ-MCMC), the latter further promoting the idea behind an experimental adaptation of the model indicator proposal densities in the main sampling stage. We note that the experienced runtimes and degeneracies in parameter extraction limit the shown examples to the extraction of a low but realistic number of signals.« less
  • We introduce a fast Markov Chain Monte Carlo (MCMC) exploration of the astrophysical parameter space using a modified version of the publicly available code Code Investigating GALaxy Emission (CIGALE). The original CIGALE builds a grid of theoretical spectral energy distribution (SED) models and fits to photometric fluxes from ultraviolet to infrared to put constraints on parameters related to both formation and evolution of galaxies. Such a grid-based method can lead to a long and challenging parameter extraction since the computation time increases exponentially with the number of parameters considered and results can be dependent on the density of sampling points,more » which must be chosen in advance for each parameter. MCMC methods, on the other hand, scale approximately linearly with the number of parameters, allowing a faster and more accurate exploration of the parameter space by using a smaller number of efficiently chosen samples. We test our MCMC version of the code CIGALE (called CIGALEMC) with simulated data. After checking the ability of the code to retrieve the input parameters used to build the mock sample, we fit theoretical SEDs to real data from the well-known and -studied Spitzer Infrared Nearby Galaxy Survey sample. We discuss constraints on the parameters and show the advantages of our MCMC sampling method in terms of accuracy of the results and optimization of CPU time.« less
  • We use the Markov Chain Monte Carlo method to investigate a global constraints on the generalized Chaplygin gas (GCG) model as the unification of dark matter and dark energy from the latest observational data: the Constitution dataset of type supernovae Ia (SNIa), the observational Hubble data (OHD), the cluster X-ray gas mass fraction, the baryon acoustic oscillation (BAO), and the cosmic microwave background (CMB) data. In a non-flat universe, the constraint results for GCG model are, Ω{sub b}h{sup 2} = 0.0235{sup +0.0021}{sub −0.0018} (1σ) {sup +0.0028}{sub −0.0022} (2σ), Ω{sub k} = 0.0035{sup +0.0172}{sub −0.0182} (1σ) {sup +0.0226}{sub −0.0204} (2σ), A{submore » s} = 0.753{sup +0.037}{sub −0.035} (1σ) {sup +0.045}{sub −0.044} (2σ), α = 0.043{sup +0.102}{sub −0.106} (1σ) {sup +0.134}{sub −0.117} (2σ), and H{sub 0} = 70.00{sup +3.25}{sub −2.92} (1σ) {sup +3.77}{sub −3.67} (2σ), which is more stringent than the previous results for constraint on GCG model parameters. Furthermore, according to the information criterion, it seems that the current observations much support ΛCDM model relative to the GCG model.« less