$callexciton

</math>Function for calculate the exciton distribution. We usually use this equation for organic material. Behavior of exciton will follow this equation. You can see the detail in Subroutine_exciton1D.

Singlet Rate Equation: $\frac{d{n}_{ex}^S}{dt}=D^S{\mathrm{\nabla }}^2{n}_{ex}^S-(k_r^S+k_{nr}^S+k_e^{S}n+k_h^{S}p+k_{TS}{n}_{ex}^T)n_{ex}^S+\alpha \frac{{\gamma }_{TS}{n_{ex}^T}^2}{2}+G_S$

Triplet Rate Equation: $\frac{d{n}_{ex}^T}{dt}=D^S{\mathrm{\nabla }}^2{n}_{ex}^T-(k_r^T+k_{nr}^T+k_e^{T}n+k_h^{T}p)n_{ex}^T-\frac{({\gamma }_{TS}+{\gamma }_{TT}){n_{ex}^T}^2}{2}+G_T$

Where　

• $D$ is diffusion coefficient.
• $\tau$ is relaxation time of exciton.
• $\gamma$ is annihilation rate constant.
• $G$ is exciton generation rate.

Format

$callexciton
n
a 4 b c d f
d kr knr gamma g

Parameter Explanation

• n : the number of tables we usually set n as 5.
• a : The type of exciton solver mode

 1: Time-dependent triplet solver 
 123: Time-dependent triplet and singlet solver (For TADF OLEDs model)
 3: Triplet Exciton Solver (For PhOLEDs model)
 6: Singlet and Triplet Exciton Solver (For TADF OLEDs model)
 4: Triplet Exciton Solver with exciton blocking boundary 

• b : Start time (For time-dependent solver)
• c : dt (For time-dependent solver)
• d : End time (For time-dependent solver)
• e : savenum (For time-dependent solver)
• D : diffusion coefficient. $(c{m}^2{s}^{-1})$
• kr : radiatvie rate constant $(s^{-1})$
• knr :non-radiative rate constant $(s^{-1})$
• gamma : quenching coefficient. $(c{m}^2{s}^{-1})$
• g : generation rate if you wanna let whole recombination rate change into exciton you should set g as 1.

Example

$callexciton
5
2e-14 20000 3000 1e-12 1
2e-14 20000 3000 1e-12 1
2e-14 20000 3000 1e-12 1
2e-14 20000 3000 1e-12 1
2e-14 20000 3000 1e-12 1

static TTA model (mode 7)

Format

$callexciton
20
7 1 1
DS DT krS knrS krT knrT kisc krisc keS khS keT khT kST gammaTS gammaTT a DrefS DrefT ES ET 

Parameter Explanation ...

Subroutine_exciton1D,