Abstract:Recent achievement of BN-graphene alloy material has enabled the potential of bandgap tuning through both sub-10 nm width control and BN concentration variation. However, its mechanics, which is necessary for prediction of stability in functional applications, is not well studied. Here, molecular dynamics simulation is performed to conduct uniaxial tensile test for BN-doped graphene nanoribbons (BN-GNR) with varying widths and BN atom fractions. Efforts are made to study the constitutive relations for the edges and the whole BN-GNR and explore the fracture mechanisms of the hybrid nanoribbons. The substantial softening effect of the edges induced by wrinkling alters the impact of BN concentration on the stiffness in the sub-20 nm regime deviating from the linear behaviour observed in the bulk case. Fracture properties are unexpectedly independent of BN concentrations unlike in the bulk and the failure behaviour is rather decided by the graphene ribbon edge structure. Here the armchair edges serve as the source of crack nucleation at an early stage leading to weakened strength and reduced stretchability, whereas zigzag edges do not promote early crack nucleation and leads to the size dependence of fracture properties.
What problem does this paper attempt to address?
The problems that this paper attempts to solve mainly focus on exploring the elastic modulus and fracture characteristics of boron nitride - doped graphene nanoribbons (BN - GNR) with different widths and h - BN atomic fractions, and studying the relationship between these properties and the edge structures. Specifically, the paper conducts uniaxial tensile tests through molecular dynamics simulations to explore the following issues:
1. **Size - dependence of elastic modulus and fracture characteristics**: Study how the elastic modulus of BN - GNR with different widths changes with the width, and whether the fracture characteristics are affected by the width.
2. **Effect of edge structure on mechanical properties**: Explore the influence of different edge structures (such as armchair - edged and zigzag - edged) on the elastic modulus and fracture behavior of BN - GNR.
3. **Effect of doping concentration on mechanical properties**: Analyze the mechanical properties of BN - GNR under different h - BN atomic fractions, especially the elastic modulus and fracture characteristics.
### Main findings
1. **Size - dependence of elastic modulus**:
- For the armchair - edged BN - GNR (BN - AGNR) doped with 20% h - BN, as the width increases from 5 nm to 20 nm, the elastic modulus increases from approximately 691 GPa to 793 GPa, and then gradually approaches the elastic modulus of the infinite BN - AGNR, which is 821 GPa.
- The elastic modulus of the zigzag - edged BN - GNR (BN - ZGNR) also shows a similar size - dependence, but its elastic modulus is approximately 20% higher than that of BN - AGNR as a whole.
2. **Edge - softening effect**:
- The edge atomic chains show significant bond contraction at 0 K, causing the bond length of the outer atomic chains to approach the length of the carbon triple bond, which may produce a hardening effect. However, at 300 K, due to thermal perturbations, the edge atomic chains form periodic wrinkles, and these wrinkles lead to a negative elastic modulus of the edge atomic chains, thus producing a softening effect.
- The softening effect of the armchair - edged is more significant, while the softening effect of the zigzag - edged is relatively weak.
3. **Edge influence on fracture characteristics**:
- The BN - GNR with armchair - edged shows lower fracture stress and strain at a finite width, mainly because cracks nucleate at the edge in the early stage and then expand inward along the zigzag direction.
- The BN - GNR with zigzag - edged shows the highest fracture stress and strain at a width of 5 nm, even higher than that of the infinite BN - GNR, but as the width increases, the fracture performance gradually decreases.
### Conclusion
The paper reveals the complex behavior of the elastic modulus and fracture characteristics of boron nitride - doped graphene nanoribbons at different widths and doping concentrations through molecular dynamics simulations. The edge structure has a significant impact on these mechanical properties, especially the softening effect of the armchair - edged and the early crack nucleation mechanism of the zigzag - edged. These findings provide an important theoretical basis for designing BN - GNR with specific mechanical properties.