REVIEW 3 major objections 3 minor 27 references
Absolute intensity normalisation of powder neutron scattering data
T0 review · 3 major / 3 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read A simple table of conversion factors turns the Rietveld scale factor from any standard refinement into an absolute intensity scale for powder neutron data, so no extra calibration measurement is needed.
desk verdict A clean practical derivation of Rietveld-scale-to-absolute conversion factors, with an accuracy claim that the worked example cannot actually support. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the absolute normalisation factor $s$, defined by $I_{\mathrm{expt}} = s I_{\mathrm{abs}}$, which converts measured counts into barn sr⁻¹ atom⁻¹. After Rietveld refinement, the ratio $I_{\mathrm{expt}}/I_{\mathrm{abs}}$ is rewritten as $\mathrm{RietveldScale}\times I_0/I_{\mathrm{abs}}$, where $I_0$ is the bare nuclear Bragg profile programmed into the refinement code. Comparing the program's $I_0$ with the physically absolute $I_{\mathrm{abs}}$ and changing variables from $Q$ to $2\theta$ or time-of-flight $t$ produces the factors in Table I; the GSAS-II variants differ because GSAS-II divides by unit-cell volume and uses centidegrees.
What would settle it
Take a polycrystalline sample whose absolute intensity has been independently calibrated with a vanadium standard, apply Table I to a fresh refinement, and compare the normalised intensities against that independent calibration; a systematic offset beyond the statistical error of the Rietveld scale factor would show that the program's internal intensity differs from the assumed bare Bragg profile.
Extended reading notes
Core claim
The central claim is that the Rietveld scale factor can be converted into an absolute normalisation factor $s$, defined by $I_{\mathrm{expt}} = s I_{\mathrm{abs}}$, where $I_{\mathrm{abs}}$ is the true nuclear Bragg intensity in units of barn sr⁻¹ atom⁻¹. Because the refined Rietveld scale factor equals $I_{\mathrm{expt}}/I_0$, one just needs to know how the program's internal $I_0$ differs from $I_{\mathrm{abs}}$. The paper works this out for FullProf and GSAS-II, obtaining factors such as $s^{\mathrm{Fullprof}}_{2\theta} = \mathrm{FullprofScale}\times 2\pi^2 N V/(45\lambda^3)$ for constant-wavelength data and analogous factors for time-of-flight data and GSAS-II; Table I collects all four. The derivation is tested on Dy$_3$Mg$_2$Sb$_3$O$_{14}$, where constant-wavelength and time-of-flight data, normalised independently, give nearly identical absolute magnetic diffuse scattering. The method is deliberately limited to polycrystalline samples with a well-refined average structure and is not intended for liquids, amorphous materials, or poorly crystallised samples.
Load-bearing premise
The only load-bearing assumption is that the Rietveld program's internal intensity calculation is exactly the bare nuclear Bragg profile described in the paper, with absorption and other overall intensity corrections already applied to the data before refinement, so that the refined scale factor really equals the ratio of measured to calculated intensity.
Editorial extensions
If this is right
- Any polycrystalline sample with a well-refined crystal structure can be placed on an absolute intensity scale with no additional measurement beyond the Rietveld refinement itself, using the appropriate factor from Table I.
- The normalised data are directly comparable across instruments and even across techniques, as the constant-wavelength and time-of-flight results for Dy3Mg2Sb3O14 show.
- Magnetic diffuse scattering is expressed in absolute units, so magnetic spectral weights and ordered moments can be quantified without the systematic uncertainty typical of vanadium-based normalisation.
- The same factors apply to energy-integrated spectroscopy data from direct-geometry instruments, allowing phonon and magnetic-excitation spectra to be normalised absolutely.
- The method fails for liquids, amorphous materials, or poorly crystallised samples where no well-refined average structure exists.
Reading between the lines
- Not stated in the paper, but the method would allow archived powder datasets to be re-normalised retrospectively if their Rietveld scale factors and cell parameters survive, since no new measurement is needed.
- Not stated in the paper, but a magnetic-structure refinement could use the same logic to put magnetic Bragg intensities on an absolute scale and extract ordered moment magnitudes without an external standard.
- Not stated in the paper, but the approach offers a cross-instrument consistency check: two instruments with correct program conventions should give the same factor $s$ for the same sample.
- Not stated in the paper, but comparing $s$ across many samples on one instrument could reveal whether a program version silently changes its internal intensity convention.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a method for placing powder neutron diffraction data on an absolute intensity scale by using the intensity scale factor obtained from a Rietveld refinement. The derivation relates the refinement scale factor to the ratio of the measured profile to the internally calculated profile, and Table I gives conversion formulae for FullProf and GSAS-II for both constant-wavelength and time-of-flight data. A worked example for Dy3Mg2Sb3O14 compares constant-wavelength data from DCS and time-of-flight data from GEM after normalization and reports agreement between the two instruments. The authors claim the accuracy is likely significantly better than the ~20% typical of other methods.
Significance. The method is attractive because it avoids the separate vanadium calibration measurement and uses the sample's own nuclear Bragg scattering as an internal standard. The derivation is transparent and the resulting Table I is directly usable by practitioners. The worked example demonstrates that the two instruments agree after the proposed normalization, which supports internal consistency. However, the absolute scale of Table I is not independently validated, and the method depends sensitively on the exact internal intensity conventions of each Rietveld program, a point the paper itself acknowledges for one data format. If the absolute scale is confirmed, the paper would be a useful practical contribution; as it stands, the evidence does not yet support the claimed accuracy.
major comments (3)
- [Worked example, Fig. 1(c)] The agreement between the DCS and GEM datasets in Fig. 1(c) validates that the two conversion factors used are internally consistent, but it does not validate their absolute scale. Since both factors are derived from the same assumption about the Rietveld program's internal intensity I0, a common multiplicative error in Table I would cancel when the two independently normalized datasets are compared. To support the statement that the accuracy is 'likely to be significantly better than the ~20% systematic uncertainty' of other methods, the authors should compare at least one normalized dataset against an independent absolute calibration, such as a vanadium-standard normalization or a well-characterized standard sample with a known absolute cross-section.
- [§2, Eq. (3) and Table I footnote] The derivation assumes that the Rietveld scale factor is exactly the ratio of the measured intensity to the program-internal I0 of Eq. (6). The footnote in Table I shows that this assumption is fragile: for some TOF data formats FullProf applies an undocumented internal factor of 1000, so the conversion factor must be rescaled by 0.001. Because other data formats or program versions could contain similar undocumented multipliers, the current presentation leaves no way for a user to know whether a given format requires such a correction. The manuscript should either provide a procedure to determine the effective internal scaling for a specific data format and software version (for example, using a standard sample), or explicitly state that the table is valid only for formats and versions for which the program documentation guarantees no additional overall multiplier.
- [§2, 'Time-of-flight diffraction', Eq. (19)] The derivation leading to Eq. (19) relies on the identification of DIFC (or dtt1) with the instrument parameter in Eq. (20), and on the specific form of the FullProf intensity expression. The paper states that the GSAS-II results differ because its intensity divides by V and uses centidegrees, but the derivation of the GSAS-II factor, especially the appearance of the constant 4500 in Eq. (13) and the corresponding TOF factor in Table I, is not shown. Given that the entire method depends on these constants, the authors should provide a derivation or a precise citation to the relevant program-manual equations for each entry in Table I, so that the factors can be independently checked.
minor comments (3)
- [Introduction] In the sentence 'since it allows for quantitative determination the magnitude of ordered magnetic moments', the word 'of' is missing between 'determination' and 'the magnitude'.
- [Discussion and Conclusions] The sentence 'The main limitation of this approach it that' should read 'The main limitation of this approach is that'.
- [Worked example, Fig. 1 caption] The caption of Fig. 1 would benefit from stating explicitly that the vertical axis in panels (a) and (b) is in arbitrary linear units, since the unit of the Rietveld scale factor is not stated in the text.
Circularity Check
No significant circularity: the normalisation factors are derived from published scattering equations and program documentation; the Rietveld scale factor is an input, not a predicted output.
full rationale
The paper's central chain is Eqs. (1)-(3), which define the absolute normalisation factor s as RietveldScale times I0/Iabs. I0 is taken from the FullProf/GSAS-II program definitions (Eq. (6) and following, with cited manuals [14] and [20]), and Iabs is computed from the standard polycrystalline Bragg intensity formula Eq. (4) (Squires [19]) by analytic change of variable. Table I then follows by algebra; no step fits a parameter to the target data and re-presents it as a prediction. The Rietveld scale factor is a fitted input used in the normalisation, not a quantity the paper claims to predict. The worked example compares two independently measured datasets and is a demonstration, not the derivation; although a common multiplicative error in Table I would not be detected by this comparison, that is a validation limitation, not circularity. Self-citations ([18], [25], [26]) are data or application references and are not load-bearing. The only caveat is program-convention dependence (e.g., the Table I footnote about FullProf internally multiplying some TOF data by 1000), which is an external-validity issue rather than a circular-reasoning issue.
Assumptions & free parameters
assumptions (4)
- standard math The powder-averaged nuclear Bragg intensity is given by Eq. (4) in barn sr^-1 atom^-1, Iabs(Q) = (2pi^2/NV) sum_G m_G |F_G|^2 / G^2 R_Q(Q-G).
- domain assumption Fullprof's internal profile intensity has the form I0(x) = sum_G m_G |F_G|^2 L_x R_x(x-x_G), with no additional overall intensity factor; GSAS-II differs by dividing by V and using centidegrees.
- domain assumption Experimental data are corrected for absorption and other overall intensity effects before Rietveld refinement, so the refined scale factor only multiplies the pure nuclear Bragg profile.
- domain assumption The structure factor F_G uses accurate nuclear scattering lengths and Debye-Waller factors for the sample.
Cite this review
Pith. "Pith review of Absolute intensity normalisation of powder neutron scattering data." pith.science (2026). https://pith.science/paper/RMDZNWUY
@misc{pith2026241118547,
author = {Pith},
title = {Pith review of: Absolute intensity normalisation of powder neutron scattering data},
year = {2026},
howpublished = {\url{https://pith.science/paper/RMDZNWUY}},
note = {Machine review of arXiv:2411.18547}
}
read the original abstract
An important property of neutron scattering data is that they can be normalised in absolute intensity units. In practice, however, such normalisation is often not performed, since it can be time-consuming and subject to systematic uncertainties. Here, a straightforward approach is presented for absolute intensity normalisation of neutron scattering data from polycrystalline samples. This approach uses the intensity scale factor obtained from a Rietveld refinement to normalise the data to the nuclear Bragg profile of the sample. Factors to convert the Rietveld scale factor into an absolute normalisation factor are tabulated for constant-wavelength and time-of-flight data refined using the popular programs Fullprof and GSAS-II. An example of the application of this method to experimental data is presented. Advantages, disadvantages, and extensions of this approach to spectroscopic data are discussed.
Figures
Reference graph
Works this paper leans on
-
[1]
B. H. Toby, T. Egami, Acta Crystallogr . A48, 336 (1992)
work page 1992
-
[2]
S. J. L. Billinge, M. G. Kanatzidis, Chem. Commun. pp. 749– 760 (2004)
work page 2004
-
[3]
D. A. Keen, R. L. McGreevy, Nature 344, 423 (1990)
work page 1990
-
[4]
M. G. Tucker, D. A. Keen, M. T. Dove, A. L. Goodwin, Q. Hui, J. Phys.: Condens. Matter 19, 335218 (2007)
work page 2007
-
[5]
K. R. A. Ziebeck, P. J. Brown, J. Phys. F: Met. Phys. 10, 2015 (1980)
work page 1980
-
[6]
O. Steinsvoll, C. F. Majkrzak, G. Shirane, J. Wicksted, Phys. Rev. B 30, 2377 (1984)
work page 1984
- [7]
-
[8]
Mourigal, et al., Nature Physics 9, 435 (2013)
M. Mourigal, et al., Nature Physics 9, 435 (2013)
work page 2013
Show all 27 references
-
[9]
G. D. Wignall, F. S. Bates, J. Appl. Crystallogr .20, 28 (1987)
1987
-
[10]
G. Xu, Z. Xu, J. M. Tranquada, Rev. Sci. Instrum. 84, 083906 (2013)
2013
-
[11]
A. K. Soper, GudrunN and GudrunX: programs for correcting raw neutron and X-ray diffraction data to differential scattering cross section (Science & Technology Facilities Council Swin- don, UK, 2011)
2011
-
[12]
H. M. Rietveld, J. Appl. Crystallogr .2, 65 (1969)
1969
-
[13]
Rodr ´ıguez-Carvajal, Physica B 192, 55 (1993)
J. Rodr ´ıguez-Carvajal, Physica B 192, 55 (1993)
1993
-
[14]
Rodr ´ıguez-Carvajal, FullProf Manual , Laboratoire L ´eon Brillouin, CEA/Saclay, 91191 Gif sur Yvette Cedex, France (2001)
J. Rodr ´ıguez-Carvajal, FullProf Manual , Laboratoire L ´eon Brillouin, CEA/Saclay, 91191 Gif sur Yvette Cedex, France (2001)
2001
-
[15]
B. H. Toby, R. B. V on Dreele, J. Appl. Crystallogr . 46, 544 (2013)
2013
-
[16]
A. A. Coelho, J. Appl. Crystallogr .51, 210 (2018)
2018
-
[17]
V . F. Sears,Neutron News 3, 26 (1992)
1992
-
[18]
J. A. M. Paddison, et al., Nat. Commun. 7, 13842 (2016)
2016
-
[19]
G. L. Squires, Introduction to the Theory of Thermal Neutron Scattering (Cambridge University Press, Cambridge, 1978)
1978
-
[20]
A. C. Larson, R. B. V on Dreele, GSAS Manual , Los Alamos National Laboratory (1985)
1985
-
[21]
Zhang, J
Y . Zhang, J. Liu, M. G. Tucker, Acta Crystallogr . A 79, 20 (2023). 5
2023
-
[22]
J. R. D. Copley, J. C. Cook, Chem. Phys. 292, 477 (2003)
2003
-
[23]
A. C. Hannon, Nucl. Instrum. Methods Phys. Res. A 551, 88 (2005)
2005
-
[24]
P. F. Peterson, D. Olds, M. T. McDonnell, K. Page, J. Appl. Crystallogr .54, 317 (2021)
2021
-
[25]
J. A. M. Paddison, et al., Cell Rep. Phys. Sci. 5 (2024)
2024
-
[26]
J. A. M. Paddison, et al., npj Quantum Mater .9, 48 (2024)
2024
-
[1000]
In such cases, this conversion factor should be replaced by FullprofScale × 0.001 × 4πNV sin θ dtt1 when applied to the original data. where the geometrical factors have been combined into the Lorentz factor for time-of-flight diffraction [21], Lt ≡ d4 sin θ , (18) where d = 2...
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.