🔬 materials science

On the Thermomechanical Properties and Fracture Patterns of the Novel Nonbenzenoid Carbon Allotrope (Biphenylene Network): A Reactive Molecular Dynamics Study

This study employs reactive molecular dynamics simulations to reveal that the novel biphenylene network (BPN) exhibits exceptional thermomechanical properties comparable to graphene, including a Young's modulus of ~1019.4 GPa and a melting point of ~4024 K, while demonstrating a unique, direction-dependent inelastic fracture process involving four distinct morphological transitions that are suppressed by the presence of nanocracks.

Original authors: M. L. Pereira Júnior, W. F. da Cunha, R. T. de Sousa Junior, G. D. Amvame Nze, D. S. Galvão, L. A. Ribeiro Júnior

Published 2026-09-09

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: M. L. Pereira Júnior, W. F. da Cunha, R. T. de Sousa Junior, G. D. Amvame Nze, D. S. Galvão, L. A. Ribeiro Júnior

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Carbon has long been the star of materials science, famous for its ability to form incredibly strong and versatile structures. For decades, the most celebrated member of this family was graphene, a single layer of carbon atoms arranged in a perfect honeycomb pattern. This material is renowned for being nearly unbreakable and an excellent conductor of electricity, leading scientists to believe it could revolutionize everything from electronics to energy storage. However, graphene has a significant limitation: it lacks a natural energy gap, which makes it difficult to use in many types of computer chips and light-based devices. This challenge has driven researchers to look for new ways to arrange carbon atoms, hoping to create materials that keep graphene's strength while gaining new electronic abilities. Recently, a new two-dimensional carbon structure called the biphenylene network was successfully created in a laboratory, offering a fresh alternative to the classic honeycomb design.

This new material looks quite different from graphene. Instead of a uniform grid of six-sided rings, it is a mosaic made of four-sided, six-sided, and eight-sided rings of carbon atoms fused together. While scientists had already studied how this new structure behaves electronically, they did not yet understand how it would hold up under physical stress or heat. To answer these questions, a team of researchers used powerful computer simulations to watch how the material reacts when pulled apart or heated up. They treated the atoms as individual particles that could break and reform their connections, allowing them to observe the exact moment the material fails. Their goal was to see if this new carbon form was as tough as graphene and to understand how tiny flaws, which are common in real-world manufacturing, would affect its strength.

The researchers began by simulating the material being stretched in different directions, much like pulling on a piece of fabric. They found that the biphenylene network behaves in a surprisingly complex way before it finally breaks. Unlike graphene, which snaps abruptly once it reaches its limit, this new material goes through a series of transformations. As it is pulled, the atoms rearrange themselves, shifting from one shape to another. The material passes through several distinct stages where the rings of atoms change their configuration, effectively softening the structure and allowing it to stretch further without breaking immediately. In some cases, the material even temporarily turns into a structure that looks very much like graphene before it eventually tears apart. This ability to change shape and absorb energy suggests that the material has a unique resilience that is different from the brittle nature of graphene.

The study also looked at what happens when the material has small cracks or missing pieces, which are inevitable in any real-world application. The researchers simulated cracks running both parallel and perpendicular to the direction of the pull. They discovered that the orientation of the crack matters greatly. When a crack runs across the direction of the pull, the material becomes weaker and breaks sooner. However, when the crack runs along the direction of the pull, the atoms near the crack can actually move closer together and reconnect, allowing the material to hold its strength better than expected. This means that the material's durability depends heavily on how the flaw is aligned relative to the force being applied.

Beyond pulling, the team also tested how the material handles extreme heat. They simulated heating the structure from room temperature up to ten thousand degrees Kelvin to see when it would melt. The results showed that the material remains solid and stable up to a very high temperature, with a melting point of 4024 Kelvin. This is remarkably close to the melting point of graphene, indicating that the new material is just as thermally stable as its famous cousin. Even as it heats up, the material undergoes structural changes, eventually breaking down into a gas-like state of individual atoms and short chains.

The simulations provided specific numbers for how stiff and strong the material is. The researchers calculated a value for its stiffness, known as Young's modulus, which ranged from about 570 to 1019 gigapascals depending on the direction of the pull and the presence of defects. These values are comparable to the stiffness of graphene, confirming that the new material is indeed a heavyweight contender in the world of carbon-based materials. The study concluded that while the biphenylene network shares the exceptional strength and heat resistance of graphene, it possesses a distinct personality. It does not simply snap under pressure; instead, it bends, reshapes, and transforms in a series of steps before failing. This behavior, combined with its ability to maintain strength even when flawed in certain ways, suggests that this new carbon allotrope could be a promising candidate for future technologies, provided engineers can learn to work with its unique way of handling stress.

Technical Summary: Thermomechanical Properties and Fracture Patterns of the Biphenylene Network (BPN)

Problem Statement
While graphene has established itself as a paradigmatic two-dimensional carbon allotrope, its application in optoelectronics is limited by a zero bandgap. Recently, a new nonbenzenoid carbon allotrope, the biphenylene network (BPN), was experimentally synthesized. BPN consists of a periodic arrangement of fused four-, six-, and eight-membered carbon rings. Although theoretical studies have explored its electronic structure and potential as a battery material, detailed investigations into its thermomechanical stability, fracture mechanisms, and the influence of lattice defects (specifically nanocracks) remain absent from the literature. This study addresses the gap in understanding the mechanical response and failure dynamics of pristine and defective BPN under uniaxial tensile loading and thermal stress.

Methodology
The authors employed fully atomistic reactive molecular dynamics (MD) simulations using the ReaxFF force field, implemented within the LAMMPS software package. This approach allows for the dynamic formation and breaking of chemical bonds, which is essential for modeling fracture mechanisms.

Models: Three structural configurations were simulated: non-defective (pristine) BPN, BPN with a horizontally aligned nanocrack (Defect-H), and BPN with a vertically aligned nanocrack (Defect-V).
Simulation Protocol: Structures were equilibrated at 300 K under an NPT ensemble. Tensile strain was applied uniaxially along both the x-direction (parallel to the horizontal crack) and y-direction (parallel to the vertical crack) at a strain rate of $1.0 \times 10^{-6} \text{ fs}^{-1}$ up to 50% strain.
Analysis: The study utilized von Mises stress ( $\sigma_{VM}$ ) per-atom calculations to identify fracture initiation sites and stress accumulation. Elastic properties (Young's modulus, fracture strain, ultimate strength) were derived from stress-strain curves. Thermal stability was assessed via heating ramp simulations from 300 K to 10,000 K to determine melting behavior.

Key Results

Mechanical Response and Fracture Dynamics:
- Distinct Inelastic Regimes: Unlike graphene, which transitions directly from an elastic to a brittle regime after a critical strain, BPN exhibits distinct intermediate inelastic stages. Under uniaxial tensile loading, non-defective BPN is converted into four distinct morphologies before fracture starts. This conversion process is dependent on the stretching direction. Some of the formed structures are mainly composed of eight-membered rings, which have different shapes in each morphology. Importantly, in the presence of nanocracks, no new morphologies are formed.
- Graphitization: In one of the morphologies formed by non-defective BPN, a graphitization process occurs, where the structure transforms into a graphene-like configuration prior to total fracture.
- Fracture Patterns: The fracture process is characterized by the formation of linear atomic chains (LACs) leading to complete separation.
- Defect Influence:
  - X-direction Stretching: The vertical defect (Defect-V) significantly weakens the lattice, reducing the fracture strain. Conversely, the horizontal defect (Defect-H) allows for bond reconstructions within the nanocrack as atoms move closer, slightly increasing the strain rate and resilience compared to the pristine structure.
  - Y-direction Stretching: The fracture strain is generally lower (approx. 7%) compared to x-direction stretching. Stress accumulation occurs in four- and six-membered rings rather than bonds parallel to the strain. Defects here act as softening agents, particularly the horizontal defect, which reduces the ultimate strength.
- Elastic Properties: The calculated Young's modulus ranges from 570.9 to 1019.4 GPa, and the ultimate strength varies between 46.9 and 64.4 GPa depending on the defect orientation and loading direction. These values are comparable to those reported for graphene.
Thermal Stability:
- Melting Point: The BPN melting point was determined to be 4024 K, identified by a pronounced peak in heat capacity ( $C_V$ ) and a discontinuity in total energy. This value is comparable to pristine graphene (4095 K) and amorphous graphene (3626 K).
- Melting Process: The melting process involves three regimes: structural integrity maintenance (300–1000 K), structural phase transitions and graphitization (1100–3900 K), and the transition to a gas-like phase (3900–5000 K). Complete melting occurs around 5000 K.

Significance and Claims
The paper claims that the Biphenylene Network shares exceptional mechanical properties with graphene, specifically regarding its Young's modulus and melting point, despite its distinct nonbenzenoid topology. The primary contribution of this work is the detailed mapping of BPN's unique fracture dynamics, which differ fundamentally from graphene by exhibiting multiple inelastic stages and structural morphological changes rather than immediate brittle failure.

The authors emphasize that the presence and orientation of nanocracks critically alter the mechanical response, either softening the structure or facilitating bond reconstructions that enhance resilience. The observation of a graphitization process prior to fracture highlights a unique structural evolution pathway for this material. The study concludes that BPN represents a promising 2D carbon material with thermomechanical stability comparable to graphene, warranting further investigation into its structural and electronic properties for potential applications.

More like this