Conserved dynamics and interface roughening in spontaneous imbibition: A phase field model

The propagation and roughening of a fluid-gas interface through a disordered medium in the case of capillary driven spontaneous imbibition is considered. The system is described by a conserved (model B) phase-field model, with the structure of the disordered medium appearing as a quenched random field . The flow of liquid into the medium is obtained by imposing a non-equilibrium boundary condition on the chemical potential, which reproduces Washburn's equation for the slowing down motion of the average interface position H. The interface is found to be superrough, with global roughness exponent , indicating anomalous scaling. The spatial extent of the roughness is determined by a length scale arising from the conservation law. The interface advances by avalanche motion, which causes temporal multiscaling and qualitatively reproduces the experimental results of Horv'ath and Stanley (Phys. Rev. E 52, 5166 (1995)) on the temporal scaling of the interface. Received 24 November 1999     

Front propagation in chaotic and noisy reaction-diffusion systems: a discrete-time map approach

We study the front propagation in reaction-diffusion systems whose reaction dynamics exhibits an unstable fixed point and chaotic or noisy behaviour. We have examined the influence of chaos and noise on the front propagation speed and on the wandering of the front around its average position. Assuming that the reaction term acts periodically in an impulsive way, the dynamical evolution of the system can be written as the convolution between a spatial propagator and a discrete-time map acting locally. This approach allows us to perform accurate numerical analysis. They reveal that in the pulled regime the front speed is basically determined by the shape of the map around the unstable fixed point, while its chaotic or noisy features play a marginal role. In contrast, in the pushed regime the presence of chaos or noise is more relevant. In particular the front speed decreases when the degree of chaoticity is increased, but it is not straightforward to derive a direct connection between the chaotic properties (e.g. the Lyapunov exponent) and the behaviour of the front. As for the fluctuations of the front position, we observe for the noisy maps that the associated mean square displacement grows in time as t ^1/2 in the pushed case and as t ^1/4 in the pulled one, in agreement with recent findings obtained for continuous models with multiplicative noise. Moreover we show that the same quantity saturates when a chaotic deterministic dynamics is considered for both pushed and pulled regimes. Received 17 July 2001     

Front Propagation Techniques to Calculate the Largest Lyapunov Exponent of Dilute Hard Disk Gases

A kinetic approach is adopted to describe the exponential growth of a small deviation of the initial phase space point, measured by the largest Lyapunov exponent, for a dilute system of hard disks, both in equilibrium and in a uniform shear flow. We derive a generalized Boltzmann equation for an extended one-particle distribution that includes deviations from the reference phase space point. The equation is valid for very low densities n, and requires an unusual expansion in powers of 1/|ln n|. It reproduces and extends results from the earlier, more heuristic clock model and may be interpreted as describing a front propagating into an unstable state. The asymptotic speed of propagation of the front is proportional to the largest Lyapunov exponent of the system. Its value may be found by applying the standard front speed selection mechanism for pulled fronts to the case at hand. For the equilibrium case, an explicit expression for the largest Lyapunov exponent is given and for sheared systems we give explicit expressions that may be evaluated numerically to obtain the shear rate dependence of the largest Lyapunov exponent.     

Suppression of antiferromagnetic correlations by quenched dipole-type impurities

The effects of quenched dipole moments on a two-dimensional Heisenberg antiferromagnet are found exactly, by applying the renormalization group to the appropriate classical non-linear sigma model. Such dipole moments represent random fields with power law correlations. At low temperatures, they also represent the long range effects of quenched random strong ferromagnetic bonds on the antiferromagnetic correlation length, , of a two-dimensional Heisenberg antiferromagnet. It is found that the antiferromagnetic long range order is destroyed for any non-zero concentration, x, of the dipolar defects, even at zero temperature. Below a line , where T is the temperature, is independent of T, and decreases exponentially with x. At higher temperatures, it decays exponentially with , with an effective stiffness constant , which decreases with increasing x/T. The latter behavior is the same as for annealed dipole moments, and we use our quenched results to interpolate between the two types of averaging for the problem of ferromagnetic bonds in an antiferromagnet. The results are used to estimate the three-dimensional Néel temperature of a lamellar system with weakly coupled planes, which decays linearly with x at small concentrations, and drops precipitously at a critical concentration. These predictions are shown to reproduce successfully several of the prominent features of experiments on slightly doped copper oxides. Received 22 October 1998     

Growth velocity and the topography of Ni-Zn binary alloy electrodeposits

We show that the electrodeposition of Ni-Zn alloys at the lowest growth velocities, v < 0.5μm/s, exclusively proceeds from an abnormal co-deposition phenomenon. The growth process in this v region greatly depends on the initial [Co²⁺] concentration of the film deposition bath. A theoretical approach of this process including the role of the saturation surface roughness of the alloy, , leads to an estimation of the transport properties of the ad-atoms involved during the deposit formation. Their surface diffusion coefficient varying between 1.76×10^-10 and 2.40×10^-8 cm^-2/s exhibits a minimal value, D _s = 2.10×10^-10 cm^-2/s located between v = 0.17 and 0.35μm/s. The spatial scaling analysis of the local roughness, σ, examined according to the power-law σ≈L ^α reveals that the resulting roughness exponents concurs with the Kardar-Parisi-Zhang dynamics including the restricted surface diffusion. Two main v regions leads to different fractal textural features of the alloy film surface. Below 0.10 μm/s, the roughness exponent obtained is α≈ 0.6, depicting a limited ad-atom mobility. Over v = 0.30μm/s, this exponent stabilises at α≈ 0.82, indicating an increase of the surface diffusion. Received 16 August 2000 and Received in final form 20 June 2001     

Universality classes of self-avoiding fixed-connectivity membranes

We present an analysis of extensive large-scale Monte Carlo simulations of self-avoiding fixed-connectivity membranes for sizes (number of faces) ranging from 512 to 17672 (triangular) plaquettes. Self-avoidance is implemented via impenetrable plaquettes. We simulate the impenetrable plaquette model in both three and four bulk dimensions. In both cases we find the membrane to be flat for all temperatures: the size exponent in three dimensions is ν = 0.95(5) (Hausdorff dimension d _H = 2.1(1)). The single flat phase appears, furthermore, to be equivalent to the large bending rigidity phase of non-self-avoiding fixed-connectivity membranes --the roughness exponent in three dimensions is ξ = 0.63(4). This suggests that there is a unique universality class for flat fixed-connectivity membranes without attractive interactions. Finally, we address some theoretical and experimental implications of our work. Received 23 June 2000 and Received in final form 25 October 2000     

Local waiting time fluctuations along a randomly pinned crack front

The propagation of an interfacial crack along a heterogeneous weak plane of a transparent Plexiglas block is followed using a high resolution fast camera. We show that the fracture front dynamics is governed by local and irregular avalanches with very large size and velocity fluctuations. We characterize the intermittent dynamics observed, i.e., the local pinnings and depinnings of the crack front by measuring the local waiting time fluctuations along the crack front during its propagation. The deduced local front line velocity distribution exhibits a power law behavior, P(v) alpha v-eta with eta=2.55+/-0.15, for velocities v larger than the average front speed . The burst size distribution is also a power law, P(S) alpha S-gamma with gamma=1.7+/-0.1. Above a characteristic length scale of disorder Ld approximately 15 microm, the avalanche clusters become anisotropic providing an estimate of the roughness exponent of the crack front line, H=0.66.     

Numerical study of local and global persistence in directed percolation

The local persistence probability P _l (t) that a site never becomes active up to time t, and the global persistence probability P _g (t) that the deviation of the global density from its mean value does not change its sign up to time t are studied in a (1+1)-dimensional directed percolation process by Monte-Carlo simulations. At criticality, starting from random initial conditions, P _l (t) decays algebraically with the exponent . The value is found to be independent of the initial density and the microscopic details of the dynamics, suggesting is an universal exponent. The global persistence exponent is found to be equal or larger than . This contrasts with previously known cases where . It is shown that in the special case of directed-bond percolation, P _l (t) can be related to a certain return probability of a directed percolation process with an active source (wet wall). Received: 15 December 1997 / Revised: 6 April 1998 / Accepted: 29 May 1998     

Block persistence

We define a block persistence probability p _l (t) as the probability that the order parameter integrated on a block of linear size l has never changed sign since the initial time in a phase-ordering process at finite temperature T<T _c . We argue that in the scaling limit of large blocks, where z is the growth exponent (), is the global (magnetization) persistence exponent and f(x) decays with the local (single spin) exponent for large x. This scaling is demonstrated at zero temperature for the diffusion equation and the large-n model, and generically it can be used to determine easily from simulations of coarsening models. We also argue that and the scaling function do not depend on temperature, leading to a definition of at finite temperature, whereas the local persistence probability decays exponentially due to thermal fluctuations. These ideas are applied to the study of persistence for conserved models. We illustrate our discussions by extensive numerical results. We also comment on the relation between this method and an alternative definition of at finite temperature recently introduced by Derrida [Phys. Rev. E 55, 3705 (1997)]. Received: 25 February 1998 / Revised: 24 July 1998 / Accepted: 27 July 1998     

Island distance in one-dimensional epitaxial growth

The typical island distance in submonolayer epitaxial growth depends on the growth conditions via an exponent . This exponent is known to depend on the substrate dimensionality, the dimension of the islands, and the size i^* of the critical nucleus for island formation. In this paper we study the dependence of on i^* in one-dimensional epitaxial growth. We derive that for and confirm this result by computer simulations. Received: 26 May 1998 / Accepted: 23 June 1998     

Log-periodic oscillations for biased diffusion of a polymer chain in a porous medium

A coarse-grained off-lattice bead-spring model is used to reveal the complex dynamics of a polymer chain in a quenched porous medium in the presence of an external field B. The behavior of the mean square displacement (MSD) of the center chain bead and that of the center of mass of the chain as a function of time is studied at different values of the barrier concentration C, the field strength B and the chain length N. In a field, important information on the way in which chains move between obstacles and overcome them is gained from the MSD vs. time analysis in the directions parallel and perpendicular to the flow. Instead of a steady approach to uniform drift-like motion at low C, for sufficiently strong field B we observe logarithmic oscillations in the effective exponents describing the time dependence of the MSD along and perpendicular to field. A common nature of this phenomenon with oscillatory behavior, observed earlier for biased diffusion of tracers on random lattices, is suggested. Received 7 August 1998     

Circular-like maps: sensitivity to the initial conditions, multifractality and nonextensivity

Dissipative one-dimensional maps may exhibit special points (e.g., chaos threshold) at which the Lyapunov exponent vanishes. Consistently, the sensitivity to the initial conditions has a power-law time dependence, instead of the usual exponential one. The associated exponent can be identified with 1/(1-q), where q characterizes the nonextensivity of a generalized entropic form currently used to extend standard, Boltzmann-Gibbs statistical mechanics in order to cover a variety of anomalous situations. It has been recently proposed (Lyra and Tsallis, Phys. Rev. Lett. 80, 53 (1998)) for such maps the scaling law , where and are the extreme values appearing in the multifractal function. We generalize herein the usual circular map by considering inflexions of arbitrary power z, and verify that the scaling law holds for a large range of z. Since, for this family of maps, the Hausdorff dimension d_f equals unity for all z in contrast with q which does depend on z, it becomes clear that d_f plays no major role in the sensitivity to the initial conditions. Received 5 February 1999     

Fluctuations, response and aging dynamics in a simple glass-forming liquid out of equilibrium

By means of molecular dynamics computer simulations we investigate the out of equilibrium relaxation dynamics of a simple glass former, a binary Lennard-Jones system, after a quench to low temperatures. We find that one-time quantities, such as the energy or the structure factor, show only a weak time dependence. By comparing the out of equilibrium structure factor with equilibrium data we find evidence that during the aging process the system remains in that part of phase space that mode-coupling theory classifies as liquid like. Two-times correlation functions show a strong time and waiting time dependence. For large and times corresponding to the early -relaxation regime the correlators approach the Edwards-Anderson value by means of a power-law in time. For large but fixed values of the relaxation dynamics in the -relaxation regime seems to be independent of the observable and temperature. The -relaxation shows a power-law dependence on time with an exponent which is independent of but depends on the observable. We find that at long times the correlation functions can be expressed as and compute the function h(t). This function is found to show a t-dependence which is a bit stronger than a logarithm and to depend on the observable considered. If the system is quenched to very low temperatures the relaxation dynamics at long times shows fast drops as a function of time. We relate these drops to relatively local rearrangements in which part of the sample relaxes its stress by a collective motion of 50-100 particles. Finally we discuss our measurements of the time dependent response function. We find that at long times the correlation functions and the response are not related by the usual fluctuation dissipation theorem but that this relation is similar to the one found for spin glasses with one step replica symmetry breaking. Received 17 May 1999     

Influence of hydrodynamics on the dynamics of a homopolymer

We present an analytical approach of the dynamics of a polymer when it is quenched from a solvent into a good or bad solvent. The dynamics is studied by means of a Langevin equation, first in the absence of hydrodynamic effect, then taking into account the hydrodynamic interactions with the solvent. The variation of the radius of gyration is studied as a function of time. In both cases, for the first stage of collapse or swelling, the evolution is described by a power law with a characteristic time proportional to N ^4/3 (N), where N is the number of monomers, without (with) hydrodynamic interactions. At larger times, scaling laws are derived for the diffusive relaxation time. Received: 10 March 1998 / Received in final form: 15 September 1998 / Accepted: 25 September 1998     

Growing mechanisms in the QKPZ equation aaand the DPD models

The roughening of interfaces moving in inhomogeneous media is investigated by numerical integration of the phenomenological stochastic differential equation proposed by Kardar, Parisi, and Zhang [Phys. Rev. Lett. 56, 889 (1986)] with quenched noise (QKPZ) [Phys. Rev. Lett. 74, 920 (1995)]. We express the evolution equations for the mean height and the roughness into two contributions: the local and the lateral one in order to compare them with the local and the lateral contributions obtained for the directed percolation depinning models (DPD) introduced independently by Tang and Leschhorn [Phys. Rev A 45, R8309 (1992)] and Buldyrev et al. [Phys. Rev A 45, R8313 (1992)]. These models are classified in the same universality class of the QKPZ although the mechanisms of growth are quite different. In the DPD models the lateral contribution is a coupled effect of the competition between the local growth and the lateral one. In these models the lateral contribution leads to an increasing of the roughness near the criticality while in the QKPZ equation this contribution always flattens the roughness. Received 7 April 2000 and Received in final form 7 March 2001     

Calculation of scattering from stretched copolymers using the tube model: a generalisation of the RPA

There are many experimental situations in which polymer chains are constrained or localised into a small region of space (e.g. melt chains confined to a “tube”, network chains pinned by crosslinks). We show that detailed consideration of the quenched variables is vital in these experiments. This paper provides a crucial link between microscopic models with localising constraints and scattering patterns by a generalisation of the Random Phase Approximation (RPA) which allows for quenched translational variables. A method is developed which deals with correlations between the quenched variables brought about by incompressiblity (for example, in a polymer melt there are correlations between tubes because of the interaction between chains). As an example, the generalised RPA is applied to models based on the Warner-Edwards picture of the tube. Theoretical results for a melt of H-shaped copolymers are compared with experimental scattering. Early results suggest that to fit the scattering we may be forced to relax one of the central assumptions of the tube model; that the tube deforms affinely, that all chains retract by the same amount or that the tube diameter does not couple to the strain. Received 26 October 1998 and Received in final form 19 March 1999     

Roughness and dynamics of a contact line of a viscous fluid on a disordered substrate

We have studied the roughness and the dynamics of the contact line of a viscous liquid on a disordered substrate. We have used photolithographic techniques to obtain a controlled disorder with a correlation length ξ = 10μm. Liquids with different viscosity were used: water and aqueous glycerol solution. We have found that the roughness W of the contact line depends neither on the viscosity nor on the velocity v of the contact line for v in the range 0.2-20μm/s. W is found to scale with the length L of the line as L ^ζ with a roughness exponent ζ = 0.51±0.03. This value is similar to the one obtained with superfluid helium. In the present experiment, we have checked that the motion of the contact line is actually overdamped, so that the phenomenological equation first proposed by Ertas and Kardar should be relevant. However, our measurement of ζ is in disagreement with the predicted value ζ = 0.39. We have also analyzed the avalanche-like motion of the contact line. We find that the size distribution does not follow a power law dependence. Received 18 April 2002     

Influence of Boron doping on microcrystalline silicon growth

Microcrystalline silicon (μc-Si:H) thin films with and without boron doping are deposited using the radio-frequency plasma-enhanced chemical vapour deposition method. The surface roughness evolutions of the silicon thin films are investigated using ex situ spectroscopic ellipsometry and an atomic force microscope. It is shown that the growth exponent β and the roughness exponent α are about 0.369 and 0.95 for the undoped thin film, respectively. Whereas, for the boron-doped μc-Si:H thin film, β increases to 0.534 and α decreases to 0.46 due to the shadowing effect.     

Short-time behavior of the kinetic spherical model with long-ranged interactions

The kinetic spherical model with long-ranged interactions and an arbitrary initial order m₀ quenched from a very high temperature to T is solved. In the short-time regime, the bulk order increases with a power law in both the critical and phase-ordering dynamics. To the latter dynamics, a power law for the relative order is found in the intermediate time-regime. The short-time scaling relations of small m₀ are generalized to an arbitrary m₀ and all the time larger than . The characteristic functions for the scaling of m₀ and for are obtained. The crossover between scaling regimes is discussed in detail. Received 17 September 1999