Conductividad térmica de nuevos alótropos carbono 2D con anillos 558: perspectivas desde dinámica molecular

Resumen

Este estudio investiga transporte térmico dos alótropos carbono bidimensionales recién propuestos —grafeno-línea octágono-pentágono (OPG-L) y grafeno-zigzag octágono-pentágono (OPG-Z)— usando simulaciones dinámica molecular inversa y equilibrio. La conductividad térmica (TC) crece constantemente con tamaño muestra; extrapolación longitud infinita produce valores 310–332 W m⁻¹ K⁻¹ OPG-L y 247–228 W m⁻¹ K⁻¹ OPG-Z. Comparado grafeno (≈ 3000–5000 W m⁻¹ K⁻¹), menor TC origina velocidades grupo fonón y caminos libres medios reducidos. Temperatura y tensión tensil suprimen sistemáticamente TC, ofreciendo ajustabilidad para aplicaciones gestión térmica.

Introducción

Alótropos carbono como diamante, nanotubos carbono y grafeno son celebrados por su conductividad térmica excepcional, con nanotubos pared simple y grafeno alcanzando varios miles W m⁻¹ K⁻¹. Trabajo teórico reciente ha revelado multitud redes carbono 2D —grafino, octagrafeno, T-grafeno y H-net— cada con motivos anillos distintos y prometedoras propiedades electrónicas o mecánicas. Su et al. introdujeron dos láminas 2D energéticamente favorables, cinéticamente estables compuestas octágonos y pentágonos (anillos 558). Estas estructuras, OPG-L y OPG-Z, poseen energías formación 0.31 y 0.34 eV átomo⁻¹ respectivamente, sustancialmente menores grafino (0.76 eV átomo⁻¹) y exhiben comportamiento electrónico anisotrópico. Sin embargo, sus características disipación térmica permanecen inexploradas, motivando presente investigación.

Modelo y Métodos

Modelos atómicos OPG-L (Fig. 1a) y OPG-Z (Fig. 1b) construidos desde cintas 558 repetitivas, quiralidad definida análogamente grafeno (brazo-silla y zigzag). LAMMPS empleado con potencial Tersoff optimizado modificado capturar interacciones C-C accuradamente. Constantes red equilibrio obtenidas (OA = 3.63 Å, OB = 9.38 Å OPG-L; OA = 6.78 Å, OB = 5.04 Å OPG-Z) concuerdan previos resultados primeros principios.

Dinámica molecular no-equilibrio inversa (rNEMD) usada calcular TC. Tras 0.25 ns equilibración Nosé-Hoover a 300 K, muestra particionada 50 láminas dirección transferencia calor. Átomo más caliente lámina fría y más frío lámina caliente intercambiaron energía cinética, generando flujo calor estable. Gradiente temperatura resultante produce TC vía ley Fourier: \[ \kappa = \frac{J}{2A\,\partial T/\partial L} \] donde factor 2 cuenta flujo calor bi-direccional y espesor tomado 0.34 nm, espaciado intercapa grafeno.

Cálculos MD equilibrio (EMD) realizados usando formalismo Green-Kubo. TC corriendo obtenida función autocorrelación flujo calor, valores tamaño-infinito extraídos promediando 1–2 ns tiempo correlación.

Resultados y Discusiones

Dependencia tamaño: TC aumenta monótonamente longitud efectiva muestra (L_eff = L/2). Para OPG-L (zigzag) y OPG-Z (zigzag) 50-nm, valores TC 125 y 94 W m⁻¹ K⁻¹ respectivamente; a 1000 nm suben a 296 y 236 W m⁻¹ K⁻¹. Ajustes lineales κ⁻¹ versus L_eff⁻¹ proveen TC tamaño-infinito extrapolado: 310–332 W m⁻¹ K⁻¹ OPG-L y 247–228 W m⁻¹ K⁻¹ OPG-Z, en excelente acuerdo resultados EMD (313–344 W m⁻¹ K⁻¹ OPG-L, 261–233 W m⁻¹ K⁻¹ OPG-Z).

Efecto temperatura: En 200–500 K, TC disminuye debido dispersión Umklapp aumentada. De 300 K a 500 K, TC cae 36–42%. Conductividades relativas siguen tendencia observada grafeno, confirmando universalidad dispersión fonón-fonón redes 2D.

Efecto tensión: Tensión tensil uniaxial aplicada dirección transferencia calor suprime monótonamente TC ambos alótropos. A 10% tensión, reducciones TC alcanzan 49% (OPG-L zigzag) y 31% (OPG-Z brazo-silla). Análisis densidad estados vibracionales revela ablandamiento fonón pronunciado modos en-plano, explicando transporte térmico disminuido.

Propiedades mecánicas: Bajo tensión uniaxial, OPG-L y OPG-Z exhiben módulos Young 538–648 GPa y 492–550 GPa respectivamente, comparables grafino pero menores grafeno prístino. Deformaciones últimas alcanzan 17.2% (OPG-L zigzag) y 8.7% (OPG-L brazo-silla), indicando comportamiento robusto pero sensible tensión.

Conclusiones

Mediante simulaciones rNEMD y EMD exhaustivas, hemos cuantificado dependencia tamaño, temperatura y tensión TC en OPG-L y OPG-Z. Sus conductividades son marcadamente menores grafeno debido velocidades fonón y caminos libres medios reducidos, yet permanecen ajustables via ingeniería tensión. Estos insights establecen OPG-L y OPG-Z como candidatos prometedores gestión térmica adaptada dispositivos micro- y nanoescala.

Abreviaturas

558

Anillos cinco-cinco-ocho miembros

OPG-L

Octágono y pentágono grafeno-línea

OPG-Z

Octágono y pentágono grafeno-zigzag

rNEMD

Dinámica molecular no-equilibrio inversa

TC

Conductividad térmica

VDOS

Densidad estados vibracionales