Condutividade Térmica de Novos Alotropos 2D de Carbono com Anéis 558: Insights de Dinâmica Molecular

Resumo

Este estudo investiga transporte térmico de dois novos alotropos bidimensionais de carbono propostos — grafeno-linha octágono-pentágono (OPG-L) e grafeno-ziguezague octágono-pentágono (OPG-Z) — usando simulações de dinâmica molecular reversa e de equilíbrio. Condutividade térmica (CT) cresce firmemente com tamanho da amostra; extrapolação para comprimento infinito rende valores de 310-332 W m⁻¹ K⁻¹ para OPG-L e 247-228 W m⁻¹ K⁻¹ para OPG-Z. Comparado com grafeno (≈3000-5000 W m⁻¹ K⁻¹), menor CT origina-se de velocidades de grupo de fonons e caminhos livres médios reduzidos. Temperatura e tensão de tração suprimem sistematicamente CT, oferecendo ajustabilidade para aplicações de gestão térmica.

Introdução

Alotropos de carbono como diamante, nanotubos de carbono e grafeno são celebrados por sua condutividade térmica excepcional, com nanotubos de parede única e grafeno atingindo vários milhares de W m⁻¹ K⁻¹. Trabalhos teóricos recentes revelaram uma série de redes 2D de carbono — grafinas, octagrafeno, T-grafeno e H-net — cada uma com motivos de anéis distintos e propriedades eletrônicas ou mecânicas promissoras. Su et al. introduziram duas folhas 2D energeticamente favoráveis, cineticamente estáveis compostas de octágonos e pentágonos (anéis 558). Estas estruturas, OPG-L e OPG-Z, possuem energias de formação de 0,31 e 0,34 eV átomo⁻¹, respectivamente, substancialmente menores que grafinas (0,76 eV átomo⁻¹) e exibem comportamento eletrônico anisotrópico. No entanto, suas características de dissipação térmica permanecem inexploradas, motivando presente investigação.

Modelo e Métodos

Modelos atômicos de OPG-L (Fig. 1a) e OPG-Z (Fig. 1b) são construídos a partir de fitas 558 repetidas, com quiralidade definida analogamente a grafeno (poltrona e ziguezague). LAMMPS foi empregado com potencial Tersoff otimizado modificado para capturar interações C-C com precisão. Constantes de rede de equilíbrio obtidas (OA = 3,63 Å, OB = 9,38 Å para OPG-L; OA = 6,78 Å, OB = 5,04 Å para OPG-Z) concordam com resultados prévios de primeiros princípios.

Dinâmica molecular não-equilíbrio reversa (rNEMD) foi usada para calcular CT. Após 0,25 ns de equilíbrio Nosé-Hoover a 300 K, amostra foi particionada em 50 lâminas ao longo da direção de transferência de calor. Átomo mais quente na lâmina fria e átomo mais frio na lâmina quente trocaram energia cinética, gerando fluxo de calor estável. Gradiente de temperatura resultante rende CT via lei de Fourier: \[ \kappa = \frac{J}{2A\,\partial T/\partial L} \] onde fator 2 conta para fluxo de calor bi-direcional e espessura tomada como 0,34 nm, espaçamento intercamada do grafeno.

Cálculos MD de equilíbrio (EMD) realizados usando formalismo Green-Kubo. CT corrente obtida da função de autocorrelação de fluxo de calor, e valores de tamanho infinito extraídos por média sobre 1-2 ns de tempo de correlação.

Resultados e Discussões

Dependência de tamanho: CT aumenta monoticamente com comprimento efetivo da amostra (L_eff = L/2). Para OPG-L (ziguezague) e OPG-Z (ziguezague) de 50 nm, valores de CT são 125 e 94 W m⁻¹ K⁻¹, respectivamente; a 1000 nm sobem para 296 e 236 W m⁻¹ K⁻¹. Ajustes lineares de κ⁻¹ versus L_eff⁻¹ fornecem CT de tamanho infinito extrapolado: 310-332 W m⁻¹ K⁻¹ para OPG-L e 247-228 W m⁻¹ K⁻¹ para OPG-Z, em excelente concordância com resultados EMD (313-344 W m⁻¹ K⁻¹ para OPG-L, 261-233 W m⁻¹ K⁻¹ para OPG-Z).

Efeito de temperatura: Em 200-500 K, CT diminui devido a espalhamento Umklapp aprimorado. De 300 K a 500 K, CT cai 36-42%. Condutividades relativas seguem tendência observada em grafeno, confirmando universalidade de espalhamento fonon-fonon em redes 2D.

Efeito de tensão: Tensão de tração uniaxial aplicada junto à direção de transferência de calor suprime monoticamente CT em ambos alotropos. A 10% tensão, reduções de CT atingem 49% (OPG-L ziguezague) e 31% (OPG-Z poltrona). Análise de densidade de estados vibracionais revela amolecimento pronunciado de fonons em modos no plano, explicando transporte térmico diminuído.

Propriedades mecânicas: Sob tensão uniaxial, OPG-L e OPG-Z exibem módulos de Young de 538-648 GPa e 492-550 GPa, respectivamente, comparáveis a grafinas mas menores que grafeno pristino. Deformações últimas atingem 17,2% (OPG-L ziguezague) e 8,7% (OPG-L poltrona), indicando comportamento robusto mas sensível a tensão.

Conclusões

Através de simulações abrangentes rNEMD e EMD, quantificamos dependência de tamanho, temperatura e tensão de CT em OPG-L e OPG-Z. Suas condutividades são marcadamente menores que grafeno devido a velocidades de fonons e caminhos livres médios reduzidos, mas permanecem ajustáveis via engenharia de tensão. Estes insights estabelecem OPG-L e OPG-Z como candidatos promissores para gestão térmica personalizada em dispositivos micro e nanoescala.

Abreviações

558

Anéis de cinco-cinco-oito membros

OPG-L

Grafeno-linha octágono e pentágono

OPG-Z

Grafeno-ziguezague octágono e pentágono

rNEMD

Dinâmica molecular não-equilíbrio reversa

TC

Condutividade térmica

VDOS

Densidade de estados vibracionais