Lattice initialisation and finite-size effects of non-equilibrium molecular dynamics simulations for heat transfer across graphene-copper interfaces

Resumen Técnico: Inicialización de la Red y Efectos de Tamaño Finito en Simulaciones NEMD de Transferencia de Calor Grafeno–Cobre

Planteamiento del Problema
Las simulaciones de Dinámica Molecular de No Equilibrio (NEMD) se utilizan ampliamente para predecir la conductancia térmica interfacial (ITC) y la resistencia de Kapitza en compuestos de cobre–grafeno, los cuales son críticos para aplicaciones de enfriamiento electrónico. Sin embargo, la fiabilidad de estas predicciones se ve a menudo comprometida por las incertidumbres en las elecciones de la configuración computacional. Específicamente, la literatura carece de una comprensión sistemática de cómo las longitudes del dominio computacional finito y, más críticamente, la inicialización de las redes atómicas (incluyendo la deformación residual y las constantes de red) influyen en las propiedades de transporte térmico predichas. Si bien se sabe que la NEMD es sensible al tamaño del dominio debido a las limitaciones del camino libre medio de los fonones, el alcance con el que las estrategias de inicialización de la red —elegidas a menudo arbitrariamente para coincidir con los parámetros de red experimentales o para minimizar la deformación— afectan la dinámica interfacial permanece poco explorado.

Metodología
Los autores emplearon simulaciones directas de NEMD utilizando el paquete LAMMPS para modelar la transferencia de calor a través de una interfaz de grafeno de capa única embebida entre dominios de cobre. El estudio utilizó una combinación de campos de fuerza: el Método de Átomo Embebido (EAM) para el cobre, el potencial de Orden de Enlace Reactivo de Interacción Adaptativa (AIREBO) para el carbono, y un potencial de Lennard-Jones para las interacciones cobre–carbono.

Para aislar fuentes específicas de incertidumbre, el estudio restringió las variaciones a la longitud del dominio ($L_z$), la estrategia de inicialización de la red y los métodos de imposición de temperatura, manteniendo constantes el campo de fuerza, el termostato y el área de la sección transversal. Se compararon dos estrategias distintas de inicialización de la red:

• Caso I: Inicialización basada en parámetros de red derivados experimentalmente ($a_{Cu} = 3.6$ Å, $a_G = 2.46$ Å), lo que resulta en un sistema con deformación residual y una capa de grafeno plana.
• Caso II: Inicialización basada en un procedimiento de relajación de la red (equilibrio NPT) para minimizar la deformación por desajuste, lo que resulta en parámetros de red ligeramente diferentes ($a_{Cu} \approx 3.61$ Å, $a_G \approx 2.41$ Å) y una capa de grafeno visiblemente arrugada.

Se probaron dos esquemas de imposición de temperatura para evaluar los artefactos de frontera:

1. $\Delta T$ Constante: Diferencia de temperatura fija a través de todo el dominio.
2. $\delta T$ Constante: Diferencia de temperatura fija específicamente a través de la interfaz, con temperaturas de frontera extrapoladas.

El post-procesamiento consistió en calcular las temperaturas locales mediante el teorema de equipartición, ajustar perfiles de temperatura lineales para derivar el flujo de calor y la resistencia de Kapitza ($R_K$), y analizar la densidad de estados vibracionales (VDOS) para comprender los espectros de fonones y su solapamiento.

Contribuciones Clave y Resultados

• Sensibilidad a la Inicialización de la Red: El hallazgo más significativo es la extrema sensibilidad de la resistencia de Kapitza a las elecciones de inicialización de la red. A pesar de que los dos casos difieren solo en unos pocos puntos porcentuales en los parámetros de red y la densidad atómica, el Caso I (mayor deformación residual, grafeno plano) produjo una ITC de aproximadamente 557 MW/m²K, mientras que el Caso II (menor deformación residual, grafeno arrugado) produjo una ITC de aproximadamente 298 MW/m²K, una diferencia de casi un factor de dos.
• Mecanismo de Resistencia: El estudio desafía la visión convencional de que el solapamiento espectral por sí solo dicta la transferencia de calor interfacial. Aunque el Caso II exhibió un ligero mayor solapamiento de modos de fonones entre el cobre y el grafeno, mostró una resistencia térmica significativamente mayor. Los autores atribuyen esto a la formación de una "capa límite de amortiguamiento" en el Caso II de menor deformación. Esta capa se caracteriza por:
  • Desorden estructural y arrugamiento de la interfaz de grafeno.
  • Mayor desorden estructural local en los átomos de cobre adyacentes (distribución más amplia de la distancia de los vecinos más cercanos).
  • Ensanchamiento espectral y aumento de la dispersión en las regiones de cobre interfacial, lo que suprime la transferencia de calor.
• Efectos de Tamaño Finito en la Conductividad Bulk vs. ITC:
  • Conductancia Interfacial (ITC): Una vez contabilizados los efectos relacionados con la deformación (variaciones de densidad), la ITC no mostró una dependencia significativa de la longitud del dominio o del método específico de imposición de temperatura ($\Delta T$ vs. $\delta T$). Esto sugiere que la ITC es robusta ante artefactos de frontera y elecciones de termostato dentro del rango investigado.
  • Conductividad del Cobre Bulk: En contraste, la conductividad térmica de los dominios de cobre adyacentes a la interfaz exhibió una clara dependencia tanto de la longitud del dominio como de la temperatura. Los dominios más pequeños mostraron una conductividad reducida debido a la dispersión en las fronteras que limita el camino libre medio de los fonones. Además, la conductividad varió entre los lados caliente y frío debido a los mecanismos de dispersión dependientes de la temperatura (por ejemplo, la dispersión Umklapp), consistente con la teoría establecida de la conductividad térmica de red.
• Análisis de los Espectros de Fonones: El estudio demostró que los desplazamientos en los espectros de fonones bulk (debidos a la deformación) son insuficientes para explicar las diferencias en la ITC. En su lugar, la dinámica interfacial local, específicamente el ensanchamiento de los modos de fonones y el aumento del desorden en las inmediaciones de la interfaz, son los factores dominantes que gobiernan la resistencia térmica.

Significancia y Reivindicaciones
El artículo sostiene que, si bien las simulaciones NEMD pueden producir valores numéricos precisos para las propiedades térmicas, la fiabilidad de los resultados interfaciales depende fuertemente del protocolo de inicialización. Los autores argumentan que el "desplazamiento dependiente de la deformación" en las constantes de red es un motor primario de la variabilidad de la ITC, superando a menudo los efectos de tamaño finito.

El estudio concluye que el solapamiento espectral es una métrica insuficiente para predecir la transferencia de calor interfacial en este sistema; el desorden estructural y dinámico local en la interfaz juega un papel más crítico. Por consiguiente, los autores afirman que para determinar de manera fiable las propiedades de transporte térmico para compuestos de grafeno–cobre, los investigadores deben controlar y reportar cuidadosamente las estrategias de inicialización de la red y la deformación residual, en lugar de confiar únicamente en la convergencia del tamaño del dominio o en los argumentos de solapamiento espectral. Los hallazgos resaltan la necesidad de una mayor exploración de la sensibilidad de la NEMD a los parámetros de entrada variados continuamente para distinguir la física interfacial genuina de los artefactos de simulación.