Probing Intrinsic Elastic Properties of Multilayer Graphene -- a New Mechanical Constant

We present measurements on in-plane Young’s modulus and the Grüneisen parameter of multilayer graphene with varying number of layers, obtained through {\it in situ} bulge tests. Accurate determination of their elastic parameters poses a significant experimental challenge due to the substantial differences in mechanical behavior between intra- and inter-layers. To address this, we develop a novel theoretical model with first-principles calculations to investigate thickness-dependent incomplete strain transfer between the layers. Our findings show that the experimentally measured elastic constants, which deviate from computed intrinsic values, fail to fully capture ideal mechanical couplings between layers. As a solution, we propose a new mechanical modulus that integrates the Grüneisen parameter and in-plane Young’s modulus, providing a more reliable representation of their mechanical properties, independent of unavoidable interlayer effects.

💡 Research Summary

The authors present a comprehensive study of the intrinsic elastic properties of multilayer graphene, focusing on the in‑plane Young’s modulus (E) and the Grüneisen parameter (γ). Using an in‑situ bulge test combined with atomic force microscopy (AFM) and Raman spectroscopy, they measured the bulge height (h) and the pressure‑dependent shift of the G‑band for samples ranging from a single layer up to twenty‑one layers. The conventional bulge‑test analysis assumes perfect clamping and a spherical bulge shape, leading to a simple relation P = k h³ where k ∝ E. However, for multilayer van‑der‑Waals (vdW) stacks, incomplete strain transfer due to interlayer sliding invalidates this assumption, especially as thickness increases.

To address this, the paper introduces the Bulge Test Model (BTM). BTM treats the supported region (r > R) as experiencing a shear stress τ that varies linearly through the stack: the bottom layer experiences both graphene‑substrate shear (τ_gs) and graphene‑graphene shear (τ_gg), while upper layers experience only τ_gg. The model derives a modified bulge height h = h₀ + h₁, where h₀ is the ideal spherical contribution and h₁ = –u_R accounts for the inward pulling of the membrane at the rim caused by shear. The dimensionless parameters ζ and ξ, functions of pressure P, shear τ, Young’s modulus E, radius R, and thickness t, control the magnitude of h₁.

Applying BTM with the density‑functional‑theory (DFT) intrinsic values (E_DFT ≈ 1075 GPa, γ_DFT independent of Nₗₐᵧ) reproduces the experimentally observed “apparent” Young’s modulus E_ap, which decreases with layer number and saturates around Nₗₐᵧ ≈ 10. The saturation arises because the effective shear stress transitions from the lower τ_gs (≈ 40 kPa, graphene‑substrate) to the higher τ_gg (≈ 1640 kPa, inter‑graphene) as more layers are added. Analytical approximations show that when P ≈ τ, the effective stiffness coefficient k_eff = P h³ drops sharply, leading to an underestimation of E if the simple model is used.

For stacks thicker than ten layers, the measured E_ap deviates further from BTM predictions. The authors attribute this to pressure‑induced radial wrinkles and local delamination near the hole edge. Wrinkles introduce circumferential compression, reducing the local modulus, while delamination lowers the effective shear τ_gs, both of which amplify the apparent softening. These phenomena mark a mechanical crossover from two‑dimensional (2D) behavior, where interlayer sliding dominates, to three‑dimensional (3D) bulk‑like response. The critical crossover thickness is identified as roughly ten graphene layers.

Recognizing that neither E nor γ alone can capture the intrinsic mechanical coupling in multilayer vdW systems, the authors propose a new mechanical constant M = E · γ. Because M combines the stiffness and anharmonic lattice response, it remains robust against interlayer sliding and can be directly extracted from the bulge‑test data. This constant provides a more reliable metric for comparing the mechanical performance of multilayer graphene and other vdW materials, facilitating the design of flexible electronics, nano‑electromechanical systems, and strain‑engineered devices.

In summary, the paper delivers (i) high‑precision experimental measurements of E and γ for multilayer graphene, (ii) a physically grounded BTM that explains the discrepancy between experiment and DFT by incorporating interlayer shear, (iii) identification of a thickness‑driven 2D‑to‑3D mechanical crossover around ten layers, and (iv) the introduction of a novel, shear‑independent mechanical constant M that unifies elastic and anharmonic properties for vdW stacks.