Anatomy of a Spin: The Information-Theoretic Structure of Classical Spin Systems

Collective organization in matter plays a significant role in its expressed physical properties. Typically, it is detected via an order parameter, appropriately defined for each given system's observed emergent patterns. Recent developments in information theory, however, suggest quantifying collective organization in a system- and phenomenon-agnostic way: decompose the system's thermodynamic entropy density into a localized entropy, that solely contained in the dynamics at a single location, and a bound entropy, that stored in space as domains, clusters, excitations, or other emergent structures. We compute this decomposition and related quantities explicitly for the nearest-neighbor Ising model on the 1D chain, the Bethe lattice with coordination number k=3, and the 2D square lattice, illustrating its generality and the functional insights it gives near and away from phase transitions. In particular, we consider the roles that different spin motifs play (in cluster bulk, cluster edges, and the like) and how these affect the dependencies between spins.


I. INTRODUCTION
Collective behavior underlies a vast array of fascinating phenomena, many of critical importance to contemporary science and technology. Unusual material propertiessuch as superconductivity, metal-insulator transitions, and heavy fermions-have been attributed to collective behavior arising from the compounded interaction of system components [1]. Collective behavior is by no means limited to complex materials, however: the behavior of financial markets [2], epileptic seizures [3] and consciousness [4], and animal flocking behavior [5] are all now also seen as examples. And, we now appreciate that it appears most prominently via phase transitions [6, and references therein]. Operationally, collective behavior is detected and quantified using correlations or, more recently, by mutual information [7], estimated either locally via pairwise component interactions or globally. Exploring spin systems familiar from statistical mechanics, here we argue that these diagnostics can be substantially refined to become more incisive tools for quantifying collective behavior. This is part of the larger endeavor of discovering ways to automatically detect collective behavior and the emergence of organization [8]. Along these lines, much effort has been invested to explore information-theoretic properties of the Ising model of statistical mechanics. Both Shaw [9] and Arnold [10] studied information in 1D spin-strips within the 2D Ising model. Feldman and Crutchfield [11][12][13] explored several generalizations of the excess entropy, a well known mutual information measure for time series, in 1D and 2D Ising models as a function of coupling strength, showing that they were sensitive to spin patterns and can be used as a generalized order parameter. In one of the more thorough-going studies to date, Lau and Grassberger [14] computed the excess entropy using adjacent rings of spins to probe criticality in the 2D cylindrical Ising model. Barnett et al. [6] tracked information flows in a kinetic Ising system using both the pairwise and global mutual informations (generalized as the total correlation [15], it appears below) and the transfer entropy [16]. Abdallah and Plumbley [17] employed an extensive version of the dual total correlation [15] on small, random spin glasses. Despite their successful application, the measures used to monitor collective phenomena in these studies were motivated information theoretically rather than physically, and so their thermodynamic relevance and structural interpretation have remained unclear.
To address this we decompose the thermodynamic entropy density for spin systems into a set of information measures, some already familiar in information theory [18], complex systems [19], and elsewhere [17]. For example, one measure that falls out naturally-the bound entropyis that part of the thermodynamic entropy accounting for entropy shared between spins. This is in contrast to monitoring collective behavior as the difference between the thermodynamic entropy and the entropy of a hypothetical distribution over uncorrelated spins-the total correla-tion mentioned above. In this way, the decomposition provides a physical basis for informational measures of collective behavior, in particular showing how the measures lead to interpretations of emergent structure in system configurations.
To make this argument and illustrate its consequences, Section II first defines notation and lays the groundwork for our decomposition. Section III then derives the decomposition. Though this decomposition is given for regular spin lattices, in principle it is meaningful for any system with a well defined thermodynamic entropy. Section IV outlines the computational methods for estimating the resulting quantities, using them to interpret emergent organization in the nearest-neighbor ferromagnetic Ising model in one dimension (Section IV A), in the Bethe lattice with coordination number k = 3 (Section IV B), and in the two-dimensional square Ising lattice (Section IV C). Having established the different types of entropy in spin systems and connections to their underlying physical structures, we conclude by suggesting applications and future directions.

II. SPIN ENTROPIES
We write σ to denote a configuration of spins on a lattice L, σ i for the particular spin state at lattice site i ∈ L, and σ \i the collection of all spin states in a configuration excluding σ i at site i. The random variable over all possible spin configurations is denoted σ, for a particular spin variable σ i , and for all spin variables but the one at site i, σ \i . As a shorthand, we write σ ∈ σ and similar to mean indexing into the event space of the random variable σ. We study the ferromagnetic spin-1/2 Ising model with nearest-neighbor interactions in the thermodynamic limit whose Hamiltonian is given by: where i, j denotes all pairs (i, j) such that the sites i and j are directly connected in L and the interaction strength J is positive. We assume the system is in equilibrium and isolated. And so, the probability of configuration σ occurring is given by the Boltzmann distribution: where Z is the partition function.
The Boltzmann entropy of a statistical mechanical system assumes the constituent degrees of freedom (here spins) are uncorrelated. To determine it, we use the isolated spin entropy: where p (↑) = (1 + m)/2 is the probability of a spin being up, p (↓) = (1 − m)/2 is the probability of a spin being down, and m = (# ↑ −# ↓)/ |L| is average magnetization in a configuration. The site index 0 was chosen arbitrarily and represents any single spin in the lattice. The system's Boltzmann entropy is then the extensive quantity H B = |L| · H [σ 0 ]. As Jaynes [20] emphasizes, though correct for an ideal gas, H B is not the empirically correct system entropy if a system develops internal correlations. In our case, if spins are correlated across the lattice, one must consider entire configurations, leading to the Gibbs entropy: and the thermodynamic entropy density h, the entropy per spin when the entire lattice is considered:

III. DECOMPOSING A SPIN'S THERMODYNAMIC ENTROPY
Since we are particularly interested here in monitoring the appearance of internal correlations, the difference between the Boltzmann and Gibbs entropies suggests itself as our first measure of a system's internal organization.
If each spin were independent, H [σ 0 ] = h. For example, this is true at T = ∞ for the Ising model and for site percolation models [24]. If there are correlations between spins, however, then h < H [σ 0 ], as just noted for the extensive quantities. Their difference is the total correlation density [15]: |L| , here is the Kullback-Leibler divergence [7] between the distribution over entire configurations and the product of isolated-spin marginal distributions. Now called the total correlation [15], this measure of internal correlation was wholly anticipated by Gibbs, as Jaynes notes [20,Eq. (4)]. Since ρ vanishes only when each spin is independent, one might consider it a measure of pattern or structure in a spin system. While this is certainly reasonable as an operational definition, our next step decomposes both h and ρ into more nuanced quantities.
Continuing, we recast the Gibbs thermodynamic entropy H [σ] as the sum of two nonnegative terms; starting from the configuration entropy in its standard statistical form Eq. (3) and manipulating it into two terms: where R [σ] and B [σ], known as the residual entropy and dual total correlation [15,17,18], are the following measures: and Similarly, we can begin with the total correlation of spins T [σ] and break it in to two terms: and again, we consider the spatial density q = Q[σ] /|L| in the thermodynamic limit.

It can be difficult to initially intuit the difference between T [σ] and B [σ]
. Both measure dependencies among spins, but against different reference ensembles. T [σ] is the difference between the Boltzmann and Gibbs entropies (Eq. (5)) and, therefore, quantifies the spin distribution's deviation from a hypothetical distribution in which each spin is independent from the others. Since R [σ] is the amount of actual independent entropy in the spin distribution, B [σ] quantifies the amount of spin-spin dependency in configurations not in reference to a hypothetical independent-spin distribution, but rather to the actually realized amount of independence. This difference is pictorially represented in Fig. 1.
Let us now interpret the physics captured by these entropy components. First, H [σ 0 ] quantifies the entropy per spin ignoring any dependencies (e.g. correlations) between spins. The thermodynamic entropy density h, however, is the entropy per spin including dependencies. Continuing in this way, from Eq. (7) r is the entropy per spin remaining after these dependencies have been factored out, meaning it is the average amount of independent entropy-the entropy per spin when we know the state of the other spins in the lattice. With this logic, b is that portion of the thermodynamic entropy density coming from dependencies and is, therefore, the average amount of dependent entropy per spin. As we shall see, in comparison to ρ, b provides a complementary, quantitatively different, and more physically grounded approach to quantifying dependencies among spins.
This completes, in effect, our decomposition of the isolated-spin entropy H [σ 0 ], shown schematically in Fig. 2. To take stock, let's back up a bit. Recall that ρ is the difference between H [σ 0 ], which ignores the dependencies between the spins, and the thermodynamic entropy density h. That is, ρ is the difference between the isolated- . First, into the Gibbs thermodynamic entropy density h and the total correlation ρ. Second, the thermodynamic entropy into the localized density r and the dual total correlation information density b, and the total correlation density ρ into b and the enigmatic information density q.
spin entropy and how much entropy (h) there actually is. Whereas, b is the difference between the thermodynamic entropy h and the amount of independent entropy r there is. In addition, b is also part of ρ [15], leaving q which therefore quantifies potential dependencies (explained shortly) not present in the thermodynamic entropy h.
Let's explore the consequences of this decomposition. Although r and b are components of the thermodynamic entropy density h-itself a measure of the irreducible perspin randomness-individually they have very different meanings that give insight into the interplay of spin thermalization and spin ordering. The first, r, is the average amount of unshared entropy, quantifying randomness remaining in a single spin given the exact configuration of the rest of the lattice. For this reason we refer to it as the localized entropy. In the limit of high and low temperature, r's behavior is intuitive. At low temperature, spins in the lattice are entirely aligned, and so there is no randomness in individual sites given the global alignment. At high temperatures, each spin is independent, and so there is full randomness in a site, even knowing its surroundings.
The other quantity, b, is the entropy remaining in a spin when the localized entropy is deducted, and so it captures the thermodynamic entropy shared between a spin and the rest of the lattice. We refer to it as the bound entropy. Its limits are similarly intuitive. At low temperatures there is no randomness to share, and so b = 0. Similarly, at high temperatures all randomness is localized, and so no entropy can be shared and bound entropy also vanishes. At intermediate temperatures, though spins are in flux, a spin's context does at least in part influence its alignment and so b > 0. We refer to the temperature at which b is maximized by T b .
Finally, we consider the enigmatic information density q. As a component of ρ it also reflects dependencies, but in a way different from b. The bound entropy captures dependencies that are also part of the thermodynamic entropy, but q captures those that are not. As we shall see in Section V, this means it is sensitive to the interior of clusters. At low temperatures, i.e., temperatures below T c , we expect q to be small, since H [σ 0 ] is small. At high temperatures, there are no dependencies between spins, and so we also expect q to be small. At intermediate temperatures, just like b, we expect q to take larger values.
The nonzero values of ρ, b, and q are all driven by Ising spin clustering. Consider for the moment randomly permuting the sites in the lattice, resulting in a site percolation model with site occupancy probability given by a simple function of the magnetization of the original Ising lattice. Then, the spin orientation at any particular site is now statistically independent of its neighbors, and so ρ = b = q = 0, and therefore H [σ 0 ] = h = r. There are, however, clusters in this permuted lattice-in fact, there may exist a giant cluster. These clusters are driven purely by the lattice topology. We therefore conclude that ρ, b, and q are each sensitive to differing aspects of how the pairwise Ising Hamiltonian affects the shape and distribution of spin clusters.
To summarize, the disorder generated at each site is measured by the Gibbs thermodynamic entropy density h. And, one portion (b) participates in spatial organization and the other portion (r) does not. Thus, b is key to monitoring the emergence of spatial structures during pattern formation dynamics.

IV. RESULTS
To explore what the entropy components reveal, we calculate them for the ferromagnetic spin-1 /2 Ising model on one-(1D) and two-dimensional (2D) square lattices and on the Bethe lattice. These results demonstrate that our intuitions regarding the behaviors of the information measures are correct, but also raise subtle questions regarding the details of exactly what they capture.
As part of this, we report analytical results for the 1D chain and the Bethe lattice and show that they match results from simulation. And, we use analytic results whenever possible for the 2D lattice. For example, we report analytical results for h, H [σ 0 ], and ρ. Then, only r needs to be estimated from simulation and that result is then combined with analytic quantities to compute b and q. To provide a comparison, estimating b purely from simulation requires accurate h and r. In this, we make use of the global Markovian property of spin lattices [26]. To obtain r we need only condition a spin on all its nearest neighbors-spins directly coupled via the Hamiltonian of Eq. (1). For h estimated purely from simulation, we use Ref. [27]'s method. In both cases, Ref. [28]'s entropy estimator H 2 was used to decrease the number of configuration samples while still providing convergence to known analytic results. Overall, when comparing and when analytic results are available, simulation results match to within standard error bars smaller than plot line widths. All simulation results were obtained using the Wolff algorithm [29]. In summary, in the reported results for those several quantities that required estimated or partial estimation from simulations, sufficient care was taken that the statistical errors fall far below the level needed to support the main conclusions.

A. 1D Ising Model
Simple Ising models have been broadly adapted to understand simple collective phenomena in fields ranging from surface physics [30,31] and biophysics [32,33] to economics [34]. The Ising model on a one-dimensional lattice (L = Z) can be fully analyzed in closed form. Unfortunately, its emergent patterns are rather constrained; for example, it does not exhibit a phase transition [35]. However, the exact solutions provide an important benchmark and so it is a necessary candidate with which to start and provides a base from which to generalize [36][37][38][39]. Analytic results were computed using a transfer matrix approach combined with the aforementioned Markovian property. For comparison, simulations were performed on a lattice of size N = 1024.
Independent of sporting a phase transition or not, the argument that b maximizes between temperature extremes still holds and Fig. 3 verifies this. The bound entropy b reaches a maximal value at T b ≈ 1.0117 J/k B and this is not (and cannot be) associated with a phase transition. At lower temperatures, spin-spin shared entropy b is the dominant contributor to the Gibbs thermodynamic entropy density h = r + b. Whereas at higher temperatures, the localized entropy r is the dominant contributor. Similarly, at low temperatures q is the dominant contributor to ρ, and b is its dominant contributor at high temperatures. Section V below isolates why b peaks where it does and to which spin-configuration features each measure is sensitive. For now, let's continue with system informational phenomenology.

B. Bethe Lattice Ising Model
What role does the underlying lattice topology play in determining the balance between local randomness and spatially shared information that equilibrium configurations achieve? Spins on the Bethe lattice, in which L is an infinite Cayley tree with coordination number k, is an ideal candidate since the absence of loops makes it possible to compute quantities analytically. Moreover, spin configurations exhibit a phase transition at critical temperature [36]: We analytically calculated the entropy decomposition for a Bethe lattice with coordination number k, the details of which appear in Appendix A. Simulation results, matching the analytic to high precision, were performed on a 1, 000 node random 3-regular graph [40], which has no boundary and is locally Bethe lattice-like. Figure 4 presents the results for a Bethe lattice with coordination number k = 3, though other ks behave similarly. Interestingly, all information measures have a discontinuity in their first derivatives, and this happens at the phase transition at T c . Furthermore, the bound entropy b is maximized there: Opposite the 1D lattice, at small T the dominant contributor to Gibbs entropy h is the local randomness r. Thus, not only does the change in lattice topology induce a phase transition, but it also inverts the informational components' contributions to thermodynamic entropy density. This, in turn, indicates a rather different underlying mechanism that supports entropy production in the low temperature regime. In the 1D case entropy is spatially extended and configurations of spins are dominated by large clusters.  Whereas, in the Bethe lattice deviations from uniformity come in the form of isolated spins differing from their surroundings. Also, unlike 1D, the spin-spin shared entropy b is the dominant contributor to ρ at both low and high temperatures.

C. 2D Ising Model
Unlike its 1D sibling, the nearest-neighbor ferromagnetic Ising model in two dimensions (L = Z 2 ) has a phase transition at a finite critical temperature T c = 2 /ln(1+ √ 2) ≈ 2.2692 J/k B . In the 2D case, although Onsager's solution lets us calculate H [σ 0 ], h, and ρ [36], we do not have an analytic form for r and so the curves for r, b, and q in Fig. 5 partly rely on estimates from simulation, as explained above. The simulations were conducted on a 128 × 128 lattice and quantities averaged over 200, 000 configuration updates. Again, the resulting standard error bars were smaller than the plotted line widths.
Phase transition aside, the behaviors of h, r, and b, seen in Fig. 5, are qualitatively similar to those in the 1D lattice. However, unlike 1D, at low temperatures r is the dominant contributor to h, similar to the Bethe lattice below its transition. Also, paralleling the Bethe lattice, b is the dominant contributor to ρ at both low and high temperatures. Like the 1D system, bound entropy is maximized within the disordered phase at T b ≈ 2.675 J/k B , above T c . Let us now turn to what configurational structures lead to the overall behaviors of these measures and how they contribute to the Gibbs entropy density h.

V. DISCUSSION
At first, it is striking that T b is not generically identical to T c . To understand why, we need to investigate to which spin-configuration motifs the various measures are sensitive. To accomplish this, we take a local approach. Each of the measures examined above is an average density, which is important when discussing the lattice as a whole and its bulk thermodynamic properties. However, each spin configuration motif contributes differently to this average. For example, an up spin surrounded by down spins contributes to the measures differently than an up spin surrounded by other up spins.

A. Motif Entropies
To quantify this, we appeal to spatial-and temporallocal forms of the averaged densities considered so far. This is possible on square lattices L = Z d , due to the existence of conditional forms of the entropy density; see Ref. [41,Thm. 2.9]. The conditional forms for each of the component quantities are configuration-weighted averages over quantities of the form log 2 (•) [15]. For example, the thermodynamic entropy density of Eq. (4) can be shown to be [41]: where h i (σ) = − log 2 p (σ i | ← − σ i ) and ← − σ i is the set of spins whose indices are "lexicographically" less than i. h i (σ) is a spatially local or pointwise measure and quantifies how much the individual spin σ i contributes to the global entropy. We can more directly understand how local motifs contribute to the average entropy by plotting the local measures at each site within a given lattice configuration. This approach roughly parallels that in, e.g., Ref. [42], but here we employ different information measures and do not assume directionality, since we average the pointwise values resulting from each of the possible orientations centered at the given spin.
In the following, we define a cluster to be a maximal set of contiguous spins oriented identically. The edge of a cluster is that set of spins in a cluster with at least one neighbor not in the cluster. The remaining spins in the cluster are its bulk. An isolated spin is a cluster consisting of a single spin. There are several immediate similarities that allow for structural interpretation. For example, motif q is positive within the bulk of a cluster and negative on its edges. Highlighting opposite features, motif r is negligible within the bulk, but positive along edges, particularly corners and isolated spins. The motif bound entropy b, however, is more nuanced in its behavior. Considering Figure 6, at r r r r r r r r r r r r r r r r r r r r r q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q lower temperatures it is sensitive to cluster edges more so than cluster bulk. At higher temperatures, however, this relationship flips and it is more sensitive to the bulk. For all temperatures, it is negative for isolated spins-clusters of size 1.
In contrast, we see that at not-too-high temperatures motif b is largely sensitive to cluster boundaries. This leads one to speculate that, as in Ref. [6], b's maximization in the disordered phase is due to complex interactions between cluster sizes and their surface area. This may also shed light on why b does in fact peak at T c on the Bethe lattice. On the square lattice, on the one hand, cluster boundaries have a tendency to smooth due to the presence of correlations flowing along short loops in the lattice topology that tend to align spins on a boundary. On the other, the Bethe lattice has no such loops, and so there is no pressure to reduce surface area. This suggests that b will generically peak at T c for systems where boundaries are not constrained by system energetics and that the energetic smoothing of cluster boundaries drives b to peak in the more "textured" disordered phase.
Interestingly, for general lattices such as the Bethe lattice, but especially for random graph topologies, local (conditional) forms of the thermodynamic entropy density (analogs of Eq. (13)) and other motif measures are unknown and may simply not exist. While techniques, such as Ref. [43]'s, exist for estimating thermodynamic entropy density for a ferromagnetic Ising model on an arbitrary graph, the interpretation of a thermodynamic entropy density in such systems is problematic as each spin site may have differing connectivity. Therefore, while global averages may exist for arbitrary topologies and may even be tractably estimated, their structural meaning is vastly more challenging.
And, this is critical to understanding thermodynamic and statistical mechanical properties of such systems and, more generally, correlations and information measures on network dynamical systems. If local information estimates, such as those used for the motif entropy analyses above, existed, then studies can be undertaken that determine, say, which nodes in a lattice or network graph contribute most to collective behavior. Failing such measures, though, one must be wary of claims to answer such questions. These concerns are touched on by the informational analyses of frustration in spin glasses by Robinson et al. [44], whose local entropies are temporally averaged, spatially local entropies in a heterogeneous system. They are unlike our motif entropies, though, which are both spatially and temporally local and, therefore, have meaning only in a homogeneous system. More study is required, however.

B. Actively Storing Information versus Communicating It
Assuming Glauber dynamics for the nearest-neighbor 2D Ising system, Barnett et al. [6] computed a global transfer entropy [16]-the average information communicated between the entire lattice and a single spin. They found that it peaks within the disordered phase at a temperature Although there are serious concerns with how the transfer entropy conflates two-spin (dyadic) and multispin (polyadic) dependencies [45], those concerns are not relevant in these Ising systems due to their known, dyadic Hamiltonian. To compare our measures with theirs, we examine the active storage of randomly generated information in spatial patterns, as measured by the ratio of b /h-how much of the localized randomness (h) is converted by a system into spin-spin (spatially) shared information (b). Figure 8 plots the ratio as a function of temperature for the 2D Ising system: the ratio peaks in the disordered phase and, surprisingly, it appears to do so at T T gl . This implies that there is a strong connection between b (storing thermal fluctuations as spatial correlation) and the potential for information communication in a system. Hopefully, the result suggests that Barnett et al.'s T T gl is not a result specific only to their choice of Glauber dynamics, but rather an intrinsic property of the 2D Ising model. It remains to be seen why spatially shared information is maximized within the disordered phase of the 2D Ising system, though.

VI. CONCLUSIONS
As noted at the beginning, even the earliest debates over entropy's statistical foundation turned on contrasting its thermodynamic components-Boltzmann's isolated-spin entropy versus the Gibbs global entropy. From a modern perspective their difference, well known to Gibbs, is a generalized mutual information that measures the degree of individual-component independence, now called the total correlation. Keying off this, we showed that the Gibbs thermodynamic entropy density naturally decomposes further into two functionally distinct informational components: one quantifying independence among constituent spins and the other, dependence. The one quantifying dependence, the bound entropy b, captures collective behavior by expressing how much of the thermodynamic entropy density is locally shared. We then demonstrated the behavior of the bound entropy and related quantities for the nearest-neighbor ferromagnetic Ising model on a variety of lattices. We found that it tends to a maximum at intermediate temperatures, though not always at the magnetic phase transition. Our analyses support our earlier hypothesis [46] that q, as the dominant component of the persistent mutual information [47], should generically be maximized at critical points. Though not detailed here, this observation holds in three dimensions and in simulations of the Potts model [48] with the number of states s = 3 on lattices with 1, 2, and 3 dimensions.
This brief phenomenological study of thermodynamic entropy components served to give a physical grounding for information measures and what they reveal in spin systems on various lattice topologies. The results suggest many avenues of potential research. One topic to explore is the behavior of b in the Potts model, which switches from exhibiting a second-order phase transition in the magnetization to a first-order transition when the number of spin states exceeds four. Another setting of particular interest is to study the behavior of b in a frustrated system, such as the antiferromagnetic Ising model on a triangular lattice, paralleling Ref. [44]. Finally, as alluded to above, although beyond the present scope of this work, a next step is to consider informational measures for spins on arbitrary graphs, with the goal of providing insight into the roles that different nodes play in information processing and storage in complex dynamical networks.
expression we derive the joint probability of the central spin and its neighbors.
Let g n (σ 0 ) → g(σ 0 ) in the limit of n → ∞ shells. Then, the spins not close to the Cayley tree leaves behave like those on the Bethe lattice with the same coordination number. In this limit, by explicitly accounting for the bonds between the central spin and its neighbors, we can write the partition function as: . . . [1,k] σ 0 σ i   i∈ [1,k] [g(σ where the product is essentially over all branches 1 step away from the central spin, and σ i is the spin to which it is anchored. From this, we obtain the joint probability of the central spin σ 0 and its k neighbors as: p(σ 0 ,σ 1 , . . . , σ k ) = exp   βJ i∈ [1,k] σ 0 σ i   i∈ [1,k] [g(σ i )] k−1 /Z .
Above T c , Eq. (A2) has only one root: x = 1, which is stable. Below T c , however, there are there are 3 roots: x 0 , 1, x 0 −1 where x 0 < 1. The stable roots are x 0 and x 0 −1 . Setting x = x 0 provides the distribution p(σ 0 , σ 1 , . . . , σ k ) where the symmetry breaking prefers spins aligned up and where setting x = x 0 −1 provides the distribution with the symmetry broken such that downward spins are preferred.