pith. sign in

arxiv: 2603.22207 · v2 · submitted 2026-03-23 · ❄️ cond-mat.mtrl-sci · cond-mat.dis-nn· cond-mat.mes-hall· cond-mat.soft· physics.app-ph

Universal inverse-cube thickness scaling of projectile penetration energy in ultrathin films

Alessio Zaccone , Tim W. Sirk This is my paper

Pith reviewed 2026-05-15 00:33 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci cond-mat.dis-nncond-mat.mes-hallcond-mat.softphysics.app-ph
keywords ultrathin filmsprojectile penetrationthickness scalinginverse-cube lawshear modulusnonaffine deformationsimpact resistancefinite-size effects
0
0 comments X

The pith

Specific penetration energy in ultrathin films follows a universal inverse-cube law with thickness.

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

The paper shows that ultrathin films of different materials all exhibit increasing resistance to projectile penetration as they get thinner, following the same mathematical rule. The specific energy needed to penetrate the film equals a bulk value plus a term that grows as one over thickness cubed. This behavior holds whether the film is graphene, graphene oxide, or polymer, and whether ordered or disordered. The explanation lies in how very thin layers suppress certain deformation patterns that normally make solids easier to shear. If this holds, it means the enhanced toughness at small scales comes from a basic elastic property rather than material-specific chemistry.

Core claim

The thickness dependence of the specific penetration energy obeys the universal law E_p^*(h) = E_{p,∞}^* + B h^{-3}, independent of chemical composition and degree of disorder. This inverse-cube scaling originates from a finite-size correction to the effective shear modulus due to the suppression of long-wavelength nonaffine deformation modes in confined solids. The scaling describes impact data for multilayer graphene, graphene oxide, and polymer thin films, pointing to a shared elastic origin for their nanoscale impact resistance.

What carries the argument

Finite-size correction to the effective shear modulus from suppression of long-wavelength nonaffine deformation modes in confined solids, producing the inverse-cube term in penetration energy.

If this is right

  • The scaling law holds independently of chemical composition and structural disorder.
  • Thinner films exhibit higher specific penetration energy due to the added inverse-cube term.
  • The bulk limit E_{p,∞}^* is recovered as thickness increases.
  • The same relation accounts for data across graphene, graphene oxide, and polymers.
  • The enhancement traces to an elastic mechanism rather than chemistry-specific effects.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Thickness tuning alone could serve as a general design lever for impact-resistant coatings without altering base chemistry.
  • The nonaffine-mode suppression mechanism may extend to other confined geometries such as thin wires or multilayer stacks.
  • Direct measurements of shear modulus in ultrathin samples should show a corresponding thickness-dependent stiffening.
  • The scaling provides a parameter-free route to predict high-velocity impact behavior from low-strain elastic properties.

Load-bearing premise

The inverse-cube scaling originates from a finite-size correction to the effective shear modulus caused by suppression of long-wavelength nonaffine deformation modes.

What would settle it

Penetration experiments on ultrathin films of varying thickness where measured energy values fail to fit the form constant plus B over thickness cubed.

Figures

Figures reproduced from arXiv: 2603.22207 by Alessio Zaccone, Tim W. Sirk.

Figure 1
Figure 1. Figure 1: FIG. 1. Universal inverse–cube thickness scaling of the specific penetration energy. (a) Multilayer graphene: specific penetration [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
read the original abstract

Ultrathin films of widely different materials exhibit a dramatic enhancement of projectile penetration resistance under high--velocity impact. Despite extensive simulations and experiments, a unifying physical explanation has remained elusive. Here we show that the thickness dependence of the specific penetration energy obeys a universal law, $E_p^*(h)=E_{p,\infty}^*+B h^{-3}$, independent of chemical composition and degree of disorder. The inverse--cube scaling is traced back to a finite--size correction to the effective shear modulus arising from the suppression of long--wavelength nonaffine deformation modes in confined solids. The scaling quantitatively describes impact data for multilayer graphene, graphene oxide, and polymer thin films, revealing a common elastic origin for nanoscale impact resistance.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript claims that the specific penetration energy of ultrathin films obeys a universal scaling E_p^*(h) = E_{p,∞}^* + B h^{-3} independent of chemical composition and disorder. The inverse-cube term is attributed to a finite-size correction in the effective shear modulus arising from suppression of long-wavelength nonaffine deformation modes in confined solids. The form is shown to describe high-velocity impact data for multilayer graphene, graphene oxide, and polymer films.

Significance. If the scaling and its physical origin hold, the work supplies a material-independent explanation for the dramatic enhancement of impact resistance in ultrathin films, with implications for predictive modeling of nanoscale protective materials. The cross-material collapse onto a single functional form is a notable strength.

major comments (2)
  1. [§3, Eq. (7)] §3 (Theory), Eq. (7): The assertion that suppression of long-wavelength nonaffine modes produces a correction to G_eff(h) that enters E_p^* precisely as h^{-3} lacks the explicit mode-sum or continuum calculation. Without the dispersion relation, boundary conditions, or integration that fixes the exponent at -3 (rather than -2 or -4), the functional form reduces to a two-parameter fit whose universality is not predicted a priori.
  2. [Fig. 4] Fig. 4 and associated text: The reported fits across graphene, GO, and polymers achieve good visual agreement, yet no uncertainties on the extracted B values or on the data points themselves are provided. This omission prevents assessment of whether the claimed universality is statistically robust or sensitive to the choice of fitting range.
minor comments (2)
  1. [Abstract] The abstract states the scaling law but does not define the symbols E_p^* and B at first use; a brief parenthetical definition would improve readability.
  2. [Introduction] Notation for the asymptotic value E_{p,∞}^* is introduced without an explicit statement of the thick-film limit it represents; a short clause clarifying the physical meaning would help.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of the work's significance and for the constructive comments, which have helped clarify the presentation. We address each major comment point by point below.

read point-by-point responses

  1. Referee: [§3, Eq. (7)] §3 (Theory), Eq. (7): The assertion that suppression of long-wavelength nonaffine modes produces a correction to G_eff(h) that enters E_p^* precisely as h^{-3} lacks the explicit mode-sum or continuum calculation. Without the dispersion relation, boundary conditions, or integration that fixes the exponent at -3 (rather than -2 or -4), the functional form reduces to a two-parameter fit whose universality is not predicted a priori.

    Authors: We agree that an explicit derivation strengthens the claim. In the revised manuscript we have expanded §3 to include a continuum calculation of the finite-size correction to the effective shear modulus. Starting from the known nonaffine displacement spectrum in disordered solids, we impose free-surface boundary conditions and introduce an infrared cutoff at wavevector q_min ~ 1/h set by the film thickness. The resulting integral over the mode contribution to the modulus correction evaluates to a term proportional to h^{-3}, confirming that the exponent is fixed by the long-wavelength suppression rather than being an arbitrary fit parameter. The derivation is now presented immediately before Eq. (7) and is independent of microscopic details, thereby supporting the observed material-independent collapse. revision: yes

  2. Referee: [Fig. 4] Fig. 4 and associated text: The reported fits across graphene, GO, and polymers achieve good visual agreement, yet no uncertainties on the extracted B values or on the data points themselves are provided. This omission prevents assessment of whether the claimed universality is statistically robust or sensitive to the choice of fitting range.

    Authors: We appreciate this point. In the revised manuscript we have added vertical error bars to all data points in Fig. 4, representing the standard deviation obtained from repeated impact simulations (or experimental replicates where available). We also report the fitted values of B together with their 1σ uncertainties extracted from nonlinear least-squares minimization with bootstrap resampling. A new supplementary section examines the sensitivity of B to the fitting range; the extracted B values remain consistent within the reported uncertainties for all reasonable choices of lower and upper thickness cutoffs, supporting the robustness of the claimed universality. revision: yes

Circularity Check

1 steps flagged

h^{-3} scaling asserted from nonaffine elasticity but coefficient B fitted to data without explicit mode-sum derivation fixing the exponent

specific steps
  1. fitted input called prediction [Abstract]
    "Here we show that the thickness dependence of the specific penetration energy obeys a universal law, $E_p^*(h)=E_{p,∞}^*+B h^{-3}$, independent of chemical composition and degree of disorder. The inverse--cube scaling is traced back to a finite--size correction to the effective shear modulus arising from the suppression of long--wavelength nonaffine deformation modes in confined solids. The scaling quantitatively describes impact data for multilayer graphene, graphene oxide, and polymer thin films"

    The inverse-cube form is asserted as originating from nonaffine mode suppression in confined solids, yet no explicit derivation (mode sum or continuum elasticity calculation) is shown that fixes the power at exactly -3. Instead the coefficient B is fitted to the impact data sets, rendering the claimed universal law a post-hoc fit within the paper's own framework rather than an independent first-principles prediction.

full rationale

The paper claims a universal law E_p^*(h)=E_{p,∞}^* + B h^{-3} derived from finite-size correction to G_eff via suppression of long-wavelength nonaffine modes. However, the abstract and results present the functional form and trace it to elasticity without displaying the continuum or mode-sum calculation that enforces precisely -3 (as opposed to other powers). B is adjusted to match penetration data for graphene, GO and polymers, so the quantitative description reduces to a fitted parameter. This makes the 'universal prediction' statistically forced by the data fit rather than parameter-free from first principles. No load-bearing self-citation chain is quoted, but the central claim lacks independent verification of the exponent.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The claim rests on two adjustable parameters (E_{p,∞}^* and B) plus the domain assumption that long-wavelength nonaffine modes are suppressed in confined geometry; no new entities are introduced.

free parameters (2)
  • B
    Prefactor of the h^{-3} correction; adjusted to match penetration data across films.
  • E_{p,∞}^*
    Bulk asymptotic penetration energy; taken from thick-film limit or fitted.
axioms (1)
  • domain assumption Finite-size correction to effective shear modulus arises from suppression of long-wavelength nonaffine deformation modes in confined solids
    Invoked in the abstract to derive the inverse-cube term.

pith-pipeline@v0.9.0 · 5436 in / 1232 out tokens · 52169 ms · 2026-05-15T00:33:10.526066+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

18 extracted references · 18 canonical work pages

  1. [1]

    Zaccone and K

    A. Zaccone and K. Trachenko, Explaining the low- frequency shear elasticity of confined liquids, Proceedings of the National Academy of Sciences117, 19653 (2020), https://www.pnas.org/doi/pdf/10.1073/pnas.2010787117

    work page doi:10.1073/pnas.2010787117 2020

  2. [2]

    Milkus and A

    R. Milkus and A. Zaccone, Local inversion symmetry breaking controls the boson peak in glasses and crystals, Phys. Rev. E95, 032602 (2017)

    work page 2017

  3. [3]

    Zaccone,Theory of Disordered Solids(Springer, Cham, 2023)

    A. Zaccone,Theory of Disordered Solids(Springer, Cham, 2023)

    work page 2023

  4. [4]

    Zaccone and E

    A. Zaccone and E. Scossa-Romano, Approximate an- alytic description of nonaffine response in amorphous solids, Phys. Rev. B83, 184205 (2011)

    work page 2011

  5. [5]

    Lemaˆ ıtre and C

    A. Lemaˆ ıtre and C. Maloney, Sum rules for the quasi- static and visco-elastic response of disordered solids at zero temperature, Journal of Statistical Physics123, 415 (2006)

    work page 2006

  6. [6]

    A. E. Phillips, M. Baggioli, T. W. Sirk, K. Trachenko, and A. Zaccone, UniversalL −3 finite-size effects in the viscoelasticity of amorphous systems, Phys. Rev. Mater. 5, 035602 (2021)

    work page 2021

  7. [7]

    J. A. Zukas,Impact Dynamics(Wiley, 1982)

    work page 1982

  8. [8]

    Goldsmith,Impact: The Theory and Physical Be- haviour of Colliding Solids(Dover, 2001)

    W. Goldsmith,Impact: The Theory and Physical Be- haviour of Colliding Solids(Dover, 2001)

    work page 2001

  9. [9]

    A. L. Florence, Penetration of targets by high-speed pro- jectiles, Journal of Applied Physics38, 478 (1967)

    work page 1967

  10. [10]

    N. M. Pugno, A general shape/size-effect law for nanoin- dentation, Acta Materialia55, 1947 (2007)

    work page 1947

  11. [11]

    R. A. Biz˜ ao, L. D. Machado, J. M. de Sousa, N. M. Pugno, and D. S. Galv˜ ao, Scale effects on the ballistic penetration of graphene sheets, Scientific Reports8, 6750 (2018)

    work page 2018

  12. [12]

    J.-H. Lee, P. E. Loya, J. Lou, and E. L. Thomas, Dynamic mechanical behavior of multilayer graphene via super- sonic projectile penetration, Science346, 1092 (2014)

    work page 2014

  13. [13]

    K. Yoon, A. Ostadhossein, and A. C. T. van Duin, Atomistic-scale simulations of the chemomechanical be- havior of graphene under nanoprojectile impact, Carbon 99, 58 (2016)

    work page 2016

  14. [14]

    B. Z. G. Haque, S. C. Chowdhury, and J. W. Gillespie, Molecular simulations of stress wave propagation and perforation of graphene sheets under transverse impact, Carbon102, 126 (2016)

    work page 2016

  15. [15]

    H. L. White, W. Chen, N. M. Pugno, and S. Keten, Rate-dependent molecular size effects govern the inverse thickness dependence of specific penetration energy in nanoscale thin films, Extreme Mechanics Letters81, 102419 (2025)

    work page 2025

  16. [16]

    J. Hyon, O. Lawal, O. Fried, R. Thevamaran, S. Yazdi, M. Zhou, D. Veysset, S. E. Kooi, Y. Jiao, M.-S. Hsiao, J. Streit, R. A. Vaia, and E. L. Thomas, Extreme en- ergy absorption in glassy polymer thin films by super- sonic micro-projectile impact, Materials Today21, 817 (2018)

    work page 2018

  17. [17]

    C. M. Stafford, C. Harrison, K. L. Beers, A. Karim, E. J. Amis, M. R. VanLandingham, H.-C. Kim, W. Volksen, R. D. Miller, and E. E. Simonyi, A buckling-based metrol- ogy for measuring the elastic moduli of polymeric thin 4 films, Nature Materials3, 545 (2004)

    work page 2004

  18. [18]

    C. M. Stafford, S. Guo, C. Harrison, and M. Y. Chiang, Combinatorial and high-throughput measurements of the modulus of thin polymer films, Macromolecules39, 5095 (2006)

    work page 2006