On the scale-invariant distribution of the diffusion coefficient for classical particles diffusing in disordered media

The scaling form of the whole distribution P(D) of the random diffusion coefficient D(x) in a model of classically diffusing particles is investigated. The renormalization group approach above the lower critical dimension d=0 is applied to the distribution P(D) using the n-replica approach. In the annealed approximation (n=1), the inverse gaussian distribution is found to be the stable one under rescaling. This identification is made based on symmetry arguments and subtle relations between this model and that of fluc- tuating interfaces studied by Wallace and Zia. The renormalization-group flow for the ratios between consecutive cumulants shows a regime of pure diffusion for small disorder, in which P(D) goes to delta(D-<D>), and a regime of strong disorder where the cumulants grow infinitely large and the diffusion process is ill defined. The boundary between these two regimes is associated with an unstable fixed-point and a subdiffusive behavior: <x**2>=Ct**(1-d/2). For the quenched case (n goes to 0) we find that unphysical operators are generated raisng doubts on the renormalizability of this model. Implications to other random systems near their lower critical dimension are discussed.