$callexciton

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.

$\frac{d{n}_{ex}}{dt}=D{\mathrm{\nabla }}^2{n}_{ex}(r)-\frac{n_{ex}(r)}{\tau }-\gamma {n_{ex}(r)}^2+G$

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 b c d f g
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 2:

• 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

Subroutine_exciton1D,