Eur. Phys. J. B (2017) 90: 136
https://doi.org/10.1140/epjb/e2017-80145-1

Regular Article

Structure-correlated diffusion anisotropy in nanoporous channel networks by Monte Carlo simulations and percolation theory

¹ Institut für Theoretische Physik, Justus-Liebig-Universität Giessen, 35392 Giessen, Germany
² Felix Bloch Institute for Solid State Physics, University of Leipzig, 04103 Leipzig, Germany

^a e-mail: valiullin@uni-leipzig.de

Received: 11 March 2017
Received in final form: 26 May 2017
Published online: 17 July 2017

Abstract

Nanoporous silicon consisting of tubular pores imbedded in a silicon matrix has found many technological applications and provides a useful model system for studying phase transitions under confinement. Recently, a model for mass transfer in these materials has been elaborated [Kondrashova et al., Sci. Rep. 7, 40207 (2017)], which assumes that adjacent channels can be connected by “bridges” (with probability p_bridge) which allows diffusion perpendicular to the channels. Along the channels, diffusion can be slowed down by “necks” which occur with probability p_neck. In this paper we use Monte-Carlo simulations to study diffusion along the channels and perpendicular to them, as a function of p_bridge and p_neck, and find remarkable correlations between the diffusivities in longitudinal and radial directions. For clarifying the diffusivity in radial direction, which is governed by the concentration of bridges, we applied percolation theory. We determine analytically how the critical concentration of bridges depends on the size of the system and show that it approaches zero in the thermodynamic limit. Our analysis suggests that the critical properties of the model, including the diffusivity in radial direction, are in the universality class of two-dimensional lattice percolation, which is confirmed by our numerical study.