No Arabic abstract
In this paper we study bond percolation on a one-dimensional chain with power-law bond probability $C/ r^{1+sigma}$, where $r$ is the distance length between distinct sites. We introduce and test an order $N$ Monte Carlo algorithm and we determine as a function of $sigma$ the critical value $C_{c}$ at which percolation occurs. The critical exponents in the range $0<sigma<1$ are reported and compared with mean-field and $varepsilon$-expansion results. Our analysis is in agreement, up to a numerical precision $approx 10^{-3}$, with the mean field result for the anomalous dimension $eta=2-sigma$, showing that there is no correction to $eta$ due to correlation effects.
In this note we study metastability phenomena for a class of long-range Ising models in one-dimension. We prove that, under suitable general conditions, the configuration -1 is the only metastable state and we estimate the mean exit time. Moreover, we illustrate the theory with two examples (exponentially and polynomially decaying interaction) and we show that the critical droplet can be macroscopic or mesoscopic, according to the value of the external magnetic field.
Understanding the low-temperature pure state structure of spin glasses remains an open problem in the field of statistical mechanics of disordered systems. Here we study Monte Carlo dynamics, performing simulations of the growth of correlations following a quench from infinite temperature to a temperature well below the spin-glass transition temperature $T_c$ for a one-dimensional Ising spin glass model with diluted long-range interactions. In this model, the probability $P_{ij}$ that an edge ${i,j}$ has nonvanishing interaction falls as a power-law with chord distance, $P_{ij}propto1/R_{ij}^{2sigma}$, and we study a range of values of $sigma$ with $1/2<sigma<1$. We consider a correlation function $C_{4}(r,t)$. A dynamic correlation length that shows power-law growth with time $xi(t)propto t^{1/z}$ can be identified in the data and, for large time $t$, $C_{4}(r,t)$ decays as a power law $r^{-alpha_d}$ with distance $r$ when $rll xi(t)$. The calculation can be interpreted in terms of the maturation metastate averaged Gibbs state, or MMAS, and the decay exponent $alpha_d$ differentiates between a trivial MMAS ($alpha_d=0$), as expected in the droplet picture of spin glasses, and a nontrivial MMAS ($alpha_d e 0$), as in the replica-symmetry-breaking (RSB) or chaotic pairs pictures. We find nonzero $alpha_d$ even in the regime $sigma >2/3$ which corresponds to short-range systems below six dimensions. For $sigma < 2/3$, the decay exponent $alpha_d$ follows the RSB prediction for the decay exponent $alpha_s = 3 - 4 sigma$ of the static metastate, consistent with a conjectured statics-dynamics relation, while it approaches $alpha_d=1-sigma$ in the regime $2/3<sigma<1$; however, it deviates from both lines in the vicinity of $sigma=2/3$.
Finite-$N$ effects unavoidably drive the long-term evolution of long-range interacting $N$-body systems. The Balescu-Lenard kinetic equation generically describes this process sourced by ${1/N}$ effects but this kinetic operator exactly vanishes by symmetry for one-dimensional homogeneous systems: such systems undergo a kinetic blocking and cannot relax as a whole at this order in ${1/N}$. It is therefore only through the much weaker ${1/N^{2}}$ effects, sourced by three-body correlations, that these systems can relax, leading to a much slower evolution. In the limit where collective effects can be neglected, but for an arbitrary pairwise interaction potential, we derive a closed and explicit kinetic equation describing this very long-term evolution. We show how this kinetic equation satisfies an $H$-theorem while conserving particle number and energy, ensuring the unavoidable relaxation of the system towards the Boltzmann equilibrium distribution. Provided that the interaction is long-range, we also show how this equation cannot suffer from further kinetic blocking, i.e., the ${1/N^{2}}$ dynamics is always effective. Finally, we illustrate how this equation quantitatively matches measurements from direct $N$-body simulations.
We consider near-critical two-dimensional statistical systems at phase coexistence on the half plane with boundary conditions leading to the formation of a droplet separating coexisting phases. General low-energy properties of two-dimensional field theories are used in order to find exact analytic results for one- and two-point correlation functions of both the energy density and order parameter fields. The subleading finite-size corrections are also computed and interpreted within an exact probabilistic picture in which interfacial fluctuations are characterized by the probability density of a Brownian excursion. The analytical results are compared against high-precision Monte Carlo simulations we performed for the specific case of the Ising model. The numerical results are found to be in good agreement with the analytic results in absence of adjustable parameters. The explicit analysis of the closed-form expression for order parameter correlations reveals the long-ranged character of interfacial correlations and their confinement within the interfacial region. The analysis of correlations is then carried out in momentum space through the notion of interface structure factor, which we extend to the case of systems bounded by a flat wall. The presence of the wall and its associated entropic repulsion leads to a specific term in the interface structure factor which we identify.
We investigate the nonequilibrium dynamics following a quench to zero temperature of the non-conserved Ising model with power-law decaying long-range interactions $propto 1/r^{d+sigma}$ in $d=2$ spatial dimensions. The zero-temperature coarsening is always of special interest among nonequilibrium processes, because often peculiar behavior is observed. We provide estimates of the nonequilibrium exponents, viz., the growth exponent $alpha$, the persistence exponent $theta$, and the fractal dimension $d_f$. It is found that the growth exponent $alphaapprox 3/4$ is independent of $sigma$ and different from $alpha=1/2$ as expected for nearest-neighbor models. In the large $sigma$ regime of the tunable interactions only the fractal dimension $d_f$ of the nearest-neighbor Ising model is recovered, while the other exponents differ significantly. For the persistence exponent $theta$ this is a direct consequence of the different growth exponents $alpha$ as can be understood from the relation $d-d_f=theta/alpha$; they just differ by the ratio of the growth exponents $approx 3/2$. This relation has been proposed for annihilation processes and later numerically tested for the $d=2$ nearest-neighbor Ising model. We confirm this relation for all $sigma$ studied, reinforcing its general validity.