Impact of dynamic electrostatic disorder on hole mobility in rubrene: a nonadiabatic molecular dynamics investigation

Resumen Técnico: Impacto del Desorden Electroestático Dinámico en la Movilidad de Huecos en el Rubreno

Planteamiento del Problema
Los cristales moleculares orgánicos de alta movilidad, como el rubreno, son fundamentales para la electrónica orgánica; sin embargo, lograr una descripción cuantitativamente predictiva de su transporte de carga sigue siendo un desafío. El transporte de carga en estos materiales ocurre en un régimen intermedio entre el salto de pequeño polarón (small-polaron hopping) y el transporte de banda, caracterizado por una deslocalización transitoria donde la función de onda del portador de carga es modulada por fluctuaciones térmicas. Si bien la dinámica molecular no adiabática de base cuántica y clásica mixta (específicamente, el Surface Hopping basado en Orbitales de Fragmentos, o FOB-SH) ofrece una vía prometedora para simular este régimen sin asumir un mecanismo de transporte específico, las simulaciones previas a gran escala de cristales apolares como el rubreno han omitido típicamente el desorden electroestático dinámico. Esta omisión fue impulsada por el costo computacional de evaluar las interacciones coulombianas de largo alcance (que típicamente requieren la sumación de Ewald) y la suposición de que tales interacciones son despreciables en sistemas apolares. Sin embargo, omitir estas interacciones elimina las contribuciones electroestáticas a la energía de reorganización y suprime las fluctuaciones de la energía de sitio, lo que puede conducir a predicciones inexactas de la deslocalización y la movilidad de los portadores.

Metodología
Este estudio integra el desorden electroestático dinámico en las simulaciones de FOB-SH para el transporte de huecos a temperatura ambiente en el rubreno. Los autores emplean el método de sumación electroestática en el espacio real de Fuerza de Desplazamiento Amortiguado (DSF, por sus siglas en inglés), combinado con un esquema eficiente de adición-sustracción, para superar los cuellos de botella computacionales asociados con la evaluación de las contribuciones electroestáticas para múltiples estados localizados de carga.

Los componentes metodológicos clave incluyen:

• Tratamiento Electroestático: El método DSF aproxima las interacciones de Coulomb de largo alcance utilizando una interacción en el espacio real amortiguada y desplazada a cero en un corte finito ($R_{cut} = 16$ Å). Un esquema de adición-sustracción descompone la energía electroestática de un estado de carga localizada en la energía del sistema neutro más términos de corrección, reduciendo la escala computacional de $O(N_{mol})$ a casi $O(1)$ respecto a un único cálculo.
• Efectos de Polarización: Para contabilizar la polarización electrónica sin el costo de campos de fuerza polarizables explícitos, se aplicó un esquema implícito basado en el escalado de carga (factor $\gamma = 0.80$) a las cargas parciales atómicas.
• Protocolo de Simulación: Las simulaciones de FOB-SH se realizaron a 300 K utilizando una supercelda de $50 \times 13 \times 1$ celdas unitarias (1300 moléculas). Las simulaciones propagaron la función de onda del hueco de forma cuántica mientras evolucionaban los núcleos de forma clásica. Los resultados se compararon contra un modelo de referencia de "carga cero" (donde las cargas parciales se establecieron en cero) y se validaron mediante cálculos de Ewald de Partícula de Malla Suave (SPME) para la energía de reorganización.
• Análisis: El estudio cuantificó la energía de reorganización ($\lambda$), el desorden de la energía de sitio ($\sigma_E$), las fluctuaciones del acoplamiento electrónico, la Razón de Participación Inversa (IPR) para medir la deslocalización de la función de onda y la movilidad del portador de carga ($\mu$) derivada del desplazamiento cuadrático medio (MSD).

Resultos Clave

1. Energía de Reorganización: La inclusión de las interacciones electroestáticas (vía DSF o SPME) aumenta la energía de reorganización para la transferencia de huecos en aproximadamente 20–40 meV (de un valor base de 152 meV a ~175–191 meV) en comparación con el modelo de carga cero. El análisis espectral revela que este aumento surge de un acoplamiento mejorado a través de tanto los movimientos de la red de baja frecuencia (<200 cm$^{-1}$) como de los modos intramolecular de mayor frecuencia. El método DSF captura la mayor parte de esta contribución electroestática, mostrando una buena concordancia con los referentes de SPME.
2. Desorden de la Energía de Sitio: La inclusión de la electroestática aumenta significativamente las fluctuaciones de la raíz cuadrada media de las energías de sitio ($\sigma_E$), pasando de 63.5 meV (carga cero) a 69.2 meV (DSF). Este aumento es cuantitativamente consistente con el incremento en la energía de reorganización, confirmando que la electroestática modifica principalmente el desorden energético diagonal, dejando los acoplamientos electrónicos fuera de la diagonal prácticamente inalterados.
3. Deslocalización de la Función de Onda: El aumento del desorden de la energía de sitio reduce la extensión espacial de la función de onda del hueco. El promedio de la Razón de Participación Inversa (IPR) disminuye de 13 (carga cero) a 9 (DSF), indicando un portador de carga más localizado.
4. Movilidad de Carga: La reducción en la deslocalización conduce a una disminución significativa en la movilidad de huecos predicha. A lo largo de la dirección de alta movilidad $a$, la movilidad cae de 34.8 cm$^2$V$^{-1}$s$^{-1}$ (carga cero) a 21.1 cm$^2$V$^{-1}$s$^{-1}$ (DSF). Este valor concuerda estrechamente con la referencia experimental de 20.2 cm$^2$V$^{-1}$s$^{-1}$. Además, la inclusión de la electroestática mejora la razón de anisotropía predicha ($\mu_a/\mu_b$), acercándola a las observaciones experimentales.

Significancia y Reivindicaciones
El artículo demuestra que el desorden electroestático dinámico es un componente esencial para las simulaciones cuantitativamente predictivas del transporte de carga en semiconductores orgánicos moleculares, incluso en aquellos compuestos por moléculas apolares como el rubreno. Al incluir explícitamente estas interacciones mediante el esquema computacionalmente eficiente de DSF, los autores muestran que los modelos previos que omitían la electrostática sobreestimaban la deslocalización y la movilidad de los portadores. El estudio establece que las fluctuaciones de la energía de sitio electroestática son el mecanismo principal responsable de la localización dinámica observada en los cristales orgánicos de alta movilidad, y su inclusión es necesaria para cerrar la brecha entre las simulaciones teóricas y la realidad experimental. El trabajo valida el método DSF como una alternativa viable y escalable a la sumación de Ewald para la dinámica no adiabática a gran escala en sistemas periódicos.