Beyond Gaussian Statistics in Polymer Melts: Statistical Masking of Persistent Local Constraints

Jos\'e A. Martins

REVIEW 2 major objections 1 minor 14 references

Beyond Gaussian Statistics in Polymer Melts: Statistical Masking of Persistent Local Constraints

T0 review · 2 major / 1 minor · reviewed 2026-06-29 · grok-4.3

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Polymer chains recover Gaussian end-to-end statistics through statistical masking by accumulating random segments, not by erasing persistent local alignments.

desk verdict The paper claims Gaussian recovery in long polymer chains comes from statistical masking by RCS segments rather than erasure of persistent ACS domains at ~35%, backed by MD fits to q-Gaussians and direct entropy ratios. read the letter →

arxiv 2605.25989 v2 pith:OVY5KS4W submitted 2026-05-25 cond-mat.soft cond-mat.stat-mechphysics.bio-ph

José A. Martins This is my paper · ORCID

classification cond-mat.softcond-mat.stat-mechphysics.bio-ph

keywords polymermeltsGaussianstatisticsq-GaussianconformationalheterogeneitypolyethylenestatisticalmaskingKuhnscaleend-to-enddistance

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Short polymer chains show clear non-Gaussian end-to-end distance statistics because of conformational heterogeneities that exist at the Kuhn scale. These heterogeneities form a mosaic of extended aligned segments and coiled segments that persists in chains of all lengths above a critical mass. Longer chains accumulate more independent random conformational sequences that progressively obscure the non-Gaussian signatures from the aligned segments. The q-Gaussian distribution captures this process, with its entropic index q rising toward 1 as the masking effect strengthens. The ratio of Tsallis entropy to Boltzmann-Gibbs entropy falls toward unity, confirming the approach to extensive statistics without loss of the underlying local constraints.

What carries the argument

The q-Gaussian distribution whose entropic index q serves as a heterogeneity index that quantifies how the accumulation of independent random conformational sequences masks the non-Gaussian effects of persistent aligned chain segments.

What would settle it

Simulation or experimental data on very long chains that either show aligned segments dropping well below 35 percent or produce end-to-end distributions that deviate from the q-Gaussian form while q equals 1.

Watch

Extended reading notes

Core claim

The paper shows that Gaussian recovery in long polyethylene chains occurs because independent random conformational sequences accumulate and mask the non-Gaussian signatures of persistent aligned chain segments, which remain at roughly 35 percent even in chains up to C500. Both unentangled and entangled chains follow q-Gaussian distributions whose entropic index q increases from 0.67 to 0.99 with chain length, tracking the growth of the masking segments. The ratio of Tsallis to Boltzmann-Gibbs entropy, computed directly from the data, decreases from 1.80 to 1.03 over the same range, establishing q as a quantitative heterogeneity index.

Load-bearing premise

Conformational heterogeneity at the Kuhn scale, with aligned segments remaining near 35 percent, persists across all chain lengths and is accurately captured by the q-Gaussian form without those heterogeneities being erased.

Editorial extensions

If this is right

End-to-end distance distributions remain accurately described by a q-Gaussian for both unentangled and entangled chains.
Aligned chain segments persist at approximately 35 percent in all chains above the critical mass.
The entropic index q increases systematically with chain length as more random conformational sequences accumulate.
The ratio of Tsallis to Boltzmann-Gibbs entropy decreases toward 1, confirming the statistical masking process.

Reading between the lines

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

Local conformational constraints may continue to influence dynamic properties such as relaxation times even in long chains where end-to-end statistics appear Gaussian.
Similar masking mechanisms could apply to other systems that exhibit persistent heterogeneities at a fixed length scale while global statistics converge.
Direct visualization or labeling of aligned versus random segments in longer chains would provide an independent test of the masking fraction.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

read the letter

The main thing to know is that this work argues the shift to Gaussian end-to-end statistics in polyethylene melts is a masking effect: random conformational sequences accumulate and hide the non-Gaussian signatures from aligned chain segments that stay present at roughly 35% even in long chains. Short chains lack those RCS pieces, so their distributions deviate more, and q in the q-Gaussian fit rises toward 1 as RCS build up. They also compute the Tsallis-to-Boltzmann entropy ratio straight from the binned data and show it approaches 1.

The simulations track segment types at the Kuhn scale across C50 to C500, which gives a concrete way to link local structure to the global statistics. Treating q as a heterogeneity index is a clean move, and computing the entropy ratio without extra fitting avoids one layer of circularity. That part of the evidence is straightforward.

The soft spot is the segment classification itself. The constant ACS fraction is the load-bearing observation for the masking claim, but the criteria for calling a segment aligned or extended could shift with chain length if global constraints change local bond statistics even slightly. The abstract does not show tests of threshold sensitivity or order-parameter independence from N, so the persistence result needs that check before the interpretation is secure. If the classification holds under variation, the masking story strengthens; if not, it weakens.

This is for people in polymer physics and soft matter who work on chain statistics and local order parameters. It engages the existing literature on non-Gaussian deviations and supplies simulation numbers rather than just a new model. The thinking is internally consistent once the classification is accepted.

I would bring it to a reading group to discuss the classification method. I would not cite it in the next year because the result is still tied to one set of polyethylene runs. It deserves peer review because the simulation evidence and the alternative mechanism are specific enough to get useful referee comments, even with the classification question open.

Referee Report

2 major / 1 minor

Summary. The manuscript uses atomistic MD simulations of polyethylene melts to argue that recovery of Gaussian end-to-end distance statistics in long chains occurs via statistical masking by accumulating random conformational sequences (RCS), rather than erasure of persistent aligned chain segments (ACS) at the Kuhn scale. ACS domains remain at ~35% fraction for all chains above a critical length; end-to-end distributions are fit to q-Gaussians with q rising from 0.67 (C50) to 0.99 (C500), and the directly computed Tsallis-to-Boltzmann entropy ratio falls from 1.80 to 1.03, positioning q as a heterogeneity index.

Significance. If the masking interpretation is substantiated, the work offers a mechanistic distinction between homogenization and statistical obscuration in polymer statistics, with potential implications for entangled dynamics and non-extensive entropy descriptions. The direct entropy-ratio computation from raw data and the systematic q trend constitute concrete, falsifiable elements.

major comments (2)

[Abstract] Abstract: The central claim that ACS domains persist at a constant ~35% fraction (distinguishing masking from erasure) rests on Kuhn-scale classification of ACS/RCS/CE whose independence from chain length is not demonstrated; if the order parameters (local alignment, extension, or relaxation times) are sensitive to global constraints that strengthen with length, the measured fraction could remain stable while underlying heterogeneities are modified.
[Results/Methods] Results/Methods: Full simulation protocols, error bars on the reported q values and ACS fractions, raw distribution data, and robustness checks on the q-Gaussian fits (e.g., alternative binning or maximum-likelihood procedures) are absent, preventing verification that the q trend and entropy ratios are not artifacts of fitting or classification choices derived from the same trajectories.

minor comments (1)

[Abstract] Notation: The definition of the entropic index q and its relation to the heterogeneity index should be stated explicitly with the functional form of the q-Gaussian used for fitting.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and constructive feedback on our manuscript. We address each of the major comments below and will incorporate the necessary revisions to strengthen the presentation.

read point-by-point responses

Referee: [Abstract] Abstract: The central claim that ACS domains persist at a constant ~35% fraction (distinguishing masking from erasure) rests on Kuhn-scale classification of ACS/RCS/CE whose independence from chain length is not demonstrated; if the order parameters (local alignment, extension, or relaxation times) are sensitive to global constraints that strengthen with length, the measured fraction could remain stable while underlying heterogeneities are modified.

Authors: We thank the referee for highlighting this important point. In the manuscript, the ACS fraction is reported as approximately 35% for all chains above the critical length based on our simulations across C50 to C500. To explicitly demonstrate the independence from chain length, we will add a new figure in the revised manuscript showing the ACS fraction versus chain length, with error bars from multiple trajectories. Regarding the potential sensitivity of order parameters to global constraints, our definitions are strictly local (Kuhn-scale alignment and relaxation), and we have checked that they remain consistent. However, we will include additional discussion and perhaps a correlation analysis to address this concern. revision: yes
Referee: [Results/Methods] Results/Methods: Full simulation protocols, error bars on the reported q values and ACS fractions, raw distribution data, and robustness checks on the q-Gaussian fits (e.g., alternative binning or maximum-likelihood procedures) are absent, preventing verification that the q trend and entropy ratios are not artifacts of fitting or classification choices derived from the same trajectories.

Authors: We agree with the referee that these details are crucial for reproducibility and verification. The full simulation protocols were summarized due to length constraints but will be expanded in the Methods section of the revised manuscript. We will also include error bars on q values and ACS fractions (computed from independent runs), provide raw end-to-end distance distribution data in supplementary materials, and add robustness checks for the q-Gaussian fits using alternative methods such as maximum likelihood estimation and different binning schemes. These additions will confirm that the trends are not artifacts. revision: yes

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity; claims rest on independent simulation observables

full rationale

The paper measures end-to-end distributions, fits q-Gaussian parameters, computes S_q/S_1 directly from binned data without fitting, and separately classifies local ACS/RCS/CE domains at the Kuhn scale. The masking interpretation follows from the observed persistence of ACS fraction (~35%) and the correlation between q and RCS accumulation, but neither the fitted q nor the domain fractions are defined in terms of each other or the final claim. No equation reduces to its input by construction, no parameter is renamed as a prediction, and no self-citation chain bears the central result. The derivation remains self-contained against the simulation data.

Assumptions & free parameters 1 free parameters · 1 assumptions · 2 invented entities

The central claim rests on fitted q parameter, assumption that q-Gaussian applies, and new segment classifications (ACS/RCS) derived from the simulations themselves with no external validation.

free parameters (1)

entropic index q = 0.67 to 0.99
Fitted parameter to q-Gaussian distributions for different chain lengths C50 to C500

assumptions (1)

domain assumption q-Gaussian function accurately describes end-to-end distance distributions
Invoked to quantify non-extensive statistics and heterogeneity

invented entities (2)

ACS (slow-relaxing extended aligned chain segments)
purpose: Represent persistent local constraints at Kuhn scale
Postulated from simulation conformational analysis with no independent evidence outside the paper
RCS (random conformational sequences)
purpose: Represent coiled segments that accumulate to mask ACS
Defined from simulation data to explain masking effect

how reviews work

0 comments

Cite this review

Pith. "Pith review of Beyond Gaussian Statistics in Polymer Melts: Statistical Masking of Persistent Local Constraints." pith.science (2026). https://pith.science/paper/OVY5KS4W

@misc{pith2026260525989,
  author       = {Pith},
  title        = {Pith review of: Beyond Gaussian Statistics in Polymer Melts: Statistical Masking of Persistent Local Constraints},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OVY5KS4W}},
  note         = {Machine review of arXiv:2605.25989}
}

abstract

Short polymer chains exhibit clear deviations from Gaussian end-to-end distance statistics, yet the molecular mechanism by which Gaussian behavior is recovered in long chains remains unestablished. Atomistic molecular dynamics simulations of polyethylene melts reveal that conformational heterogeneity persists at the Kuhn scale across all chain lengths, consisting of a mosaic of slow-relaxing, extended aligned chain segments (ACS) and coiled segments -- random conformational sequences (RCS) and chain ends (CE). We show that the end-to-end distance distributions for both unentangled and entangled chains are accurately described by a $q$-Gaussian function, with the entropic index $q$ increasing systematically from $0.67$ (C50) to $0.99$ (C500). This evolution tracks the emergence and accumulation of RCS segments, which are absent in short chains, establishing $q$ as a quantitative ``heterogeneity index''. The $q < 1$ values are a signature of non-extensive statistics, with the ratio of Tsallis to Boltzmann-Gibbs entropy ($S_q/S_1$), computed directly from simulation data without fitting, decreasing from $1.80$ (C50) to $1.03$ (C500). Crucially, we demonstrate that Gaussian recovery does not result from the erasure of Kuhn-scale heterogeneities, as ACS domains persist in all chain lengths above the critical mass ($\approx 35\%$). Instead, the transition to Gaussian statistics is a statistical masking effect, where the accumulation of independent RCS segments progressively obscures the non-Gaussian signatures of the persistent ACS domains.

Figures

Figures reproduced from arXiv: 2605.25989 by the authors.

**Figure 1.** FIG. 1. Data histogram for the end-to-end distances of the different polymer systems and fits with a Gaussian, Eq. (1), and [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗

**Figure 2.** FIG. 2. Linearity test in [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗

**Figure 3.** FIG. 3. Tail concentration of the last data quartile ( [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗

Figures from the paper (5 more)

**Figure 4.** Figure 4: FIG. 4. Entropy non-extensivity quantified by the ratio [PITH_FULL_IMAGE:figures/full_fig_p011_4.png]

**Figure 5.** Figure 5: FIG. 5. Segment length distributions in terms of united atoms (UA) for ACS (red) and RCS (blue) in all four systems. C50 has no RCS [PITH_FULL_IMAGE:figures/full_fig_p012_5.png]

**Figure 6.** Figure 6: FIG. 6. End-to-end distance distributions [PITH_FULL_IMAGE:figures/full_fig_p012_6.png]

**Figure 7.** Figure 7: FIG. 7. Correlation between the average end-to-end distance of the chain and the number of ACS and RCS segments present in each chain. [PITH_FULL_IMAGE:figures/full_fig_p013_7.png]

**Figure 8.** Figure 8: FIG. 8. Time-weighted mean entropy ratio [PITH_FULL_IMAGE:figures/full_fig_p014_8.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

14 extracted references

[1]

42 This constraint defines the domain of the function forq<1

Theq-exponential function In nonextensive statistical mechanics Tsallis introduced a generalized exponential function, theq-exponential [exp q(u)] as a generalization of [exp(u)].36,38,39 It is defined for a single argumentuas expq(u) = [1+ (1−q)u] 1 1−q + ,(4) where the operator[·] + ensures that the argument remains pos- itive, guaranteeing real values ...
[2]

Forq<1 the distribution has compact sup- port, while forq>1 it exhibits long tails (Appendix A and supplementary material, Sec

Radialq-Gaussian probability density in 3D The radialq-Gaussian probability density (with units of L−1) for the scalar end-to-end distanceRis defined by Pq(R) =C qR2 expq −βqR2 ,(8) whereC q is the normalization constant (unitsL −3) ensuringR ∞ 0 Pq(R)dR=1. Forq<1 the distribution has compact sup- port, while forq>1 it exhibits long tails (Appendix A and ...
[3]

As shown in the Appendix (Sec

Second moment Since the moments of interest remain finite for the range ofqvalues relevant to the polymer melts studied in this work (q<5/3 in three dimensions), we employ the standard expec- tation value ⟨Rk⟩q = Z ∞ 0 RkPq(R)dR,(9) avoiding the use of escort distributions. As shown in the Appendix (Sec. A 4) and supplementary material (Sec. S1.2), for bo...
[4]

Inserting this definition in Eq

Mapping of the q-Gaussian second moment onto an equivalent random walk To relate thisq-Gaussian mean-square end-to-end distance to an effective random walk, we define βq = 3 7−5q 1 Nq l2q ,(11) that equalsβ g whenq=1. Inserting this definition in Eq. (10), the second moment may be expressed as a random walk ofNq steps, each with lengthl q ⟨R2⟩q =N ql2 q ....
[5]

renor- malized

Theq-logarithm and linearization In the Gaussian case (q=1), as shown in Eq. (3), the plot of the logarithm of the volumetric densityρ g(R)∝P g(R)/R2 versusR 2 is linear. For theq-Gaussian, this linearity is pre- served only when utilizing theq-logarithm. Theq-logarithm is the inverse function of theq-exponential [ln q[expq(u)] =u] and is defined fory>0 a...
[6]

Boltzmann–Gibbs entropy (S1) The classical extensive entropy is defined using the volu- metric probability densityρ(R)(units ofL −3) integrated over the 3D volume. Expressed in terms of the radial probability densityP(R)(units ofL −1), it is: S1 =−k B Z ∞ 0 P(R)ln ρ(R) ρ0 dR+const.,(16) whereρ(R) =P(R)/(4πR 2)andρ 0 =1 Å −3 is a reference density introduc...
[7]

experimental

Tsallisq-entropy (S q) To capture possible non-extensive behavior due the ex- istence of different types and relative populations of Kuhn segments, we use the generalized Tsallis entropy. 36,43 In the discrete implementation, the dimensionless probabilityp k is raised to the powerq Sq =k B 1− ∑k pq k q−1 .(18) Asq→1,S q converges toS 1. A value ofq̸=1 sig...
[8]

Every chaini(i=1,

Space and time (Sp&Ti) averaging In this conventional approach (typical of steady-state molecular dynamics sampling), all sampled chain configura- tions in space and time are combined into a single ensemble. Every chaini(i=1, . . . ,N ch), at every framet, contributes an end-to-end distanceR i,t to the overall distribution. The space- time probability dis...

Show all 14 references

[9]

variables

Frame–by–frame (FrbyFr) averaging In this approach, each simulation frame is treated as an independent, quasi-equilibrium instantaneous configuration. The molecular dynamics thermostat maintains thermal equi- librium even as conformational variables fluctuate. In each frame th...
[10]

Forq<1 the distribution has compact support (supplementary material, Sec

Normalization of the partition functionZ q The radialq-Gaussian probability density in three dimen- sions may be written as Pq(R) = 4πR2 Zq expq(−βqR2),(A1) whereZ q ensures R ∞ 0 Pq(R)dR=1. Forq<1 the distribution has compact support (supplementary material, Sec. S1.1) and th...
[11]

(A4) and (A5) in the definition ofP q(R), Eq

Normalized 3D Radialq-Gaussian Replacing Eqs. (A4) and (A5) in the definition ofP q(R), Eq. (A1), yields the normalizedq-Gaussian probability den- sity 17 P(q<1) q (R) =4πR 2 (1−q) 3/2 βq π 3/2 Γ 5 2 + 1 1−q Γ 1+ 1 1−q 1−(1−q)β qR2 1 1−q (A6a) P(q>1) q (R) =4πR 2 (q−1) 3/2 βq ...
[12]

To verify that the prefactor in Eqs

Gaussian limitq→1 Forq→1 one has exp q(−x)→exp(−x). To verify that the prefactor in Eqs. (A6a) and (A6b) converge to the same Gaussian distribution, we defineν=1/(1−q)and use the asymptotic ratio of Gamma functions,Γ(ν+a)/Γ(ν+b)≈ νa−b forν→∞. For theq<1casea=5/2,b=1 and (1−q) ...
[13]

The energy landscapes and motions of proteins,

Second moment forq<1 The second moment is defined as ⟨R2⟩q = I2 I0 = R ∞ 0 4πR4[1−(1−q)β qR2 ]1/(1−q)dRR ∞ 0 4πR2[1−(1−q)β qR2 ]1/(1−q)dR .(A7) TheI 0 was evaluated above and is given by Eq. (A4) or (A5), depending on theqvalue. To evaluateI 2 we follow the same steps as in Eq...

2003
[14]

Gromacs 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation,

pp. 79–98. 44B. Hess, C. Kutzner, D. Spoel, and E. Lindahl, “Gromacs 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation,” J. Chem. Theory Comput4, 435–447 (2008). 45B. Hess, H. Bekker, H. Berendsen, and J. Fraaije, “Lincs: A linear con- strain...

2008

Pith tools

Reviewed June 29, 2026 · model on record in the stance chip above.

[1] [1]

42 This constraint defines the domain of the function forq<1

Theq-exponential function In nonextensive statistical mechanics Tsallis introduced a generalized exponential function, theq-exponential [exp q(u)] as a generalization of [exp(u)].36,38,39 It is defined for a single argumentuas expq(u) = [1+ (1−q)u] 1 1−q + ,(4) where the operator[·] + ensures that the argument remains pos- itive, guaranteeing real values ...

[2] [2]

Forq<1 the distribution has compact sup- port, while forq>1 it exhibits long tails (Appendix A and supplementary material, Sec

Radialq-Gaussian probability density in 3D The radialq-Gaussian probability density (with units of L−1) for the scalar end-to-end distanceRis defined by Pq(R) =C qR2 expq −βqR2 ,(8) whereC q is the normalization constant (unitsL −3) ensuringR ∞ 0 Pq(R)dR=1. Forq<1 the distribution has compact sup- port, while forq>1 it exhibits long tails (Appendix A and ...

[3] [3]

As shown in the Appendix (Sec

Second moment Since the moments of interest remain finite for the range ofqvalues relevant to the polymer melts studied in this work (q<5/3 in three dimensions), we employ the standard expec- tation value ⟨Rk⟩q = Z ∞ 0 RkPq(R)dR,(9) avoiding the use of escort distributions. As shown in the Appendix (Sec. A 4) and supplementary material (Sec. S1.2), for bo...

[4] [4]

Inserting this definition in Eq

Mapping of the q-Gaussian second moment onto an equivalent random walk To relate thisq-Gaussian mean-square end-to-end distance to an effective random walk, we define βq = 3 7−5q 1 Nq l2q ,(11) that equalsβ g whenq=1. Inserting this definition in Eq. (10), the second moment may be expressed as a random walk ofNq steps, each with lengthl q ⟨R2⟩q =N ql2 q ....

[5] [5]

renor- malized

Theq-logarithm and linearization In the Gaussian case (q=1), as shown in Eq. (3), the plot of the logarithm of the volumetric densityρ g(R)∝P g(R)/R2 versusR 2 is linear. For theq-Gaussian, this linearity is pre- served only when utilizing theq-logarithm. Theq-logarithm is the inverse function of theq-exponential [ln q[expq(u)] =u] and is defined fory>0 a...

[6] [6]

Boltzmann–Gibbs entropy (S1) The classical extensive entropy is defined using the volu- metric probability densityρ(R)(units ofL −3) integrated over the 3D volume. Expressed in terms of the radial probability densityP(R)(units ofL −1), it is: S1 =−k B Z ∞ 0 P(R)ln ρ(R) ρ0 dR+const.,(16) whereρ(R) =P(R)/(4πR 2)andρ 0 =1 Å −3 is a reference density introduc...

[7] [7]

experimental

Tsallisq-entropy (S q) To capture possible non-extensive behavior due the ex- istence of different types and relative populations of Kuhn segments, we use the generalized Tsallis entropy. 36,43 In the discrete implementation, the dimensionless probabilityp k is raised to the powerq Sq =k B 1− ∑k pq k q−1 .(18) Asq→1,S q converges toS 1. A value ofq̸=1 sig...

[8] [8]

Every chaini(i=1,

Space and time (Sp&Ti) averaging In this conventional approach (typical of steady-state molecular dynamics sampling), all sampled chain configura- tions in space and time are combined into a single ensemble. Every chaini(i=1, . . . ,N ch), at every framet, contributes an end-to-end distanceR i,t to the overall distribution. The space- time probability dis...