The random propagation of molecules in a fluid medium is characterized by the spontaneous diffusion law as well as the interaction between the environment and molecules. In this paper, we embody the anomalous diffusion theory for modeling and analysis in molecular communication. We employ H-diffusion to model a non-Fickian behavior of molecules in diffusive channels. H-diffusion enables us to model anomalous diffusion as the subordinate relationship between self-similar parent and directing processes and their corresponding probability density functions with two H-variates in a unified fashion. In addition, we introduce standard H-diffusion to make a bridge of normal diffusion across wellknown anomalous diffusions such as space-time fractional diffusion, Erdélyi-Kober fractional diffusion, grey Brownian motion, fractional Brownian motion, and Brownian motion. We then characterize the statistical properties of uncertainty of the random propagation time of a molecule governed by Hdiffusion laws by introducing a general class of molecular noise-called H-noise. Since H-noise can be an algebraic tailed process, we provide a concept of H-noise power using finite logarithm moments based on zero-order statistics. Finally, we develop a unifying framework for error probability analysis in a timing-based molecular communication system with a concept of signal-to-noise power ratio. ). 2 3 for timing and amplitude modulations. This work further extended to a connectivity problem with a random time constraint in a one-dimensional nanonetwork, where the random locations of molecules at the initial time are modeled by poisson point process [20]. The Cox process has been considered in [21] to capture the dynamic variation of the molecule concentration arising from the mobility of anomalously diffusive molecules, and the spatial ordering of the molecular communication performance has been characterized in terms of the error rate in the presence of interfering molecules. However, there is no comprehensive study on the modeling of anomalous diffusion channels in the context of molecular communication. Various anomalous diffusion processes typically can be modeled numerous ways including continuous random walk (CTRW), generalized diffusion equation, generalized master equation, fractional Brownian motion, and fractional kinetic equation (fractional diffusion equation) [15], [23]-[25]. 3 In particular, the CTRW simply describes diffusion of molecules in the medium with arbitrary distributions of jump lengths and waiting times. 4 In addition, the combination of a stochastic operational time-a directing process-and the self-similar parent process is equivalent to the subordination integral mechanism for the product of two random variables in the context of subordinated processes [28]-[31]. 5 This subordination law generates the solution of the fractional diffusion equation in purely analytical ways using the machinery offered by convolution properties of the Mellin transform. 6In this paper, we embody anomalous diffusion according to th...