# Comment on: Kinetic Roughening in Slow Combustion of Paper

PACS numbers: 64.60.Ht, 05.40.+j, 05.70.Ln

In a recent Letter, Maunuksela et al. present experiments on the combustion of paper . They observe the roughening of the burning front and study the scaling properties of the correlation function of the height of the burning front. Maunuksela et al. find good agreement between the measured exponents characterizing the front roughening and the predictions of the Kardar-Parisi-Zhang (KPZ) equation , $`\chi _{KPZ}=1/2`$ ($`\chi `$ is the roughness exponent which characterizes the scaling of the saturated height-height correlation function $`G(r)`$ with distance $`G(r)r^{2\chi }`$). Reference  also comments on the results of earlier experiments on paper burning by Zhang et al.  which measured $`\chi 0.71`$, but offers no explanation for the difference in the measured value of the exponent.

Here, we show that the results of Maunuksela et al.  and of Zhang et al.  may be both understood under the framework of interface motion in disordered media .

For many experimental cases, the dominant source of noise in the dynamics is the disorder in the medium which is not time dependent. The universality classes for interface motion in disordered media have been identified and the values of the exponents are known, especially for $`(1+1)`$ dimensions . One of the universality classes for interface motion in disordered media  can be described by the directed percolation depinning (DPD) model . For a driving force $`F`$ smaller than a critical value $`F_c`$, which defines the depinning transition, the interface eventually stops, and the scaling of $`G(r)`$ is characterized by an exponent $`\chi _{DPD}0.63`$ . If, on the other hand, $`F>F_c`$, the interface moves with an average velocity $`v(FF_c)^\theta `$ (where $`\theta `$ is a critical exponent), and $`G(r)`$ presents two scaling regimes in the steady-state. For length-scales smaller than the correlation length given by the disorder of the medium $`\xi (FF_c)^\nu _{}`$ (where $`\nu _{}`$ is a critical exponent), $`G(r)`$ scales with an exponent $`\chi _{DPD}^{}0.75`$, while for length-scales larger than $`\xi `$, $`G(r)`$ scales with the KPZ exponent, $`\chi _{KPZ}=1/2`$ .

To check that the results reported by Maunuksela et al. are consistent with the above theory, we digitize the data reported in  and plot it in Fig. 1 along with the predictions of the DPD model (see Refs. ). It is visually apparent that the experimental and theoretical data sets have the same two scaling regimes as discussed in Ref. . Furthermore, the same theory may explain why Zhang et al. found no crossover to a KPZ dominated regime. Zhang’s experiments were performed near the depinning transition, which implies that $`\xi `$ was nearly as large as the system size. Thus, the quenched disorder dominates the scaling of $`G(r)`$ leading to the observation of only one scaling regime with a roughness exponent $`\chi \chi _{DPD}^{}`$. On the other hand, in the experiments reported by Maunuksela et al., a shorter correlation length $`\xi `$ is found, so that quenched disorder dominates only for small length scales where an exponent with values close to $`\chi _{DPD}^{}`$ was measured. For length scales larger than $`\xi `$, time-dependent disorder dominates, and the KPZ exponent $`\chi _{KPZ}`$ was measured, in agreement with the theory . These results suggest that the DPD universality class may provide a compelling explanation of the experimental findings of Refs. .