Researcher profile

Stefan Boettcher

Stefan Boettcher contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
15topics
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

17 published item(s)

preprint2026arXiv

Instability cascades in crumpling mylar sheets follow a log-Poisson statistic

The process of aging following a hard quench into a glassy state is characterized universally, for a wide class of materials, by logarithmic evolution of state variables and a power-law decay of two-time correlation functions that collapse only for the ratio of those times. This stands in stark contrast with relaxation in equilibrium materials, where time-translational invariance holds. It is by now widely recognized that these aging processes, which ever so slowly relax a complex disordered material after a quench, are facilitated by activated events. Yet, theories often cited to describe such a non-equilibrium process can be shown to miss pertinent aspects that are inherent to many experiments. A case in point are recent experiments on crumpling sheets of mylar loaded by a weight whose acoustic emissions are measured while the material buckles. Using extensive simulations to generate long time-series of such buckling events, we show that crumpling is a log-Poisson process activated by increasingly rare record-sized fluctuations in a slowly stiffening material characterized by a logarithmically growing length-scale. Crumpling thus adds to a range of glassy materials exhibiting the log-Poisson property, which can be used to discriminate between theories.

preprint2022arXiv

Density Fluctuations in Granular Piles Traversing the Glass Transition: A Grain-Scale Characterization of the Transition via the Internal Energy

The transition into a glassy state of the ensemble of static, mechanically stable configurations of a tapped granular pile is explored using extensive molecular dynamics simulations. We show that different horizontal sub-regions ("layers") along the height of the pile traverse this transition in a similar manner but at distinct tap intensities. We supplement the conventional approach based purely on properties of the static configurations with investigations of the grain-scale dynamics by which the tap energy is transmitted throughout the pile. We find that the effective energy that particles dissipate is a function of each particle's location in the pile and, moreover, that its value plays a distinctive role in the transformation between configurations. This internal energy provides a "temperature-like" parameter that allows us to align the transition into the glassy state for all layers, as well as different annealing schedules, at a critical value.

preprint2022arXiv

Inability of a graph neural network heuristic to outperform greedy algorithms in solving combinatorial optimization problems like Max-Cut

In Nature Machine Intelligence 4, 367 (2022), Schuetz et al provide a scheme to employ graph neural networks (GNN) as a heuristic to solve a variety of classical, NP-hard combinatorial optimization problems. It describes how the network is trained on sample instances and the resulting GNN heuristic is evaluated applying widely used techniques to determine its ability to succeed. Clearly, the idea of harnessing the powerful abilities of such networks to ``learn'' the intricacies of complex, multimodal energy landscapes in such a hands-off approach seems enticing. And based on the observed performance, the heuristic promises to be highly scalable, with a computational cost linear in the input size $n$, although there is likely a significant overhead in the pre-factor due to the GNN itself. However, closer inspection shows that the reported results for this GNN are only minutely better than those for gradient descent and get outperformed by a greedy algorithm, for example, for Max-Cut. The discussion also highlights what I believe are some common misconceptions in the evaluations of heuristics.

preprint2022arXiv

Transport and Localization in Quantum Walks on a Random Hierarchy of Barriers

We study transport within a spatially heterogeneous one-dimensional quantum walk with a combination of hierarchical and random barriers. Recent renormalization group calculations for a spatially disordered quantum walk with a regular hierarchy of barriers alone have shown a gradual decrease in transport but no localization for increasing (but finite) barrier sizes. In turn, it is well-known that extensive random disorder in the spatial barriers is sufficient to localize a quantum walk on the line. Here we show that adding only a sparse (sub-extensive) amount of randomness to a hierarchy of barriers is sufficient to induce localization such that transport ceases. Our numerical results suggest the existence of a localization transition for a combination of both, the strength of the regular barrier hierarchy at large enough randomness as well as the increasing randomness at sufficiently strong barriers in the hierarchy.

preprint2021arXiv

Analysis of landscape hierarchy during coarsening and aging in Ising spin glasses

We use record dynamics (RD), a coarse-grained description of the ubiquitous relaxation phenomenology known as "aging", as a diagnostic tool to find universal features that distinguish between the energy landscapes of Ising spin models and the ferromagnet. According to RD, a non-equilibrium system after a quench relies on fluctuations that randomly generate a sequence of irreversible record-sized events (quakes or avalanches) that allow the system to escape ever-higher barriers of meta-stable states within a complex, hierarchical energy landscape. Once these record events allow the system to overcome such barriers, the system relaxes by tumbling into the following meta-stable state that is marginally more stable. Within this framework, a clear distinction can be drawn between the coarsening dynamics of an Ising ferromagnet and the aging of the spin glass, which are often put in the same category. To that end, we interpolate between the spin glass and ferromagnet by varying the admixture $p$ of ferromagnetic over anti-ferromagnetic bonds from the glassy state (at 50% each) to wherever clear ferromagnetic behavior emerges. The accumulation of record events grows logarithmic with time in the glassy regime, with a sharp transition at a specific admixture into the ferromagnetic regime where such activations saturate quickly. We show this effect both for the Edwards-Anderson model on a cubic lattice as well as the Sherrington-Kirkpatrick (mean-field) spin glass. While this transition coincides with a previously observed zero-temperature equilibrium transition in the former, that transition has not yet been described for the latter.

preprint2020arXiv

Ground State Properties of the Diluted Sherrington-Kirkpatrick Spin Glass

We present a numerical study of ground states of the dilute versions of the Sherrington-Kirkpatrick (SK) mean-field spin glass. In contrast to so-called "sparse" mean-field spin glasses that have been studied widely on random networks of finite (average or regular) degree, the networks studied here are randomly bond-diluted to an overall density $p$, such that the average degree diverges as $\sim pN$ with the system size $N$. Ground-state energies are obtained with high accuracy for random instances over a wide range of fixed $p$. Since this is a NP-hard combinatorial problem, we employ the Extremal Optimization heuristic to that end. We find that the exponent describing the finite-size corrections, $ω$, varies continuously with $p$, a somewhat surprising result, as one would not expect that gradual bond-dilution would change the $T=0$ universality class of a statistical model. For $p\to1$, the familiar result of $ω(p=1)\approx\frac{2}{3}$ for SK is obtained.

preprint2020arXiv

Quantum Ultra-Walks: Walks on a Line with Hierarchical Spatial Heterogeneity

We discuss the model of a one-dimensional, discrete-time walk on a line with spatial heterogeneity in the form of a variable set of ultrametric barriers. Inspired by the homogeneous quantum walk on a line, we develop a formalism by which the classical ultrametric random walk as well as the quantum walk can be treated in parallel by using a &#34;coined&#34; walk with internal degrees of freedom. For the random walk, this amounts to a $2^{\rm nd}$-order Markov process with a \emph{stochastic} coin, better known as an (anti-)persistent walk. When this coin varies spatially in the hierarchical manner of &#34;ultradiffusion,&#34; it reproduces the well-known results of that model. The exact analysis employed for obtaining the walk dimension $d_{w}$, based on the real-space renormalization group (RG), proceeds virtually identical for the corresponding quantum walk with a $unitary$ coin. However, while the classical walk remains robustly diffusive ($d_{w}=\frac{1}{2}$) for a wide range of barrier heights, unitarity provides for a quantum walk dimension $d_{w}$ that varies continuously, for even the smallest amount of heterogeneity, from ballistic spreading ($d_{w}=1$) in the homogeneous limit to confinement ($d_{w}=\infty$) for diverging barriers. Yet for any $d_{w}<\infty$ the quantum ultra-walk never appears to localize.

preprint2020arXiv

Universal Features of Annealing and Aging in Compaction of Granular Piles

We explore the compaction dynamics of a granular pile after a hard quench from a liquid into the glassy regime. First, we establish that the otherwise athermal granular pile during tapping exhibits annealing behavior comparable to glassy polymer or colloidal systems. Like those other systems, the pile undergoes a glass transition and &#34;freezes&#34; into different non-equilibrium glassy states at low agitation for different annealing speeds, starting from the same initial equilibrium state at high agitation. Then, we quench the system instantaneously from the highly-agitated state to below the glass transition regime to study the ensuing aging dynamics. In this classical aging protocol, the density increases (i.e., the potential energy of the pile decreases) logarithmically over several decades in time. Instead of system-wide, thermodynamic measures, here we identify the intermittent, irreversible events (&#34;quakes&#34;) that actually drive the glassy relaxation process. We find that the event rate decelerates hyperbolically, which explains the observed increase in density when the integrated contribution to the downward displacements is evaluated. We argue that such a hyperbolically decelerating event rate is consistent with a log-Poisson process, also found as a universal feature of aging in many thermal glasses.

preprint2019arXiv

Analysis of the Relation between Quadratic Unconstrained Binary Optimization (QUBO) and the Spin Glass Ground-State Problem

We analyze the transformation of QUBO from its conventional Boolean presentation into an equivalent spin glass problem with coupled $\pm1$ spin variables exposed to a site-dependent external field. We find that in a widely used testbed for QUBO these fields tend to be rather large compared to the typical coupling and many spins in each optimal configurations simply align with the fields irrespective of their constraints. Thereby, the testbed instances tend to exhibit large redundancies - seemingly independent variables which contribute little to the hardness of the problem, however. We demonstrate various consequences of this insight, for QUBO solvers as well as for heuristics developed for finding spin glass ground states. To this end, we implement the Extremal Optimization (EO) heuristic, in a new adaptation for the QUBO problem. We also propose a novel way to assess the quality of heuristics for increasing problem sizes based on asymptotic scaling.

preprint2013arXiv

Accumulation of beneficial mutations in one dimension

When beneficial mutations are relatively common, competition between multiple unfixed mutations can reduce the rate of fixation in well-mixed asexual populations. We introduce a one dimensional model with a steady accumulation of beneficial mutations. We find a transition between periodic selection and multiple-mutation regimes. In the multiple-mutation regime, the increase of fitness along the lattice bears a striking similarity to surface growth phenomena, with power law growth and saturation of the interface width. We also find significant differences compared to the well-mixed model. In our lattice model, the transition between regimes happens at a much lower mutation rate due to slower fixation times in one dimension. Also the rate of fixation is reduced with increasing mutation rate due to the more intense competition, and it saturates with large population size.

preprint2012arXiv

Finite-Size Corrections for Ground States of Edwards-Anderson Spin Glasses

Extensive computations of ground state energies of the Edwards-Anderson spin glass on bond-diluted, hypercubic lattices are conducted in dimensions d=3,..,7. Results are presented for bond-densities exactly at the percolation threshold, p=p_c, and deep within the glassy regime, p>p_c, where finding ground-states becomes a hard combinatorial problem. Finite-size corrections of the form 1/N^w are shown to be consistent throughout with the prediction w=1-y/d, where y refers to the &#34;stiffness&#34; exponent that controls the formation of domain wall excitations at low temperatures. At p=p_c, an extrapolation for $d\to\infty$ appears to match our mean-field results for these corrections. In the glassy phase, w does not approach the value of 2/3 for large d predicted from simulations of the Sherrington-Kirkpatrick spin glass. However, the value of w reached at the upper critical dimension does match certain mean-field spin glass models on sparse random networks of regular degree called Bethe lattices.

preprint2010arXiv

Simulations of Ground State Fluctuations in Mean-Field Ising Spin Glasses

The scaling of fluctuations in the distribution of ground-state energies or costs with the system size N for Ising spin glasses is considered using an extensive set of simulations with the Extremal Optimization heuristic across a range of different models on sparse and dense graphs. These models exhibit very diverse behaviors, and an asymptotic extrapolation is often complicated by higher-order corrections. The clearest picture, in fact, emerges from the study of graph-bipartitioning, a combinatorial optimization problem closely related to spin glasses. Aside from two-spin interactions with discrete bonds, we also consider problems with Gaussian bonds and three-spin interactions, which behave differently to a significant degree.

preprint2009arXiv

Conjecture on the maximum cut and bisection width in random regular graphs

Asymptotic properties of random regular graphs are object of extensive study in mathematics. In this note we argue, based on theory of spin glasses, that in random regular graphs the maximum cut size asymptotically equals the number of edges in the graph minus the minimum bisection size. Maximum cut and minimal bisection are two famous NP-complete problems with no known general relation between them, hence our conjecture is a surprising property of random regular graphs. We further support the conjecture with numerical simulations. A rigorous proof of this relation is obviously a challenge.

preprint2008arXiv

The Peculiar Phase Structure of Random Graph Bisection

The mincut graph bisection problem involves partitioning the n vertices of a graph into disjoint subsets, each containing exactly n/2 vertices, while minimizing the number of &#34;cut&#34; edges with an endpoint in each subset. When considered over sparse random graphs, the phase structure of the graph bisection problem displays certain familiar properties, but also some surprises. It is known that when the mean degree is below the critical value of 2 log 2, the cutsize is zero with high probability. We study how the minimum cutsize increases with mean degree above this critical threshold, finding a new analytical upper bound that improves considerably upon previous bounds. Combined with recent results on expander graphs, our bound suggests the unusual scenario that random graph bisection is replica symmetric up to and beyond the critical threshold, with a replica symmetry breaking transition possibly taking place above the threshold. An intriguing algorithmic consequence is that although the problem is NP-hard, we can find near-optimal cutsizes (whose ratio to the optimal value approaches 1 asymptotically) in polynomial time for typical instances near the phase transition.

preprint2005arXiv

Extremal Optimization for Sherrington-Kirkpatrick Spin Glasses

Extremal Optimization (EO), a new local search heuristic, is used to approximate ground states of the mean-field spin glass model introduced by Sherrington and Kirkpatrick. The implementation extends the applicability of EO to systems with highly connected variables. Approximate ground states of sufficient accuracy and with statistical significance are obtained for systems with more than N=1000 variables using $\pm J$ bonds. The data reproduces the well-known Parisi solution for the average ground state energy of the model to about 0.01%, providing a high degree of confidence in the heuristic. The results support to less than 1% accuracy rational values of $ω=2/3$ for the finite-size correction exponent, and of $ρ=3/4$ for the fluctuation exponent of the ground state energies, neither one of which has been obtained analytically yet. The probability density function for ground state energies is highly skewed and identical within numerical error to the one found for Gaussian bonds. But comparison with infinite-range models of finite connectivity shows that the skewness is connectivity-dependent.

preprint1995arXiv

SPHERICALLY SYMMETRIC RANDOM WALKS II. DIMENSIONALLY DEPENDENT CRITICAL BEHAVIOR

A recently developed model of random walks on a $D$-dimensional hyperspherical lattice, where $D$ is {\sl not} restricted to integer values, is extended to include the possibility of creating and annihilating random walkers. Steady-state distributions of random walkers are obtained for all dimensions $D>0$ by solving a discrete eigenvalue problem. These distributions exhibit dimensionally dependent critical behavior as a function of the birth rate. This remarkably simple model exhibits a second-order phase transition with a nontrivial critical exponent for all dimensions $D>0$.

preprint1995arXiv

SPHERICALLY SYMMETRIC RANDOM WALKS III. POLYMER ADSORPTION AT A HYPERSPHERICAL BOUNDARY

A recently developed model of random walks on a $D$-dimensional hyperspherical lattice, where $D$ is {\sl not} restricted to integer values, is used to study polymer growth near a $D$-dimensional attractive hyperspherical boundary. The model determines the fraction $P(κ)$ of the polymer adsorbed on this boundary as a function of the attractive potential $κ$ for all values of $D$. The adsorption fraction $P(κ)$ exhibits a second-order phase transition with a nontrivial scaling coefficient for $0<D<4$, $D\neq 2$, and exhibits a first-order phase transition for $D>4$. At $D=4$ there is a tricritical point with logarithmic scaling. This model reproduces earlier results for $D=1$ and $D=2$, where $P(κ)$ scales linearly and exponentially, respectively. A crossover transition that depends on the radius of the adsorbing boundary is found.