Elastic straining of free-standing monolayer graphene

Since graphene was isolated by mechanical exfoliation in 2004, it has been pursued for electronic and optoelectronic applications due to its near-perfect two-dimensional crystal structure and exceptional mechanical, optical, thermal, and electronic properties. Theory predicts a Young’s modulus of ~1 ± 0.1 TPa and ideal tensile strength of 100–130 GPa, with maximum strains to failure up to ~13–19% (armchair) and ~20–26% (zigzag), with nonlinear elasticity at large strain. In large-area CVD-grown graphene, defects (point defects, grain boundaries, edge defects) are expected to reduce tensile strength, but direct experimental evidence under uniform uniaxial tension has been lacking. Mechanical straining of graphene is also a powerful route to modulate lattice, phonons, reactivity, magnetism, and electronic/optoelectronic properties; however, significant effects generally require large elastic strain, especially for monolayer graphene with zero intrinsic bandgap. Experimentally, tensile strain to failure has often been limited to ~1% due to defects, far below the range needed for robust strain engineering. This motivates a method to impose well-controlled, sample-wide uniaxial elastic strain substantially above 1% on large-area, crack-free monolayer graphene. In this work, the authors develop controlled transfer and shaping of free-standing CVD monolayer graphene and perform in situ tensile tests in SEM using a push-to-pull device, aiming to directly measure elastic properties, stretchability, and fracture behavior of single-crystalline monolayer graphene under uniaxial loading.

Prior experimental studies on graphene’s fracture used pre-cracked samples (monolayer, bilayer, multilayer) in SEM/TEM, observing brittle failure with breaking stresses well below ideal strength, which characterizes fracture toughness but not the realistic strength of crack-free monolayers relevant to devices. Indirect nanoindentation on suspended membranes established an intrinsic Young’s modulus ~1 TPa and strength up to ~130 GPa in pristine regions, and ~15% strength reduction at grain boundaries; however, indentation probes local, non-uniform stress fields and is highly sensitive to local atomic structure and geometry, limiting its relevance to large-area performance. Burst tests on suspended membranes reported wide distributions of burst pressures attributed to wrinkles and defects. Overall, because engineering strength/strain to failure in brittle materials are governed by the weakest locations, a direct, uniaxial tensile test on large-area, crack-free graphene is required for quantitative assessment. Beyond mechanics, extensive literature has shown that elastic strain in graphene can tune electronic, chemical, and optical properties, and strain also plays a role in emergent phenomena in twisted bilayer graphene, underscoring the need to experimentally realize large, uniform, reversible elastic strains.

Sample preparation and transfer: CVD-grown monolayer graphene on copper (Sigma Aldrich), coated with PMMA (<100 nm), was released by Cu etching, rinsed in DI water, and fished onto a push-to-pull (PTP) micromechanical device. After overnight drying to ensure adhesion, PMMA was removed in acetone using a critical point dryer to protect the suspended region over the device gap. The graphene covering the device provided robust clamping by van der Waals adhesion, and the suspended region was patterned into a ribbon by focused ion beam (FIB) using small current to minimize damage (typical sample: width ~3.4 μm, gauge length ~3.0 μm). Raman spectroscopy (514 nm, 1.7 mW) verified monolayer, high-quality graphene (I2D/IG ≈ 3, negligible D peak). TEM at 80 kV imaged edges (single dark line) and SAED confirmed monolayer and single crystallinity across the ribbon via multi-point diffraction.

In situ SEM tensile testing: The PTP device with graphene was mounted on a Hysitron pico-indenter (PI85/P185) inside an FEI Quanta 450 SEM (5–20 kV). Displacement-controlled loading actuated the PTP device while recording load–displacement and SEM videos. Initial cycles pre-stretched draped membranes to define zero-strain. Subsequent cycles applied increasing displacement with full unloading to assess elastic recoverability. The load–displacement response initially reflected device stiffness until the graphene was fully tightened; thereafter, the slope captured combined device+sample stiffness. The tensile stiffness of graphene was obtained by subtracting the device stiffness (measured post-fracture) from the pre-fracture slope. Tensile strain was measured directly from SEM images of gauge length change. To interpret stiffness as Young’s modulus, the stress/strain state was assumed predominantly uniaxial stress (σz=0) over most of the suspended region due to free transverse contraction, validated by FEM; end-clamped regions enforce uniaxial strain locally. Using nominal thickness t=0.335 nm, gauge length l, and width W, 2D and 3D moduli were evaluated from stiffness k (e.g., E2D≈k in uniaxial stress assumption; alternative relations and Poisson’s ratio ν=0.169 considered for completeness). Tensile stress was computed as σ=E3D·ε.

Molecular dynamics (MD) simulations: Atomistic models of free-standing graphene ribbons (length 20 nm, width 10 nm) with periodic boundary along loading and free boundaries along width/out-of-plane were constructed. Edge defects (bare and H-terminated) were created based on ion-irradiation damage motifs. Simulations used LAMMPS with the AIREBO potential; cutoff was set to 0.2 nm to avoid spurious strengthening post-failure. Simulations began at 0 K; temperature rise during loading was ~10 K, making thermal effects negligible. Uniform tensile strain was applied by deforming the simulation box at 1 ns−1; lower rates were checked for consistency. For crack nucleation/propagation in clamped samples, two rows of aromatic rings at ends were fixed or displaced at constant rate. Atomic virial stresses were mapped to 2D meshes for local stress analysis and averaged to obtain tensile stress. Simulations also examined effects of in-plane misalignment and stress localization near clamping ends.

• Monolayer, single-crystalline CVD graphene ribbons exhibited fully recoverable, sample-wide elastic strains of approximately 2.3%, 3.6%, and 4.7% in cyclic loading, with linear load–displacement and complete unloading recovery.
• In monotonic tests to failure, peak engineering strain reached ~5.8% with purely elastic behavior up to brittle fracture.
• From pre- and post-fracture slopes, sample tensile stiffness ~350 N/m was extracted (example: pre-fracture ~460 N/m; device ~110 N/m). Under a uniaxial stress assumption, E2D ≈ 309 N/m, comparable to ab initio predictions (~348 N/m). Using t=0.335 nm, E3D ≈ 920 GPa. Across multiple samples, Young’s modulus ranged ~900–1000 GPa, near theoretical values for pristine graphene.
• The corresponding engineering tensile strength was ~53 GPa for the ~5.8% fracture strain case; representative strengths across samples were ~50–60 GPa.
• Fracture was brittle and initiated near edges close to clamping ends; post-failure load–displacement reverted to device stiffness.
• MD simulations showed that FIB-like edge defects reduce tensile strength; ratios of ideal lattice strength σ0 to measured engineering strengths σm for defective edges match experimental reductions (σ0/σm in the same range). Peak local principal stress near defects (σp) is not a reliable fracture criterion due to anisotropic fracture toughness; average principal stress (σa) correlates better with failure.
• Crack paths observed by post-mortem TEM included armchair and zigzag edges, consistent with anisotropic Griffith-type criteria depending on loading direction relative to lattice.
• In-plane misalignment between clamps had minor effect on measured strength/strain to failure per simulations.
• The demonstrated elastic stretchability and near-ideal stiffness enable realistic strain-engineering of graphene in devices.

The study directly addresses the open question of graphene’s sample-wide elastic stretchability and realistic strength under uniaxial tension. By overcoming gripping and handling through an optimized transfer, clamping, and shaping protocol on a PTP platform, the authors demonstrate large, uniform, reversible elastic strains up to ~5–6% and near-theoretical Young’s modulus (~0.9–1.0 TPa) in single-crystalline CVD monolayer graphene. The measured engineering strengths of ~50–60 GPa, while below the ideal strength (100–130 GPa), are shown via MD to be limited primarily by edge defects introduced during FIB cutting; thus, intrinsic lattice strength in the membrane interior remains high. Fracture initiates near edges at clamping ends due to combined stress concentration and boundary constraints, and proceeds along armchair/zigzag directions consistent with anisotropic fracture mechanics. These results significantly surpass the ~1% strain often reported for finite graphene sheets and validate that high-quality, crack-free, single-crystal CVD graphene can sustain large, uniform elastic strains, thereby enabling practical strain engineering of band structure and properties. The findings also highlight that careful control of edge states and sample preparation could further improve engineering strength towards intrinsic limits, and that device-relevant large-area mechanical performance can be reliably quantified by direct tensile testing rather than local probes.

The authors developed a robust, SEM-based in situ tensile testing strategy for free-standing, single-crystalline CVD monolayer graphene using a push-to-pull microdevice and optimized transfer/clamping. They achieved near-ideal tensile stiffness and Young’s modulus (~900–1000 GPa) with fully recoverable sample-wide elastic strains up to ~6% and representative engineering tensile strengths of ~50–60 GPa. Brittle fracture initiates from edge defects near clamping ends; MD simulations corroborate that FIB-induced edge damage reduces strength while the membrane interior retains near-pristine properties. These results reaffirm graphene’s high mechanical resilience beyond local probing and pave the way for dynamically strain-tuned electronics/optoelectronics and ultra-strong composites. Future work should focus on minimizing edge damage (e.g., alternative cutting methods), optimizing clamping geometries to reduce stress concentrations, exploring higher-aspect-ratio specimens, and extending the approach to other 2D materials and heterostructures for elastic strain engineering.

• Edge defects introduced by FIB cutting reduce engineering strength relative to the ideal lattice; strength is thus limited by preparation-induced edge damage. Alternative patterning may yield higher strengths.
• Interpretation of stiffness as Young’s modulus requires assumptions about the stress/strain state; while uniaxial stress is argued and FEM-validated for most of the suspended region, end-clamped areas experience different constraints, and the stress/strain field is not perfectly uniform.
• A nominal thickness of 0.335 nm is used to convert 2D to 3D modulus and stress, which is a convention and may affect absolute 3D values.
• Possible small in-plane misalignment between clamping ends can introduce non-uniform displacement fields; simulations suggest minor impact, but not eliminated experimentally.
• Tests are performed in SEM under high vacuum and at room temperature; environmental effects (humidity, temperature) and dynamic loading were not explored.
• MD simulations were conducted near 0 K with AIREBO potential and finite sizes (20×10 nm), which may not capture all thermally activated or long-length-scale phenomena.
• Results pertain to high-quality, single-crystalline CVD monolayer graphene; generalization to polycrystalline or defect-rich films may require additional study.