Researcher profile

David B. Wilson

David B. Wilson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
19works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2011arXiv

A sharp threshold for minimum bounded-depth and bounded-diameter spanning trees and Steiner trees in random networks

In the complete graph on n vertices, when each edge has a weight which is an exponential random variable, Frieze proved that the minimum spanning tree has weight tending to zeta(3)=1/1^3+1/2^3+1/3^3+... as n goes to infinity. We consider spanning trees constrained to have depth bounded by k from a specified root. We prove that if k > log_2 log n+omega(1), where omega(1) is any function going to infinity with n, then the minimum bounded-depth spanning tree still has weight tending to zeta(3) as n -> infinity, and that if k < log_2 log n, then the weight is doubly-exponentially large in log_2 log n - k. It is NP-hard to find the minimum bounded-depth spanning tree, but when k < log_2 log n - omega(1), a simple greedy algorithm is asymptotically optimal, and when k > log_2 log n+omega(1), an algorithm which makes small changes to the minimum (unbounded depth) spanning tree is asymptotically optimal. We prove similar results for minimum bounded-depth Steiner trees, where the tree must connect a specified set of m vertices, and may or may not include other vertices. In particular, when m = const * n, if k > log_2 log n+omega(1), the minimum bounded-depth Steiner tree on the complete graph has asymptotically the same weight as the minimum Steiner tree, and if 1 <= k <= log_2 log n-omega(1), the weight tends to (1-2^{-k}) sqrt{8m/n} [sqrt{2mn}/2^k]^{1/(2^k-1)} in both expectation and probability. The same results hold for minimum bounded-diameter Steiner trees when the diameter bound is 2k; when the diameter bound is increased from 2k to 2k+1, the minimum Steiner tree weight is reduced by a factor of 2^{1/(2^k-1)}.

preprint2011arXiv

Double-dimer pairings and skew Young diagrams

We study the number of tilings of skew Young diagrams by ribbon tiles shaped like Dyck paths, in which the tiles are &#34;vertically decreasing&#34;. We use these quantities to compute pairing probabilities in the double-dimer model: Given a planar bipartite graph $G$ with special vertices, called nodes, on the outer face, the double-dimer model is formed by the superposition of a uniformly random dimer configuration (perfect matching) of $G$ together with a random dimer configuration of the graph formed from $G$ by deleting the nodes. The double-dimer configuration consists of loops, doubled edges, and chains that start and end at the boundary nodes. We are interested in how the chains connect the nodes. An interesting special case is when the graph is $ε(\Z\times\N)$ and the nodes are at evenly spaced locations on the boundary $\R$ as the grid spacing $ε\to 0$.

preprint2011arXiv

XOR-Ising Loops and the Gaussian Free Field

We find by simulation that the interfaces in the exclusive-or (XOR) of two independent 2D Ising spin configurations at the critical temperature form an ensemble of loops that have the same distribution as the contour lines of the Gaussian free field, but with the heights of the contours spaced \sqrt{2} times as far apart as they are for the conformal loop ensemble CLE_4 or the double dimer model on the square lattice. For domains with boundary, various natural boundary conditions for the two Ising models correspond to certain boundary heights for the Gaussian free field.

preprint2010arXiv

Driving sandpiles to criticality and beyond

A popular theory of self-organized criticality relates driven dissipative systems to systems with conservation. This theory predicts that the stationary density of the abelian sandpile model equals the threshold density of the fixed-energy sandpile. We refute this prediction for a wide variety of underlying graphs, including the square grid. Driven dissipative sandpiles continue to evolve even after reaching criticality. This result casts doubt on the validity of using fixed-energy sandpiles to explore the critical behavior of the abelian sandpile model at stationarity.

preprint2010arXiv

Event-chain Monte Carlo algorithms for hard-sphere systems

In this paper we present the event-chain algorithms, which are fast Markov-chain Monte Carlo methods for hard spheres and related systems. In a single move of these rejection-free methods, an arbitrarily long chain of particles is displaced, and long-range coherent motion can be induced. Numerical simulations show that event-chain algorithms clearly outperform the conventional Metropolis method. Irreversible versions of the algorithms, which violate detailed balance, improve the speed of the method even further. We also compare our method with a recent implementations of the molecular-dynamics algorithm.

preprint2010arXiv

The approach to criticality in sandpiles

A popular theory of self-organized criticality relates the critical behavior of driven dissipative systems to that of systems with conservation. In particular, this theory predicts that the stationary density of the abelian sandpile model should be equal to the threshold density of the corresponding fixed-energy sandpile. This &#34;density conjecture&#34; has been proved for the underlying graph Z. We show (by simulation or by proof) that the density conjecture is false when the underlying graph is any of Z^2, the complete graph K_n, the Cayley tree, the ladder graph, the bracelet graph, or the flower graph. Driven dissipative sandpiles continue to evolve even after a constant fraction of the sand has been lost at the sink. These results cast doubt on the validity of using fixed-energy sandpiles to explore the critical behavior of the abelian sandpile model at stationarity.

preprint2010arXiv

Winding angle variance of Fortuin-Kasteleyn contours

The variance in the winding number of various random fractal curves, including the self-avoiding walk, the loop-erased random walk, contours of FK clusters, and stochastic Loewner evolution, have been studied by numerous researchers. Usually the focus has been on the winding at the endpoints. We measure the variance in winding number at typical points along the curve. More generally, we study the winding at points where k strands come together, and some adjacent strands may be conditioned not to hit each other. The measured values are consistent with an interesting conjecture.

preprint2009arXiv

Boundary Partitions in Trees and Dimers

Given a finite planar graph, a grove is a spanning forest in which every component tree contains one or more of a specified set of vertices (called nodes) on the outer face. For the uniform measure on groves, we compute the probabilities of the different possible node connections in a grove. These probabilities only depend on boundary measurements of the graph and not on the actual graph structure, i.e., the probabilities can be expressed as functions of the pairwise electrical resistances between the nodes, or equivalently, as functions of the Dirichlet-to-Neumann operator (or response matrix) on the nodes. These formulae can be likened to generalizations (for spanning forests) of Cardy&#39;s percolation crossing probabilities, and generalize Kirchhoff&#39;s formula for the electrical resistance. Remarkably, when appropriately normalized, the connection probabilities are in fact integer-coefficient polynomials in the matrix entries, where the coefficients have a natural algebraic interpretation and can be computed combinatorially. A similar phenomenon holds in the so-called double-dimer model: connection probabilities of boundary nodes are polynomial functions of certain boundary measurements, and as formal polynomials, they are specializations of the grove polynomials. Upon taking scaling limits, we show that the double-dimer connection probabilities coincide with those of the contour lines in the Gaussian free field with certain natural boundary conditions. These results have direct application to connection probabilities for multiple-strand SLE_2, SLE_8, and SLE_4.

preprint2008arXiv

Combinatorics of Tripartite Boundary Connections for Trees and Dimers

A grove is a spanning forest of a planar graph in which every component tree contains at least one of a special subset of vertices on the outer face called nodes. For the natural probability measure on groves, we compute various connection probabilities for the nodes in a random grove. In particular, for &#34;tripartite&#34; pairings of the nodes, the probability can be computed as a Pfaffian in the entries of the Dirichlet-to-Neumann matrix (discrete Hilbert transform) of the graph. These formulas generalize the determinant formulas given by Curtis, Ingerman, and Morrow, and by Fomin, for parallel pairings. These Pfaffian formulas are used to give exact expressions for reconstruction: reconstructing the conductances of a planar graph from boundary measurements. We prove similar theorems for the double-dimer model on bipartite planar graphs.

preprint2008arXiv

The electrical response matrix of a regular 2n-gon

Consider a unit-resistive plate in the shape of a regular polygon with 2n sides, in which even-numbered sides are wired to electrodes and odd-numbered sides are insulated. The response matrix, or Dirichlet-to-Neumann map, allows one to compute the currents flowing through the electrodes when they are held at specified voltages. We show that the entries of the response matrix of the regular 2n-gon are given by the differences of cotangents of evenly spaced angles, and we describe some connections with the limiting distributions of certain random spanning forests.

preprint2008arXiv

Tug-of-war and the infinity Laplacian

We prove that every bounded Lipschitz function F on a subset Y of a length space X admits a tautest extension to X, i.e., a unique Lipschitz extension u for which Lip_U u = Lip_{boundary of U} u for all open subsets U of X that do not intersect Y. This was previously known only for bounded domains R^n, in which case u is infinity harmonic, that is, a viscosity solution to Delta_infty u = 0. We also prove the first general uniqueness results for Delta_infty u = g on bounded subsets of R^n (when g is uniformly continuous and bounded away from zero), and analogous results for bounded length spaces. The proofs rely on a new game-theoretic description of u. Let u^epsilon(x) be the value of the following two-player zero-sum game, called tug-of-war: fix x_0=x \in X minus Y. At the kth turn, the players toss a coin and the winner chooses an x_k with d(x_k, x_{k-1})< epsilon. The game ends when x_k is in Y, and player one&#39;s payoff is F(x_k) - (epsilon^2/2) sum_{i=0}^{k-1} g(x_i) We show that the u^εconverge uniformly to u as epsilon tends to zero. Even for bounded domains in R^n, the game theoretic description of infinity-harmonic functions yields new intuition and estimates; for instance, we prove power law bounds for infinity-harmonic functions in the unit disk with boundary values supported in a delta-neighborhood of a Cantor set on the unit circle.

preprint2008arXiv

Two-player Knock &#39;em Down

We analyze the two-player game of Knock &#39;em Down, asymptotically as the number of tokens to be knocked down becomes large. Optimal play requires mixed strategies with deviations of order sqrt(n) from the naive law-of-large numbers allocation. Upon rescaling by sqrt(n) and sending n to infinity, we show that optimal play&#39;s random deviations always have bounded support and have marginal distributions that are absolutely continuous with respect to Lebesgue measure.

preprint2004arXiv

Balanced Boolean functions that can be evaluated so that every input bit is unlikely to be read

A Boolean function of n bits is balanced if it takes the value 1 with probability 1/2. We exhibit a balanced Boolean function with a randomized evaluation procedure (with probability 0 of making a mistake) so that on uniformly random inputs, no input bit is read with probability more than Theta(n^{-1/2} sqrt{log n}). We give a balanced monotone Boolean function for which the corresponding probability is Theta(n^{-1/3} log n). We then show that for any randomized algorithm for evaluating a balanced Boolean function, when the input bits are uniformly random, there is some input bit that is read with probability at least Theta(n^{-1/2}). For balanced monotone Boolean functions, there is some input bit that is read with probability at least Theta(n^{-1/3}).

preprint2002arXiv

On the critical exponents of random k-SAT

There has been much recent interest in the satisfiability of random Boolean formulas. A random k-SAT formula is the conjunction of m random clauses, each of which is the disjunction of k literals (a variable or its negation). It is known that when the number of variables n is large, there is a sharp transition from satisfiability to unsatisfiability; in the case of 2-SAT this happens when m/n --> 1, for 3-SAT the critical ratio is thought to be m/n ~ 4.2. The sharpness of this transition is characterized by a critical exponent, sometimes called ν=ν_k (the smaller the value of νthe sharper the transition). Experiments have suggested that ν_3 = 1.5+-0.1, ν_4 = 1.25+-0.05, ν_5=1.1+-0.05, ν_6 = 1.05+-0.05, and heuristics have suggested that ν_k --> 1 as k --> infinity. We give here a simple proof that each of these exponents is at least 2 (provided the exponent is well-defined). This result holds for each of the three standard ensembles of random k-SAT formulas: m clauses selected uniformly at random without replacement, m clauses selected uniformly at random with replacement, and each clause selected with probability p independent of the other clauses. We also obtain similar results for q-colorability and the appearance of a q-core in a random graph.

preprint2001arXiv

The Scaling Window of the 2-SAT Transition

We consider the random 2-satisfiability problem, in which each instance is a formula that is the conjunction of m clauses of the form (x or y), chosen uniformly at random from among all 2-clauses on n Boolean variables and their negations. As m and n tend to infinity in the ratio m/n --> alpha, the problem is known to have a phase transition at alpha_c = 1, below which the probability that the formula is satisfiable tends to one and above which it tends to zero. We determine the finite-size scaling about this transition, namely the scaling of the maximal window W(n,delta) = (alpha_-(n,delta),alpha_+(n,delta)) such that the probability of satisfiability is greater than 1-delta for alpha < alpha_- and is less than delta for alpha > alpha_+. We show that W(n,delta)=(1-Theta(n^{-1/3}),1+Theta(n^{-1/3})), where the constants implicit in Theta depend on delta. We also determine the rates at which the probability of satisfiability approaches one and zero at the boundaries of the window. Namely, for m=(1+epsilon)n, where epsilon may depend on n as long as |epsilon| is sufficiently small and |epsilon|*n^(1/3) is sufficiently large, we show that the probability of satisfiability decays like exp(-Theta(n*epsilon^3)) above the window, and goes to one like 1-Theta(1/(n*|epsilon|^3)) below the window. We prove these results by defining an order parameter for the transition and establishing its scaling behavior in n both inside and outside the window. Using this order parameter, we prove that the 2-SAT phase transition is continuous with an order parameter critical exponent of 1. We also determine the values of two other critical exponents, showing that the exponents of 2-SAT are identical to those of the random graph.

preprint1999arXiv

How to Couple from the Past Using a Read-Once Source of Randomness

We give a new method for generating perfectly random samples from the stationary distribution of a Markov chain. The method is related to coupling from the past (CFTP), but only runs the Markov chain forwards in time, and never restarts it at previous times in the past. The method is also related to an idea known as PASTA (Poisson arrivals see time averages) in the operations research literature. Because the new algorithm can be run using a read-once stream of randomness, we call it read-once CFTP. The memory and time requirements of read-once CFTP are on par with the requirements of the usual form of CFTP, and for a variety of applications the requirements may be noticeably less. Some perfect sampling algorithms for point processes are based on an extension of CFTP known as coupling into and from the past; for completeness, we give a read-once version of coupling into and from the past, but it remains unpractical. For these point process applications, we give an alternative coupling method with which read-once CFTP may be efficiently used.

preprint1999arXiv

Scaling Limits for Minimal and Random Spanning Trees in Two Dimensions

A general formulation is presented for continuum scaling limits of stochastic spanning trees. A spanning tree is expressed in this limit through a consistent collection of subtrees, which includes a tree for every finite set of endpoints in $\R^d$. Tightness of the distribution, as $δ\to 0$, is established for the following two-dimensional examples: the uniformly random spanning tree on $δ\Z^2$, the minimal spanning tree on $δ\Z^2$ (with random edge lengths), and the Euclidean minimal spanning tree on a Poisson process of points in $\R^2$ with density $δ^{-2}$. In each case, sample trees are proven to have the following properties, with probability one with respect to any of the limiting measures: i) there is a single route to infinity (as was known for $δ> 0$), ii) the tree branches are given by curves which are regular in the sense of Hölder continuity, iii) the branches are also rough, in the sense that their Hausdorff dimension exceeds one, iv) there is a random dense subset of $\R^2$, of dimension strictly between one and two, on the complement of which (and only there) the spanning subtrees are unique with continuous dependence on the endpoints, v) branching occurs at countably many points in $\R^2$, and vi) the branching numbers are uniformly bounded. The results include tightness for the loop erased random walk (LERW) in two dimensions. The proofs proceed through the derivation of scale-invariant power bounds on the probabilities of repeated crossings of annuli.