Stretched exponential relaxation in the Coulomb glass

The relaxation of the specific heat and the entropy to their equilibrium values is investigated numerically for the three-dimensional Coulomb glass at very low temperatures. The long time relaxation follows a stretched exponential function, f (t) = f ₀exp - (t/τ)^β , with the exponent β increasing with the temperature. The relaxation time diverges as an Arrhenius law when T→ 0. Received 24 May 2001 and Received in final form 12 September 2001     

Multifractal properties of aperiodic Ising model on hierarchical lattices: role of the geometric fluctuations

The role of the geometric fluctuations on the multifractal properties of the local magnetization of aperiodic ferromagnetic Ising models on hierarchical lattices is investigated. The geometric fluctuations are introduced by generalized Fibonacci sequences. The local magnetization is evaluated via an exact recurrent procedure encompassing real space renormalization group decimation. The symmetries of the local magnetization patterns induced by the aperiodic couplings is found to be strongly (weakly) different, with respect to the ones of the corresponding homogeneous systems, when the geometric fluctuations are relevant (irrelevant) to change the critical properties of the system. At the criticality, the measure defined by the local magnetization is found to exhibit a non-trivial F(α) spectra being shifted to higher values of α when relevant geometric fluctuations are considered. The critical exponents are found to be related with some special points of the F(α) function and agree with previous results obtained by the quite distinct transfer matrix approach. Received 2 April 2001 and Received in final form 14 August 2001     

Finite-size scaling in systems with long-range interaction

The finite-size critical properties of the (n) vector ϕ⁴ model, with long-range interaction decaying algebraically with the interparticle distance r like r ^{-d - σ}, are investigated. The system is confined to a finite geometry subject to periodic boundary condition. Special attention is paid to the finite-size correction to the bulk susceptibility above the critical temperature T _c. We show that this correction has a power-law nature in the case of pure long-range interaction i.e. 0 < σ < 2 and it turns out to be exponential in case of short-range interaction i.e.σ = 2. The results are valid for arbitrary dimension d, between the lower ( d _< = σ) and the upper ( d _> = 2σ) critical dimensions. Received 2 July 2001 and Received in final form 4 Septembre 2001     

Core percolation in random graphs: a critical phenomena analysis

We study both numerically and analytically what happens to a random graph of average connectivity α when its leaves and their neighbors are removed iteratively up to the point when no leaf remains. The remnant is made of isolated vertices plus an induced subgraph we call the core. In the thermodynamic limit of an infinite random graph, we compute analytically the dynamics of leaf removal, the number of isolated vertices and the number of vertices and edges in the core. We show that a second order phase transition occurs at α = e = 2.718 ... : below the transition, the core is small but above the transition, it occupies a finite fraction of the initial graph. The finite size scaling properties are then studied numerically in detail in the critical region, and we propose a consistent set of critical exponents, which does not coincide with the set of standard percolation exponents for this model. We clarify several aspects in combinatorial optimization and spectral properties of the adjacency matrix of random graphs. Received 31 January 2001 and Received in final form 26 June 2001     

Random quantum magnets with broad disorder distribution

We study the critical behavior of Ising quantum magnets with broadly distributed random couplings (J), such that P(ln J) ∼ | ln J|^{-1 - α}, α > 1, for large | ln J| (Lévy flight statistics). For sufficiently broad distributions, α < , the critical behavior is controlled by a line of fixed points, where the critical exponents vary with the Lévy index, α. In one dimension, with = 2, we obtained several exact results through a mapping to surviving Riemann walks. In two dimensions the varying critical exponents have been calculated by a numerical implementation of the Ma-Dasgupta-Hu renormalization group method leading to ≈ 4.5. Thus in the region 2 < α < , where the central limit theorem holds for | ln J| the broadness of the distribution is relevant for the 2d quantum Ising model. Received 6 December 2000 and Received in final form 22 January 2001     

Dynamics of a single ion in a perturbed Penning trap. Sextupolar perturbation

In the frame work of classical mechanics, we study the nonlinear dynamics of a single ion trapped in a Penning trap perturbed by an electrostatic sextupolar perturbation. The perturbation is caused by a deformation in the configuration of the electrodes. By using a Hamiltonian formulation, we obtain that the system is governed by three parameters: the z-component of the canonical angular momentum P _φ - which is a constant of the motion because the perturbation we assume is axial-symmetric -, the parameter δ that determines the ratio between the axial and the cyclotron frequencies, and the parameter a which indicates how far from the ideal design the electrodes are. We study the case P _φ = 0. By means of surfaces of section, we show that the phase space structure is made of three fundamental families of orbits: arch, loop and box orbits. The coexistence of these kinds of orbits depends on the parameter δ. The escape is also explained on the basis of the shape of the potential energy surface as well as of the phase space structure. Received 6 September 2001 / Received in final form 19 March 2002 Published online 28 June 2002     

Failure time and critical behaviour of fracture precursors in heterogeneous materials

The acoustic emission of fracture precursors, and the failure time of samples of heterogeneous materials (wood, fiberglass) are studied as a function of the load features and geometry. It is shown that in these materials the failure time is predicted with a good accuracy by a model of microcrack nucleation proposed by Pomeau. We find that the time interval δt between events (precursors) and the energy ɛ are power law distributed and that the exponents of these power laws depend on the load history and on the material. In contrast, the cumulated acoustic energy E presents a critical divergency near the breaking time τ which is E∼ . The positive exponent γ is independent, within error bars, on all the experimental parameters. Received 31 July 2001 and Received in final form 17 December 2001     

Renormalization group treatment of the scaling properties of finite systems with subleading long-range interaction

The finite size behavior of the susceptibility, Binder cumulant and some even moments of the magnetization of a fully finite O(n) cubic system of size L are analyzed and the corresponding scaling functions are derived within a field-theoretic ɛ-expansion scheme under periodic boundary conditions. We suppose a van der Waals type long-range interaction falling apart with the distance r as r ^{- (d + σ)}, where 2 < σ < 4, which does not change the short-range critical exponents of the system. Despite that the system belongs to the short-range universality class it is shown that above the bulk critical temperature T _c the finite-size corrections decay in a power-in-L, and not in an exponential-in-L law, which is normally believed to be a characteristic feature for such systems. Received 8 August 2001     

Universality in diffusion front growth dynamics

We have studied the scaling properties of diffusion fronts by numerical calculations based on the mean field approach in the context of a lattice gas model, performed in a triangular lattice. We find that the height-height correlation function scales with time t and length l as C(l, t) ≈l ^α f (t/l ^α/β) with α = 0.62±0.01 and β = 0.39±0.02. These exponent values are identical to those characterising the roughness of the diffusion fronts evolving through a square lattice [1,2], thus confirming their universality. Received 14 November 2001 / Received in final form 20 April 2002 Published online 31 July 2002     

Dynamical properties of the slithering-snake algorithm: A numerical test of the activated-reptation hypothesis

Correlations in the motion of reptating polymers in a melt are investigated by means of Monte Carlo simulations of the three-dimensional slithering-snake version of the bond-fluctuation model. Surprisingly, the slithering-snake dynamics becomes inconsistent with classical reptation predictions at high chain overlap (created either by chain length N or by the volume fraction φ of occupied lattice sites), where the relaxation times increase much faster than expected. This is due to the anomalous curvilinear diffusion in a finite time window whose upper bound (N) is set by the density of chain ends φ/N. Density fluctuations created by passing chain ends allow a reference polymer to break out of the local cage of immobile obstacles created by neighboring chains. The dynamics of dense solutions of “snakes” at t ≪ is identical to that of a benchmark system where all chains but one are frozen. We demonstrate that the subdiffusive dynamical regime is caused by the slow creeping of a chain out of its correlation hole. Our results are in good qualitative agreement with the activated-reptation scheme proposed recently by Semenov and Rubinstein (Eur. Phys. J. B, 1 (1998) 87). Additionally, we briefly comment on the relevance of local relaxation pathways within a slithering-snake scheme. Our preliminary results suggest that a judicious choice of the ratio of local to slithering-snake moves is crucial to equilibrate a melt of long chains efficiently. Received: 18 December 2002 / Accepted: 3 April 2003 / Published online: 12 May 2003 RID="a" ID="a"e-mail: jwittmer@dpm.univ-lyon1.fr RID="b" ID="b"Current address: University of Illinois at Urbana-Champaign.     

Short-range plasma model for intermediate spectral statistics

We propose a plasma model for spectral statistics displaying level repulsion without long-range spectral rigidity, i.e. statistics intermediate between random matrix and Poisson statistics similar to the ones found numerically at the critical point of the Anderson metal-insulator transition in disordered systems and in certain dynamical systems. The model emerges from Dysons one-dimensional gas corresponding to the eigenvalue distribution of the classical random matrix ensembles by restricting the logarithmic pair interaction to a finite number k of nearest neighbors. We calculate analytically the spacing distributions and the two-level statistics. In particular we show that the number variance has the asymptotic form Σ²(L) ∼χL for large L and the nearest-neighbor distribution decreases exponentially when s→∞, P(s) ∼ exp(- Λs) with Λ = 1/χ = kβ + 1, where β is the inverse temperature of the gas (β = 1, 2 and 4 for the orthogonal, unitary and symplectic symmetry class respectively). In the simplest case of k = β = 1, the model leads to the so-called Semi-Poisson statistics characterized by particular simple correlation functions e.g. P(s) = 4s exp(- 2s). Furthermore we investigate the spectral statistics of several pseudointegrable quantum billiards numerically and compare them to the Semi-Poisson statistics. Received 13 September 2000     

Zeroth and second laws of thermodynamics simultaneously questioned in the quantum microworld

A new and rather trivial model is suggested with mechanism that implies simultaneous violation of the zeroth and the second laws of thermodynamics. Mathematically rigorous quantum theory reduces to a trivial application of the Golden rule formula. It yields exciton on-energy-shell diffusion caused by bath-nonassisted excitation hopping between tails of different exciton site levels ε₁ < ε₂ broadened by bath-assisted finite life-time effects. The elastic character of the hopping implies 1 ↔ 2-symmetric transfer rate W. Thus the net diffusion exciton flow W(P ₁ - P ₂) and also, as argued, the net energy flow are possible due to different near-to-equilibrium exciton populations P ₁ > P ₂. As the sites are provided with two different baths, the population imbalance and the flows survive even for slightly different local bath temperatures T ₁ < T ₂ < T ₁ε₂/ε₁. Thus spontaneous exciton and also energy flows against temperature step become possible, in contradiction with the Clausius form of the second law. Violations of both the laws disappear in the high-temperature, i.e. classical limit Received 16 May 2001 and Received in final form 20 September 2001     

Energy landscapes, lowest gaps, and susceptibility of elastic manifolds at zero temperature

We study the effect of an external field on (1 + 1) and (2 + 1) dimensional elastic manifolds, at zero temperature and with random bond disorder. Due to the glassy energy landscape the configuration of a manifold changes often in abrupt, “first order”-type of large jumps when the field is applied. First the scaling behavior of the energy gap between the global energy minimum and the next lowest minimum of the manifold is considered, by employing exact ground state calculations and an extreme statistics argument. The scaling has a logarithmic prefactor originating from the number of the minima in the landscape, and reads ΔE ₁∼L ^θ[ln(L _z L ^{- ζ})]^-1/2, where ζ is the roughness exponent and θ is the energy fluctuation exponent of the manifold, L is the linear size of the manifold, and L_z is the system height. The gap scaling is extended to the case of a finite external field and yields for the susceptibility of the manifolds ∼L ^{2D + 1 - θ}[(1 - ζ)ln(L)]^1/2. We also present a mean field argument for the finite size scaling of the first jump field, h ₁∼L ^{d - θ}. The implications to wetting in random systems, to finite-temperature behavior and the relation to Kardar-Parisi-Zhang non-equilibrium surface growth are discussed. Received December 2000 and Received in final form April 2001     

Bulk singularities at critical end points: a field-theory analysis

A class of continuum models with a critical end point is considered whose Hamiltonian [φ,ψ] involves two densities: a primary order-parameter field, φ, and a secondary (noncritical) one, ψ. Field-theoretic methods (renormalization group results in conjunction with functional methods) are used to give a systematic derivation of singularities occurring at critical end points. Specifically, the thermal singularity ∼ | t|^{2 - α} of the first-order line on which the disordered or ordered phase coexists with the noncritical spectator phase, and the coexistence singularity ∼ | t|^{1 - α} or ∼ | t|^β of the secondary density <ψ> are derived. It is clarified how the renormalization group (RG) scenario found in position-space RG calculations, in which the critical end point and the critical line are mapped onto two separate fixed points _CEP ^* and _λ ^*, translates into field theory. The critical RG eigenexponents of _CEP ^* and _λ ^* are shown to match. _CEP ^* is demonstrated to have a discontinuity eigenperturbation (with eigenvalue y = d), tangent to the unstable trajectory that emanates from _CEP ^* and leads to _λ ^*. The nature and origin of this eigenperturbation as well as the role redundant operators play are elucidated. The results validate that the critical behavior at the end point is the same as on the critical line. Received 18 January 2001     

Universality of a dynamical percolative approach to 1/f γ noise

A dynamical percolative model explaining the universality of 1/ f ^γ noise is reported. Exponents γ ranging from 0 to 2 are obtained under the hypothesis that noise originates from random switching events between two ON-OFF states in elemental parts (switchers) of a physical system. The usual noise behaviour with γ very close to 1 in an arbitrarily wide frequency range is obtained assuming a statistical distribution of switcher relaxation time τ proportional to τ ^-1 , as in McWhorter's model. The impact of these results with respect to recent self-organised criticality models is discussed. Received 6 November 2000 and Received in final form 22 May 2001     

Dynamical effects of a one-dimensional multibarrier potential of finite range

We discuss the properties of a large number N of one-dimensional (bounded) locally periodic potential barriers in a finite interval. We show that the transmission coefficient, the scattering cross section σ, and the resonances of σ depend sensitively upon the ratio of the total spacing to the total barrier width. We also show that a time dependent wave packet passing through the system of potential barriers rapidly spreads and deforms, a criterion suggested by Zaslavsky for chaotic behaviour. Computing the spectrum by imposing (large) periodic boundary conditions we find a Wigner type distribution. We investigate also the S-matrix poles; many resonances occur for certain values of the relative spacing between the barriers in the potential. Received 1st August 2001 and Received in final form 18 November 2001     

Cold ignition of combustion-like waves of cryo-chemical reactions in solids

Fronts of weakly exothermal chemical reaction may propagate in solids at very low temperatures ( 4K≤T≤77K) thanks to a quite unusual mechanism, involving a feedback between the heat produced by the reaction and the disruption of the solid matrix. In this class of phenomena, the reaction may be induced by mechanical constraints, without a large elevation of temperature. On the basis of a simple phenomenological model, we investigate ignition of a propagating front by initially (i) disrupting a localized zone of the solid matrix, or by (ii) introducing a temperature jump, leading to a thermal shock with strong temperature gradients. In particular, we show that reaction can be initiated by disrupting only a very small fraction of the sample. Applications to the problem of initiation of solid explosives by friction or shocks is briefly discussed. Received 26 January 2001 and Received in final form 3 May 2001     

Momentum transport and torque scaling in Taylor-Couette flow from an analogy with turbulent convection

We generalize an analogy between rotating and stratified shear flows. This analogy is summarized in Table 1. We use this analogy in the unstable case (centrifugally unstable flow vs. convection) to compute the torque in Taylor-Couette configuration, as a function of the Reynolds number. At low Reynolds numbers, when most of the dissipation comes from the mean flow, we predict that the non-dimensional torque G = T/ν² L, where L is the cylinder length, scales with Reynolds number R and gap width η, G = 1.46η^3/2(1 - η)^-7/4 R ^3/2. At larger Reynolds number, velocity fluctuations become non-negligible in the dissipation. In these regimes, there is no exact power law dependence the torque versus Reynolds. Instead, we obtain logarithmic corrections to the classical ultra-hard (exponent 2) regimes: G = 0.50 . These predictions are found to be in excellent agreement with avail-able experimental data. Predictions for scaling of velocity fluctuations are also provided. Received 7 June 2001 and Received in final form 7 December 2001     

Structure of gels and aggregates of disk-like colloids

The static structure factor (S(q)) of dispersions and gels of disk-like mineral colloids (Laponite) was investigated using time- and ensemble-averaged light scattering. The evolution of S(q) in time after increasing the ionic strength of well-dispersed Laponite suspensions shows that Laponite aggregates and forms fractal clusters. The structure of the aggregates does not depend on the ionic strength, but the rate of growth increases very strongly with the ionic strength. At concentrations below about 3 g/l (0.12% v/v) the aggregates sediment while at higher concentrations space-filling gels are formed. The gels are homogeneous on length scales larger than the correlation length which decreases strongly with decreasing ionic strength and increasing concentration. However, the local structure is the same, independent of the concentration and the ionic strength. Received 6 August 2000 and Received in final form 16 March 2001     

Directed spiral site percolation on the square lattice

A new site percolation model, directed spiral percolation (DSP), under both directional and rotational (spiral) constraints is studied numerically on the square lattice. The critical percolation threshold p _c ≈ 0.655 is found between the directed and spiral percolation thresholds. Infinite percolation clusters are fractals of dimension d _f ≈ 1.733. The clusters generated are anisotropic. Due to the rotational constraint, the cluster growth is deviated from that expected due to the directional constraint. Connectivity lengths, one along the elongation of the cluster and the other perpendicular to it, diverge as p → p _c with different critical exponents. The clusters are less anisotropic than the directed percolation clusters. Different moments of the cluster size distribution P _s(p) show power law behaviour with | p - p _c| in the critical regime with appropriate critical exponents. The values of the critical exponents are estimated and found to be very different from those obtained in other percolation models. The proposed DSP model thus belongs to a new universality class. A scaling theory has been developed for the cluster related quantities. The critical exponents satisfy the scaling relations including the hyperscaling which is violated in directed percolation. A reasonable data collapse is observed in favour of the assumed scaling function form of P _s(p). The results obtained are in good agreement with other model calculations. Received 10 November 2002 / Received in final form 20 February 2003 Published online 23 May 2003 RID="a" ID="a"e-mail: santra@iitg.ernet.in