A necessidade de estados de armadilha em modelos de dispositivos orgânicos

Esta seção explica por que estados de armadilha devem ser considerados ao simular materiais desordenados, como misturas polímero:fulereno, sistemas de pequenas moléculas ou semicondutores amorfos. Sem um tratamento explícito de aprisionamento e liberação, qualquer modelo de dispositivo deixará de capturar a física do transporte de carga e da recombinação. O uso do formalismo completo de recombinação e aprisionamento de Shockley–Read–Hall (SRH) é, portanto, essencial para obter resultados fisicamente significativos.

1. A estrutura física e energética de materiais desordenados

Semicondutores inorgânicos tradicionais, como Si cristalino ou GaAs, são ao mesmo tempo altamente ordenados e extremamente puros — frequentemente alcançando pureza “nove noves” (99,9999999%). Semicondutores orgânicos, em contraste, raramente excedem 99,9% de pureza, tornando-os cerca de um milhão de vezes mais ricos em defeitos do que seus equivalentes inorgânicos. Estruturalmente, a diferença é igualmente marcante: semicondutores inorgânicos formam redes cristalinas regulares, como bolinhas de gude cuidadosamente organizadas em um tabuleiro de paciência (??). O próprio silício adota a estrutura cúbica do diamante, uma rede quase perfeita (??). Materiais orgânicos, em contraste, são sistemas moleculares “moles”, com polímeros emaranhados que lembram um prato de espaguete à bolonhesa — o espaguete representando as cadeias poliméricas e o molho, as pequenas moléculas (??). Simulações de misturas orgânicas confirmam essa imagem, mostrando empacotamento polimérico semelhante a espaguete com derivados de fulereno intercalados por toda a estrutura (??).

Essas diferenças estruturais dão origem a paisagens energéticas muito diferentes. Em semicondutores cristalinos, elétrons e lacunas se movem livremente em bandas de condução e valência bem definidas, experimentando apenas resistência modesta sob um campo aplicado. Esse transporte do tipo banda é ilustrado em ??. Em semicondutores orgânicos desordenados, porém, impurezas e desordem estrutural introduzem uma distribuição densa de estados de armadilha localizados dentro do gap de banda. Em vez de se propagarem livremente por estados estendidos, os portadores devem saltar termicamente entre armadilhas. Esse transporte dominado por armadilhas é mostrado esquematicamente em ??.

A implicação é clara: enquanto estados de armadilha podem frequentemente ser negligenciados em semicondutores ordenados, eles dominam a física de sistemas desordenados. Qualquer modelo realista de dispositivo deve, portanto, incluir uma descrição detalhada das distribuições de armadilhas e da cinética de Shockley–Read–Hall. OghmaNano faz exatamente isso, tornando possível simular sistemas desordenados como células solares orgânicas, OFETs e perovskitas com precisão física.

Para entender por que isso importa, passamos agora para a relação entre densidade de portadores e nível de Fermi, examinando alguns exemplos concretos.

2. Por que estados de armadilha são importantes para modelagem de dispositivos (sem matemática)

Uma característica central de semicondutores orgânicos e outros semicondutores desordenados é que a densidade de portadores é uma função forte tanto da tensão aplicada quanto da intensidade de iluminação. À medida que a polarização ou a intensidade luminosa aumenta, mais carga é injetada ou fotogerada no dispositivo. Como esses materiais possuem um grande número de estados de armadilha dentro do gap de banda, os portadores primeiro preenchem essas armadilhas antes de ocupar estados estendidos. Esse preenchimento de armadilhas significa que mesmo pequenas mudanças na tensão podem produzir grandes mudanças na densidade de portadores livres.

Isso importa porque a recombinação no dispositivo depende diretamente das densidades de portadores. A forma geral da taxa de recombinação é

\[ R = k_r \, n(V)\, p(V), \]

onde \(k_r\) é a constante de recombinação, e \(n(V)\) e \(p(V)\) são as densidades de elétrons e lacunas dependentes da tensão. Se a forma funcional de \(n(V)\) (e \(p(V)\)) estiver errada porque os estados de armadilha foram negligenciados, então a taxa de recombinação também estará errada. Isso leva diretamente a previsões incorretas da tensão de circuito aberto (\(V_{OC}\)) e de outras características importantes do dispositivo.

A mobilidade é afetada de maneira semelhante. Em sistemas desordenados, a mobilidade efetiva do portador depende do balanço entre portadores livres e portadores aprisionados. Uma expressão simples é

\[ \mu_e(n) = \frac{\mu_e^0 \, n_{\text{free}}}{n_{\text{free}} + n_{\text{trap}}}, \]

onde \(\mu_e^0\) é a mobilidade intrínseca do elétron, \(n_{\text{free}}\) é a densidade de portadores móveis e \(n_{\text{trap}}\) é a densidade de portadores aprisionados. Se a relação densidade–tensão estiver errada, então a dependência prevista mobilidade–tensão também estará errada. Juntos, esses erros significam que a curva J–V simulada não corresponderá ao experimento, mesmo que os parâmetros de recombinação ou mobilidade sejam razoáveis em outros aspectos.

3. Por que estados de armadilha são importantes para modelagem de dispositivos (com matemática)

Para descrever corretamente as densidades de portadores, é necessário levar em conta a densidade de estados (DoS) subjacente. A Figura ?? esquematiza a DoS correspondente às estruturas de banda ordenadas e desordenadas mostradas anteriormente em ?? e ??. Em um semicondutor ordenado, a DoS possui uma borda de banda abrupta (banda parabólica), e a população de Fermi–Dirac fica acima da borda da banda de condução. Em um semicondutor desordenado, porém, a DoS exibe uma cauda de estados de armadilha localizados que se estende profundamente dentro do gap de banda (comumente modelada com caudas exponenciais ou gaussianas). A consequência é que os portadores ocupam distribuições fundamentalmente diferentes nos dois casos.

Formalmente, a densidade total de elétrons é dada pela integral da DoS ponderada pela ocupação de Fermi–Dirac:

\[ n(E_f,T) = \int_{E_{\min}}^{\infty} \rho(E)\, f(E,E_f,T)\, dE, \]

onde \(E_f\) é o quasi-nível de Fermi, \(\rho(E)\) é a DoS e \(f(E,E_f,T)\) é a distribuição de Fermi–Dirac. Em materiais ordenados, \(\rho(E)\) é abrupta, de modo que apenas estados acima da borda da banda de condução contribuem. Em materiais desordenados, \(\rho(E)\) inclui extensas caudas de armadilha, de modo que os portadores podem ser armazenados próximos ao nível de Fermi, produzindo densidades de carga uma a duas ordens de grandeza maiores do que em um cristal ordenado sob a mesma polarização. Isso é observado diretamente em experimentos de extração de carga.

A conclusão é direta: obter corretamente a dependência da densidade de portadores não é opcional. Sem um modelo de estados de armadilha, tanto a recombinação quanto a mobilidade estarão erradas em função da tensão, e curvas J–V realistas não poderão ser reproduzidas. OghmaNano inclui esses estados de armadilha explicitamente, permitindo modelagem precisa de dispositivos desordenados, como misturas PM6:Y6 e P3HT:PCBM, bem como semicondutores mais ordenados quando as armadilhas são desativadas.

4. Por que você não deve usar recombinação de Langevin em modelos de dispositivos

A taxa clássica de recombinação de Langevin é definida como

\[ R_{\text{free}} = q \, k_r \, \frac{\mu_e + \mu_h}{2 \, \epsilon_0 \epsilon_r} \, n p , \]

onde \(R_{\text{free}}\) é a taxa de recombinação, \(k_r\) é o fator empírico de redução de Langevin, \(\mu_e\) e \(\mu_h\) são as mobilidades de elétrons e lacunas, \(n\) e \(p\) são as densidades de portadores, e \(\epsilon_0 \epsilon_r\) é a permissividade dielétrica. À primeira vista isso parece fisicamente razoável: assume-se que a recombinação ocorre sempre que um elétron e uma lacuna se aproximam o suficiente sob movimento browniano para sentir o campo coulombiano um do outro. Essa imagem é apropriada para portadores perfeitamente livres em líquidos simples ou condutores iônicos. Mas em fotovoltaicos orgânicos (OPV) e outros semicondutores desordenados, as hipóteses por trás do modelo de Langevin deixam de valer.

Por que ele falha? Estudos experimentais — especialmente do início da década de 2010 — mostraram rapidamente que a recombinação de Langevin não podia reproduzir de forma autoconsistente tanto as curvas J–V no escuro quanto sob iluminação. O modelo superestimava sistematicamente as taxas de recombinação, frequentemente por ordens de magnitude. Para tornar os ajustes possíveis, pesquisadores introduziram um “fator de redução de Langevin” \(k_r\), por vezes tão pequeno quanto 10⁻³. Embora conveniente, esse ajuste foi na verdade uma admissão de que o próprio mecanismo não era válido nesses sistemas.

Os problemas ficam claros quando analisamos mais de perto a equação:

• Dependência errada da densidade de portadores. A taxa escala estritamente como \(np\), ou seja, segunda ordem. Em experimentos, a recombinação em OPVs frequentemente segue leis de potência não inteiras, como \((np)^{1.5}\), refletindo dinâmica limitada por armadilhas e a distribuição de estados localizados. A forma de Langevin simplesmente não consegue capturar isso.
• Tratamento incorreto da mobilidade. A equação assume mobilidades fixas \(\mu_e\) e \(\mu_h\). Na realidade, a mobilidade em semicondutores desordenados depende fortemente da densidade de portadores (veja a discussão acima). Incorporar a dependência errada da mobilidade distorce ainda mais a taxa de recombinação prevista.
• Ignora armadilhas e desordem. O mecanismo de Langevin trata todos os portadores como livres. Em OPVs e materiais semelhantes, a maioria dos portadores passa tempo aprisionada em estados localizados e só intermitentemente contribui para a condução. A recombinação, portanto, depende da paisagem de armadilhas, não apenas do movimento de portadores livres.

Em conjunto, essas questões significam que a recombinação de Langevin é, na melhor das hipóteses, uma aproximação grosseira e, na pior, uma aproximação enganosa. Mesmo com um fator de redução ajustado \(k_r\), ela não consegue capturar a física correta da recombinação assistida por armadilhas e a dependência da mobilidade com a tensão. Usar recombinação de Langevin em modelos de dispositivos é, portanto, como forçar um pino quadrado em um buraco redondo: pode lhe dar um número, mas não lhe dará resultados fisicamente significativos.

5. Como se pode fazer a recombinação de Langevin “funcionar” em modelos de dispositivos

Os problemas centrais da recombinação clássica de Langevin são sua dependência incorreta da densidade de portadores e a necessidade de um fator de redução arbitrário. Uma forma pela qual pesquisadores tentaram fazer a recombinação de Langevin “funcionar” foi introduzindo uma dependência da densidade de portadores nas próprias mobilidades, como em:

\[ R_{\text{free}} = q \, k_r \, \frac{\alpha \mu_e(n) + \beta \mu_h(n)}{2 \, \epsilon_0 \epsilon_r} \, n_{\text{tot}} p_{\text{tot}} , \]

Aqui, define-se uma mobility edge: portadores acima da mobility edge contribuem para a condução, enquanto aqueles abaixo são considerados aprisionados. As mobilidades médias podem então ser expressas como

\[ \mu_e(n) = \frac{\mu_e^0 \, n_{\text{free}}}{n_{\text{free}} + n_{\text{trap}}}, \qquad \mu_h(p) = \frac{\mu_h^0 \, p_{\text{free}}}{p_{\text{free}} + p_{\text{trap}}}. \]

Se a densidade de portadores livres for muito menor que a densidade de portadores aprisionados, isso leva a uma taxa efetiva de recombinação

\[ R(n,p) = q \, k_r \, \frac{\alpha \mu_e^0 \, n_{\text{free}} p_{\text{trap}} + \beta \mu_h^0 \, p_{\text{free}} n_{\text{trap}}} {2 \, \epsilon_0 \epsilon_r}. \]

Dessa forma, a recombinação de Langevin é efetivamente reinterpretada em termos de interações entre portadores livres e aprisionados (\(n_{\text{free}}p_{\text{trap}}\) e \(p_{\text{free}}n_{\text{trap}}\)). Isso é, em essência, equivalente à imagem de recombinação Shockley–Read–Hall (SRH): portadores livres recombinando com portadores aprisionados.

Embora essa abordagem funcione razoavelmente bem em regime estacionário, ela se baseia em uma hipótese forte: que todos os portadores em uma dada posição compartilham um único quasi-nível de Fermi, isto é, estão em equilíbrio local com velocidade de termalização infinita. Isso pode ser plausível sob condições de regime estacionário, quando os portadores têm tempo para se equilibrar, mas essa hipótese falha no domínio do tempo. Na realidade, a distribuição densa de estados de armadilha em semicondutores orgânicos torna improvável que os portadores possam agir como um único gás equilibrado. Em contraste, o formalismo SRH evita essa suposição e, portanto, é uma descrição fisicamente mais sólida da recombinação e do aprisionamento em materiais desordenados.