In the past years, organic electrochemical transistors (OECTs) have emerged as potential transducers inapplications that require the conversion of ion fluxes to electronic current. For the rational optimization and understanding of the fundamentals of OECTs and OECT-based applications, however, it is essential to have theoretical models capable to predict experimental data. Within drift-diffusion models,ion flux from the electrolyte into the organic semiconducting layer is considered to take place due to the action of an electrical field, but also because of diffusion processes generated by the concentration gradients. The governing equation of the drift-diffusion model is the Nernst-Planck equation.Thus, in this project, a numerical approach is followed in order to solve the Nernst-Planck equation in one dimension,and model the ion migration from the electrolyte to the semiconductor. To evaluate the accuracy of the implementation, standard boundary conditions used in the literature to solve analytically the drift-diffusion equations were considered.In doing so,the numerical results were in good agreement with the analytical solutions,achieving maximum errors in the order of 1%. Aiming to a better representation of OECTs, closed boundary conditions are considered. Here, the temporal evolution of the concentration profiles showed a convergence to an exponential steady state distribution, which is in good agreement with the result expected theoretically. A further situation investigated was the consideration of a non-uniform electric field acting on the system, assumed to be finite in the electrolyte regionandzerointhe semiconductor.This Consideration impacts principally the temporal evolution of the concentration in each region. In order to consider the distinct compositions in electrolyte and semiconductors, different values of diffusion coefficients were introduced for each region. This extension has visible impacts in the time that the system needs to achieve the steady state. Moreover, the introduction of the chemical potential gradient as the driving force of diffusion leaded to significant variations in the results obtained with the model. Here,the so-called uphill diffusion reported in the literature was observed. With the numerical approach,it was possible to consider different types of pulsed gate voltages, which allowed to simulate the charge and discharge processes of OECTs. For all cases, oscillatory curves similar to experimental measurements were obtained. Therefore, the numerical approach allowed to go beyond the analytical description, and develop an extensive investigation of the impact that different considerations have in the results.

Nos últimos anos, os transistores eletroquímicos orgânicos (OECTs) surgiram como transdutores em aplicações que requerem a conversão de fluxos iônicos em corrente eletrônica. Para a otimização e compreensão dos fundamentos das OECTs e suas aplicações, é essencial ter modelos teóricos capazes de prever dados experimentais. Nos modelos de drift-diffusion, o fluxo de íons do eletrólito para a camada semicondutora orgânica é considerado como tendo lugar devido à ação de um campo elétrico, mas também devido aos processos de difusão gerados pelos gradientes de concentração. A equação governante do modelo de drift-diffusion é a equação de Nernst-Planck. Assim, neste projeto, uma abordagem numérica é seguida para resolver a equação de Nernst-Planck em uma dimensão e modelar a migração iônica do eletrólito para o semicondutor. Para avaliar a precisão da implementação, foram consideradas as condições de contorno padrão utilizadas na literatura para resolver analiticamente as equações de drift-diffusion. Ao fazê-lo, os resultados numéricos obtidos mostraram boa concordância com as soluções analíticas, alcançando erros na ordem de 1% .Visando uma melhor representação das OECTs, são consideradas as condições de contorno fechado. Aqui, a evolução temporal dos perfis de concentração mostraram uma convergência para uma distribuição exponencial no estado estacionário, que está em boa concordância com o resultado esperado teoricamente. Adicionalmente, foi considerado um campo elétrico não uniforme atuando sobre o sistema, assumido como finito na região do eletrólito e zero no semicondutor. Esta consideração afeta principalmente a evolução temporal da concentração em cada região. Como objetivo de considerar as composições distintas em eletrólitos e semicondutores, foram introduzidos diferentes valores de coeficientes de difusão para cada região. Esta extensão tem impactos visíveis no tempo que o sistema precisa para alcançar o estado estacionário. Além disso, a introdução do gradiente do potencial químico como a força responsável da difusão dos íons levou a variações significativas nos resultados obtidos como modelo. Aqui, a chamada difusão scendente relatada na literatura foi observada. Com a abordagem numérica, foi possível considerar diferentes tipos de tensões pulsadas no eletrodo da porta, o que permitiu simular os processos de carga e descarga dos OECTs. Para todos os casos, foram obtidas curvas oscilatórias semelhantes às medidas experimentais. Portanto, a abordagem numérica permitiu ir além da descrição analítica, e desenvolver uma extensa investigação do impacto que diferentes considerações têm nos resultados.