Tutorial de Simulação Eletro-Térmica: Autoaquecimento no OghmaNano (Parte A)

1. Introdução

O autoaquecimento é uma característica definidora de dispositivos eletrônicos e optoeletrônicos em operação. Sempre que a corrente flui, a potência elétrica é dissipada e convertida em calor, elevando a temperatura local do dispositivo. Esse aumento de temperatura retroalimenta diretamente transporte, recombinação e aprisionamento por meio da forte dependência da temperatura dos parâmetros do material. A simulação eletro-térmica é o tratamento autoconsistente desse problema elétrico–térmico acoplado.

Na prática, o autoaquecimento surge de múltiplos mecanismos físicos que atuam simultaneamente. Aquecimento por transporte de portadores converte energia elétrica em calor à medida que portadores de carga se movem através de energias de borda de banda espacialmente variáveis, abrangendo tanto aquecimento Joule em regiões homogêneas quanto aquecimento ou resfriamento Peltier em interfaces. Recombinação de portadores transfere energia eletrônica diretamente para a rede cristalina, enquanto calor adicional é gerado por perdas parasíticas em série e shunt. Capturar com precisão o comportamento do dispositivo sob polarização requer portanto que todos os termos de geração de calor sejam avaliados autoconsistentemente a partir da solução elétrica e retroalimentados por meio de um modelo de difusão térmica.

Uma simulação eletro-térmica consiste portanto em duas resoluções acopladas, refletindo o fato de que os problemas elétrico e térmico normalmente operam em escalas físicas de comprimento muito diferentes:

• uma resolução de transporte elétrico totalmente acoplado (drift–diffusion, Poisson, recombinação e armadilhas), geralmente confinada à região eletricamente ativa do dispositivo, e
• uma resolução de difusão térmica para o campo de temperatura da rede cristalina, que pode se estender muito além da região ativa para contatos, substratos e dissipadores de calor.

A estratégia de acoplamento implementada no OghmaNano é ilustrada em ??. Para uma determinada tensão aplicada, as equações elétricas são resolvidas usando o campo de temperatura atual. Os termos de geração de calor são então avaliados e passados ao solucionador térmico, que atualiza a temperatura da rede. Essa iteração externa continua até que os resíduos elétrico e térmico satisfaçam ambos os critérios de convergência.

Neste tutorial, você executará uma simulação eletro-térmica totalmente acoplada, identificará os mecanismos ativos de aquecimento e examinará como a temperatura da rede evolui com a tensão aplicada. Isso estabelece o fluxo de trabalho eletro-térmico que será estendido na Parte B para malhas térmicas detalhadas, condições de contorno e análise espacial da distribuição de calor.

2. Abrir o exemplo térmico

Na janela principal do OghmaNano, clique em Nova simulação. Na biblioteca de dispositivos, dê duplo clique em Simulação térmica (o pequeno ícone de vela no canto inferior esquerdo), mostrado em ??. Isso carrega um projeto de exemplo eletro-térmico completo com o modelo térmico e o acoplamento do solucionador já configurados. Após abrir o exemplo, a interface principal de simulação é exibida (??). O dispositivo é uma pilha de diodos multicamada com injeção, transporte e recombinação, escolhida porque produz densidades de corrente realistas e autoaquecimento mensurável sob polarização. Embora os materiais aqui correspondam a uma arquitetura orgânica, o objetivo do exemplo é demonstrar o fluxo de trabalho eletro-térmico em si, e não física específica de materiais.

Antes de executar a simulação, clique em Parâmetros elétricos na janela principal para abrir o editor de parâmetros elétricos (??). Essa visualização lista os parâmetros de transporte, recombinação e armadilhas usados pelo solucionador elétrico, que juntos definem os perfis de densidade de corrente e recombinação que geram calor durante a simulação eletro-térmica. Nesta etapa, nenhum parâmetro precisa ser modificado.

Neste exemplo, os parâmetros elétricos foram obtidos ajustando o modelo eletro-térmico acoplado a medições experimentais relatadas em In Situ Visualization and Quantification of Electrical Self-Heating in Conjugated Polymer Diodes Using Raman Spectroscopy (S. Maity, C. Ramanan, F. Ariese, R. C. I. MacKenzie, E. von Hauff, Adv. Electron. Mater., 2021, 2101208; https://doi.org/10.1002/aelm.202101208 ). Como resultado, vários parâmetros aparecem com múltiplas casas decimais; isso reflete o ponto final numérico do procedimento de ajuste, e não uma implicação de precisão absoluta dos parâmetros. Tal ajuste não é necessário para modelagem eletro-térmica geral: para a maioria dos dispositivos, parâmetros padrão da literatura ou valores nominais razoáveis são suficientes para obter insight fisicamente significativo sobre mecanismos de aquecimento, aumento de temperatura e acoplamento eletro-térmico, sendo o ajuste de parâmetros necessário apenas quando é exigida concordância quantitativa com um experimento específico.

3. Examinando a faixa térmica (fontes de calor habilitadas)

Antes de executar a simulação, navegue até a faixa Térmica na janela principal. O ícone de vela fornece um indicador simples de estado: uma vela acesa indica que o modelo térmico está habilitado, enquanto uma vela apagada indica que a simulação será executada a temperatura fixa. Os dois estados são mostrados em ?? e ??. Quando o modelo térmico está habilitado, múltiplos mecanismos de geração de calor podem ser ativados. Cada termo habilitado contribui para a fonte volumétrica total de calor usada pelo solucionador térmico. Neste exemplo, os seguintes mecanismos estão incluídos:

• Aquecimento por drift (rotulado Joule heating na faixa), incluindo efeitos termoelétricos
• Aquecimento parasítico, decorrente de perdas resistivas em série e shunt
• Aquecimento por recombinação, devido à transferência de energia dos portadores para a rede
• Aquecimento óptico, disponível quando a absorção óptica é convertida em calor

A faixa também fornece acesso a Parâmetros térmicos, Malha térmica, e Condições de contorno. Esses controles são importantes porque o problema térmico normalmente é definido em um domínio físico muito maior do que o problema elétrico. Seu papel é discutido em detalhe na Seção 6 e na Parte B.

3.1 Aquecimento por drift (contribuições Joule e Peltier)

A fonte de calor relacionada ao drift usada no OghmaNano é escrita como:

\[ Q_{\mathrm{drift}} = \mathbf{J}_n \cdot \nabla E_C + \mathbf{J}_p \cdot \nabla E_V \]

Essa forma captura tanto o aquecimento Joule resistivo convencional quanto efeitos termoelétricos (semelhantes ao efeito Peltier) associados a variações espaciais na energia dos portadores. Em regiões dominadas por transporte resistivo, o termo é positivo e corresponde à geração de calor. Próximo a interfaces ou em regiões de forte curvatura de banda, contribuições locais negativas podem ocorrer, correspondendo a resfriamento ou aquecimento da rede induzido por portadores. A presença de regiões positivas e negativas em gráficos de aquecimento por drift é portanto esperada e reflete o balanço energético local, e não ruído numérico.

3.2 Aquecimento parasítico (perdas em série e shunt)

Além de processos intrínsecos de transporte e recombinação, dispositivos práticos dissipam energia em elementos resistivos de série e shunt. No modelo eletro-térmico, essas perdas são representadas diretamente como um termo volumétrico de geração de calor:

\[ Q_{\mathrm{parasitic}} = \frac{I^2 R_s + \dfrac{V^2}{R_{sh}}}{V_{\mathrm{dev}}} \]

Aqui, \(V_{\mathrm{dev}}\) é o volume sobre o qual a dissipação parasítica é distribuída. Como a localização microscópica das perdas parasíticas geralmente não é resolvida pelo modelo elétrico, essa contribuição é tratada como uma fonte de calor espacialmente distribuída, garantindo conservação global de energia sem introduzir suposições sobre a localização de pontos quentes não resolvidos.

3.3 Aquecimento por recombinação

A recombinação de portadores transfere energia eletrônica para a rede cristalina. Neste modelo, o termo correspondente de geração de calor é escrito como:

\[ Q_{\mathrm{rec}} = \left(\langle w_n \rangle + E_g + \langle w_p \rangle \right)\,R \]

Aqui, \(R\) é a taxa de recombinação, enquanto \(\langle w_n \rangle\) e \(\langle w_p \rangle\) representam contribuições médias de energia dos portadores relativas às bordas de banda. O gap de banda \(E_g\) define a escala dominante de energia. Como a recombinação é frequentemente espacialmente localizada, o aquecimento por recombinação geralmente exibe uma distribuição espacial diferente daquela associada ao aquecimento por drift, tornando útil analisar essas contribuições separadamente.

Outros controles térmicos disponíveis na faixa são examinados em detalhe na Parte B. Para a Parte A, o objetivo é simplesmente identificar quais fontes de calor estão habilitadas e compreender sua interpretação física antes de executar a simulação.

4. Executar a simulação eletro-térmica

Execute a simulação pressionando o botão azul Executar simulação na janela principal (ou pressione F9). Após a execução ser concluída, todos os resultados são gravados no diretório de saída mostrado em ??.

A execução eletro-térmica produz as saídas elétricas padrão, como a curva corrente–tensão (JV), juntamente com saídas térmicas adicionais. Em particular, juntamente com a característica JV (??), é gerada uma saída de temperatura da rede dependente da tensão (??), representando a temperatura autoconsistente do dispositivo sob polarização.

Em uma simulação eletro-térmica, a curva JV não é mais avaliada a uma temperatura fixa da rede cristalina. Em vez disso, a temperatura é determinada autoconsistentemente em cada ponto de polarização por meio da resolução acoplada elétrico–térmica. Uma vez que um dispositivo entra em condução, a potência dissipada \(P \sim IV\) pode aumentar rapidamente, e o aumento de temperatura resultante pode modificar de forma mensurável parâmetros de transporte e recombinação, levando a mudanças observáveis na inclinação e curvatura da curva JV.