Order-to-chaos transition in the hardness of random Boolean satisfiability problems

TL;DR

This paper investigates the emergence of transient chaos in the dynamics of an analog algorithm solving random Boolean satisfiability problems, revealing a phase transition in problem hardness related to constraint density.

Contribution

It introduces a dynamical systems approach to analyze the hardness transition in random $k$-SAT problems, identifying a second-order phase transition and metastable basins.

Findings

Hardness varies widely and follows a lognormal distribution.

A second-order phase transition occurs at a critical constraint density $oldsymbol{\alpha_{ ext{ extless} ext{ extgreater}}}$ in 3-SAT.

Transient chaos is characterized by an exponential-algebraic dependence of the escape rate on parameters.

Abstract

Transient chaos is an ubiquitous phenomenon characterizing the dynamics of phase space trajectories evolving towards a steady state attractor in physical systems as diverse as fluids, chemical reactions and condensed matter systems. Here we show that transient chaos also appears in the dynamics of certain efficient algorithms searching for solutions of constraint satisfaction problems that include scheduling, circuit design, routing, database problems or even Sudoku. In particular, we present a study of the emergence of hardness in Boolean satisfiability ($k$-SAT), a canonical class of constraint satisfaction problems, by using an analog deterministic algorithm based on a system of ordinary differential equations. Problem hardness is defined through the escape rate $\kappa$, an invariant measure of transient chaos of the dynamical system corresponding to the analog algorithm, and it expresses the rate at which the trajectory approaches a solution.We show that for a given density of constraints and fixed number of Boolean variables $N$, the hardness of formulas in random $k$-SAT ensembles has a wide variation, approximable by a lognormal distribution. We also show that when increasing the density of constraints $\alpha$, hardness appears through a second-order phase transition at $\alpha_{\chi}$ in the random 3-SAT ensemble where dynamical trajectories become transiently chaotic. A similar behavior is found in 4-SAT as well, however, such transition does not occur for 2-SAT. This behavior also implies a novel type of transient chaos in which the escape rate has an exponential-algebraic dependence on the critical parameter $\kappa \sim N^{B|\alpha - \alpha_{\chi}|^{1-\gamma}}$ with $0< \gamma < 1$. We demonstrate that the transition is generated by the appearance of metastable basins in the solution space as the density of constraints $\alpha$ is increased.