{{Short description|Model in statistical mechanics}} {{no footnotes|date=June 2024}} thumb|Phase diagram of random subcube model.
In statistical mechanics, the '''random-subcube model (RSM)''' is an exactly solvable model that reproduces key properties of hard constraint satisfaction problems (CSPs) and optimization problems, such as geometrical organization of solutions, the effects of frozen variables, and the limitations of various algorithms like decimation schemes.
The RSM consists of a set of ''N'' binary variables, where solutions are defined as points in a hypercube. The model introduces clusters, which are random subcubes of the hypercube, representing groups of solutions sharing specific characteristics. As the density of constraints increases, the solution space undergoes a series of phase transitions similar to those observed in CSPs like random k-satisfiability (k-SAT) and random k-coloring (k-COL). These transitions include clustering, condensation, and ultimately the unsatisfiable phase where no solutions exist.
The RSM is equivalent to these real CSPs in the limit of large constraint size. Notably, it reproduces the cluster size distribution and freezing properties of k-SAT and k-COL in the large-k limit. This is similar to how the random energy model is the large-p limit of the p-spin glass model.
== Setup ==
=== Subcubes ===
There are <math>N</math> particles. Each particle can be in one of two states <math>-1, +1</math>.
The state space <math>\{-1, +1\}^N</math> has <math>2^N</math> states. Not all are available. Only those satisfying the constraints are allowed.
Each constraint is a subset <math>A_i</math> of the state space. Each <math>A_i</math> is a "subcube", structured like <math display=block>A_i = \prod_{j \in 1:N} A_{ij}</math> where each <math>A_{ij}</math> can be one of <math>\{-1\}, \{+1\}, \{-1, +1\}</math>.
The available states is the union of these subsets: <math display="block">S = \cup_i A_i </math>
=== Random subcube model ===
Each random subcube model is defined by two parameters <math>\alpha, p \in (0, 1)</math>.
To generate a random subcube <math>A_i</math>, sample its components <math>A_{ij}</math> IID according to <math display=block> \begin{aligned} Pr(A_{ij} &= \{-1\}) &= p/2 \\ Pr(A_{ij} &= \{+1\}) &= p/2 \\ Pr(A_{ij} &= \{-1, +1\}) &= 1-p \end{aligned} </math>
Now sample <math>2^{(1-\alpha)N}</math> random subcubes, and union them together.
=== Entropies ===
The entropy density of the <math>r</math>-th cluster in bits is <math display=block>s_r := \frac 1N \log_2 |A_r|</math>
The entropy density of the system in bits is <math display=block>s := \frac 1N \log_2 |\cup_r A_r|</math>
== Phase structure ==
=== Cluster sizes and numbers ===
Let <math>n(s)</math> be the number of clusters with entropy density <math>s</math>, then it is binomially distributed, thus <math display="block"> \begin{aligned} E[n(s)] &= 2^{(1-\alpha)N} P \to 2^{N\Sigma(s) + o(N)} \\ Var[n(s)] &= 2^{(1-\alpha)N} P(1-P) \\ \frac{Var[n(s)]}{E[n(s)]^2} &\to 2^{-N\Sigma(s)} \end{aligned} </math> where <math display=block> \begin{aligned} P &:= \binom{N}{sN}p^{(1-s)N}(1-p)^{sN}, \\ \Sigma(s) &:= 1-\alpha - D_{KL}(s \| 1-p) \\ D_{KL}(s \| 1-p) &:= s\log_2\frac{s}{1-p} + (1-s) \log_2\frac{1-s}{p} \end{aligned} </math>
By the Chebyshev inequality, if <math>\Sigma > 0</math>, then <math>n(s)</math> concentrates to its mean value. Otherwise, since <math>E[n(s)] \to 0</math>, <math>n(s)</math> also concentrates to <math>0</math> by the Markov inequality.
Thus, <math display="block">n(s) \to \begin{cases} 2^{N\Sigma(s) + o(N)} \quad &\text{if }\Sigma(s) > 0\\ 0 \quad &\text{if }\Sigma(s) < 0 \end{cases}</math> almost surely as <math>N\to\infty</math>.
When <math>\Sigma = 0</math> exactly, the two forces exactly balance each other out, and <math>n(s)</math> does not collapse, but instead converges in distribution to the Poisson distribution <math>Poisson(1)</math> by the law of small numbers.
=== Liquid phase === For each state, the number of clusters it is in is also binomially distributed, with expectation<math display="block">2^{(1-\alpha)N}(1-p/2)^N = 2^{N(\log_2(2-p) - \alpha)}</math>
So if <math>\alpha < \log_2(2-p)</math>, then it concentrates to <math>2^{N(\log_2(2-p) - \alpha)}</math>, and so each state is in an exponential number of clusters.
Indeed, in that case, the probability that ''all'' states are allowed is<math display="block">[1-[1-(1 - p/2)^N]^{2^{(1-\alpha) N}}]^{2^N}\sim e^{-e^{-2^{N(\log_2(2-p) - \alpha)} + N\ln 2}} \to 1</math>
Thus almost surely, all states are allowed, and the entropy density is 1 bit per particle.
=== Clustered phase ===
If <math>\alpha > \alpha_d := \log_2(2-p)</math>, then it concentrates to zero exponentially, and so most states are not in any cluster. Those that do are exponentially unlikely to be in more than one. Thus, we find that almost all states are in zero clusters, and of those in at least one cluster, almost all are in just one cluster. The state space is thus roughly speaking the disjoint union of the clusters.
Almost surely, there are <math>n(s) = 2^{N\Sigma(s)}</math> clusters of size <math>2^{Ns}</math>, therefore, the state space is dominated by clusters with optimal entropy density <math>s^* = \arg \max_s (\Sigma (s) + s)</math>.
Thus, in the clustered phase, the state space is almost entirely partitioned among <math>2^{N\Sigma(s^*)}</math> clusters of size <math>2^{Ns^*}</math> each. Roughly, the state space looks like exponentially many equally-sized clusters.
=== Condensation phase ===
Another phase transition occurs when <math>\Sigma(s^*) = 0</math>, that is,<math display="block">\alpha = \alpha_c := \frac{p}{(2-p)}+\log _2(2-p)</math>When <math>\alpha > \alpha_c</math>, the optimal entropy density becomes unreachable, as there almost surely exists zero clusters with entropy density <math>s^*</math>. Instead, the state space is dominated by clusters with entropy close to <math>s_c</math>, the larger solution to <math>\Sigma(s_c) = 0</math>.
Near <math>s_c</math>, the contribution of clusters with entropy density <math>s = s_c - \delta</math> to the total state space is <math display="block">\underbrace{2^{Ns}}_{\text{size of clusters}} \times \underbrace{2^{N\Sigma(s)}}_{\text{number of clusters}} = 2^{N(s + \Sigma(s))} = 2^{N(s_c - \delta - \Sigma'(s_c)\delta)}</math> At large <math>N</math>, the possible entropy densities are <math>s_c, s_c - 1/N, s_c - 2/N, \dots </math>. The contribution of each is <math display="block">2^{Ns_c}, 2^{Ns_c}2^{-(1+\Sigma'(s_c))}, 2^{Ns_c}2^{-2(1+\Sigma'(s_c))}, \dots</math>
We can tabulate them as follows:
{| class="wikitable" |+ ! <math>s</math> ! <math>s_c</math> ! <math>s_c - 1/N</math> ! <math>s_c - 2/N</math> |- | #clusters | <math>1</math> | <math>2^{\Sigma'(s_c)}</math> | <math>2^{2\Sigma'(s_c)}</math> |- | contributes | <math>\frac{1}{1 - 2^{-( 1 + \Sigma'(s_c))}}</math> | <math>\frac{2^{-( 1 + \Sigma'(s_c))}}{1 - 2^{-( 1 + \Sigma'(s_c))}}</math> | <math>\frac{2^{-2( 1 + \Sigma'(s_c))}}{1 - 2^{-( 1 + \Sigma'(s_c))}}</math> |}
Thus, we see that for any <math>\epsilon > 0</math>, at <math>N \to \infty</math> limit, over <math>1-\epsilon</math> of the total state space is covered by only a finite number of clusters. The state space looks partitioned into clusters with exponentially decaying sizes. This is the condensation phase.
=== Unsatisfiable phase === When <math>\alpha > 1</math>, the number of clusters is zero, so there are no states.
== Extensions == {{Stub section|date=June 2024}} The RSM can be extended to include energy landscapes, allowing for the study of glassy behavior, temperature chaos, and the dynamic transition.
== See also ==
==References==
* {{Cite journal |last1=Mora |first1=Thierry |last2=Zdeborová |first2=Lenka |date=2008-06-01 |title=Random Subcubes as a Toy Model for Constraint Satisfaction Problems |url=https://doi.org/10.1007/s10955-008-9543-x |journal=Journal of Statistical Physics |language=en |volume=131 |issue=6 |pages=1121–1138 |doi=10.1007/s10955-008-9543-x |issn=1572-9613|arxiv=0710.3804 |bibcode=2008JSP...131.1121M }} * {{cite journal |last=Zdeborová |first=Lenka |title=Statistical Physics of Hard Optimization Problems |date=2008-06-25 |arxiv=0806.4112 |bibcode=2008PhDT.......107Z |journal=Acta Physica Slovaca |volume=59 |issue=3 |pages=169–303}} * {{Cite book |last1=Mézard |first1=Marc |url= |title=Information, Physics, and Computation |last2=Montanari |first2=Andrea |date=2009-01-22 |publisher=Oxford University Press |isbn=978-0-19-154719-5 |language=en}} {{Reflist}}
Category:Statistical mechanics