Modeling Compressive Instability in Two-Dimensional Ti2COx MXenes

Technical Summary: Modeling Compressive Instability in Two-Dimensional Ti2COx MXenes

Problem Statement
While the mechanical properties of MXenes under tensile loading have been extensively characterized, their behavior under compressive and shear loads remains largely unexplored. In practical applications such as flexible electronics, wearable devices, and MXene-polymer composites, nanosheets frequently encounter compressive stresses due to bending, matrix shrinkage, or external loading. This compressive environment can lead to buckling, wrinkling, or structural failure. Current literature lacks a fundamental understanding of the compressive instability mechanisms of MXenes, particularly regarding the effects of surface terminations, atomic defects, and size-dependent phenomena that classical continuum mechanics may fail to capture.

Methodology
The study employs reactive molecular dynamics (MD) simulations using the Large Scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) and the ReaxFF force field (parameterized in [36]) to investigate the compressive and post-buckling behavior of Ti2C and Ti2CO2 nanosheets.

• Simulation Setup: Nanosheets with in-plane dimensions of 14×14 nm² were modeled at 1 K to isolate intrinsic mechanical responses. Simulations utilized a 0.2 fs timestep.
• Loading Conditions: The study examined uniaxial and biaxial compression, as well as shear loading, along both armchair and zigzag crystallographic directions.
• Variables Investigated: The research systematically analyzed the effects of:
  • Strain rates (ranging from $0.725 \times 10^6$ to $2.85 \times 10^8$ s⁻¹).
  • Defect types: Isolated point defects (Ti and C vacancies), random vacancy clusters, and specific Ti vacancy clusters.
  • Lateral confinement pressure (modeled via harmonic springs/Winkler foundation).
  • Surface termination (comparing pristine Ti2C with oxygen-terminated Ti2CO2).
• Validation and Modeling: The ReaxFF potential was validated against Density Functional Theory (DFT) and experimental data for structural parameters and Young's modulus. To address discrepancies between atomistic results and classical theory, a stress-driven nonlocal continuum model was implemented and calibrated against the MD data.

Key Results

1. Strain Rate and Anisotropy: Higher strain rates resulted in increased buckling stresses and strains, accompanied by more unstable post-buckling behavior characterized by larger stress oscillations. Ti2C nanosheets exhibited greater resistance to buckling along the armchair direction compared to the zigzag direction.
2. Defect Influence:
   • Isolated Point Defects: Up to a 3% defect fraction, isolated vacancies primarily reduced the critical buckling stress and increased post-buckling stress oscillations but did not significantly alter the buckling strain or global mode shapes. Ti vacancies had a more pronounced weakening effect than C vacancies.
   • Vacancy Clusters: Clusters of Ti vacancies (specifically groups of 4–5 vacancies) fundamentally altered the deformation mechanism. Instead of classical thin-plate buckling, the nanosheets exhibited localized, wavy out-of-plane deflections driven by severe localized compliance around the defect clusters.
3. Surface Termination and Confinement:
   • Oxygen Termination: The presence of –O surface groups substantially increased the buckling stress (from ~1 GPa to ~3.5 GPa) and reduced the anisotropy between loading directions. However, Ti2CO2 nanosheets were found to fracture at large compressive strains, whereas pristine Ti2C retained structural integrity up to strains exceeding 0.35.
   • Lateral Confinement: Applying lateral confinement pressure (simulating an elastic matrix) significantly increased both the buckling stress and strain while dampening post-buckling oscillations.
4. Buckling Modes: Under uniaxial compression, the nanosheets exhibited plate-like first buckling modes. Under biaxial compression, a symmetric, dome-like buckling shape was observed. Shear loading produced elliptical deflection modes.
5. Continuum vs. Atomistic: Classical Euler-Bernoulli beam theory significantly underestimated the critical buckling strains of the nanosheets. A stress-driven nonlocal continuum formulation, calibrated with a characteristic length parameter of 5 Å, successfully reproduced the size-dependent buckling strains observed in the MD simulations.

Significance and Claims
The paper claims to fill a critical gap in the literature by providing the first systematic atomistic investigation into the compressive and shear instability of MXene nanosheets. The study demonstrates that:

• Classical continuum mechanics is insufficient for predicting the buckling behavior of 2D MXenes due to size effects, necessitating nonlocal formulations for accurate modeling.
• The mechanical stability of MXenes is highly sensitive to the specific nature of defects (isolated vs. clustered) and surface chemistry.
• Lateral confinement and surface termination are effective levers for enhancing compressive stability.

The authors conclude that these findings provide a fundamental understanding of MXene morphological transformations under compression. This knowledge is presented as essential for designing robust MXene-based composites and flexible electronics, and it is suggested that these insights may stimulate future research into MXene-based nanotubes, nanoscrolls, and folded architectures. The study emphasizes that the observed trends are mechanistic insights expected to be qualitatively transferable to lower strain-rate conditions, despite the high strain rates inherent to MD simulations.