Dynamical Freezing of Relaxation to Equilibrium.

We provide evidence of an extremely slow thermalization occurring in the discrete nonlinear Schrödinger model. At variance with many similar processes encountered in statistical mechanics-typically ascribed to the presence of (free) energy barriers-here the slowness has a purely dynamical origin: it is due to the presence of an adiabatic invariant, which freezes the dynamics of a tall breather. Consequently, relaxation proceeds via rare events, where energy is suddenly released towards the background. We conjecture that this exponentially slow relaxation is a key ingredient contributing to the nonergodic behavior recently observed in the negative-temperature region of the discrete nonlinear Schrödinger equation.

Lattice models of random advection and diffusion and their statistics.

We study in detail a one-dimensional lattice model of a continuum, conserved field (mass) that is transferred deterministically between neighboring random sites. The model belongs to a wider class of lattice models capturing the joint effect of random advection and diffusion and encompassing as specific cases some models studied in the literature, such as those of Kang-Redner, Kipnis-Marchioro-Presutti, Takayasu-Taguchi, etc. The motivation for our setup comes from a straightforward interpretation of the advection of particles in one-dimensional turbulence, but it is also related to a problem of synchronization of dynamical systems driven by common noise. For finite lattices, we study both the coalescence of an initially spread field (interpreted as roughening), and the statistical steady-state properties. We distinguish two main size-dependent regimes, depending on the strength of the diffusion term and on the lattice size. Using numerical simulations and a mean-field approach, we study the statistics of the field. For weak diffusion, we unveil a characteristic hierarchical structure of the field. We also connect the model and the iterated function systems concept.

Nabiximols oromucosal spray in patients with multiple sclerosis-related bladder dysfunction: A prospective study.

BACKGROUND: Spasticity and urinary disturbances can profoundly impact the daily lives of persons with multiple sclerosis (pwMS). Cannabis has been associated with improvement in sphincteric disturbances. To our knowledge, few studies have evaluated the effect of nabiximols oromucosal spray (Sativex®) on urinary disturbances by instrumental methods. OBJECTIVES: This longitudinal study was conducted to assess the effect of nabiximols oromucosal spray on urinary disturbances by clinical and urodynamic evaluation in pwMS. MATERIALS AND METHODS: Neurological, spasticity, and quality of life (QoL) assessments were performed before (T0), and at one (T1) and six (T6) months after the start of nabiximols treatment. At these same time points, patients were assessed for urinary disturbances by the International Prostatic Symptoms Score (IPSS) and a urodynamic test evaluating maximum detrusor pressure (Pdet), bladder filling capacity (CCmax), uninhibited detrusor contractions (UDC), bladder volume at first desire (BVFD), post-void residual volume (PVR) and voluntary abdominal pressure (PA). RESULTS: Of 31 pwMS enrolled in the study, 25 reached T1 and 18 reached T6. Mean IPSS total score, its subscores, and IPSS QoL decreased significantly from T0 to T6 (p = 0.000), with no differences according to sex, age, MS type, disease duration and disability at baseline. Pdet improved significantly from T0 to T6 (p = 0.0171), and CCmax changed only marginally (p = 0.0494); results were similar in patient subgroups naïve to or previously exposed to urological treatment. All patients with overactive bladder showed improvement in their urodynamic assessment based on significant reduction of Pdet (p = 0.0138). In patients with mainly hypotonic bladder, mean Pdet decreased from T0 to T6 without reaching statistical significance; most urodynamic parameters showed a trend to improve. Mean numerical scale scores for MS spasticity, and for spasms, pain and tremors, decreased significantly from T0 to T6. The mean 'physical health composite' score of the MS Quality of Life-54 questionnaire increased significantly from T0 to T6 (p = 0.0126). DISCUSSION AND CONCLUSION: Our data suggest that nabiximols has an appreciable effect on ameliorating subjective perception of urinary disturbances and appears to have a positive effect on objective urodynamic parameters, particularly in patients with hyperactive bladder.

Coarsening scenarios in unstable crystal growth.

Crystal surfaces may undergo thermodynamical as well as kinetic, out-of-equilibrium instabilities. We consider the case of mound and pyramid formation, a common phenomenon in crystal growth and a long-standing problem in the field of pattern formation and coarsening dynamics. We are finally able to attack the problem analytically and get rigorous results. Three dynamical scenarios are possible: perpetual coarsening, interrupted coarsening, and no coarsening. In the perpetual coarsening scenario, mound size increases in time as L~t(n), where the coarsening exponent is n=1/3 when faceting occurs, otherwise n=1/4.

Disorder strength and field-driven ground state domain formation in artificial spin ice: experiment, simulation, and theory.

Quenched disorder affects how nonequilibrium systems respond to driving. In the context of artificial spin ice, an athermal system comprised of geometrically frustrated classical Ising spins with a twofold degenerate ground state, we give experimental and numerical evidence of how such disorder washes out edge effects and provide an estimate of disorder strength in the experimental system. We prove analytically that a sequence of applied fields with fixed amplitude is unable to drive the system to its ground state from a saturated state. These results should be relevant for other systems where disorder does not change the nature of the ground state.

Condensation induced by coupled transport processes.

Several lattice models display a condensation transition in real space when the density of a suitable order parameter exceeds a critical value. We consider one of such models with two conservation laws, in a onedimensional open setup where the system is attached to two external reservoirs. Both reservoirs impose subcritical boundary conditions at the chain ends. When such boundary conditions are equal, the system is in equilibrium below the condensation threshold and no condensate can appear. Instead, when the system is kept out of equilibrium, localization may arise in an internal portion of the lattice. We discuss the origin of this phenomenon, the relevance of the number of conservation laws, and the effect of the pinning of the condensate on the dynamics of the out-of-equilibrium state.

Diversity enabling equilibration: disorder and the ground state in artificial spin ice.

We report a novel approach to the question of whether and how the ground state can be achieved in square artificial spin ices where frustration is incomplete. We identify two sources of randomness that affect the approach to ground state: quenched disorder in the island response to fields and randomness in the sequence of driving fields. Numerical simulations show that quenched disorder can lead to final states with lower energy, and randomness in the sequence of driving fields always lowers the final energy attained by the system. We use a network picture to understand these two effects: disorder in island responses creates new dynamical pathways, and a random sequence of driving fields allows more pathways to be followed.

Finite-size localization scenarios in condensation transitions.

We consider the phenomenon of condensation of a globally conserved quantity H=∑_{i=1}^{N}ε_{i} distributed on N sites, occurring when the density h=H/N exceeds a critical density h_{c}. We numerically study the dependence of the participation ratio Y_{2}=〈ε_{i}^{2}〉/(Nh^{2}) on the size N of the system and on the control parameter δ=(h-h_{c}), for various models: (i) a model with two conservation laws, derived from the discrete nonlinear Schrödinger equation; (ii) the continuous version of the zero-range process class, for different forms of the function f(ε) defining the factorized steady state. Our results show that various localization scenarios may appear for finite N and close to the transition point. These scenarios are characterized by the presence or the absence of a minimum of Y_{2} when plotted against N and by an exponent γ≥2 defined through the relation N^{*}≃δ^{-γ}, where N^{*} separates the delocalized region (N≪N^{*}, Y_{2} vanishes with increasing N) from the localized region (N≫N^{*}, Y_{2} is approximately constant). We finally compare our results with the structure of the condensate obtained through the single-site marginal distribution.

Kinetics of the two-dimensional long-range Ising model at low temperatures.

We study the low-temperature domain growth kinetics of the two-dimensional Ising model with long-range coupling J(r)∼r^{-(d+σ)}, where d=2 is the dimensionality. According to the Bray-Rutenberg predictions, the exponent σ controls the algebraic growth in time of the characteristic domain size L(t), L(t)∼t^{1/z}, with growth exponent z=1+σ for σ<1 and z=2 for σ>1. These results hold for quenches to a nonzero temperature T>0 below the critical temperature T_{c}. We show that, in the case of quenches to T=0, due to the long-range interactions, the interfaces experience a drift which makes the dynamics of the system peculiar. More precisely, we find that in this case the growth exponent takes the value z=4/3, independently of σ, showing that it is a universal quantity. We support our claim by means of extended Monte Carlo simulations and analytical arguments for single domains.

Vertex dynamics in finite two-dimensional square spin ices.

Local magnetic ordering in artificial spin ices is discussed from the point of view of how geometrical frustration controls dynamics and the approach to steady state. We discuss the possibility of using a particle picture based on vertex configurations to interpret the time evolution of magnetic configurations. Analysis of possible vertex processes allows us to anticipate different behaviors for open and closed edges and the existence of different field regimes. Numerical simulations confirm these results and also demonstrate the importance of correlations and long-range interactions in understanding particle population evolution. We also show that a mean-field model of vertex dynamics gives important insights into finite size effects.

Quasideterministic dynamics, memory effects, and lack of self-averaging in the relaxation of quenched ferromagnets.

We discuss the interplay between the degree of dynamical stochasticity, memory persistence, and violation of the self-averaging property in the aging kinetics of quenched ferromagnets. We show that, in general, the longest possible memory effects, which correspond to the slowest possible temporal decay of the correlation function, are accompanied by the largest possible violation of self-averaging and a quasideterministic descent into the ergodic components. This phenomenon is observed in different systems, such as the Ising model with long-range interactions, including the mean-field, and the short-range random-field Ising model.

Modified Kuramoto-Sivashinsky equation: stability of stationary solutions and the consequent dynamics.

We study the effect of a higher-order nonlinearity in the standard Kuramoto-Sivashinsky equation: partial differentialxG(Hx). We find that the stability of steady states depends on dv/dq , the derivative of the interface velocity on the wave vector q of the steady state. If the standard nonlinearity vanishes, coarsening is possible, in principle, only if G is an odd function of Hx. In this case, the equation falls in the category of the generalized Cahn-Hilliard equation, whose dynamical behavior was recently studied by the same authors. Alternatively, if G is an even function of Hx, we show that steady-state solutions are not permissible.

Nonlinear dynamics in one dimension: a criterion for coarsening and its temporal law.

We develop a general criterion about coarsening for some classes of nonlinear evolution equations describing one-dimensional pattern-forming systems. This criterion allows one to discriminate between the situation where a coarsening process takes place and the one where the wavelength is fixed in the course of time. An intermediate scenario may occur, namely "interrupted coarsening." The power of the criterion on which a brief account has been given [Politi and Misbah, Phys. Rev. Lett. 92, 090601 (2004)], and which we extend here to more general equations, lies in the fact that the statement about the occurrence of coarsening, or selection of a length scale, can be made by only inspecting the behavior of the branch of steady state periodic solutions. The criterion states that coarsening occurs if lambda'(A)>0 while a length scale selection prevails if lambda'(A)<0, where lambda is the wavelength of the pattern and A is the amplitude of the profile (prime refers to differentiation). This criterion is established thanks to the analysis of the phase diffusion equation of the pattern. We connect the phase diffusion coefficient D(lambda) (which carries a kinetic information) to lambda'(A), which refers to a pure steady state property. The relationship between kinetics and the behavior of the branch of steady state solutions is established fully analytically. Another important and new result which emerges here is that the exploitation of the phase diffusion coefficient enables us to determine in a rather straightforward manner the dynamical coarsening exponent. Our calculation, based on the idea that |D(lambda)| approximately lambda2/t, is exemplified on several nonlinear equations, showing that the exact exponent is captured. We are not aware of another method that so systematically provides the coarsening exponent. Contrary to many situations where the one-dimensional character has proven essential for the derivation of the coarsening exponent, this idea can be used, in principle, at any dimension. Some speculations about the extension of the present results are outlined.

Transition to coarsening for confined one-dimensional interfaces with bending rigidity.

We discuss the nonlinear dynamics and fluctuations of interfaces with bending rigidity under the competing attractions of two walls with arbitrary permeabilities. This system mimics the dynamics of confined membranes. We use a two-dimensional hydrodynamic model, where membranes are effectively one-dimensional objects. In a previous work [T. Le Goff et al., Phys. Rev. E 90, 032114 (2014)], we have shown that this model predicts frozen states caused by bending rigidity-induced oscillatory interactions between kinks (or domain walls). We here demonstrate that in the presence of tension, potential asymmetry, or thermal noise, there is a finite threshold above which frozen states disappear, and perpetual coarsening is restored. Depending on the driving force, the transition to coarsening exhibits different scenarios. First, for membranes under tension, small tensions can only lead to transient coarsening or partial disordering, while above a finite threshold, membrane oscillations disappear and perpetual coarsening is found. Second, potential asymmetry is relevant in the nonconserved case only, i.e., for permeable walls, where it induces a drift force on the kinks, leading to a fast coarsening process via kink-antikink annihilation. However, below some threshold, the drift force can be balanced by the oscillatory interactions between kinks, and frozen adhesion patches can still be observed. Finally, at long times, noise restores coarsening with standard exponents depending on the permeability of the walls. However, the typical time for the appearance of coarsening exhibits an Arrhenius form. As a consequence, a finite noise amplitude is needed in order to observe coarsening in observable time.

Process of irreversible nucleation in multilayer growth. I. Failure of the mean-field approach.

The formation of stable dimers on top of terraces during epitaxial growth is investigated in detail. In this paper we focus on mean-field theory, the standard approach to study nucleation. Such theory is shown to be unsuitable for the present problem, because it is equivalent to considering adatoms as independent diffusing particles. This leads to an overestimate of the correct nucleation rate by a factor N, which has a direct physical meaning: on average, a visited lattice site is visited N times by a diffusing adatom. The dependence of N on the size of the terrace and on the strength of step-edge barriers is derived from well-known results for random walks. The spatial distribution of nucleation events is shown to be different from the mean-field prediction, for the same physical reason. In the following paper we develop an exact treatment of the problem.

Process of irreversible nucleation in multilayer growth. II. Exact results in one and two dimensions.

We study irreversible dimer nucleation on top of terraces during epitaxial growth in one and two dimensions, for all values of the step-edge barrier. The problem is solved exactly by transforming it into a first passage problem for a random walker in a higher-dimensional space. The spatial distribution of nucleation events is shown to differ markedly from the mean-field estimate except in the limit of very weak step-edge barriers. The nucleation rate is computed exactly, including numerical prefactors.

Universality classes for unstable crystal growth.

Universality has been a key concept for the classification of equilibrium critical phenomena, allowing associations among different physical processes and models. When dealing with nonequilibrium problems, however, the distinction in universality classes is not as clear and few are the examples, such as phase separation and kinetic roughening, for which universality has allowed to classify results in a general spirit. Here we focus on an out-of-equilibrium case, unstable crystal growth, lying in between phase ordering and pattern formation. We consider a well-established 2+1-dimensional family of continuum nonlinear equations for the local height h(x,t) of a crystal surface having the general form ∂_{t}h(x,t)=-∇·[j(∇h)+∇(∇^{2}h)]: j(∇h) is an arbitrary function, which is linear for small ∇h, and whose structure expresses instabilities which lead to the formation of pyramidlike structures of planar size L and height H. Our task is the choice and calculation of the quantities that can operate as critical exponents, together with the discussion of what is relevant or not to the definition of our universality class. These aims are achieved by means of a perturbative, multiscale analysis of our model, leading to phase diffusion equations whose diffusion coefficients encapsulate all relevant information on dynamics. We identify two critical exponents: (i) the coarsening exponent, n, controlling the increase in time of the typical size of the pattern, L∼t^{n}; (ii) the exponent ß, controlling the increase in time of the typical slope of the pattern, M∼t^{ß}, where M≈H/L. Our study reveals that there are only two different universality classes, according to the presence (n=1/3, ß=0) or the absence (n=1/4, ß>0) of faceting. The symmetry of the pattern, as well as the symmetry of the surface mass current j(∇h) and its precise functional form, is irrelevant. Our analysis seems to support the idea that also space dimensionality is irrelevant.

Frozen states and order-disorder transition in the dynamics of confined membranes.

The adhesion dynamics of a membrane confined between two permeable walls is studied using a two-dimensional hydrodynamic model. The membrane morphology decomposes into adhesion patches on the upper and the lower walls and obeys a nonlinear evolution equation that resembles that of phase-separation dynamics, which is known to lead to coarsening, i.e., to the endless growth of the adhesion patches. However, due to the membrane bending rigidity, the system evolves toward a frozen state without coarsening. This frozen state exhibits an order-disorder transition when increasing the permeability of the walls.

Large-scale patterns of genetic variation in a female-biased dispersing passerine: the importance of sex-based analyses.

Dispersal affects the distribution, dynamics and genetic structure of natural populations, and can be significantly different between sexes. However, literature records dealing with the dispersal of migratory birds are scarce, as migratory behaviour can notably complicate the study of dispersal. We used the barn swallow Hirundo rustica as model taxon to investigate patterns of genetic variability in males and in females of a migratory species showing sex-biased dispersal. We collected blood samples (n = 186) over the period 2006 to 2011 from adults (H. r. rustica subspecies) nesting in the same breeding site at either high (Ireland, Germany and Russia) or low (Spain, Italy and Cyprus) latitude across Europe. We amplified the Chromo Helicase DNA gene in all birds in order to warrant a sex-balanced sample size (92 males, 94 females). We investigated both uniparental (mitochondrial ND2 gene) and biparental (microsatellite DNA: 10 loci) genetic systems. The mtDNA provided evidence for demographic expansion yet no significant partition of the genetic variability was disclosed. Nevertheless, a comparatively distant Russian population investigated in another study, whose sequences were included in the present dataset, significantly diverged from all other ones. Different to previous studies, microsatellites highlighted remarkable genetic structure among the studied populations, and pointed to the occurrence of differences between male and female barn swallows. We produced evidence for non-random patterns of gene flow among barn swallow populations probably mediated by female natal dispersal, and we found significant variability in the philopatry of males of different populations. Our data emphasize the importance of taking into account the sex of sampled individuals in order to obtain reliable inferences on species characterized by different patterns of dispersal between males and females.

Coarsening dynamics in one dimension: the phase diffusion equation and its numerical implementation.

Many nonlinear partial differential equations (PDEs) display a coarsening dynamics, i.e., an emerging pattern whose typical length scale L increases with time. The so-called coarsening exponent n characterizes the time dependence of the scale of the pattern, L(t)≈t(n), and coarsening dynamics can be described by a diffusion equation for the phase of the pattern. By means of a multiscale analysis we are able to find the analytical expression of such diffusion equations. Here, we propose a recipe to implement numerically the determination of D(λ), the phase diffusion coefficient, as a function of the wavelength λ of the base steady state u(0)(x). D carries all information about coarsening dynamics and, through the relation |D(L)|=/~L(2)/t, it allows us to determine the coarsening exponent. The main conceptual message is that the coarsening exponent is determined without solving a time-dependent equation, but only by inspecting the periodic steady-state solutions. This provides a much faster strategy than a orward time-dependent calculation. We discuss our method for several different PDEs, both conserved and not conserved.