GraphCombEx: A Software Tool for Exploration of Combinatorial Optimisation Properties of Large Graphs

We present a prototype of a software tool for exploration of multiple combinatorial optimisation problems in large real-world and synthetic complex networks. Our tool, called GraphCombEx (an acronym of Graph Combinatorial Explorer), provides a unified framework for scalable computation and presentation of high-quality suboptimal solutions and bounds for a number of widely studied combinatorial optimisation problems. Efficient representation and applicability to large-scale graphs and complex networks are particularly considered in its design. The problems currently supported include maximum clique, graph colouring, maximum independent set, minimum vertex clique covering, minimum dominating set, as well as the longest simple cycle problem. Suboptimal solutions and intervals for optimal objective values are estimated using scalable heuristics. The tool is designed with extensibility in mind, with the view of further problems and both new fast and high-performance heuristics to be added in the future. GraphCombEx has already been successfully used as a support tool in a number of recent research studies using combinatorial optimisation to analyse complex networks, indicating its promise as a research software tool

On kk-Core Percolation in Four Dimensions

The $k$-core percolation on the Bethe lattice has been proposed as a simple model of the jamming transition because of its hybrid first-order/second-order nature. We investigate numerically $k$-core percolation on the four-dimensional regular lattice. For $k=4$ the presence of a discontinuous transition is clearly established but its nature is strictly first order. In particular, the $k$-core density displays no singular behavior before the jump and its correlation length remains finite. For $k=3$ the transition is continuous

Remarks on Bootstrap Percolation in Metric Networks

We examine bootstrap percolation in d-dimensional, directed metric graphs in the context of recent measurements of firing dynamics in 2D neuronal cultures. There are two regimes, depending on the graph size N. Large metric graphs are ignited by the occurrence of critical nuclei, which initially occupy an infinitesimal fraction, f_* -> 0, of the graph and then explode throughout a finite fraction. Smaller metric graphs are effectively random in the sense that their ignition requires the initial ignition of a finite, unlocalized fraction of the graph, f_* >0. The crossover between the two regimes is at a size N_* which scales exponentially with the connectivity range \lambda like_* \sim \exp\lambda^d. The neuronal cultures are finite metric graphs of size N \simeq 10^5-10^6, which, for the parameters of the experiment, is effectively random since N<< N_*. This explains the seeming contradiction in the observed finite f_* in these cultures. Finally, we discuss the dynamics of the firing front

Facilitated spin models on Bethe lattice: bootstrap percolation, mode-coupling transition and glassy dynamics

We show that facilitated spin models of cooperative dynamics introduced by Fredrickson and Andersen display on Bethe lattices a glassy behaviour similar to the one predicted by the mode-coupling theory of supercooled liquids and the dynamical theory of mean-field disordered systems. At low temperature such cooperative models show a two-step relaxation and their equilibration time diverges at a finite temperature according to a power-law. The geometric nature of the dynamical arrest corresponds to a bootstrap percolation process which leads to a phase space organization similar to the one of mean-field disordered systems. The relaxation dynamics after a subcritical quench exhibits aging and converges asymptotically to the threshold states that appear at the bootstrap percolation transition.Comment: 7 pages, 6 figures, minor changes, final version to appear in Europhys. Let

Boson-exchange parquet solver for dual fermions

We present and implement a parquet approximation within the dual-fermion formalism based on a partial bosonization of the dual vertex function which substantially reduces the computational cost of the calculation. The method relies on splitting the vertex exactly into single-boson exchange contributions and a residual four-fermion vertex, which physically embody, respectively, long- and short-range spatial correlations. After recasting the parquet equations in terms of the residual vertex, these are solved using the truncated-unity method of Eckhardt et al. [Phys. Rev. B 101, 155104 (2020)2469-995010.1103/PhysRevB.101.155104], which allows for a rapid convergence with the number of form factors in different regimes. While our numerical treatment of the parquet equations can be restricted to only a few Matsubara frequencies, reminiscent of Astretsov et al. [Phys. Rev. B 101, 075109 (2020)2469-995010.1103/PhysRevB.101.075109], the one- and two-particle spectral information is fully retained. In applications to the two-dimensional Hubbard model the method agrees quantitatively with a stochastic summation of diagrams over a wide range of parameters

An Unusual Clinical Presentation of Merkel Cell Carcinoma: A Case Report

Introduction. Merkel cell carcinoma is a rare, aggressive neuroendocrine cell carcinoma arising in the sun-exposed skin of elderly patients. Most of these tumors are located in the dermis. An unusual clinical presentation of such a tumor in the subcutis, if not biopsied, may be easily mistaken as a benign lesion. Case Presentation. An 83-year-old white woman presented with a several-month history of a painless 7 mm subcutaneous mass that was initially thought to be a lipoma. A conservative follow-up was planned. At the insistence of the patient, an excisional biopsy of the mass was performed revealing a subcutaneous Merkel cell carcinoma. The tumor cells stained positively for CK 20, chromogranin, and synaptophysin. No other primary or metastatic tumors found after a thorough work-up. The patient was treated with local irradiation. She remains disease free at her six-month follow-up visit. Conclusion. When a new growth is encountered in the sun-exposed skin of elderly patients, a biopsy is warranted even if the lesion clinically appears benign

Hysteresis in the Random Field Ising Model and Bootstrap Percolation

We study hysteresis in the random-field Ising model with an asymmetric distribution of quenched fields, in the limit of low disorder in two and three dimensions. We relate the spin flip process to bootstrap percolation, and show that the characteristic length for self-averaging $L^*$ increases as $exp(exp (J/\Delta))$ in 2d, and as $exp(exp(exp(J/\Delta)))$ in 3d, for disorder strength $\Delta$ much less than the exchange coupling J. For system size $1 << L < L^*$, the coercive field $h_{coer}$ varies as $2J - \Delta \ln \ln L$ for the square lattice, and as $2J - \Delta \ln \ln \ln L$ on the cubic lattice. Its limiting value is 0 for L tending to infinity, both for square and cubic lattices. For lattices with coordination number 3, the limiting magnetization shows no jump, and $h_{coer}$ tends to J.Comment: 4 pages, 4 figure