Source author record

Sergio Caracciolo

Sergio Caracciolo 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

34works
20topics
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

34 published item(s)

preprint2023arXiv

Simulated annealing, optimization, searching for ground states

The chapter starts with a historical summary of first attempts to optimize the spin glass Hamiltonian, comparing it to recent results on searching largest cliques in random graphs. Exact algorithms to find ground states in generic spin glass models are then explored in Section 1.2, while Section 1.3 is dedicated to the bidimensional case where polynomial algorithms exist and allow for the study of much larger systems. Finally Section 1.4 presents a summary of results for the assignment problem where the finite size corrections for the ground state can be studied in great detail.

preprint2020arXiv

The Dyck bound in the concave 1-dimensional random assignment model

We consider models of assignment for random $N$ blue points and $N$ red points on an interval of length $2N$, in which the cost for connecting a blue point in $x$ to a red point in $y$ is the concave function $|x-y|^p$, for $0<p<1$. Contrarily to the convex case $p>1$, where the optimal matching is trivially determined, here the optimization is non-trivial. The purpose of this paper is to introduce a special configuration, that we call the \emph{Dyck matching}, and to study its statistical properties. We compute exactly the average cost, in the asymptotic limit of large $N$, together with the first subleading correction. The scaling is remarkable: it is of order $N$ for $p<\frac{1}{2}$, order $N \ln N$ for $p=\frac{1}{2}$, and $N^{\frac{1}{2}+p}$ for $p>\frac{1}{2}$, and it is universal for a wide class of models. We conjecture that the average cost of the Dyck matching has the same scaling in $N$ as the cost of the optimal matching, and we produce numerical data in support of this conjecture. We hope to produce a proof of this claim in future work.

preprint2015arXiv

Deterministic Abelian Sandpile and square-triangle tilings

The Abelian Sandpile Model, seen as a deterministic lattice automaton, on two-dimensional periodic graphs generates complex regular patterns displaying (fractal) self-similarity. In particular, on a variety of lattices and initial conditions, at all sizes, there appears what we call an exact Sierpinski structure: the volume is filled with periodic patterns, glued together along straight lines, with the topology of a triangular Sierpinski gasket. Various lattices (square, hexagonal, kagome,...), initial conditions, and toppling rules show Sierpinski structures which are apparently unrelated and involve different mechanisms. As will be shown elsewhere, all these structures fall under one roof, and are in fact different projections of a unique mechanism pertinent to a family of deterministic surfaces in a 4-dimensional lattice. This short note gives a description of this surface, and of the combinatorics associated to its construction.

preprint2015arXiv

Quadratic stochastic Euclidean bipartite matching problem

We propose a new approach for the study of the quadratic stochastic Euclidean bipartite matching problem between two sets of $N$ points each, $N\gg 1$. The points are supposed independently randomly generated on a domain $Ω\subset\mathbb R^d$ with a given distribution $ρ(\mathbf x)$ on $Ω$. In particular, we derive a general expression for the correlation function and for the average optimal cost of the optimal matching. A previous ansatz for the matching problem on the flat hypertorus is obtained as particular case.

preprint2015arXiv

Scaling hypothesis for the Euclidean bipartite matching problem II. Correlation functions

We analyze the random Euclidean bipartite matching problem on the hypertorus in $d$ dimensions with quadratic cost and we derive the two--point correlation function for the optimal matching, using a proper ansatz introduced by Caracciolo et al. to evaluate the average optimal matching cost. We consider both the grid--Poisson matching problem and the Poisson--Poisson matching problem. We also show that the correlation function is strictly related to the Green's function of the Laplace operator on the hypertorus.

preprint2014arXiv

Correlation function for the Grid-Poisson Euclidean matching on a line and on a circle

We compute the two-point correlation function for spin configurations which are obtained by solving the Euclidean matching problem, for one family of points on a grid, and the second family chosen uniformly at random, when the cost depends on a power $p$ of the Euclidean distance. We provide the analytic solution in the thermodynamic limit, in a number of cases ($p>1$ open b.c.\ and $p=2$ periodic b.c., both at criticality), and analyse numerically other parts of the phase diagram.

preprint2014arXiv

On the one dimensional Euclidean matching problem: exact solutions, correlation functions and universality

We discuss the equivalence relation between the Euclidean bipartite matching problem on the line and on the circumference and the Brownian bridge process on the same domains. The equivalence allows us to compute the correlation function and the optimal cost of the original combinatoric problem in the thermodynamic limit; moreover, we solve also the minimax problem on the line and on the circumference. The properties of the average cost and correlation functions are discussed.

preprint2014arXiv

Replica Symmetry Breaking in Cold Atoms and Spin Glasses

We consider a system composed by N atoms trapped within a multimode cavity, whose theoretical description is captured by a disordered multimode Dicke model. We show that in the resonant, zero field limit the system exactly realizes the Sherrington-Kirkpatrick model. Upon a redefinition of the temperature, the same dynamics is realized in the dispersive, strong field limit. This regime also gives access to spin-glass observables which can be used to detect Replica Symmetry Breaking.

preprint2014arXiv

Scaling hypothesis for the Euclidean bipartite matching problem

We propose a simple yet very predictive form, based on a Poisson's equation, for the functional dependence of the cost from the density of points in the Euclidean bipartite matching problem. This leads, for quadratic costs, to the analytic prediction of the large $N$ limit of the average cost in dimension $d=1,2$ and of the subleading correction in higher dimension. A non-trivial scaling exponent, $γ_d=\frac{d-2}{d}$, which differs from the monopartite's one, is found for the subleading correction. We argue that the same scaling holds true for a generic cost exponent in dimension $d>2$.

preprint2013arXiv

Growth and form of melanoma cell colonies

We study the statistical properties of melanoma cell colonies grown in vitro by analyzing the results of crystal violet assays at different concentrations of initial plated cells and for different growth times. The distribution of colony sizes is described well by a continuous time branching process. To characterize the shape fluctuations of the colonies, we compute the distribution of eccentricities. The experimental results are compared with numerical results for models of random division of elastic cells, showing that experimental results are best reproduced by restricting cell division to the outer rim of the colony. Our results serve to illustrate the wealth of information that can be extracted by a standard experimental method such as the crystal violet assay.

preprint2013arXiv

Noncommutative determinants, Cauchy-Binet formulae, and Capelli-type identities II. Grassmann and quantum oscillator algebra representation

We prove that, for $X$, $Y$, $A$ and $B$ matrices with entries in a non-commutative ring such that $[X_{ij},Y_{k\ell}]=-A_{i\ell} B_{kj}$, satisfying suitable commutation relations (in particular, $X$ is a Manin matrix), the following identity holds: $ \mathrm{coldet} X \mathrm{coldet} Y = < 0 | \mathrm{coldet} (a A + X (I-a^{\dagger} B)^{-1} Y) |0 > $. Furthermore, if also $Y$ is a Manin matrix, $ \mathrm{coldet} X \mathrm{coldet} Y =\int \mathcal{D}(ψ, ψ^{\dagger}) \exp [ \sum_{k \geq 0} \frac{1}{k+1} (ψ^{\dagger} A ψ)^{k} (ψ^{\dagger} X B^k Y ψ) ] $. Notations: $ < 0 |$, $| 0 >$, are respectively the bra and the ket of the ground state, $a^{\dagger}$ and $a$ the creation and annihilation operators of a quantum harmonic oscillator, while $ψ^{\dagger}_i$ and $ψ_i$ are Grassmann variables in a Berezin integral. These results should be seen as a generalization of the classical Cauchy-Binet formula, in which $A$ and $B$ are null matrices, and of the non-commutative generalization, the Capelli identity, in which $A$ and $B$ are identity matrices and $[X_{ij},X_{k\ell}]=[Y_{ij},Y_{k\ell}]=0$.

preprint2012arXiv

Algebraic/combinatorial proofs of Cayley-type identities for derivatives of determinants and pfaffians

The classic Cayley identity states that \det(\partial) (\det X)^s = s(s+1)...(s+n-1) (\det X)^{s-1} where X=(x_{ij}) is an n-by-n matrix of indeterminates and \partial=(\partial/\partial x_{ij}) is the corresponding matrix of partial derivatives. In this paper we present straightforward combinatorial proofs of a variety of Cayley-type identities, both old and new. The most powerful of these proofs employ Grassmann algebra (= exterior algebra) and Grassmann-Berezin integration. Among the new identities proven here are a pair of "diagonal-parametrized" Cayley identities, a pair of "Laplacian-parametrized" Cayley identities, and the "product-parametrized" and "border-parametrized" rectangular Cayley identities.

preprint2012arXiv

Exact integration of height probabilities in the Abelian Sandpile Model

The height probabilities for the recurrent configurations in the Abelian Sandpile Model on the square lattice have analytic expressions, in terms of multidimensional quadratures. At first, these quantities have been evaluated numerically with high accuracy, and conjectured to be certain cubic rational-coefficient polynomials in 1/pi. Later their values have been determined by different methods. We revert to the direct derivation of these probabilities, by computing analytically the corresponding integrals. Yet another time, we confirm the predictions on the probabilities, and thus, as a corollary, the conjecture on the average height.

preprint2012arXiv

Transfer matrix for Kogut-Susskind fermions in the spin basis

In the absence of interaction it is well known that the Kogut-Susskind regularizations of fermions in the spin and flavor basis are equivalent to each other. In this paper we clarify the difference between the two formulations in the presence of interaction with gauge fields. We then derive an explicit expression of the transfer matrix in the spin basis by a unitary transformation on that one in the flavor basis which is known. The essential key ingredient is the explicit construction of the fermion Fock space for variables which live on blocks. Therefore the transfer matrix generates time translations of two lattice units.

preprint2011arXiv

Chiral symmetry breaking and quark confinement in the nilpotency expansion of QCD

We apply to lattice QCD a bosonization method previously developed in which dynamical bosons are generated by time-dependent Bogoliubov transformations. The transformed action can be studied by an expansion in the inverse of the nilpotency index, which is the number of fermionic states in the structure function of composite bosons. When this number diverges the model is solved by the saddle point method which has a variational interpretation. We give a stationary covariant solution for a background matter field whose fluctuations describe mesons. In the saddle point approximations live fermionic quasiparticles with quark quantum numbers which are confined, in the sense that they propagate only in pointlike color singlets. Conditions for chiral symmetry breaking are determined, to be studied numerically, and a derivation of mesons-nucleons action is outlined.

preprint2011arXiv

Diquarks in the nilpotency expansion of QCD and their role at finite chemical potential

We assume that the most important quark correlations are pairwise at all baryon densities. We introduce correlated pairs by means of Bogoliubov transformations which are functions of time and spatial gauge fields, in the formalism of the transfer matrix with lattice regularization. The dependence on time and gauge fields allows us to enforce gauge invariance and other symmetries term by term in the transformed quantities. The resulting action should be suitable for the description of multiquark mesons and baryons as states of a quark and a diquark. We derive the quark contribution to the free energy at finite chemical potential in a certain approximation. Its expression cannot be evaluated analytically, but it has a definite sign.

preprint2011arXiv

Multiple and inverse topplings in the Abelian Sandpile Model

The Abelian Sandpile Model is a cellular automaton whose discrete dynamics reaches an out-of-equilibrium steady state resembling avalanches in piles of sand. The fundamental moves defining the dynamics are encoded by the toppling rules. The transition monoid corresponding to this dynamics in the set of stable configurations is abelian, a property which seems at the basis of our understanding of the model. By including also antitoppling rules, we introduce and investigate a larger monoid, which is not abelian anymore. We prove a number of algebraic properties of this monoid, and describe their practical implications on the emerging structures of the model.

preprint2010arXiv

Absence of sign problem in the (saddle point approximation of the) nilpotency expansion of QCD at finite chemical potential

We have developed a method to derive the (approximate) quark contribution to the fermion free energy of QCD on a lattice, at finite temperature and chemical potential, with Kogut-Susskind fermions in the flavor basis. We show here the expression at zero temperature. This result has been obtained at the lowest order of the nilpotency expansion. At this order the well known "sign problem" does not arise and the quark contribution to the action can be used as a statistical weight in the Monte Carlo simulations.

preprint2010arXiv

Conservation laws for strings in the Abelian Sandpile Model

The Abelian Sandpile generates complex and beautiful patterns and seems to display allometry. On the plane, beyond patches, patterns periodic in both dimensions, we remark the presence of structures periodic in one dimension, that we call strings. We classify completely their constituents in terms of their principal periodic vector k, that we call momentum. We derive a simple relation between the momentum of a string and its density of particles, E, which is reminiscent of a dispersion relation, E=k^2. Strings interact: they can merge and split and within these processes momentum is conserved. We reveal the role of the modular group SL(2,Z) behind these laws.

preprint2010arXiv

Geometrical Properties of Two-Dimensional Interacting Self-Avoiding Walks at the Theta-Point

We perform a Monte Carlo simulation of two-dimensional N-step interacting self-avoiding walks at the theta point, with lengths up to N=3200. We compute the critical exponents, verifying the Coulomb-gas predictions, the theta-point temperature T_theta = 1.4986(11), and several invariant size ratios. Then, we focus on the geometrical features of the walks, computing the instantaneous shape ratios, the average asphericity, and the end-to-end distribution function. For the latter quantity, we verify in detail the theoretical predictions for its small- and large-distance behavior.

preprint2009arXiv

Noncommutative determinants, Cauchy-Binet formulae, and Capelli-type identities. I. Generalizations of the Capelli and Turnbull identities

We prove, by simple manipulation of commutators, two noncommutative generalizations of the Cauchy-Binet formula for the determinant of a product. As special cases we obtain elementary proofs of the Capelli identity from classical invariant theory and of Turnbull's Capelli-type identities for symmetric and antisymmetric matrices.

preprint2008arXiv

Bogoliubov transformations and fermion condensates in lattice field theories

We apply generalized Bogoliubov transformations to the transfer matrix of relativistic field theories regularized on a lattice. We derive the conditions these transformations must satisfy to factorize the transfer matrix into two terms which propagate fermions and antifermions separately, and we solve the relative equations under some conditions. We relate these equations to the saddle point approximation of a recent bosonization method and to the Foldy-Wouthuysen transformations which separate positive from negative energy states in the Dirac Hamiltonian.

preprint2007arXiv

Composite boson dominance in relativistic field theories

We apply a new bosonization technique to relativistic field theories of fermions whose partition function is dominated by bosonic composites, and derive the effective action for these bosons. The derivation respects all symmetries, including gauge invariance, with the exception of Euclidean invariance which must be checked a posteriori. We use a lattice regularization which should make applications to gauge theories easier. We test the method on a fermion field theory with quartic interaction in the limit when the number of flavours N_f is large, and show that it reproduces the exact results in the bosonic sector, namely condensation of a compositeboson with the right mass which breaks the discrete chiral invariance of the model. Moreover we determine the structure function of the condensed composite, whose spatial part turns out to be identical to that of the Cooper pairs of the BCS model of superconductivity.

preprint2005arXiv

Correction-to-scaling exponents for two-dimensional self-avoiding walks

We study the correction-to-scaling exponents for the two-dimensional self-avoiding walk, using a combination of series-extrapolation and Monte Carlo methods. We enumerate all self-avoiding walks up to 59 steps on the square lattice, and up to 40 steps on the triangular lattice, measuring the mean-square end-to-end distance, the mean-square radius of gyration and the mean-square distance of a monomer from the endpoints. The complete endpoint distribution is also calculated for self-avoiding walks up to 32 steps (square) and up to 22 steps (triangular). We also generate self-avoiding walks on the square lattice by Monte Carlo, using the pivot algorithm, obtaining the mean-square radii to ~0.01% accuracy up to N = 4000. We give compelling evidence that the first non-analytic correction term for two-dimensional self-avoiding walks is Delta_1 = 3/2. We compute several moments of the endpoint distribution function, finding good agreement with the field-theoretic predictions. Finally, we study a particular invariant ratio that can be shown, by conformal-field-theory arguments, to vanish asymptotically, and we find the cancellation of the leading analytic correction.

preprint2004arXiv

Two-Dimensional Heisenberg Model with Nonlinear Interactions: 1/N Corrections

We investigate a two-dimensional classical $-vector model with a generic nearest-neighbor interaction $W(\bsigma_i\cdot \bsigma_j)$ in the large-N limit, focusing on the finite-temperature transition point at which energy-energy correlations become critical. We show that this transition belongs to the Ising universality class. However, the width of the region in which Ising behavior is observed scales as $1/N^{3/2}$ along the magnetic direction and as 1/N in the thermal direction; outside a crossover to mean-field behavior occurs. This explains why only mean-field behavior is observed for $N=\infty$

preprint2003arXiv

Comment on ``Dynamic behavior of anisotropic non-equilibrium driving lattice gases''

In a recent Letter Albano and Saracco study the dynamic critical behavior of some anisotropic driven lattice gases by Monte Carlo (MC) simulations. In this Comment we point out that the Ansätze they use to relate the measured scaling exponents with the critical exponents analytically computed within different field-theoretical approaches do not take properly into account the strongly anisotropic nature of the phase transition, by implicitly assuming $z = z_{\bot} = z_{\parallel}$. As a consequence, at variance with the claims by the authors, their MC data are not conclusive to determine which one of the field theories proposed in the literature correctly describes the universal properties of the phase transition in these lattice gases.

preprint2002arXiv

Two-Dimensional Heisenberg Model with Nonlinear Interactions

We investigate a two-dimensional classical $N$-vector model with a nonlinear interaction $(1 + \bsigma_i\cdot \bsigma_j)^p$ in the large-N limit. As observed for N=3 by Blöte {\em et al.} [Phys. Rev. Lett. {\bf 88}, 047203 (2002)], we find a first-order transition for $p>p_c$ and no finite-temperature phase transitions for $p < p_c$. For $p>p_c$, both phases have short-range order, the correlation length showing a finite discontinuity at the transition. For $p=p_c$, there is a peculiar transition, where the spin-spin correlation length is finite while the energy-energy correlation length diverges.

preprint2000arXiv

Finite-Size Correlation Length and Violations of Finite-Size Scaling

We address the problem of the definition of the finite-volume correlation length. First, we study the large-N limit of the N-vector model, and we show the existence of several constraints on the definition if regularity of the finite-size scaling functions and correct anomalous behaviour above the upper critical dimension are required. Then, we study in detail a model in which the zero mode in prohibited. Such a model is a generalization of the fixed-magnetization Ising model which is equivalent to the lattice gas. Also in this case, we find that the finite-volume correlation length must satisfy appropriate constraints in order to obtain regular finite-size scaling functions, and, above the upper critical dimension, an anomalous scaling behaviour. The large-N results are confirmed by a one-loop calculation in the lattice $ϕ^4$ theory.

preprint2000arXiv

Operator Product Expansion on the Lattice: a Numerical Test in the Two-Dimensional Non-Linear Sigma-Model

We consider the short-distance behaviour of the product of the Noether O(N) currents in the lattice nonlinear sigma-model. We compare the numerical results with the predictions of the operator product expansion, using one-loop perturbative renormalization-group improved Wilson coefficients. We find that, even on quite small lattices (m a \approx 1/6), the perturbative operator product expansion describes that data with an error of 5-10% in a large window 2a \ltapprox x \ltapprox m^{-1}. We present a detailed discussion of the possible systematic errors.

preprint1998arXiv

Corrections to Finite-Size Scaling in the Lattice N-Vector Model for Infinite N

We compute the corrections to finite-size scaling for the N-vector model on the square lattice in the large-N limit. We find that corrections behave as log L/L^2. For tree-level improved hamiltonians corrections behave as 1/L^2. In general l-loop improvement is expected to reduce this behaviour to 1/(L^2 \log^l L). We show that the finite-size-scaling and the perturbative limit do not commute in the calculation of the corrections to finite-size scaling. We present also a detailed study of the corrections for the RP^N-model.

preprint1993arXiv

A General Limitation on Monte Carlo Algorithms of Metropolis Type

We prove that for any Monte Carlo algorithm of Metropolis type, the autocorrelation time of a suitable ``energy''-like observable is bounded below by a multiple of the corresponding ``specific heat''. This bound does not depend on whether the proposed moves are local or non-local; it depends only on the distance between the desired probability distribution $π$ and the probability distribution $π^{(0)}$ for which the proposal matrix satisfies detailed balance. We show, with several examples, that this result is particularly powerful when applied to non-local algorithms.