🔬 materials science

Irida-Graphene Phonon Thermal Transport via Non-equilibrium Molecular Dynamics Simulations

This study utilizes non-equilibrium molecular dynamics simulations to reveal that the newly proposed 2D carbon allotrope, Irida-Graphene, exhibits an intrinsic room-temperature thermal conductivity of approximately 215 W/mK—significantly lower than pristine graphene due to phonon scattering from its porous 3-6-8 ring structure and reduced group velocities—while maintaining isotropic thermal transport with notable size effects.

Original authors: Isaac M. Felix, Raphael M. Tromer, Leonardo D. Machado, Douglas S. Galvão, Luiz A. Ribeiro, Marcelo L. Pereira

Published 2026-09-09

📖 3 min read☕ Coffee break read

CC BY 4.0

Original authors: Isaac M. Felix, Raphael M. Tromer, Leonardo D. Machado, Douglas S. Galvão, Luiz A. Ribeiro, Marcelo L. Pereira

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

Heat is a constant challenge in the world of tiny machines. As electronic devices shrink to the size of a grain of sand, the heat they generate becomes harder to manage, often threatening to melt the very components they are meant to power. To solve this, scientists look to the realm of two-dimensional materials—sheets of atoms so thin they are essentially flat. Among these, a single layer of carbon atoms arranged in a honeycomb pattern, known as graphene, has long been the star. It conducts heat with incredible efficiency, acting like a superhighway for thermal energy. However, researchers are constantly searching for new variations of these materials that might offer different properties, perhaps trading some of that extreme speed for other useful characteristics, or simply finding new ways to control how heat moves through a solid.

In a recent study, a team of researchers turned their attention to a newly proposed carbon structure called Irida-Graphene. Unlike the perfect honeycomb of graphene, this material is built from a mix of triangles, hexagons, and octagons, creating a sheet with a porous, open design. The scientists wanted to understand how heat travels through this unique lattice. Using powerful computer simulations that mimic the movement of atoms, they tracked how thermal energy flows across the material. They found that while Irida-Graphene is still a good conductor of heat, it is significantly slower at the job than its famous cousin. At room temperature, the material conducts heat at a rate of approximately 215 watts per meter per kelvin. This is a substantial drop from the roughly 1,200 watts per meter per kelvin observed in pristine graphene when measured with the same methods.

The reason for this slowdown lies in the material's own architecture. In a perfect sheet like graphene, heat-carrying vibrations, known as phonons, can travel long distances without interruption. In Irida-Graphene, the mix of different ring shapes creates a bumpy landscape. As these vibrations move through the material, they constantly bump into the structural irregularities of the porous design. This scattering effect slows them down, much like a runner navigating a path filled with obstacles rather than a clear track. The researchers also observed that the speed at which these vibrations travel is lower in Irida-Graphene than in graphene, further contributing to the reduced ability to move heat.

The study also revealed that the size of the material matters greatly. When the sheet is very small, heat moves in a direct, unimpeded way. As the sheet grows longer, the heat begins to scatter more, shifting into a different mode of transport. This transition happens over a specific distance, which the team calculated to be around 55 nanometers. Despite being less conductive than graphene, the material still performs better than many other two-dimensional substances, such as certain forms of boron or molybdenum disulfide. This suggests that Irida-Graphene could still find a home in thermal management systems where a moderate, controlled flow of heat is more desirable than the extreme speed of graphene, offering a new tool for engineers designing the next generation of compact electronics.

Technical Summary: Irida-Graphene Phonon Thermal Transport via Non-equilibrium Molecular Dynamics Simulations

Problem Statement
The thermal management of nanomaterials is critical for the advancement of nanoelectronics and thermoelectric devices. While pristine graphene is renowned for its exceptional thermal conductivity, the discovery and characterization of new two-dimensional (2D) carbon allotropes with diverse topologies are essential for expanding the material library for specific applications. Recently, a new 2D carbon allotrope named Irida-Graphene (Irida-G) was proposed, characterized by a flat sheet topology composed of 3-, 6-, and 8-membered carbon rings. Despite its structural stability and metallic properties, the thermal transport characteristics of Irida-G remained unexplored. Understanding these properties is necessary to determine its viability for thermal management applications and to understand how its unique porous structure influences phonon scattering compared to standard graphene.

Methodology
The authors employed classical reactive molecular dynamics (MD) simulations using the LAMMPS software package. Interatomic forces were modeled using the second-generation Reactive Empirical Bond Order (REBO) potential, which was validated for Irida-G by computing phonon dispersions via the GULP software; the absence of imaginary frequencies confirmed the structural stability of the model.

To investigate thermal transport, the Reverse Non-Equilibrium Molecular Dynamics (RNEMD) method was utilized. The simulation setup involved:

System Configuration: Periodic boundary conditions in the $xy$-plane and free boundaries in the $z$ -direction. Systems ranged in size from 13,824 to over 1 million atoms, with lengths varying from 10 nm to 800 nm.
Heat Flux Induction: The system was partitioned into slabs. Kinetic energy was exchanged between "hot" and "cold" regions to establish a steady-state temperature gradient.
Analysis: Thermal conductivity ( $\kappa$ ) was calculated using Fourier's law based on the time-averaged heat flux and temperature gradient. The study also analyzed the vibrational density of states (VDOS) via the velocity autocorrelation function and calculated phonon group velocities from dispersion relations to elucidate the mechanisms of heat transport.

Key Contributions and Results

Intrinsic Thermal Conductivity: The simulations revealed that Irida-G possesses an intrinsic thermal conductivity ( $\kappa_\infty$ ) of approximately 130 W/mK at room temperature when calculated directly from the REBO potential. However, acknowledging that classical force fields often underestimate experimental values (as seen with graphene, where the REBO potential yields ~1200 W/mK versus experimental values of 1500–5000 W/mK), the authors applied a semi-empirical correction method. This correction, scaling the result based on graphene's known deviation, suggests a "true" intrinsic thermal conductivity for Irida-G of approximately 215 W/mK.
Size Effects and Transport Regimes: The study identified three distinct heat transport regimes dependent on system length ( $L$ $L$ ):
- Ballistic Regime ( $L < 30$ nm): Thermal conductivity scales linearly with length.
- Ballistic-Diffusive Transition ( $30 < L < 100$ nm): A transition region where dependence on length decreases.
- Diffusive Regime ( $L > 100$ nm): Thermal conductivity becomes largely independent of length.
  The effective phonon mean free path was estimated to be approximately 55 nm.
Mechanism of Reduction: The thermal conductivity of Irida-G is significantly lower than that of pristine graphene (approx. 1200 W/mK in the same simulation framework). This reduction is attributed to enhanced phonon scattering caused by the material's porous structure (the arrangement of triangles and octagons).
Phonon Dynamics:
- Group Velocity: The average phonon group velocity for Irida-G was calculated at 410 m/s, significantly lower than graphene's 780 m/s. Since lattice thermal conductivity is proportional to the square of the group velocity, this difference is a primary factor in the reduced conductivity.
- Vibrational Spectrum: The VDOS analysis indicated that flexural (ZA) phonons, particularly at low frequencies (below 10 THz), are the primary heat carriers. These modes exhibit longer wavelengths and fewer scatterings, similar to trends observed in other 2D materials.
Isotropy: The thermal transport properties were found to be practically isotropic, with negligible differences between the $x$ and $y$ directions.

Significance and Claims
The paper claims that while Irida-G exhibits lower thermal conductivity than pristine graphene, its value (~215 W/mK) remains competitive with other known 2D materials such as hexagonal boron nitride (h-BN), PHA-graphene, and biphenylene networks. The study establishes that the unique topological arrangement of Irida-G introduces significant phonon scattering, effectively tuning its thermal properties. The authors conclude that Irida-G holds potential for specific thermal management applications where moderate thermal conductivity is required, sharing characteristics with other emerging 2D materials. The work provides a foundational understanding of how structural porosity in carbon allotropes influences thermal transport, bridging the gap between theoretical proposals and functional material characterization.

More like this