Going after the k-SAT Threshold

Amin Coja-Oghlan, Konstantinos Panagiotou

Introduction

Since the early 2000s physicists have developed a sophisticated but highly non-rigorous technique called the “cavity method” for the study of random constraint satisfaction problems. This method allowed them to put forward a very detailed conjectured picture according to which various phase transitions affect both computational and structural properties of random CSPs. In addition, the cavity method has inspired new message passing algorithms called Belief/Survey Propagation guided decimation. Over the past few years there has been significant progress in turning bits and pieces of the physics picture into rigorous theorems. Examples include results on the interpolation method or the geometry of the solution space and their algorithmic implications .

In spite of this progress, substantial gaps remain. Perhaps most importantly, in most random CSPs the threshold for the existence of solutions is not known precisely. In the relatively simple case of the random kk-NAESAT (“Not-All-Equal-Satisfiability”) problem the difference between the best current lower and upper bounds is as tiny as 2−Ω(k)2^{-\Omega(k)} . By contrast, in random graph kk-coloring, a problem already studied by Erdős and Rényi in the 1960s, the best current bounds differ by Θ(ln⁡k)\Theta(\ln k) . Hence, the difference is unbounded in terms of the number of colors. Even worse, in random kk-SAT the gap is as big as Θ(k)\Theta(k) . Yet random kk-SAT is probably the single most important example of a random CSP, not least due to the great amount of experimental and algorithmic work conducted on it (e.g., ).

The reason for the large gap in random kk-SAT is that the satisfiability problem lacks a certain symmetry property. This property is vital to the current rigorous proof methods, particularly the second moment method, on which most of the previous work is based (e.g., ). More precisely, in random graph coloring the different colors all play the exact same role: for any proper coloring of a graph, another proper coloring can be obtained by simply permuting the color classes (e.g., color all red vertices blue and vice versa). Similarly, in kk-NAESAT, where the requirement is that in each clause at least one literal must be true and at least one false, the binary inverse of any NAE-solution is a NAE-solution as well. By contrast, in kk-SAT there is an inherent asymmetry between the Boolean values ‘true’ and ‘false’.

As has been noticed in prior work , the second moment method is fundamentally ill-posed to deal with such asymmetries. Roughly speaking, the second moment method is based on the assumption that in a random CSP instance, two randomly chosen solutions are perfectly uncorrelated. But in random kk-SAT, this is simply not the case. Indeed, suppose that a variable xx appears much more often positively than negatively throughout the formula. Then it seems reasonable to expect that most satisfying assignments set xx to ‘true’, thereby satisfying all clauses where xx appears positively. More generally, define the majority vote σmaj\sigma_{maj} to be the assignment that sets variable xx to true if it appears more often positively than negatively, and to false otherwise. Then we expect that the satisfying assignments of a random formula “gravitate toward” σmaj\sigma_{maj}. Unfortunately, the correlations among satisfying assignments induced by this drift toward σmaj\sigma_{maj} doom the second moment method. Previously this issue was sidestepped by symmetrizing the problem artificially . But this inevitably leaves a Θ(k)\Theta(k) gap.

The main contribution of the present work is a new asymmetric second moment method that enables us to tackle this problem head on. A key feature of this method is that we harness the Belief Propagation calculation from physics, called the “replica symmetric case” of the cavity method in physics jargon. We are going to employ Belief Propagation directly as an “educated guess” in the design the random variable upon which our proof is based in order to quantify how much a typical satisfying assignment leans toward σmaj\sigma_{maj}.

This is in contrast to most prior work on the subject, where individual statements hypothesized on the basis of physics arguments were proved via completely different methods (with the notable exception of the interpolation technique ). Hence, we view the present work as a pivotal step in the long-term effort of providing a rigorous foundation for the physicists’ cavity method. In fact, the general approach developed here does not hinge on particular properties of the kk-SAT problem, and thus we expect that the technique will extend to other asymmetric problems as well. Examples include not only other random CSPs that are asymmetric per se, but also instances of random problems that arise at intermediate steps of message passing algorithms such as Belief/Survey Propagation guided decimation, even if the initial problem is symmetric. In particular, we believe that getting a handle on asymmetric problems is a necessary step to analyze such message passing algorithms accurately.

To state our results precisely, we let k≥3,n>0k\geq 3,n>0 be integers and we let V={x1,…,xn}V=\left\{{x_{1},\ldots,x_{n}}\right\} be a set of nn Boolean variables. Further, let \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}=\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{k}(n,m) denote a Boolean formula with mm clauses of length kk over the variables VV chosen uniformly at random among all (2n)km(2n)^{km} such formulas. Let r=m/nr=m/n denote the density. We say that an event occurs with high probability (‘w.h.p.’) if its probability tends to 11 as n→∞n\rightarrow\infty.

where ok(1)o_{k}(1) hides a term that tends to 00 for large kk. The best prior lower bound is due to Achlioptas and Peres , who used a “symmetric” second moment argument to show

The bounds (1) and (2) leave an additive gap of k⋅ln⁡22+12+ok(1)k\cdot\frac{\ln 2}{2}+\frac{1}{2}+o_{k}(1), i.e., the gap is unbounded in terms of kk.

There is εk=ok(1)\varepsilon_{k}=o_{k}(1) such that

Apart from the quantitative improvement, the main point of this paper is that we manage to solve the problem of asymmetry in random CSPs for the first time. To explain this point, we start by discussing what we mean by asymmetry and how it derails the second moment method. That this is so was already intuited in . In the next section, we are going to verify and elaborate on those discussions.

Asymmetry and the second moment method

The second moment method. In general, the second moment method works as follows. Suppose that Z=Z(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}) is a non-negative random variable such that Z>0Z>0 only if Φ\textstyle\Phi is satisfiable. Moreover, suppose that for some density r>0r>0 there is a number C=C(k)>0C=C(k)>0 that may depend on kk but not on nn such that

Hence, we “just” need to find a random variable that satisfies (5). Let \mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}) denote the set of satisfying assignments; then certainly Z=\left|{\mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}})}\right| is the most obvious choice. However, this “vanilla” second moment argument turns out to fail spectacularly. We need to understand why.

Asymmetry and the majority vote. The origin of the problem is that kk-SAT is asymmetric in the following sense. Suppose that all we know about the random formula Φ\textstyle\Phi is for each variable xx the number dxd_{x} of times that xx appears as a positive literal in the formula, and the number d¬xd_{\neg x} of negative occurrences. Then our best stab at constructing a satisfying assignment seems to be the “majority vote” assigment σmaj\sigma_{maj} where we set xx to true if dx>d¬xd_{x}>d_{\neg x} and to false otherwise. Indeed, by maximizing the total number of true literal occurrences, of which a satisfying assignment must put one in every clause, σmaj\sigma_{maj} also maximizes the probability of being satisfiable.

Our proof of Theorem 1.1 allows us to formalize this observation, thereby verifying a conjecture from . Let \mboxdist(⋅,⋅)\mbox{dist}(\cdot,\cdot) denote the Hamming distance.

Hence, the average Hamming distance of \sigma\in\mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}) from σmaj\sigma_{maj} is strictly smaller than n/2n/2, i.e., the set \mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}) is “skewed toward” σmaj\sigma_{maj} w.h.p.

This asymmetry dooms the second moment method. To see why, let

Let us highlight this tradeoff, as it is characteristic of the kind of trouble that asymmetry causes. For ξ>0\xi>0 independent of nn but sufficiently small it turns out that for a certain constant c>0c>0,

That is, the probability is exponentially small but, like in the Chernoff bound, the exponent is a quadratic function of ξ\xi. By comparison, increasing the majority weight by ξ\xi boosts the expected number of satisfying assignments by a linear exponential factor: there is c′>0c^{\prime}>0 such that

For any k≥3k\geq 3 and r>2k/kr>2^{k}/k we have

In summary, the drift toward σmaj\sigma_{maj} and the resulting fluctuations of the majority weight induce a tremendous source of variance, derailing the “vanilla” second moment argument.

A quick fix? We saw that to make an asymmetric second moment argument work, we need to rule out fluctuations of the majority weight. A sensible way of implementing this is by actually fixing the entire vector \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}=(d_{x},d_{\neg x})_{x\in V} that counts the positively/negatively occurrences of each variable. More precisely, given a non-negative integer vector \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}=(d_{x},d_{\neg x})_{x\in V} with ∑x∈Vdx+d¬x=km\sum_{x\in V}d_{x}+d_{\neg x}=km let \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}} denote a uniformly random kk-CNF in which each variable xx appears dxd_{x} times positively and d¬xd_{\neg x} times negatively. Then we can split the generation of a random formula Φ\textstyle\Phi into two steps:

First, choose the occurrence vector d\textstyle d randomly from the “correct” distribution D\textstyle D.

Then, choose a random formula \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}.

The “correct” D\textstyle D is as follows. Let {\mathchoice{\mbox{\boldmath\displaystyle e}}{\mbox{\boldmath\textstyle e}}{\mbox{\boldmath\scriptstyle e}}{\mbox{\boldmath\scriptscriptstyle e}}}=(e_{x},e_{\neg x})_{x\in V} be a family of independent Poisson variables with mean kr/2kr/2 each. Moreover, let E{\cal E} be the event that ∑x∈Vex+e¬x=km\sum_{x\in V}e_{x}+e_{\neg x}=km. Let D\textstyle D be the conditional distribution of e\textstyle e given E{\cal E}. Then standard arguments show that the outcome of first choosing d\textstyle d and then \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}} is exactly the uniformly random Φ\textstyle\Phi.

The point of generating Φ\textstyle\Phi in two steps as above is that given the outcome d\textstyle d of the first step, the majority weight is fixed. Hence, if we could show that given a “typical” d\textstyle d, the second moment succeeds for |\mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}})| we would obtain a lower bound on rk−SATr_{k-SAT}. Unfortunately, matters are not so simple.

The explanation for this is that even if we fix d\textstyle d, various other types of fluctuations remain, turning \left|{\mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}})}\right| into a “lottery”. For instance, even given d\textstyle d the number of clauses that are unsatisfied under σmaj\sigma_{maj} fluctuates. Hence, the inherent asymmetry of kk-SAT puts not only the majority weight but also various other parameters on a slippery slope. What we need is a way of controlling all these fluctuations simultaneously. We will present our solution in Section 5.

Catching the kk-SAT threshold? Before we come to that, let us discuss what it would take to eliminate the (small but non-zero) gap left by Theorem 1.1, i.e., how far we are from “catching” the kk-SAT threshold. The physicists’ cavity method comes in two installments. The (relatively speaking) simpler “replica symmetric” version is based on Belief Propagation. Theorem 1.1 provides a rigorous proof of the best possible bound on the kk-SAT threshold that can be obtained from this version of the cavity method (up to possibly the precise error term εk\varepsilon_{k}) .

In we managed to prove rigorously that the 1RSB prediction for the random kk-NAESAT threshold is correct (up to an additive 2−Ω(k)2^{-\Omega(k)}). However, depends heavily on the fact that kk-NAESAT is symmetric. While it would be very interesting to combine the merits of the present paper with those of , this appears to be quite challenging. Thus, putting the 1RSB calculation for random kk-SAT on a rigorous foundation remains an important open problem. That said, we believe that any such attempt would need to build upon the techniques developed in this paper.

Related work

Also in random kk-XORSAT (random linear equations mod 2) the threshold for the existence of solutions is known precisely . The proof relies on computing the second moment of the number of solutions (after the instance has been stripped down to a suitable core). In contrast to random kk-SAT, the random kk-XORSAT problem is symmetric (cf. Remark 5.5 below), albeit in a more subtle way than kk-NAESAT.

Other problems where the second moment method succeeds are symmetric as well. Pioneering the use of the second moment method in random CSPs, Achlioptas and Moore computed the random kk-NAESAT threshold within an additive 1/21/2. By enhancing this argument with insights from physics this gap can be narrowed to a mere 2−Ω(k)2^{-\Omega(k)} . Moreover, the best current bounds on the random (hyper)graph kk-colorability thresholds are based on “vanilla” second moment arguments as well . In summary, in all the previous second moment arguments, the issue of asymmetry either did not appear at all by the nature of the problem , or it was sidestepped .

The best current algorithms for random kk-SAT find satisfying assignments w.h.p. for densities up to 1.817⋅2k/k1.817\cdot 2^{k}/k (better for small kk) resp. 2kln⁡(k)/k2^{k}\ln(k)/k (better for large kk) , a factor of Θ(k/ln⁡k)\Theta(k/\ln k) below the satisfiability threshold. By comparison, the Lovász Local Lemma and its algorithmic version succeed up to r=Θ(2k/k2)r=\Theta(2^{k}/k^{2}) .

Apart from experimental work , very little is known about the physics-inspired message passing algorithms (“Belief/Survey Propagation guided decimation”) . The most basic variant of Belief Propagation guided decimation is known to fail w.h.p. on random formulas if r>c⋅2k/kr>c\cdot 2^{k}/k for some constant c>0c>0 . However, it is conceivable that Survey Propagation and/or other variants of Belief Propagation perform better.

Preliminaries

We shall make repeated use of the following local limit theorem for the sums of independent random variables, see and .

where ζ\zeta and ξ\xi are the solutions to the equations

From this we can rather easily derive the following well-known statement about the rate function of the binomial distribution.

If p,qp,q remain fixed as n→∞n\rightarrow\infty, then

The following form of the chain rule will prove useful.

Let g:Ra→Rbg:\mathbf{R}^{a}\rightarrow\mathbf{R}^{b} and f:Rb→Rf:\mathbf{R}^{b}\rightarrow\mathbf{R} be of class C2C^{2}, i.e, with continuous second derivatives. Then for any x0∈Rax_{0}\in\mathbf{R}^{a} and with y0=g(x0)y_{0}=g(x_{0}) we have for any i,j∈[a]i,j\in\left[{a}\right]

Finally, we need the following version of the inverse function theorem that states under which conditions a given system of equations can be solved around a specific point uu. A detailed exposition can be found in .

Let U⊂RhU\subset\mathbf{R}^{h} be open and let f∈C1(U)f\in C^{1}(U). Assume that u∈Uu\in U and λ>0\lambda>0 are such that

Then for each y∈Rhy\in\mathbf{R}^{h} such that ∥y−f(u)∥≤λ/2\left\|{y-f(u)}\right\|\leq\lambda/2 there is precisely one x∈Rhx\in\mathbf{R}^{h} such that ∥x−u∥≤r\left\|{x-u}\right\|\leq r and f(x)=yf(x)=y. Furthermore, the inverse map f−1f^{-1} is C1C^{1} on {x∈Rh:∥x−u∥2<λ}\left\{{x\in\mathbf{R}^{h}:\left\|{x-u}\right\|_{2}<\lambda}\right\}, and Df−1(x)=(Df(x))−1Df^{-1}(x)=(Df(x))^{-1} on this set.

to denote the fact that ∥ξ−η∥∞≤O(1/n)\left\|{\xi-\eta}\right\|_{\infty}\leq O(1/n).

for some sequence εk=ok(1)\varepsilon_{k}=o_{k}(1) that tends to 00 sufficiently slowly.

The random variable

Our goal is to make the second moment method work for a random variable that counts “asymmetric” satisfying assignments. In this section, we develop this random variable. The starting point, and the key ingredient, is simply a map p:Z→[0,1]p:\mathbf{Z}\rightarrow\left[{0,1}\right]. For the sake of clarity, we start by setting up the framework for generic maps pp; below we will use the Belief Propagation formalism to pick the “optimal” pp.

The idea is that pp prescribes how strongly the assignments that we work with lean toward the majority vote. Informally speaking, we are going to work with assignments such that a variable xx that occurs dxd_{x} times positively and d¬xd_{\neg x} times negatively has a p(dx−d¬x)p(d_{x}-d_{\neg x}) chance of being set to ‘true’. Before we give a formal definition, we need to fix the number of times that each variable appears positively or negatively.

Fixing the majority weight. As we saw in Section 2, in order to make the second moment argument work, we need to rule out fluctuations of the majority weight. To achieve this, we follow the strategy outlined in Section 2. That is, we are going to work with formulas \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}} with a given vector \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}=(d_{x},d_{\neg x})_{x\in V} of occurrence counts, where each variable xx appears precisely dxd_{x} times positively and d¬xd_{\neg x} times negatively. As in Section 2, we let D\textstyle D denote the (conditional Poisson) distribution over sequences d\textstyle d such that first choosing d\textstyle d from D\textstyle D and then generating \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}} is equivalent to choosing a kk-CNF Φ\textstyle\Phi uniformly at random.

Fixing the marginals. Now, fix one such vector d\textstyle d. Then the map p:Z→[0,1]p:\mathbf{Z}\rightarrow\left[{0,1}\right] induces a map p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}} from the set L={x,¬x:x∈V}L=\left\{{x,\neg x:x\in V}\right\} of literals to [0,1]\left[{0,1}\right] in the natural way: we let

The idea is that, given d\textstyle d, we should set variable xx to ‘true’ with probability p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}(x).

To formalize this, we call p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}(l) the p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}-type of the literal ll. Let \mathcal{T}=\mathcal{T}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}=\left\{{p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}(l):l\in L}\right\} be the set of all possible p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}-types. We say that σ:V→{0,1}\sigma:V\rightarrow\left\{{0,1}\right\} has p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}-marginals if for any type t\in\mathcal{T}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}} we have

i.e., among all occurrences of literals of type tt, a tt fraction is true under σ\sigma. This definition captures the above idea that variable xx has a p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}(x) chance of being ‘true’.

Given \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}} there is a simple way of generating the random formula \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}}}. Namely, create dld_{l} clones of each literal ll, and put all the clones of a given p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}-type on a pile. Then the formula \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}}} is simply the result of matching the clones on the type tt pile randomly to all the clauses where a literal of type tt is required.

As in those papers, the problem admits a remarkably simple solution: let us call an assignment σ\sigma good in \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}}} if

Together with Paley-Zygmund (5), Theorem 5.1 shows that with d\textstyle d chosen from D\textstyle D and m\textstyle m chosen from \mathchoice{\mbox{\boldmath\displaystyle M}}{\mbox{\boldmath\textstyle M}}{\mbox{\boldmath\scriptstyle M}}{\mbox{\boldmath\scriptscriptstyle M}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}} w.h.p.

Guessing the marginals. For a set ∅≠S⊂{0,1}V\emptyset\neq S\subset\left\{{0,1}\right\}^{V} and a variable xx we define the SS-marginal of xx as

The definition of ‘p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}-judicious’ is guided by the idea that p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}(x) should prescribe the marginal of xx in the set of all p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}-judicious satisfying assignments. Hence, in order to make the set of p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}-judicious assignments as good an approximation of the entire set of satisfying assignments as possible, we better pick pp so that p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}(x) is a good approximation to the actual marginal \mu_{\mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}})}(x) of xx in the set of all satisfying assignments. The problem is that, because of the asymmetry of the kk-SAT problem, these marginals are highly non-trivial quantities. Indeed, on general formulas Φ\Phi the marginals μS(Φ)(x)\mu_{\mathcal{S}(\Phi)}(x) are #P\#P-hard to compute.

We observe that (20) is in line with the notion that \mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}) is “skewed toward” σmaj\sigma_{maj}. Indeed, the conjecture quantifies how much so. Motivated by Conjecture 5.2, we define

Under the distribution D\textstyle D, the random variables dx,d¬xd_{x},d_{\neg x} are asymptotically independent Poisson with mean kr/2kr/2 (cf. Section 2). Therefore,

Belief Propagation actually leads to a stronger prediction than Conjecture 5.2. Namely, it yields a conjecture for \mu_{\mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}})}(x) up to an additive error then tends to 00 as n→∞n\rightarrow\infty. However, (a) this stronger conjecture is not in explicit form, and (b) it does not only depend on dx,d¬xd_{x},d_{\neg x}, but also on various other parameters. In any case, even a more accurate prediction would not yield a better constant than 32ln⁡2\frac{3}{2}\ln 2 in Theorem 1.1.

2 Typical degree sequences

We need to collect a few basic properties of the sequence d\textstyle d chosen from D\textstyle D. Let us call a sequence \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}=(d_{l})_{l\in L} of non-negative integers such that ∑l∈Ldl=km\sum_{l\in L}d_{l}=km a signed degree sequence. For a kk-CNF Φ\Phi let \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}\left({\Phi}\right)=(d_{l}\left({\Phi}\right))_{l\in L} denote the vector whose entry dl(Φ)d_{l}\left({\Phi}\right) is equal to the number of times that literal ll occurrs in Φ\Phi. Then \mathchoice{\mbox{\boldmath\displaystyle D}}{\mbox{\boldmath\textstyle D}}{\mbox{\boldmath\scriptstyle D}}{\mbox{\boldmath\scriptscriptstyle D}}=\mathchoice{\mbox{\boldmath\displaystyle D}}{\mbox{\boldmath\textstyle D}}{\mbox{\boldmath\scriptstyle D}}{\mbox{\boldmath\scriptscriptstyle D}}_{k}(n,m) is just the distribution of the signed degree sequence \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}).

Let d\textstyle d be a signed degree sequence. A kk-CNF Φ\Phi over VV is d\textstyle d-compatible if \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}(\Phi)=\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}. Thus,

is a uniformly random d\textstyle d-compatible kk-CNF.

For d\textstyle d chosen from D\textstyle D the following statements hold w.h.p.

∑x∈V(dx−d¬x)2∼km.\sum_{x\in V}(d_{x}-d_{\neg x})^{2}\sim km.

Proof. We use the following description of the distribution D\textstyle D. Let \mathchoice{\mbox{\boldmath\displaystyle e}}{\mbox{\boldmath\textstyle e}}{\mbox{\boldmath\scriptstyle e}}{\mbox{\boldmath\scriptscriptstyle e}}=(e_{l})_{l\in L} be a family of indepedent Po(kr/2){\rm Po}(kr/2) variables. Moreover, let E{\cal E} be the event that ∑l∈Lel=km\sum_{l\in L}e_{l}=km. It is well known that e\textstyle e given E{\cal E} has distribution D\textstyle D. Furthermore, a simple calculation based on Stirling’s formula yields

Furthermore, as ex,e¬xe_{x},e_{\neg x} are independent for any x∈Vx\in V, we have

Because e^l≤ln⁡2n\hat{e}_{l}\leq\ln^{2}n and the random variables {(e^x−e^¬x)2}x∈V\left\{{(\hat{e}_{x}-\hat{e}_{\neg x})^{2}}\right\}_{x\in V} are mutually independent, Azuma’s inequality yields

thereby proving the first claim. The second claim follows from the first by means of the Cauchy-Schwarz inequality: w.h.p.

Finally, the third assertion is immediate from the second. □\Box

Let d\textstyle d be chosen from D\textstyle D. Then w.h.p. the following is true.

Proof. We use the alternative description of D\textstyle D from the proof of Lemma 5.7. That is, \mathchoice{\mbox{\boldmath\displaystyle e}}{\mbox{\boldmath\textstyle e}}{\mbox{\boldmath\scriptstyle e}}{\mbox{\boldmath\scriptscriptstyle e}}=(e_{l})_{l\in L} is a family of indepedent Po(kr/2){\rm Po}(kr/2) variables, and E{\cal E} is the event that ∑l∈Lel=km\sum_{l\in L}e_{l}=km. Let λ=kr/2\lambda=kr/2. For any fixed set S⊂LS\subset L the random variable XS=∑l∈SelX_{S}=\sum_{l\in S}e_{l} has distribution Po(∣S∣λ){\rm Po}(|S|\lambda) (because the sum of two independent Poisson variables is Poisson). Therefore, letting μ=10∣S∣max⁡{kr,ln⁡(n/∣S∣)}\mu=10|S|\max\left\{{kr,\ln(n/|S|)}\right\}, we obtain from Stirling’s formula

For 1≤s≤2n1\leq s\leq 2n let X_{s}=\sum_{S:\left|{S}\right|=s}\mathchoice{\mbox{\boldmath\displaystyle 1}}{\mbox{\boldmath\textstyle 1}}{\mbox{\boldmath\scriptstyle 1}}{\mbox{\boldmath\scriptscriptstyle 1}}_{X_{S}>\mu}. Then (27) yields

because μ≥10sln⁡(n/s)\mu\geq 10s\ln(n/s). Thus, the first claim follows from (22) and the union bound.

To prove the second claim, we use Lemma 5.6. For S⊂LS\subset L we let YSY_{S} be the total number of occurrences of literals from SS in Φ\textstyle\Phi. Then YSY_{S} has distribution Bin(km,∣S∣/2n){\rm Bin}(km,|S|/2n) with mean ∣S∣kr/2|S|kr/2. By the Chernoff bound,

Hence, letting Y_{s}=\sum_{S:|S|=s}\mathchoice{\mbox{\boldmath\displaystyle 1}}{\mbox{\boldmath\textstyle 1}}{\mbox{\boldmath\scriptstyle 1}}{\mbox{\boldmath\scriptscriptstyle 1}}_{Y_{S}<kr|S|/3}, we get from (28) for s≥n2−0.8ks\geq n2^{-0.8k}

For any t∈Tt\in\mathcal{T} we let n(t)n(t) be the number of variables x∈Vx\in V such that p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}(x)=t.

Let d\textstyle d be chosen from D\textstyle D. Then w.h.p. for any type t∈Tt\in\mathcal{T} we have

Thus, the assertion follows by combining (22), (29) and (30). □\Box

For each t∈Tt\in\mathcal{T} we let π(t)\pi(t) denote the fraction of literal occurrences of pp-type tt, i.e.,

Since ∣L∣=O(1)\left|{\mathcal{L}}\right|=O(1) as n→∞n\rightarrow\infty by the construction of pp, the assertion follows from (31) and the union bound. □\Box

The first moment

Let ρ>32ln⁡2\rho>\frac{3}{2}\ln 2 be such that r=2kln⁡2−ρr=2^{k}\ln 2-\rho.

W.h.p. \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}} are such that

Let χ(z)=−zln⁡z−(1−z)ln⁡(1−z)\chi(z)=-z\ln z-(1-z)\ln(1-z) denote the entropy function. Then w.h.p. d\textstyle d is such that

Taylor expanding χ(z)\chi(z) around z=1/2z=1/2 and plugging in the definition (21) of pp, we obtain that w.h.p. d\textstyle d is such that

As a next step, we compute the probability of \sigma\in\mathcal{S}_{p}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}}}) for \sigma\in\mathcal{H}_{p}(\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}).

W.h.p. d\textstyle d, m\textstyle m are such that for any \sigma\in\mathcal{H}_{p}(\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}),

Let us defer the proof of Lemma 6.3, which is the core of the first moment computation, for a little while. Combining (32)–(34), we see that w.h.p. over the choice of \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}} we have

W.h.p. over the choice of \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}} we have

The proof of Lemma 6.4 is based on arguments developed in for analyzing the geometry of the set of satisfying assignments. Combining (35) and Lemma 6.4 yields Proposition 6.1.

2 Proof of Lemma 6.3

For any σ∈{0,1}V\sigma\in\left\{{0,1}\right\}^{V} with pp-marginals we have

The proof of Proposition 6.5 consists of two steps. We defer the proof of the following lemma to Section 6.3.

With the assumptions of Proposition 6.5 and with δ,δ′\delta,\delta^{\prime} defined by

Proof of Proposition 6.5. Let Δ=100k2kln⁡k\Delta=100k2^{k}\ln k and let δ,δ′\delta,\delta^{\prime} be as in Lemma 6.6. Using the alternative description of the distribution D\textstyle D from the proof of Lemma 5.7 and applying Azuma’s inequality, one can easily verify that w.h.p.

Similarly, invoking (36) once more, we see that w.h.p.

w.h.p. Thus, Proposition 6.5 is a direct consequence of Lemmas 5.9 and 6.6 and (37), (38). □\Box

3 Proof of Lemma 6.6

We begin by determining the number σ∈{0,1}V\sigma\in\left\{{0,1}\right\}^{V} with pp-marginals. The following is an easy consequence of Lemma 6.2.

W.h.p. for d\textstyle d chosen from D\textstyle D we have

Proof. This follows from Lemma 6.2 by Taylor expanding χ(⋅)\chi(\cdot) around 12\frac{1}{2}. □\Box

For any d\textstyle d-compatible formula \Phi\in\Gamma_{\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}}} we can define a map

Proof. Since the total number ∣L∣\left|{\mathcal{L}}\right| of clause types is bounded, the assertion follows from a repeated application of Lemma 4.1 (the local limit theorem). □\Box

From this point on we fix q\textstyle q as in Lemma 6.12.

where we used the approximation ln⁡(1+x)=x−12x2+O(x3)\ln(1+x)=x-\frac{1}{2}x^{2}+O(x^{3}). Thus, Lemma 5.10 yields

In the second moment calculation we will need to know that

Let δ,δ′>0\delta,\delta^{\prime}>0 be such that

Proof of Lemma 6.6. Lemma 6.6 is a direct consequence of Corollaries 6.7, 6.9, 6.11 and 6.15. □\Box

4 Proof of Lemma 6.4

Assume that m\textstyle m is feasible. Let Z\mathcal{Z} denote the number of good pp-satisfying assignments.

The proof of Proposition 6.16 is based on three lemmas.

Let d\textstyle d be chosen from D\textstyle D and let m\textstyle m be chosen from \mathchoice{\mbox{\boldmath\displaystyle M}}{\mbox{\boldmath\textstyle M}}{\mbox{\boldmath\scriptstyle M}}{\mbox{\boldmath\scriptscriptstyle M}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}.

Proof. This follows from a similar application of Markov’s inequality as in the proof of Lemma 5.6. □\Box

With the assumptions of Proposition 6.16 the random variable

The proof of Lemma 6.18 can be found in Section 6.5. Moreover, in Section 6.6 we prove the following.

Suppose that r≤2kln⁡2r\leq 2^{k}\ln 2. Let ξ=k2−k/2\xi=k2^{-k/2}. Let Z′′Z^{\prime\prime} be the number of pairs (\sigma,\tau)\in\mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}})^{2} such that

Finally, Proposition 6.16 follows immediately from Lemmas 6.17, 6.18 and 6.19.

5 Proof of Lemma 6.18

Let Φ\Phi be a kk-CNF and σ∈S(Φ)\sigma\in\mathcal{S}(\Phi). We say that a variable xx is ξ\xi-rigid in (Φ,σ)(\Phi,\sigma) if for any τ∈S(Φ)\tau\in\mathcal{S}(\Phi) with τ(x)≠σ(x)\tau(x)\neq\sigma(x) we have \mboxdist(σ,τ)≥ξn\mbox{dist}(\sigma,\tau)\geq\xi n. Let λ=kr/(2k−1)\lambda=kr/(2^{k}-1).

Proof. Fix an assignment σ∈{0,1}V\sigma\in\left\{{0,1}\right\}^{V}, say \sigma=\mathchoice{\mbox{\boldmath\displaystyle 1}}{\mbox{\boldmath\textstyle 1}}{\mbox{\boldmath\scriptstyle 1}}{\mbox{\boldmath\scriptscriptstyle 1}}. Then the number of clauses supported by each x∈Vx\in V is asymptotically Poisson with mean λ\lambda. Let Ex{\cal E}_{x} be the event that xx supports no more than 12 clauses. Then

The events (Ex)x∈V({\cal E}_{x})_{x\in V} are negatively correlated. Therefore, the total number XX of variables x∈Vx\in V for which Ex{\cal E}_{x} occurs is stochastically dominated by a binomial variable Bin(n,12k122−k){\rm Bin}(n,\frac{1}{2}k^{12}2^{-k}). Hence, the first assertion follows from Chernoff bounds.

Let us call a set S⊂VS\subset V self-contained if each variable in SS supports at least ten clauses that consist of variables in SS only. There is a simple process that yields a (possibly empty) self-contained set SS.

For each variable xx that supports at least one clause, choose such a clause CxC_{x} randomly.

Let RR be the set of all variables that support at least 12 clauses.

While there is a variable x∈Rx\in R that supports fewer than ten clauses \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{i}\neq C_{x} that consist of variables of RR only, remove xx from RR.

The clauses CxC_{x} will play a special role later.

Proof. Let σ∈{0,1}V\sigma\in\left\{{0,1}\right\}^{V} be an assignment, say \sigma=\mathchoice{\mbox{\boldmath\displaystyle 1}}{\mbox{\boldmath\textstyle 1}}{\mbox{\boldmath\scriptstyle 1}}{\mbox{\boldmath\scriptscriptstyle 1}}. Let QQ be the set of all variables that support fewer than 12 clauses. By Lemma 6.20 we may condition on ∣Q∣≤k122−kn|Q|\leq k^{12}2^{-k}n. Assume that ∣R∣≤(1−k15/2k)n|R|\leq(1-k^{15}/2^{k})n. Then there exists a set S⊂V∖(R∪Q)S\subset V\setminus(R\cup Q) of size 12k15n/2k≤S≤k15n/2k\frac{1}{2}k^{15}n/2^{k}\leq S\leq k^{15}n/2^{k} such that each variable in SS supports ten clauses that contain another variable from S∪QS\cup Q. With s=∣S∣/ns=|S|/n the probability of this event is bounded by

Hence, the expected number of set SS for which the aforementioned event occurs is bounded by

Let us call a variable xx is attached if xx supports a clause whose other k−1k-1 variables belong to RR.

Proof. Let F=V∖RF=V\setminus R. By Proposition 6.21 we may assume that ∣F∣≤nk15/2k|F|\leq nk^{15}/2^{k}. Therefore, for each of the “special” clause CxC_{x} that we reserved for each xx that supports at least one clause the probability of containing a variable from F∖{x}F\setminus\left\{{x}\right\} is bounded by

Furthermore, these events are independent (because the clauses CxC_{x} were disregarded in the construction of RR). Hence, the number of variables x∉Rx\not\in R that support at least one clause but that are not attached is dominated by Bin(∣F∣,3k162k){\rm Bin}(|F|,\frac{3k^{16}}{2^{k}}). The assertion thus follows from Chernoff bounds. □\Box

Let us call S⊂VS\subset V dense if each variable in SS supports at least ten clauses and at most 2k2k clauses such that at least ten of them feature another variable from SS.

For d\textstyle d chosen from D\textstyle D, m\textstyle m chosen from \mathchoice{\mbox{\boldmath\displaystyle M}}{\mbox{\boldmath\textstyle M}}{\mbox{\boldmath\scriptstyle M}}{\mbox{\boldmath\scriptscriptstyle M}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}} and any σ∈{0,1}V\sigma\in\left\{{0,1}\right\}^{V} the following holds w.h.p. Let A\mathcal{A} be the event that σ\sigma is a pp-satisfying assignment of \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}}}. Then

Due to negative correlation, in total we obtain

(The factor (2k10)∣S∣{{2k}\choose{10}}^{|S|} accounts for the number of ways to choose 1010 out of the at most 2k2k clauses that each variable in SS supports.)

For 0<s≤1/k50<s\leq 1/k^{5} let XsX_{s} be the number of sets SS of size ∣S∣=sn|S|=sn for which D(S)\mathcal{D}(S) occurs. Then

There are two cases to consider. First, if s≤ln⁡(n)/ns\leq\ln(n)/n, then the term in the brackets is clearly o(1)o(1). Second, if s≥ln⁡(n)/ns\geq\ln(n)/n, then we have the following bound. Since s≤smax⁡=2−0.99ks\leq s_{\max}=2^{-0.99k} and as x↦x9ln⁡10xx\mapsto x^{9}\ln^{10}x is monotonically increasing for x<0.1x<0.1, we have

Hence, the entire bracket is bounded by 2−k/22^{-k/2}. Summing over all possible ss and using Markov’s inequality completes the proof. □\Box

Let us call a variable x∈Vx\in V ξ\xi-rigid in σ∈S(Φ)\sigma\in\mathcal{S}(\Phi) if for any τ∈S(Φ)\tau\in\mathcal{S}(\Phi) with τ(x)≠σ(x)\tau(x)\neq\sigma(x) we have \mboxdist(σ,τ)≥ξn\mbox{dist}(\sigma,\tau)\geq\xi n.

W.h.p. for d\textstyle d chosen from D\textstyle D and for m\textstyle m chosen from \mathchoice{\mbox{\boldmath\displaystyle M}}{\mbox{\boldmath\textstyle M}}{\mbox{\boldmath\scriptstyle M}}{\mbox{\boldmath\scriptscriptstyle M}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}} the following is true. Let σ∈{0,1}V\sigma\in\left\{{0,1}\right\}^{V} and let A\mathcal{A} be the event that σ\sigma is a pp-satisfying assignment of \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}}}. Moreover, let YY be the number of variables that are not 2−0.99k2^{-0.99k}-rigid. Then

Proof. Let ξ=2−0.99k\xi=2^{-0.99k}. We condition on the event A\mathcal{A}. Consider a variable zz that is either attached or in RR. Let \tau\in\mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}}}) be such that τ(z)≠σ(x)\tau(z)\neq\sigma(x) and \mboxdist(σ,τ)<n/20.99k\mbox{dist}(\sigma,\tau)<n/2^{0.99k}. Because zz is attached or in RR, the set

is non-empty. Moreover, Δ\Delta is dense by the construction of RR. Thus, Lemma 6.23 shows that \mboxdist(σ,τ)≥∣Δ∣≥n/20.99k\mbox{dist}(\sigma,\tau)\geq\left|{\Delta}\right|\geq n/2^{0.99k} w.h.p. Hence, w.h.p. all zz that are either attached or in RR are ξ\xi-rigid.

Further, let R{\mathcal{R}} be the event that

no more than (1+1/k2)2−kn(1+1/k^{2})2^{-k}n variables support no clause at all and

at most n/(k22k)n/(k^{2}2^{k}) variables x∉Rx\not\in R that support at least one clause are not attached

Then Lemma 6.20 and Corollary 6.22 imply together with Proposition 6.5 that

Hence, the total number of vertices that either do not support a clause or that are not attached is bounded by (1+2/k2)2−kn(1+2/k^{2})2^{-k}n □\Box

Proof of Lemma 6.18. Suppose that r=2kln⁡2−cr=2^{k}\ln 2-c. By Proposition 6.5 we have

w.h.p. Now, assume that in \sigma\in\mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}}}) all but at most (1+2k−2)2−kn(1+2k^{-2})2^{-k}n variables are ξ\xi-rigid with ξ=2−0.99k\xi=2^{-0.99k}. If \tau\in\mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}}}) is such that \mboxdist(σ,τ)≤ξn\mbox{dist}(\sigma,\tau)\leq\xi n, then σ,τ\sigma,\tau agree on all ξ\xi-rigid variables of σ\sigma. Hence,

As c−ln⁡22+ok(1)>(1+ok(1))ln⁡2c-\frac{\ln 2}{2}+o_{k}(1)>(1+o_{k}(1))\ln 2 for c>32ln⁡2+εc>\frac{3}{2}\ln 2+\varepsilon and kk large enough, the assertion follows. □\Box

6 Proof of Lemma 6.19

By Markov’s inequality, it suffices to bound the expected number of paris (\sigma,\tau)\in\mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}) at the given Hamming distances. More precisely, let Zx\mathcal{Z}_{x} be the number of pairs (\sigma,\tau)\in\mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}) such that \mboxdist(σ,τ)/n=x\mbox{dist}(\sigma,\tau)/n=x. Let h(x)=−xln⁡x−(1−x)ln⁡(1−x)h(x)=-x\ln x-(1-x)\ln(1-x) and set

If k2−k≤x≤k−2k2^{-k}\leq x\leq k^{-2}, then 1−ln⁡x−k+k2x≤1−ln⁡k+1<01-\ln x-k+k^{2}x\leq 1-\ln k+1<0. Moreover, if k−2≤x≤(2k)−1k^{-2}\leq x\leq(2k)^{-1}, then 1−ln⁡x−k+k2x≤1+2ln⁡k−34k<01-\ln x-k+k^{2}x\leq 1+2\ln k-\frac{3}{4}k<0.

The last expression is negative for x<0.05x<0.05 (and kk not too small).

Hence, for 0.01≤x<12−k−20.01\leq x<\frac{1}{2}-k^{-2} we have h′(x)+q′(x)>0h^{\prime}(x)+q^{\prime}(x)>0. Thus, h(x)+q(x)+ln⁡2h(x)+q(x)+\ln 2 is monotonically increasing in this interval. Now, let x=12−εx=\frac{1}{2}-\varepsilon for k−2≤ε≤k2−k/2k^{-2}\leq\varepsilon\leq k2^{-k/2}. Then

The function h(x)h(x) satisfies h(1−y)=h(y)h(1-y)=h(y) for 0<y<1/20<y<1/2. Furthermore, q(x)q(x) is monotonically decreasing. Therefore, for any x≥12+k2−k/2x\geq\frac{1}{2}+k2^{-k/2} we have

In each case we have ln⁡2+h(x)+q(x)<0\ln 2+h(x)+q(x)<0. Thus, the assertion follows from (47) and Markov’s inequality.

The second moment

The overlap of two assignments σ,τ∈{0,1}V\sigma,\tau\in\left\{{0,1}\right\}^{V} is the vector

In words, O(σ,τ)\mathcal{O}(\sigma,\tau) captures the fraction of occurrences of literals of each type tt that are true under both σ,τ\sigma,\tau. Since \mathchoice{\mbox{\boldmath\displaystyle\sigma}}{\mbox{\boldmath\textstyle\sigma}}{\mbox{\boldmath\scriptstyle\sigma}}{\mbox{\boldmath\scriptscriptstyle\sigma}},\mathchoice{\mbox{\boldmath\displaystyle\tau}}{\mbox{\boldmath\textstyle\tau}}{\mbox{\boldmath\scriptstyle\tau}}{\mbox{\boldmath\scriptscriptstyle\tau}} are independent and have pp-marginals, we have

Set O∗=[t2]t∈T\mathcal{O}^{*}=\left[{t^{2}}\right]_{t\in\mathcal{T}}.

Let Z′′Z^{\prime\prime} be the number of pairs (σ,τ)(\sigma,\tau) of p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}-judicious satisfying assignments of \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}}} such that

Moreover, let Z′Z^{\prime} be the number of pairs (σ,τ)(\sigma,\tau) of p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}-judicious satisfying assignments of \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}}} such that

The proof of Proposition 7.1 can be found in Section 7.2. Let ZZ denote the number of pp-satisfying assignments of \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}. Furthermore, let Z\mathcal{Z} signify the number of good pp-satisfying assignments of \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}. In Section 8 we are going to establish the following.

Proof. Let YY be the number of pairs (σ,τ)(\sigma,\tau) of good pp-satisfying assignments of \mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}}} such that

Combining (50) with Proposition 7.1 and 7.2, we obtain for d\textstyle d chosen from D\textstyle D w.h.p.

The second part of Theorem 5.1 follows directly from Corollary 7.3.

2 Proof of Proposition 7.1

We begin by relating the overlap to the Hamming distance.

W.h.p. \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}} are such that for all pairs σ,τ∈{0,1}V\sigma,\tau\in\left\{{0,1}\right\}^{V} satisfying (48) we have

W.h.p. \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}} are such that for any σ,τ∈{0,1}V\sigma,\tau\in\left\{{0,1}\right\}^{V} that satisfy (48) and that have pp-marginals we have

Let (\hat{\mathchoice{\mbox{\boldmath\displaystyle\sigma}}{\mbox{\boldmath\textstyle\sigma}}{\mbox{\boldmath\scriptstyle\sigma}}{\mbox{\boldmath\scriptscriptstyle\sigma}}},\hat{\mathchoice{\mbox{\boldmath\displaystyle\tau}}{\mbox{\boldmath\textstyle\tau}}{\mbox{\boldmath\scriptstyle\tau}}{\mbox{\boldmath\scriptscriptstyle\tau}}}) denote a random pair chosen from this distribution.

Proposition 6.5 and Lemma 7.4 ensure that w.h.p. d\textstyle d is such that

In addition, let SS be the event that \max_{j\in\left[{k}\right]}\hat{\mathchoice{\mbox{\boldmath\displaystyle\sigma}}{\mbox{\boldmath\textstyle\sigma}}{\mbox{\boldmath\scriptstyle\sigma}}{\mbox{\boldmath\scriptscriptstyle\sigma}}}_{ij}=\max_{j\in\left[{k}\right]}\hat{\mathchoice{\mbox{\boldmath\displaystyle\tau}}{\mbox{\boldmath\textstyle\tau}}{\mbox{\boldmath\scriptstyle\tau}}{\mbox{\boldmath\scriptscriptstyle\tau}}}_{ij} for all i∈[m]i\in\left[{m}\right]. We claim that

Indeed, any dd-compatible formula Φ\Phi induces a pair (σ^∣Φ,τ^∣Φ)∈Ω^(\hat{\sigma}|_{\Phi},\hat{\tau}|_{\Phi})\in\hat{\Omega} defined by σ^ij∣Φ=σ(Φij)\hat{\sigma}_{ij}|_{\Phi}=\sigma(\Phi_{ij}), τ^ij∣Φ=τ(Φij)\hat{\tau}_{ij}|_{\Phi}=\tau(\Phi_{ij}). Clearly, the distribution of the random pair (\hat{\sigma}|_{\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}},\hat{\tau}|_{\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}}) is identical to the distribution of (\hat{\mathchoice{\mbox{\boldmath\displaystyle\sigma}}{\mbox{\boldmath\textstyle\sigma}}{\mbox{\boldmath\scriptstyle\sigma}}{\mbox{\boldmath\scriptscriptstyle\sigma}}},\hat{\mathchoice{\mbox{\boldmath\displaystyle\tau}}{\mbox{\boldmath\textstyle\tau}}{\mbox{\boldmath\scriptstyle\tau}}{\mbox{\boldmath\scriptscriptstyle\tau}}}) given BB.

Due to independence, the probability of the event SS is easy to compute. Indeed, with q=q10+q11q=q^{10}+q^{11} inclusion/exclusion yields

Using (52) and simplifying completes the proof. □\Box

Let λ>2−k\lambda>2^{-k} and t∈Tt\in\mathcal{T}. For d\textstyle d chosen from D\textstyle D the following is true w.h.p. Let H′′\mathcal{H}^{\prime\prime} be the set of all pairs σ,τ∈{0,1}V\sigma,\tau\in\left\{{0,1}\right\}^{V} such that ∣Ot(σ,τ)−1/4∣>λ.|\mathcal{O}_{t}(\sigma,\tau)-1/4|>\lambda. Then

If σ′,τ′′,σ′′,τ′′∈{0,1}V\sigma^{\prime},\tau^{\prime\prime},\sigma^{\prime\prime},\tau^{\prime\prime}\in\left\{{0,1}\right\}^{V} are such that there is a literal l0l_{0} with T(l0)=t\mathcal{T}(l_{0})=t such that σ′′(l)=σ′(l),τ′′(l)=τ′(l)\sigma^{\prime\prime}(l)=\sigma^{\prime}(l),\tau^{\prime\prime}(l)=\tau^{\prime}(l) for all l∉{l0,¬l0}l\not\in\left\{{l_{0},\neg l_{0}}\right\}, then

Therefore, by Azuma’s inequality for any λ>0\lambda>0 we have

where the last step follows from part 2 of Proposition 6.5. □\Box

Proof of Proposition 7.1. Let H′′H^{\prime\prime} be the set of pairs (σ,τ)(\sigma,\tau) such that

σ,τ\sigma,\tau satisfy (48) and have pp-marginals, and

∥O(σ,τ)−O∗∥∞>ξ\left\|{\mathcal{O}(\sigma,\tau)-\mathcal{O}^{*}}\right\|_{\infty}>\xi.

Then by Lemma 7.6 and the second part of Proposition 6.5 w.h.p. (over the choice of d\textstyle d) we have

Furthermore, by Lemma 7.5 w.h.p. (again over the choice of d\textstyle d) we have

Combining (54) and (55), we obtain that w.h.p. d\textstyle d is such that

Therefore, the definition of the distribution \mathchoice{\mbox{\boldmath\displaystyle M}}{\mbox{\boldmath\textstyle M}}{\mbox{\boldmath\scriptstyle M}}{\mbox{\boldmath\scriptscriptstyle M}}_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}} entails that w.h.p. d\textstyle d is such that

Thus, the assertion follows from Markov’s inequality. □\Box

Proof of Proposition 7.2

We keep the notation and the assumptions of Section 7.

For two assignments σ,τ\sigma,\tau and a formula Φ\Phi with signed degree distribution d\textstyle d we define a matrix

In Section 9 we are going to prove the following.

For any assignment σ\sigma with pp-marginals we have

Recall that O∗=(t2)t∈T\mathcal{O}^{*}=(t^{2})_{t\in\mathcal{T}}. In Section 8.3 we will prove the following.

W.h.p. \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}} are such that the following holds. For any O=(Ot)t∈T\mathcal{O}=(\mathcal{O}_{t})_{t\in\mathcal{T}} such that \left\|{\mathcal{O}-\frac{1}{4}\mathchoice{\mbox{\boldmath\displaystyle 1}}{\mbox{\boldmath\textstyle 1}}{\mbox{\boldmath\scriptstyle 1}}{\mbox{\boldmath\scriptscriptstyle 1}}}\right\|_{\infty}\leq 2\xi we have

In Section 8.4 we will show the following.

There exists a constant η>0\eta>0 such that w.h.p. \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}} are such that the following holds. For all O\mathcal{O} with \left\|{\mathcal{O}-\frac{1}{4}\mathchoice{\mbox{\boldmath\displaystyle 1}}{\mbox{\boldmath\textstyle 1}}{\mbox{\boldmath\scriptstyle 1}}{\mbox{\boldmath\scriptscriptstyle 1}}}\right\|_{\infty}\leq 2\xi we have

Recall that n(t)n(t) is the number of variables of type t∈Tt\in\mathcal{T}. In Section 10 we are going to prove the following.

W.h.p. \mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}},\mathchoice{\mbox{\boldmath\displaystyle m}}{\mbox{\boldmath\textstyle m}}{\mbox{\boldmath\scriptstyle m}}{\mbox{\boldmath\scriptscriptstyle m}} are such that the following holds. For all vectors λ=(λt)t∈T\lambda=(\lambda_{t})_{t\in\mathcal{T}} with ∥λ∥∞≤1/8\left\|{\lambda}\right\|_{\infty}\leq 1/8 we have

Proof of Proposition 7.2. Suppose that O∈[0,1]T\mathcal{O}\in\left[{0,1}\right]^{\mathcal{T}} satisfies ∥O−O∗∥∞≤ξ\left\|{\mathcal{O}-\mathcal{O}^{*}}\right\|_{\infty}\leq\xi. By Proposition 8.4 w.h.p.

For an assignment σ\sigma with pp-marginals let

Then by part 3 of Proposition 8.1, Corollary 8.4 and Corollary 6.9 we have

For a vector λ=(λt)t∈T\lambda=(\lambda_{t})_{t\in\mathcal{T}} let

Moreover, for c=c(k)>0c=c(k)>0 a sufficiently large number let Λ=cnZ≥0T\Lambda=\frac{c}{\sqrt{n}}\mathbf{Z}_{\geq 0}^{\mathcal{T}} be the positive T\mathcal{T}-dimensional grid scaled by a factor of c/nc/\sqrt{n}. In addition, let hh be the number of assignments σ\sigma with pp-marginals. Then by Proposition 8.5 and (59) there is a number ζ=ζ(k)>0\zeta=\zeta(k)>0 such that

It will be convenient to work with a different probability space. Namely, let Ω^\hat{\Omega} be the set of all pairs (σ^,τ^)(\hat{\sigma},\hat{\tau}) of 0/10/1 vectors

2 Proof of Proposition 8.2

the remaining entries of q\textstyle q are determined by (60)–(62). Then the following is immediate from the construction.

The last step follows from the local limit theorem for the multinomial distribution because

3 Proof of Corollary 8.3

observe that Pˉ\bar{P} depends on O\mathcal{O} but not on the specific choice of ω\omega. Then by Propositions 8.1 and 8.2

Summing over all possible overlap matrices ω\omega of assignments with pp-marginals, we get

4 Proof of Proposition 8.4

We claim that there exist numbers 0<ck<ck′0<c_{k}<c_{k}^{\prime} (independent of ω\omega) such that w.h.p. d\textstyle d is such that

Once more, the conditional probability that this random variable equals its expectation lies in [ck,3n−1/2,ck,4n−1/2]\left[{c_{k,3}n^{-1/2},c_{k,4}n^{-1/2}}\right] for certain ck,4>ck,3>0c_{k,4}>c_{k,3}>0. Hence, setting ck=ck,1ck,3c_{k}=c_{k,1}c_{k,3} and ck′=ck,2ck,4c_{k}^{\prime}=c_{k,2}c_{k,4}, we obtain (63).

Proof of Proposition 8.1

We keep the notation and the assumptions of Section 7.

In Section 9.2 we will establish the following.

For any ω\omega in the domain of P\mathcal{P} we have

Proof. This follows directly from Propositions 9.2 and 9.3 and Taylor’s formula. □\Box

Finally, in Section 9.6 we will show of Proposition 8.1 follows from Proposition 9.1 and Corollary 9.4.

2 Proof of Proposition 9.1

There exists a vector q\textstyle q such that

Proof. This follows from applying the inverse function theorem in a similar way as in the proof of Lemma 6.10. □\Box

In the rest of this section, we fix q\textstyle q as in Lemma 9.9.

By inclusion/exclusion, we obtain from (70) and (71) that

Due to (70) and (72), a repeated application of Lemma 4.1 (the local limit theorem) yields

Invoking Lemma 9.8 and using the large deviations principle for the binomial distribution (Lemma 4.2), we can easily determine the unconditional probability of B′∩C′B^{\prime}\cap C^{\prime}: we have

Finally, Proposition 9.1 follows from Fact 9.6 and Lemma 9.10.

3 Proof of Proposition 9.2

Proof. Equation (73) from the proof of Lemma 9.10 shows that

Furthermore, applying Azuma’s inequality just as in the previous paragraph, we find that

Let H(ω)H(\omega) be the number of pairs (σ^,τ^)∈Ω^(\hat{\sigma},\hat{\tau})\in\hat{\Omega} such (σ^,τ^)∈B′(\hat{\sigma},\hat{\tau})\in B^{\prime} and \hat{\mathchoice{\mbox{\boldmath\displaystyle\omega}}{\mbox{\boldmath\textstyle\omega}}{\mbox{\boldmath\scriptstyle\omega}}{\mbox{\boldmath\scriptscriptstyle\omega}}}(\hat{\sigma},\hat{\tau})=\omega. We claim that

This can be verified either by representing H(ω)H(\omega) as a product of binomial coefficients and applying Stirlings formula or, alternatively, by using (75). Indeed, assume that (81) is false. Then for small enough ε>0\varepsilon>0 there is δ>0\delta>0 such that for some ω′\omega^{\prime} with ∥ω′−ω∗∥∞∼ε\left\|{\omega^{\prime}-\omega^{*}}\right\|_{\infty}\sim\varepsilon we have

(with both ε,δ\varepsilon,\delta possibly dependent on kk but not on nn). Letting

which is a contradiction. Hence, (81) follows.

with ε,δ\varepsilon,\delta independent of nn. Let

4 Proof of Proposition 9.3

for all i,j,h∈[k]i,j,h\in\left[{k}\right]. Therefore, the assertion follows from Lemma 4.3 (the chain rule) and Lemma 9.12. □\Box

Let ε>0\varepsilon>0. We say that Ψ∈C2((0,1)2,R)\Psi\in C^{2}((0,1)^{2},\mathbf{R}) is ε\varepsilon-tame on Y⊂(0,1)2\mathcal{Y}\subset(0,1)^{2} if the following conditions hold:

For all y∈(0,1)y\in(0,1) we have Ψ(y,y)=0\Psi(y,y)=0.

On Y\mathcal{Y} we have ∣∑i=12∂2Ψ∂zi∂zj∣≤ε\left|{\sum_{i=1}^{2}\frac{\partial^{2}\Psi}{\partial z_{i}\partial z_{j}}}\right|\leq\varepsilon for any j=1,2j=1,2.

On Y\mathcal{Y} we have ∣∑i,j=12∂2Ψ∂zi∂zj∣≤ε2\left|{\sum_{i,j=1}^{2}\frac{\partial^{2}\Psi}{\partial z_{i}\partial z_{j}}}\right|\leq\varepsilon^{2}.

On Y\mathcal{Y} we have ∣∂2Ψ∂zi∂zj∣≤100|\frac{\partial^{2}\Psi}{\partial z_{i}\partial z_{j}}|\leq 100 for any i,j=1,2i,j=1,2.

Let f:(0,1)k→R2f:(0,1)^{k}\rightarrow\mathbf{R}^{2}, (z1,…,zk)↦(f1(z1,…,zk),f2(z1,…,zk))(z_{1},\ldots,z_{k})\mapsto(f_{1}(z_{1},\ldots,z_{k}),f_{2}(z_{1},\ldots,z_{k})) be a C2C^{2}-function. We say that ff is ε\varepsilon-benign on W\mathcal{W} if the following statements are true on W\mathcal{W}:

∣∂f1∂z1−∂f2∂z1∣<ε\left|{\frac{\partial f_{1}}{\partial z_{1}}-\frac{\partial f_{2}}{\partial z_{1}}}\right|<\varepsilon.

∣∂fi∂zj∣<ε\left|{\frac{\partial f_{i}}{\partial z_{j}}}\right|<\varepsilon for any 1<j≤k1<j\leq k and i=1,2i=1,2 and ∣∂fi∂z1∣≤100\left|{\frac{\partial f_{i}}{\partial z_{1}}}\right|\leq 100.

∣∂2fi∂zh∂zj∣<ε\left|{\frac{\partial^{2}f_{i}}{\partial z_{h}\partial z_{j}}}\right|<\varepsilon for any ii and (h,j)≠(1,1)(h,j)\neq(1,1).

∣∂2f1∂z12−∂2f2∂z12∣<ε\left|{\frac{\partial^{2}f_{1}}{\partial z_{1}^{2}}-\frac{\partial^{2}f_{2}}{\partial z_{1}^{2}}}\right|<\varepsilon and ∣∂2f1∂z12∣≤100|\frac{\partial^{2}f_{1}}{\partial z_{1}^{2}}|\leq 100.

There is an absolute constant C>0C>0 such that the following is true. Assume that ff is ε\varepsilon-benign on W\mathcal{W} and that Ψ\Psi is ε\varepsilon-tame on f(W)f(\mathcal{W}). Then on W\mathcal{W} we have

Proof. By Lemma 4.3 (the chain rule), we have

Since by T4 and Taylor’s formula we have ∂Ψ∂yh=Ok(ε)\frac{\partial\Psi}{\partial y_{h}}=O_{k}(\varepsilon), B3 implies that for (i,j)≠(1,1)(i,j)\neq(1,1)

Furthermore, as ∂Ψ∂yh=Ok(ε)\frac{\partial\Psi}{\partial y_{h}}=O_{k}(\varepsilon), B4 yields

the last step follows from T2 and Taylor’s formula.

To deal with the second sum, we consider four cases.

By B2 we have ∂fa∂zi∂fb∂zj≤Ok(ε2)\frac{\partial f_{a}}{\partial z_{i}}\frac{\partial f_{b}}{\partial z_{j}}\leq O_{k}(\varepsilon^{2}), and thus

Hence, in all cases we obtain a bound of Ok(ε2)O_{k}(\varepsilon^{2}). □\Box

Proof. It is straightforward to work out the differentials of ψ\psi: we have

Differentiating once more with respect to y1y_{1}, we get

Therefore, at y1=y2+εy_{1}=y_{2}+\varepsilon the second derivatives work out to be

Hence, ψ\psi is tame. Furthermore, differentiating (y1,y2)↦(1−y2)ψ(y1,y2)(y_{1},y_{2})\mapsto(1-y_{2})\psi(y_{1},y_{2}) yields

Hence, the fact that (1−y2)ψ(y1,y2)(1-y_{2})\psi(y_{1},y_{2}) is ε\varepsilon-tame follows from the fact that ψ\psi is. □\Box

With \mathchoice{\mbox{\boldmath\displaystyle q}}{\mbox{\boldmath\textstyle q}}{\mbox{\boldmath\scriptstyle q}}{\mbox{\boldmath\scriptscriptstyle q}}=\mathchoice{\mbox{\boldmath\displaystyle q}}{\mbox{\boldmath\textstyle q}}{\mbox{\boldmath\scriptstyle q}}{\mbox{\boldmath\scriptscriptstyle q}}(\omega) the functions

Finally, Proposition 9.3 follows directly from Lemmas 9.13, 9.14, 9.15 and 9.16.

5 Proof of Lemma 9.12

Proceeding to the second derivative, we highlight the following (folklore) fact.

Proof. This is a simple consequence of Cramer’s rule. Indeed, let Aij′A_{ij}^{\prime} be the matrix obtained from AA by omitting row ii and column jj. Then

Thus, we need to differentiate det⁡Ats′\det A_{ts}^{\prime} and det⁡A\det A. For any i≠ji\neq j we have

Similarly, for i≠ji\neq j and s≠ts\neq t we have

Thus, the assertion follows from the quotient rule. □\Box

for any h,i,j∈[k]h,i,j\in\left[{k}\right]. Thus,

6 Completing the proof of Proposition 8.1

Therefore, the third assertion follows from Remark 6.14.

Enumeration of Assignments with pp-Marginals

In this section we will prove Lemma 6.2 and Proposition 8.5. Before we present the actual details we will introduce an appropriate framework, which will enable us to perform the enumeration of assignments with pp-marginals, and pairs of such assignments with a given overlap.

In Section 5 we said that an assignment σ∈{0,1}V\sigma\in\left\{{0,1}\right\}^{V} has p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}-marginals if for any type t∈Tt\in\mathcal{T} we have

In words, the fraction of literal occurrences of type tt that are true under σ\sigma equals p(t)p(t) up to an error of O(1/n)O(1/n). However, due to technical reasons and because it simplifies some of our calculations significantly, we will actually work with a slightly refined definition. Let us say that a signature (s,d+,d−)(s,d^{+},d^{-}) is good, if d+,d−<3kr/4d^{+},d^{-}<3kr/4 and 0<(d+−d)2≤100k2kln⁡k0<(d^{+}-d)^{2}\leq 100k2^{k}\ln k. Instead of requiring that the fraction of literal occurrences of type tt equals p(t)p(t), we require that this is true for every good signature. That is, we say that an assignment σ∈{0,1}V\sigma\in\left\{{0,1}\right\}^{V} has p_{\mathchoice{\mbox{\boldmath\displaystyle d}}{\mbox{\boldmath\textstyle d}}{\mbox{\boldmath\scriptstyle d}}{\mbox{\boldmath\scriptscriptstyle d}}}-marginals if for any good s∈Ts\in T

and moreover, that fraction of literal occurrences of all other variables is 1/21/2, i.e.,

We are going to prove Lemma 6.2 and Proposition 8.5 with this modified definition. It is easily checked that this modification does not affect any of the arguments in the previous sections.

Let t1/2t_{1/2} be the type such that p(t1/2)=1/2p(t_{1/2})=1/2. Since Lt1/2=L¬t1/2L_{t_{1/2}}=L_{\neg t_{1/2}} it follows that in this special case

With the above notation, an assignment σ\sigma has pp-marginals if and only if

The next proposition is the first step towards the estimation of the total number of assignments with pp-marginals, c.f. Lemma 6.2. We denote by H(x)=−xln⁡x−(1−x)ln⁡(1−x)H(x)=-x\ln x-(1-x)\ln(1-x) the entropy of xx, and with [zn]f(z)[z^{n}]f(z) the nn-th coefficient in the Taylor series expansion of an analytic function ff around 0.

W.h.p. d\bf d chosen from D\bf D has the following property. There is a constant C>0C>0 such that if we denote by S\cal S the set of signatures s∈Ts\in T with the property p(s)>1/2p(s)>1/2, then

Proof. First of all, note that if for an assignment σ\sigma and a signature s∈Ts\in T with p(s)>1/2p(s)>1/2 we have ws(σ)=π(s)kmw_{s}(\sigma)=\pi(s)km, then the fraction of variables in VsV_{s} that are set to true is p(s)p(s). Thus, the fraction of variables set to false is 1−p(s)=p(¬s)1-p(s)=p(\neg s), and we infer that

Consequently, for any such ss the number of partial assignments σs:Vs→{0,1}\sigma_{s}:V_{s}\to\{0,1\}, with the property that the fraction of satisfied variables is p(s)p(s) is

Since w.h.p. d\bf d is such that ∣Vs∣=(1+o(1))αsn|V_{s}|=(1+o(1))\alpha_{s}n for some αs=αs(k)\alpha_{s}=\alpha_{s}(k), this provides the exponential terms in (93).

It remains to bound the number of partial assignments σ′:Vt1/2→{0,1}\sigma^{\prime}:V_{t_{1/2}}\to\{0,1\} such that wt1/2=12π(t1/2)kmw_{t_{1/2}}=\frac{1}{2}\pi(t_{1/2})km. Define the generating function

By definition, the sought quantity is [zπ(t1/2)km/2]F(z)[z^{\pi(t_{1/2})km/2}]F(z). Moreover, the definition of F(z)F(z) and (92) imply that

Lemma 6.2 follows immediately from the next statement, which is shown in Section 10.1.

W.h.p. d\bf d chosen from D\bf D has the following property. There is a constant C=C(k)>0C=C(k)>0 such that is we write N=∣Vt1/2∣N=|V_{t_{1/2}}|, then

We proceed with the proof of Proposition 8.5, i.e., we want to enumerate pairs of assignments with pp-marginals that have a specific overlap. Let s∈Ts\in T be a signature. For any σ,τ∈{0,1}n\sigma,\tau\in\{0,1\}^{n} denote the by the ss-overlap os(σ,τ)o_{s}(\sigma,\tau) the number of literal occurrences that are satisfied in both σ\sigma and τ\tau, where we consider only literals of signature ss, i.e.,

Similarly, for any type t∈Tt\in\cal T we denote by ot(σ,τ)o_{t}(\sigma,\tau) the number of satisfied literal occurrences in both σ\sigma and τ\tau, where only literals of type tt are considered. Note that ot(σ,τ)=O(σ,τ)tπ(t)kmo_{t}(\sigma,\tau)={\cal O}(\sigma,\tau)_{t}\pi(t)km, where O\cal O is defined in Section 7.1. For the special case t=t1/2t=t_{1/2} it follows

Let us begin with a simple observation. Let s∈ts\in t such that p(s)>1/2p(s)>1/2, and let σ,τ\sigma,\tau be two assignments with pp-marginals. Note that if ws(σ,τ)=(1+δ)p(s)2π(s)kmw_{s}(\sigma,\tau)=(1+\delta)p(s)^{2}\pi(s)km, for some δ≥−1\delta\geq-1, then the fraction of variables in VsV_{s} that are set to true in σ\sigma and τ\tau is (1+δ)p(s)2(1+\delta)p(s)^{2}. Consequently, the number of variables that are set to false in both assignments is (1−p(s))∣Vs∣−(p(s)∣Vs∣−(1+δ)p(s)2∣Vs∣)(1-p(s))|V_{s}|-(p(s)|V_{s}|-(1+\delta)p(s)^{2}|V_{s}|), and therefore

In words, the overlap in ss determines the overlap in ¬s\neg s. However, note that the s′s^{\prime}-overlap, for any s′≠s,¬ss^{\prime}\neq s,\neg s, is not affected by the quantities ws(σ,τ)w_{s}(\sigma,\tau) and w¬s(σ,τ)w_{\neg s}(\sigma,\tau).

Let t∈Tt\in\cal T be a type. With the previous observation at hand we are able to estimate the number of pairs of pp-satisfying assignments with a given tt- and ¬t\neg t-overlap. The proof can be found in Section 10.2.

There is a c>0c>0 such that the following is true. Let ε,ε′>0\varepsilon,\varepsilon^{\prime}>0. Let t∈Tt\in\cal T be a type such that p(t)≠1/2p(t)\neq 1/2. Denote by Ht,¬t2(ε,ε′){\cal H}^{2}_{t,\neg t}(\varepsilon,\varepsilon^{\prime}) the set of pairs σ\sigma, τ\tau of assignments with pp-marginals, such that

What remains is to enumerate pairs of pp-satisfying assignments with a given t1/2t_{1/2}-overlap. The next proposition provides this number as the coefficient of an appropriately defined generating function.

Let ε∈(−1/4,1/4)\varepsilon\in(-1/4,1/4). Let H1/22(ε){\cal H}^{2}_{1/2}(\varepsilon) denote the set of pairs σ′,τ′\sigma^{\prime},\tau^{\prime} of assignments to the variables in Vt1/2V_{t_{1/2}} such that

Then H1/22(ε)=[(xy)π(t1/2)km/2 u(1/4+ε)π(t1/2)km]F(x,y,u){\cal H}^{2}_{1/2}(\varepsilon)=[(xy)^{\pi(t_{1/2})km/2}\,u^{(1/4+\varepsilon)\pi(t_{1/2})km}]F(x,y,u), where

Proof. Assign to a pair of assignments σ′,τ′\sigma^{\prime},\tau^{\prime} to the variables in Vt1/2V_{t_{1/2}} the weight xwt1/2(σ′) ywt1/2(τ′) uot1/2(σ′,τ′)x^{w_{t_{1/2}}(\sigma^{\prime})}\,y^{w_{t_{1/2}}(\tau^{\prime})}\,u^{o_{t_{1/2}}(\sigma^{\prime},\tau^{\prime})}. Then, by using (92) and (94)

Summing this expression up yields the claimed statement. □\Box

The next statement provides the asymptotic value of the sought coefficients of F(x,y,u)F(x,y,u) from the previous proposition. The proof can be found in Section 10.3.

W.h.p. d\bf d chosen from D\bf D has the following property. There is a constant C=C(k,ε)>0C=C(k,\varepsilon)>0 such that if we write N=∣Vt1/2∣N=|V_{t_{1/2}}| and M=π(t1/2)kmM=\pi(t_{1/2})km, then

and ρ\rho is the solution to the equation

In order to complete the proof of Proposition 8.5 we will estimate the exponential term in the previous statement as a function of ε\varepsilon. Note that if ε=0\varepsilon=0, then clearly ρ=1\rho=1 and E=4NE=4^{N}. Let ∣ε∣<1/100|\varepsilon|<1/100. We begin with providing bounds for the value of ρ\rho from Equation (96). Let fg(ρ)=g/(2+2ρg)f_{g}(\rho)=g/(2+2\rho^{g}), where g≥3g\geq 3. Then fg(1)=g/4,fg′(1)=−g2/8f_{g}(1)=g/4,f_{g}^{\prime}(1)=-g^{2}/8 and

Note that if 0≤ρ≤10\leq\rho\leq 1, then, with room to spare, ∣fg′′(ρ)∣≤g3|f_{g}^{\prime\prime}(\rho)|\leq g^{3}. Moreover, if ρ>1\rho>1, then we may estimate fg′′f_{g}^{\prime\prime} as follows:

Let us write ρ=1+δ\rho=1+\delta. Taylor’s theorem then implies that ∣fg(ρ)−(g/4−g2δ/8)∣≤g2δ2|f_{g}(\rho)-(g/4-g^{2}\delta/8)|\leq g^{2}\delta^{2}. By writing gv=dv+d¬vg_{v}=d_{v}+d_{\neg v} and recalling that M=∑v∈Vt1/2gvM=\sum_{v\in V_{t_{1/2}}}g_{v} we infer from (96)

In view of these inequalities we might expect that whenever ε\varepsilon is not too large, then δ≈−ε8M/S2\delta\approx-\varepsilon{8M}/{S_{2}}. This can be made precise as follows. By solving the quadratic equations explicitly we infer that δ\delta satisfies

Note that d\bf d is such that w.h.p. S2=Θ(krM)S_{2}=\Theta(krM) and S3=Θ((kr)2M)S_{3}=\Theta((kr)^{2}M). Thus, for sufficiently large kk

The square-root with the minus sign can be estimated analogously. We infer that

With the approximate value of ρ\rho at hand we can proceed with estimating the exponential term in (95). First of, we rearrange terms to obtain

Regarding the last term involving the product in (98), we bound it by the following probabilistic considerations. Note that

Let (Sv)v∈Vt1/2(S_{v})_{v\in V_{t_{1/2}}} be a family of independent random variables, which are uniformly distributed in {−1,+1}\{-1,+1\}. Then the last expression in the previous display is equal to the expected value of ρ{−1/2∑vsvgv}\rho\left\{-1/2\sum_{v}s_{v}g_{v}\right\}. We obtain

Note that since either ρt/2≥1\rho^{t/2}\geq 1 or ρ−t/2≥1\rho^{-t/2}\geq 1 we may assume without loss of generality that ρ≥1\rho\geq 1. The advantage of the above formulation is that we can estimate rather easily the probability for a large deviation of the sum S=∑vSvgvS=\sum_{v}S_{v}g_{v}. Indeed, if we change the value of any SvS_{v} to obtain a new sum S′S^{\prime}, then ∣S−S′∣=2gv|S-S^{\prime}|=2g_{v}. By applying Azuma-Hoeffding we obtain

Thus, by using (97) and noting that ε≤0\varepsilon\leq 0 due to our assumption ρ≥1\rho\geq 1 we obtain the bound

Since the exponent is convex in tt, it can easily be seen that it is maximized at t=4M∣ε∣t=4M|\varepsilon|, where its value equals

Thus, μ=O(N)e8ε2M2S2\mu=O(\sqrt{N})e^{8\varepsilon^{2}\frac{M^{2}}{S_{2}}}, and by combining (98) and (99) we infer that E≤Ne−8ε2M2S2E\leq\sqrt{N}e^{-8\varepsilon^{2}\frac{M^{2}}{S_{2}}}. But since S2=Θ(krM)S_{2}=\Theta(krM) and M=Θ(krN)M=\Theta(krN), this is at most Ne−cε2N\sqrt{N}e^{-c\varepsilon^{2}N}, for some c>0c>0.

Proposition 8.5 then follows immediately from Propositions 10.3-10.5, and the (aforementioned) observation that the tt- and t′t^{\prime}-overlap of \mathchoice{\mbox{\boldmath\displaystyle\sigma}}{\mbox{\boldmath\textstyle\sigma}}{\mbox{\boldmath\scriptstyle\sigma}}{\mbox{\boldmath\scriptscriptstyle\sigma}},\mathchoice{\mbox{\boldmath\displaystyle\tau}}{\mbox{\boldmath\textstyle\tau}}{\mbox{\boldmath\scriptstyle\tau}}{\mbox{\boldmath\scriptscriptstyle\tau}} are independent for t≠t′,¬tt\neq t^{\prime},\neg t.

Set M=π(t1/2)kmM=\pi(t_{1/2})km. By the virtue of Cauchy’s integral formula we obtain

Since FF is analytic in C\mathbf{C}, CC can be any curve enclosing the origin. To estimate the integral we will use the saddle point method, which is commonly used to determine the asymptotic behavior of integrals that involve a large parameter, and are simultaneously subject to huge variations. For an excellent overview and numerous applications we refer the reader to .

The main idea is to choose CC such that the integrand ’peaks’ at a unique point, so that the main contribution to the integral comes from a small neighborhood of this maximum. We choose CC to be the unit circle centered at the origin, i.e., C={eiθ:−π<θ<π}C=\{e^{i\theta}:-\pi<\theta<\pi\}. Moreover, let θ0=θ0(n)=N−2/5\theta_{0}=\theta_{0}(n)=N^{-2/5}, and write C0={eiθ:∣θ∣≤θ0(n)}C_{0}=\{e^{i\theta}:|\theta|\leq\theta_{0}(n)\} for the restriction of CC to the segment with ∣θ∣≤θ0(n)|\theta|\leq\theta_{0}(n). Then we may write I=I0+I1I=I_{0}+I_{1}, where

By changing variables, the first integral becomes

Moreover, by using the trivial bound for complex integrals and the fact ∣z∣=1|z|=1 on CC we obtain

Our subsequent proof strategy is as follows. We will first compute the asymptotic value of the integral over the ’central region’; in particular, we show that

for an appropriate c>0c>0. Then, by using (101) we show that I1=o(I0)I_{1}=o(I_{0}). The two statements combined yield then immediately the conclusion of the proposition.

We proceed with showing (102). Recall that ∣θ∣≤θ0=N−2/5|\theta|\leq\theta_{0}=N^{-2/5}, and note that for any d,d′d,d^{\prime}, by applying Taylor’s Theorem

Observe that d\bf d is w.h.p. such that Sj=(1+o(1))cjNS_{j}=(1+o(1))c_{j}N for some cj=cj(k)>0c_{j}=c_{j}(k)>0, where 2≤j≤92\leq j\leq 9. Using (100) we infer that the integrand satisfies

This proves (102). To complete the proof we will show that sup⁡z∈C∖C0∣F(z)∣\sup_{z\in C\setminus C_{0}}\left|F(z)\right| is asymptotically negligible compared to I0I_{0}. First, for any v∈Vt1/2v\in V_{t_{1/2}}

Let us collect some basic properties of fvf_{v}. Note that if dv=d¬vd_{v}=d_{\neg v}, then fv(θ)=2f_{v}(\theta)=2 for any −π<θ<π-\pi<\theta<\pi. Otherwise, ff is maximized for any

For a pair (d+,d−)∈N2(d_{+},d_{-})\in\mathbf{N}^{2} let Vd+,d−⊆Vt1/2V_{d_{+},d_{-}}\subseteq V_{t_{1/2}} denote the set of variables vv such that dv=d+d_{v}=d_{+} and d¬v=d−d_{\neg v}=d_{-}, and write Nd+,d−=∣Vd+,d−∣N_{d_{+},d_{-}}=|V_{d_{+},d_{-}}|. Then,

Note that ∑s=(d+,d−)Ns=N\sum_{s=(d_{+},d_{-})}N_{s}=N. Thus, ∣F(eiθ)∣≤2N|F(e^{i\theta})|\leq 2^{N} for all θ\theta. However, this bound is achieved only if all factors are maximized simultaneously. We will argue in the sequel that if ∣θ∣∈(θ0,π)|\theta|\in(\theta_{0},\pi), then a linear (in NN) fraction of the factors is ≤2−O(N−4/5)\leq 2-O(N^{-4/5}). It follows for some α>0\alpha>0 that

To see the claim, consider the specific pair (d+′,d−′)=(kr,kr−1)(d^{\prime}_{+},d^{\prime}_{-})=(kr,kr-1), and note that if kk is sufficiently large, then kr−1>kr/2+10k2kln⁡kkr-1>kr/2+10\sqrt{k2^{k}\ln k}. So, indeed Vd+,d−⊆Vt1/2V_{d_{+},d_{-}}\subseteq V_{t_{1/2}}. Furthermore, d\bf d is such that w.h.p. there is a constant α=α(k)>0\alpha=\alpha(k)>0 such that Nd+,d−≥αNN_{d_{+},d_{-}}\geq\alpha N. It follows that for all variables v∈Vd+′,d−′v\in V_{d^{\prime}_{+},d^{\prime}_{-}}

It can easily be verified that fvf_{v} is monotone increasing for −π<θ<0-\pi<\theta<0 and decreasing for 0<θ<π0<\theta<\pi. Thus, for any ∣θ∣∈(θ0,π)|\theta|\in(\theta_{0},\pi) we have fv(θ)≤max⁡{fv(θ0),fv(−θ0)}f_{v}(\theta)\leq\max\{f_{v}(\theta_{0}),f_{v}(-\theta_{0})\}. By using the Taylor series expansion of the cosine and the square root we obtain that

We conclude that fv(θ)≤2−O(n−4/5)f_{v}(\theta)\leq 2-O(n^{-4/5}) for at least αN\alpha N variables vv, and the proof is completed.

2 Proof of Proposition 10.3

We will exploit a concentration inequality due to McDiarmid . We present it here in a simplified form that is appropriate for our purpose. Given a finite non-empty set BB, we denote by Sym(B)Sym(B) the set of all ∣B∣!|B|! permutations of the elements of BB. Let B1,…,BNB_{1},\dots,B_{N} be a family of finite non-empty sets, and denote by Ω=Sym(B1)×⋯×Sym(BN)\Omega=Sym(B_{1})\times\dots\times Sym(B_{N}). Moreover, let \mathchoice{\mbox{\boldmath\displaystyle\pi}}{\mbox{\boldmath\textstyle\pi}}{\mbox{\boldmath\scriptstyle\pi}}{\mbox{\boldmath\scriptscriptstyle\pi}}=(\pi_{1},\dots,\pi_{N}) be a family of independent random permutations, where πi\pi_{i} is drawn uniformly from Sym(Bi)Sym(B_{i}).

Let cc and rr be positive constants. Suppose that h:Ω→R+h:\Omega\to\mathbf{R}_{+} is such that for any π∈Ω\pi\in\Omega the following conditions are satisfied.

If π′\pi^{\prime} can be obtained from π\pi by swapping two elements, then ∣h(π)−h(π′)∣≤c|h(\pi)-h(\pi^{\prime})|\leq c.

If h(π)≥sh(\pi)\geq s, then there is a set of at most rsrs coordinates such that h(π′)≥sh(\pi^{\prime})\geq s for any π′∈Ω\pi^{\prime}\in\Omega that agrees with π\pi on these coordinates.

Let Z=h(\mathchoice{\mbox{\boldmath\displaystyle\pi}}{\mbox{\boldmath\textstyle\pi}}{\mbox{\boldmath\scriptstyle\pi}}{\mbox{\boldmath\scriptscriptstyle\pi}}) and let mm be the median of ZZ. Then, for any t>0t>0

Let us proceed with the proof of Proposition 10.3. We will assume without loss of generality that tt is such that p(t)>1/2p(t)>1/2. We will abbreviate p=p(t)p=p(t), q=p(¬t)q=p(\neg t). Let σ\sigma be an arbitrary assignment with pp-marginals. Moreover, denote by τ\textstyle\tau an assignment that is obtained by selecting for any signature s∈ts\in t uniformly at random p∣Vs∣p|V_{s}| variables from VsV_{s} and setting them to true, and setting all other variables in V∖VtV\setminus V_{t} arbitrarily so that τ\textstyle\tau has pp-marginals. Equivalently, we may generate τ\textstyle\tau by permuting the variables in VsV_{s} randomly, and setting the first p∣Vs∣p|V_{s}| variables in that permutation to true, for all s∈ts\in t. With this notation we obtain

The latter probability can be estimated with Theorem 10.6. Indeed, note that

if τ,τ′\tau,\tau^{\prime} have pp-marginals and can be obtained by swapping the truth assignment of two variables, then

if wt(σ,τ)≥sw_{t}(\sigma,\tau)\geq s, then there is a set SS of ≤s/min⁡v∈Vtdv≤2s/kr\leq s/\min_{v\in V_{t}}d_{v}\leq 2s/kr variables that are set to true, and any τ′\tau^{\prime} with pp-marginals that sets all variables is SS to true satisfies wt(σ,τ′)≥sw_{t}(\sigma,\tau^{\prime})\geq s.

Exactly the same argument, where we interchange the roles of tt and ¬t\neg t, shows that also

3 Proof of Proposition 10.5

Set M=π(t1/2)kmM=\pi(t_{1/2})km. By applying Cauchy’s integral formula we obtain

The function FF is analytic in C3\mathbf{C}^{3}, implying that C1,C2,CoC_{1},C_{2},C_{o} can be any curves enclosing the origin. We choose

where ρ\rho is the solution to the Equation (96). Some remarks are in place here. The choice of the integration paths may seem arbitrary at this point. Note, however, that FF is symmetric with respect to xx and yy, and thus it is natural to assume similar integration curves for them. Moreover, the choice of ρ\rho is guided by the general principles of the saddle-point method and is such that the integrand has a unique maximum at (θ,φ,ψ)=(0,0,0)(\theta,\varphi,\psi)=(0,0,0). Indeed, as we will show subsequently, the integrand is around (0,0,0)(0,0,0) of elliptic type; this allows us to reduce the estimation of the main terms to the evaluation of a 3-dimensional Gaussian integral.

Denote by C\cal C the restriction of the circles C1,C2,CoC_{1},C_{2},C_{o} to a small region around the origin, i.e.,

and I1I_{1} is the integral over (C1×C2×Co)∖C(C_{1}\times C_{2}\times C_{o})\setminus\cal C. By changing variables we obtain

Regarding I1I_{1}, we will use the trivial bound

We begin with estimating I0I_{0} by providing an appropriate asymptotic expansion of it for points around the origin. First of all, note that for any v∈Vt1/2v\in V_{t_{1/2}} we have hv(0,0,0)=2+2ρdv+d¬vh_{v}(0,0,0)=2+2\rho^{d_{v}+d_{\neg v}} and thus

The second derivatives at (0,0,0)(0,0,0) are given by

Furthermore, the mixed second derivatives are

We will also need crude bounds for the third-order derivatives in order to establish an accurate approximation for HH around the origin. Note that hvh_{v} linearly exponential in θ,φ,ψ\theta,\varphi,\psi and dv,d¬vd_{v},d_{\neg v}. Thus, every time we take a derivative with respect to some variable, the norm of each single term in the expression of hvh_{v} can increase by at most mv=max⁡{dv,d¬v}m_{v}=\max\{d_{v},d_{\neg v}\}. Thus, uniformly for (θ,φ,ψ)∈[−N2/5,N2/5](\theta,\varphi,\psi)\in[-N^{2/5},N^{2/5}] we have that

By using the uniform estimate 1+x=ex−x2/2+Θ(x3)1+x=e^{x-x^{2}/2+\Theta(x^{3})}, where we set 1+x=hv(θ,φ,ψ)/hv(0,0,0)1+x=h_{v}(\theta,\varphi,\psi)/h_{v}(0,0,0) we infer that

Finally, since (θ,φ,ψ)∈[−N2/5,N2/5](\theta,\varphi,\psi)\in[-N^{2/5},N^{2/5}] the error term is of order at most (dv+d¬v)3N−6/5(d_{v}+d_{\neg v})^{3}N^{-6/5}. In order to obtain an approximation for HH we form the product over all v∈Vt1/2v\in V_{t_{1/2}}. Observe that the (linear in the variables) exponential factor e−i(θ+φ)M/2−iψ(1/4+ε)Me^{-i(\theta+\varphi)M/2-i\psi(1/4+\varepsilon)M} cancels exactly with the first order terms in (105). By abbreviating

we obtain uniformly for any (θ,φ,ψ)∈[−N−2/5,N−2/5]3(\theta,\varphi,\psi)\in[-N^{-2/5},N^{-2/5}]^{3}

Observe that d\bf d is such that w.h.p. all quantities S.,.S_{.,.} and S3S_{3} are linear in NN. Thus, we are left with computing

In order to compute this integral we rescale each variable with N−1/2N^{-1/2}. By writing s.,.s_{.,.} for S.,./NS_{.,.}/N we obtain

A termwise comparison and elementary algebraic manipulations yield that

Thus, the squares can be completed and the integral in the above expression equals a constant depending on the family s.,.s_{.,.}; this shows that asymptotically I1I_{1} is proportional to N−3/2⋅EN^{-3/2}\cdot E.

In order to complete the proof we will use (104) to show that I1I_{1} is asymptotically negligible compared to I0I_{0}. Recall the definition of HH from (103). It follows that the absolute value of HH is given by

Let us abbreviate Dv=dv−d¬vD_{v}=d_{v}-d_{\neg v}. A lengthy calculation, which can be performed easily with the help of MAPLE, yields that

Note that we can get an upper bound for fvf_{v} if we replace all terms involving a cosine by one; this implies that ∣H∣≤ρ−(1−4ε)M/2∏v(2+2ρdv+d¬v)=E|H|\leq\rho^{-(1-4\varepsilon)M/2}\prod_{v}(2+2\rho^{d_{v}+d_{\neg v}})=E. Moreover, the bound is achieved only if all factors are maximized simultaneously, and this happens for example when we choose (θ,φ,ψ)=(0,0,0)(\theta,\varphi,\psi)=(0,0,0). We will argue in the sequel that if (θ,φ,ψ)∈(C1×C2×Co)∖C(\theta,\varphi,\psi)\in(C_{1}\times C_{2}\times C_{o})\setminus\cal C, i.e., at least one of the variables θ,φ,ψ\theta,\varphi,\psi is assigned a value not lying in [−N−2/5,N−2/5][-N^{-2/5},N^{-2/5}], then there is a subset of variables V′⊂Vt1/2V^{\prime}\subset V_{t_{1/2}} such that ∣V′∣≥αN|V^{\prime}|\geq\alpha N for some α>0\alpha>0 and for all v∈V′v\in V^{\prime} it holds fv(0,0,0)≤fv(0,0,0)−O(N−4/5)f_{v}(0,0,0)\leq f_{v}(0,0,0)-O(N^{-4/5}). Indeed, if this is true, then

Since ρ\rho is bounded and d\bf d is such that w.h.p. dv+d¬v=o(log⁡n)d_{v}+d_{\neg v}=o(\log n), it follows that ∣H∣|H| smaller that EE by an exponential factor, which shows with (104) that I1=o(I0)I_{1}=o(I_{0}).

To see that a set V′V^{\prime} with the desired properties exists, let us assume that at least one of θ,φ,ψ\theta,\varphi,\psi is in absolute value at least N−2/5N^{-2/5}. For a pair (d+,d−)∈N2(d_{+},d_{-})\in\mathbf{N}^{2} let Vd+,d−⊆Vt1/2V_{d_{+},d_{-}}\subseteq V_{t_{1/2}} denote the set of variables vv such that dv=d+d_{v}=d_{+} and d¬v=d−d_{\neg v}=d_{-}, and write Nd+,d−=∣Vd+,d−∣N_{d_{+},d_{-}}=|V_{d_{+},d_{-}}|. Consider the specific pair (d+,d−)=(kr,kr−1)(d_{+},d_{-})=(kr,kr-1), and note that for all such variables we have Dv=1D_{v}=1. Furthermore, d\bf d is such that w.h.p. there is a constant β=β(k)>0\beta=\beta(k)>0 such that Nd+,d−≥βNN_{d_{+},d_{-}}\geq\beta N. Then we may assume that

as otherwise there is nothing to show. This impliesthat the arguments of all cosines appearing in the expression of fvf_{v} are close to multiples of 2π2\pi, and in particular,

this follows directly from the series expansion of the cosine around integer multiples of 2π2\pi, which lack a linear term. Next, consider the pair (d+′,d−′)=(kr,kr−2)(d^{\prime}_{+},d^{\prime}_{-})=(kr,kr-2); again d\bf d is such that w.h.p. there is a constant β′=β′(k)>0\beta^{\prime}=\beta^{\prime}(k)>0 such that Nd+′,d−′≥β′NN_{d^{\prime}_{+},d^{\prime}_{-}}\geq\beta^{\prime}N. Note that for these variables we have Dv=2D_{v}=2. Then, as previously, we may also assume that

But then, by the same argument as above, ∣2φ+d+′ψ∣=O(N−2/5) ( mod  2π)|2\varphi+d^{\prime}_{+}\psi|=O(N^{-2/5})~(\bmod~2\pi). Since d+=d+′d_{+}=d^{\prime}_{+} and, by assumption, ∣φ∣<π|\varphi|<\pi, by combining this with the third term in (106), we infer that ∣φ∣=O(N−2/5)|\varphi|=O(N^{-2/5}). In turn, together with the second term in (106), this implies that also ∣θ∣=O(N−2/5)|\theta|=O(N^{-2/5}). Finally, the fact ∣θ+φ+ψ∣=O(N−2/5) ( mod  2π)|\theta+\varphi+\psi|=O(N^{-2/5})~(\bmod~2\pi) from (106) then also implies that ∣δ∣=O(N−2/5)|\delta|=O(N^{-2/5}). Everything together yields that (θ,φ,ψ)∈C(\theta,\varphi,\psi)\in\cal C, a contradiction.

Proof of Corollary 2.2

As a direct consequence of our second moment argument, the Paley-Zygmund inequality, and a concentration result on the number of satisfying assignments from we obtain the following.

There is a number δ=δ(k)>0\delta=\delta(k)>0 such that

for all x∈Vx\in V. Furthermore, if we let E{\cal E} be the event that ∑l∈Lel=km\sum_{l\in L}e_{l}=km, then \mathchoice{\mbox{\boldmath\displaystyle e}}{\mbox{\boldmath\textstyle e}}{\mbox{\boldmath\scriptstyle e}}{\mbox{\boldmath\scriptscriptstyle e}}=(e_{l})_{l\in L} given E{\cal E} has the same distribution as d\textstyle d. Moreover,

Finally, the assertion follows from (108) and (109). □\Box

Proof of Lemma 2.3

The expected majority weight in Φ\textstyle\Phi is easily computed. In Φ\textstyle\Phi, for each xx the numbers dxd_{x}, d¬xd_{\neg x} of positive/negative occurrences are asymptotically independently Poisson with mean kr/2kr/2. Therefore, for any d=Θ(kr)d=\Theta(kr) we obtain

By comparison, given that, say, the all-true assignment is satisfying, the number dxd_{x} of positive occurrences has distribution Po((1+1/(2k−1))kr/2){\rm Po}((1+1/(2^{k}-1))kr/2), while d¬xd_{\neg x} has distribution Po((1−1/(2k−1))kr/2){\rm Po}((1-1/(2^{k}-1))kr/2). The normal approximation to the Poisson distribution yields for d=Θ(kr)d=\Theta(kr),

for a certain constant c>0c>0. Consequently,

Both with and without conditioning on \mathchoice{\mbox{\boldmath\displaystyle 1}}{\mbox{\boldmath\textstyle 1}}{\mbox{\boldmath\scriptstyle 1}}{\mbox{\boldmath\scriptscriptstyle 1}}\in\mathcal{S}(\mathchoice{\mbox{\boldmath\displaystyle\Phi}}{\mbox{\boldmath\textstyle\Phi}}{\mbox{\boldmath\scriptstyle\Phi}}{\mbox{\boldmath\scriptscriptstyle\Phi}}), wmajw_{maj} enjoys the following Lipschitz property: changing one single clause can alter the value of wmajw_{maj} by at most k/(km)=1/(rn)k/(km)=1/(rn). Therefore, Azuma’s inequality yields

In effect, for a certain constant ζ>0\zeta>0 we have

Combining (112) and (113) with a simple counting argument yields Lemma 2 from the extended abstract.

Acknowledgment. The first author thanks Dimitris Achlioptas for helpful discussions on the second moment method. We also thank Charilaos Efthymiou for helpful comments that have led to an improved presentation.

References