Researcher profile

A. K. Hartmann

A. K. Hartmann contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
17works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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)

preprint2014arXiv

Analysis of the phase transition in the $2D$ Ising ferromagnet using a Lempel-Ziv string parsing scheme and black-box data-compression utilities

In this work we consider information-theoretical observables to analyze short symbolic sequences, comprising time-series that represent the orientation of a single spin in a $2D$ Ising ferromagnet on a square lattice of size $L^2=128^2$, for different system temperatures $T$. The latter were chosen from an interval enclosing the critical point $T_{\rm c}$ of the model. At small temperatures the sequences are thus very regular, at high temperatures they are maximally random. In the vicinity of the critical point, nontrivial, long-range correlations appear. Here, we implement estimators for the entropy rate, excess entropy (i.e. "complexity") and multi-information. First, we implement a Lempel-Ziv string parsing scheme, providing seemingly elaborate entropy rate and multi-information estimates and an approximate estimator for the excess entropy. Furthermore, we apply easy-to-use black-box data compression utilities, providing approximate estimators only. For comparison and to yield results for benchmarking purposes we implement the information-theoretic observables also based on the well-established M-block Shannon entropy, which is more tedious to apply compared to the the first two "algorithmic" entropy estimation procedures. To test how well one can exploit the potential of such data compression techniques, we aim at detecting the critical point of the $2D$ Ising ferromagnet. Among the above observables, the multi-information, which is known to exhibit an isolated peak at the critical point, is very easy to replicate by means of both efficient algorithmic entropy estimation procedures. Finally, we assess how good the various algorithmic entropy estimates compare to the more conventional block entropy estimates and illustrate a simple modification that yields enhanced results.

preprint2013arXiv

Biased and greedy random walks on two-dimensional lattices with quenched randomness: the "greedy" ant within a disordered environment

The principle characteristics of biased greedy random walks (BGRWs) on two-dimensional lattices with real-valued quenched disorder on the lattice edges are studied. Here, the disorder allows for negative edge-weights. In previous studies, considering the negative-weight percolation (NWP) problem, this was shown to change the universality class of the existing, static percolation transition. In the presented study, four different types of BGRWs and an algorithm based on the ant colony optimization (ACO) heuristic were considered. Regarding the BGRWs, the precise configurations of the lattice walks constructed during the numerical simulations were influenced by two parameters: a disorder parameter rho that controls the amount of negative edge weights on the lattice and a bias strength B that governs the drift of the walkers along a certain lattice direction. Here, the pivotal observable is the probability that, after termination, a lattice walk exhibits a total negative weight, which is here considered as percolating. The behavior of this observable as function of rho for different bias strengths B is put under scrutiny. Upon tuning rho, the probability to find such a feasible lattice walk increases from zero to one. This is the key feature of the percolation transition in the NWP model. Here, we address the question how well the transition point rho_c, resulting from numerically exact and "static" simulations in terms of the NWP model can be resolved using simple dynamic algorithms that have only local information available, one of the basic questions in the physics of glassy systems.

preprint2013arXiv

Typical and large-deviation properties of minimum-energy paths on disordered hierarchical lattices

We perform numerical simulations to study the optimal path problem on disordered hierarchical graphs with effective dimension d=2.32. Therein, edge energies are drawn from a disorder distribution that allows for positive and negative energies. This induces a behavior which is fundamentally different from the case where all energies are positive, only. Upon changing the subtleties of the distribution, the scaling of the minimum energy path length exhibits a transition from self-affine to self-similar. We analyze the precise scaling of the path length and the associated ground-state energy fluctuations in the vincinity of the disorder critical point, using a decimation procedure for huge graphs. Further, using an importance sampling procedure in the disorder we compute the negative-energy tails of the ground-state energy distribution up to 12 standard deviations away from its mean. We find that the asymptotic behavior of the negative-energy tail is in agreement with a Tracy-Widom distribution. Further, the characteristic scaling of the tail can be related to the ground-state energy flucutations, similar as for the directed polymer in a random medium.

preprint2012arXiv

A computational mechanics approach to estimate entropy and (approximate) complexity for the dynamics of the 2D Ising Ferromagnet

We present a numerical analysis of the entropy rate and statistical complexity related to the spin flip dynamics of the 2D Ising Ferromagnet at different temperatures T. We follow an information theoretic approach and test three different entropy estimation algorithms to asses entropy rate and statistical complexity of binary sequences. The latter are obtained by monitoring the orientation of a single spin on a square lattice of side-length L=256 at a given temperature parameter over time. The different entropy estimation procedures are based on the M-block Shannon entropy (a well established method that yields results for benchmarking purposes), non-sequential recursive pair substitution (providing an elaborate and an approximate estimator) and a convenient data compression algorithm contained in the zlib-library (providing an approximate estimator only). We propose an approximate measure of statistical complexity that emphasizes on correlations within the sequence and which is easy to implement, even by means of black-box data compression algorithms. Regarding the 2D Ising Ferromagnet simulated using Metropolis dynamics and for binary sequences of finite length, the proposed approximate complexity measure is peaked close to the critical temperature. For the approximate estimators, a finite-size scaling analysis reveals that the peak approaches the critical temperature as the sequence length increases. Results obtained using different spin-flip dynamics are briefly discussed. The suggested complexity measure can be extended to non-binary sequences in a straightforward manner.

preprint2012arXiv

Analysis of the loop length distribution for the negative weight percolation problem in dimensions d=2 through 6

We consider the negative weight percolation (NWP) problem on hypercubic lattice graphs with fully periodic boundary conditions in all relevant dimensions from d=2 to the upper critical dimension d=6. The problem exhibits edge weights drawn from disorder distributions that allow for weights of either sign. We are interested in in the full ensemble of loops with negative weight, i.e. non-trivial (system spanning) loops as well as topologically trivial ("small") loops. The NWP phenomenon refers to the disorder driven proliferation of system spanning loops of total negative weight. While previous studies where focused on the latter loops, we here put under scrutiny the ensemble of small loops. Our aim is to characterize -using this extensive and exhaustive numerical study- the loop length distribution of the small loops right at and below the critical point of the hypercubic setups by means of two independent critical exponents. These can further be related to the results of previous finite-size scaling analyses carried out for the system spanning loops. For the numerical simulations we employed a mapping of the NWP model to a combinatorial optimization problem that can be solved exactly by using sophisticated matching algorithms. This allowed us to study here numerically exact very large systems with high statistics.

preprint2012arXiv

Information theoretic approach to ground-state phase transitions for two and three-dimensional frustrated spin systems

The information theoretic observables entropy (a measure of disorder), excess entropy (a measure of complexity) and multi information are used to analyze ground-state spin configurations for disordered and frustrated model systems in 2D and 3D. For both model systems, ground-state spin configurations can be obtained in polynomial time via exact combinatorial optimization algorithms, which allowed us to study large systems with high numerical accuracy. Both model systems exhibit a continuous transition from an ordered to a disordered ground state as a model parameter is varied. By using the above information theoretic observables it is possible to detect changes in the spatial structure of the ground states as the critical point is approached. It is further possible to quantify the scaling behavior of the information theoretic observables in the vicinity of the critical point. For both model systems considered, the estimates of critical properties for the ground-state phase transitions are in good agreement with existing results reported in the literature.

preprint2012arXiv

Is negative-weight percolation compatible with SLE?

We study numerically the geometrical properties of minimally weighted paths that appear in the negative-weight percolation (NWP) model on two-dimensional lattices assuming a combination of periodic and free boundary conditions (BCs). Each realization of the disorder consists of a random fraction 1-rho of bonds with unit strength and a fraction rho of bond strengths drawn from a Gaussian distribution with zero mean and unit width. For each such sample, the path is forced to span the lattice along the direction with the free BCs. The path and a set of negatively weighted loops form a ground state (GS). A ground state on such a lattice can be determined performing a non-trivial transformation of the original graph and applying sophisticated matching algorithms. Here we examine whether the geometrical properties of the paths are in accordance with predictions of Schramm-Loewner evolution (SLE). Measuring the fractal dimension and reviewing Schramm's left passage formula indicates that the paths cannot be described in terms of SLE.

preprint2011arXiv

Mean-field behavior of the negative-weight percolation model on random regular graphs

We investigate both analytically and numerically the ensemble of minimum-weight loops and paths in the negative-weight percolation model on random graphs with fixed connectivity and bimodal weight distribution. This allows us to study the mean-field behavior of this model. The analytical study is based on a conjectured equivalence with the problem of self-avoiding walks in a random medium. The numerical study is based on a mapping to a standard minimum-weight matching problem for which fast algorithms exist. Both approaches yield results which are in agreement, on the location of the phase transition, on the value of critical exponents, and on the absence of any sizeable indications of a glass phase. By these results, the previously conjectured upper critical dimension of d_u=6 is confirmed.

preprint2011arXiv

Optimal Vertex Cover for the Small-World Hanoi Networks

The vertex-cover problem on the Hanoi networks HN3 and HN5 is analyzed with an exact renormalization group and parallel-tempering Monte Carlo simulations. The grand canonical partition function of the equivalent hard-core repulsive lattice-gas problem is recast first as an Ising-like canonical partition function, which allows for a closed set of renormalization group equations. The flow of these equations is analyzed for the limit of infinite chemical potential, at which the vertex-cover problem is attained. The relevant fixed point and its neighborhood are analyzed, and non-trivial results are obtained both, for the coverage as well as for the ground state entropy density, which indicates the complex structure of the solution space. Using special hierarchy-dependent operators in the renormalization group and Monte-Carlo simulations, structural details of optimal configurations are revealed. These studies indicate that the optimal coverages (or packings) are not related by a simple symmetry. Using a clustering analysis of the solutions obtained in the Monte Carlo simulations, a complex solution space structure is revealed for each system size. Nevertheless, in the thermodynamic limit, the solution landscape is dominated by one huge set of very similar solutions.

preprint2010arXiv

A dedicated algorithm for calculating ground states for the triangular random bond Ising model

In the presented article we present an algorithm for the computation of ground state spin configurations for the 2d random bond Ising model on planar triangular lattice graphs. Therefore, it is explained how the respective ground state problem can be mapped to an auxiliary minimum-weight perfect matching problem, solvable in polynomial time. Consequently, the ground state properties as well as minimum-energy domain wall (MEDW) excitations for very large 2d systems, e.g. lattice graphs with up to N=384x384 spins, can be analyzed very fast. Here, we investigate the critical behavior of the corresponding T=0 ferromagnet to spin-glass transition, signaled by a breakdown of the magnetization, using finite-size scaling analyses of the magnetization and MEDW excitation energy and we contrast our numerical results with previous simulations and presumably exact results.

preprint2010arXiv

Configurational statistics of densely and fully packed loops in the negative-weight percolation model

By means of numerical simulations we investigate the configurational properties of densely and fully packed configurations of loops in the negative-weight percolation (NWP) model. In the presented study we consider 2d square, 2d honeycomb, 3d simple cubic and 4d hypercubic lattice graphs, where edge weights are drawn from a Gaussian distribution. For a given realization of the disorder we then compute a configuration of loops, such that the configurational energy, given by the sum of all individual loop weights, is minimized. For this purpose, we employ a mapping of the NWP model to the "minimum-weight perfect matching problem" that can be solved exactly by using sophisticated polynomial-time matching algorithms. We characterize the loops via observables similar to those used in percolation studies and perform finite-size scaling analyses, up to side length L=256 in 2d, L=48 in 3d and L=20 in 4d (for which we study only some observables), in order to estimate geometric exponents that characterize the configurations of densely and fully packed loops. One major result is that the loops behave like uncorrelated random walks from dimension d=3 on, in contrast to the previously studied behavior at the percolation threshold, where random-walk behavior is obtained for d>=6.

preprint2010arXiv

Large-deviation properties of largest component for random graphs

Distributions of the size of the largest component, in particular the large-deviation tail, are studied numerically for two graph ensembles, for Erdoes-Renyi random graphs with finite connectivity and for two-dimensional bond percolation. Probabilities as small as 10^-180 are accessed using an artificial finite-temperature (Boltzmann) ensemble. The distributions for the Erdoes-Renyi ensemble agree well with previously obtained analytical results. The results for the percolation problem, where no analytical results are available, are qualitatively similar, but the shapes of the distributions are somehow different and the finite-size corrections are sometimes much larger. Furthermore, for both problems, a first-order phase transition at low temperatures T within the artificial ensemble is found in the percolating regime, respectively.

preprint2010arXiv

Minimum-Free-Energy Distribution of RNA Secondary Structures: Entropic and Thermodynamic Properties of Rare Events

We study the distribution of the minimum free energy (MFE) for the Turner model of pseudoknot free RNA secondary structures over ensembles of random RNA sequences. In particular, we are interested in those rare and intermediate events of unexpected low MFEs. Generalized ensemble Markov-chain Monte Carlo methods allow us to explore the rare-event tail of the MFE distribution down to probabilities like $10^{-70}$ and to study the relationship between the sequence entropy and structural properties for sequence ensembles with fixed MFEs. Entropic and structural properties of those ensembles are compared with natural RNA of the same reduced MFE (z-score).

preprint2010arXiv

Numerical Solution-Space Analysis of Satisfiability Problems

The solution-space structure of the 3-Satisfiability Problem (3-SAT) is studied as a function of the control parameter alpha (ratio of number of clauses to the number of variables) using numerical simulations. For this purpose, one has to sample the solution space with uniform weight. It is shown here that standard stochastic local-search (SLS) algorithms like "ASAT" and "MCMCMC" (also known as "parallel tempering") exhibit a sampling bias. Nevertheless, unbiased samples of solutions can be obtained using the "ballistic-networking approach", which is introduced here. It is a generalization of "ballistic search" methods and yields also a cluster structure of the solution space. As application, solutions of 3-SAT instances are generated using ASAT plus ballistic networking. The numerical results are compatible with a previous analytic prediction of a simple solution-space structure for small values of alpha and a transition to a clustered phase at alpha_c ~ 3.86, where the solution space breaks up into several non-negligible clusters. Furthermore, in the thermodynamic limit there are, for values of alpha close to the SATUNSAT transition alpha_s ~ 4.267, always clusters without any frozen variables. This may explain why some SLS algorithms are able to solve very large 3-SAT instances close to the SAT-UNSAT transition.

preprint2010arXiv

The upper critical dimension of the negative-weight percolation problem

By means of numerical simulations we investigate the geometric properties of loops on hypercubic lattice graphs in dimensions d=2 through 7, where edge weights are drawn from a distribution that allows for positive and negative weights. We are interested in the appearance of system-spanning loops of total negative weight. The resulting negative-weight percolation (NWP) problem is fundamentally different from conventional percolation, as we have seen in previous studies of this model for the 2d case. Here, we characterize the transition for hypercubic systems, where the aim of the present study is to get a grip on the upper critical dimension d_u of the NWP problem. For the numerical simulations we employ a mapping of the NWP model to a combinatorial optimization problem that can be solved exactly by using sophisticated matching algorithms. We characterize the loops via observables similar to those in percolation theory and perform finite-size scaling analyses, e.g. 3d hypercubic systems with side length up to L=56 sites, in order to estimate the critical properties of the NWP phenomenon. We find our numerical results consistent with an upper critical dimension d_u=6 for the NWP problem.

preprint2009arXiv

Low-energy excitations in the three-dimensional random-field Ising model

The random-field Ising model (RFIM), one of the basic models for quenched disorder, can be studied numerically with the help of efficient ground-state algorithms. In this study, we extend these algorithm by various methods in order to analyze low-energy excitations for the three-dimensional RFIM with Gaussian distributed disorder that appear in the form of clusters of connected spins. We analyze several properties of these clusters. Our results support the validity of the droplet-model description for the RFIM.

preprint1995arXiv

Cluster-Exact Approximation of Spin Glass Groundstates

We present an algorithm which calculates groundstates of Ising spin glasses approximately. It works by randomly selecting clusters of spins which exhibit no frustrations. The spins which were not selected, contribute to the local fields of the selected spins. For the spin--cluster a groundstate is exactly calaculated by using graphtheoretical methods. The other spins remain unchanged. This procedure is repeated many times resulting in a state with low energy. The total time complexity of this scheme is approximately cubic. We estimate that the groundstate energy density of the infinite system for the +/- J model is -1.400 +/- 0.005 (2d) and -1.766 +/- 0.002 (3d). The distribution of overlaps for selected systems is calculated in order to characterize the algorithm.