Modeling future cliff-front waves during sea level rise and implications for coastal cliff retreat rates.

It is often assumed that future coastal cliff retreat rates will accelerate as global sea level rises, but few studies have investigated how SLR (sea level rise) might change cliff-front wave dynamics. Using a new simple numerical model, this study simulates the number and type (breaking, broken, or unbroken) of cliff-front waves under future SLR scenarios. Previous research shows breaking waves deliver more energy to cliffs than broken waves, and unbroken waves generate minimal impact. Here, we investigated six cliff-platform profiles from three regions (USA, New Zealand, and UK) with varied tidal ranges and wave climates. Model inputs included 2013-2100 hindcast/forecast incident wave height and tidal water level, and three future SLR scenarios. Results show the number of both cliff-front breaking and broken waves generally increase for a high-elevation (relative to tide) cliff-platform junction. In contrast, breaking/broken wave occurrence decrease by 38-92% for a near-horizontal shore platform with a low-elevation cliff-platform junction under a high SRL scenario, leading to high (96-97%) unbroken wave occurrence. Overall, results suggest the response of cliff-front waves to future SLR is complex and depends on shore platform geometries and SLR scenarios, indicating that future cliff retreat rates may not homogeneously accelerate under SLR.

Metastate analysis of the ground states of two-dimensional Ising spin glasses.

Using an efficient polynomial-time ground-state algorithm we investigate the Ising spin glass state at zero temperature in two dimensions. For large sizes, we show that the spin state in a central region is independent of the interactions far away, indicating a "single-state" picture, presumably the droplet model. Surprisingly, a single power law describes corrections to this result down to the smallest sizes studied.

Nontrivial maturation metastate-average state in a one-dimensional long-range Ising spin glass: Above and below the upper critical range.

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}∝1/R_{ij}^{2σ}, and we study a range of values of σ with 1/2<σ<1. We consider a correlation function C_{4}(r,t). A dynamic correlation length that shows power-law growth with time ξ(t)∝t^{1/z} can be identified in the data and, for large time t, C_{4}(r,t) decays as a power law r^{-α_{d}} with distance r when r≪ξ(t). The calculation can be interpreted in terms of the maturation metastate averaged Gibbs state, or MMAS, and the decay exponent α_{d} differentiates between a trivial MMAS (α_{d}=0), as expected in the droplet picture of spin glasses, and a nontrivial MMAS (α_{d}≠0), as in the replica-symmetry-breaking (RSB) or chaotic pairs pictures. We find nonzero α_{d} even in the regime σ>2/3 which corresponds to short-range systems below six dimensions. For σ<2/3, the decay exponent α_{d} follows the RSB prediction for the decay exponent α_{s}=3-4σ of the static metastate, consistent with a conjectured statics-dynamics relation, while it approaches α_{d}=1-σ in the regime 2/3<σ<1; however, it deviates from both lines in the vicinity of σ=2/3.

Ordering behavior of the two-dimensional Ising spin glass with long-range correlated disorder.

The standard short-range two-dimensional Ising spin glass is numerically well accessible, in particular, because there are polynomial-time ground-state algorithms. On the other hand, in contrast to higher dimensional spin glasses, it does not exhibit a rich behavior, i.e., no ordered phase at finite temperature. Here, we investigate whether long-range correlated bonds change this behavior. This would still keep the model numerically well accessible while exhibiting a more interesting behavior. The bonds are drawn from a Gaussian distribution with a two-point correlation for bonds at distance r that decays as (1+r^{2})^{-a/2}, a≥0. We study numerically with exact algorithms the ground-state and domain-wall excitations. Our results indicate that the inclusion of bond correlations still does not lead to a spin-glass order at any finite temperature. A further analysis reveals that bond correlations have a strong effect at local length scales, inducing ferro- and antiferromagnetic domains into the system. The length scale of ferro- and antiferromagnetic order diverges exponentially as the correlation exponent approaches a critical value, a→a_{crit}=0. Thus, our results suggest that the system becomes a ferro- or antiferromagnet only in the limit a→0.

Percolation of Fortuin-Kasteleyn clusters for the random-bond Ising model.

We apply generalizations of the Swendson-Wang and Wolff cluster algorithms, which are based on the construction of Fortuin-Kasteleyn clusters, to the three-dimensional ±1 random-bond Ising model. The behavior of the model is determined by the temperature T and the concentration p of negative (antiferromagnetic) bonds. The ground state is ferromagnetic for 0≤p0, our data suggest that the percolation transition is universal, irrespective of whether the ground state exhibits ferromagnetic or spin-glass order, and is in the universality class of standard percolation. This shows that correlations in the bond occupancy of the Fortuin-Kasteleyn clusters are irrelevant, except for p=0 where the clusters are strictly tied to Ising correlations so the percolation transition is in the Ising universality class.

Critical and Griffiths-McCoy singularities in quantum Ising spin glasses on d-dimensional hypercubic lattices: A series expansion study.

We study the ±J transverse-field Ising spin-glass model at zero temperature on d-dimensional hypercubic lattices and in the Sherrington-Kirkpatrick (SK) model, by series expansions around the strong-field limit. In the SK model and in high dimensions our calculated critical properties are in excellent agreement with the exact mean-field results, surprisingly even down to dimension d=6, which is below the upper critical dimension of d=8. In contrast, at lower dimensions we find a rich singular behavior consisting of critical and Griffiths-McCoy singularities. The divergence of the equal-time structure factor allows us to locate the critical coupling where the correlation length diverges, implying the onset of a thermodynamic phase transition. We find that the spin-glass susceptibility as well as various power moments of the local susceptibility become singular in the paramagnetic phase before the critical point. Griffiths-McCoy singularities are very strong in two dimensions but decrease rapidly as the dimension increases. We present evidence that high enough powers of the local susceptibility may become singular at the pure-system critical point.

Stability of the quantum Sherrington-Kirkpatrick spin glass model.

I study in detail the quantum Sherrington-Kirkpatrick (SK) model, i.e., the infinite-range Ising spin glass in a transverse field, by solving numerically the effective one-dimensional model that the quantum SK model can be mapped to in the thermodynamic limit. I find that the replica symmetric solution is unstable down to zero temperature, in contrast to some previous claims, and so there is not only a line of transitions in the (longitudinal) field-temperature plane (the de Almeida-Thouless, AT, line) where replica symmetry is broken, but also a quantum de Almeida-Thouless (QuAT) line in the transverse field-longitudinal field plane at T=0. If the QuAT line also occurs in models with short-range interactions its presence might affect the performance of quantum annealers when solving spin glass-type problems with a bias (i.e., magnetic field).

de Almeida-Thouless instability in short-range Ising spin glasses.

We use high-temperature series expansions to study the ±J Ising spin glass in a magnetic field in d-dimensional hypercubic lattices for d=5-8 and in the infinite-range Sherrington-Kirkpatrick (SK) model. The expansions are obtained in the variable w=tanh^{2}J/T for arbitrary values of u=tanh^{2}h/T complete to order w^{10}. We find that the scaling dimension Δ associated with the ordering-field h^{2} equals 2 in the SK model and for d≥6. However, in agreement with the work of Fisher and Sompolinsky [Phys. Rev. Lett. 54, 1063 (1985)PRLTAO0031-900710.1103/PhysRevLett.54.1063], there is a violation of scaling in a finite field, leading to an anomalous h-T dependence of the de Almeida-Thouless (AT) [J. Phys. A 11, 983 (1978)JPHAC50305-447010.1088/0305-4470/11/5/028] line in high dimensions, whereas scaling is restored as d→6. Within the convergence of our series analysis, we present evidence supporting an AT line in d≥6. In d=5, the exponents γ and Δ are substantially larger than mean-field values, but we do not see clear evidence for the AT line in d=5.

Publisher's Note: Critical and Griffiths-McCoy singularities in quantum Ising spin glasses on d-dimensional hypercubic lattices: A series expansion study [Phys. Rev. E 96, 022139 (2017)].

This corrects the article DOI: 10.1103/PhysRevE.96.022139.

Spin glass behavior in a random Coulomb antiferromagnet.

We study spin glass behavior in a random Ising Coulomb antiferromagnet in two and three dimensions using Monte Carlo simulations. In two dimensions, we find a transition at zero temperature with critical exponents consistent with those of the Edwards-Anderson model, though with large uncertainties. In three dimensions, evidence for a finite-temperature transition, as occurs in the Edwards-Anderson model, is rather weak. This may indicate that the sizes are too small to probe the asymptotic critical behavior, or possibly that the universality class is different from that of the Edwards-Anderson model and has a lower critical dimension equal to three.

Universal dynamic scaling in three-dimensional Ising spin glasses.

We use a nonequilibrium Monte Carlo simulation method and dynamical scaling to study the phase transition in three-dimensional Ising spin glasses. The transition point is repeatedly approached at finite velocity v (temperature change versus time) in Monte Carlo simulations starting at a high temperature. This approach has the advantage that the equilibrium limit does not have to be strictly reached for a scaling analysis to yield critical exponents. For the dynamic exponent we obtain z=5.85(9) for bimodal couplings distribution and z=6.00(10) for the Gaussian case. Assuming universal dynamic scaling, we combine the two results and obtain z=5.93±0.07 for generic 3D Ising spin glasses.

Finite-size scaling above the upper critical dimension.

We present a unified view of finite-size scaling (FSS) in dimension d above the upper critical dimension, for both free and periodic boundary conditions. We find that the modified FSS proposed some time ago to allow for violation of hyperscaling due to a dangerous irrelevant variable applies only to k=0 fluctuations, and "standard" FSS applies to k≠0 fluctuations. Hence the exponent η describing power-law decay of correlations at criticality is unambiguously η=0. With free boundary conditions, the finite-size "shift" is greater than the rounding. Nonetheless, using T-T(L), where T(L) is the finite-size pseudocritical temperature, rather than T-T(c), as the scaling variable, the data do collapse onto a scaling form that includes the behavior both at T(L), where the susceptibility χ diverges like L(d/2), and at the bulk T(c), where it diverges like L(2). These claims are supported by large-scale simulations on the five-dimensional Ising model.

Spin glasses in the nonextensive regime.

Spin systems with long-range interactions are "nonextensive" if the strength of the interactions falls off sufficiently slowly with distance. It has been conjectured for ferromagnets and, more recently, for spin glasses that, everywhere in the nonextensive regime, the free energy is exactly equal to that for the infinite range model in which the characteristic strength of the interaction is independent of distance. In this paper we present the results of Monte Carlo simulations of the one-dimensional long-range spin glasses in the nonextensive regime. Using finite-size scaling, our results for the transition temperatures are consistent with this prediction. We also propose and provide numerical evidence for an analogous result for dilutedlong-range spin glasses in which the coordination number is finite, namely, that the transition temperature throughout the nonextensive regime is equal to that of the infinite-range model known as the Viana-Bray model.

Complexity of several constraint-satisfaction problems using the heuristic classical algorithm WalkSAT.

We determine the complexity of several constraint-satisfaction problems using the heuristic algorithm WalkSAT. At large sizes N, the complexity increases exponentially with N in all cases. Perhaps surprisingly, out of all the models studied, the hardest for WalkSAT is the one for which there is a polynomial time algorithm.

Exponential complexity of the quantum adiabatic algorithm for certain satisfiability problems.

We determine the complexity of several constraint satisfaction problems using the quantum adiabatic algorithm in its simplest implementation. We do so by studying the size dependence of the gap to the first excited state of "typical" instances. We find that, at large sizes N, the complexity increases exponentially for all models that we study. We also compare our results against the complexity of the analogous classical algorithm WalkSAT and show that the harder the problem is for the classical algorithm, the harder it is also for the quantum adiabatic algorithm.

de Almeida-Thouless line in vector spin glasses.

We consider the infinite-range spin glass in which the spins have m>1 components (a vector spin glass). Applying a magnetic field which is random in direction, there is a de Almeida-Thouless (AT) line below which the "replica symmetric" solution is unstable, just as for the Ising (m=1) case. We calculate the location of this AT line for Gaussian random fields for arbitrary m and verify our results by numerical simulations for m=3 .

First-order phase transition in the quantum adiabatic algorithm.

We simulate the quantum adiabatic algorithm (QAA) for the exact cover problem for sizes up to N=256 using quantum Monte Carlo simulations incorporating parallel tempering. At large N, we find that some instances have a discontinuous (first-order) quantum phase transition during the evolution of the QAA. This fraction increases with increasing N and may tend to 1 for N-->infinity.

Continuum and lattice heat currents for oscillator chains.

We show that two commonly used definitions for the heat current give different results-through the Kubo formula-for the heat conductivity of oscillator chains. The difference exists for finite chains, and is expected to be important more generally for small structures. For a chain of N particles that are tethered at the ends, the ratio of the heat conductivities calculated with the two currents differs from unity by O(1/N). For a chain held at constant pressure, the difference from unity decays more slowly, and is consistent with O(1/Neta) with 1>eta>0.5.

Study of the de Almeida-Thouless line using power-law diluted one-dimensional Ising spin glasses.

We test for the existence of a spin-glass phase transition, the de Almeida-Thouless line, in an externally applied (random) magnetic field by performing Monte Carlo simulations on a power-law diluted one-dimensional Ising spin glass for very large system sizes. We find that a de Almeida-Thouless line occurs only in the mean-field regime, which corresponds, for a short-range spin glass, to dimension d larger than 6.