Abstract:: Because charge carriers of many organic semiconductors (OSCs) exhibit fractional drift diffusion (Fr-DD) transport properties, the need to develop a Fr-DD model solver becomes more apparent. However, the current research on solving the governing equations of the Fr-DD model is practically nonexistent. In this paper, an iterative solver with high precision is developed to solve both the transient and steady-state Fr-DD model for organic semiconductor devices. The Fr-DD model is composed of two fractional-order carriers (i.e., electrons and holes) continuity equations coupled with Poisson’s equation. By treating the current density as constants within each pair of consecutive grid nodes, a linear Caputo’s fractional-order ordinary differential equation (FrODE) can be produced, and its analytic solution gives an approximation to the carrier concentration. The convergence of the solver is guaranteed by implementing a successive over-relaxation (SOR) mechanism on each loop of Gummel’s iteration. Based on our derivations, it can be shown that the Scharfetter–Gummel discretization method is essentially a special case of our scheme. In addition, the consistency and convergence of the two core algorithms are proved, with three numerical examples designed to demonstrate the accuracy and computational performance of this solver. Finally, we validate the Fr-DD model for a steady-state organic field effect transistor (OFET) by fitting the simulated transconductance and output curves to the experimental data. This work aimed to develop a Fr-DD model solver for simulating the anomalous dynamics of OSC devices. Two algorithms based on a novel discretization scheme and successively over-relaxed Gummel’s iteration are proposed here to solve the transient and steady-state Fr-DD model equations. This study has identified the consistency of the two algorithms by showing that the truncation error from the discretized divergence of current density functions will vanish with the spatial step size to a positive fractional order of 1 − β . The convergence analysis reveals that the Gummel mapping is a contraction mapping if we consider the successive over-relaxation mechanism in the Gummel’s iteration, which thus completes the proof of convergence. Three numerical examples, including one benchmark example and two others constructed from the perspective of engineering applications are employed to demonstrate the algorithms’ accuracy and computational performance. It is found in the first example that altering α and β can impact the spatial convergence order but only varying α will affect the convergence order in time. The increase rate of CPU time is less than the shrinking rate of temporal step size and lower than the growth rate of spatial grid points. These findings suggest that our solver has high precision and fast computational speed as it limits the computational error to a predefined satisfactory level (from ∼ 10 − 4 to ∼ 10 − 6 ) at a relatively small expense of CPU time (from ∼ 20 s to ∼ 100 s). The results reported in two numerical examples reveal the prediction and characterization of the transient-state and steady-state dynamics for any type of organic semiconductor device. Finally, we provide experimental verification for the fractional drift diffusion model of a DNTT-OFET. We stipulate that this is the first to date exploration of the Fr-DD model solver laying the groundwork for future research into fractional drift diffusion modeling of flexible organic electronics.

Accurate Description of Coulombic Interactions in Organic Field‐Effect Transistors Enabled by Efficient 3D Poisson's Equation Solver with Mixed Boundary Conditions

A Transfer Method for High‐Mobility, Bias‐Stable, and Flexible Organic Field‐Effect Transistors

Designing Contact Independent High‐Performance Low‐Cost Flexible Electronics

Semiconductor/dielectric Interface in Organic Field-Effect Transistors: Charge Transport, Interfacial Effects, and Perspectives with 2D Molecular Crystals

Interface Collaborative Strategy for High Mobility Organic Single‐Crystal Field‐Effect Transistors with Ideal Current–Voltage Curves

A Fractional Drift Diffusion Model for Organic Semiconductor Devices

Charge Distribution and Contact Resistance Model for Coplanar Organic Field-Effect Transistors

An efficient model algorithm for two-dimensional field-effect transistors

Precise Extraction of Charge Carrier Mobility for Organic Transistors

A Generalized Quantum Chemical Approach for Modeling of Single Molecular Field Effect Transistors

Recent Efforts in Understanding and Improving the Nonideal Behaviors of Organic Field‐Effect Transistors

Pursuing Polymer Dielectric Interfacial Effect in Organic Transistors for Photosensing Performance Optimization.

Refined drift-diffusion model for the simulation of charge transport across layer interfaces in organic semiconductor devices

Multi-scale Calculation of the Electric Properties of Organic-Based Devices from the Molecular Structure

De Novo Calculation of the Charge Carrier Mobility in Amorphous Small Molecule Organic Semiconductors

Interfacial Flat-Lying Molecular Monolayers for Performance Enhancement in Organic Field-Effect Transistors

Contact Resistance in Organic Field‐Effect Transistors: Conquering the Barrier

On the Nature of Charge-Injecting Contacts in Organic Field-Effect Transistors

Toward Efficient Charge Transport of Polymer-Based Organic Field-Effect Transistors: Molecular Design, Processing, and Functional Utilization

Charge-Trapping-Induced Non-Ideal Behaviors in Organic Field-Effect Transistors.