Experimental approaches to universal outofequilibrium scaling laws: turbulent liquid crystal and other developments
Abstract
This is a brief survey of recent experimental studies on outofequilibrium scaling laws, focusing on two prominent situations where nontrivial universality classes have been identified theoretically: absorbingstate phase transitions and growing interfaces. First the article summarizes main results obtained for electricallydriven turbulent liquid crystal, which exhibited the scaling laws for the directed percolation class at the transition between two turbulent regimes, and those for the KardarParisiZhang class in the supercritical phase where one turbulent regime invades the other. Other experimental investigations on these universality classes and related situations are then overviewed and discussed. Some remarks on analyses of these scaling laws are also given from the practical viewpoints.
1 Introduction
Physics of critical phenomena and other scaleinvariant problems like the Brownian motion is a landmark in statistical mechanics, which has provided us with simple and unified views on diverse phenomena in nature [1, 2]. Scale invariance refers to absence of characteristic scales. It is known to arise, typically, when a system is undergoing a secondorder phase transition. A prototypical example is the Ising model. While patches or islands of up/down spins with finite characteristic lengths are formed in the ferro and paramagnetic phases, diverges at the critical temperature, leading to fractal spin configuration where fluctuations take place in all scales, so that scale invariance arises. Descriptions of such systems should then also be scaleinvariant, wiping out most microscopic details like the choice of material or the mechanism of interactions as one moves toward larger scales of interest. Roughly, this is how universality arises in such scaleinvariant phenomena, often in terms of power laws characterized by universal scaling (critical) exponents – the idea formulated and established notably by Wilson by means of renormalization group theory [3]. These scaling laws usually depend solely on global properties of the systems, such as spatial dimensions, symmetry, and conservation laws, which group various systems into a few universality classes that share the same critical exponents and scaling functions. This is now deeply understood for equilibrium systems by theoretical frameworks like renormalization group and conformal field theory [1, 2], and, importantly, firmly underpinned by ample experimental evidence: for example, the Isingclass criticality has been observed in real experiments of magnets, liquidvapor systems and alloys, to name but a few [2].
It is then natural to expect similar universality
for systems driven out of equilibrium,
since the concept does not seem to necessitate being at equilibrium.
Indeed, theoretical developments have adapted renormalization group
and fieldtheoretic techniques to outofequilibrium systems,
identifying a number of universality classes
that have no counterpart in equilibrium, hence truly out of equilibrium
[4, 5, 6, 7].
Two noteworthy classes in this context are the directed percolation (DP) class,
governing continuous phase transitions into an absorbing state
[4, 5, 6],
and the KardarParisiZhang (KPZ) class for generic nonlinear growth processes
[4, 8, 7, 9, 10, 11]
In the author’s view, main difficulties in studying such outofequilibrium scaling laws experimentally stem from a chief characteristic of outofequilibrium systems: formation of dissipative structure. For instance, the RayleighBénard convection forms convective rolls and other types of structure, which behave as single dynamical units when they start to oscillate or fluctuate [12, 13]. The effective system size is then given by the aspect ratio of the system, typically in the order of tens or a hundred at most for fluid systems, which is much smaller than the Avogadro number. This makes truly macroscopic behavior difficult to access. Similarly, since each degree of freedom has a relatively large scale, thermal noise is in most cases negligible. Stochastic processes then result from other sources of effective noise, such as chaos or external heterogeneity, which however tend to introduce complex and possibly longrange noise correlations, as well as quenched disorder. It is well known that these factors can affect or even destroy the scaling laws that one would expect in their absence [5, 6, 7]. In addition, the presence of ordered structure can easily introduce effectively longrange interactions to the system, which may also have a decisive impact on the scaling laws.
To overcome all these difficulties, in a series of papers [14, 15, 16, 17, 18] the author and his coauthors have focused on the electroconvection of nematic liquid crystal [19, 20]. This is a system much studied in the context of pattern formation and is argued to be analogous to the RayleighBénard convection, except that liquid crystal molecules are driven by an electric field through the CarrHelfrich instability [19]. Thanks to its efficient driving mechanism, the electroconvection is usually studied in a very thin container of liquid crystal, which implies a very large system size. Moreover, working with turbulent regimes called the dynamic scattering modes 1 and 2 (DSM1 and DSM2) [19, 20], which are actually regimes of spatiotemporal chaos with exponentially decaying correlation in space and time, one can safely eliminate the effect of longrange correlations. Finally, the source of noise in this system is intrinsic turbulent fluctuations, which overwhelm uncontrolled heterogeneities of the experimental setup. In this simple yet nontrivial system for studying outofequilibrium scaling laws, the author and coworkers have investigated a phase transition between the DSM1 and DSM2 states [14, 15] as well as growth processes in the DSM2 phase, where stable DSM2 clusters invade the metastable DSM1 phase [16, 17, 18].
The article is organized as follows. First, section 2 summarizes main results obtained for the DSM1DSM2 transition and for the DSM2 growth processes. As will be explained there, the former turned out to exhibit the DPclass criticality [14, 15], while the latter shows the scaling laws of the KPZ class [16, 17, 18]. Section 3 overviews other experimental studies related to these two problems, with more weight on recent active developments. In the hope of facilitating these directions of research, section 4 describes some general remarks and practical hints on analyses of these scaling laws, followed by concluding remarks in section 5. Note, however, that this article is not intended to be a comprehensive review of these matters: there are many experimental studies that the author was not able to cover in this article, especially on the case where they are related in some indirect manner, or they fit better to different theoretical contexts (e.g., cases with longrange interactions, quenched disorder, etc.). The author also reminds the readers that the article deals with scaling laws associated with absorbingstate transitions and growth processes, leaving aside other interesting outofequilibrium scaling laws, in particular those for fully developed turbulence [21].
2 Main results of the turbulent liquid crystal experiment
This section describes main experimental results obtained for the liquidcrystal turbulence. Specifically, section 2.1 deals with the DPclass critical behavior found at the DSM1DSM2 transition [14, 15], and section 2.2 shows the KPZclass interface fluctuations identified in the DSM2 growth processes [16, 17, 18]. For complete descriptions, the readers are referred to the original papers, in particular [15] and [18].
The standard experimental setup for both studies is briefly outlined in figure 1(a). A thin container of inner dimensions 16 mm 16 mm 12 m was filled with nematic liquid crystal (4methoxybenzylidene)4butylaniline (MBBA) doped with tetrabutylammonium bromide. The inner surfaces of the container were treated so that liquid crystal molecules are aligned in a desired direction on the surfaces, here along an axis parallel to the surfaces, unless otherwise specified. One can apply an alternating electric field perpendicularly to the sample, using transparent electrodes coated on the surfaces. Here, a low frequency of the field was chosen, specifically 250 Hz, in order to work within the conductive regime of the electroconvection [19, 20]. The convection was observed through the transmitted light, emitted from lightemitting diodes and captured by a chargecoupled device camera. The temperature of the sample was kept constant with high precision throughout each experiment, by means of a thermocontroller and heatinsulating walls (the latter not shown in figure 1(a)). Lastly, in the presented setup, one can shoot ultraviolet laser pulses to the liquid crystal sample, which are used to trigger local nucleation of the DSM2 state in the DSM1 region [15].
In the lowfrequency region chosen here, one can observe a variety of patterns such as rolls and grids as the amplitude of the applied voltage is increased, finally reaching the regimes of spatiotemporal chaos, DSM1 and DSM2 [19, 20], separated by the critical voltage . The difference between DSM1 and DSM2 lies in their density of topological defects, or specifically, disclinations, in the liquidcrystal orientation (figure 1(b)). The DSM2 state consists of a high density of disclinations [22, 23], which leads, in macroscopic scales, to a lower light transmittance and to the loss of nematic anisotropy. Spontaneous nucleation of the DSM2 state is therefore prevented by a high energy barrier for formation of topological defects. While for such events do occur in realistic time scales [22, 23], followed by domain growth of the stable DSM2 phase, for voltages near the threshold, spontaneous nucleation of the DSM2 state becomes a practically unobservable event. The dynamics is then characterized by spatiotemporal intermittency of the DSM2 state, where dark DSM2 patches randomly migrate amid the bright DSM1 region (figure 1(c) and video 1), sometimes splitting, sometimes disappearing. This stochastic process results from random transport, stretching, and shrinkage of disclinations, due to local chaotic flow of the electroconvection. Note also that, for the DSM2 state to be formed and maintained, constituent disclinations should break surface anchoring of the liquidcrystal alignment [24] by reaching both the top and bottom surfaces of the container. Therefore, despite the quasitwodimensional geometry of the container, the DSM1DSM2 configuration is defined in a purely twodimensional space. This also implies that the minimum size of a single DSM2 domain is in the order of the cell thickness , specifically [22], which would give a cutoff length scale for coarsegrained descriptions of the DSM1DSM2 dynamics at stake.
2.1 Directed percolation class governing the DSM1DSM2 transition [14, 15]
Here we overview the spatiotemporal intermittency described above, which characterizes the DSM1DSM2 transition. The two states can be distinguished by the light intensity (see video 1), but for an automatic binarization one should also take into account that the DSM1DSM2 configuration evolves more slowly than local intensity fluctuations and that DSM2 domains cannot be smaller than [22] (see [15] for details). Figure 1(d) shows results of such binarization, indicating how DSM2 patches evolve in space and time in the steady state, at some voltages slightly larger than . This shows that, for larger , more space is occupied by DSM2. In contrast, for lower than the threshold, all DSM2 patches eventually disappear. Therefore, one can choose the DSM2 area fraction, denoted by , as an order parameter for this transition.
Figure 2(a) shows the timeaveraged order parameter against the applied voltage in the steady state, clearly indicating that the transition is continuous. We therefore assume the usual scaling form , or, instead, an asymptotically equivalent one
(1) 
following convention for the electroconvection, which takes into account its driving force proportional to [19, 20]. The data indeed showed a clear power law over 3 digits, providing a precise estimate of the critical voltage at and the critical exponent , as shown in the inset of figure 2(a). Here, the numbers in the parentheses indicate the range of errors in the last digit. The estimate for therefore rules out, e.g., the Landau meanfield exponent , but instead agrees with the exponent for the DP class in dimensions, [25, 26, 5, 6]. Similarly, the correlation length and time can be measured from DSM1 interval distributions in the steady state [15]. The length or the duration separating two neighboring DSM2 domains shows powerlaw distribution at criticality, and , while for it is cut off by an exponential tail as and for large and . From this, and were found to show powerlaw divergence at the critical point (figures 2(c,d)) as follows:
(2) 
The data yielded and , both in agreement with the DPclass values and . Moreover, the values of the exponents and for the powerlaw distributions also supported the DPclass scenario, and [15], to be compared with and . This also confirms the following scaling relations:
(3) 
One can also study dynamic critical behavior by, e.g.,
suddenly decreasing the applied voltage from a value deep in the DSM2 phase
to a target value near the critical point.
Figure 2(b) shows relaxation of the order parameter
after such critical quenching, which indicates powerlaw decay
with at criticality
(red bold curve)
(4) 
with and a universal scaling function . This gives and correctly accounts for the experimental result. Moreover, the scaling form in equation (4) implies that data for different voltages should collapse on the single universal function , when is plotted against . This is indeed the case as shown in the inset of figure 2(b), where the data in the main panel overlap onto the DPclass universal scaling function obtained independently by a numerical simulation. The result here substantiates, therefore, both the scaling relation and the universality in the scaling function .
It is also worthwhile to study another aspect of dynamic critical behavior, by shooting laser pulses to create a single DSM2 nucleus in the fully DSM1 state and track the resultant DSM2 cluster. Power laws at criticality were then identified for the probability that the cluster survives until time , as well as for the mean volume and the mean squared radius of the cluster, with the corresponding exponents of the DP universality class [15]. Actually, that survival probability is known to serve as another order parameter for the DPclass and related phase transitions, whose exponent is generally different from . For the DP class, however, a special timereversal symmetry called the rapidity symmetry leads to [5, 6]. The liquidcrystal experiment therefore underpins the presence of this nontrivial symmetry in a real system. In such ways, the three sets of experiments for the steady state, critical quench, and critical spreading have revealed a total of 12 critical exponents, 5 scaling functions, and 8 scaling relations, in full agreement with those characterizing the dimensional DP class [14, 15].
Theoretically, the appearance of the DPclass criticality in this experiment can be accounted for on the basis of Pomeau’s general argument on spatiotemporal intermittency [27] and Janssen and Grassberger’s DP conjecture [28, 29, 5, 6]. Pomeau’s argument [27] is illustrated in figure 1(e): consider, for the sake of simplicity, a twostate lattice model with active or inactive local states, which correspond to the DSM2 and DSM1 states in the experiment, respectively. An active site can activate neighbors and/or relax into the inactive state, both stochastically at constant rates. This reproduces, at least qualitatively, spatiotemporal intermittency like the one observed in the liquidcrystal experiment (figure 1(c) and video 1) and is nothing but a percolation process directed in the time forward direction. Indeed, the model explained here, called the contact process, is recognized as a prototypical model exhibiting the DPclass transition [5, 6]. Of course it is an oversimplified model for the electroconvection, but then the DP conjecture implies that such differences do not change the critical behavior. Specifically, the DSM1DSM2 transition can be regarded as a transition into an absorbing state, i.e., a global state that the system stays in forever once it enters, since, in practice, the DSM1 state admits no spontaneous nucleation of the DSM2 state. The DP conjecture then states, in short, that phase transitions into a single absorbing state generically belong to the DP universality class, unless the system has additional symmetry, a conservation law, longrange interactions, or quenched disorder [28, 29, 5, 6]. This is true of the DSM1DSM2 transition as argued in section 1 and discussed in more detail in [15]. As a result, the DPclass criticality arises in such a real, microscopically complex system as well.
2.2 KardarParisiZhang class for the DSM2 growing interfaces [16, 17, 18]
So far we have reviewed critical behavior at the DSM1DSM2 transition, where a single cluster of the DSM2 state does not grow or shrink in an obvious way. This is because, in terms of Pomeau’s scenario in figure 1(e), the activation rate is essentially balanced with the recession rate. By contrast, when the applied voltage is sufficiently raised, the activation rate overwhelms the recession rate, leading to sustained growth of the DSM2 cluster bordered by a fluctuating interface (figures 3(a,b)). Microscopically, while disclinations inside the DSM2 cluster are maintained at a more or less constant density, on the border, strong chaotic flow continuously elongates, multiplies, and randomly transport disclinations, leading to the random local expansion of the DSM2 cluster in macroscopic scales. One can easily observe such growing clusters using spontaneous nucleations at high applied voltages, regardless of the chosen alignment of liquid crystal molecules. Under planar alignment, however, elliptic clusters are produced because of the anisotropy in the horizontal plane [22, 23]. A simple and controlled experiment can be performed, instead, by measuring isotropic clusters under homeotropic alignment, where molecules are aligned perpendicularly to the surfaces, and by triggering nucleation by laser pulses instead of relying on spontaneous nucleation. The use of laser also allows changing the cluster shape. Figures 3(a,b) show growing DSM2 clusters produced in such a way by point and lineshaped laser pulses, exhibiting circular and flat interfaces, respectively. One can therefore investigate statistical properties of these growing interfaces, as well as consequences of the different geometries.
Figures 3(a,b) indicate that the interfaces roughen in the course of time with fractallike intricate structure. This suggests scale invariance of the process, which may therefore be described by another set of universal scaling laws. This turned out to be indeed the case [16, 17, 18], characterized by the KPZ universality class, known as the simplest generic class for such stochastic growth processes due to shortrange interactions [8, 7, 9, 10, 11]. To characterize the interface fluctuations, one defines the local height as in figures 3(a,b), as a function of the lateral coordinate and time . Roughness is then quantified by the standard deviation of , measured within a strip of lateral length and averaged over space and ensemble, denoted by . For scaleinvariant growth, this quantity exhibits the following power laws called the FamilyVicsek scaling [30]:
(5) 
with a scaling function ,
scaling exponents , and ,
and a crossover length .
It was indeed verified for the DSM2 growth
for both circular and flat interfaces
(6) 
with two constant parameters and , a rescaled random variable and a rescaled coordinate .
Now we are ready to probe finer statistical properties
of interface fluctuations, namely their distribution and correlation,
through the rescaled variable .
First, the skewness and the kurtosis were computed from the cumulants,
defined by and ,
respectively
(7) 
with and . These analytical results, together with the experimental support, imply that the TracyWidom distributions are universal characteristics of the dimensional KPZ class. In particular, the single KPZ class splits into at least two universality subclasses separating the circular (or curved) and flat interfaces.
Correlation can also be studied experimentally. For the spatial correlation, analytical studies have shown that it is described by time correlation of the particular stochastic processes called the Airy and Airy processes, for the curved and flat interfaces, respectively [35, 36, 37]. Specifically, the twopoint spatial correlation
(8) 
is given by
(9) 
with , the Airy process ( or ), and certain coefficient linked to parameters in the KPZ equation, given by in the case of isotropic growth [18]. This result for solvable models was also reproduced experimentally as shown in figure 4(a), again with finitetime corrections as expected. This actually implies that the difference between the circular and flat interfaces is not necessarily merely quantitative, because it is theoretically known that the Airy correlation for the circular interfaces decreases as , while the Airy correlation for the flat interfaces decays faster than exponentially [42].
In contrast to the spatial correlation, time correlation remains intractable by analytical studies, leaving room for unique contributions by experiments and other empirical or phenomenological approaches. The liquidcrystal experiment showed that it is also different between the two cases, qualitatively, as shown in figures 4(b,c) for the twotime correlation function
(10) 
For the flat interfaces (figure 4(b)), correlation decays toward zero, following the natural scaling ansatz that takes into account aging effect in the scaleinvariant growth, with (inset). By contrast, the results for the circular interfaces (figure 4(c)) suggest that the same ansatz does not work as well (as rescaled data do not perfectly overlap) and, moreover, the correlation decays as in the rescaled unit, hence remains strictly positive in the original scale. These claims are not definitive since was varied by less than a digit experimentally, but the same result was obtained numerically [43], and the latter claim is also supported by phenomenological theory due to Singha [44]. Similarly, the flat and circular interfaces were found to be different in many other quantities, notably in persistence properties, as summarized in table 1 of [18].
According to analytical developments [34, 35, 36, 37],
the curved (circular) and flat subclasses are not the only ones
constituting the KPZ class.
Another important subclass is that for the stationary interfaces,
which still grow but without changing their statistical properties.
To study this subclass theoretically, the initial condition is set
to be a stationary interface itself,
e.g., a flat interface in the limit ,
which is equivalent to the onedimensional Brownian motion
for the dimensional KPZ class
[7, 9, 10, 11].
Then the height difference from the initial time was analytically shown
to obey Baik and Rains’ distribution
[45, 34].
This subclass is, however, unable to be reached experimentally or numerically,
unless a stationary interface is artificially prepared
as an initial condition
3 Other experimental systems
While the previous section dealt entirely with the liquidcrystal experiment, this is by no means an only experimental system to study absorbingstate transitions or growing interfaces. On the contrary, universality implies, at least in principle, that the same macroscopic descriptions apply to a vast variety of phenomena and situations, as is indeed established on the numerical side. The present section overviews other experimental systems related to absorbingstate transitions and to growing interfaces, especially those concerning the DP class and the KPZ class.
3.1 Experiments on absorbingstate transitions
(1+1)D System  

annular RayleighBénard [12]  
annular RayleighBénard [13]  
linear RayleighBénard [13]  
depinning in paper wetting [47]  
viscous fingering [48]  
fluid vortices [49]  
TaylorDean [50]  
TaylorCouette [51]    
granular flow [52, 53]  
torsional Couette [54]  
ferrofluidic spikes [55]  
lateral heat convection [56]  
viscoelastic TaylorCouette [57]  
DP [58]  
(2+1)D System  
liquid columns [59]  
DSM1DSM2 [14, 15]  
DP [25, 26] 
Number in parentheses is the range of error given by the authors of each article.
The estimates for and were given for the two different axes in the plane, parallel (upper raw) and perpendicular (lower raw) to the initial molecular alignment.
Table 1 presents a list of experiments
– most probably not a complete list –
performed in situations where one might expect
the DPclass criticality to arise.
This shows that the values of experimentally measured critical exponents
are rather widely distributed instead of being universal,
as already discussed in section 1,
in a strong contrast with the status of numerical investigations
[5, 6].
One notices in particular that,
even if some of the exponents agree with the DPclass values,
others may be significantly different
(see, e.g., the values for the linear RayleighBénard convection
[13] or for the ferrofluidic spikes
[55] in table 1).
Therefore, for experimental investigations on absorbingstate transitions,
it is essential to determine a complete set
of independent critical exponents,
such as
or for the DP class
Spatiotemporal intermittency of ferrofluidic spikes.
One of the most thorough experiments
on onedimensional spatiotemporal intermittency was carried out
by Rupp et al[55] using magnetic fluids.
When a horizontal layer of magnetic fluid
is subjected to a strong enough vertical magnetic field,
the Rosensweig instability leads to formation of liquid spikes.
Using a sharp edge of a cylindrical electromagnet
to produce an inhomogeneous magnetic field,
Rupp et altrapped ferrofluidic spikes on this edge,
which were regularly aligned at a constant interval.
Then they added an alternating component to the magnetic field,
which induced spatiotemporal chaos or spatiotemporal intermittency
of the spikes, depending on the amplitude of the alternating component.
They analyzed the steady state of this spatiotemporal intermittency
for different amplitudes of the alternating component,
determining five critical exponents,
,
as listed in table 1.
Interestingly, while those for the order parameter ()
and for the spatial correlation ( and )
indicated the values for the dimensional DP class,
those for the temporal correlation ( and )
were found significantly different from the DP class
Onset of pipeflow turbulence.
The onset of a sustained turbulent state
in a flow along a long pipe is an outstanding issue in fluid mechanics,
which dates back to the seminal work by Reynolds [60]
in 1883.
Since the laminar flow in this setup, called the HagenPoiseuille flow,
is known to be linearly stable for all Reynolds numbers,
the problem consists in persistence of a turbulent state
induced by a finite disturbance to the system.
For relatively low Reynolds numbers ,
this turbulent state was shown to have a rather localized structure
along the stream, hence called the turbulent puff
[61, 62].
Recently, Hof and coworkers have performed a series of experiments
[62, 63, 64, 65]
using carefully designed pipes of unprecedented lengths
(typically and maximum ,
which amount to 3750 and 7500 times the pipe diameter, respectively).
They found that the survival probability of a turbulent puff decays
exponentially in time, suggesting a memoryless stochastic decay process,
with the mean life time that grows very rapidly with
but does not diverge near the apparent critical point
[63, 64].
At first glance, this result may imply
the absence of a welldefined critical Reynolds number.
However, their subsequent work [65] showed that
a puff can also split into two as a result of another stochastic event,
whose typical time grows rapidly with decreasing .
These results were further confirmed by another experimental group
[66]
as well as by numerical simulations
[63, 67, 65].
In view of all these results,
Hof’s group located the critical point at ,
at which the decaying rate and the splitting rate are balanced
[65]
Apart from determining the critical point, Hof et al’s scenario based on the decaying and splitting of turbulent puffs is reminiscent of Pomeau’s argument [27] and the contact process (figure 1(e)). This leads to the simplest scenario where the critical dynamics of the pipeflow turbulence is governed by the DP universality class, if one puts aside the rarity of the DPclass criticality in experiments (table 1). Therefore, a conclusion should be made after direct analysis of its critical behavior, but the time scale needed for a single decay or splitting event at the critical point is desperately long to do so by experiments or by direct numerical simulations. A way to circumvent this difficulty is studying whether common factors that affect the DP class are present or not in pipeflow turbulence. As an example of this direction of research, Hof’s group measured interval distributions between two turbulent puffs both experimentally [68] and numerically [69] and found no indication of longrange correlation, though they were unable to reach the steady state near the critical point, where one would expect powerlaw distribution of laminar intervals characterized by the exponent . It could also be informative to study quenching dynamics after a sudden decrease of the Reynolds number [68, 70] in the context of the critical quench experiment for absorbingstate transitions (see section 2.1). In a different line, Hof et al’s approach for measuring the life time and the splitting time seems to provide a useful way to revisit other hydrodynamic systems exhibiting spatiotemporal intermittency. Hof’s group indeed applied it to the plane Couette flow successfully [71]. This system is particularly interesting, because they found, numerically, powerlaw distributions of laminar intervals near the critical point [71].
Reversibleirreversible transitions in colloidal suspensions.
Another interesting and nontrivial example of absorbingstate transitions was found in transitions between reversible and irreversible motions of colloids under a periodic shear. Pine et al[72] studied colloidal motion in a viscous fluid confined in concentric cylinders (a Couette cell), the inner one being rotated back and forth slowly. In the absence of colloids, such a flow at low Reynolds numbers is described by the Stokes equation, and hence reversible with respect to a periodic displacement of the boundary. This is also true of the system mixed with nonBrownian colloids, i.e., particles large enough to neglect their Brownian motion, when they are shared periodically at small shear strain amplitudes; however, Pine et alfound that the particle motion becomes irreversible and effectively diffusive when the strain amplitude exceeds a certain threshold [72]. The origin of this irreversible motion was argued to be in chaos due to hydrodynamic interactions of nearby particles, but the connection between this microscopic chaos and the macroscopic irreversibility remained to be shown.
A reasonable explanation for this reversibleirreversible transition was later given by Corté et al[73] by means of a simple model (figure 5). Consider a twodimensional space and assume a finite interaction range of particles. These particles are distributed randomly at the initial time. When the system is sheared along the axis with the strain amplitude , each particle is displaced by (figure 5(b)). This may cause nearby particles to encounter each other, which results in random displacement of the particle positions after a cycle of the periodic shear, while those which did not experience collision return to their original positions (figure 5(c)). Then, for small , the particles eventually find such a configuration that no pair of particles collides any more, whereas such a state is never reached for large enough . This is the reversibleirreversible transition. It can be regarded as a transition to an absorbing state, since a collisionless configuration implies no further rearrangement of particles. Here, however, the transition is coupled with a conserved field because of the fixed particle number. Therefore, one expects another universality class to arise, called the conserved DP class or the Manna class [5, 6, 74], instead of the DP class. Indeed, Corté et al’s model as well as the colloidal experiment showed indications of an an absorbingstate transition, though they did not obtain a quantitative agreement with the conserved DP class [73] ( and from twodimensional simulations [73], to be compared with and for the conserved DP class [6, 74]; and for the experiment in three dimensions [73], to be compared with and for the conserved DP class [6, 74]). Moreover, Corté et alreported an interesting finding that a rheological quantity, namely the elastic component of the complex modulus, , behaves similarly to the order parameter [73], though the reason for this is yet to be clarified. After this work, Menon and Ramaswamy [75] introduced a lattice version of Corté et al’s model and found a quantitative agreement with the conserved DP class both for two and three dimensions. They also argued that longrange hydrodynamic interactions may explain the different critical exponents observed in the experiment.
Some interesting studies have been carried out subsequently. Franceschini et al[76] used rods instead of spherical particles, and revealed two different transitions as was increased, the first one for the formation of nematic alignment and the second one being the reversibleirreversible transition. Interestingly, in this setup, the latter transition was nicely accounted for by the conserved DP class, with the correct values of and for three dimensions [76]. From different perspectives, Mangan, Reichhardt, and Olson Reichhardt argued an analogous transition in driven systems with quenched disorder, such as vortices in superconductors (or other types of particles) under periodic driving [77] and plastic depinning of trapped particles [78], on the basis of their numerical simulations. Indeed, Okuma et al[79, 80] experimentally identified the predicted transitions using vortices driven by a Corbino disk made of an amorphous superconductor (a)MoGe. In their experiments, the reversibleirreversible transition under periodic driving gave and the depinning transition under constant driving yielded . Finally, it would be worth mentioning that the existence of the conserved DP class was recently challenged by numerical studies due to Basu et al[81, 82]. They claimed that common models for the conserved DP class actually showed the DPclass exponents, when they used a “natural” initial condition that takes into account correlated configuration of inactive particles in the stationary state, hence the previously reported exponents for the conservedDP class were numerical artifacts. This is still controversial as no analytic argument is given yet, especially on the different fieldtheoretic descriptions of the two classes [5, 6, 74, 83, 84]. In any case, though, it is true that the two universality classes have very similar sets of the critical exponents as well as the scaling functions [74], usually indistinguishable in experiments. In this context, studying effect of walls, originally proposed for numerical simulations [84], may also be useful for experiments.
3.2 Experiments on growing interfaces
Now we turn our attention to interface fluctuations due to local growth processes. Because of their ubiquity and importance, a great number of experiments have been carried out for a wide variety of systems, such as fluid flow in porous media, paper wetting, growing bacterial colony, to name but a few [7, 9, 10, 11, 85]. However, these studies have reported a wide range of the exponent values, mostly significantly different from those for the KPZ class; for onedimensional interfaces, the KPZclass values are and , while experiments tend to indicate larger values, typically and [7, 9, 10, 11] (see, e.g., table 11.1 in [7]). For some of the experiments, the deviation from the KPZ class can be understood as consequences of quenched disorder [7, 86, 47, 87, 88, 89, 90, 91, 92], longrange correlation [93], and/or powerlaw noise distribution [94, 95], which are common factors that intervene in typical experimental systems. In particular, imbibition in porous media has been studied in depth along these lines, for which the presence of such effects and consequences in scaling behavior are now rather clear in theory and experiments [85].
Among those factors affecting the scaling laws,
effect of quenched disorder has been particularly investigated
and seems to be relevant in many experiments.
First of all, in the presence of quenched noise,
interfaces may be pinned or unpinned, according to the driving strength.
While scaling behavior is expected to remain unchanged
when interfaces are driven strongly,
a few different universality classes may arise
near the depinning transition,
depending on the anisotropy and the nonlinearity
of the growth process at stake [89].
If the nonlinearity has its origin in kinematics
and hence vanishes at the transition,
as expected for simple isotropic growth processes,
the critical dynamics will be governed
by the quenched EdwardsWilkinson class
[89, 87, 90]
On the other hand, there have been reported a few experiments showing the KPZclass exponents, the number growing rather rapidly in recent years. These experiments are overviewed in the following subsections, apart from the liquidcrystal experiment already described in section 2.2. Note that there are also a few indirect realizations of the KPZclass scaling laws, specifically, in rupture lines [96, 97] and in crystal facets [98], which are not dealt with in this article.
Colony growth of bacteria and other cells.
Growth of bacterial colonies on agar is one of the systems well studied
in the context of pattern formation [99],
showing a variety of patterns like dendric and concentric ones.
One can also observe compact colonies for certain concentrations
of agar and nutrients, which can therefore be regarded as realizations
of growing interfaces.
Vicsek et al[100] studied such interfaces
using Escherichia coli and Bacillus subtilis
and reported and , respectively.
This latter result was later reproduced
by Wakita et al[101]
with a similar exponent value .
However, Wakita et alalso used a mutant strain of B. subtilis
that does not secrete surfactant,
in order to study another type of compact growing colonies
exhibiting less rough interfaces.
It is in this case that their data were found to be consistent
with the KPZclass exponents and
[101]
Recently, Huergo et alperformed analogous experiments for eukaryotic cells, namely Vero cells [103, 104] and cancer cells HeLa [105], which form colonies directly on a Petri dish. Then they found the two KPZclass exponents and for all the above cases, which may be surprising in view of the usually nonKPZ behavior of the bacterial colonies and of more complex characters of eukaryotic cells. Moreover, for Vero cells, Huergo et alrealized both flat and circular interfaces, which are particularly intriguing from the viewpoint of the recent developments on the geometrydependent universality in the onedimensional KPZ class (see section 2.2). Therefore, performing analyses developed for the liquidcrystal turbulence [18] would be an important attempt for testing the robustness of this particular universality in such very different experimental systems.
Slow combustion of paper.
Another interesting experimental example of the KPZ class was reported on slow, flameless combustion of paper, almost concurrently with Wakita et al’s work on bacteria. While an earlier experiment on slow combustion indicated [106], Timonen’s group found the KPZclass exponents by a series of highly controlled experiments [107, 108]: they prepared paper with uniform distribution of oxidizer KNO, and realized uniform heat transfer, by controlling air flow and by compensating heat loss at paper boundaries by means of heater filaments. In addition to the KPZclass exponents, they also determined values of the coefficients in the KPZ equation [108, 109] and amplitude ratios in the saturated regime [108], i.e., the longtime regime where the interface roughness does not grow any more because of a finite size. Then they found, in terms of the quantity called the universal coupling constant, for three different sets of samples [108]. This quantity is now understood to be , with the random variable of the BaikRains distribution for the stationary interfaces [34]. Their results therefore indicate that the saturated regime of finitesize systems may also be characterized by the BaikRains distribution, though this statement should be confirmed by measuring other cumulants from the experimental and numerical sides, as well as by analytical studies on finitesize systems (see footnote 4 on page 4).
In subsequent work, Timonen and coworkers measured firstpassage properties of interface fluctuations, namely the persistence probability [110], i.e., the probability that fluctuations remain positive (or negative) over given durations in time (or lengths in space). Their results for the saturated regime showed theoretically predicted power laws both in space and time. In contrast, for the growing regime (flat interfaces), their temporal persistence probability did not show asymmetry between positive and negative fluctuations [110], as opposed to earlier numerical work [111] and to the liquidcrystal experiment [18]. Similarly, the spatial persistence probability was found to decay algebraically for the papercombustion experiment, in contrast to an exponential decay found in the liquidcrystal experiment [18] which was later supported by an analytical study [112]. This apparent discrepancy may be due to a rather narrow time window for the growing regime in the papercombustion experiment, as the authors discussed [110], but this needs to be clarified by further investigations.
Finally, Timonen’s group also studied the onepoint distribution. They found that the experimentally obtained histograms for the growing and saturated regimes were reasonably fitted by the GOE TracyWidom and BaikRains distributions, respectively, by adjusting the mean and the variance of these theoretical distributions to the experimental data [113]. However, when the mean and the variance are adjusted as free parameters, or, equivalently, normalized to zero and unity, respectively, the difference between these distributions becomes barely visible only below (see figure 6 and discussions in section 4.2), which is obscured by statistical and experimental errors in the data presented in [113]. They also provided values of the skewness for the growing and saturated regimes, specifically 0.33 and 0.32, respectively, which however lay between the skewness of the GOE TracyWidom distribution, 0.2935, and that of the BaikRains distribution, 0.359, without favoring one of them. To avoid the loss of information by normalizing the distributions, one can measure the two constant parameters and in equation (6) and rescale the height variable by (section 4.2). This can be compared with the theoretical distributions without the need of ad hoc fitting, which are readily distinguishable in this case. Anyhow, the series of the papercombustion experiments by Timonen’s group provides a beautiful example where a variety of theoretical achievements on the KPZ class can be tested in depth in experiments.
Particle deposition on coffee ring.
Very recently, Yunker et al[114] found an experimental realization of the KPZclass scaling laws associated with the coffee ring effect. The coffee ring effect [115] refers to formation of a ring stain after evaporation of a suspension droplet – the phenomenon that one finds when a droplet of coffee falls onto a glass table and dries out. This pattern results from deposition of particles driven by capillary flow toward the edge of the droplet [115]. Yunker et al[114] studied this deposition process with elongated particles and found three distinct regimes of the kinetic roughening by varying the aspect ratio of the particles. For spheres, i.e., unit aspect ratio, they found , which is explained by a spatially uncorrelated Poissonlike deposition process, . For highly elongated particles with aspect ratios , they obtained and and compared with the quenched KPZ class. In between, for the particles with aspect ratios 1.1 to 1.2, they revealed and , in agreement with the KPZ class. Moreover, they also measured the skewness and the kurtosis as functions of time, which were found to converge to the values of the GUE TracyWidom distribution within the range of error. This confirms the characteristic distribution for the curvedinterface subclass, as we have seen in section 2.2. From a microscopic viewpoint, particles with larger aspect ratios deform the airwater interface more strongly, leading to longer ranges of interparticle attraction [114]. This gives an intuitive explanation on why the Poissonlike and the KPZclass growth processes arise for unit and intermediate aspect ratios, respectively. Apart from the scientific importance, this experiment tells us that the profound universal characteristics of the KPZ class may even underlie common phenomena that we experience in our daily lives.
Fronts of chemical waves in disordered media
Another very recent experiment showing the KPZclass exponents was reported
on fronts of chemical reactions in porous media.
Atis et al[116, 117] carried out experiments
on the iodate arsenous acid autocatalytic reaction,
,
between two vertical acrylic glass plates,
placed parallel to each other
and filled with bidisperse glass beads in between.
The cell was filled with reactant solution
and the chemical reaction was initiated from the bottom end,
by immersing it in a reservoir of reagent containing autocatalyst.
Without external flow, the reaction forms a smooth horizontal front,
propagating upward with a stationary concentration profile.
However, in the presence of flow
introduced by injection or suction of the solution,
the smooth front profile is destabilized
by disordered flow in the porous medium, forming selfaffine interfaces.
Atis et althen found four different regimes
by varying the flow rate:
(1) supportive (upward) flow, downstream (upward) front propagation (SD),
(2) weak adverse (downward) flow, upstream (upward) front propagation (AU),
(3) moderate adverse (downward) flow, eventually static (frozen) front (AS),
and (4) strong adverse (downward) flow,
downstream (downward) front propagation (AD)
[116], which were also reproduced
by their lattice Boltzmann simulations [118].
They experimentally measured the two scaling exponents and
for the four regimes and found that those for the SD and AD regimes
are consistent with the KPZclass exponents and ,
while for the AS regime
4 Remarks on analyses of scaling laws in experiments
As already stressed, absorbingstate transitions and growing interfaces have been extensively studied by numerical simulations, but not by as many experiments, especially if it concerns universality. As a result, analyses used by numerical studies are often employed for experimental work. Such methods are of course valid in experiments, but not necessarily optimal, because experiments usually have stronger restrictions on the measurable quantities and on the number of available independent experimental runs. In addition, secondary or uncontrolled factors may intervene in the phenomenon of interest. In view of this situation, the present section describes some general and practical remarks on analyses of the scaling laws for absorbingstate transitions and for growing interfaces, focusing in particular on how to test the critical behavior of the DP/conservedDP class for the former, and the universal distributions of the dimensional KPZ class for the latter.
4.1 Experimental tests of the DP/conservedDP class
Numerical studies over many years have established
the following three main approaches to characterize critical behavior
of absorbingstate transitions:
steadystate measurement, critical quenching, and critical spreading,
as we have already seen briefly in section 2.1
(see reviews [5, 6]
for detailed protocols).
For the critical quenching,
one measures relaxation of the order parameter
from a fully active initial condition,
which can typically be achieved by a sudden change
in the value of the control parameter
(the applied voltage for the electroconvection).
For the critical spreading,
one tracks the fate of a single active patch
and measures the probability
that the resulting active cluster survives until time ,
together with its volume and squared radius.
First of all, as remarked in section 3.1,
the DPclass criticality should be tested at least by three independent
critical exponents, hopefully more.
To this end, it is always worth trying all the above three approaches,
whenever possible.
Data collapse as in figure 2(b) is also helpful;
comparing the constructed scaling function
with numerical data for the DP class
That having said, those three protocols have strengths and weaknesses. Numerically, the critical spreading protocol is widely regarded as the best method, since one only needs to evolve a few number of active sites for most of runs and the measurement is completely free from finitesize effects [5, 6]. However, it is far less useful in experiments, because one does not benefit from the former point and will be faced with poor statistics for the survival probability. One also needs to develop a method to introduce a single active patch into the given experimental system. Concerning the steadystate measurement, while most experimental studies have relied on it, it is known to be problematic, mainly because, for finitesize systems, the system eventually experiences an unfortunate moment at which all active patches disappear at the same time. In this sense, the system always ends up with an absorbing state, even if it is apparently in the active phase. Of course, this does not occur in practice if the system size is large enough, and if the system is not too close to the critical point. Conversely, as far as an absorbingstate transition is concerned, one is expected to encounter such an event near the critical point; otherwise spontaneous nucleation of active patches may be present at a nonnegligible rate, wiping out the critical behavior expected for absorbingstate transitions. In contrast, the critical quench protocol seems to be possible and useful in many experiments. This allows us to compare the scaling function in addition to the critical exponents, as exemplified in figure 2(b). Finally, for some experiments, one is obliged to probe critical behavior through only a few observables, which are not necessarily appropriate quantities to characterize critical phenomena. In this case, one may use hysteresis to extract an estimate of the exponent [119]: if spontaneous nucleation rate is tiny but nonzero, which is a realistic assumption for experiments, ramping the control parameter up and down through the critical point produces hysteresis loops, whose average size scales as with the ramping rate . Because one can choose any observable to measure hysteresis associated with the critical behavior (as far as such hysteresis is visible), this method may be useful in many experimental systems.
4.2 Experimental tests of the KPZclass universal fluctuations
Concerning the scaleinvariant fluctuations of growing interfaces, although the methods to measure the scaling exponents through the FamilyVicsek scaling (5) are now quite standard and have been applied to many experiments [7, 9, 10, 11], experimental results are often more complicated, because of, among others, different kinds of crossover effects. First, the FamilyVicsek scaling of course arises for lengths and times much larger than characteristic scales of microscopic interactions. Even if this condition is met and the asymptotic growth is indeed governed by the KPZ class, the nonlinear term is negligible at the early times because the interfaces are not yet sufficiently rough; therefore one expects crossover from the linear growth regime, called the EdwardsWilkinson regime ( and in one dimension), to the KPZ regime () [11]. This can take place quite late, if nonlinearity in the growth process is weak. To add, in the presence of quenched disorder, as is often the case in experiments, one has additional crossover effects as discussed in the beginning of section 3.2. The presence of a wall may also introduce another crossover effect [120]. One needs to overcome all these crossover effects to infer asymptotic scaling behavior, which is difficult in many experiments because of the rather limited ranges of accessible scales.
Despite this difficulty, we have seen in section 3.2 that some recent experiments have clearly identified the scaling exponents of the dimensional KPZ class. Then, in view of the recent analytical developments discussed in section 2.2, a central question to ask is whether the universality beyond the scaling exponents, at the level of the distribution and correlation functions, is still valid for those experimental systems. Since only few experiments have tested this universality yet, here we discuss a general procedure for this, assuming that the KPZclass exponents are already identified in a given experimental (or numerical) system.
First of all, one should be aware that the theoretical distribution functions discussed in this context, namely the TracyWidom distributions and the BaikRains distribution, have only quantitative differences in their functional forms. While they can readily be distinguished in their original definition of the stochastic variable (figure 6(a)), the difference is quite subtle when they are arbitrarily shifted and scaled (figure 6(b)). When they are normalized to have zero mean and unit variance, for example, the difference is barely visible only beneath approximately in the probability density, which is usually obscured by statistical and experimental errors in most experiments. Therefore, it is essential to rescale the height variable appropriately according to equation (6), by using the two constant parameters and , instead of fitting the theoretical distribution functions by adjusting the mean and the variance as free parameters.
The two parameters and can be determined as follows:

Estimate the asymptotic growth rate by measuring as a function of time. Since with a constant , one may plot against and make a linear extrapolation to read the intercept.

Determine the amplitude of the interface fluctuations. This can be obtained by measuring the secondorder cumulant . Note that one needs to set here the value of . This can be chosen arbitrary, e.g., , but one may choose the variance of the compared theoretical distribution to facilitate the comparison. This choice does not bias the results [17]. Alternatively, Alves et al[121] proposed estimating from the firstorder cumulant by , on the basis of their numerical simulations. While they found it more accurate than that measured from the secondorder cumulant for some of the models they studied [121], this method did not work better for the liquidcrystal experiment. The precision of the two estimates is determined by finitetime corrections in the corresponding cumulants, which are systemdependent, and so is the accuracy of the two methods. Note also that the parameter may be estimated from coefficients of the KPZ equation [16, 18, 121], but recent numerical and theoretical work showed that this relation may be changed in the presence of correlated noise (even shortranged one) [122], which is unavoidable in experimental systems. Therefore, it would be safer to estimate directly from the cumulants of height fluctuations, as explained here.
Once the two parameters and are determined, one can make a histogram of the rescaled height to infer the distribution of the random variable . This can be directly compared with probability density functions as plotted in figure 6(a). In addition, it is always worth plotting time series of the differences in the thorder cumulants, . This not only provides a quantitative measure of the agreement with the theoretical distribution, but also yields finitetime corrections to the asymptotic distribution, whose universal aspect is not fully understood yet [18, 121, 123]. Finally, using the rescaled variables, one can investigate many other statistical properties related to spacetime correlation. Experimental (or numerical) investigations of time correlation are particularly important, as it still remains unsolved by analytical means.
5 Concluding remarks
Physics of scaleinvariant phenomena driven far from equilibrium needs further investigations in many aspects, especially when one compares with their equilibrium counterparts. As one cannot rely on the celebrated equilibrium statistical mechanics, connection between microscopic evolution and resulting macroscopic dynamics is a very delicate problem for outofequilibrium systems. This makes it much more difficult to discriminate between relevant and irrelevant ingredients that determine the universality class. In this respect, experimental investigations play a unique role in testing the robustness of the outofequilibrium universality under real situations, in the presence of nontrivial interactions, quenched disorder (more or less), selfgenerated outofequilibrium noise, and so forth. Experiments thereby complement numerical studies, which can control ingredients of the system as desired, albeit sometimes oversimplified compared with real systems.
As overviewed in this article, recent years have marked considerable progress in experimental studies on outofequilibrium scaling laws for absorbingstate transitions and growing interfaces, with a growing number of experiments showing direct link or evidence for the DP class and the KPZ class. Having more experimental examples will not only reinforce powerful universality of the outofequilibrium scaling laws; by comparing with other experimental systems, one may elucidate what are the relevant parameters that separate experimental systems following and not following the simple universality scenario in practical situations. Let us also recall that both the DP class and the KPZ class correspond to the most fundamental situation without any additional symmetry or conservation laws. Firm experimental grounds for these two classes would therefore help exploring other universality classes in the presence of extra symmetry or conservation in experimental systems, as partly discussed on the conserved DP class in section 3.1. Finally, as stressed in this article, the dimensional KPZ class provides a unique situation for studies of outofequilibrium scaling laws, which allows direct comparison of exact theoretical solutions and quantitative experimental results, at the level of detailed statistical properties such as the distribution and correlation functions. While theoretical studies can deal with situations hard to realize experimentally, experimentalists can measure statistical properties that are analytically intractable. Of course numerical studies are also indispensable in many of these approaches. The author believes that the current situation is the beginning of such cooperative developments among theory, simulations and experiments, which will deepen our understanding on outofequilibrium scaling laws and possibly provide novel perspectives on them.
References
Footnotes
 Note that the KPZ problem can be mapped to an equilibrium model under disordered environment, namely the directed polymer in random medium [7, 9, 10, 11]. Moreover, in one dimension, the steady state of the KPZ class satisfies a certain form of the fluctuationdissipation theorem [8, 7, 9, 10, 11]. These facts, however, do not indicate that the KPZ class describes equilibrium systems in the usual way universality classes do, so are not contradictory to the outofequilibrium nature of the KPZ class or corresponding growth processes.
 The critical voltage in the quench experiment was determined from figure 2(b), independently of the estimate from the steadystate experiment. This is to avoid influence of a slight, uncontrolled shift in the experimental conditions during the days that separated the two sets of experiments [15].
 However, for circular interfaces, in general, spatial average in the definition of results in a slight bias in the estimate of from the FamilyVicsek scaling [18]. One should use instead the standard deviation defined by ensemble average, .
 Definitions of the cumulants are given, up to the fourth order, by , , and , with ensemble average and .
 The GUE and GOE random matrices and the corresponding TracyWidom distributions are defined as follows. First, consider an complex Hermitian matrix (hence ) with random matrix elements. Specifically, the real and imaginary parts are independently drawn from the Gaussian distribution, satisfying for all , , and for . This ensemble of random matrices is called GUE. Similarly, GOE is defined by replacing complex Hermitian matrices with real symmetric matrices. Then, whether GUE or GOE, the matrix has real random eigenvalues, the largest of which scales as for large , in the normalization adopted here. The distribution of this random variable in the limit is called the TracyWidom distribution. Analytic expressions are available for GOE, GUE, and their symplectic counterpart [32, 33]. In passing, it is suggestive to compare the above expression for and equation (6) for the growing interfaces.
 Note however that, for the sake of simplicity, in this article is defined by multiplying to the standard definition in random matrix theory, as suggested by analytical results for solvable growth models [34, 35, 36, 37].
 See [35, 36, 37] for reviews on these analytical developments.
 It is not clear, at least strictly, how the properties of the stationary subclass appear in the steady state reached within a finitesize system, since one cannot avoid influence from the boundary in the asymptotic limit. As a matter of fact, no exact solution is obtained yet for systems with finite sizes.
 Although one needs in general four independent exponents to characterize absorbingstate transitions, for the DP class the number reduces to three thanks to the rapidity symmetry , as briefly discussed in section 2.1.
 The authors of [55] cited an erroneous value of for the DP class, which led them to conclude that four of the five exponents agreed with the DP class.
 Strictly, the critical point is not exactly given by the balance of the two local transition rates. However, for the pipeflow turbulence, both the life time and the splitting time grow superexponentially fast near the critical point [64, 67, 65], so that the difference from the true critical point would be smaller than the given range of error [65].
 The EdwardsWilkinson equation refers to the linear continuum equation obtained by setting in equation (7). If one replaces the usual noise term by quenched one, , one obtaines the quenched EdwardsWilkinson equation. One can similarly define the quenched KPZ equation.
 Wakita et aldid not provide an estimate of , but they presented data for the roughness growth [101], which are consistent with albeit the narrow range of the scaling region.
 In the AS regime, front profiles are eventually frozen, forming a sawtoothlike pattern [116]. At a given time before the complete formation of the final pattern, some parts of the front are already frozen, while others are still evolving. The scaling exponents were measured only from these evolving parts of the front [117].
 Numerical data for some of the (2+1)dimensional DPclass scaling functions are available from the author upon request; namely, those for the order parameter relaxation in critical quenching (figure 2(b) inset), as well as those for the survival probability, the volume expansion, and the squared radius of the active cluster in critical spreading.
References
 Goldenfeld N 1992 Lectures on Phase Transitions and the Renormalization Group Frontiers in Physics (Colorado: Westview Press)
 Henkel M 1999 Conformal Invariance and Critical Phenomena (Berlin: Springer)
 Wilson K G 1975 Rev. Mod. Phys. 47 773–840
 Täuber U C (to be published) Critical Dynamics: A Field Theory Approach to Equilibrium and NonEquilibrium Scaling Behavior (Cambridge: Cambridge Univ. Press)
 Hinrichsen H 2000 Adv. Phys. 49 815–958
 Henkel M, Hinrichsen H and Lübeck S 2009 NonEquilibrium Phase Transitions: Volume 1: Absorbing Phase Transitions 1st ed Theoretical and Mathematical Physics (Dordrecht: Springer)
 Barabási A L and Stanley H E 1995 Fractal Concepts in Surface Growth (Cambridge: Cambridge Univ. Press)
 Kardar M, Parisi G and Zhang Y C 1986 Phys. Rev. Lett. 56 889–892
 Meakin P 1993 Phys. Rep. 235 189–289
 HalpinHealy T and Zhang Y C 1995 Phys. Rep. 254 215–414
 Krug J 1997 Adv. Phys. 46 139–282
 Ciliberto S and Bigazzi P 1988 Phys. Rev. Lett. 60 286–289
 Daviaud F, Bonetti M and Dubois M 1990 Phys. Rev. A 42 3388–3399
 Takeuchi K A, Kuroda M, Chaté H and Sano M 2007 Phys. Rev. Lett. 99 234503
 Takeuchi K A, Kuroda M, Chaté H and Sano M 2009 Phys. Rev. E 80 051116
 Takeuchi K A and Sano M 2010 Phys. Rev. Lett. 104 230601
 Takeuchi K A, Sano M, Sasamoto T and Spohn H 2011 Sci. Rep. 1 34
 Takeuchi K A and Sano M 2012 J. Stat. Phys. 147 853–890
 de Gennes P G and Prost J 1995 The Physics of Liquid Crystals 2nd ed (International Series of Monographs on Physics vol 83) (New York: Oxford Univ. Press)
 Kai S and Zimmermann W 1989 Prog. Theor. Phys. Suppl. 99 458–492
 Frisch U 1995 Turbulence: The Legacy of A. N. Kolmogorov (Cambridge: Cambridge Univ. Press)
 Kai S, Zimmermann W, Andoh M and Chizumi N 1989 J. Phys. Soc. Jpn. 58 3449–3452
 Kai S, Zimmermann W, Andoh M and Chizumi N 1990 Phys. Rev. Lett. 64 1111–1114
 Fazio V S U and Komitov L 1999 Europhys. Lett. 46 38–42
 Grassberger P and Zhang Y C 1996 Physica A 224 169–179
 Voigt C A and Ziff R M 1997 Phys. Rev. E 56 R6241–R6244
 Pomeau Y 1986 Physica D 23 3–11
 Janssen H K 1981 Z. Phys. B 42 151–154
 Grassberger P 1982 Z. Phys. B 47 365–374
 Family F and Vicsek T 1985 J. Phys. A 18 L75–L81
 Mehta M L 2004 Random Matrices 3rd ed (Pure and Applied Mathematics vol 142) (San Diego: Elsevier)
 Tracy C A and Widom H 1994 Commun. Math. Phys. 159 151–174
 Tracy C A and Widom H 1996 Commun. Math. Phys. 177 727–754
 Prähofer M and Spohn H 2000 Phys. Rev. Lett. 84 4882–4885
 Kriecherbauer T and Krug J 2010 J. Phys. A 43 403001
 Sasamoto T and Spohn H 2010 J. Stat. Mech. 2010 P11013
 Corwin I 2012 Random Matrices Theory Appl. 1 1130001
 Johansson K 2000 Commun. Math. Phys. 209 437–476
 Sasamoto T and Spohn H 2010 Phys. Rev. Lett. 104 230602
 Amir G, Corwin I and Quastel J 2011 Commun. Pure Appl. Math. 64 466–537
 Calabrese P and Le Doussal P 2011 Phys. Rev. Lett. 106 250603
 Bornemann F, Ferrari P L and Prähofer M 2008 J. Stat. Phys. 133 405–415
 Takeuchi K A 2012 J. Stat. Mech. 2012 P05007
 Singha S B 2005 J. Stat. Mech. 2005 P08006
 Baik J and Rains E M 2000 J. Stat. Phys. 100 523–541
 Takeuchi K A 2013 Phys. Rev. Lett. 110 210604
 Buldyrev S V, Barabási A L, Caserta F, Havlin S, Stanley H E and Vicsek T 1992 Phys. Rev. A 45 R8313–R8316
 Michalland S, Rabaud M and Couder Y 1993 Europhys. Lett. 22 17–22
 Willaime H, Cardoso O and Tabeling P 1993 Phys. Rev. E 48 288–295
 Degen M M, Mutabazi I and David Andereck C 1996 Phys. Rev. E 53 3495–3504
 Colovas P W and Andereck C D 1997 Phys. Rev. E 55 2736–2741
 Daerr A and Douady S 1999 Nature 399 241–243
 Hinrichsen H, JiménezDalmaroni A, Rozov Y and Domany E 1999 Phys. Rev. Lett. 83 4999–5002
 Cros A and Le Gal P 2002 Phys. Fluids 14 3755–3765
 Rupp P, Richter R and Rehberg I 2003 Phys. Rev. E 67 036209
 Lepiller V, Prigent A, Dumouchel F and Mutabazi I 2007 Phys. Fluids 19 054101
 Latrache N, Crumeyrolle O and Mutabazi I 2012 Phys. Rev. E 86 056305
 Jensen I 1999 J. Phys. A 32 5233–5249
 Pirat C, Naso A, Meunier J L, Maïssa P and Mathis C 2005 Phys. Rev. Lett. 94 134502
 Reynolds O 1883 Phil. Trans. R. Soc. London 174 935–982
 Wygnanski I J and Champagne F H 1973 J. Fluid. Mech. 59 281–335
 Hof B, van Doorne C W H, Westerweel J, Nieuwstadt F T M, Faisst H, Eckhardt B, Wedin H, Kerswell R R and Waleffe F 2004 Science 305 1594–1598
 Hof B, Westerweel J, Schneider T M and Eckhardt B 2006 Nature 443 59–62
 Hof B, de Lozar A, Kuik D J and Westerweel J 2008 Phys. Rev. Lett. 101 214501
 Avila K, Moxey D, de Lozar A, Avila M, Barkley D and Hof B 2011 Science 333 192–196
 Kuik D J, Poelma C and Westerweel J 2010 J. Fluid. Mech. 645 529–539
 Avila M, Willis A P and Hof B 2010 J. Fluid. Mech. 646 127–136
 Samanta D, Lozar A D and Hof B 2011 J. Fluid. Mech. 681 193–204
 Avila M and Hof B 2013 Phys. Rev. E 87 063012
 Moxey D and Barkley D 2010 Proc. Natl. Acad. Sci. USA 107 8091–8096
 Shi L, Avila M and Hof B 2013 Phys. Rev. Lett. 110 204502
 Pine D J, Gollub J P, Brady J F and Leshansky A M 2005 Nature 438 997–1000
 Corté L, Chaikin P M, Gollub J P and Pine D J 2008 Nat. Phys. 4 420–424
 Lübeck S 2004 Int. J. Mod. Phys. B 18 3977–4118
 Menon G I and Ramaswamy S 2009 Phys. Rev. E 79 061108
 Franceschini A, Filippidi E, Guazzelli E and Pine D J 2011 Phys. Rev. Lett. 107 250603
 Mangan N, Reichhardt C and Reichhardt C J O 2008 Phys. Rev. Lett. 100 187002
 Reichhardt C and Reichhardt C J O 2009 Phys. Rev. Lett. 103 168301
 Okuma S, Tsugawa Y and Motohashi A 2011 Phys. Rev. B 83 012503
 Okuma S, Kawamura Y and Tsugawa Y 2012 J. Phys. Soc. Jpn. 81 114718
 Basu M, Basu U, Bondyopadhyay S, Mohanty P K and Hinrichsen H 2012 Phys. Rev. Lett. 109 015702
 Basu U, Basu M and Mohanty P K 2013 Eur. Phys. J. B 86 1–7
 Vespignani A, Dickman R, Muñoz M A and Zapperi S 1998 Phys. Rev. Lett. 81 5676–5679
 Bonachela J A and Muñoz M A 2007 Physica A 384 89–93
 Alava M, Dubé M and Rost M 2004 Adv. Phys. 53 83–175
 Tang L H and Leschhorn H 1992 Phys. Rev. A 45 R8309–R8312
 Amaral L A N, Barabási A L and Stanley H E 1994 Phys. Rev. Lett. 73 62–65
 Amaral L A N, Barabási A L, Buldyrev S V, Harrington S T, Havlin S, SadrLahijany R and Stanley H E 1995 Phys. Rev. E 51 4655–4673
 Tang L H, Kardar M and Dhar D 1995 Phys. Rev. Lett. 74 920–923
 Leschhorn H, Nattermann T, Stepanow S and Tang L H 1997 Ann. Phys. (Berlin) 509 1–34
 Csahók Z, Honda K and Vicsek T 1993 J. Phys. A 26 L171–L178
 Csahók Z, Honda K, Somfai E, Vicsek M and Vicsek T 1993 Physica A 200 136–154
 Medina E, Hwa T, Kardar M and Zhang Y C 1989 Phys. Rev. A 39 3053–3075
 Zhang Y C 1990 J. Phys. (Paris) 51 2129–2134
 Horváth V K, Family F and Vicsek T 1991 Phys. Rev. Lett. 67 3207–3210
 Kertész J, Horváth V K and Weber F 1993 Fractals 1 67–74
 Engøy T, Måløy K J, Hansen A and Roux S 1994 Phys. Rev. Lett. 73 834–837
 Degawa M, Stasevich T J, Cullen W G, Pimpinelli A, Einstein T L and Williams E D 2006 Phys. Rev. Lett. 97 080601
 Matsushita M, Hiramatsu F, Kobayashi N, Ozawa T, Yamazaki Y and Matsuyama T 2004 Biofilms 1 305–317
 Vicsek T, Cserző M and Horváth V K 1990 Physica A 167 315–321
 Wakita J i, Itoh H, Matsuyama T and Matsushita M 1997 J. Phys. Soc. Jpn. 66 67–72
 Hallatschek O, Hersen P, Ramanathan S and Nelson D R 2007 Proc. Natl. Acad. Sci. USA 104 19926–19930
 Huergo M A C, Pasquale M A, Bolzán A E, Arvia A J and González P H 2010 Phys. Rev. E 82 031903
 Huergo M A C, Pasquale M A, González P H, Bolzán A E and Arvia A J 2011 Phys. Rev. E 84 021917
 Huergo M A C, Pasquale M A, González P H, Bolzán A E and Arvia A J 2012 Phys. Rev. E 85 011918
 Zhang J, Zhang Y C, Alstrøm P and Levinsen M T 1992 Physica A 189 383–389
 Maunuksela J, Myllys M, Kähkönen O P, Timonen J, Provatas N, Alava M J and AlaNissila T 1997 Phys. Rev. Lett. 79 1515–1518
 Myllys M, Maunuksela J, Alava M, AlaNissila T, Merikoski J and Timonen J 2001 Phys. Rev. E 64 036101
 Maunuksela J, Myllys M, Merikoski J, Timonen J, Kärkkäinen T, Welling M and Wijngaarden R 2003 Eur. Phys. J. B 33 193–202
 Merikoski J, Maunuksela J, Myllys M, Timonen J and Alava M J 2003 Phys. Rev. Lett. 90 024501
 Kallabis H and Krug J 1999 Europhys. Lett. 45 20–25
 Ferrari P L and Frings R 2013 J. Stat. Mech. 2013 P02001
 Miettinen L, Myllys M, Merikoski J and Timonen J 2005 Eur. Phys. J. B 46 55–60
 Yunker P J, Lohr M A, Still T, Borodin A, Durian D J and Yodh A G 2013 Phys. Rev. Lett. 110 035501
 Deegan R D, Bakajin O, Dupont T F, Huber G, Nagel S R and Witten T A 1997 Nature 389 827–829
 Atis S, Saha S, Auradou H, Salin D and Talon L 2013 Phys. Rev. Lett. 110 148301
 Atis S, Awadhesh K, Talon L and Salin D (in preparation)
 Saha S, Atis S, Salin D and Talon L 2013 Europhys. Lett. 101 38003
 Takeuchi K A 2008 Phys. Rev. E 77 030103(R)
 Allegra N, Fortin J Y and Henkel M 2013 arXiv 1309.1634
 Alves S G, Oliveira T J and Ferreira S C 2013 J. Stat. Mech. 2013 P05007
 Agoritsas E, Bustingorry S, Lecomte V, Schehr G and Giamarchi T 2012 Phys. Rev. E 86 031144
 Ferrari P L and Frings R 2011 J. Stat. Phys. 144 1123–1150
 Bornemann F 2010 Math. Comput. 79 871–915