Jiao Shi*ab,
Jia-Long Zhanga,
Jia-Xing Jia and
Bo Song*c
aCollege of Water Resources and Architectural Engineering, Northwest A&F University, Yangling 712100, China. E-mail: shijiaoau@163.com
bState Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian 116024, China
cInstitute of Advanced Structure Technology, Beijing Institute of Technology, Beijing, 100081, China. E-mail: songbo943@163.com
First published on 4th January 2021
Developing novel graphene-based materials with unique mechanical properties is of significance to meet the requirements in new applications. The pristine graphene shows a brittle fracture when the stretching strain on it exceeds the critical value. Further, it fails to bear the external load. Herein, to enhance the ductility of the pristine graphene, we proposed a corrugated sandwich carbon network based on few-layered graphene, in which the two surface layers are bonded with several corrugated core layers via benzene molecules. The effects of factors such as the geometry, temperature, and strain rate on the ductility of the carbon network were evaluated using the uniaxial tension tests by molecular dynamics simulations. Results show that the new carbon material has more than one peak fracture strain in stretching. The second peak fracture strain is proportional to the length difference between the surface layers and core layers. Hence, the carbon network has a tunable ductility, which suggests a flexible design of such novel materials in a nanostructure/nanodevice with large deformation.
At the nanoscale, the motion and/or deformation of an object is usually realized by a probe.15,21,22 The large-deformed object should have high ductility under an external load. The literature indicates that when a pristine graphene ribbon is under uniaxial tension along the armchair and zigzag directions, its critical fracture strains are about 14% and 21%, respectively, at room temperature.19,23 Once the stretching strain exceeds the corresponding critical value, it breaks in brittle style. As a result, the low ductility hampers the pristine graphene's application in a nanostructure/nanodevice with large deformation.
For enhancing the graphene ductility, this study proposes a new sandwich-like graphene/benzene-based carbon network (Fig. 1), which mimics the sandwich structures.24–26 In the new carbon network, the two surface layers and core layers are bonded with benzene molecules. The tensile properties of the new carbon network are evaluated by molecular dynamics simulations.
Model | Nc | Np | Lp (Å) | hc (Å) | h (Å) | Surface layer | Core layer | RL layer |
---|---|---|---|---|---|---|---|---|
P3C2 | 2 | 3 | 44.28 | 8.0 | 33.04 | 1080C | 1200C | 111.1% |
P3C3 | 3 | 3 | 44.28 | 8.0 | 46.72 | 1080C | 1200C | 111.1% |
P4C2 | 2 | 4 | 31.98 | 8.0 | 33.04 | 1040C | 1200C | 115.4% |
P4C3 | 3 | 4 | 31.98 | 8.0 | 46.72 | 1040C | 1200C | 115.4% |
P5C2 | 2 | 5 | 24.60 | 8.0 | 33.04 | 1000C | 1200C | 120.0% |
P5C3 | 3 | 5 | 24.60 | 8.0 | 46.72 | 1000C | 1200C | 120.0% |
P6C2 | 2 | 6 | 19.68 | 8.0 | 33.04 | 960C | 1200C | 125.0% |
P6C3 | 3 | 6 | 19.68 | 8.0 | 46.72 | 960C | 1200C | 125.0% |
(1) |
For a deep understanding of the carbon network's tunable ductility, the uniaxial tension tests are conducted with constant strain rate loading. Before stretching, the initial configuration shown in Fig. 1 is reshaped by minimizing its potential energy. Then, the system is relaxed in NPT (constant number of atoms, constant pressure and temperature of the system) ensemble for 10 ps. The environment temperature is controlled by the Nose–Hoover thermostat.38,39 After that, a uniaxial tensile loading with a specified strain rate is applied along the x-direction without any pressure in the other two directions. The time step for the integration of atoms' motion equations is set as 1 fs.
Moreover, the uniaxial tensile tests for a mono-layer pristine graphene ribbon along the same direction are conducted for comparison. The size of the pristine graphene ribbon is 131.6 Å × 19.9 Å, which is similar to that of the surface layers in the P3C3 model. For mitigating the size effect, periodic boundaries are adopted in the stretching plane. Critical strains of PGN under different temperatures obtained here (see the second column in Table 2) agree well with the corresponding values given by Shi et al.19 and Zhao et al.,23 which also confirm the validity of the simulation methods used here.
T/K | References | PGN | ●Peak 1 | ↓Δσ/GPa | ♦Trough | ▼Peak 2 | ||||
---|---|---|---|---|---|---|---|---|---|---|
εc | εc | 1st σ/GPa | 1st εc | σ/GPa | εTr | 2nd σ/GPa | Δε | 2nd εc | ||
8 | 26.9% (ref. 19) | 27.0% | 33.98 | 19.1% | 22.79 | 11.19 | 20.4% | 22.70 | 8.2% | 28.6% |
50 | 25.0% | 32.69 | 19.1% | 20.85 | 11.84 | 18.7% | 20.61 | 7.8% | 26.5% | |
100 | 23.7% (ref. 19) | 23.8% | 31.62 | 18.5% | 21.10 | 10.52 | 18.8% | 20.03 | 7.0% | 25.8% |
200 | 23.1% | 29.46 | 17.4% | 19.62 | 9.84 | 17.8% | 19.38 | 7.4% | 25.2% | |
300 | 20.4% (ref. 19)/20.0% (ref. 23) | 21.7% | 26.22 | 16.4% | 18.31 | 7.91 | 16.9% | 18.12 | 7.7% | 24.6% |
500 | 18.5% | 19.66 | 13.6% | 14.49 | 5.17 | 15.6% | 15.95 | 7.6% | 23.2% |
VPE = (P(t) − P(t0))/N, | (2) |
The stretching strength of a sample will be expressed by the critical values of x-strain, and the effects of the major factors, such as temperature, strain rate/loading speed, Np and Nc, on the strength will be demonstrated. Accordingly, the stress along the stretching direction, referred to as x-stress, will be recorded at each x-strain (the engineering strain rate times steps in the x-direction). Actually, the x-stress is the virial stress40 considering 3.4 Å as the thickness of a surface layer. The normal stress in the x-direction in the virial stress tensor is defined as
(3) |
When stretching under the given engineering strain rate έ, the length of the sample can be calculated by
L(t) = L(t0) × [1 + έ(t − t0)], | (4) |
In Fig. 2, the VPE of the four layers (i.e., top, core 1, core 2 and bottom) in the P3C2 sample at a temperature between 8 K and 500 K is illustrated. The results indicate several characteristics of the sample under uniaxial tension.
Fig. 2 Histories of VPE per atom of the P3C2 sample under the strain rate of 0.1% per ps at different temperatures. Note that 1st strain = 1st εc and 2nd strain = 2nd εc. |
First, the two surface layers have approximately the same critical strain, i.e., 1st critical strain labeled as 1st εc. Similarly, the two core layers have the same 2nd critical strain (and labeled as 2nd εc). The two critical strains demonstrate the sample's stretching strength. Details on ductility will be given below.
Second, temperature influences the tensile strength significantly. For example, at a temperature below 100 K, the two surface layers have a critical strain of ∼19%. However, the value of 1st εc reduces to 16.4% at 300 K or 13.6% at 500 K, which is lower than that of the pristine graphene under the same conditions (see Tables 2 or S1 in ESI†). This reduction is due to the sp3 bonds that connect surface layers and core layers via benzene molecules.
Third, the VPE curve of a surface layer contains two stages, i.e., continuous increasing and a sharp drop after the strain exceeds the 1st εc. However, the VPE curve of the core layer has one more stage, i.e., VPE keeps unchanged near zero with a small x-strain, which demonstrates that the core layer has no elastic deformation in this stage.
Fourth, for the sample with x-strain between 11.5% and 20%, the VPE of the core layers decreases first and then increases untill approaching 1st εc. The VPE decreases due to two facts. One is the core layers being attracted by and attaching to the surface layer (Fig. 3a). The other is the core layer attached during deformation (snapshot in Fig. 3a with x-strain of 12.3%). In this period, the core buckles41–44 and the carbon network shrinks along the y-direction according to the snapshot in Fig. 3a. After that, the core layers work together with the two surface layers against the stretching load, which means that the system is in the co-working state. For example, at 8 K, the carbon network is kept undamaged before x-strain approaching 19%. However, at 19.1% (1st εc), the surface layer breaks (Fig. 3b). The core layers break after the x-strain reaches 28.6% (2nd εc) (Fig. 3c and Table 2), which is greater than that of the pristine graphene under the same conditions (detailed data can also be found in Table S1 in ESI†). It means that the sandwich carbon network's ductility is enhanced when considering the bond breakage in surface layers as “plastic deformation”.
Fig. 4 shows the stress–strain curves of the P3C2 sample under different stretching conditions. In Fig. 4a, the curves differ slightly before x-strain approaching 11.5%, and the (first) slope/modulus is ∼130 GPa (Fig. 4b). When x-strain is between 12.5% and 19.1%, the increment of stress is ∼15.3 GPa, which means the (second) modulus is ∼232 GPa. The second modulus of the sample is much higher than the first modulus due to two facts, which are the four layers work together in the co-working period, and the sample shrinks in y-direction leading to a lower volume of the sample (see eqn (3)). For the sample with the x-strain between 20% and 28.5%, the stress increases 10.5 GPa, which indicates that the (third) modulus of the sample is ∼124 GPa. It is lower than the first modulus only because the simulation box's volume becomes higher when the surface layer breaks (Fig. 3c).
Fig. 4 Stress–strain curves of the P3C2 sample under tension, (a) at different temperatures with strain rate of 0.1% per ps, and (b) at different loading speeds with temperature of 300 K. |
From eqn (3), we know that the magnitude of stress is inversely proportional to the simulation box's volume. To reflect the stretching strength of the carbon network, we measure the strength using the two critical strains. As shown in Table 2, 1st εc decreases with the increase in temperature, and 2nd εc changes from 28.6% at 8 K to 23.2% at 500 K. However, after breaking of surface layers, the increment of strain is not sensitive to temperature.
Fig. 4b shows the stress–strain curves to different loading speeds. The curves are overlapped when x-strain is between 0 and 11.5%, or between 13.2% and 14.6%, or between 18% and 24.6%, which means that the loading speeds slightly influence the three moduli. When x-strain is in the interval (11.5%, 13.2%), noticeable differences exist among the curves since the core layers have different buckling process at different loading speeds, when in the interval of (16.4%, 17.4%), curves differ, as the surface layers have different breakage processes.
(b) P3C3 model
Fig. 5a gives the VPE curves of the five layers in the P3C3 model under the strain rate of 0.1% per ps at 300 K. The two surface layers have the same values of VPE before breakage happens. The top layer breaks at x-strain = 16.7%, while the bottom layer breaks when the x-strain exceeds 17.4%. The three core layers have approximately the same VPE in tension before breakage happening at x-strain = 25.8% (Fig. 5b). The co-working region of x-strain is between 12.3% and 16.6% (highlighted with a gray background in Fig. 5a).
Fig. 5 VPE per atom of the P3C3 sample under the strain rate of 0.1% per ps at 300 K. (a) VPE curves, (b) snapshots. |
The same phenomenon, i.e., VPE of the core layers decreases soon after the x-strain > 11.9%, as shown in Fig. 5a. According to the snapshots in Fig. 5b, the sample shrinks sharply between x-strain = 11.9% and 12.3%. The core structure buckles,41–44 and the three core layers attach. Between 16.6% and 17.4% of x-strain, VPE of the core layers has fluctuations caused by the co-working of the surface layers and the core 1 and core 3 layers in deformation (similar to 25.7%-snapshot in Fig. 5b).
Fig. 6 shows the stress–strain curves of the P3C3 sample considering the effects of temperature and loading rate. It demonstrates that the critical strains are sensitive to temperature but insensitive to the loading rate. The results listed in Table 3 indicate that the critical strains are slightly different from those obtained from the P3C2 model at the same temperature below 500 K. Hence, the critical strain is independent of the number of core layers.
Fig. 6 Stress–strain curves of the P3C3 sample under tension. (a) At different temperatures at the strain rate of 0.1% per ps, (b) at different loading speeds with a temperature of 300 K. |
T | PGN | ●Peak 1 | ↓Δσ/GPa | ♦Trough | ▼Peak 2 | ||||
---|---|---|---|---|---|---|---|---|---|
εc | 1st σ/GPa | 1st εc | σ/GPa | εTr | 2nd σ/GPa | Δε | 2nd εc | ||
8 K | 27.0% | 28.2 | 19.1% | 15.5 | 12.6 | 19.7% | 24.7 | 9.2% | 28.9% |
50 K | 25.0% | 28.0 | 19.0% | 15.3 | 12.7 | 19.4% | 23.1 | 7.8% | 26.9% |
100 K | 23.8% | 27.0 | 19.2% | 14.7 | 12.2 | 19.7% | 22.1 | 7.0% | 26.7% |
200 K | 23.1% | 25.3 | 18.1% | 14.5 | 11.1 | 18.6% | 20.6 | 7.4% | 25.4% |
300 K | 21.7% | 22.4 | 16.6% | 13.5 | 9.0 | 17.5% | 21.0 | 7.7% | 25.7% |
500 K | 18.5% | 14.2 | 12.9% | 11.8 | 2.4 | 13.3% | 17.2 | 7.6% | 23.3% |
However, compared to P3C2 at the same temperature, P3C3 has a lower stress level due to the larger volume with the same critical strain according to eqn (3). Comparing the values of a sharp jump of stress (↓σ) at the first critical strains in Tables 2 and 3, P3C3 has a lower stress decline than P3C2, i.e., the stress fluctuation of P3C3 is weaker. It implies that more core layers will reduce stress fluctuation, which is significant for the application of the carbon network in bearing displacement load.
We also calculated the three moduli of P3C3 according to its stress–strain curves. It shows that the first modulus is only 103 GPa (Fig. 6b), which is less than 130 GPa of P3C2. In this stage, only surface layers bear the external load. The two samples have the same surface layers, but P3C3 has a larger volume leading to lower stress levels. Hence, P3C3 has a lower 1st modulus. The second and third moduli are 217 GPa and 128 GPa, respectively, which are slightly different from those calculated for P3C2. In this stage, only the core layers are under stretching, and the sample volume is approximately proportional to the number of core layers. Therefore, the 2nd or the 3rd modulus of the two samples differ slightly. It is also the major reason for providing no more discussion on the samples with more core layers in this study.
Fig. 7 Stress–strain curves of the samples under loading speed of 0.1% per ps at different temperatures. (a) P4C2 and P4C3, (b) P5C2 and P5C3, and (c) P6C2 and P6C3. |
Another interesting phenomenon in Fig. 7 is that more sudden jumps in the stress–strain curves happen after the surface layers break (i.e., x-strain > 1st εc), when the sample has more periodic cells. The reason is that the surface breakage first happens at one benzene-bonded line in the surface layer (either the top or bottom). However, the broken surface layers and the core layers are still connected via the rest lines of benzene molecules (Fig. 8b and 9b). It also indicates the ductility of the sandwich carbon network is improved when the local breakages are considered as the plastic deformation of the sample. During further stretching, the co-working segments in the surface layers break one by one (Fig. 8b and 9b), and each breakage leads to a sudden drop of stress. In further relaxing the system, VPE (Fig. 8a and 9a) decreases due to self-attraction between core layers.
When comparing the two samples with the same Np but different Nc, we find that the difference between critical strains is negligible at the same temperature (Fig. 7 and Table S1 in the ESI†). It further verifies that the number of core layers slightly influences the critical strains.
When observing the sudden change of the configurations in Fig. 8b, it can be found that the sample has no buckling of core layers before the surface layer breaks (comparing snapshots to x-strain = 14.9% and 15.8%). The reason is that the core layers are longer than the surface layers. The core layer is still not bearing the external deformation when the surface layers break entirely. Without heavy elongation, the core layers cannot get closer to each other and have no chance to attach. After the breakage of surface layers, the core layer shrinks quickly in the transverse direction to bear the stretching deformation.
Fig. 10 indicates that the 2nd critical strain is proportional to Np but decreases with an increase in temperature. More importantly, the 2nd critical strains of the sandwich carbon network are greater than those of pristine graphene under the same conditions, particularly with higher Np in the system. For example, when RLayer = 125.0% (Np = 6), the critical strain at room temperature (300 K) can reach 39.7%, which is almost twice as much as that for pristine graphene (21.7%).
Alternatively, compared to that of pristine graphene, the ductility of the sandwich carbon network is highly improved. Since the 2nd critical strain increases linearly with Np, a new carbon network with higher ductility can be designed with specific Np. This will benefit its application in bearing large deformation.
First, due to the length difference between the surface and core layers, the carbon network under uniaxial tensile load contains two critical strains. The surface layers break at the 1st critical strain, while the core layers break at the higher 2nd critical strain. The 2nd critical strain represents the ductility of the carbon network.
Second, when the difference between the surface and core layers' lengths is less than 20%, the stress–strain curve of the carbon network has three different moduli. The 1st modulus appears when only the surface layers bear the load. The 2nd modulus is the highest among the three moduli due to the co-working of all the layers. The 3rd modulus reflects the stiffness of the core layers.
Third, the critical strains of the sandwich carbon network decrease with temperature, but slightly depend on either loading speeds or the number of core layers.
Finally, the 2nd critical strain in the new carbon network is higher than that of the pristine graphene and increases with an increase in Np. Alternatively, compared to that of the pristine graphene, the ductility of the new carbon network is improved and tunable.
Footnote |
† Electronic supplementary information (ESI) available: Detailed information of critical strains and stresses for all the models adopted in this study. See DOI: 10.1039/d0ra09094h |
This journal is © The Royal Society of Chemistry 2021 |