Pith. Sign in

REVIEW

Global Analysis of Nuclear Parton Distribution Functions using Markov Chain Monte Carlo

High Energy Physics - Phenomenology · Submitted 27 Oct 2026 · 271 references

Pith's verdict: First genuine MCMC-based nPDF global fit with a clever Pb-only single-nucleus setup, but the central multimodality claim rests on chain ergodicity that the thesis itself has not established.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on October 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 4 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.37254 (version 1)
Pith IDDGVGDPOB, Pith's permanent record of this paper and its review
CitationNasim Derakhshanian. "Global Analysis of Nuclear Parton Distribution Functions using Markov Chain Monte Carlo." arXiv:2610.37254 (2026). https://pith.science/paper/DGVGDPOB
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

By sampling the nuclear-PDF posterior directly with MCMC instead of expanding the fit statistic around a single minimum, this thesis finds non-Gaussian, multimodal parameter structure in the valence sector that Hessian uncertainty bands…

The reading

Global fits of nuclear parton distributions have long estimated parameter uncertainties with the Hessian method, which expands the fit statistic quadratically around one best-fit point and thereby assumes a Gaussian posterior. This thesis argues that for nuclear PDFs, where data are sparse and parameter correlations are strong, that assumption breaks down, and that direct sampling of the posterior with Markov Chain Monte Carlo (MCMC) reveals structure the Hessian cannot see: skewed margins, flat directions, and secondary minima, most visibly in the valence-quark parameters. To make the case it performs the first global single-nucleus fit of lead PDFs from lead data alone, with the standard mass-number dependence frozen at $A=208$, and a complementary multi-nuclei fit that lets lighter nuclei constrain the lead distributions. The practical consequence is that Hessian-based nuclear-PDF error bands inherit the Gaussian assumption; if the posterior is genuinely multimodal, those bands misstate the uncertainty in exactly the region where new nuclear data are arriving.

What carries the argument

The load-bearing machinery is the adaptive Metropolis-Hastings sampler applied to a ten-parameter nuclear-PDF ansatz (an nCTEQ-style functional form with standard analytic $A$-dependence), with the posterior proportional to $\exp(-\chi^2/2)$ built from a covariance-matrix likelihood. The algorithm adapts its proposal covariance during burn-in and then runs fixed-kernel chains; convergence is assessed through trace plots, cumulative means, and $\Gamma$-method autocorrelation times, with thinning chosen from those times. Two data setups carry the scientific comparison: a Pb-only fit ($A$ frozen at 208, lead data only) that isolates the intrinsic posterior structure, and a Multi-nuclei fit in which the $A$-dependence links lead to lighter nuclei. The observed object that distinguishes the Bayesian picture from the Hessian picture is the appearance of secondary minima in the valence sector, which a single quadratic minimum cannot represent.

Extended reading notes

Core claim

Within the nCTEQ analysis framework, the thesis constructs a Bayesian posterior $P(\theta|D)\propto\exp(-\chi^2(\theta)/2)$ over the ten parameters of its nuclear-PDF ansatz and samples it with an adaptive Metropolis-Hastings algorithm. The posterior of the Pb-only fit is reported as markedly non-Gaussian, with the valence parameters $a_{dv3}$ and $a_{uv3}$ showing secondary minima beside the main minimum, so that a Hessian expansion around the single best point misrepresents the uncertainty. Comparing fits with and without lighter-nucleus data, the thesis finds that multi-nuclei constraints shift and narrow the quark distributions while the gluon stays comparatively stable. The central scientific claim is that the valence-sector posterior is genuinely multimodal and strongly correlated, making MCMC-based sampling, rather than the Hessian approximation, the appropriate uncertainty tool for nuclear PDF fits of this kind.

Load-bearing premise

The results assume the MCMC chains have converged and visited every relevant mode of the posterior; convergence is checked only empirically, and the paper's own trace plots show parameter drift and leftover hints of multimodality, so the reported secondary minima could in principle be sampling artifacts rather than true features of the likelihood.

Editorial extensions

If this is right

  • If the posterior is genuinely multimodal, Hessian-built 68% uncertainty bands understate the spread of valence-quark nuclear PDFs, and every observable computed from those bands inherits the error.
  • MCMC sampling becomes a necessary diagnostic, not an optional cross-check, for nPDF fits whose parameter spaces are strongly correlated and weakly constrained.
  • The Pb-only fit establishes a controlled single-nucleus benchmark: the first lead-PDF extraction constrained by lead data alone, against which future single-nucleus fits can be compared.
  • Including lighter nuclei changes the extracted quark distributions but leaves the gluon stable, implying that current datasets constrain nuclear gluons only weakly in the fitted kinematic range.
  • The comparison maps where the Hessian quadratic approximation is adequate (approximately Gaussian, single-mode regions) and where it is not, giving future analyses a picture of the unreliable regions.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same multimodality check should be applied to proton global fits, whose parameter spaces are similarly small and correlated; the companion proton MCMC study cited in the thesis is the natural place to look for valence-sector secondary minima there.
  • A nested-sampling or thermodynamic-integration run on the same likelihood would supply the Bayesian evidence and the relative posterior weight of each minimum, upgrading 'secondary minima exist' into a quantitative statement about how much posterior mass each mode carries.
  • If the gluon's stability under added $A$-dependence survives new data, it suggests the quark-sector shifts reflect genuinely improved quark constraints, while gluon insensitivity is a statement about current data coverage (mainly DIS and Drell-Yan) rather than about nuclear gluon physics itself.
Original abstract

Global analyses of nuclear parton distribution functions (nPDFs) have traditionally used the Hessian method for uncertainty estimation, relying on a local quadratic approximation of the chi-squared function. This Gaussian approximation can be inadequate when experimental constraints are limited and parameter correlations are strong, as it cannot fully capture asymmetric uncertainties, flat directions, or multiple minima. This thesis develops a Bayesian framework for nPDF determination using Markov Chain Monte Carlo (MCMC) techniques within the nCTEQ analysis framework. An adaptive Metropolis-Hastings algorithm enables direct sampling of the posterior distribution, revealing pronounced non-Gaussian behavior and secondary minima, particularly in the valence distributions. The analysis presents the first global single-nucleus determination of lead (Pb) PDFs based exclusively on lead data, complemented by a multi-nuclei analysis incorporating different nuclear targets through a standard analytic dependence on the mass number. Including lighter nuclei significantly constrains quark distributions and changes the extracted lead PDFs, while the gluon distribution remains comparatively stable. A comparison with Hessian-based results highlights the limitations of the Gaussian approximation in regions with strong nonlinearity or multimodality. Overall, this work establishes Bayesian MCMC inference as a flexible framework for nPDF determination, providing a more comprehensive characterization of parameter uncertainties and the underlying parameter space.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-10-28: 1 request, 6,329 input tokens.

What the abstract itself states, as probabilities: claims first 0.89 · states limitation 0.28 · human participants 0.01 · controlled comparison 0.47 · randomized assignment 0.01 · states quantitative result 0.07 · compares with prior results 0.97 · mentions code or data availability 0.03.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

Global analyses of nuclear parton distribution functions (nPDFs) have traditionally used the Hessian method for uncertainty estimation, relying on a local quadratic approximation of the chi-squared function. This Gaussian approximation can be inadequate when experimental constraints are limited and parameter correlations are strong, as it cannot fully capture asymmetric uncertainties, flat directions, or multiple minima. This thesis develops a Bayesian framework for nPDF determination using Markov Chain Monte Carlo (MCMC) techniques within the nCTEQ analysis framework. An adaptive Metropolis-Hastings algorithm enables direct sampling of the posterior distribution, revealing pronounced non-Gaussian behavior and secondary minima, particularly in the valence distributions. The analysis presents the first global single-nucleus determination of lead (Pb) PDFs based exclusively on lead data, complemented by a multi-nuclei analysis incorporating different nuclear targets through a standard analytic dependence on the mass number. Including lighter nuclei significantly constrains quark distributions and changes the extracted lead PDFs, while the gluon distribution remains comparatively stable. A comparison with Hessian-based results highlights the limitations of the Gaussian approximation in regions with strong nonlinearity or multimodality. Overall, this work establishes Bayesian MCMC inference as a flexible framework for nPDF determination, providing a more comprehensive characterization of parameter uncertainties and the underlying parameter space.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Global analyses of nuclear parton distribution functions (nPDFs) have traditionally used the Hessian method for uncertainty estimation, relying on a local quadratic approximation of the chi-squared function. background 1.00 0.05 0.02 asserted 1.00
s1 This Gaussian approximation can be inadequate when experimental constraints are limited and parameter correlations are strong, as it cannot fully capture asymmetric uncertainties, flat directions, or multiple minima. background 0.57 0.07 0.05 asserted 0.91
s2 This thesis develops a Bayesian framework for nPDF determination using Markov Chain Monte Carlo (MCMC) techniques within the nCTEQ analysis framework. contribution 0.88 0.92 0.08 asserted 1.00
s3 An adaptive Metropolis-Hastings algorithm enables direct sampling of the posterior distribution, revealing pronounced non-Gaussian behavior and secondary minima, particularly in the valence distributions. result 0.71 0.95 0.11 asserted 1.00
s4 The analysis presents the first global single-nucleus determination of lead (Pb) PDFs based exclusively on lead data, complemented by a multi-nuclei analysis incorporating different nuclear targets through a standard analytic dependence on the mass number. contribution 0.89 0.95 0.08 asserted 1.00
s5 Including lighter nuclei significantly constrains quark distributions and changes the extracted lead PDFs, while the gluon distribution remains comparatively stable. result 0.99 0.93 0.07 asserted 1.00
s6 A comparison with Hessian-based results highlights the limitations of the Gaussian approximation in regions with strong nonlinearity or multimodality. interpretation 0.58 0.84 0.09 asserted 0.99
s7 Overall, this work establishes Bayesian MCMC inference as a flexible framework for nPDF determination, providing a more comprehensive characterization of parameter uncertainties and the underlying parameter space. interpretation 0.80 0.87 0.89 asserted 1.00

What kind of paper the classifier reports this is: method or tool 0.92 · evidence from analysis of existing data 0.75.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 2. 1: Classification of elementary particles according to the Standard Model of particle
Figure 3. 1: Nuclear correction ratio F Fe 2 /F D 2 , the figure taken from [126]. 1. Shadowing (x ≲ 0.1): At very low x, the nuclear structure function is suppressed relative to that of free nucleons, yielding RA(x) < 1. This suppression, known as shadowing [29, 30, 31], is believed to arise from coherent multiple scattering and the overlap of partons from different nucleons. In this regime, gluon recombination and coherence effects reduce the effective parton densities, especially for gluons and sea 4This expectation holds exactly only for isoscalar nuclei, where the number of protons equals
Figure 6. 2: Scan of the χ 2 function along the parameter a uv 3 in the Pb-only analysis. The parameter a uv 3 is varied while all other parameters are fixed at their values corresponding to the MCMC global minimum of the full fit. The profile exhibits a nearly constant χ 2 over a broad interval, indicating weak sensitivity of the likelihood to this parameter, together with a sharp, localized increase around a uv 3 = 0 that strongly disfavors this region. volume of the parameter space, the Markov chain can efficiently traverse it and explore the extended, nearly degenerate likelihood regi
Figure 7. 13: Pb-only fit: Comparison of 68% uncertainty bands for the Pb valence (uv, dv), light-sea (u¯+d¯), and gluon (g) distributions obtained using three MCMC-based estimators: cumulative-χ 2 (red), Percentile (blue), and MCSE (orange). For each flavor, the upper panel shows the absolute distributions, while the lower panel displays ratios to the cumulative-χ 2 best-fit prediction. The differences are more pronounced in the uncertainty bands. While the cumulative-χ 2 method consistently yields the widest (most conservative) uncertainties, the Percentile and MCSE bands overlap very closely

Showing 4 of the paper's 7 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:2610.37254 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.