Source author record

Jon Machta

Jon Machta appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

10 published item(s)

preprint2015arXiv

Chaos in spin glasses revealed through thermal boundary conditions

We study the fragility of spin glasses to small temperature perturbations numerically using population annealing Monte Carlo. We apply thermal boundary conditions to a three-dimensional Edwards-Anderson Ising spin glass. In thermal boundary conditions all eight combinations of periodic versus antiperiodic boundary conditions in the three spatial directions are present, each appearing in the ensemble with its respective statistical weight determined by its free energy. We show that temperature chaos is revealed in the statistics of crossings in the free energy for different boundary conditions. By studying the energy difference between boundary conditions at free-energy crossings, we determine the domain-wall fractal dimension. Similarly, by studying the number of crossings, we determine the chaos exponent. Our results also show that computational hardness in spin glasses and the presence of chaos are closely related.

preprint2012arXiv

Computational Study of a Multistep Height Model

An equilibrium random surface multistep height model proposed in [Abraham and Newman, EPL, 86, 16002 (2009)] is studied using a variant of the worm algorithm. In one limit, the model reduces to the two-dimensional Ising model in the height representation. When the Ising model constraint of single height steps is relaxed, the critical temperature and critical exponents are continuously varying functions of the parameter controlling height steps larger than one. Numerical estimates of the critical exponents can be mapped via a single parameter-- the Coulomb gas coupling-- to the exponents of the O(n) loop model on the honeycomb lattice with n <= 1.

preprint2012arXiv

Packing Squares in a Torus

The densest packings of N unit squares in a torus are studied using analytical methods as well as simulated annealing. A rich array of dense packing solutions are found: density-one packings when N is the sum of two square integers; a family of "gapped bricklayer" Bravais lattice solutions with density N/(N+1); and some surprising non-Bravais lattice configurations, including lattices of holes as well as a configuration for N=23 in which not all squares share the same orientation. The entropy of some of these configurations and the frequency and orientation of density-one solutions as N goes to infinity are discussed.

preprint2011arXiv

Monte Carlo Methods for Rough Free Energy Landscapes: Population Annealing and Parallel Tempering

Parallel tempering and population annealing are both effective methods for simulating equilibrium systems with rough free energy landscapes. Parallel tempering, also known as replica exchange Monte Carlo, is a Markov chain Monte Carlo method while population annealing is a sequential Monte Carlo method. Both methods overcome the exponential slowing associated with high free energy barriers. The convergence properties and efficiency of the two methods are compared. For large systems, population annealing initially converges to equilibrium more rapidly than parallel tempering for the same amount of computational work. However, parallel tempering converges exponentially and population annealing inversely in the computational work so that ultimately parallel tempering approaches equilibrium more rapidly than population annealing.

preprint2011arXiv

Natural Complexity, Computational Complexity and Depth

Depth is a complexity measure for natural systems of the kind studied in statistical physics and is defined in terms of computational complexity. Depth quantifies the length of the shortest parallel computation required to construct a typical system state or history starting from simple initial conditions. The properties of depth are discussed and it is compared to other complexity measures. Depth can only be large for systems with embedded computation.

preprint2011arXiv

Parallel Complexity of Random Boolean Circuits

Random instances of feedforward Boolean circuits are studied both analytically and numerically. Evaluating these circuits is known to be a P-complete problem and thus, in the worst case, believed to be impossible to perform, even given a massively parallel computer, in time much less than the depth of the circuit. Nonetheless, it is found that for some ensembles of random circuits, saturation to a fixed truth value occurs rapidly so that evaluation of the circuit can be accomplished in much less parallel time than the depth of the circuit. For other ensembles saturation does not occur and circuit evaluation is apparently hard. In particular, for some random circuits composed of connectives with five or more inputs, the number of true outputs at each level is a chaotic sequence. Finally, while the average case complexity depends on the choice of ensemble, it is shown that for all ensembles it is possible to simultaneously construct a typical circuit together with its solution in polylogarithmic parallel time.

preprint2010arXiv

Population Annealing with Weighted Averages: A Monte Carlo Method for Rough Free Energy Landscapes

The population annealing algorithm introduced by Hukushima and Iba is described. Population annealing combines simulated annealing and Boltzmann weighted differential reproduction within a population of replicas to sample equilibrium states. Population annealing gives direct access to the free energy. It is shown that unbiased measurements of observables can be obtained by weighted averages over many runs with weight factors related to the free energy estimate from the run. Population annealing is well suited to parallelization and may be a useful alternative to parallel tempering for systems with rough free energy landscapes such as spin glasses. The method is demonstrated for spin glasses.

preprint2009arXiv

Strengths and Weaknesses of Parallel Tempering

Parallel tempering, also known as replica exchange Monte Carlo, is studied in the context of two simple free energy landscapes. The first is a double well potential defined by two macrostates separated by a barrier. The second is a `golf course' potential defined by microstates having two possible energies with exponentially more high energy states than low energy states. The equilibration time for replica exchange is analyzed for both systems. For the double well system, parallel tempering with a number of replicas that scales as the square root of the barrier height yields exponential speedup of the equilibration time. On the other hand, replica exchange yields only marginal speed-up for the golf course system. For the double well system, the free energy difference between the two wells has a large effect on the equilibration time. Nearly degenerate wells equilibrate much more slowly than strongly asymmetric wells. It is proposed that this difference in equilibration time may lead to a bias in measuring overlaps in spin glasses. These examples illustrate the strengths and weaknesses of replica exchange and may serve as a guide for understanding and improving the method in various applications.

preprint2000arXiv

Critical dynamics of two-replica cluster algorithms

The dynamic critical behavior of the two-replica cluster algorithm is studied. Several versions of the algorithm are applied to the two-dimensional, square lattice Ising model with a staggered field. The dynamic exponent for the full algorithm is found to be less than 0.4. It is found that odd translations of one replica with respect to the other together with global flips are essential for obtaining a small value of the dynamic exponent.