Replica symmetry broken states of some glass models.

We have studied in detail the M-p balanced spin-glass model, especially the case p=4. These types of model have relevance to structural glasses. The models possess two kinds of broken replica states; those with one-step replica symmetry breaking (1RSB) and those with full replica symmetry breaking (FRSB). To determine which arises requires studying the Landau expansion to quintic order. There are nine quintic-order coefficients, and five quartic-order coefficients, whose values we determine for this model. We show that it is only for 2≤M<2.4714⋯ that the transition at mean-field level is to a state with FRSB, while for larger M values there is either a continuous transition to a state with 1RSB (when M≤3) or a discontinuous transition for M>3. The Gardner transition from a 1RSB state at low temperatures to a state with FRSB also requires the Landau expansion to be taken to quintic order. Our result for the form of FRSB in the Gardner phase is similar to that found when 2≤M<2.4714⋯, but differs from that given in the early paper of Gross et al. [Phys. Rev. Lett. 55, 304 (1985)0031-900710.1103/PhysRevLett.55.304]. Finally we discuss the effects of fluctuations on our mean-field solutions using the scheme of Höller and Read [Phys. Rev. E 101, 042114 (2020)2470-004510.1103/PhysRevE.101.042114] and argue that such fluctuations will remove both the continuous 1RSB transition and discontinuous 1RSB transitions when 8>d≥6 leaving just the FRSB continuous transition. We suggest values for M and p which might be used in simulations to confirm whether fluctuation corrections do indeed remove the 1RSB transitions.

Study of the de Almeida-Thouless line in the one-dimensional diluted power-law XY spin glass.

We study the de Almeida-Thouless (AT) line in the one-dimensional power-law diluted XY spin-glass model, in which the probability that two spins separated by a distance r interact with each other, decays as 1/r^{2σ}. Tuning the exponent σ is equivalent to changing the space dimension of a short-range model. We develop a heat bath algorithm to equilibrate XY spins; using this in conjunction with the standard parallel tempering and overrelaxation sweeps, we carry out large-scale Monte Carlo simulations. For σ=0.6, which is in the mean-field regime above six dimensions-it is similar to being in 10 dimensions-we find clear evidence for an AT line. For σ=0.75 and σ=0.85, which are in the non-mean-field regime and similar to four and three dimensions, respectively, our data is like that found in previous studies of the Ising and Heisenberg spin glasses when reducing the temperature at fixed field. For σ=0.75, there is evidence from finite-size-scaling studies for an AT transition but for σ=0.85, the evidence for a transition is nonexistent. We have also studied these systems at fixed temperature varying the field and discovered that at both σ=0.75 and at σ=0.85 there is evidence of an AT transition! Confusingly, the correlation length and spin-glass susceptibility as a function of the field are both entirely consistent with the predictions of the droplet picture and hence the nonexistence of an AT line. In the usual finite-size critical point scaling studies used to provide evidence for an AT transition, there is seemingly good evidence for an AT line at σ=0.75 for small values of the system size N, which is strengthening as N is increased, but for N>2048 the trend changes and the evidence then weakens as N is further increased. We have also studied with fewer bond realizations the system at σ=0.70, which is the analog of a system with short-range interactions just below six dimensions, and found that it is similar in its behavior to the system at σ=0.75 but with larger finite-size corrections. The evidence from our simulations points to the complete absence of the AT line in dimensions outside the mean-field region and to the correctness of the droplet picture. Previous simulations which suggested there was an AT line can be attributed to the consequences of studying systems which are just too small. The collapse of our data to the droplet scaling form is poor for σ=0.75 and to some extent also for σ=0.85, when the correlation length becomes of the order of the length of the system, due to the existence of excitations which only cost a free energy of O(1), just as envisaged in the TNT picture of the ordered state of spin glasses. However, for the case of σ=0.85 we can provide evidence that for larger system sizes, droplet scaling will prevail even when the correlation length is comparable to the system size.

Free-energy barriers in the Sherrington-Kirkpatrick model.

The free-energy landscape of the Sherrington-Kirkpatrick (SK) Ising spin glass is simple in the framework of the Thouless-Anderson-Palmer (TAP) equations as each solution (which are minima of the free energy) has associated with it a nearby index-one saddle point. The free-energy barrier to escape the minimum is just the difference between the saddle point free energy and that at its associated minimum. This difference is calculated for the states with free energies f>f_{c}. It is very small for these states, decreasing as 1/N^{2}, where N is the number of spins in the system. These states are not marginally stable. We argue that such small barriers are why numerical studies never find these states when N is large. Instead, the states that are found are those that have marginal stability. For them the barriers are at least of O(1). f_{c} is the free energy per spin below which the states develop broken replica-symmetry-like overlaps with each other. In the regime f

Droplet-scaling versus replica symmetry breaking debate in spin glasses revisited.

Simulational studies of spin glasses since the early 2010s have focused on the so-called replicon exponent α as a means of determining whether the low-temperature phase of spin glasses is described by the replica symmetry breaking picture of Parisi or by the droplet-scaling picture. On the latter picture, it should be zero, but we shall argue that it will only be zero for systems of linear dimension L>L^{*}. The crossover length L^{*} may be of the order of hundreds of lattice spacings in three dimensions and approach infinity in six dimensions. We use the droplet-scaling picture to show that the apparent nonzero value of α when L

Marginally jammed states of hard disks in a one-dimensional channel.

We have studied a class of marginally jammed states in a system of hard disks confined in a narrow channel-a quasi-one-dimensional system-whose exponents are not those predicted by theories valid in the infinite dimensional limit. The exponent γ which describes the distribution of small gaps takes the value 1 rather than the infinite dimensional value 0.41269⋯. Our work shows that there exist jammed states not found within the tiling approach of Ashwin and Bowles. The most dense of these marginal states is an unusual state of matter that is asymptotically crystalline.

Possible instability of one-step replica symmetry breaking in p-spin Ising models outside mean-field theory.

The fully connected Ising p-spin model has for p>2 a discontinuous phase transition from the paramagnetic phase to a stable state with one-step replica symmetry breaking (1RSB). However, simulations in three dimension do not look like these mean-field results and have features more like those which would arise with full replica symmetry breaking (FRSB). To help understand how this might come about we have studied in the fully connected p-spin model the state of two-step replica symmetry breaking (2RSB). It has a free energy degenerate with that of 1RSB, but the weight of the additional peak in P(q) vanishes. We expect that the state with full replica symmetry breaking (FRSB) is also degenerate with that of 1RSB. We suggest that finite-size effects will give a nonvanishing weight to the FRSB features, as also will fluctuations about the mean-field solution. Our conclusion is that outside the fully connected model in the thermodynamic limit, FRSB is to be expected rather than 1RSB.

Realizable solutions of the Thouless-Anderson-Palmer equations.

We show that the only solutions of the Thouless-Anderson-Palmer (TAP) equations for the Sherrington-Kirkpatrick model of Ising spin glasses which can be found by iteration are those whose free energy lies on the border between replica-symmetric and broken-replica-symmetric states, when the number of spins N is large. Convergence to this same borderline also happens in quenches from a high-temperature initial state to a locally stable state where each spin is parallel to its local field; both are examples of self-organized criticality. At this borderline the band of eigenvalues of the Hessian associated with a solution extends to zero, so the states reached have marginal stability. We have also investigated the factors which determine the free-energy difference between a stationary solution corresponding to a saddle point and its associated minimum, which is the barrier which has to be surmounted to escape from the vicinity of a TAP minimum or pure state.

Absence of Hyperuniformity in Amorphous Hard-Sphere Packings of Nonvanishing Complexity.

We relate the structure factor S(k→0) in a system of jammed hard spheres of number density ρ to its complexity per particle Σ(ρ) by the formula S(k→0)=-1/[ρ^{2}Σ^{″}(ρ)+2ρΣ^{'}(ρ)]. We have verified this formula for the case of jammed disks in a narrow channel, for which it is possible to find Σ(ρ) and S(k) analytically. Hyperuniformity, which is the vanishing of S(k→0), will therefore not occur if the complexity is nonzero. An example is given of a jammed state of hard disks in a narrow channel which is hyperuniform when generated by dynamical rules that produce a nonextensive complexity.

Gardner Transition in Physical Dimensions.

The Gardner transition is the transition that at mean-field level separates a stable glass phase from a marginally stable phase. This transition has similarities with the de Almeida-Thouless transition of spin glasses. We have studied a well-understood problem, that of disks moving in a narrow channel, which shows many features usually associated with the Gardner transition. We show that some of these features are artifacts that arise when a disk escapes its local cage during the quench to higher densities. There is evidence that the Gardner transition becomes an avoided transition, in that the correlation length becomes quite large, of order 15 particle diameters, even in our quasi-one-dimensional system.

Fractal dimension of interfaces in Edwards-Anderson spin glasses for up to six space dimensions.

The fractal dimension of domain walls produced by changing the boundary conditions from periodic to antiperiodic in one spatial direction is studied using both the strong-disorder renormalization group algorithm and the greedy algorithm for the Edwards-Anderson Ising spin-glass model for up to six space dimensions. We find that for five or fewer space dimensions, the fractal dimension is lower than the space dimension. This means that interfaces are not space filling, thus implying that replica symmetry breaking is absent in space dimensions fewer than six. However, the fractal dimension approaches the space dimension in six dimensions, indicating that replica symmetry breaking occurs above six dimensions. In two space dimensions, the strong-disorder renormalization group results for the fractal dimension are in good agreement with essentially exact numerical results, but the small difference is significant. We discuss the origin of this close agreement. For the greedy algorithm there is analytical expectation that the fractal dimension is equal to the space dimension in six dimensions and our numerical results are consistent with this expectation.

Multicritical Point on the de Almeida-Thouless Line in Spin Glasses in d>6 Dimensions.

The de Almeida-Thouless (AT) line in Ising spin glasses is the phase boundary in the temperature T and magnetic field h plane below which replica symmetry is broken. Using perturbative renormalization group (RG) methods, we show that, when the dimension d of space is just above six, there is a multicritical point (MCP) on the AT line, which separates a low-field regime, in which the critical exponents have mean-field values, from a high-field regime, where the RG flows run away to infinite coupling strength; as d approaches six from above, the MCP approaches the zero-field critical point exponentially in 1/(d-6). Thus, on the AT line, perturbation theory for the critical properties breaks down at a sufficiently large magnetic field even above 6 dimensions, as well as for all nonzero fields when d≤6, as was known previously. We calculate the exponents at the MCP to first order in ϵ=d-6>0. The fate of the MCP as d increases from just above six to infinity is not known.

Fractal Dimension of Interfaces in Edwards-Anderson and Long-range Ising Spin Glasses: Determining the Applicability of Different Theoretical Descriptions.

The fractal dimension of excitations in glassy systems gives information on the critical dimension at which the droplet picture of spin glasses changes to a description based on replica symmetry breaking where the interfaces are space filling. Here, the fractal dimension of domain-wall interfaces is studied using the strong-disorder renormalization group method pioneered by Monthus [Fractals 23, 1550042 (2015)FRACEG0218-348X10.1142/S0218348X15500425] both for the Edwards-Anderson spin-glass model in up to 8 space dimensions, as well as for the one-dimensional long-ranged Ising spin-glass with power-law interactions. Analyzing the fractal dimension of domain walls, we find that replica symmetry is broken in high-enough space dimensions. Because our results for high-dimensional hypercubic lattices are limited by their small size, we have also studied the behavior of the one-dimensional long-range Ising spin-glass with power-law interactions. For the regime where the power of the decay of the spin-spin interactions with their separation distance corresponds to 6 and higher effective space dimensions, we find again the broken replica symmetry result of space filling excitations. This is not the case for smaller effective space dimensions. These results show that the dimensionality of the spin glass determines which theoretical description is appropriate. Our results will also be of relevance to the Gardner transition of structural glasses.

Metastable minima of the Heisenberg spin glass in a random magnetic field.

We have studied zero-temperature metastable minima in classical m-vector component spin glasses in the presence of m-component random fields for two models, the Sherrington-Kirkpatrick (SK) model and the Viana-Bray (VB) model. For the SK model we have calculated analytically its complexity (the log of the number of minima) for both the annealed case where one averages the number of minima before taking the log and the quenched case where one averages the complexity itself, both for fields above and below the de Almeida-Thouless (AT) field, which is finite for m>2. We have done numerical quenches starting from a random initial state (infinite temperature state) by putting spins parallel to their local fields until there is no further decrease of the energy and found that in zero field it always produces minima that have zero overlap with each other. For the m=2 and m=3 cases in the SK model the final energy reached in the quench is very close to the energy E_{c} at which the overlap of the states would acquire replica symmetry-breaking features. These minima have marginal stability and will have long-range correlations between them. In the SK limit we have analytically studied the density of states ρ(λ) of the Hessian matrix in the annealed approximation. Despite the fact that in the presence of a random field there are no continuous symmetries, the spectrum extends down to zero with the usual sqrt[λ] form for the density of states for fields below the AT field. However, when the random field is larger than the AT field, there is a gap in the spectrum, which closes up as the AT field is approached. The VB model behaves differently and seems rather similar to studies of the three-dimensional Heisenberg spin glass in a random vector field.

Interface free-energy exponent in the one-dimensional Ising spin glass with long-range interactions in both the droplet and broken replica symmetry regions.

The one-dimensional Ising spin-glass model with power-law long-range interactions is a useful proxy model for studying spin glasses in higher space dimensions and for finding the dimension at which the spin-glass state changes from having broken replica symmetry to that of droplet behavior. To this end we have calculated the exponent that describes the difference in free energy between periodic and antiperiodic boundary conditions. Numerical work is done to support some of the assumptions made in the calculations and to determine the behavior of the interface free-energy exponent of the power law of the interactions. Our numerical results for the interface free-energy exponent are badly affected by finite-size problems.

Glasslike behavior of a hard-disk fluid confined to a narrow channel.

Disks moving in a narrow channel have many features in common with the glassy behavior of hard spheres in three dimensions. In this paper we study the caging behavior of the disks that sets in at characteristic packing fraction ϕ(d). Four-point overlap functions similar to those studied when investigating dynamical heterogeneities have been determined from event-driven molecular dynamics simulations and the time-dependent dynamical length scale has been extracted from them. The dynamical length scale increases with time and, on the equilibration time scale, it is proportional to the static length scale associated with the zigzag ordering in the system, which grows rapidly above ϕ(d). The structural features responsible for the onset of caging and the glassy behavior are easy to identify as they show up in the structure factor, which we have determined exactly from the transfer-matrix approach.

Finite-size critical scaling in Ising spin glasses in the mean-field regime.

We study in Ising spin glasses the finite-size effects near the spin-glass transition in zero field and at the de Almeida-Thouless transition in a field by Monte Carlo methods and by analytical approximations. In zero field, the finite-size scaling function associated with the spin-glass susceptibility of the Sherrington-Kirkpatrick mean-field spin-glass model is of the same form as that of one-dimensional spin-glass models with power-law long-range interactions in the regime where they can be a proxy for the Edwards-Anderson short-range spin-glass model above the upper critical dimension. We also calculate a simple analytical approximation for the spin-glass susceptibility crossover function. The behavior of the spin-glass susceptibility near the de Almeida-Thouless transition line has also been studied, but here we have only been able to obtain analytically its behavior in the asymptotic limit above and below the transition. We have also simulated the one-dimensional system in a field in the non-mean-field regime to illustrate that when the Imry-Ma droplet length scale exceeds the system size one can then be erroneously lead to conclude that there is a de Almeida-Thouless transition even though it is absent.

Understanding the ideal glass transition: lessons from an equilibrium study of hard disks in a channel.

We use an exact transfer-matrix approach to compute the equilibrium properties of a system of hard disks of diameter σ confined to a two-dimensional channel of width 1.95σ at constant longitudinal applied force. At this channel width, which is sufficient for next-nearest-neighbor disks to interact, the system is known to have a great many jammed states. Our calculations show that the longitudinal force (pressure) extrapolates to infinity at a well-defined packing fraction ϕ(K) that is less than the maximum possible ϕ(max), the latter corresponding to a buckled crystal. In this quasi-one-dimensional problem there is no question of there being any real divergence of the pressure at ϕ(K). We give arguments that this avoided phase transition is a structural feature, the remnant in our narrow channel system of the hexatic to crystal transition, but that it has the phenomenology of the (avoided) ideal glass transition. We identify a length scale ξ̃(3) as our equivalent of the penetration length for amorphous order: In the channel system, it reaches a maximum value of around 15σ at ϕ(K), which is larger than the penetration lengths that have been reported for three-dimensional systems. It is argued that the α-relaxation time would appear on extrapolation to diverge in a Vogel-Fulcher manner as the packing fraction approaches ϕ(K).

Sample convection in liquid-state NMR: why it is always with us, and what we can do about it.

Many NMR experiments on liquids suffer if the sample convects. This is particularly true for applications, such as the measurement of diffusion, that rely on spatial labelling of spins. It is widely assumed that, in most well-conducted experiments with stable temperature regulation, samples do not convect. Unfortunately this is not the case. It is shown here that typical NMR samples show measurable convective flow for all but a very narrow range of temperatures; convection is seen both above and below this range, which can be as small as a degree or so for a mobile solvent such as chloroform. This convection is driven by both vertical and horizontal temperature gradients. Measurements of convection velocity are presented for a range of samples, sample tubes, probes, and temperatures. Both decreasing sample tube inner diameter and changing sample tube material from glass to sapphire can slow convection markedly, with sapphire tubes being particularly effective. Such tubes are likely to be particularly helpful for accurate measurement of diffusion by NMR.

Dealing with correlated choices: how a spin-glass model can help political parties select their policies.

Starting from preferences on N proposed policies obtained via questionnaires from a sample of the electorate, an Ising spin-glass model in a field can be constructed from which a political party could find the subset of the proposed policies which would maximize its appeal, form a coherent choice in the eyes of the electorate, and have maximum overlap with the party's existing policies. We illustrate the application of the procedure by simulations of a spin glass in a random field on scale-free networks.

Boolean decision problems with competing interactions on scale-free networks: equilibrium and nonequilibrium behavior in an external bias.

We study the equilibrium and nonequilibrium properties of Boolean decision problems with competing interactions on scale-free networks in an external bias (magnetic field). Previous studies at zero field have shown a remarkable equilibrium stability of Boolean variables (Ising spins) with competing interactions (spin glasses) on scale-free networks. When the exponent that describes the power-law decay of the connectivity of the network is strictly larger than 3, the system undergoes a spin-glass transition. However, when the exponent is equal to or less than 3, the glass phase is stable for all temperatures. First, we perform finite-temperature Monte Carlo simulations in a field to test the robustness of the spin-glass phase and show that the system has a spin-glass phase in a field, i.e., exhibits a de Almeida-Thouless line. Furthermore, we study avalanche distributions when the system is driven by a field at zero temperature to test if the system displays self-organized criticality. Numerical results suggest that avalanches (damage) can spread across the whole system with nonzero probability when the decay exponent of the interaction degree is less than or equal to 2, i.e., that Boolean decision problems on scale-free networks with competing interactions can be fragile when not in thermal equilibrium.