Fixed Point Properties of the Ising Ferromagnet on the Hanoi Networks

The Ising model with ferromagnetic couplings on the Hanoi networks is analyzed with an exact renormalization group. In particular, the fixed-points are determined and the renormalization-group flow for certain initial conditions is analyzed. Hanoi networks combine a one-dimensional lattice structure with a hierarchy of small-world bonds to create a mix of geometric and mean-field properties. Generically, the small-world bonds result in non-universal behavior, i.e. fixed points and scaling exponents that depend on temperature and the initial choice of coupling strengths. It is shown that a diversity of different behaviors can be observed with seemingly small changes in the structure of the networks. Defining interpolating families of such networks, we find tunable transitions between regimes with power-law and certain essential singularities in the critical scaling of the correlation length, similar to the so-called inverted Berezinskii-Kosterlitz-Thouless transition previously observed only in scale-free or dense networks.Comment: 13 pages, revtex, 12 fig. incl.; fixed confusing labels, published version. For related publications, see http://www.physics.emory.edu/faculty/boettcher

Extremal Optimization at the Phase Transition of the 3-Coloring Problem

We investigate the phase transition of the 3-coloring problem on random graphs, using the extremal optimization heuristic. 3-coloring is among the hardest combinatorial optimization problems and is closely related to a 3-state anti-ferromagnetic Potts model. Like many other such optimization problems, it has been shown to exhibit a phase transition in its ground state behavior under variation of a system parameter: the graph's mean vertex degree. This phase transition is often associated with the instances of highest complexity. We use extremal optimization to measure the ground state cost and the ``backbone'', an order parameter related to ground state overlap, averaged over a large number of instances near the transition for random graphs of size $n$ up to 512. For graphs up to this size, benchmarks show that extremal optimization reaches ground states and explores a sufficient number of them to give the correct backbone value after about $O(n^{3.5})$ update steps. Finite size scaling gives a critical mean degree value $\alpha_{\rm c}=4.703(28)$. Furthermore, the exploration of the degenerate ground states indicates that the backbone order parameter, measuring the constrainedness of the problem, exhibits a first-order phase transition.Comment: RevTex4, 8 pages, 4 postscript figures, related information available at http://www.physics.emory.edu/faculty/boettcher

Renormalization Group for Critical Phenomena in Complex Networks

We discuss the behavior of statistical models on a novel class of complex “Hanoi” networks. Such modeling is often the cornerstone for the understanding of many dynamical processes in complex networks. Hanoi networks are special because they integrate small-world hierarchies common to many social and economical structures with the inevitable geometry of the real world these structures exist in. In addition, their design allows exact results to be obtained with the venerable renormalization group (RG). Our treatment will provide a detailed, pedagogical introduction to RG. In particular, we will study the Ising model with RG, for which the fixed points are determined and the RG flow is analyzed. We show that the small-world bonds result in non-universal behavior. It is shown that a diversity of different behaviors can be observed with seemingly small changes in the structure of hierarchical networks generally, and we provide a general theory to describe our findings

Extremal Optimization for Graph Partitioning

Extremal optimization is a new general-purpose method for approximating solutions to hard optimization problems. We study the method in detail by way of the NP-hard graph partitioning problem. We discuss the scaling behavior of extremal optimization, focusing on the convergence of the average run as a function of runtime and system size. The method has a single free parameter, which we determine numerically and justify using a simple argument. Our numerical results demonstrate that on random graphs, extremal optimization maintains consistent accuracy for increasing system sizes, with an approximation error decreasing over runtime roughly as a power law t^(-0.4). On geometrically structured graphs, the scaling of results from the average run suggests that these are far from optimal, with large fluctuations between individual trials. But when only the best runs are considered, results consistent with theoretical arguments are recovered.Comment: 34 pages, RevTex4, 1 table and 20 ps-figures included, related papers available at http://www.physics.emory.edu/faculty/boettcher

Sarma phase in relativistic and non-relativistic systems

We investigate the stability of the Sarma phase in two-component fermion systems in three spatial dimensions. For this purpose we compare strongly-correlated systems with either relativistic or non-relativistic dispersion relation: relativistic quarks and mesons at finite isospin density and spin-imbalanced ultracold Fermi gases. Using a Functional Renormalization Group approach, we resolve fluctuation effects onto the corresponding phase diagrams beyond the mean-field approximation. We find that fluctuations induce a second order phase transition at zero temperature, and thus a Sarma phase, in the relativistic setup for large isospin chemical potential. This motivates the investigation of the cold atoms setup with comparable mean-field phase structure, where the Sarma phase could then be realized in experiment. However, for the non-relativistic system we find the stability region of the Sarma phase to be smaller than the one predicted from mean-field theory. It is limited to the BEC side of the phase diagram, and the unitary Fermi gas does not support a Sarma phase at zero temperature. Finally, we propose an ultracold quantum gas with four fermion species that has a good chance to realize a zero-temperature Sarma phase.Comment: version published in Phys.Lett.B; 10 pages, 5 figure

Continuous extremal optimization for Lennard-Jones Clusters

In this paper, we explore a general-purpose heuristic algorithm for finding high-quality solutions to continuous optimization problems. The method, called continuous extremal optimization(CEO), can be considered as an extension of extremal optimization(EO) and is consisted of two components, one is with responsibility for global searching and the other is with responsibility for local searching. With only one adjustable parameter, the CEO's performance proves competitive with more elaborate stochastic optimization procedures. We demonstrate it on a well known continuous optimization problem: the Lennerd-Jones clusters optimization problem.Comment: 5 pages and 3 figure

Ordinary Percolation with Discontinuous Transitions

Percolation on a one-dimensional lattice and fractals such as the Sierpinski gasket is typically considered to be trivial because they percolate only at full bond density. By dressing up such lattices with small-world bonds, a novel percolation transition with explosive cluster growth can emerge at a nontrivial critical point. There, the usual order parameter, describing the probability of any node to be part of the largest cluster, jumps instantly to a finite value. Here, we provide a simple example of this transition in form of a small-world network consisting of a one-dimensional lattice combined with a hierarchy of long-range bonds that reveals many features of the transition in a mathematically rigorous manner.Comment: RevTex, 5 pages, 4 eps-figs, and Mathematica Notebook as Supplement included. Final version, with several corrections and improvements. For related work, see http://www.physics.emory.edu/faculty/boettcher

Condensation transition in a model with attractive particles and non-local hops

We study a one dimensional nonequilibrium lattice model with competing features of particle attraction and non-local hops. The system is similar to a zero range process (ZRP) with attractive particles but the particles can make both local and non-local hops. The length of the non-local hop is dependent on the occupancy of the chosen site and its probability is given by the parameter $p$. Our numerical results show that the system undergoes a phase transition from a condensate phase to a homogeneous density phase as $p$ is increased beyond a critical value $p_c$. A mean-field approximation does not predict a phase transition and describes only the condensate phase. We provide heuristic arguments for understanding the numerical results.Comment: 11 Pages, 6 Figures. Published in Journal of Statistical Mechanics: Theory and Experimen