Elastic, plastic, and fracture mechanisms in graphene materials

Anisotropic Fracture of Graphene Revealed by Surface Steps on Graphite

The anisotropic fracture toughness G(θ) is an intrinsic feature of graphene and is fundamental for fabrication, functioning, and robustness of graphene-based devices. However, existing results show significant discrepancies on the anisotropic factor, i.e., the ratio between zigzag (ZZ) and armchair (AC) directions, G_{ZZ}/G_{AC}, both qualitatively and quantitatively. Here, we investigate the anisotropic fracture of graphene by atomic steps on cleaved graphite surfaces. Depending on the relation between the peeling direction and local lattice orientation, two categories of steps with different structures and behaviors are observed. In one category are straight steps well aligned with local ZZ directions, while in the other are steps consisting of nanoscale ZZ and AC segments. Combined with an analysis on fracture mechanics, the microscale morphology of steps and statistics of their directions provides a measurement on the anisotropic factor of G_{ZZ}/G_{AC}=0.971, suggesting that the ZZ direction has a slightly lower fracture toughness. The results provide an experimental benchmark for the widely scattered existing results, and offer constraints on future models of graphene fracture.

Uniaxially Strained Graphene: Structural Characteristics and G-Mode Splitting

The potential use of graphene in various strain engineering applications requires an accurate characterization of its properties when the material is under different mechanical loads. In this work, we present the strain dependence of the geometrical characteristics at the atomic level and the Raman active G-band evolution in a uniaxially strained graphene monolayer, using density functional theory methods as well as molecular dynamics atomistic simulations for strains that extend up to the structural failure. The bond length and bond angle variations with strain, applied either along the zigzag or along the armchair direction, are discussed and analytical relations describing this dependence are provided. The G-mode splitting with strain, as obtained by first principles' methods, is also presented. While for small strains, up to around 1%, the G-band splitting is symmetrical in the two perpendicular directions of tension considered here, this is no longer the case for larger values of strains where the splitting appears to be larger for strains along the zigzag direction. Further, a crossing is observed between the lower frequency split G-mode component and the out-of-plane optical mode at the Γ point for large uniaxial strains (>20%) along the zigzag direction.

Atomistic Mechanics of Torn Back Folded Edges of Triangular Voids in Monolayer WS2

Triangular nanovoids in 2D materials transition metal dichalcogenides have vertex points that cause stress concentration and lead to sharp crack propagation and failure. Here, the atomistic mechanics of back folding around triangular nanovoids in monolayer WS2 sheets is examined. Combining atomic-resolution images from annular dark-field scanning transmission electron microscopy with reactive molecular modelling, it is revealed that the folding edge formation has statistical preferences under geometric conditions based on the orientation mismatch. It is further investigated how loading directions and strong interlayer friction, interplay with WS2 lattice's crack preference, govern the deformation and fracture pattern around folding edges. These results provide fundamental insights into the combination of fracture and folding in flexible monolayer crystals and the resultant Moiré lattices.

Width Dependent Elastic Properties of Graphene Nanoribbons

The mechanical response of graphene nanoribbons under uniaxial tension, as well as its dependence on the nanoribbon width, is presented by means of numerical simulations. Both armchair and zigzag edged graphene nanoribbons are considered. We discuss results obtained through two different theoretical approaches, viz. density functional methods and molecular dynamics atomistic simulations using empirical force fields especially designed to describe interactions within graphene sheets. Apart from the stress-strain curves, we calculate several elastic parameters, such as the Young’s modulus, the third-order elastic modulus, the intrinsic strength, the fracture strain, and the Poisson’s ratio versus strain, presenting their variation with the width of the nanoribbon.

Graphene Oxides Derivatives Prepared by an Electrochemical Approach: Correlation between Structure and Properties

Graphene oxide (GO) can be defined as a single monolayer of graphite with oxygen-containing functionalities such as epoxides, alcohols, and carboxylic acids. It is an interesting alternative to graphene for many applications due to its exceptional properties and feasibility of functionalization. In this study, electrochemically exfoliated graphene oxides (EGOs) with different amounts of surface groups, hence level of oxidation, were prepared by an electrochemical two-stage approach using graphite as raw material. A complete characterization of the EGOs was carried out in order to correlate their surface topography, interlayer spacing, defect content, and specific surface area (SSA) with their electrical, thermal, and mechanical properties. It has been found that the SSA has a direct relationship with the d-spacing. The EGOs electrical resistance decreases with increasing SSA while rises with increasing the D/G band intensity ratio in the Raman spectra, hence the defect content. Their thermal stability under both nitrogen and dry air atmospheres depends on both their oxidation level and defect content. Their macroscopic mechanical properties, namely the Young’s modulus and tensile strength, are influenced by the defect content, while no correlation was found with their SSA or interlayer spacing. Young moduli values as high as 54 GPa have been measured, which corroborates that the developed method preserves the integrity of the graphene flakes. Understanding the structure-property relationships in these materials is useful for the design of modified GOs with controllable morphologies and properties for a wide range of applications in electrical/electronic devices.

Chaotic dynamics of graphene and graphene nanoribbons

We study the chaotic dynamics of graphene structures, considering both a periodic, defect free, graphene sheet and graphene nanoribbons (GNRs) of various widths. By numerically calculating the maximum Lyapunov exponent, we quantify the chaoticity for a spectrum of energies in both systems. We find that for all cases, the chaotic strength increases with the energy density and that the onset of chaos in graphene is slow, becoming evident after more than 10⁴ natural oscillations of the system. For the GNRs, we also investigate the impact of the width and chirality (armchair or zigzag edges) on their chaotic behavior. Our results suggest that due to the free edges, the chaoticity of GNRs is stronger than the periodic graphene sheet and decreases by increasing width, tending asymptotically to the bulk value. In addition, the chaotic strength of armchair GNRs is higher than a zigzag ribbon of the same width. Furthermore, we show that the composition of ¹²C and ¹³C carbon isotopes in graphene has a minor impact on its chaotic strength.

Electronic and optical properties of strained graphene and other strained 2D materials: a review

This review presents the state of the art in strain and ripple-induced effects on the electronic and optical properties of graphene. It starts by providing the crystallographic description of mechanical deformations, as well as the diffraction pattern for different kinds of representative deformation fields. Then, the focus turns to the unique elastic properties of graphene, and to how strain is produced. Thereafter, various theoretical approaches used to study the electronic properties of strained graphene are examined, discussing the advantages of each. These approaches provide a platform to describe exotic properties, such as a fractal spectrum related with quasicrystals, a mixed Dirac-Schrödinger behavior, emergent gravity, topological insulator states, in molecular graphene and other 2D discrete lattices. The physical consequences of strain on the optical properties are reviewed next, with a focus on the Raman spectrum. At the same time, recent advances to tune the optical conductivity of graphene by strain engineering are given, which open new paths in device applications. Finally, a brief review of strain effects in multilayered graphene and other promising 2D materials like silicene and materials based on other group-IV elements, phosphorene, dichalcogenide- and monochalcogenide-monolayers is presented, with a brief discussion of interplays among strain, thermal effects, and illumination in the latter material family.

Theory of electrical conductivity and dielectric permittivity of highly aligned graphene-based nanocomposites

Highly aligned graphene-based nanocomposites are of great interest due to their excellent electrical properties along the aligned direction. Graphene fillers in these composites are not necessarily perfectly aligned, but their orientations are highly confined to a certain angle, [Formula: see text] with 90° giving rise to the randomly oriented state and 0° to the perfectly aligned one. Recent experiments have shown that electrical conductivity and dielectric permittivity of highly aligned graphene-polymer nanocomposites are strongly dependent on this distribution angle, but at present no theory seems to exist to address this issue. In this work we present a new effective-medium theory that is derived from the underlying physical process including the effects of graphene orientation, filler loading, aspect ratio, percolation threshold, interfacial tunneling, and Maxwell-Wagner-Sillars polarization, to determine these two properties. The theory is formulated in the context of preferred orientational average. We highlight this new theory with an application to rGO/epoxy nanocomposites, and demonstrate that the calculated in-plane and out-of-plane conductivity and permittivity are in agreement with the experimental data as the range of graphene orientations changes from the randomly oriented to the highly aligned state. We also show that the percolation thresholds of highly aligned graphene nanocomposites are in general different along the planar and the normal directions, but they converge into a single one when the statistical distribution of graphene fillers is spherically symmetric.

Low-energy theory for strained graphene: an approach up to second-order in the strain tensor

An analytical study of low-energy electronic excited states in uniformly strained graphene is carried out up to second-order in the strain tensor. We report a new effective Dirac Hamiltonian with an anisotropic Fermi velocity tensor, which reveals the graphene trigonal symmetry being absent in first-order low-energy theories. In particular, we demonstrate the dependence of the Dirac-cone elliptical deformation on the stretching direction with respect to graphene lattice orientation. We further analytically calculate the optical conductivity tensor of strained graphene and its transmittance for a linearly polarized light with normal incidence. Finally, the obtained analytical expression of the Dirac point shift allows a better determination and understanding of pseudomagnetic fields induced by nonuniform strains.

Numerical Investigation of the Fracture Mechanism of Defective Graphene Sheets

Despite the unique occurrences of structural defects in graphene synthesis, the fracture mechanism of a defective graphene sheet has not been fully understood due to the complexities of the defects. In this study, the fracture mechanism of the monolayer graphene with four common types of defects (single vacancy defect, divacancy defect, Stone–Wales defect and line vacancy defect) were investigated systematically for mechanical loading along armchair and zigzag directions, by using the finite element method. The results demonstrated that all four types of defects could cause significant fracture strength loss in graphene sheet compared with the pristine one. In addition, the results revealed that the stress concentration occurred at the carbon–carbon bonds along the same direction as the displacement loading due to the deficiency or twist of carbon–carbon bonds, resulting in the breaking of the initial crack point in the graphene sheet. The fracture of the graphene sheet was developed following the direction of the breaking of carbon–carbon bonds, which was opposite to that of the displacement loading.