There is a large variety of quantum and classical systems in which the quenched disorder plays a dominant r\^ole over quantum, thermal, or stochastic fluctuations : these systems display strong spatial heterogeneities, and many averaged observables are actually governed by rare regions. A unifying approach to treat the dynamical and/or static singularities of these systems has emerged recently, following the pioneering RG idea by Ma and Dasgupta and the detailed analysis by Fisher who showed that the Ma-Dasgupta RG rules yield asymptotic exact results if the broadness of the disorder grows indefinitely at large scales. Here we report these new developments by starting with an introduction of the main ingredients of the strong disorder RG method. We describe the basic properties of infinite disorder fixed points, which are realized at critical points, and of strong disorder fixed points, which control the singular behaviors in the Griffiths-phases. We then review in detail applications of the RG method to various disordered models, either (i) quantum models, such as random spin chains, ladders and higher dimensional spin systems, or (ii) classical models, such as diffusion in a random potential, equilibrium at low temperature and coarsening dynamics of classical random spin chains, trap models, delocalization transition of a random polymer from an interface, driven lattice gases and reaction diffusion models in the presence of quenched disorder. For several one-dimensional systems, the Ma-Dasgupta RG rules yields very detailed analytical results, whereas for other, mainly higher dimensional problems, the RG rules have to be implemented numerically. If available, the strong disorder RG results are compared with another, exact or numerical calculations.Comment: review article, 195 pages, 36 figures; final version to be published in Physics Report

The effect of quenched disorder on nonequilibrium phase transitions in the directed percolation universality class is studied by a strong disorder renormalization group approach and by density matrix renormalization group calculations. We show that for sufficiently strong disorder the critical behavior is controlled by a strong disorder fixed point and in one dimension the critical exponents are conjectured to be exact: beta=(3-sqrt[5])/2 and nu( perpendicular )=2. For disorder strengths outside the attractive region of this fixed point, disorder dependent critical exponents are detected. Existing numerical results in two dimensions can be interpreted within a similar scenario.

The phase transitions and critical properties of two types of inhomogeneous systems are reviewed. In one case, the local critical behaviour results from the particular shape of the system. Here scale-invariant forms like wedges or cones are considered as well as general parabolic shapes. In the other case the system contains defects, either narrow ones in the form of lines or stars, or extended ones where the couplings deviate from their bulk values according to power laws. In each case the perturbation may be irrelevant, marginal or relevant. In the marginal case one finds local exponents which depend on a parameter. In the relevant case unusual stretched exponential behaviour and/or local first order transitions appear. The discussion combines mean field theory, scaling considerations, conformal transformations and perturbation theory. A number of examples are Ising models for which exact results can be obtained. Some walks and polymer problems are considered, too.