🔬 physics

Tuning the Interfacial Thermal Conductance of Graphene–C 3 B Heterostructures via Interfacial Engineering: A Molecular Dynamics Study

This molecular dynamics study reveals that interfacial thermal conductance in graphene–C3B heterostructures is maximized at sharp interfaces and significantly suppressed by interfacial mixing due to enhanced phonon scattering and localized vibrational modes, while the heterostructure length has a relatively minor impact.

Original authors: Mohammad Zakeri, Mohammad Moein Kharrazi, Erfan Derakhshan, Ehsan Taheran, Ali Rajabpour

Published 2026-07-07

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Mohammad Zakeri, Mohammad Moein Kharrazi, Erfan Derakhshan, Ehsan Taheran, Ali Rajabpour

Original paper licensed under CC BY 4.0 (https://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

Imagine you have two different types of high-speed train tracks laid side-by-side. One track is made of pure, smooth graphene (a super-strong, single-layer carbon material), and the other is made of a newer material called C3B (a mix of carbon and boron). When you want to move heat energy from one track to the other, the place where they meet—the "interface"—is critical.

This paper is like a detailed engineering report on what happens at that meeting point. The researchers used powerful computer simulations to act as a "molecular microscope," watching how heat travels when they mess with the boundary between these two materials.

Here is the breakdown of their findings using simple analogies:

1. The "Clean Cut" vs. The "Messy Blend"

The researchers tested different ways to join these two tracks.

The Sharp Interface (0% mixing): Imagine the two tracks meeting with a perfectly clean, straight cut. No atoms from one side are on the other side.
The Mixed Interface: Imagine taking a blender and mixing the atoms from both tracks together in the middle section. They tested mixing ratios ranging from a tiny bit of blending to a complete 50/50 mash-up.

2. The Surprising "U-Shape" Result

The most important discovery is that the amount of heat flowing across the junction didn't just go down steadily as they mixed more atoms. Instead, it behaved like a U-shaped curve:

At the top (Clean Cut): Heat flows very well. The tracks are distinct, and the connection is efficient.
In the middle (40–60% mixing): This is the "traffic jam." When they mixed the atoms partially, heat flow dropped to its lowest point. It's like trying to drive a car through a construction zone where the road surface is constantly changing from asphalt to gravel to dirt. The heat energy (the car) gets confused, scatters, and gets stuck.
At the other top (100% mixing): Surprisingly, if they mixed the atoms completely so the whole interface became a uniform blend of the new material, the heat flow started to recover. It wasn't as good as the clean cut, but it was much better than the messy middle ground.

The Analogy: Think of it like a conversation between two people speaking different languages.

Sharp Interface: They use a perfect translator. Communication is fast.
Partial Mixing: They try to speak a broken, mixed-up version of both languages. It's chaotic, full of misunderstandings, and communication slows to a crawl.
Full Mixing: They both agree to speak a new, third language together. It's not as fast as the perfect translator, but it's much smoother than the chaotic broken version.

3. Why Does the "Messy Middle" Fail?

The paper explains this using phonons. In simple terms, phonons are tiny packets of sound or vibration that carry heat through solid materials.

In a clean interface, these vibration packets glide across easily.
In the "messy middle" (40–60% mixing), the random mix of Carbon and Boron atoms creates a bumpy, disordered landscape. The vibration packets hit these bumps, scatter in random directions, or get trapped in little pockets near the interface. The paper calls this "phonon scattering." It's like a pinball machine where the ball (heat) gets stuck bouncing around the middle instead of reaching the other side.

4. Does the Pattern of Mixing Matter?

The researchers also asked: "Does it matter how we mix them? Is it better to have a gradient (slowly changing from one side to the other) or a random mix?"

The Answer: Not really. As long as the total amount of mixing is the same, the specific pattern (whether it's a smooth gradient or a random scatter) didn't change the heat flow much.
The Takeaway: It's the total chaos (the disorder) that kills the heat flow, not the specific shape of that chaos.

5. Length Doesn't Matter Much

They also checked if making the tracks longer changed the results. They found that while length has a tiny effect, the mixing ratio is the boss. Whether the tracks are short or long, a messy middle section will always block heat much more than a clean or fully blended section.

Summary

The main lesson from this study is that for these specific materials (Graphene and C3B), keeping the boundary clean and sharp is the best way to move heat. If you must mix them, don't do it halfway; either keep them separate or blend them completely. The "in-between" state creates a traffic jam for heat that is very hard to overcome.

Technical Summary: Tuning the Interfacial Thermal Conductance of Graphene–C $_3$ B Heterostructures via Interfacial Engineering

Problem Statement
Two-dimensional (2D) heterostructures, particularly those combining graphene with semiconducting materials like C $_3$ B, are promising candidates for next-generation nanoelectronic and energy devices. However, practical fabrication inevitably introduces an interfacial region where atomic intermixing occurs between the adjacent layers. While the impact of interfaces on thermal transport is recognized, the specific role of interfacial atomic mixing in graphene–C $_3$ B systems remains insufficiently understood. Specifically, the coupled effects of mixing ratio, the spatial distribution of mixed atoms, and finite-length effects on Interfacial Thermal Conductance (ITC) have not been systematically quantified. This study addresses this gap to determine how interfacial disorder governs phonon-mediated heat transport in these carbon-based heterostructures.

Methodology
The authors employed Nonequilibrium Molecular Dynamics (NEMD) simulations using the LAMMPS package with the Tersoff potential to model heat transport across graphene–C $_3$ B heterostructures.

System Configuration: Simulations utilized a 4 nm-thick interfacial region with varying atomic mixing ratios ( $m$ ) defined as the fraction of C $_3$ B atoms within the graphene lattice. Ratios ranged from 0% (sharp interface) to 100% (fully mixed).
Thermal Transport Calculation: A steady-state heat flux was induced by coupling the ends of the structure to hot (320 K) and cold (280 K) reservoirs. ITC was calculated based on the temperature discontinuity ( $\Delta T$ ) across the interface and the heat flux ( $q$ ).
Finite-Size Correction: To isolate intrinsic interfacial properties, simulations were performed for heterostructure lengths ( $L$ ) of 30, 60, 90, 120, and 150 nm. The results were extrapolated to the length-independent limit ( $h_\infty$ ) using a linear relation between inverse conductance and inverse length.
Mechanistic Analysis: To elucidate transport mechanisms, the study utilized Phonon Density of States (PDOS) analysis to examine vibrational modes and Spectral Heat Current (SHC) analysis (via GPUMD) to resolve heat flux by angular frequency.

Key Contributions and Results

Non-Monotonic Dependence on Mixing Ratio: The study reveals a pronounced non-monotonic relationship between ITC and the mixing ratio. The ITC is highest for a sharp interface ( $m=0\%$ ). As mixing increases, ITC decreases significantly, reaching a minimum at intermediate mixing ratios of approximately 40–60%. Surprisingly, further increasing the mixing toward 100% results in a partial recovery of ITC, creating a U-shaped trend.
Dominance of Disorder over Spatial Distribution: While the presence of mixing drastically reduces ITC compared to a sharp interface, the specific spatial arrangement of mixed atoms (e.g., uniform vs. gradient distributions) within the fixed 4 nm thickness has only a marginal effect on the thermal conductance. The total degree of atomic disorder is the dominant factor, rather than the morphology of the mixing.
Length Independence: The effect of heterostructure length on ITC is relatively weak compared to the impact of interfacial mixing. Extrapolation confirms that the observed trends are intrinsic to the interface rather than artifacts of finite-size effects.
Mechanistic Insights via PDOS and SHC:
- Phonon Scattering and Localization: PDOS analysis indicates that intermediate mixing ratios induce pronounced phonon peak broadening and the emergence of interface-localized vibrational modes. These features, particularly in the low-frequency range, act as scattering centers that trap energy and impede heat flow.
- Spectral Mismatch: The sharp interface and the fully mixed interface ( $m=100\%$ ) exhibit better vibrational spectral matching, particularly in the low-frequency regime, facilitating more efficient phonon transmission.
- Frequency-Resolved Suppression: SHC analysis demonstrates that partially mixed interfaces suppress heat flux primarily across the mid- and high-frequency ranges (short-wavelength phonons). Long-wavelength acoustic phonons remain relatively robust against interfacial mixing, but the overall reduction in transport is driven by the disruption of short-wavelength channels.

Significance and Claims
The paper claims that for graphene-based heterostructures, preserving a sharp interface is more favorable for efficient thermal transport than introducing partial atomic mixing. The findings suggest that interfacial disorder plays a dominant role in limiting phonon transmission, primarily by enhancing scattering and inducing localized modes that obstruct heat transfer.

The study provides practical design guidelines for thermal management in 2D heterostructures, emphasizing that the overall degree of atomic disorder is a more critical parameter than the specific spatial gradient of mixing. By clarifying the non-monotonic behavior of ITC and the underlying phonon mechanisms, the work offers a clearer physical understanding of how interfacial engineering can be used to tune thermal properties, highlighting that intermediate mixing levels are particularly detrimental to thermal performance.