# The following file is created by the region editor, based on 'donor_acceptor.sit'.
# This mimics a situation like in [1]
#
# [1]  K. Maturová, S. S. van Bavel, M. M. Wienk, R. A. J. Janssen and M. Kemerink,
#      "Morphological Device Model for Organic Bulk Heterojunction Solar Cells" 
#      Nano Lett., 2009, 9 (8), pp 3032-3037
#      http://pubs.acs.org/doi/abs/10.1021/nl901511a
#
import donor_acceptor.txt

# Here, looking for the equilibrium works best if we don't initialize
# n and p using the sim2/grid/init command, but leave the initial value
# at zero.
sim2/rundirect

# Create the JV curve
sim2/iv VBI 0 1.0 100 jv.txt no