ShaneLogic/drift-diffusion-modeling

Drift-Diffusion Modeling Strategy

Apply this strategy when selecting a mathematical model to simulate Perovskite Solar Cell (PSC) behavior, especially when you need to balance computational tractability with physical interpretability.

When to Use

• Selecting a mathematical model for PSC simulation
• Designing charge transport models
• Addressing numerical stiffness in simulations with multiple timescales
• Implementing ion vacancy and electronic charge carrier motion

Procedure

1. Evaluate Modeling Approaches

Compare available methods for PSC simulation:

• DFT: High computational cost, limited to few atoms, short timescales
• Equivalent Circuit: Easy to solve, difficult to connect to device physics
• Drift-Diffusion: Intermediate computational cost, tractable, interpretable physics

2. Select Drift-Diffusion Model

Choose drift-diffusion when it meets your requirements:

• Incorporates parameters from DFT calculations
• Accounts for both ion vacancy motion and electronic charge carriers
• Provides interpretable connection to device physics

3. Address Numerical Challenges

Implement solutions for inherent stiffness issues:

Spatial Stiffness:

• Short Debye lengths create boundary layers with rapid solution variation
• Use non-uniform grids to handle boundary layer resolution

Temporal Stiffness:

• Large disparity in timescales between ion vacancies and charge carriers
• Implement adaptive time stepping to handle timescale differences

Realistic Conditions:

• Account for realistic ion densities (up to 10^19 cm^-3)
• These high densities exacerbate numerical stiffness

4. Apply Single-Layer Assumption

Simplify the model when appropriate:

• Treat highly doped transport layers as "quasi-metals"
• Assume uniform electric potential in transport layers
• Restrict physics calculations to the perovskite layer only

Key Concepts

• Debye Length: Scale over which charge screens electric field
• Stiffness: Numerical property of systems with widely varying timescales or lengthscales
• Ion Density Threshold: Up to 10^19 cm^-3 for realistic conditions