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

Resumo Técnico: Impacto da Desordem Eletrostática Dinâmica na Mobilidade de Lacunas em Rubreno

Definição do Problema
Cristais moleculares orgânicos de alta mobilidade, como o rubreno, são críticos para a eletrônica orgânica, embora alcançar uma descrição quantitativamente preditiva de seu transporte de carga permaneça um desafio. O transporte de carga nesses materiais ocorre em um regime intermediário entre o salto de pequeno-pólaron (small-polaron hopping) e o transporte de banda, caracterizado por uma deslocalização transitória onde a função de onda do portador de carga é modulada por flutuações térmicas. Embora a dinâmica molecular não adiabática de mistura quântico-clássica (especificamente o Fragment Orbital-based Surface Hopping, ou FOB-SH) ofereça uma rota promissora para simular esse regime sem assumir um mecanismo de transporte específico, simulações anteriores de grande escala de cristais apolares como o rubreno tipicamente negligenciaram a desordem eletrostática dinâmica. Essa omissão foi impulsionada pelo custo computacional de avaliar interações coulombianas de longo alcance (tipicamente exigindo a soma de Ewald) e pela suposição de que tais interações seriam desprezíveis em sistemas apolares. No entanto, negligenciar essas interações remove as contribuições eletrostáticas para a energia de reorganização e suprime as flutuações de energia de sítio, podendo levar a previsões imprecisas de deslocalização de portadores e mobilidade.

Metodologia
Este estudo integra a desordem eletrostática dinâmica em simulações de FOB-SH para o transporte de lacunas em temperatura ambiente no rubreno. Os autores empregam o método de soma eletrostática em espaço real Damped Shifted-Force (DSF), combinado com um esquema eficiente de adição-subtração, para superar os gargalos computacionais associados à avaliação de contribuições eletrostáticas para múltiplos estados localizados em carga.

Principais componentes metodológicos incluem:

• Tratamento Eletrostático: O método DSF aproxima as interações de Coulomb de longo alcance usando uma interação de espaço real amortecida que é deslocada para zero em um corte finito ($R_{cut} = 16$ Å). Um esquema de adição-subtração decompõe a energia eletrostática de um estado localizado em carga na energia do sistema neutro mais termos de correção, reduzindo a escala computacional de $O(N_{mol})$ para quase $O(1)$ em relação a um único cálculo.
• Efeitos de Polarização: Para contabilizar a polarização eletrônica sem o custo de campos de força polarizáveis explícitos, um esquema implícito baseado em escala de carga (fator $\gamma = 0,80$) foi aplicado às cargas parciais atômicas.
• Protocolo de Simulação: Simulações de FOB-SH foram realizadas a 300 K usando uma supercela de $50 \times 13 \times 1$ células unitárias (1300 moléculas). As simulações propagaram a função de onda da lacuna quanticamente enquanto evoluíam os núcleos classicamente. Os resultados foram comparados contra um modelo de "carga zero" (onde as cargas parciais foram definidas como zero) e validados contra cálculos de Smooth Particle-Mesh Ewald (SPME) para a energia de reorganização.
• Análise: O estudo quantificou a energia de reorganização ($\lambda$), a desordem de energia de sítio ($\sigma_E$), as flutuações de acoplamento eletrônico, a Razão de Participação Inversa (IPR) para medir a deslocalização da função de onda e a mobilidade do portador de carga ($\mu$) derivada do deslocamento quadrático médio (MSD).

Resultos Principais

1. Energia de Reorganização: A inclusão das interações eletrostáticas (via DSF ou SPME) aumenta a energia de reorganização para a transferência de lacunas em aproximadamente 20–40 meV (de uma linha de base de 152 meV para ~175–191 meV) em comparação ao modelo de carga zero. A análise espectral revela que este aumento surge do aumento do acoplamento através de ambos os movimentos de rede de baixa frequência (<200 cm$^{-1}$) e modos intramoleculares de alta frequência. O método DSF captura a maior parte desta contribuição eletrostática, mostrando boa concordância com os parâmetros de referência SPME.
2. Desordem de Energia de Sítio: A inclusão da eletrostática aumenta significativamente as flutuações de raiz quadrada média das energias de sítio ($\sigma_E$), subindo de 63,5 meV (carga zero) para 69,2 meV (DSF). Este aumento é quantitativamente consistente com o aumento da energia de reorganização, confirmando que a eletrostática modifica primariamente a desordem energética diagonal, deixando os acoplamentos eletrônicos off-diagonal amplamente inalterados.
3. Deslocalização da Função de Onda: A desordem de energia de sítio aumentada reduz a extensão espacial da função de onda da lacuna. A Razão de Participação Inversa (IPR) média diminui de 13 (carga zero) para 9 (DSF), indicando um portador de carga mais localizado.
4. Mobilidade de Carga: A redução na deslocalização leva a uma diminuição significativa na mobilidade de lacunas prevista. Ao longo da direção de alta mobilidade $a$, a mobilidade cai de 34,8 cm$^2$V$^{-1}$s$^{-1}$ (carga zero) para 21,1 cm$^2$V$^{-1}$s$^{-1}$ (DSF). Este valor está em estreita concordância com a referência experimental de 20,2 cm$^2$V$^{-1}$s$^{-1}$. Além disso, a inclusão da eletrostática melhora a razão de anisotropia prevista ($\mu_a/\mu_b$), aproximando-a das observações experimentais.

Significância e Alegações
O artigo demonstra que a desordem eletrostática dinâmica é um componente essencial para simulações quantitativamente preditivas de transporte de carga em semicondutores orgânicos moleculares, mesmo aqueles compostos por moléculas apolares como o rubreno. Ao incluir explicitamente essas interações via esquema computacionalmente eficiente DSF, os autores mostram que modelos anteriores que negligenciaram a eletrostática superestimaram a deslocalização e a mobilidade dos portadores. O estudo estabelece que as flutuações de energia de sítio eletrostática são o principal mecanismo responsável pela deslocalização dinâmica observada em cristais orgânicos de alta mobilidade, e sua inclusão é necessária para unir a lacuna entre as simulações teóricas e a realidade experimental. O trabalho valida o método DSF como uma alternativa viável e escalável à soma de Ewald para grandes escalas de dinâmica não adiabática em sistemas periódicos.