---
title: 'NEVOD-DECOR: Cosmic-Ray Detector'
url: https://www.emergentmind.com/topics/nevod-decor
type: topic
---

# NEVOD-DECOR: Cosmic-Ray Detector

Searching arXiv for NEVOD-DECOR papers to ground the article in the latest and directly relevant literature.
Tool unavailable in this environment. Proceeding with the provided arXiv records and details.
NEVOD-DECOR is a surface-based cosmic-ray installation at MEPhI (Moscow) built around the NEVOD Cherenkov water calorimeter and the DECOR coordinate-tracking detector. In the published literature, it functions as a hybrid system for the study of inclined muon bundles in extensive air showers, local muon density spectra (LMDS), calorimetric energy deposits in water, near-horizontal and upward-going albedo muons, and large-scale cosmic-ray anisotropy; later analyses extend the reconstructed primary-energy coverage to approximately \(2\,\mathrm{PeV}\)–\(3\,\mathrm{EeV}\) [2508.05274].

## 1. Configuration of the experimental complex

The NEVOD complex is centered on a Cherenkov water detector at ground level in the MEPhI campus. The water calorimeter has an active volume of \(2000\,\mathrm{m^3}\), corresponding to a \(26\,\mathrm{m}\times 9\,\mathrm{m}\times 9\,\mathrm{m}\) tank, and is instrumented with 91 quasi-spherical measuring modules arranged on 25 vertical strings. Each quasi-spherical module contains six FEU-200 photomultipliers oriented along the coordinate axes, providing near-isotropic sensitivity to Cherenkov light. The two-dynode readout extends the dynamic range from 1 to \(10^5\) photoelectrons [2208.05926].

Around the water volume, DECOR is deployed as eight vertical supermodules mounted in the surrounding galleries. In the NEVOD descriptions used for LMDS analyses, each supermodule comprises eight vertical planes of streamer-tube chambers with orthogonal strip readout and an effective area of \(8.4\,\mathrm{m^2}\) per supermodule; the full DECOR acceptance for muon bundles is approximately \(70\,\mathrm{m^2\,sr}\), with angular resolution better than \(0.8^\circ\). Other technical descriptions give each supermodule as \(3.1\times 2.7\,\mathrm{m^2}\), again with eight planes and spatial and angular resolutions at the \(\lesssim 1\,\mathrm{cm}\) and \(\lesssim 1^\circ\) level, respectively [1701.06759].

Several auxiliary subsystems have been integrated into the complex. The Calibration Telescope System (CTS) consists of two horizontal planes of 40 plastic scintillation counters each, arranged in a chessboard pattern over an area \(8\times 10\,\mathrm{m^2}\) above and below the water volume; after the 2013 upgrade, its bottom plane also recorded the muon component of inclined air showers with an effective angular acceptance peaked at \(\theta_{\mathrm{eff}}\simeq 29^\circ\) for muon bundles [1701.06759]. The broader NEVOD infrastructure also includes the NEVOD-EAS cluster-type array for air-shower size, axis, and direction reconstruction, and the PRISMA-32 en-detector array for simultaneous measurements of hadronic and electromagnetic air-shower components [1612.09441], [1512.02042].

A recurring design feature is the use of shielding by the water tank itself. In selected azimuth sectors, six of the eight DECOR supermodules are shadowed by the NEVOD water volume, raising the muon detection threshold to approximately \(2\,\mathrm{GeV}\) and suppressing electromagnetic and hadronic contamination in inclined-bundle measurements [2208.05926].

## 2. Detection principles and core observables

In DECOR, a muon bundle is defined as a group of spatially and temporally correlated muon tracks reconstructed in one or more supermodules. Bundle multiplicity \(m\) is the number of muon tracks in the event. In analyses of inclined bundles, DECOR reconstructs individual tracks in each supermodule, imposes a parallelism cut within a \(5^\circ\) cone, and determines the bundle direction \((\theta,\phi)\) from a global straight-line fit to the track segments [1701.06759].

The calorimetric observable is supplied by NEVOD. The detector records the total Cherenkov-light yield, written either as \(\Sigma=\sum(\mathrm{p.e.})\) or as \(E N_{\mathrm{pe}}=\sum_i N_{\mathrm{phe},i}\), across all photomultipliers. Since muon energy loss in water follows
\[
\frac{dE}{dX}\simeq a+bE_\mu,
\]
the total light yield is, to first order, proportional to the sum of muon energies and to the total energy deposit of the bundle in water [2208.05926], [1612.09175].

The principal density observable is the local muon density. In the 2017 LMDS formulation it is defined as
\[
\rho_\mu=\frac{m}{A_{\mathrm{eff}}(\theta)},
\]
where \(m\) is the reconstructed multiplicity and \(A_{\mathrm{eff}}(\theta)\) is the effective detection area projected into the shower front. Later NEVOD-DECOR papers write the same quantity as
\[
D=\frac{m}{S_{\mathrm{det}}(\theta,\phi)},
\]
or, in the bias-corrected form used for energy-deposit studies,
\[
D=\frac{m-B}{S_{\mathrm{det}}},
\]
where \(B\) accounts for the bias from Poisson fluctuations and the integral slope \(\beta\) of the density spectrum [1701.06759], [2208.05926], [1612.09175].

For calorimetric studies of bundles, the specific deposit is introduced to remove the trivial scaling with multiplicity and acceptance. Two equivalent notations appear in the literature:
\[
\frac{\Sigma}{D}
\quad\text{and}\quad
\epsilon\equiv \frac{E N_{\mathrm{pe}}}{D}.
\]
These observables are used as proxies for the average muon energy in the bundle after detector-response calibration with Geant4 [2208.05926], [1612.09175].

CTS implements a different but related counting principle. There, muon bundles are recorded as coincidences of hit scintillation counters in the bottom plane, and the hit-counter multiplicity \(k\) is used as a proxy for bundle size. In the CTS analyses, the effective bundle angle \(\theta_{\mathrm{eff}}\approx 29^\circ\) is taken from the measured muon-angular distribution proportional to \(\cos^{4.5}\theta\) [1701.06759].

## 3. Local muon density spectra methodology

The LMDS method is the central analysis framework of NEVOD-DECOR. At fixed zenith angle, the differential local muon density spectrum is assumed to follow a power law,
\[
\frac{dN}{d\rho_\mu}(\theta)=C(\theta)\,\rho_\mu^{-\gamma(\theta)},
\]
with spectral slope \(\gamma(\theta)\) obtained from fits of \(\ln[dN/d\rho_\mu]\) versus \(\ln\rho_\mu\) in a specified density interval [1701.06759].

In the later notation, the measured quantity is the event flux per unit density, per unit solid angle, and per unit time,
\[
\frac{dF(D,\theta)}{dD}.
\]
A phenomenological reference form,
\[
\frac{dF_0}{dD}=C\,D^{-(\beta+1)}\cos^\alpha\theta,
\]
with \(C=7\times 10^{-4}\,\mathrm{m^2\,s^{-1}\,sr^{-1}}\), \(\beta=2.1\), and \(\alpha=4.7\), is fitted to the joint \((m,\theta)\) distribution by maximum likelihood. The observed LMDS in each \((\Delta D,\Delta\theta)\) cell is then obtained from the ratio \(N_{\mathrm{obs}}/N_{\mathrm{exp}}\) multiplied by the reference spectrum, where \(N_{\mathrm{exp}}\) comes from a full Monte Carlo simulation of detector response including Poisson fluctuations and track-masking effects [2508.05274].

The theoretical LMDS is constructed by convolving a primary spectrum with simulated muon lateral distribution functions from CORSIKA. In one form used for bundle-intensity modeling,
\[
\frac{dF(D,\theta)}{dD}=\int_{E_0} J(E_0)\,P(D\mid E_0,\theta)\,dE_0,
\]
where \(P(D\mid E_0,\theta)\,dD\) is the probability that a shower of primary energy \(E_0\) and zenith \(\theta\) produces a local density between \(D\) and \(D+dD\) at the observation point. In the continuum formulation of the 2025 analysis, the integral LMDS is written as
\[
F(\ge D)=\iint N(\ge E_0(r,D))\,dS,
\]
with \(E_0\) determined from the muon lateral distribution function \(\rho(E_0,r)\) by solving \(\rho(E_0,r)=D\) [2208.05926], [2508.05274].

The mapping from local density to primary energy is simulation-based. In the 2017 reconstruction, full CORSIKA (\(v6.900\)) simulations of the shower muon component for pure proton and iron primaries, using the QGSJET model, yield an approximate calibration curve
\[
\langle E_0\rangle(\rho_\mu,\theta)\simeq A(\theta)\,\rho_\mu^{\,\alpha(\theta)},
\]
with \(\alpha(\theta)\approx 1.0\)–\(1.1\). In the anisotropy analysis, a more explicit empirical estimator is used:
\[
\log_{10}E_{\mathrm{est}}
=
7.03
+
1.07\,\log_{10}(D)
+
3.80\,\log_{10}(\sec\theta),
\]
with energies in GeV and \(D\) in \(\mathrm{m^{-2}}\) [1701.06759], [2210.13916].

The simulation backbone has evolved with time. The published analyses employ CORSIKA \(v6.900\), \(v7.5\), \(7.69\), and \(7.6900\); hadronic-model comparisons include QGSJET, QGSJET-II-04, SIBYLL-2.3, SIBYLL-2.3c, EPOS-LHC, and later EPOS LHC-R, QGSJET-III-01, and SIBYLL-2.3e. The 2025 study also uses two-dimensional muon LDFs to account for geomagnetic-induced asymmetries [1701.06759], [1612.09175], [2208.05926], [2508.05274].

## 4. Spectral features: the knee, the second knee, and composition-sensitive trends

One of the best-known NEVOD-DECOR results is the observation of a steepening in the LMDS around a primary energy of about \(10^{17}\,\mathrm{eV}\). In the 2017 analysis, DECOR reconstructed differential LMDS for effective zenith angles \(57^\circ\), \(64^\circ\), and \(72^\circ\), while CTS measured the near-vertical region at \(\theta_{\mathrm{eff}}\approx 29^\circ\). In both subsystems, two power-law fits were performed in intervals corresponding to \(10^{16}\)–\(10^{17}\,\mathrm{eV}\) and \(>10^{17}\,\mathrm{eV}\), and in both cases the slope increased above the characteristic density mapped to \(10^{17}\,\mathrm{eV}\) [1701.06759].

| Measurement | Low-\(E_0\) integral slope | High-\(E_0\) integral slope / change |
|---|---:|---:|
| DECOR combined | \(2.107 \pm 0.014\) | \(2.303 \pm 0.054\), \(\Delta\beta=0.196 \pm 0.056\) |
| CTS | \(2.057 \pm 0.009\) | \(2.310 \pm 0.120\), \(\Delta B=0.25 \pm 0.12\) |

For DECOR, the change in slope was reported as \(\Delta\beta=0.196\pm 0.056\), corresponding to \(>3.5\sigma\). For CTS, \(\Delta B=0.25\pm 0.12\), corresponding to \(>2.1\sigma\), with additional systematic uncertainties from the 1.3 correction factor for apparent density enhancement and from the assumed angular distribution contributing \(\lesssim 0.05\) to \(\Delta B\) [1701.06759].

The paper interprets this “second knee” as evidence that the muon-density spectrum steepens near \(E_0\approx 10^{17}\,\mathrm{eV}\), indicating a change in the primary cosmic-ray composition and/or interaction characteristics above approximately \(100\,\mathrm{PeV}\). It explicitly notes that, in rigidity-dependent acceleration models, such a feature may correspond to the cutoff of the heaviest Galactic components or to the onset of extragalactic protons. The same paper states that the concordance of DECOR and CTS measurements across \(15^\circ\)–\(75^\circ\) zenith angles and \(10^{15}\)–\(10^{18}\,\mathrm{eV}\) provides robust experimental evidence for this feature in the muon sector, complementary to electron-size measurements by KASCADE-Grande, Tunka-133, and IceTop [1701.06759].

The longer-baseline 2012–2023 analysis extends these composition-sensitive inferences. Using LMDS from \(40^\circ\) to \(85^\circ\), NEVOD-DECOR reconstructs an all-particle spectrum from \(2\times 10^{15}\) to \(3\times 10^{18}\,\mathrm{eV}\), recovering the knee at \(E_k\approx (3\)–\(4)\times 10^{15}\,\mathrm{eV}\), a slight flattening identified as a “second knee” near \(\sim 10^{17}\,\mathrm{eV}\), and points approaching the ankle region near \(\sim 3\times 10^{18}\,\mathrm{eV}\). In the same analysis, the WHISP-defined \(z\)-parameter rises from \(z\sim 0.2\)–0.4 at \(10^{15}\,\mathrm{eV}\) to \(z\to 0.9\)–1.0 near \(10^{18}\,\mathrm{eV}\) under the stated modeling assumptions, which the paper describes as an extremely heavy composition [2508.05274].

## 5. Muon-bundle calorimetry and the “muon puzzle”

A distinct but closely related NEVOD-DECOR result concerns the discrepancy between observed muon-bundle intensities and the predictions of widely used hadronic interaction models. In the 2022 inclined-bundle analysis, data from May 2012 to March 2021 were selected with \(m\ge 5\), yielding \(99.6\,\mathrm{k}\) events for \(\theta\ge 55^\circ\) over \(58.3\times 10^3\,\mathrm{h}\) and \(30.4\,\mathrm{k}\) events for \(40^\circ\le \theta<55^\circ\) over \(6.3\times 10^3\,\mathrm{h}\). LMDS were measured in nine zenith intervals from \(40^\circ\) to \(90^\circ\) and compared with CORSIKA \(7.69\) calculations using QGSJET-II-04, SIBYLL-2.3c, and EPOS-LHC for pure proton and pure iron primaries [2208.05926].

At moderate angles, corresponding to \(E_0\sim 10^{16}\,\mathrm{eV}\), the data follow the light-composition proton curves. At the largest angles, corresponding to \(E_0\sim 10^{18}\,\mathrm{eV}\), the data lie close to or exceed the iron curves in all three interaction models. The same analysis maps the measurements to the WHISP \(z\)-scale,
\[
z=
\frac{\ln N_\mu^{\mathrm{obs}}-\ln N_\mu^{p,\mathrm{sim}}}
{\ln N_\mu^{Fe,\mathrm{sim}}-\ln N_\mu^{p,\mathrm{sim}}},
\]
and reports that \(z\) rises steadily above \(10^{17}\,\mathrm{eV}\), reaching \(z\simeq 0.8\)–1.0 around \(10^{18}\,\mathrm{eV}\) in all three models. Under those assumptions, the required composition is described as “extremely heavy” [2208.05926].

The contradiction with fluorescence measurements of \(\langle X_{\max}\rangle\), which favor light primaries with \(A\simeq 1\)–4 at \(10^{18}\,\mathrm{eV}\), is identified as the “muon puzzle.” The NEVOD-DECOR papers formulate the problem in explicitly model-comparison terms: either modern interaction models under-produce high-energy muons, or additional processes become important at laboratory energies \(\ge 100\,\mathrm{TeV}\) [2208.05926].

NEVOD-DECOR addresses this discrepancy not only through counting but also through calorimetry. Using the Geant4 \(v10.7\) detector model, tuned on single near-horizontal muons with \(E\approx 100\,\mathrm{GeV}\), the response to artificial bundles of known energy and density is established. With the normalized specific deposit, the average muon energy is extracted as
\[
\langle E_\mu\rangle
=
100\,\mathrm{GeV}\times
\frac{(\Sigma/D)_{\mathrm{exp}}}{(\Sigma/D)_{\mathrm{MC@100\,GeV}}}.
\]
The measured \(\langle E_\mu\rangle\) at fixed density increases from approximately \(70\,\mathrm{GeV}\) at \(\theta\approx 40^\circ\) to approximately \(200\)–\(250\,\mathrm{GeV}\) near \(\theta\approx 85^\circ\). For \(65^\circ\le \theta\le 75^\circ\), the quoted values are approximately \(60\,\mathrm{GeV}\) at \(D\simeq 0.02\,\mathrm{m^{-2}}\), \(90\,\mathrm{GeV}\) at \(0.1\,\mathrm{m^{-2}}\), \(120\,\mathrm{GeV}\) at \(0.5\,\mathrm{m^{-2}}\), and \(220\,\mathrm{GeV}\) at \(1.0\,\mathrm{m^{-2}}\). Above \(E_0\gtrsim 10^{17}\,\mathrm{eV}\), the measured \(\langle E_\mu\rangle\) exceeds model predictions by \(2.6\sigma\)–\(3.8\sigma\) in the last four points [2208.05926].

An earlier energy-deposit study had already reported a milder precursor of the same tendency. For \(D>1\,\mathrm{m^{-2}}\), corresponding to \(E_0\gtrsim 10^{17}\,\mathrm{eV}\), the specific deposit \(\epsilon\) showed an approximately \(10\)–20% upward deviation from CORSIKA \(v7.5\) plus SIBYLL-2.3 expectations, although the last density bins had limited statistics (\(N_{\mathrm{events}}=9\)–49) [1612.09175]. This suggests continuity between the earlier bundle-calorimetry anomaly and the later high-statistics formulation of the muon puzzle.

## 6. Complementary programs and subsystem extensions

NEVOD-DECOR has also been used for measurements of albedo muons, defined as atmospheric muons scattered in the ground into the upper hemisphere after traversing tens of meters of soil. These were studied in the zenith range \(85^\circ\le \theta\le 95^\circ\) using coincident triggers in two oppositely mounted DECOR supermodules. Two data sets, \(10\,548\,\mathrm{h}\) in 2002–2004 and \(19\,843\,\mathrm{h}\) in 2011–2015, were combined for a total live time of about \(30\,000\,\mathrm{h}\). Across both runs, \(5.46\times 10^6\) near-horizontal muons were recorded, of which 5,717 were identified as albedo muons in the interval \(\theta=91^\circ\)–\(95^\circ\). Comparison with Monte Carlo transport through soil showed that the point-nucleus Molière model overestimates the flux by up to approximately 50% near \(\theta=90^\circ\), whereas the finite-nucleus model tracks the data within statistical errors across the full \(85^\circ\)–\(95^\circ\) range [1701.02241].

The complex has also been adapted for anisotropy studies using muon bundles as directional tracers of primary cosmic rays. For 2012–2022, about 14 million events in the local zenith-angle range \(15^\circ\)–\(75^\circ\) were accumulated. The published method corrects the bundle rate for meteorological modulation through
\[
F(t)=F_0\,[H_{500}(t)/H_0]^n,
\]
with \(n\approx -1.68\pm 0.02\) and \(H_0\approx 5530\,\mathrm{m}\). After correction, the run-to-run rms spread in bundle rate is reduced from approximately 5–7% to approximately 1–2%. Rayleigh fits in right ascension then give a first-harmonic amplitude \(a=(1.01\pm 0.38)\times 10^{-3}\) at \(\langle E\rangle\approx 1\,\mathrm{PeV}\) and \(a=(0.95\pm 0.69)\times 10^{-3}\) at \(\langle E\rangle\approx 5\,\mathrm{PeV}\), with phases broadly consistent with the Galactic-center direction within the stated uncertainties [2210.13916].

Two associated detector developments broaden the physical reach of the NEVOD-DECOR site. PRISMA-32, constructed in February 2012 on the fourth floor, uses 32 ZnS(Ag)+\(^6\)LiF scintillation en-detectors to measure prompt electromagnetic signals and delayed thermal-neutron captures. Its thermal-neutron lateral distribution is fitted by a double exponential with scales \(r_1=1.0\pm 0.2\,\mathrm{m}\) and \(r_2=11.0\pm 1.5\,\mathrm{m}\), while the mean neutron arrival-time distribution is described by two exponentials with \(\tau_1=0.6\pm 0.1\,\mathrm{ms}\) and \(\tau_2=4.0\pm 0.5\,\mathrm{ms}\) [1512.02042]. NEVOD-EAS, deployed in 2015–2016 as a cluster array over approximately \(10^4\,\mathrm{m^2}\), reconstructs shower size, core position, and arrival direction from scintillation detector stations; the central part achieved an arrival-direction resolution of \(\lesssim 2^\circ\), preliminary core-position resolution of order \(10\,\mathrm{m}\), and shower-size resolution of order 10% [1612.09441].

In the 2025 synopsis, these developments are linked to further expansion of the inclined-bundle program. The paper states that future upgrades, specifically the TREK drift-chamber tracker with approximately \(250\,\mathrm{m^2}\) area and \(3\,\mathrm{mm}\) two-track resolution, together with enhanced NEVOD water calorimetry, are intended to extend the energy range, improve muon counting precision, and enable measurement of muon energy deposits in support of ongoing work on the muon puzzle and extensive-air-shower physics [2508.05274].

Source: https://www.emergentmind.com/topics/nevod-decor