Abstract
A methodology is proposed that is suitable for efficient simulation of continuous-time Markov chains that are nearly-completely decomposable. For such Markov chains the effort to adequately explore the state space via Crude Monte Carlo (CMC) simulation can be extremely large. The purpose of this paper is to provide a fast alternative to the standard CMC algorithm, which we call Aggregate Monte Carlo (AMC). The idea of the AMC algorithm is to reduce the jumping back and forth of the Markov chain in small subregions of the state space. We accomplish this by aggregating such problem regions into single states. We discuss two methods to identify collections of states where the Markov chain may become ‘trapped’: the stochastic watershed segmentation from image analysis, and a graph-theoretic decomposition method. As a motivating application, we consider the problem of estimating the charge carrier mobility of disordered organic semiconductors, which contain low-energy regions in which the charge carrier can quickly become stuck. It is shown that the AMC estimator for the charge carrier mobility reduces computational costs by several orders of magnitude compared to the CMC estimator.
Similar content being viewed by others
References
Angulo J, Jeulin D (2007) Stochastic watershed segmentation. In: Proceedings of the 8th international symposium on mathematical morphology, pp 265–276
Bässler H (1993) Charge transport in disordered organic photoconductors. Phys Status Solidi, B 175(15):15–56
Baumeier B, May F, Lennartz C, Andrienko D (2012) Challenges for in silico design of organic semiconductors. J Mater Chem 22:10971–10976
Beucher S, Lantuéjoul C (1979) Use of watersheds in contour detection. In: Proceedings of the international workshop on image processing, real-time edge and motion detection/estimation, pp 2.1–2.12
Beucher S, Meyer F (1993) The morphological approach to segmentation: the watershed transformation. In: Mathematical morphology in image processing, vol 34, pp 433–481
Brereton T, Stenzel O, Baumeier B, Schmidt V, Kroese DP (2012) Efficient simulation of charge transport in deep-trap media. In: Laroque C, Himmelspach J, Pasupathy R, Rose O, Uhrmacher AM (eds) Proceedings of the 2012 winter simulation conference
Choi H, Szyld DB (1996) Application of threshold partitioning of sparse matrices to Markov chains. In: Proceedings of the IEEE international computer performance and dependability symposium, IPDS’96. IEEE Computer Society Press, pp 158–165
Conway A, Georganas ND (1982) Queueing networks—exact computational algorithms. MIT Press, Cambridge
Courtois P (1977) Decomposability. ACM monograph series. Academic Press, New York
Diestel R (2005) Graph theory. Springer, Berlin
Evans S (1996) Analysing system behaviour on different time scales. In: Kelly F, Zachary S, Ziedins I (eds) Stochastic networks—theory and applications. Clarendon Press, Oxford, pp 231–246
Faessel M, Jeulin D (2010) Segmentation of 3d microtomographic images of granular materials with the stochastic watershed. J Microsc 239(1):17–31
Jansen A (2012) An introduction to kinetic Monte Carlo simulations of surface reactions. Springer, Berlin
Pasveer W, Cottaar J, Tanase C, Coehoorn R, Bobbert P, Blom P, de Leeuw D, Michels M (2005) Unified description of charge-carrier mobilities in disordered semiconducting polymers. Phys Rev Lett 94:206601-1–206601-4
Rühle V, Lukyanov A, May F, Schrader M, Vehoff T, Kirkpatrick J, Baumeier B, Andrienko D (2011) Microscopic simulations of charge transport in disordered organic semiconductors. J Chem Theory Comput 7(10):3335–3345
Schönherr G, Bässler H, Silver M (1981) Dispersive hopping transport via sites having a Gaussian distribution of energies. Philos Mag B 44(1):47–61
Silverman B (1986) Density estimation for statistics and data analysis. Chapman and Hall, London
Simon H, Ando A (1961) Aggregation of variables in dynamic systems. Econometrica 29(2):111–138
Somoza A, Ortuño M (2005) Monte Carlo method for relaxation in electron glasses. Phys Rev B 72:224202-1–224202-6
Tessler N, Preezant Y, Rappaport N, Roichman Y (2009) Charge transport in disordered organic materials and its relevance to thin-film devices: a tutorial review. Adv Mater 21:2741–4761
Tse D, Gallager R, Tsitsiklis J (1995) Statistical multiplexing of multiple time-scale Markov streams. IEEE J Sel Areas Commun 13(6):1028–1038
van der Holst J, van Oost F, Coehoorn R, Bobbert P (2011) Monte Carlo study of charge transport in organic sandwich-type single-carrier devices: effects of Coulomb interactions. Phys Rev B 83:085206-1–085206-13
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Brereton, T., Stenzel, O., Baumeier, B. et al. Efficient Simulation of Markov Chains Using Segmentation. Methodol Comput Appl Probab 16, 465–484 (2014). https://doi.org/10.1007/s11009-013-9327-x
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11009-013-9327-x
Keywords
- Markov chain
- Nearly-completely decomposable
- Monte Carlo
- Segmentation
- Watershed
- Graph-theoretic decomposition
- Electron transport
- Mobility
- Organic semiconductor