Critical behavior of charge-density waves below threshold: Numerical and scaling analysis

Shapiro steps and chaos in the Frenkel-Kontorova model with substrate lateral vibration

Numerical simulations are used to examine the dynamics of the dc-driven Frenkel-Kontorova model with an oscillation substrate subjected to lateral periodic excitations in overdamped and underdamped cases, respectively. The results reveal that the system exhibits frequency locking and chaotic behaviors due to the fact that the lateral vibration of the substrate potential introduces an additional frequency and degree of freedom into the system. In the overdamped case, we show that the appearance of subharmonic Shapiro steps can be attributed to the deformation of the substrate potential or inertia. The characteristics of the steps are significantly affected by the amplitude and frequency of the lateral vibration. When the vibration frequency is relatively high, the change of the width of the first harmonic Shapiro step with increasing amplitude exhibits a Bessel-type oscillation, but the oscillation deviates from the Bessel curve at lower frequencies. In the amplitude dependence, although the oscillatory behavior of the critical depinning force at the high frequency is anomalous, local maxima (minima) of the first step width correspond to local minima (maxima) of the critical depinning force, and the largest Lyapunov exponent obtained in the pinned state represents a mirror relationship of the critical depinning force. In contrast to the overdamped system, the underdamped one exhibits both subharmonic Shapiro steps and chaotic behaviors. We show the increased inertia of the latter system plays an important role in suppressing the emergence of subharmonic steps, which is opposite to the result of the former. When the dc force changes, chaos appears not only between adjacent subharmonic Shapiro steps but also on some specific steps where chaos should be avoided. The variation regularity of the first step width and the critical depinning force is thus annihilated in vibration amplitude and frequency dependence. However, superlubricity can be achieved by careful adjustment of vibration amplitude and frequency. The findings can serve as a theoretical guideline for technological applications such as device building and voltage standards.

Criticality in sheared, disordered solids. I. Rate effects in stress and diffusion

Rate effects in sheared disordered solids are studied using molecular dynamics simulations of binary Lennard-Jones glasses in two and three dimensions. In the quasistatic (QS) regime, systems exhibit critical behavior: the magnitudes of avalanches are power-law distributed with a maximum cutoff that diverges with increasing system size L. With increasing rate, systems move away from the critical yielding point and the average flow stress rises as a power of the strain rate with exponent 1/β, the Herschel-Bulkley exponent. Finite-size scaling collapses of the stress are used to measure β as well as the exponent ν which characterizes the divergence of the correlation length. The stress and kinetic energy per particle experience fluctuations with strain that scale as L^{-d/2}. As the largest avalanche in a system scales as L^{α}, this implies α<d/2. The diffusion rate of particles diverges as a power of decreasing rate before saturating in the QS regime. A scaling theory for the diffusion is derived using the QS avalanche rate distribution and generalized to the finite strain rate regime. This theory is used to collapse curves for different system sizes and confirm β/ν.

Anisotropic avalanches and critical depinning of three-dimensional magnetic domain walls

Simulations with more than 10^{12} spins are used to study the motion of a domain wall driven through a three-dimensional random-field Ising magnet (RFIM) by an external field H. The interface advances in a series of avalanches whose size diverges at a critical external field H_{c}. Finite-size scaling is applied to determine critical exponents and test scaling relations. Growth is intrinsically anisotropic with the height of an avalanche normal to the interface ℓ_{⊥} scaling as the width along the interface ℓ_{∥} to a power χ=0.85±0.01. The total interface roughness is consistent with self-affine scaling with a roughness exponent ζ≈χ that is much larger than values found previously for the RFIM and related models that explicitly break orientational symmetry by requiring the interface to be single-valued. Because the RFIM maintains orientational symmetry, the interface develops overhangs that may surround unfavorable regions to create uninvaded bubbles. Overhangs complicate measures of the roughness exponent but decrease in importance with increasing system size.

Depinning Transition of Charge-Density Waves: Mapping onto O(n) Symmetric ϕ^{4} Theory with n→-2 and Loop-Erased Random Walks

Driven periodic elastic systems such as charge-density waves (CDWs) pinned by impurities show a nontrivial, glassy dynamical critical behavior. Their proper theoretical description requires the functional renormalization group. We show that their critical behavior close to the depinning transition is related to a much simpler model, O(n) symmetric ϕ^{4} theory in the unusual limit of n→-2. We demonstrate that both theories yield identical results to four-loop order and give both a perturbative and a nonperturbative proof of their equivalence. As we show, both theories can be used to describe loop-erased random walks (LERWs), the trace of a random walk where loops are erased as soon as they are formed. Remarkably, two famous models of non-self-intersecting random walks, self-avoiding walks and LERWs, can both be mapped onto ϕ^{4} theory, taken with formally n=0 and n→-2 components. This mapping allows us to compute the dynamic critical exponent of CDWs at the depinning transition and the fractal dimension of LERWs in d=3 with unprecedented accuracy, z(d=3)=1.6243±0.001, in excellent agreement with the estimate z=1.62400±0.00005 of numerical simulations.

Inertial effects in the dc+ac driven underdamped Frenkel-Kontorova model: Subharmonic steps, chaos, and hysteresis

The effects of inertial terms on the dynamics of the dc+ac driven Frenkel-Kontorova model were examined. As the mass of particles was varied, the response of the system to the driving forces and appearance of the Shapiro steps were analyzed in detail. Unlike in the overdamped case, the increase of mass led to the appearance of the whole series of subharmonic steps in the staircase of the average velocity as a function of average driving force in any commensurate structure. At certain values of parameters, the subharmonic steps became separated by chaotic windows while the whole structure retained scaling similar to the original staircase. The mass of the particles also determined their sensitivity to the forces governing their dynamics. Depending on their mass, they were found to exhibit three types of dynamics, from dynamical mode-locking with chaotic windows, through to a typical dc response, to essentially a free-particle response. Examination of this dynamics in both the upforce and downforce directions showed that the system may not only exhibit hysteresis, but also that large Shapiro steps may appear in the downforce direction, even in cases for which no dynamical mode-locking occurred in the upforce direction.

Localization of soft modes at the depinning transition

We characterize the soft modes of the dynamical matrix at the depinning transition, and compare the matrix with the properties of the Anderson model (and long-range generalizations). The density of states at the edge of the spectrum displays a universal linear tail, different from the Lifshitz tails. The eigenvectors are instead very similar in the two matrix ensembles. We focus on the ground state (soft mode), which represents the epicenter of avalanche instabilities. We expect it to be localized in all finite dimensions, and make a clear connection between its localization length and the Larkin length of the depinning model. In the fully connected model, we show that the weak-strong pinning transition coincides with a peculiar localization transition of the ground state.

Devil's staircase and the absence of chaos in the dc- and ac-driven overdamped Frenkel-Kontorova model

The devil's staircase structure arising from the complete mode locking of an entirely nonchaotic system, the overdamped dc+ac driven Frenkel-Kontorova model with deformable substrate potential, was observed. Even though no chaos was found, a hierarchical ordering of the Shapiro steps was made possible through the use of a previously introduced continued fraction formula. The absence of chaos, deduced here from Lyapunov exponent analyses, can be attributed to the overdamped character and the Middleton no-passing rule. A comparative analysis of a one-dimensional stack of Josephson junctions confirmed the disappearance of chaos with increasing dissipation. Other common dynamic features were also identified through this comparison. A detailed analysis of the amplitude dependence of the Shapiro steps revealed that only for the case of a purely sinusoidal substrate potential did the relative sizes of the steps follow a Farey sequence. For nonsinusoidal (deformed) potentials, the symmetry of the Stern-Brocot tree, depicting all members of particular Farey sequence, was seen to be increasingly broken, with certain steps being more prominent and their relative sizes not following the Farey rule.

Scaling ansatz for the jamming transition

We propose a Widom-like scaling ansatz for the critical jamming transition. Our ansatz for the elastic energy shows that the scaling of the energy, compressive strain, shear strain, system size, pressure, shear stress, bulk modulus, and shear modulus are all related to each other via scaling relations, with only three independent scaling exponents. We extract the values of these exponents from already known numerical or theoretical results, and we numerically verify the resulting predictions of the scaling theory for the energy and residual shear stress. We also derive a scaling relation between pressure and residual shear stress that yields insight into why the shear and bulk moduli scale differently. Our theory shows that the jamming transition exhibits an emergent scale invariance, setting the stage for the potential development of a renormalization group theory for jamming.

Pinning Susceptibility: The Effect of Dilute, Quenched Disorder on Jamming

We study the effect of dilute pinning on the jamming transition. Pinning reduces the average contact number needed to jam unpinned particles and shifts the jamming threshold to lower densities, leading to a pinning susceptibility, χ_{p}. Our main results are that this susceptibility obeys scaling form and diverges in the thermodynamic limit as χ_{p}∝|ϕ-ϕ_{c}^{∞}|^{-γ_{p}} where ϕ_{c}^{∞} is the jamming threshold in the absence of pins. Finite-size scaling arguments yield these values with associated statistical (systematic) errors γ_{p}=1.018±0.026(0.291) in d=2 and γ_{p}=1.534±0.120(0.822) in d=3. Logarithmic corrections raise the exponent in d=2 to close to the d=3 value, although the systematic errors are very large.

Reversibility and criticality in amorphous solids

The physical processes governing the onset of yield, where a material changes its shape permanently under external deformation, are not yet understood for amorphous solids that are intrinsically disordered. Here, using molecular dynamics simulations and mean-field theory, we show that at a critical strain amplitude the sizes of clusters of atoms undergoing cooperative rearrangements of displacements (avalanches) diverges. We compare this non-equilibrium critical behaviour to the prevailing concept of a 'front depinning' transition that has been used to describe steady-state avalanche behaviour in different materials. We explain why a depinning-like process can result in a transition from periodic to chaotic behaviour and why chaotic motion is not possible in pinned systems. These findings suggest that, at least for highly jammed amorphous systems, the irreversibility transition may be a side effect of depinning that occurs in systems where the disorder is not quenched.

Avalanche dynamics of elastic interfaces

Slowly driven elastic interfaces, such as domain walls in dirty magnets, contact lines wetting a nonhomogeneous substrate, or cracks in brittle disordered material proceed via intermittent motion, called avalanches. Here we develop a field-theoretic treatment to calculate, from first principles, the space-time statistics of instantaneous velocities within an avalanche. For elastic interfaces at (or above) their (internal) upper critical dimension d≥d(uc) (d(uc)=2,4 respectively for long-ranged and short-ranged elasticity) we show that the field theory for the center of mass reduces to the motion of a point particle in a random-force landscape, which is itself a random walk [Alessandro, Beatrice, Bertotti, and Montorsi (ABBM) model]. Furthermore, the full spatial dependence of the velocity correlations is described by the Brownian-force model (BFM) where each point of the interface sees an independent Brownian-force landscape. Both ABBM and BFM can be solved exactly in any dimension d (for monotonous driving) by summing tree graphs, equivalent to solving a (nonlinear) instanton equation. We focus on the limit of slow uniform driving. This tree approximation is the mean-field theory (MFT) for realistic interfaces in short-ranged disorder, up to the renormalization of two parameters at d=d(uc). We calculate a number of observables of direct experimental interest: Both for the center of mass, and for a given Fourier mode q, we obtain various correlations and probability distribution functions (PDF's) of the velocity inside an avalanche, as well as the avalanche shape and its fluctuations (second shape). Within MFT we find that velocity correlations at nonzero q are asymmetric under time reversal. Next we calculate, beyond MFT, i.e., including loop corrections, the one-time PDF of the center-of-mass velocity u[over ·] for dimension d<d(uc). The singularity at small velocity P(u[over ·])~1/u[over ·](a) is substantially reduced from a=1 (MFT) to a=1-2/9(4-d)+... (short-ranged elasticity) and a=1-4/9(2-d)+... (long-ranged elasticity). We show how the dynamical theory recovers the avalanche-size distribution, and how the instanton relates to the response to an infinitesimal step in the force.

Dimer motion on a periodic substrate: spontaneous symmetry breaking and absolute negative mobility

We consider two coupled particles moving along a periodic substrate potential with negligible inertia effects (overdamped limit). Even when the particles are identical and the substrate spatially symmetric, a sinusoidal external driving of appropriate amplitude and frequency may lead to spontaneous symmetry breaking in the form of a permanent directed motion of the dimer. Thermal noise restores ergodicity and thus zero net velocity, but entails arbitrarily fast diffusion of the dimer for sufficiently weak noise. Moreover, upon application of a static bias force, the dimer exhibits a motion opposite to that force (absolute negative mobility). The key requirement for all these effects is a nonconvex interaction potential of the two particles.

Cutting ice: nanowire regelation

Even below its normal melting temperature, ice melts when subjected to high pressure and refreezes once the pressure is lifted. A classic demonstration of this regelation phenomenon is the passing of a thin wire through a block of ice when sufficient force is exerted. Here we present a molecular-dynamics study of a nanowire cutting through ice to unravel the molecular level mechanisms responsible for regelation. In particular, we show that the transition from a stationary to a moving wire due to increased driving force changes from symmetric and continuous to asymmetric and discontinuous as a hydrophilic wire is replaced by a hydrophobic one. This is explained at the molecular level in terms of the wetting properties of the wire.

Statistics of static avalanches in a random pinning landscape

We study the minimum-energy configuration of a d -dimensional elastic interface in a random potential tied to a harmonic spring. As a function of the spring position, the center of mass of the interface changes in discrete jumps, also called shocks or "static avalanches." We obtain analytically the distribution of avalanche sizes and its cumulants within an =4-d expansion from a tree and one-loop resummation using functional renormalization. This is compared with exact numerical minimizations of interface energies for random-field disorder in d=2,3 . Connections to dynamic avalanches are mentioned.

Size distributions of shocks and static avalanches from the functional renormalization group

Interfaces pinned by quenched disorder are often used to model jerky self-organized critical motion. We study static avalanches, or shocks, defined here as jumps between distinct global minima upon changing an external field. We show how the full statistics of these jumps is encoded in the functional-renormalization-group fixed-point functions. This allows us to obtain the size distribution P(S) of static avalanches in an expansion in the internal dimension d of the interface. Near and above d=4 this yields the mean-field distribution P(S) approximately S;{-3/2}e;{-S4S_{m}} , where S_{m} is a large-scale cutoff, in some cases calculable. Resumming all one-loop contributions, we find P(S) approximately S;{-tau}exp(C(SS_{m});{1/2}-B/4(S/S_{m});{delta}) , where B , C , delta , and tau are obtained to first order in =4-d . Our result is consistent to O() with the relation tau=tau_{zeta}:=2-2/d+zeta , where zeta is the static roughness exponent, often conjectured to hold at depinning. Our calculation applies to all static universality classes, including random-bond, random-field, and random-periodic disorders. Extended to long-range elastic systems, it yields a different size distribution for the case of contact-line elasticity, with an exponent compatible with tau=2-1/d+zeta to O(=2-d) . We discuss consequences for avalanches at depinning and for sandpile models, relations to Burgers turbulence and the possibility that the relation tau=tau_{zeta} be violated to higher loop order. Finally, we show that the avalanche-size distribution on a hyperplane of codimension one is in mean field (valid close to and above d=4 ) given by P(S) approximately K_{13}(S)S , where K is the Bessel- K function, thus tau_{hyperplane}=4/3 .

93Nb NMR spin echo spectroscopy in single crystal NbSe3

We report electric field induced phase displacements of the charge density wave (CDW) in a single crystal of NbSe3 using 93Nb NMR spin-echo spectroscopy. CDW polarizations in the pinned state induced by unipolar and bipolar pulses are linear and reversible up to at least E = (0.96)ET. The polarizations have a broad distribution extending up to phase angles of order 60 degrees for electric fields close to threshold. No evidence for polarizations in excess of a CDW wavelength or for a divergence in polarization near ET are observed. The results are consistent with elastic depinning models, provided that the critical regime expected in large systems is not observable.

Spin precession and avalanches

In many magnetic materials, spin dynamics at short times are dominated by precessional motion as damping is relatively small. We describe how avalanches evolve under these conditions. The growth front is spread out over a large region and consists of rapidly fluctuating spins often above the ferromagnetic transition temperature. In the limit of no damping the system will transition to an ergodic state if the initial instability is large enough, but otherwise can die out. This dynamic nucleation phenomenon is analyzed theoretically and the implications for real materials are discussed.

Dynamics below the depinning threshold in disordered elastic systems

We study the steady-state low-temperature dynamics of an elastic line in a disordered medium below the depinning threshold. Analogously to the equilibrium dynamics, in the limit T-->0, the steady state is dominated by a single configuration which is occupied with probability 1. We develop an exact algorithm to target this dominant configuration and to analyze its geometrical properties as a function of the driving force. The roughness exponent of the line at large scales is identical to the one at depinning. No length scale diverges in the steady-state regime as the depinning threshold is approached from below. We do find a divergent length, but it is associated only with the transient relaxation between metastable states.

Universal statistics of the critical depinning force of elastic systems in random media

We study the rescaled probability distribution of the critical depinning force of an elastic system in a random medium. We put in evidence the underlying connection between the critical properties of the depinning transition and the extreme value statistics of correlated variables. The distribution is Gaussian for all periodic systems, while in the case of random manifolds there exists a family of universal functions ranging from the Gaussian to the Gumbel distribution. Both of these scenarios are a priori experimentally accessible in finite, macroscopic, disordered elastic systems.

Return to return point memory

We describe a new class of systems exhibiting return point memory (RPM), different from those discussed before in the context of ferromagnets. We show numerically that one-dimensional random Ising antiferromagnets have exact RPM when evolving from a large field, but not when started at finite field, unlike the ferromagnetic case. This implies that the standard approach to understanding ferromagnetic RPM will fail for this case. We also demonstrate RPM with a set of variables that keeps track of spin flips at each site. Conventional RPM for the spins is a projection of this result, suggesting that spin flip variables might be a more fundamental representation of the dynamics. We also present a mapping that embeds the antiferromagnetic chain in a two-dimensional ferromagnet, and prove RPM for spin-exchange dynamics in the interior of the chain with this mapping.

Crackling noise

Crackling noise arises when a system responds to changing external conditions through discrete, impulsive events spanning a broad range of sizes. A wide variety of physical systems exhibiting crackling noise have been studied, from earthquakes on faults to paper crumpling. Because these systems exhibit regular behaviour over a huge range of sizes, their behaviour is likely to be independent of microscopic and macroscopic details, and progress can be made by the use of simple models. The fact that these models and real systems can share the same behaviour on many scales is called universality. We illustrate these ideas by using results for our model of crackling noise in magnets, explaining the use of the renormalization group and scaling collapses, and we highlight some continuing challenges in this still-evolving field.

Mirror symmetry breaking through an internal degree of freedom leading to directional motion

We analyze here the minimal conditions for directional motion (net flow in phase space) of a molecular motor placed on a mirror-symmetric environment and driven by a center-symmetric and time-periodic force field. The complete characterization of the deterministic limit of the dissipative dynamics of several realizations of this minimal model reveals a complex structure in the phase diagram in parameter space, with intertwined regions of pinning (closed orbits) and directional motion. This demonstrates that the mirror symmetry breaking, which is needed for directional motion to occur, can operate through an internal degree of freedom coupled to the translational one.