Wärmeleitfähigkeit neuartiger 2D-Kohlenstoff-Allotrope mit 558-Ringen: Einblicke aus Molekulardynamik

Zusammenfassung

Diese Studie untersucht den Wärmetransport zweier neu vorgeschlagener zweidimensionaler Kohlenstoff-Allotrope – Octagon-Pentagon-Graphen-Line (OPG-L) und Octagon-Pentagon-Graphen-Zigzag (OPG-Z) – mittels Reverse- und Equilibrium-Molekulardynamik-Simulationen. Die Wärmeleitfähigkeit (TC) wächst stetig mit Probengröße; Extrapolation auf unendliche Länge ergibt Werte von 310-332 W m⁻¹ K⁻¹ für OPG-L und 247-228 W m⁻¹ K⁻¹ für OPG-Z. Im Vergleich zu Graphen (≈ 3000-5000 W m⁻¹ K⁻¹) stammt die niedrigere TC von reduzierten Phonon-Gruppengeschwindigkeiten und mittleren freien Wegen. Temperatur und Zugdehnung unterdrücken TC systematisch, was Tunbarkeit für Thermal-Management-Anwendungen bietet.

Einführung

Kohlenstoff-Allotrope wie Diamant, Kohlenstoffnanoröhren und Graphen sind für ihre außergewöhnliche Wärmeleitfähigkeit bekannt, wobei einwandige Nanoröhren und Graphen mehrere tausend W m⁻¹ K⁻¹ erreichen. Jüngste theoretische Arbeiten haben eine Vielzahl 2D-Kohlenstoff-Netzwerke enthüllt – Graphyne, Octagraphen, T-Graphen und H-Netz – jeweils mit eindeutigen Ringmotiven und vielversprechenden elektronischen oder mechanischen Eigenschaften. Su et al. führten zwei energetisch günstige, kinetisch stabile 2D-Blätter ein, bestehend aus Okta- und Pentagone (558 Ringe). Diese Strukturen, OPG-L und OPG-Z, besitzen Bildungsenergien von 0,31 und 0,34 eV Atom⁻¹, jeweils, substanziel niedriger als Graphyne (0,76 eV Atom⁻¹) und zeigen anisotropes elektronisches Verhalten. Ihre thermischen Dissipationseigenschaften bleiben jedoch unerforscht, was die vorliegende Untersuchung motiviert.

Modell und Methoden

Die atomaren Modelle von OPG-L (Abb. 1a) und OPG-Z (Abb. 1b) werden aus wiederholenden 558-Bändern aufgebaut, mit Chiralität analog zu Graphen (Armchair und Zigzag). LAMMPS wurde mit einem modifizierten optimierten Tersoff-Potenzial eingesetzt, um C-C-Wechselwirkungen genau zu erfassen. Die erhaltenen Gleichgewichts-Gitterkonstanten (OA = 3,63 Å, OB = 9,38 Å für OPG-L; OA = 6,78 Å, OB = 5,04 Å für OPG-Z) stimmen mit vorherigen First-Principles-Ergebnissen überein.

Reverse Non-Equilibrium Molecular Dynamics (rNEMD) wurde verwendet, um TC zu berechnen. Nach 0,25 ns Nosé-Hoover-Gleichgewicht bei 300 K wurde das Sample in 50 Scheiben entlang der Wärmeübertragungsrichtung partitioniert. Das heißeste Atom im kalten Slab und das kälteste Atom im heißen Slab tauschten kinetische Energie aus, erzeugend einen stationären Wärmefluss. Der resultierende Temperaturgradient liefert TC über Fouriersches Gesetz: κ = J / (2A ∂T/∂L), wo der Faktor 2 bidirektionalen Wärmefluss berücksichtigt und die Dicke als 0,34 nm, der Graphen-Zwischenschichtabstand, genommen wird.

Equilibrium MD (EMD) Berechnungen wurden mit dem Green-Kubo-Formalismus durchgeführt. Laufende TC wurde aus der Wärmefluss-Autokorrelationsfunktion erhalten, und unendliche-Größen-Werte wurden durch Mittelung über 1-2 ns Korrelationszeit extrahiert.

Ergebnisse und Diskussionen

Größenabhängigkeit: TC steigt monoton mit effektiver Probenlänge (L_eff = L/2). Für 50-nm OPG-L (Zigzag) und OPG-Z (Zigzag) betragen TC-Werte 125 und 94 W m⁻¹ K⁻¹, jeweils; bei 1000 nm steigen sie auf 296 und 236 W m⁻¹ K⁻¹. Lineare Fits von κ⁻¹ versus L_eff⁻¹ liefern extrapolierte unendliche-Größen-TC: 310-332 W m⁻¹ K⁻¹ für OPG-L und 247-228 W m⁻¹ K⁻¹ für OPG-Z, in exzellenter Übereinstimmung mit EMD-Ergebnissen (313-344 W m⁻¹ K⁻¹ für OPG-L, 261-233 W m⁻¹ K⁻¹ für OPG-Z).

Temperatureffekt: Über 200-500 K nimmt TC aufgrund verstärkter Umklapp-Streuung ab. Von 300 K auf 500 K fällt TC um 36-42 %. Die relativen Leitfähigkeiten folgen dem in Graphen beobachteten Trend, was die Universalität von Phonon-Phonon-Streuung in 2D-Gittern bestätigt.

Dehnungseffekt: Einachsige Zugdehnung entlang der Wärmeübertragungsrichtung unterdrückt TC in beiden Allotropen monoton. Bei 10 % Dehnung erreichen TC-Reduktionen 49 % (OPG-L Zigzag) und 31 % (OPG-Z Armchair). Vibrationsdichte-Zustands-Analyse zeigt deutliche Phonon-Erweichung in den In-Plane-Moden, was für den verminderten Wärmetransport verantwortlich ist.

Mechanische Eigenschaften: Unter einachsiger Spannung zeigen OPG-L und OPG-Z Young-Moduln von 538-648 GPa und 492-550 GPa, jeweils, vergleichbar mit Graphyne aber niedriger als makelloses Graphen. Ultimative Dehnungen erreichen 17,2 % (OPG-L Zigzag) und 8,7 % (OPG-L Armchair), was robustes, doch dehnungsempfindliches Verhalten anzeigt.

Schlussfolgerungen

Durch umfassende rNEMD und EMD Simulationen haben wir die Größen-, Temperatur- und Dehnungsabhängigkeit von TC in OPG-L und OPG-Z quantifiziert. Ihre Leitfähigkeiten sind markant niedriger als Graphens aufgrund reduzierter Phonon-Geschwindigkeiten und mittlerer freier Wege, bleiben aber durch Dehnungsengineering einstellbar. Diese Einsichten etablieren OPG-L und OPG-Z als vielversprechende Kandidaten für maßgeschneidertes Thermal-Management in Mikro- und Nanoskalen-Geräten.

Abkürzungen

558

Fünf-fünf-acht-gliedrige Ringe

OPG-L

Okta- und Pentagone Graphen-Line

OPG-Z

Okta- und Pentagone Graphen-Zigzag

rNEMD

Reverse Non-Equilibrium Molecular Dynamics

TC

Wärmeleitfähigkeit

VDOS

Vibrationsdichte der Zustände