Massively parallel kinetic Monte Carlo simulations of charge carrier transport in organic semiconductors

A parallel, lattice based Kinetic Monte Carlo simulation is developed that runs on a GPGPU board and includes Coulomb like particle-particle interactions. The performance of this computationally expensive problem is improved by modifying the interaction potential due to nearby particle moves, instead of fully recalculating it. This modification is achieved by adding dipole correction terms that represent the particle move. Exact evaluation of these terms is guaranteed by representing all interactions as 32-bit floating numbers, where only the integers between - 2 22 and 222 are used. We validate our method by modelling the charge transport in disordered organic semiconductors, including Coulomb interactions between charges. Performance is mainly governed by the particle density in the simulation volume, and improves for increasing densities. Our method allows calculations on large volumes including particle-particle interactions, which is important in the field of organic semiconductors.

[1]  Michael A. Gibson,et al.  Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels , 2000 .

[2]  Adrian Neagu,et al.  Kinetic Monte Carlo and cellular particle dynamics simulations of multicellular systems. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[3]  P. Clancy,et al.  Kinetic Monte Carlo simulation of the growth of polycrystalline Cu films , 2001 .

[4]  Anna Köhler,et al.  Charge transport in organic semiconductors. , 2012, Topics in current chemistry.

[5]  Kristen A Fichthorn,et al.  Synchronous relaxation algorithm for parallel kinetic Monte Carlo simulations of thin film growth. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[6]  Richard W. Vuduc,et al.  A massively parallel adaptive fast-multipole method on heterogeneous architectures , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.

[7]  P. Smereka,et al.  An energy localization principle and its application to fast kinetic Monte Carlo simulation of heteroepitaxial growth , 2009 .

[8]  Enrique Martínez,et al.  Synchronous parallel kinetic Monte Carlo for continuum diffusion-reaction systems , 2008, J. Comput. Phys..

[9]  Boris D. Lubachevsky,et al.  Efficient parallel simulations of dynamic Ising spin systems , 1988 .

[10]  A. B. Bortz,et al.  A new algorithm for Monte Carlo simulation of Ising spin systems , 1975 .

[11]  Roger Smith,et al.  Long time scale evolution of radiation-induced defects in Er 2 O 3 , 2011 .

[12]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

[13]  Guido Raos,et al.  Methodological assessment of kinetic Monte Carlo simulations of organic photovoltaic devices: the treatment of electrostatic interactions. , 2010, The Journal of chemical physics.

[14]  Pa Peter Bobbert,et al.  Calculating charge-carrier mobilities in disordered semiconducting polymers: Mean field and beyond , 2006 .

[15]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[16]  Lemi Demeyu,et al.  Monte Carlo simulations of charge carrier mobility in semiconducting polymer field-effect transistors , 2007 .

[17]  Michael Engel,et al.  Massively parallel Monte Carlo for many-particle simulations on GPUs , 2012, J. Comput. Phys..

[18]  Ignacio Martin-Bragado,et al.  MMonCa: An Object Kinetic Monte Carlo simulator for damage irradiation evolution and defect diffusion , 2013, Comput. Phys. Commun..

[19]  H. Bässler.  Charge Transport in Disordered Organic Photoconductors a Monte Carlo Simulation Study , 1993 .

[20]  N. Greenham,et al.  Effect of charge trapping on geminate recombination and polymer solar cell performance. , 2010, Nano letters.

[21]  Enrique Martínez,et al.  Billion-atom synchronous parallel kinetic Monte Carlo simulations of critical 3D Ising systems , 2010, J. Comput. Phys..

[22]  M. Pritzker,et al.  Kinetic Monte Carlo simulation of electrodeposition using the embedded-atom method , 2013, 1307.5586.

[23]  Jean-Luc Brédas,et al.  Charge transport in organic semiconductors. , 2007, Chemical reviews.

[24]  R. Coehoorn,et al.  Monte Carlo study of charge transport in organic sandwich-type single-carrier devices: Effects of Coulomb interactions , 2011 .

[25]  Michail Stamatakis,et al.  First-Principles-Based Kinetic Monte Carlo Simulation of the Structure Sensitivity of the Water–Gas Shift Reaction on Platinum Surfaces , 2011 .

[26]  Babak Sadigh,et al.  MECHANISM OF BORON DIFFUSION IN SILICON : AN AB INITIO AND KINETIC MONTE CARLO STUDY , 1999 .

[27]  Richard L. Martin,et al.  Molecular geometry fluctuations and field-dependent mobility in conjugated polymers , 2001 .

[28]  L. Greengard,et al.  Regular Article: A Fast Adaptive Multipole Algorithm in Three Dimensions , 1999 .

[29]  L. Berthier,et al.  Clustering and heterogeneous dynamics in a kinetic Monte Carlo model of self-propelled hard disks. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[30]  W. H. Weinberg,et al.  Theoretical foundations of dynamical Monte Carlo simulations , 1991 .

[31]  Jacques G. Amar,et al.  Semirigorous synchronous sublattice algorithm for parallel kinetic Monte Carlo simulations of thin film growth , 2005 .

[32]  Jacques G. Amar,et al.  Synchronous relaxation algorithm for parallel kinetic Monte Carlo , 2004 .

[33]  R. Ziff,et al.  Kinetic phase transitions in an irreversible surface-reaction model. , 1986, Physical review letters.

[34]  Denis Andrienko,et al.  Extracting nondispersive charge carrier mobilities of organic semiconductors from simulations of small systems , 2010 .

[35]  Alain Barbu,et al.  Multiscale modelling of defect kinetics in irradiated iron , 2004 .

[36]  Notker Rösch,et al.  Ethylene conversion to ethylidyne on Pd(111) and Pt(111): A first-principles-based kinetic Monte Carlo study , 2012 .

[37]  Xiaoyuan Hou,et al.  Carrier density dependence of mobility in organic solids: A Monte Carlo simulation , 2007 .

[38]  Steven J. Plimpton,et al.  Crossing the Mesoscale No-Man's Land via Parallel Kinetic Monte Carlo , 2009 .

[39]  P. Smereka,et al.  Kinetic Monte Carlo simulation of heteroepitaxial growth: Wetting layers, quantum dots, capping, and nanorings , 2012 .

[40]  N Le Novère,et al.  Conformational spread in a ring of proteins: a stochastic approach to allostery. , 2001, Journal of molecular biology.

[41]  Petr Plechác,et al.  Hierarchical fractional-step approximations and parallel kinetic Monte Carlo algorithms , 2011, J. Comput. Phys..

[42]  P. Blom,et al.  Unified description of charge-carrier mobilities in disordered semiconducting polymers. , 2005, Physical review letters.

[43]  E. Abrahams,et al.  Impurity Conduction at Low Concentrations , 1960 .

[44]  A. P. J Jansen,et al.  An Introduction to Kinetic Monte Carlo Simulations of Surface Reactions , 2012 .