Dynamic exponent in extremal models of pinning

The depinning transition of a front moving in a time-independent random potential is studied. The temporal development of the overall roughness w(L,t) of an initially flat front, , is the classical means to have access to the dynamic exponent. However, in the case of front propagation in quenched disorder via extremal dynamics, we show that the initial increase in front roughness implies an extra dependence over the system size which comes from the fact that the activity is essentially localized in a narrow region of space. We propose an analytic expression for the exponent and confirm this for different models (crack front propagation, Edwards-Wilkinson model in a quenched noise etc.). Received 27 August 1999     

Spinodal decomposition of an ABv model alloy: Patterns at unstable surfaces

We develop mean-field kinetic equations for a lattice gas model of a binary alloy with vacancies (ABv model) in which diffusion takes place by a vacancy mechanism. These equations are applied to the study of phase separation of finite portions of an unstable mixture immersed in a stable vapor. Due to a larger mobility of surface atoms, the most unstable modes of spinodal decomposition are localized at the vapor-mixture interface. Simulations show checkerboard-like structures at the surface or surface-directed spinodal waves. We determine the growth rates of bulk and surface modes by a linear stability analysis and deduce the relation between the parameters of the model and the structure and length scale of the surface patterns. The thickness of the surface patterns is related to the concentration fluctuations in the initial state. Received 28 October 1998     

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     

DMRG studies on interchain coupling model for quasi-one-dimensional organic magnet

We introduce a solid-on-solid growth process which evolves by random deposition of dimers, surface diffusion, and evaporation of monomers from the edges of plateaus. It is shown that the model exhibits a robust transition from a smooth to a rough phase. The roughening transition is driven by an absorbing phase transition at the bottom layer of the interface, which displays the same type of critical behavior as the pair contact process with diffusion 2A↦3A, 2A↦. Received 14 October 2002 Published online 14 February 2003 RID="a" ID="a"e-mail: Haye.Hinrichsen@physik.uni-wuppertal.de     

Disorder driven roughening transitions of elastic manifolds and periodic elastic media

The simultaneous effect of both disorder and crystal-lattice pinning on the equilibrium behavior of oriented elastic objects is studied using scaling arguments and a functional renormalization group technique. Our analysis applies to elastic manifolds, e.g., interfaces, as well as to periodic elastic media, e.g., charge-density waves or flux-line lattices. The competition between both pinning mechanisms leads to a continuous, disorder driven roughening transition between a flat state where the mean relative displacement saturates on large scales and a rough state with diverging relative displacement. The transition can be approached by changing the impurity concentration or, indirectly, by tuning the temperature since the pinning strengths of the random and crystal potential have in general a different temperature dependence. For D dimensional elastic manifolds interacting with either random-field or random-bond disorder a transition exists for 2<D<4, and the critical exponents are obtained to lowest order in . At the transition, the manifolds show a superuniversal logarithmic roughness. Dipolar interactions render lattice effects relevant also in the physical case of D=2. For periodic elastic media, a roughening transition exists only if the ratio p of the periodicities of the medium and the crystal lattice exceeds the critical value . For p<p _c the medium is always flat. Critical exponents are calculated in a double expansion in and and fulfill the scaling relations of random field models. Received 28 August 1998     

Dynamical linked cluster expansions with applications to disordered systems

Dynamical linked cluster expansions are linked cluster expansions with hopping parameter terms endowed with their own dynamics. We discuss physical applications to systems with annealed and quenched disorder. Examples are the bond-diluted Ising model and the Sherrington-Kirkpatrick spin glass. We derive the rules and identify the full set of graphs that contribute to the series in the quenched case. This way it becomes possible to avoid the vague extrapolation from positive integer n to n = 0, that usually goes along with an application of the replica trick. Received 13 December 2001 Published online 25 June 2002     

Adsorption of a Gaussian random copolymer chain at an interface

We consider the adsorption of an isolated, Gaussian, random, and quenched copolymer chain at an interface. We first propose a simple analytical method to obtain the adsorption/depletion transition, by averaging over the disorder the partition function instead of the free energy. The adsorption thresholds obtained by previous authors at a solid/liquid and at a liquid/liquid interface for multicopolymer chains can be rederived using this method. We also compare the adsorption thresholds obtained for bimodal and for Gaussian disorder; they only agree for small disorder. We focus on the specific case of an ideally flat asymmetric liquid/liquid interface, and consider the situation where the chain is composed of monomers of two different chemical species A and B. The replica method is developed for this case. We show that the Hartree approximation, coupled to a replica symmetry assumption, leads to the same adsorption thresholds as obtained from our general method. In order to describe the properties of the adsorbed (or depleted) chain, we develop a new approximation for long chains, within the framework of the replica theory. In most cases, the behavior of a random copolymer chain can be mapped onto that of a homopolymer chain at an asymmetric attractive interface. The values of the effective adsorption energy are different for a random and a periodic copolymer chain. Finally, we consider the case of uncorrelated annealed disorder. The behavior of an annealed chain can be mapped onto that of a homopolymer chain at an asymmetric non attractive interface; hence, an annealed chain cannot adsorb at an asymmetric interface. Received 21 January 1999     

Scaling,propagation, and kinetic roughening of flame fronts in random media

We introduce a model of two coupled reaction-diffusion equations to describe the dynamics and propagation of flame fronts in random media. The model incorporates heat diffusion, its dissipation, and its production through coupling to the background reactant density. We first show analytically and numerically that there is a finite critical value of the background density below which the front associated with the temperature field stops propagating. The critical exponents associated with this transition are shown to be consistent with meanfield theory of percolation. Second, we study the kinetic roughening associated with a moving planar flame front above the critical density. By numerically calculating the time-dependent width and equal-time height correlation function of the front, we demonstrate that the roughening process belongs to the universality class of the Kardar-Parisi-Zhang interface equation. Finally, we show how this interface equation can be analytically derived from our model in the limit of almost uniform background density.     

Adatom mobility for the solid-on-solid model

The relaxation of a crystal surface through surface diffusion is studied within the solid-on-solid model. Two types of (conserved) dynamics are considered. ForArrhenius dynamics we show that the relevant transport coefficient, the adatom mobility, has a simple analytic form: It is independent of orientation, and depends exponentially on the inverse temperature, for any surface dimensionalityd. Together with the expression for the orientation-dependent stiffness this completely determines the macroscopic evolution equation for the surface. The predictions of the macroscopic theory are checked against simulations of profile evolution and roughening ind=1. For one-dimensionalMetropolis dynamics we provide an upper bound on the adatom mobility and obtain numerical estimates of its actual value, which indicate a nontrivial orientation dependence in this case. An alternative derivation of the macroscopic dynamics directly from the master equation is presented and discussed in relation to previous approximate work.Dedicated to Professor H. Wagner on the occasion of his 60th birthday     

Molecular photodynamics involved in multi-photon excitation fluorescence microscopy

We present a simple model for calculating the fluorescence generated by the multi-photon excitation (MPE) of molecules in solution. The model takes into account internal molecular dynamics such as ground-state depletion due to inter-system crossing (ISC), as well as external molecular dynamics associated with diffusion into and out of an excitation volume confined in 3-dimensions. Internal and external molecular dynamics are combined by using a technique of linearization of a modified diffusion equation which takes into account the possibility of concentration depletion due to photobleaching. In addition, we discuss the phenomenon of pulse saturation which effectively limits the molecular excitation rate constant in the case of short pulsed excitation. Our results are specifically applied in the context of fluorescence autocorrelation functions and single-molecule detection. In the latter case, we discuss some consequences of high-order multi-photon photobleaching. Finally, we include three appendices to rigorously define the temporal and spatial profiles of an arbitrary excitation beam, and also to discuss some properties of an exact evaluation of concentration depletion due to photobleaching. Received: 9 March 1998 / Accepted: 20 April 1998     

Reptational dynamics of a polymer chain is stable against kinematic disorder

We study the diffusive motion of a (non-selfinteracting) chain through a quenched random environment, constructed such that it influences only the local dynamics but not the equilibrium configuration of the chain. Our Monte Carlo results show that this type of disorder, which we call kinematic, does not ruin reptation. This is in sharp contrast to disorder including also entropic traps and it supports the view that reptation prevails in melts, where in contrast to a gel entropic trapping is absent. Our data show the characteristic features of reptation, irrespective of the dilution or randomness of the kinematic obstacles. Our Monte Carlo results are in quantitative agreement with our recent detailed analytical evaluation of the reptation model (J. Stat. Phys. 90, 1325 (1998)). The analysis suggests that we effectively see reptation of a “blob”-chain, where the size of the blob rapidly increases with decreasing obstacle concentration. Received 30 November 1998 and Received in final form 18 December 1998     

Kinetic lattice models of disorder

We study a class of stochastic Ising (or interacting particle) systems that exhibit a spatial distribution of impurities that change with time. It may model, for instance, steady nonequilibrium conditions of the kind that may be induced by diffusion in some disordered materials. Different assumptions for the degree of coupling between the spin and the impurity configurations are considered. Two interesting well-defined limits for impurities that behave autonomously are (i) the standard (i.e., quenched) bond-diluted, random-field, random-exchange, and spin-glass Ising models, and (ii) kinetic variations of these standard cases in which conflicting kinetics simulate fast and random diffusion of impurities. A generalization of the Mattis model with disorder that describes a crossover from the equilibrium case (i) to the nonequilibrium case (ii) and the microscopic structure of a generalized heat bath are explicitly worked out as specific realizations of our class of models. We sketch a simple classification of transition rates for the time evolution of the spin configuration based on the critical behavior that is exhibited by the models in case (ii). The latter are shown to have an exact solution for any lattice dimension for some special choice of rates.     

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     

On dewetting dynamics of solid films of hydrogen isotopes and its influence on tritium ß spectroscopy

The dewetting dynamics of solid films of hydrogen isotopes, quench-condensed on a graphite substrate, was measured at various temperatures below desorption by observing the stray light from the film. A schematic model describing the dewetting process by surface diffusion is presented, which agrees qualitatively with our data. The activation energies of different hydrogen isotopes for surface diffusion were determined. The time constant for dewetting of a quench-condensed film at the working temperature of 1.86 K of the mainz neutrino mass experiment was extrapolated. Received 30 December 1999     

Minority opinion spreading in random geometry

The dynamics of spreading of the minority opinion in public debates (a reform proposal, a behavior change, a military retaliation) is studied using a diffusion reaction model. People move by discrete step on a landscape of random geometry shaped by social life (offices, houses, bars, and restaurants). A perfect world is considered with no advantage to the minority. A one person-one argument principle is applied to determine locally individual mind changes. In case of equality, a collective doubt is evoked which in turn favors the Status Quo. Starting from a large in favor of the proposal initial majority, repeated random size local discussions are found to drive the majority reversal along the minority hostile view. Total opinion refusal is completed within few days. Recent national collective issues are revisited. The model may apply to rumor and fear propagation. Received 23 January 2002     

Stochastic dynamics in a two-level model of disorder: comparison of mean-field and exact solutions

Stochastic dynamics in the presence of quenched disorder (e.g., diffusion in a random medium) is generally treated in a suitable mean-field or effective medium approximation. While numerical simulations may help determine the accuracy of such approximations in specific models, there are relatively few instances in which analytic solutions are possible, to enable a precise comparison to be made with the mean-field results. We consider in this paper a simple but general model of quenched disorder in which a system variablex jumps stochastically between two valuesx _a andx _b . However, in each level there occurs with a certain probability a branch (or internal) state into which the system may fall, and from which a jump to the other level is possible only after a return to the original (or ‘active’) state. Four different configurations of the states of the system are thus possible, and the transitions between the states are governed by Markovian transition probabilities. The moments ofx and its autocorrelation function are computed in each case, and then configuration-averaged over the four realizations. This represents the exact solution. Next, a mean-field theory of the dynamics is developed: this turns out to involve an effective waiting-time density at each of the two levels that is non-exponential in time, so that the mean-field dynamics is a non-Markovian alternating renewal process. The moments and autocorrelation ofx are again computed, and compared with the exact solutions. The extent of the differences at both short and long times is elucidated, and a numerical comparison is presented for the case of maximal disorder.     

Description of far-from-equilibrium processes by mean-field lattice gas models

Mean-field kinetic equations are a valuable tool to study the atomic dynamics and spin dynamics of simple lattice gas and Ising models. They can be derived from the microscopic master equation of the system and contain analytical expressions for kinetic coefficients and thermodynamic quantities which are usually introduced phenomenologically. We review several methods to obtain such equations, and discuss applications to the dynamics of order–disorder transitions, spinodal decomposition, and dendritic growth in the isothermal or chemical model. In the case of dendritic growth we show that the mean-field kinetic equations are equivalent to standard continuum equations for this problem and derive expressions for macroscopic quantities, e.g. the surface tension and kinetic coefficients, as functions of the microscopic order parameters. In spinodal decomposition, we focus our attention on the vacancy mechanism, which is a more faithful picture of diffusion in solids than the more widely examined exchange mechanism. We study the interfaces between an unstable mixture and a stable ‘vapour’ phase, and analyse surface modes that lead to specific surface patterns. For order–disorder transitions, studied in the framework of a repulsive two-sublattice model, we derive sets of coupled equations for the mean concentration (a conserved quantity) and for the occupational difference between the two sublattices emerging from the symmetry breaking due to ordering (non-conserved order parameter). These equations are applied to transport in the presence of ordered domains. Finally, we discuss the possibilities of improving the simple mean-field approximation by density functional theories and various forms of the dynamic pair approximation, including the path-probability method.     

Classical diffusion in strong random media

We study classical diffusion of particles in random media. Although many of our results are general, we focus on the case of an ion in a three-dimensional medium with random, quenched charge centers obeying bulk charge neutrality. Within a functional-integral framework, we calculate the effective diffusion coefficients by first-order and second-order self-consistent perturbation theory (with a Gaussian reference in both cases). We also carry out a one-loop order momentum space renormalization group calculation. The self-consistent methods are complicated numerically and fail beyond intermediate disorder strengths. In contrast, the renormalization group calculation gives an analytical result that appears valid even to high disorder strengths. The methodology, generally applicable to a quantitative calculation of effective diffusion coefficients in disordered media, resolves deficiencies in self-consistent perturbation theory approaches to this class of problems.     

Domain growth in the three-dimensional dilute Ising model

We investigate the kinetics of domain growth in the three-dimensional Ising model with quenched random site dilution, using Monte Carlo simulation technique. A crossover from the power law growth regime to a much slower growth observed in our simulation is interpreted through the roughening of the interfaces by the quenched impurities. The results are also compared with the corresponding results in two dimensions.     

Damage in fiber bundle models

We introduce a continuous damage fiber bundle model and compare its behavior with that of dry fiber bundles. Several interesting constitutive behaviors, such as plasticity, are found in this model depending on the value of the damage parameter and on the form of the disorder distribution. We compare the constitutive behavior of global load transfer models, obtained analytically, with local load transfer models numerical simulations. The evolution of the damage is studied analyzing the cluster statistics for dry and continous damage fiber bundles. Finally, it is shown that quenched random thresholds enhance damage localization. Received 27 March 2000