Showing 1–44 of 44 results for author: Angelini, C

Bethe $M$-layer construction for the percolation problem

Authors: Maria Chiara Angelini, Saverio Palazzi, Tommaso Rizzo, Marco Tarzia

Abstract: The major difference between percolation and other phase transition models is the absence of an Hamiltonian and of a partition function. For this reason it is not straightforward to identify the corresponding field theory to be used as starting point of Renormalization Group computations. Indeed, it could be identified with the field theory of $n+1$ states Potts model in the limit of $n \to 0$ onl… ▽ More The major difference between percolation and other phase transition models is the absence of an Hamiltonian and of a partition function. For this reason it is not straightforward to identify the corresponding field theory to be used as starting point of Renormalization Group computations. Indeed, it could be identified with the field theory of $n+1$ states Potts model in the limit of $n \to 0$ only by means of the mapping discovered by Kasteleyn and Fortuin for bond percolation. In this paper we show that it is possible to recover the epsilon expansion for critical exponents in finite dimension directly using the $M$-layer expansion, without the need to perform any analytical continuation. Moreover, we also show explicitly that the critical exponents for site and bond percolation are the same. This computation provides a reference for applications of the $M$-layer method to systems where the underlying field theory is unknown or disputed. △ Less

Submitted 29 June, 2024; originally announced July 2024.

Bethe $M$-layer construction on the Ising model

Authors: Saverio Palazzi, Maria Chiara Angelini, Giorgio Parisi, Tommaso Rizzo

Abstract: In statistical physics, one of the standard methods to study second order phase transitions is the renormalization group that usually leads to an expansion around the corresponding fully connected solution. Unfortunately, often in disordered models, some important finite dimensional second-order phase transitions are qualitatively different or absent in the corresponding fully connected model: in… ▽ More In statistical physics, one of the standard methods to study second order phase transitions is the renormalization group that usually leads to an expansion around the corresponding fully connected solution. Unfortunately, often in disordered models, some important finite dimensional second-order phase transitions are qualitatively different or absent in the corresponding fully connected model: in such cases the standard expansion fails. Recently, a new method, the $M$-layer one, has been introduced that performs an expansion around a different soluble mean field model: the Bethe lattice one. This new method has been already used to compute the upper critical dimension $D_U$ of different disordered systems such as the Random Field Ising model or the Spin glass model with field. If then one wants to go beyond and construct an expansion around $D_U$ to understand how critical quantities get renormalized, the actual computation of all the numerical factors is needed. This next step has still not been performed, being technically more involved. In this paper we perform this computation for the ferromagnetic Ising model without quenched disorder, in finite dimensions: we show that, at one-loop order inside the $M$-layer approach, we recover the continuum quartic field theory and we are able to identify the coupling constant $g$ and the other parameters of the theory, as a function of macroscopic and microscopic details of the model such as the lattice spacing, the physical lattice dimension and the temperature. This is a fundamental step that will help in applying in the future the same techniques to more complicated systems, for which the standard field theoretical approach is impracticable. △ Less

Submitted 22 March, 2024; v1 submitted 2 March, 2024; originally announced March 2024.

A network-constrain Weibull AFT model for biomarkers discovery

Authors: Claudia Angelini, Daniela De Canditiis, Italia De Feis, Antonella Iuliano

Abstract: We propose AFTNet, a novel network-constraint survival analysis method based on the Weibull accelerated failure time (AFT) model solved by a penalized likelihood approach for variable selection and estimation. When using the log-linear representation, the inference problem becomes a structured sparse regression problem for which we explicitly incorporate the correlation patterns among predictors u… ▽ More We propose AFTNet, a novel network-constraint survival analysis method based on the Weibull accelerated failure time (AFT) model solved by a penalized likelihood approach for variable selection and estimation. When using the log-linear representation, the inference problem becomes a structured sparse regression problem for which we explicitly incorporate the correlation patterns among predictors using a double penalty that promotes both sparsity and grouping effect. Moreover, we establish the theoretical consistency for the AFTNet estimator and present an efficient iterative computational algorithm based on the proximal gradient descent method. Finally, we evaluate AFTNet performance both on synthetic and real data examples. △ Less

Submitted 28 February, 2024; originally announced February 2024.

A Hierarchical Spatial Transformer for Massive Point Samples in Continuous Space

Authors: Wenchong He, Zhe Jiang, Tingsong Xiao, Zelin Xu, Shigang Chen, Ronald Fick, Miles Medina, Christine Angelini

Abstract: Transformers are widely used deep learning architectures. Existing transformers are mostly designed for sequences (texts or time series), images or videos, and graphs. This paper proposes a novel transformer model for massive (up to a million) point samples in continuous space. Such data are ubiquitous in environment sciences (e.g., sensor observations), numerical simulations (e.g., particle-laden… ▽ More Transformers are widely used deep learning architectures. Existing transformers are mostly designed for sequences (texts or time series), images or videos, and graphs. This paper proposes a novel transformer model for massive (up to a million) point samples in continuous space. Such data are ubiquitous in environment sciences (e.g., sensor observations), numerical simulations (e.g., particle-laden flow, astrophysics), and location-based services (e.g., POIs and trajectories). However, designing a transformer for massive spatial points is non-trivial due to several challenges, including implicit long-range and multi-scale dependency on irregular points in continuous space, a non-uniform point distribution, the potential high computational costs of calculating all-pair attention across massive points, and the risks of over-confident predictions due to varying point density. To address these challenges, we propose a new hierarchical spatial transformer model, which includes multi-resolution representation learning within a quad-tree hierarchy and efficient spatial attention via coarse approximation. We also design an uncertainty quantification branch to estimate prediction confidence related to input feature noise and point sparsity. We provide a theoretical analysis of computational time complexity and memory costs. Extensive experiments on both real-world and synthetic datasets show that our method outperforms multiple baselines in prediction accuracy and our model can scale up to one million points on one NVIDIA A100 GPU. The code is available at \url{https://github.com/spatialdatasciencegroup/HST}. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: Accepted in NeurIPS 2023

Stochastic Gradient Descent-like relaxation is equivalent to Metropolis dynamics in discrete optimization and inference problems

Authors: Maria Chiara Angelini, Angelo Giorgio Cavaliere, Raffaele Marino, Federico Ricci-Tersenghi

Abstract: Is Stochastic Gradient Descent (SGD) substantially different from Metropolis Monte Carlo dynamics? This is a fundamental question at the time of understanding the most used training algorithm in the field of Machine Learning, but it received no answer until now. Here we show that in discrete optimization and inference problems, the dynamics of an SGD-like algorithm resemble very closely that of Me… ▽ More Is Stochastic Gradient Descent (SGD) substantially different from Metropolis Monte Carlo dynamics? This is a fundamental question at the time of understanding the most used training algorithm in the field of Machine Learning, but it received no answer until now. Here we show that in discrete optimization and inference problems, the dynamics of an SGD-like algorithm resemble very closely that of Metropolis Monte Carlo with a properly chosen temperature, which depends on the mini-batch size. This quantitative matching holds both at equilibrium and in the out-of-equilibrium regime, despite the two algorithms having fundamental differences (e.g.\ SGD does not satisfy detailed balance). Such equivalence allows us to use results about performances and limits of Monte Carlo algorithms to optimize the mini-batch size in the SGD-like algorithm and make it efficient at recovering the signal in hard inference problems. △ Less

Submitted 30 May, 2024; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: 19 pages, 9 figures

Journal ref: Scientific Reports 14, 11638 (2024)

Variational Density Propagation Continual Learning

Authors: Christopher Angelini, Nidhal Bouaynaya, Ghulam Rasool

Abstract: Deep Neural Networks (DNNs) deployed to the real world are regularly subject to out-of-distribution (OoD) data, various types of noise, and shifting conceptual objectives. This paper proposes a framework for adapting to data distribution drift modeled by benchmark Continual Learning datasets. We develop and evaluate a method of Continual Learning that leverages uncertainty quantification from Baye… ▽ More Deep Neural Networks (DNNs) deployed to the real world are regularly subject to out-of-distribution (OoD) data, various types of noise, and shifting conceptual objectives. This paper proposes a framework for adapting to data distribution drift modeled by benchmark Continual Learning datasets. We develop and evaluate a method of Continual Learning that leverages uncertainty quantification from Bayesian Inference to mitigate catastrophic forgetting. We expand on previous approaches by removing the need for Monte Carlo sampling of the model weights to sample the predictive distribution. We optimize a closed-form Evidence Lower Bound (ELBO) objective approximating the predictive distribution by propagating the first two moments of a distribution, i.e. mean and covariance, through all network layers. Catastrophic forgetting is mitigated by using the closed-form ELBO to approximate the Minimum Description Length (MDL) Principle, inherently penalizing changes in the model likelihood by minimizing the KL Divergence between the variational posterior for the current task and the previous task's variational posterior acting as the prior. Leveraging the approximation of the MDL principle, we aim to initially learn a sparse variational posterior and then minimize additional model complexity learned for subsequent tasks. Our approach is evaluated for the task incremental learning scenario using density propagated versions of fully-connected and convolutional neural networks across multiple sequential benchmark datasets with varying task sequence lengths. Ultimately, this procedure produces a minimally complex network over a series of tasks mitigating catastrophic forgetting. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 6 pages, 13th Int'l Symposium on Image and Signal Processing and Analysis

Real-Space Renormalization group for spin glasses

Authors: Maria Chiara Angelini

Abstract: While in the fully-connected limit the solution of the spin-glass model is known, with the existence of a complex transition on a critical line in the temperature-external field phase diagram, in finite dimensions we don't know if a transition is present or not and, if present, if its nature is the same as the one in infinite dimensions. This work contains a review of the real-space Renormalizatio… ▽ More While in the fully-connected limit the solution of the spin-glass model is known, with the existence of a complex transition on a critical line in the temperature-external field phase diagram, in finite dimensions we don't know if a transition is present or not and, if present, if its nature is the same as the one in infinite dimensions. This work contains a review of the real-space Renormalization Group methods that have been applied to Spin Glasses, highlighting both their point of strength and weakness in characterizing the finite dimensional behaviour of the model. △ Less

Submitted 10 February, 2023; originally announced February 2023.

Comments: To appear as a contribution to the edited volume "Spin Glass Theory and Far Beyond - Replica Symmetry Breaking after 40 Years", P. Charbonneau, E. Marinari, M. Mezard, G. Parisi, F. Ricci-Tersenghi, G. Sicuro and F. Zamponi (World Scientific)

The Ising spin glass on random graphs at zero temperature: not all spins are glassy in the glassy phase

Authors: Gianmarco Perrupato, Maria Chiara Angelini, Giorgio Parisi, Federico Ricci-Tersenghi, Tommaso Rizzo

Abstract: We investigate the replica symmetry broken (RSB) phase of spin glass (SG) models in a random field defined on Bethe lattices at zero temperature. From the properties of the RSB solution we deduce a closed equation for the extreme values of the cavity fields. This equation turns out not to depend on the parameters defining the RSB, and it predicts that the spontaneous RSB does not take place homoge… ▽ More We investigate the replica symmetry broken (RSB) phase of spin glass (SG) models in a random field defined on Bethe lattices at zero temperature. From the properties of the RSB solution we deduce a closed equation for the extreme values of the cavity fields. This equation turns out not to depend on the parameters defining the RSB, and it predicts that the spontaneous RSB does not take place homogeneously on the whole system. Indeed, there exist spins having the same effective local field in all local ground states, exactly as in the replica symmetric (RS) phase, while the spontaneous RSB manifests only on the remaining spins, whose fraction vanishes at criticality. The characterization in terms of spins having fixed or fluctuating local fields can be extended also to the random field Ising model (RFIM), in which case the fluctuating spins are the only responsible for the spontaneous magnetization in the ferromagnetic phase. Close to criticality we are able to connect the statistics of the local fields acting on the spins in the RSB phase with the correlation functions measured in the paramagnetic phase. Identifying the two types of spins on given instances of SG and RFIM, we show that they participate very differently to avalanches produced by flipping a single spin. From the scaling of the number of spins inducing RSB effects close to the critical point and using the $M$-layer expansion we estimate the upper critical dimension $D_U \geq 8$ for SG. △ Less

Submitted 10 November, 2022; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: 8 pages, 6 figures

Journal ref: Phys. Rev. B 106, 174202 (2022)

Modern graph neural networks do worse than classical greedy algorithms in solving combinatorial optimization problems like maximum independent set

Authors: Maria Chiara Angelini, Federico Ricci-Tersenghi

Abstract: The recent work ``Combinatorial Optimization with Physics-Inspired Graph Neural Networks'' [Nat Mach Intell 4 (2022) 367] introduces a physics-inspired unsupervised Graph Neural Network (GNN) to solve combinatorial optimization problems on sparse graphs. To test the performances of these GNNs, the authors of the work show numerical results for two fundamental problems: maximum cut and maximum inde… ▽ More The recent work ``Combinatorial Optimization with Physics-Inspired Graph Neural Networks'' [Nat Mach Intell 4 (2022) 367] introduces a physics-inspired unsupervised Graph Neural Network (GNN) to solve combinatorial optimization problems on sparse graphs. To test the performances of these GNNs, the authors of the work show numerical results for two fundamental problems: maximum cut and maximum independent set (MIS). They conclude that "the graph neural network optimizer performs on par or outperforms existing solvers, with the ability to scale beyond the state of the art to problems with millions of variables." In this comment, we show that a simple greedy algorithm, running in almost linear time, can find solutions for the MIS problem of much better quality than the GNN. The greedy algorithm is faster by a factor of $10^4$ with respect to the GNN for problems with a million variables. We do not see any good reason for solving the MIS with these GNN, as well as for using a sledgehammer to crack nuts. In general, many claims of superiority of neural networks in solving combinatorial problems are at risk of being not solid enough, since we lack standard benchmarks based on really hard problems. We propose one of such hard benchmarks, and we hope to see future neural network optimizers tested on these problems before any claim of superiority is made. △ Less

Submitted 2 January, 2023; v1 submitted 27 June, 2022; originally announced June 2022.

Comments: Comment to "Combinatorial Optimization with Physics-Inspired Graph Neural Networks'' [Nat Mach Intell 4 (2022) 367] https://www.nature.com/articles/s42256-022-00468-6

Journal ref: Nature Machine Intelligence (2022)

Limits and performances of algorithms based on simulated annealing in solving sparse hard inference problems

Authors: Maria Chiara Angelini, Federico Ricci-Tersenghi

Abstract: The planted coloring problem is a prototypical inference problem for which thresholds for Bayes optimal algorithms, like Belief Propagation (BP), can be computed analytically. In this paper, we analyze the limits and performances of the Simulated Annealing (SA), a Monte Carlo-based algorithm that is more general and robust than BP, and thus of broader applicability. We show that SA is sub-optimal… ▽ More The planted coloring problem is a prototypical inference problem for which thresholds for Bayes optimal algorithms, like Belief Propagation (BP), can be computed analytically. In this paper, we analyze the limits and performances of the Simulated Annealing (SA), a Monte Carlo-based algorithm that is more general and robust than BP, and thus of broader applicability. We show that SA is sub-optimal in the recovery of the planted solution because it gets attracted by glassy states that, instead, do not influence the BP algorithm. At variance with previous conjectures, we propose an analytic estimation for the SA algorithmic threshold by comparing the spinodal point of the paramagnetic phase and the dynamical critical temperature. This is a fundamental connection between thermodynamical phase transitions and out of equilibrium behavior of Glauber dynamics. We also study an improved version of SA, called replicated SA (RSA), where several weakly coupled replicas are cooled down together. We show numerical evidence that the algorithmic threshold for the RSA coincides with the Bayes optimal one. Finally, we develop an approximated analytical theory explaining the optimal performances of RSA and predicting the location of the transition towards the planted solution in the limit of a very large number of replicas. Our results for RSA support the idea that mismatching the parameters in the prior with respect to those of the generative model may produce an algorithm that is optimal and very robust. △ Less

Submitted 28 June, 2023; v1 submitted 9 June, 2022; originally announced June 2022.

Journal ref: Physical Review X 13.2 (2023): 021011

Mismatching as a tool to enhance algorithmic performances of Monte Carlo methods for the planted clique model

Authors: Maria Chiara Angelini, Paolo Fachin, Simone de Feo

Abstract: Over-parametrization was a crucial ingredient for recent developments in inference and machine-learning fields. However a good theory explaining this success is still lacking. In this paper we study a very simple case of mismatched over-parametrized algorithm applied to one of the most studied inference problem: the planted clique problem. We analyze a Monte Carlo (MC) algorithm in the same class… ▽ More Over-parametrization was a crucial ingredient for recent developments in inference and machine-learning fields. However a good theory explaining this success is still lacking. In this paper we study a very simple case of mismatched over-parametrized algorithm applied to one of the most studied inference problem: the planted clique problem. We analyze a Monte Carlo (MC) algorithm in the same class of the famous Jerrum algorithm. We show how this MC algorithm is in general suboptimal for the recovery of the planted clique. We show however how to enhance its performances by adding a (mismatched) parameter: the temperature; we numerically find that this over-parametrized version of the algorithm can reach the supposed algorithmic threshold for the planted clique problem. △ Less

Submitted 19 November, 2021; v1 submitted 10 June, 2021; originally announced June 2021.

Journal ref: J. Stat. Mech. (2021) 113406

Unexpected upper critical dimension for spin glass models in a field predicted by the loop expansion around the Bethe solution at zero temperature

Authors: Maria Chiara Angelini, Carlo Lucibello, Giorgio Parisi, Gianmarco Perrupato, Federico Ricci-Tersenghi, Tommaso Rizzo

Abstract: The spin-glass transition in a field in finite dimension is analyzed directly at zero temperature using a perturbative loop expansion around the Bethe lattice solution. The loop expansion is generated by the $M$-layer construction whose first diagrams are evaluated numerically and analytically. The generalized Ginzburg criterion reveals that the upper critical dimension below which mean-field theo… ▽ More The spin-glass transition in a field in finite dimension is analyzed directly at zero temperature using a perturbative loop expansion around the Bethe lattice solution. The loop expansion is generated by the $M$-layer construction whose first diagrams are evaluated numerically and analytically. The generalized Ginzburg criterion reveals that the upper critical dimension below which mean-field theory fails is $D_U \le 8$, at variance with the classical result $D_U = 6$ yielded by finite-temperature replica field theory. Our expansion around the Bethe lattice has two crucial differences with respect to the classical one. The finite connectivity $z$ of the lattice is directly included from the beginning in the Bethe lattice, while in the classical computation the finite connectivity is obtained through an expansion in $1/z$. Moreover, if one is interested in the zero temperature ($T = 0$) transition, one can directly expand around the $T = 0$ Bethe transition. The expansion directly at $T = 0$ is not possible in the classical framework because the fully connected spin glass does not have a transition at $T = 0$, being in the broken phase for any value of the external field. △ Less

Submitted 16 March, 2022; v1 submitted 31 March, 2021; originally announced March 2021.

Comments: 14 pages, 4 figures

Journal ref: Phys. Rev. Lett. 128, 075702 (2022)

On the problem of modeling the boat wake climate; the Florida intracoastal waterway

Authors: C. Forlini, R. Qayyum, M. Malej, M. -A. Y. -H. Lam, F. Shi, C. Angelini, A. Sheremet

Abstract: The impact of boat traffic on the health of coastal ecosystems is a multi-scale process: from minutes (individual wakes) to days (tidal modulation of sediment transport), to seasons and years (traffic is seasonal). A considerable numerical effort, notwithstanding the value of a boat-by-boat numerical modeling approach, is questionable, because of the practical impossibility of specifying the exact… ▽ More The impact of boat traffic on the health of coastal ecosystems is a multi-scale process: from minutes (individual wakes) to days (tidal modulation of sediment transport), to seasons and years (traffic is seasonal). A considerable numerical effort, notwithstanding the value of a boat-by-boat numerical modeling approach, is questionable, because of the practical impossibility of specifying the exact type and navigation characteristics for every boat comprising the traffic at any given time. Here, we propose a statistical-mechanics description of the traffic using a joint probability density of the wake population in some characteristic parameter space. We attempt to answer two basic questions: (1) what is the relevant parameter space and (2) how should a numerical model be tested for a wake population? We describe the linear and nonlinear characteristics of wakes observed in the Florida Intracoastal Waters. Adopting provisionally a two-dimensional parameter space (depth- and length-based Froude numbers) we conduct numerical simulations using the open-source FUNWAVE-TVD Boussinesq model. The model performance is excellent for weakly-dispersive, completely specified wakes (e.g., the analytical linear wakes), and also for the range of Froude numbers observed in the field, or for large container ships generating relatively long waves. The model is challenged by the short waves generated by small, slow boats. However, simulations suggest that the problem is confined to the deeper water domain and linear evolution. Nonlinear wake shoaling, essential for modeling wake-induced sediment transport and wake impact on the environment, is described well. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Solving the fully-connected spherical $p$-spin model with the cavity method: equivalence with the replica results

Authors: Giacomo Gradenigo, Maria Chiara Angelini, Luca Leuzzi, Federico Ricci-Tersenghi

Abstract: The spherical $p$-spin is a fundamental model for glassy physics, thanks to its analytic solution achievable via the replica method. Unfortunately the replica method has some drawbacks: it is very hard to apply to diluted models and the assumptions beyond it are not immediately clear. Both drawbacks can be overcome by the use of the cavity method, which, however, needs to be applied with care to s… ▽ More The spherical $p$-spin is a fundamental model for glassy physics, thanks to its analytic solution achievable via the replica method. Unfortunately the replica method has some drawbacks: it is very hard to apply to diluted models and the assumptions beyond it are not immediately clear. Both drawbacks can be overcome by the use of the cavity method, which, however, needs to be applied with care to spherical models. Here we show how to write the cavity equations for spherical $p$-spin models on complete graphs, both in the Replica Symmetric (RS) ansatz (corresponding to Belief Propagation) and in the 1-step Replica Symmetry Breaking (1RSB) ansatz (corresponding to Survey Propagation). The cavity equations can be solved by a Gaussian (RS) and multivariate Gaussian (1RSB) ansatz for the distribution of the cavity fields. We compute the free energy in both ansatzes and check that the results are identical to the replica computation, predicting a phase transition to a 1RSB phase at low temperatures. The advantages of solving the model with the cavity method are many. The physical meaning of any ansatz for the cavity marginals is very clear. The cavity method works directly with the distribution of local quantities, which allows to generalize the method to dilute graphs. What we are presenting here is the first step towards the solution of the diluted version of the spherical $p$-spin model, which is a fundamental model in the theory of random lasers and interesting $per~se$ as an easier-to-simulate version of the classical fully-connected $p$-spin model. △ Less

Submitted 5 February, 2021; v1 submitted 3 July, 2020; originally announced July 2020.

Comments: 26 pages

Journal ref: J. Stat. Mech. 113302 (2020)

City Planning with Augmented Reality

Authors: Catherine Angelini, Adam S. Williams, Mathew Kress, Edgar Ramos Vieira, Newton D'Souza, Naphtali D. Rishe, Joseph Medina, Francisco R. Ortega

Abstract: We present an early study designed to analyze how city planning and the health of senior citizens can benefit from the use of augmented reality (AR) using Microsoft's HoloLens. We also explore whether AR and VR can be used to help city planners receive real-time feedback from citizens, such as the elderly, on virtual plans, allowing for informed decisions to be made before any construction begins. We present an early study designed to analyze how city planning and the health of senior citizens can benefit from the use of augmented reality (AR) using Microsoft's HoloLens. We also explore whether AR and VR can be used to help city planners receive real-time feedback from citizens, such as the elderly, on virtual plans, allowing for informed decisions to be made before any construction begins. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: This was accepted by Graphics Interface 2019 (GI '2019) HCI track as a poster paper

Comment on "Real-space renormalization-group methods for hierarchical spin glasses"

Authors: Maria Chiara Angelini, Giorgio Parisi, Federico Ricci-Tersenghi

Abstract: In the paper [Angelini M C, Parisi G, and Ricci-Tersenghi F, Ensemble renormalization group for disordered systems, Phys. Rev. B 87 134201 (2013)] we introduced a real-space renormalization group called Ensemble Renormalization Group (ERG) and we applied it to the Edwards-Anderson model, obtaining estimates for the critical exponents in good agreement with those from Monte Carlo simulations. Recen… ▽ More In the paper [Angelini M C, Parisi G, and Ricci-Tersenghi F, Ensemble renormalization group for disordered systems, Phys. Rev. B 87 134201 (2013)] we introduced a real-space renormalization group called Ensemble Renormalization Group (ERG) and we applied it to the Edwards-Anderson model, obtaining estimates for the critical exponents in good agreement with those from Monte Carlo simulations. Recently the paper [Castellana M, Real-space renormalization-group methods for hierarchical spin glasses, J. Phys. A: Math. Theor. 52 445002 (2019)] re-examined the ERG method from a different perspective, concluding that the previous results were wrong, and claiming that the ERG method predicts trivially wrong critical exponents. In this comment we explain why the conclusions reached by Castellana are wrong, as they are based on a misinterpretation of finite-size effects. We conclude that the ERG method remains a good RG method to obtain critical exponents in strongly disordered models (if properly used). △ Less

Submitted 6 November, 2019; originally announced November 2019.

Comments: Comment on arXiv:1910.04521

New loop expansion for the Random Magnetic Field Ising Ferromagnets at zero temperature

Authors: Maria Chiara Angelini, Carlo Lucibello, Giorgio Parisi, Federico Ricci-Tersenghi, Tommaso Rizzo

Abstract: We apply to the Random Field Ising Model at zero temperature (T= 0) the perturbative loop expansion around the Bethe solution. A comparison with the standard epsilon-expansion is made, highlighting the key differences that make the new expansion much more appropriate to correctly describe strongly disordered systems, especially those controlled by a T = 0 RG fixed point. This new loop expansion pr… ▽ More We apply to the Random Field Ising Model at zero temperature (T= 0) the perturbative loop expansion around the Bethe solution. A comparison with the standard epsilon-expansion is made, highlighting the key differences that make the new expansion much more appropriate to correctly describe strongly disordered systems, especially those controlled by a T = 0 RG fixed point. This new loop expansion produces an effective theory with cubic vertices. We compute the one-loop corrections due to cubic vertices, finding new terms that are absent in the epsilon-expansion. However, these new terms are subdominant with respect to the standard, supersymmetric ones, therefore dimensional reduction is still valid at this order of the loop expansion. △ Less

Submitted 11 June, 2019; originally announced June 2019.

Comments: 6 pages + Supporting Information. Submitted to PNAS

Journal ref: PNAS 117, 2268-2274 (2020)

Monte Carlo algorithms are very effective in finding the largest independent set in sparse random graphs

Authors: Maria Chiara Angelini, Federico Ricci-Tersenghi

Abstract: The effectiveness of stochastic algorithms based on Monte Carlo dynamics in solving hard optimization problems is mostly unknown. Beyond the basic statement that at a dynamical phase transition the ergodicity breaks and a Monte Carlo dynamics cannot sample correctly the probability distribution in times linear in the system size, there are almost no predictions nor intuitions on the behavior of th… ▽ More The effectiveness of stochastic algorithms based on Monte Carlo dynamics in solving hard optimization problems is mostly unknown. Beyond the basic statement that at a dynamical phase transition the ergodicity breaks and a Monte Carlo dynamics cannot sample correctly the probability distribution in times linear in the system size, there are almost no predictions nor intuitions on the behavior of this class of stochastic dynamics. The situation is particularly intricate because, when using a Monte Carlo based algorithm as an optimization algorithm, one is usually interested in the out of equilibrium behavior which is very hard to analyse. Here we focus on the use of Parallel Tempering in the search for the largest independent set in a sparse random graph, showing that it can find solutions well beyond the dynamical threshold. Comparison with state-of-the-art message passing algorithms reveals that parallel tempering is definitely the algorithm performing best, although a theory explaining its behavior is still lacking. △ Less

Submitted 3 April, 2019; originally announced April 2019.

Comments: 14 pages, 12 figures

Journal ref: Phys. Rev. E 100, 013302 (2019)

Parallel Tempering for the planted clique problem

Authors: Maria Chiara Angelini

Abstract: The theoretical information threshold for the planted clique problem is $2\log_2(N)$, however no polynomial algorithm is known to recover a planted clique of size $O(N^{1/2-ε})$, $ε>0$. In this paper we will apply a standard method for the analysis of disordered models, the Parallel-Tempering (PT) algorithm, to the clique problem, showing numerically that its time-scaling in the hard region is ind… ▽ More The theoretical information threshold for the planted clique problem is $2\log_2(N)$, however no polynomial algorithm is known to recover a planted clique of size $O(N^{1/2-ε})$, $ε>0$. In this paper we will apply a standard method for the analysis of disordered models, the Parallel-Tempering (PT) algorithm, to the clique problem, showing numerically that its time-scaling in the hard region is indeed polynomial for the analyzed sizes. We also apply PT to a different but connected model, the Sparse Planted Independent Set problem. In this situation thresholds should be sharper and finite size corrections should be less important. Also in this case PT shows a polynomial scaling in the hard region for the recovery. △ Less

Submitted 8 January, 2019; v1 submitted 16 February, 2018; originally announced February 2018.

Comments: 12 pages, 5 figures

Journal ref: J. Stat. Mech. (2018) 073404

One-loop topological expansion for spin glasses in the large connectivity limit

Authors: Maria Chiara Angelini, Giorgio Parisi, Federico Ricci-Tersenghi

Abstract: We apply for the first time a new one-loop topological expansion around the Bethe solution to the spin-glass model with field in the high connectivity limit, following the methodological scheme proposed in a recent work. The results are completely equivalent to the well known ones, found by standard field theoretical expansion around the fully connected model (Bray and Roberts 1980, and following… ▽ More We apply for the first time a new one-loop topological expansion around the Bethe solution to the spin-glass model with field in the high connectivity limit, following the methodological scheme proposed in a recent work. The results are completely equivalent to the well known ones, found by standard field theoretical expansion around the fully connected model (Bray and Roberts 1980, and following works). However this method has the advantage that the starting point is the original Hamiltonian of the model, with no need to define an associated field theory, nor to know the initial values of the couplings, and the computations have a clear and simple physical meaning. Moreover this new method can also be applied in the case of zero temperature, when the Bethe model has a transition in field, contrary to the fully connected model that is always in the spin glass phase. Sharing with finite dimensional model the finite connectivity properties, the Bethe lattice is clearly a better starting point for an expansion with respect to the fully connected model. The present work is a first step towards the generalization of this new expansion to more difficult and interesting cases as the zero-temperature limit, where the expansion could lead to different results with respect to the standard one. △ Less

Submitted 31 January, 2018; v1 submitted 31 August, 2017; originally announced August 2017.

Comments: 8 pages, 1 figure

Journal ref: Europhys. Lett. 121, 27001 (2018)

Loop expansion around the Bethe approximation through the $M$-layer construction

Authors: Ada Altieri, Maria Chiara Angelini, Carlo Lucibello, Giorgio Parisi, Federico Ricci-Tersenghi, Tommaso Rizzo

Abstract: For every physical model defined on a generic graph or factor graph, the Bethe $M$-layer construction allows building a different model for which the Bethe approximation is exact in the large $M$ limit and it coincides with the original model for $M=1$. The $1/M$ perturbative series is then expressed by a diagrammatic loop expansion in terms of so-called fat-diagrams. Our motivation is to study so… ▽ More For every physical model defined on a generic graph or factor graph, the Bethe $M$-layer construction allows building a different model for which the Bethe approximation is exact in the large $M$ limit and it coincides with the original model for $M=1$. The $1/M$ perturbative series is then expressed by a diagrammatic loop expansion in terms of so-called fat-diagrams. Our motivation is to study some important second-order phase transitions that do exist on the Bethe lattice but are either qualitatively different or absent in the corresponding fully connected case. In this case the standard approach based on a perturbative expansion around the naive mean field theory (essentially a fully connected model) fails. On physical grounds, we expect that when the construction is applied to a lattice in finite dimension there is a small region of the external parameters close to the Bethe critical point where strong deviations from mean-field behavior will be observed. In this region, the $1/M$ expansion for the corrections diverges and it can be the starting point for determining the correct non-mean-field critical exponents using renormalization group arguments. In the end, we will show that the critical series for the generic observable can be expressed as a sum of Feynman diagrams with the same numerical prefactors of field theories. However, the contribution of a given diagram is not evaluated associating Gaussian propagators to its lines as in field theories: one has to consider the graph as a portion of the original lattice, replacing the internal lines with appropriate one-dimensional chains, and attaching to the internal points the appropriate number of infinite-size Bethe trees to restore the correct local connectivity of the original model. △ Less

Submitted 19 October, 2023; v1 submitted 26 July, 2017; originally announced July 2017.

Journal ref: J. Stat. Mech. (2017) 113303

Real Space Migdal-Kadanoff Renormalisation of Glassy Systems: Recent Results and a Critical Assessment

Authors: Maria Chiara Angelini, Giulio Biroli

Abstract: In this manuscript, in honour of L. Kadanoff, we present recent progress obtained in the description of finite dimensional glassy systems thanks to the Migdal-Kadanoff renormalisation group (MK-RG). We provide a critical assessment of the method, in particular discuss its limitation in describing situations in which an infinite number of pure states might be present, and analyse the MK-RG flow in… ▽ More In this manuscript, in honour of L. Kadanoff, we present recent progress obtained in the description of finite dimensional glassy systems thanks to the Migdal-Kadanoff renormalisation group (MK-RG). We provide a critical assessment of the method, in particular discuss its limitation in describing situations in which an infinite number of pure states might be present, and analyse the MK-RG flow in the limit of infinite dimensions. MK-RG predicts that the spin-glass transition in a field and the glass transition are governed by zero-temperature fixed points of the renormalization group flow. This implies a typical energy scale that grows, approaching the transition, as a power of the correlation length, thus leading to enormously large time-scales as expected from experiments and simulations. These fixed points exist only in dimensions larger than $d_L>3$ but they nevertheless influence the RG flow below it, in particular in three dimensions. MK-RG thus predicts a similar behavior for spin-glasses in a field and models of glasses and relates it to the presence of avoided critical points. △ Less

Submitted 10 February, 2017; originally announced February 2017.

Comments: Journal of Statistical Physics, Special Issue: Dedicated to the Memory of Leo Kadanoff

Journal ref: Journal of Statistical Physics 167 (3-4), 476-498 (2017)

A Macroscopic Mathematical Model For Cell Migration Assays Using A Real-Time Cell Analysis

Authors: Ezio Di Costanzo, Vincenzo Ingangi, Claudia Angelini, Maria Francesca Carfora, Maria Vincenza Carriero, Roberto Natalini

Abstract: Experiments of cell migration and chemotaxis assays have been classically performed in the so-called Boyden Chambers. A recent technology, xCELLigence Real Time Cell Analysis, is now allowing to monitor the cell migration in real time. This technology measures impedance changes caused by the gradual increase of electrode surface occupation by cells during the course of time and provide a Cell Inde… ▽ More Experiments of cell migration and chemotaxis assays have been classically performed in the so-called Boyden Chambers. A recent technology, xCELLigence Real Time Cell Analysis, is now allowing to monitor the cell migration in real time. This technology measures impedance changes caused by the gradual increase of electrode surface occupation by cells during the course of time and provide a Cell Index which is proportional to cellular morphology, spreading, ruffling and adhesion quality as well as cell number. In this paper we propose a macroscopic mathematical model, based on \emph{advection-reaction-diffusion} partial differential equations, describing the cell migration assay using the real-time technology. We carried out numerical simulations to compare simulated model dynamics with data of observed biological experiments on three different cell lines and in two experimental settings: absence of chemotactic signals (basal migration) and presence of a chemoattractant. Overall we conclude that our minimal mathematical model is able to describe the phenomenon in the real time scale and numerical results show a good agreement with the experimental evidences. △ Less

Submitted 5 July, 2016; originally announced July 2016.

MSC Class: 92B05; 92C17; 35K99

Journal ref: PLoS ONE 11(9): e0162553, 2016

Real Space Renormalization Group Theory of Disordered Models of Glasses

Authors: Maria Chiara Angelini, Giulio Biroli

Abstract: We develop a real space renormalisation group analysis of disordered models of glasses, in particular of the spin models at the origin of the Random First Order Transition theory. We find three fixed points respectively associated to the liquid state, to the critical behavior and to the glass state. The latter two are zero-temperature ones; this provides a natural explanation of the growth of effe… ▽ More We develop a real space renormalisation group analysis of disordered models of glasses, in particular of the spin models at the origin of the Random First Order Transition theory. We find three fixed points respectively associated to the liquid state, to the critical behavior and to the glass state. The latter two are zero-temperature ones; this provides a natural explanation of the growth of effective activation energy scale and the concomitant huge increase of relaxation time approaching the glass transition. The lower critical dimension depends on the nature of the interacting degrees of freedom and is higher than three for all models. This does not prevent three dimensional systems from being glassy. Indeed, we find that their renormalisation group flow is affected by the fixed points existing in higher dimension and in consequence is non-trivial. Within our theoretical framework the glass transition results to be an avoided phase transition. △ Less

Submitted 13 April, 2016; originally announced April 2016.

Comments: 6 pages, 3 figures

Spectral Detection on Sparse Hypergraphs

Authors: Maria Chiara Angelini, Francesco Caltagirone, Florent Krzakala, Lenka Zdeborová

Abstract: We consider the problem of the assignment of nodes into communities from a set of hyperedges, where every hyperedge is a noisy observation of the community assignment of the adjacent nodes. We focus in particular on the sparse regime where the number of edges is of the same order as the number of vertices. We propose a spectral method based on a generalization of the non-backtracking Hashimoto mat… ▽ More We consider the problem of the assignment of nodes into communities from a set of hyperedges, where every hyperedge is a noisy observation of the community assignment of the adjacent nodes. We focus in particular on the sparse regime where the number of edges is of the same order as the number of vertices. We propose a spectral method based on a generalization of the non-backtracking Hashimoto matrix into hypergraphs. We analyze its performance on a planted generative model and compare it with other spectral methods and with Bayesian belief propagation (which was conjectured to be asymptotically optimal for this model). We conclude that the proposed spectral method detects communities whenever belief propagation does, while having the important advantages to be simpler, entirely nonparametric, and to be able to learn the rule according to which the hyperedges were generated without prior information. △ Less

Submitted 15 July, 2015; originally announced July 2015.

Comments: 8 pages, 5 figures

Journal ref: 2015 53rd Annual Allerton Conference on Communication, Control, and Computing, pages 66 - 73, IEEE

First Observation of CP Violation in B0->D(*)CP h0 Decays by a Combined Time-Dependent Analysis of BaBar and Belle Data

Authors: The BaBar, Belle Collaborations, :, A. Abdesselam, I. Adachi, A. Adametz, T. Adye, H. Ahmed, H. Aihara, S. Akar, M. S. Alam, J. Albert, S. Al Said, R. Andreassen, C. Angelini, F. Anulli, K. Arinstein, N. Arnaud, D. M. Asner, D. Aston, V. Aulchenko, T. Aushev, R. Ayad, V. Babu, I. Badhrees , et al. (450 additional authors not shown)

Abstract: We report a measurement of the time-dependent CP asymmetry of B0->D(*)CP h0 decays, where the light neutral hadron h0 is a pi0, eta or omega meson, and the neutral D meson is reconstructed in the CP eigenstates K+ K-, K0S pi0 or K0S omega. The measurement is performed combining the final data samples collected at the Y(4S) resonance by the BaBar and Belle experiments at the asymmetric-energy B fac… ▽ More We report a measurement of the time-dependent CP asymmetry of B0->D(*)CP h0 decays, where the light neutral hadron h0 is a pi0, eta or omega meson, and the neutral D meson is reconstructed in the CP eigenstates K+ K-, K0S pi0 or K0S omega. The measurement is performed combining the final data samples collected at the Y(4S) resonance by the BaBar and Belle experiments at the asymmetric-energy B factories PEP-II at SLAC and KEKB at KEK, respectively. The data samples contain ( 471 +/- 3 ) x 10^6 BB pairs recorded by the BaBar detector and ( 772 +/- 11 ) x 10^6, BB pairs recorded by the Belle detector. We measure the CP asymmetry parameters -eta_f S = +0.66 +/- 0.10 (stat.) +/- 0.06 (syst.) and C = -0.02 +/- 0.07 (stat.) +/- 0.03 (syst.). These results correspond to the first observation of CP violation in B0->D(*)CP h0 decays. The hypothesis of no mixing-induced CP violation is excluded in these decays at the level of 5.4 standard deviations. △ Less

Submitted 20 July, 2015; v1 submitted 15 May, 2015; originally announced May 2015.

Comments: 9 pages, 2 figures, submitted to Physical Review Letters

Spin Glass in a Field: a New Zero-Temperature Fixed Point in Finite Dimensions

Authors: Maria Chiara Angelini, Giulio Biroli

Abstract: By using real space renormalisation group (RG) methods we show that spin-glasses in a field display a new kind of transition in high dimensions. The corresponding critical properties and the spin-glass phase are governed by two non-perturbative zero temperature fixed points of the RG flow. We compute the critical exponents, discuss the RG flow and its relevance for three dimensional systems. The n… ▽ More By using real space renormalisation group (RG) methods we show that spin-glasses in a field display a new kind of transition in high dimensions. The corresponding critical properties and the spin-glass phase are governed by two non-perturbative zero temperature fixed points of the RG flow. We compute the critical exponents, discuss the RG flow and its relevance for three dimensional systems. The new spin-glass phase we discovered has unusual properties, which are intermediate between the ones conjectured by droplet and full replica symmetry breaking theories. These results provide a new perspective on the long-standing debate about the behaviour of spin-glasses in a field. △ Less

Submitted 16 July, 2015; v1 submitted 3 September, 2014; originally announced September 2014.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. Lett. 114, 095701 (2015)

The Super-Potts glass: a new disordered model for glass-forming liquids

Authors: Maria Chiara Angelini, Giulio Biroli

Abstract: We introduce a new disordered system, the Super-Potts model, which is a more frustrated version of the Potts glass. Its elementary degrees of freedom are variables that can take M values and are coupled via pair-wise interactions. Its exact solution on a completely connected lattice demonstrates that for large enough M it belongs to the class of mean-field systems solved by a one step replica symm… ▽ More We introduce a new disordered system, the Super-Potts model, which is a more frustrated version of the Potts glass. Its elementary degrees of freedom are variables that can take M values and are coupled via pair-wise interactions. Its exact solution on a completely connected lattice demonstrates that for large enough M it belongs to the class of mean-field systems solved by a one step replica symmetry breaking Ansatz. Numerical simulations by the parallel tempering technique show that in three dimensions it displays a phenomenological behaviour similar to the one of glass-forming liquids. The Super-Potts glass is therefore the first long-sought disordered model allowing one to perform extensive and detailed studies of the Random First Order Transition in finite dimensions. We also discuss its behaviour for small values of M, which is similar to the one of spin-glasses in a field. △ Less

Submitted 14 January, 2015; v1 submitted 28 July, 2014; originally announced July 2014.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. B 90, 220201(R) (2014)

The Physics of the B Factories

Authors: A. J. Bevan, B. Golob, Th. Mannel, S. Prell, B. D. Yabsley, K. Abe, H. Aihara, F. Anulli, N. Arnaud, T. Aushev, M. Beneke, J. Beringer, F. Bianchi, I. I. Bigi, M. Bona, N. Brambilla, J. B rodzicka, P. Chang, M. J. Charles, C. H. Cheng, H. -Y. Cheng, R. Chistov, P. Colangelo, J. P. Coleman, A. Drutskoy , et al. (2009 additional authors not shown)

Abstract: This work is on the Physics of the B Factories. Part A of this book contains a brief description of the SLAC and KEK B Factories as well as their detectors, BaBar and Belle, and data taking related issues. Part B discusses tools and methods used by the experiments in order to obtain results. The results themselves can be found in Part C. Please note that version 3 on the archive is the auxiliary… ▽ More This work is on the Physics of the B Factories. Part A of this book contains a brief description of the SLAC and KEK B Factories as well as their detectors, BaBar and Belle, and data taking related issues. Part B discusses tools and methods used by the experiments in order to obtain results. The results themselves can be found in Part C. Please note that version 3 on the archive is the auxiliary version of the Physics of the B Factories book. This uses the notation alpha, beta, gamma for the angles of the Unitarity Triangle. The nominal version uses the notation phi_1, phi_2 and phi_3. Please cite this work as Eur. Phys. J. C74 (2014) 3026. △ Less

Submitted 31 October, 2015; v1 submitted 24 June, 2014; originally announced June 2014.

Comments: 928 pages, version 3 (arXiv:1406.6311v3) corresponds to the alpha, beta, gamma version of the book, the other versions use the phi1, phi2, phi3 notation

Report number: SLAC-PUB-15968, KEK Preprint 2014-3

Journal ref: Eur. Phys. J. C74 (2014) 3026

Relations between Short Range and Long Range Ising models

Authors: Maria Chiara Angelini, Giorgio Parisi, Federico Ricci-Tersenghi

Abstract: We perform a numerical study of the long range (LR) ferromagnetic Ising model with power law decaying interactions ($J \propto r^{-d-σ}$) both on a one-dimensional chain ($d=1$) and on a square lattice ($d=2$). We use advanced cluster algorithms to avoid the critical slowing down. We first check the validity of the relation connecting the critical behavior of the LR model with parameters $(d,σ)$ t… ▽ More We perform a numerical study of the long range (LR) ferromagnetic Ising model with power law decaying interactions ($J \propto r^{-d-σ}$) both on a one-dimensional chain ($d=1$) and on a square lattice ($d=2$). We use advanced cluster algorithms to avoid the critical slowing down. We first check the validity of the relation connecting the critical behavior of the LR model with parameters $(d,σ)$ to that of a short range (SR) model in an equivalent dimension $D$. We then study the critical behavior of the $d=2$ LR model close to the lower critical $σ$, uncovering that the spatial correlation function decays with two different power laws: the effect of the subdominant power law is much stronger than finite size effects and actually makes the estimate of critical exponents very subtle. By including this subdominant power law, the numerical data are consistent with the standard renormalization group (RG) prediction by Sak, thus making not necessary (and unlikely, according to Occam's razor) the recent proposal by Picco of having a new set of RG fixed points, in addition to the mean-field one and the SR one. △ Less

Submitted 30 January, 2014; v1 submitted 27 January, 2014; originally announced January 2014.

Comments: 25 pages, 8 figures

Journal ref: Phys. Rev. E 89, 062120 (2014)

Evidence for the decay B0 --> omega omega and search for B0 --> omega phi

Authors: The BABAR Collaboration, J. P. Lees, V. Poireau, V. Tisserand, E. Grauges, A. Palano, G. Eigen, B. Stugu, D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, M. J. Lee, G. Lynch, H. Koch, T. Schroeder, C. Hearty, T. S. Mattison, J. A. McKenna, R. Y. So, A. Khan, V. E. Blinov, A. R. Buzykaev, V. P. Druzhinin, V. B. Golubev, E. A. Kravchenko , et al. (312 additional authors not shown)

Abstract: We describe searches for B meson decays to the charmless vector-vector final states omega omega and omega phi with 471 x 10^6 B Bbar pairs produced in e+ e- annihilation at sqrt(s) = 10.58 GeV using the BABAR detector at the PEP-II collider at the SLAC National Accelerator Laboratory. We measure the branching fraction B(B0 --> omega omega) = (1.2 +- 0.3 +0.3-0.2) x 10^-6, where the first uncertain… ▽ More We describe searches for B meson decays to the charmless vector-vector final states omega omega and omega phi with 471 x 10^6 B Bbar pairs produced in e+ e- annihilation at sqrt(s) = 10.58 GeV using the BABAR detector at the PEP-II collider at the SLAC National Accelerator Laboratory. We measure the branching fraction B(B0 --> omega omega) = (1.2 +- 0.3 +0.3-0.2) x 10^-6, where the first uncertainty is statistical and the second is systematic, corresponding to a significance of 4.4 standard deviations. We also determine the upper limit B(B0 --> omega phi) < 0.7 x 10^-6 at 90% confidence level. These measurements provide the first evidence for the decay B0 --> omega omega, and an improvement of the upper limit for the decay B0 --> omega phi. △ Less

Submitted 6 March, 2014; v1 submitted 29 November, 2013; originally announced December 2013.

Comments: 7 pages, 2 postscript figures

Report number: BABAR-PUB-13/018, SLAC-PUB-15846

Journal ref: Phys. Rev. D 89, 051101 (2014)

Study of the decay $\bar{B}^{0}\rightarrowΛ_{c}^{+}\bar{p}π^{+}π^{-}$ and its intermediate states

Authors: The Babar Collaboration, J. P. Lees, V. Poireau, V. Tisserand, E. Grauges, A. Palano, G. Eigen, B. Stugu, D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, G. Lynch, H. Koch, T. Schroeder, D. J. Asgeirsson, C. Hearty, T. S. Mattison, J. A. McKenna, R. Y. So, A. Khan, V. E. Blinov, A. R. Buzykaev, V. P. Druzhinin, V. B. Golubev, E. A. Kravchenko , et al. (330 additional authors not shown)

Abstract: We study the decay $\bar{B}^{0}\rightarrowΛ_{c}^{+}\bar{p}π^{+}π^{-}$, reconstructing the Λ_{c}^{+} baryon in the $p K^{-}π^{+}$ mode, using a data sample of $467\times 10^{6}$ $B\bar{B}$ pairs collected with the BaBar detector at the PEP-2 storage rings at SLAC. We measure branching fractions for decays with intermediate $Σ_{c}$ baryons to be… ▽ More We study the decay $\bar{B}^{0}\rightarrowΛ_{c}^{+}\bar{p}π^{+}π^{-}$, reconstructing the Λ_{c}^{+} baryon in the $p K^{-}π^{+}$ mode, using a data sample of $467\times 10^{6}$ $B\bar{B}$ pairs collected with the BaBar detector at the PEP-2 storage rings at SLAC. We measure branching fractions for decays with intermediate $Σ_{c}$ baryons to be ${\cal B}[\bar{B}^{0}\rightarrowΣ_{c}(2455)^{++}\bar{p}π^{-}]=(21.3 \pm 1.0 \pm 1.0 \pm 5.5) \times 10^{-5}$, ${\cal B}[\bar{B}^{0}\rightarrowΣ_{c}(2520)^{++}\bar{p}π^{-}]=(11.5\pm 1.0 \pm 0.5 \pm 3.0)\times 10^{-5}$, ${\cal B}[\bar{B}^{0}\rightarrowΣ_{c}(2455)^{0}\bar{p}π^{+}]=(9.1 \pm 0.7 \pm 0.4 \pm 2.4)\times10^{-5}$, and ${\cal B}[\bar{B}^{0}\rightarrowΣ_{c}(2520)^{0}\bar{p}π^{+}]= (2.2 \pm 0.7 \pm 0.1\pm 0.6) \times 10^{-5}$, where the uncertainties are statistical, systematic, and due to the uncertainty on the $Λ_{c}^{+}\rightarrow\proton\Kmπ^{+}$ branching fraction, respectively. For decays without $Σ_{c}(2455)$ or $Σ_{c}(2520)$ resonances, we measure ${\cal B}[\bar{B}^{0}\rightarrowΛ_{c}^{+}\bar{p}π^{+}π^{-}]_{\mathrm{non-Σ_{c}}}=(79 \pm 4 \pm 4 \pm 20)\times10^{-5}$. The total branching fraction is determined to be ${\cal B}[\bar{B}^{0}\rightarrowΛ_{c}^{+}\bar{p}π^{+}π^{-}]_{\mathrm{total}}=(123 \pm 5 \pm 7 \pm 32)\times10^{-5}$. We examine multibody mass combinations in the resonant three-particle $Σ_{c}\bar{p}π$ final states and in the four-particle $Λ_{c}^{+}\bar{p}π^{+}π^{-}$ final state, and observe different characteristics for the $\bar{p}π$ combination in neutral versus doubly-charged $Σ_{c}$ decays. △ Less

Submitted 28 March, 2013; v1 submitted 1 February, 2013; originally announced February 2013.

Comments: 20 pages, 30 figures

Report number: Babar-PUB-12/028, SLAC-PUB-15363

Search for direct CP-violation in singly-Cabibbo suppressed D+- --> K+ K- pi+- decays

Authors: BaBar Collaboration, J. P. Lees, V. Poireau, V. Tisserand, J. Garra Tico, E. Grauges, A. Palano, G. Eigen, B. Stugu, D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, G. Lynch, H. Koch, T. Schroeder, D. J. Asgeirsson, C. Hearty, T. S. Mattison, J. A. McKenna, R. Y. So, A. Khan, V. E. Blinov, A. R. Buzykaev, V. P. Druzhinin, V. B. Golubev , et al. (338 additional authors not shown)

Abstract: We report on a search for direct CP asymmetry in the singly Cabibbo-suppressed decay D+- --> K+ K- pi+- using a data sample of 476 fb-1 accumulated with the BaBar detector running at and just below the Y(4S) resonance. The CP-violating decay rate asymmetry A_CP is determined to be (0.35 +- 0.30 +- 0.15)%. Model-dependent and model-independent Dalitz plot analysis techniques are used to search for… ▽ More We report on a search for direct CP asymmetry in the singly Cabibbo-suppressed decay D+- --> K+ K- pi+- using a data sample of 476 fb-1 accumulated with the BaBar detector running at and just below the Y(4S) resonance. The CP-violating decay rate asymmetry A_CP is determined to be (0.35 +- 0.30 +- 0.15)%. Model-dependent and model-independent Dalitz plot analysis techniques are used to search for CP-violating asymmetries in the various intermediate states. △ Less

Submitted 21 February, 2013; v1 submitted 8 December, 2012; originally announced December 2012.

Comments: 14 pages, 6 eps figures, to be submitted to PRD

Report number: BABAR-PUB-12/014, SLAC-PUB-15077

Branching fraction and form-factor shape measurements of exclusive charmless semileptonic B decays, and determination of |V_{ub}|

Authors: BABAR Collaboration, J. P. Lees, V. Poireau, V. Tisserand, J. Garra Tico, E. Grauges, A. Palano, G. Eigen, B. Stugu, D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, G. Lynch, H. Koch, T. Schroeder, D. J. Asgeirsson, C. Hearty, T. S. Mattison, J. A. McKenna, R. Y. So, A. Khan, V. E. Blinov, A. R. Buzykaev, V. P. Druzhinin, V. B. Golubev , et al. (340 additional authors not shown)

Abstract: We report the results of a study of the exclusive charmless semileptonic decays, B^0 --> pi^- l^+ nu, B^+ --> pi^0 l^+ nu, B^+ --> omega l^+ nu, B^+ --> eta l^+ nu and B^+ --> eta^' l^+ nu, (l = e or mu) undertaken with approximately 462x10^6 B\bar{B} pairs collected at the Upsilon(4S) resonance with the BABAR detector. The analysis uses events in which the signal B decays are reconstructed with a… ▽ More We report the results of a study of the exclusive charmless semileptonic decays, B^0 --> pi^- l^+ nu, B^+ --> pi^0 l^+ nu, B^+ --> omega l^+ nu, B^+ --> eta l^+ nu and B^+ --> eta^' l^+ nu, (l = e or mu) undertaken with approximately 462x10^6 B\bar{B} pairs collected at the Upsilon(4S) resonance with the BABAR detector. The analysis uses events in which the signal B decays are reconstructed with a loose neutrino reconstruction technique. We obtain partial branching fractions in several bins of q^2, the square of the momentum transferred to the lepton-neutrino pair, for B^0 --> pi^- l^+ nu, B^+ --> pi^0 l^+ nu, B^+ --> omega l^+ nu and B^+ --> eta l^+ nu. From these distributions, we extract the form-factor shapes f_+(q^2) and the total branching fractions BF(B^0 --> pi^- l^+ nu) = (1.45 +/- 0.04_{stat} +/- 0.06_{syst})x10^-4 (combined pi^- and pi^0 decay channels assuming isospin symmetry), BF(B^+ --> omega l^+ nu) = (1.19 +/- 0.16_{stat} +/- 0.09_{syst})x10^-4 and BF(B^+ --> eta l^+ nu) = (0.38 +/- 0.05_{stat} +/- 0.05_{syst})x10^-4. We also measure BF(B^+ --> eta^' l^+ nu) = (0.24 +/- 0.08_{stat} +/- 0.03_{syst})x10^-4. We obtain values for the magnitude of the CKM matrix element V_{ub} by direct comparison with three different QCD calculations in restricted q^2 ranges of B --> pi l^+ nu decays. From a simultaneous fit to the experimental data over the full q^2 range and the FNAL/MILC lattice QCD predictions, we obtain |V_{ub}| = (3.25 +/- 0.31)x10^-3, where the error is the combined experimental and theoretical uncertainty. △ Less

Submitted 6 November, 2012; v1 submitted 6 August, 2012; originally announced August 2012.

Comments: 35 pages, 14 figures, submitted to PRD

Report number: BABAR-PUB12/015; SLAC-PUB-15208

Journal ref: Phys.Rev.D86:092004,2012

Observation of Time Reversal Violation in the B0 Meson System

Authors: The BABAR Collaboration, J. P. Lees, V. Poireau, V. Tisserand, J. Garra Tico, E. Grauges, A. Palanoab, G. Eigen, B. Stugu, D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, G. Lynch, H. Koch, T. Schroeder, D. J. Asgeirsson, C. Hearty, T. S. Mattison, J. A. McKenna, R. Y. So, A. Khan, V. E. Blinov, A. R. Buzykaev, V. P. Druzhinin, V. B. Golubev , et al. (343 additional authors not shown)

Abstract: Although CP violation in the B meson system has been well established by the B factories, there has been no direct observation of time reversal violation. The decays of entangled neutral B mesons into definite flavor states ($B^0$ or $\bar{B}^0$), and $J/ψK_S^0$ or $c\bar{c} K_S^0$ final states (referred to as $B_+$ or $B_-$), allow comparisons between the probabilities of four pairs of T-conjugat… ▽ More Although CP violation in the B meson system has been well established by the B factories, there has been no direct observation of time reversal violation. The decays of entangled neutral B mesons into definite flavor states ($B^0$ or $\bar{B}^0$), and $J/ψK_S^0$ or $c\bar{c} K_S^0$ final states (referred to as $B_+$ or $B_-$), allow comparisons between the probabilities of four pairs of T-conjugated transitions, for example, $\bar{B}^0 \rightarrow B_-$ and $B_- \rightarrow \bar{B}^0$, as a function of the time difference between the two B decays. Using 468 million $B\bar{B}$ pairs produced in $Υ(4S)$ decays collected by the BABAR detector at SLAC, we measure T-violating parameters in the time evolution of neutral B mesons, yielding $ΔS_T^+ = -1.37 \pm 0.14 (stat.) \pm 0.06 (syst.)$ and $ΔS_T^- = 1.17 \pm 0.18 (stat.) \pm 0.11 (syst.)$. These nonzero results represent the first direct observation of T violation through the exchange of initial and final states in transitions that can only be connected by a T-symmetry transformation. △ Less

Submitted 2 January, 2013; v1 submitted 24 July, 2012; originally announced July 2012.

Comments: 7 pages, 3 postscript figures, 1 table, appendix with 4 pages containing supplementary material. Accepted by Phys. Rev. Lett

Report number: BABAR-PUB-12/011; SLAC-PUB-15192

Journal ref: Phys. Rev. Lett. 109, 211801 (2012)

Compressed sensing with sparse, structured matrices

Authors: Maria Chiara Angelini, Federico Ricci-Tersenghi, Yoshiyuki Kabashima

Abstract: In the context of the compressed sensing problem, we propose a new ensemble of sparse random matrices which allow one (i) to acquire and compress a ρ0-sparse signal of length N in a time linear in N and (ii) to perfectly recover the original signal, compressed at a rate α, by using a message passing algorithm (Expectation Maximization Belief Propagation) that runs in a time linear in N. In the lar… ▽ More In the context of the compressed sensing problem, we propose a new ensemble of sparse random matrices which allow one (i) to acquire and compress a ρ0-sparse signal of length N in a time linear in N and (ii) to perfectly recover the original signal, compressed at a rate α, by using a message passing algorithm (Expectation Maximization Belief Propagation) that runs in a time linear in N. In the large N limit, the scheme proposed here closely approaches the theoretical bound ρ0 = α, and so it is both optimal and efficient (linear time complexity). More generally, we show that several ensembles of dense random matrices can be converted into ensembles of sparse random matrices, having the same thresholds, but much lower computational complexity. △ Less

Submitted 12 April, 2013; v1 submitted 12 July, 2012; originally announced July 2012.

Comments: 7 pages, 6 figures

Journal ref: Proc. Fiftieth Annual Allerton Conference, p. 808 (2012)

Initial-State Radiation Measurement of the e+e- -> pi+pi-pi+pi- Cross Section

Authors: J. P. Lees, V. Poireau, V. Tisserand, J. Garra Tico, E. Grauges, M. Martinelli, D. A. Milanes, A. Palano, M. Pappagallo, G. Eigen, B. Stugu, D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, G. Lynch, H. Koch, T. Schroeder, D. J. Asgeirsson, C. Hearty, T. S. Mattison, J. A. McKenna, A. Khan, V. E. Blinov, A. R. Buzykaev, V. P. Druzhinin , et al. (362 additional authors not shown)

Abstract: We study the process e+e- -> pi+pi-pi+pi-gamma, with a photon emitted from the initial-state electron or positron, using 454.3 fb^-1 of data collected with the BABAR detector at SLAC, corresponding to approximately 260,000 signal events. We use these data to extract the non-radiative sigma(e+e- ->pi+pi-pi+pi-) cross section in the energy range from 0.6 to 4.5 Gev. The total uncertainty of the cros… ▽ More We study the process e+e- -> pi+pi-pi+pi-gamma, with a photon emitted from the initial-state electron or positron, using 454.3 fb^-1 of data collected with the BABAR detector at SLAC, corresponding to approximately 260,000 signal events. We use these data to extract the non-radiative sigma(e+e- ->pi+pi-pi+pi-) cross section in the energy range from 0.6 to 4.5 Gev. The total uncertainty of the cross section measurement in the peak region is less than 3%, higher in precision than the corresponding results obtained from energy scan data. △ Less

Submitted 27 August, 2012; v1 submitted 26 January, 2012; originally announced January 2012.

Report number: SLAC-PUB-14857

Journal ref: PRD 85, 112009 (2012)

Ensemble renormalization group for disordered systems

Authors: Maria Chiara Angelini, Giorgio Parisi, Federico Ricci-Tersenghi

Abstract: We propose and study a renormalization group transformation that can be used also for models with strong quenched disorder, like spin glasses. The method is based on a mapping between disorder distributions, chosen such as to keep some physical properties (e.g., the ratio of correlations averaged over the ensemble) invariant under the transformation. We validate this ensemble renormalization group… ▽ More We propose and study a renormalization group transformation that can be used also for models with strong quenched disorder, like spin glasses. The method is based on a mapping between disorder distributions, chosen such as to keep some physical properties (e.g., the ratio of correlations averaged over the ensemble) invariant under the transformation. We validate this ensemble renormalization group by applying it to the hierarchical model (both the diluted ferromagnetic version and the spin glass version), finding results in agreement with Monte Carlo simulations. △ Less

Submitted 28 April, 2013; v1 submitted 29 November, 2011; originally announced November 2011.

Comments: 7 pages, 10 figures

Journal ref: Phy. Rev. B 87, 134201 (2013)

A Measurement of the Semileptonic Branching Fraction of the B_s Meson

Authors: The BaBar Collaboration, J. P. Lees, V. Poireau, V. Tisserand, J. Garra Tico, E. Grauges, M. Martinelli, D. A. Milanes, A. Palano, M. Pappagallo, G. Eigen, B. Stugu, D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, G. Lynch, H. Koch, T. Schroeder, D. J. Asgeirsson, C. Hearty, T. S. Mattison, J. A. McKenna, A. Khan, V. E. Blinov, A. R. Buzykaev , et al. (362 additional authors not shown)

Abstract: We report a measurement of the inclusive semileptonic branching fraction of the B_s meson using data collected with the BaBar detector in the center-of-mass (CM) energy region above the Upsilon(4S) resonance. We use the inclusive yield of phi mesons and the phi yield in association with a high-momentum lepton to perform a simultaneous measurement of the semileptonic branching fraction and the pr… ▽ More We report a measurement of the inclusive semileptonic branching fraction of the B_s meson using data collected with the BaBar detector in the center-of-mass (CM) energy region above the Upsilon(4S) resonance. We use the inclusive yield of phi mesons and the phi yield in association with a high-momentum lepton to perform a simultaneous measurement of the semileptonic branching fraction and the production rate of B_s mesons relative to all B mesons as a function of CM energy. The inclusive semileptonic branching fraction of the B_s meson is determined to be B(B_s to l nu X)=9.5 (+2.5/-2.0)(stat)(+1.1/-1.9)(syst)%, where l indicates the average of e and mu. △ Less

Submitted 10 January, 2012; v1 submitted 25 October, 2011; originally announced October 2011.

Comments: 8 pages, 7 eps figures

Report number: BaBar-Pub-11/021; SLAC-PUB-14653

Journal ref: Phys. Rev. D 85, 011101(R) (2012)

Search for hadronic decays of a light Higgs boson in the radiative decay Upsilon --> gamma A0

Authors: The BABAR Collaboration, J. P. Lees, V. Poireau, V. Tisserand, J. Garra Tico, E. Grauges, M. Martinelli, D. A. Milanes, A. Palano, M. Pappagallo, G. Eigen, B. Stugu, D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, G. Lynch, H. Koch, T. Schroeder, D. J. Asgeirsson, C. Hearty, T. S. Mattison, J. A. McKenna, R. Y. So, A. Khan, V. E. Blinov , et al. (364 additional authors not shown)

Abstract: We search for hadronic decays of a light Higgs boson (A0) produced in radiative decays of an Upsilon(2S) or Upsilon(3S) meson, Upsilon --> gamma A0. The data have been recorded by the BABAR experiment at the Upsilon(3S) and Upsilon(2S) center of mass energies, and include (121.3 \pm 1.2) x 10^6 Upsilon(3S) and (98.3 \pm 0.9) x 10^6 Upsilon(2S) mesons. No significant signal is observed. We set 90%… ▽ More We search for hadronic decays of a light Higgs boson (A0) produced in radiative decays of an Upsilon(2S) or Upsilon(3S) meson, Upsilon --> gamma A0. The data have been recorded by the BABAR experiment at the Upsilon(3S) and Upsilon(2S) center of mass energies, and include (121.3 \pm 1.2) x 10^6 Upsilon(3S) and (98.3 \pm 0.9) x 10^6 Upsilon(2S) mesons. No significant signal is observed. We set 90% confidence level upper limits on the product branching fractions B(Upsilon(nS)-->gamma A0) x B(A0-->hadrons) (n=2 or 3) that range from 1 x 10^{-6} for an A0 mass of 0.3 GeV/c^2 to 8 x 10^{-5} at 7 GeV/c^2. △ Less

Submitted 17 August, 2011; originally announced August 2011.

Comments: 7 pages, 5 postscript figures, submitted to Physical Review Letters

Report number: Report-no: BABAR-PUB-11/019, SLAC-PUB-14541

Observation of the baryonic B decay B0bar --> Lambda_c^+ anti-Lambda K-

Authors: The BABAR Collaboration, J. P. Lees, V. Poireau, V. Tisserand, J. Garra Tico, E. Grauges, M. Martinelli, D. A. Milanes, A. Palano, M. Pappagallo, G. Eigen, B. Stugu, L. Sun, D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, G. Lynch, H. Koch, T. Schroeder, D. J. Asgeirsson, C. Hearty, T. S. Mattison, J. A. McKenna, A. Khan, V. E. Blinov , et al. (363 additional authors not shown)

Abstract: We report the observation of the baryonic B decay B0bar --> Lambda_c^+ anti-Lambda K- with a significance larger than 7 standard deviations based on 471x10^6$ BBbar pairs collected with the BABAR detector at the PEP-II storage ring at SLAC. We measure the branching fraction for the decay B0bar --> Lambda_c^+ anti-Lambda K- to be (3.8 \pm 0.8_{stat} \pm 0.2_{sys} \pm 1.0_{Lambda_c^+})x10^{-5}. The… ▽ More We report the observation of the baryonic B decay B0bar --> Lambda_c^+ anti-Lambda K- with a significance larger than 7 standard deviations based on 471x10^6$ BBbar pairs collected with the BABAR detector at the PEP-II storage ring at SLAC. We measure the branching fraction for the decay B0bar --> Lambda_c^+ anti-Lambda K- to be (3.8 \pm 0.8_{stat} \pm 0.2_{sys} \pm 1.0_{Lambda_c^+})x10^{-5}. The uncertainties are statistical, systematic, and due to the uncertainty in the Lambda_c^+ branching fraction. We find that the Lambda_c^+ K^- invariant mass distribution shows an enhancement above 3.5 GeV/c^2. △ Less

Submitted 16 August, 2011; originally announced August 2011.

Comments: 7 pages, 7 figures, submitted to Phys. Rev. D (Rapid Communications)

Study of radiative bottomonium transitions using converted photons

Authors: The BABAR Collaboration, J. P. Lees, V. Poireau, E. Prencipe, V. Tisserand, J. GarraTico, E. Grauges, M. Martinelli, D. A. Milanes, A. Palano, M. Pappagallo, G. Eigen, B. Stugu, L. Sun, D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, G. Lynch, H. Koch, T. Schroeder, D. J. Asgeirsson, C. Hearty, T. S. Mattison, J. A. McKenna, A. Khan , et al. (367 additional authors not shown)

Abstract: We use 111+/-1 million Upsilon(3S) and 89+/-1 million Upsilon(2S) events recorded by the BaBar detector at the PEP-II B-factory at SLAC to perform a study of radiative transitions between bottomonium states using photons that have been converted to e+e- pairs by the detector material. We observe Upsilon(3S) -> gamma chi_b0,2(1P) decay, make precise measurements of the branching fractions for chi_b… ▽ More We use 111+/-1 million Upsilon(3S) and 89+/-1 million Upsilon(2S) events recorded by the BaBar detector at the PEP-II B-factory at SLAC to perform a study of radiative transitions between bottomonium states using photons that have been converted to e+e- pairs by the detector material. We observe Upsilon(3S) -> gamma chi_b0,2(1P) decay, make precise measurements of the branching fractions for chi_b1,2(1P,2P) -> gamma Upsilon(1S) and chi_b1,2(2P) -> gamma Upsilon(2S) decays, and search for radiative decay to the eta_b(1S) and eta_b(2S) states. △ Less

Submitted 21 December, 2011; v1 submitted 27 April, 2011; originally announced April 2011.

Comments: 19 pages, 8 figures, published by Phys. Rev. D

Report number: BABAR-PUB-11/002; SLAC-PUB-14446

Journal ref: Phys. Rev. D 84, 072002 (2011)

Evidence for the h_b(1P) meson in the decay Upsilon(3S) --> pi0 h_b(1P)

Authors: The BABAR Collaboration, J. P. Lees, V. Poireau, E. Prencipe, V. Tisserand, J. Garra Tico, E. Grauges, M. Martinelli, D. A. Milanes, A. Palano, M. Pappagallo, G. Eigen, B. Stugu, L. Sun, D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, G. Lynch, I. L. Osipenkov, H. Koch, T. Schroeder, D. J. Asgeirsson, C. Hearty, T. S. Mattison, J. A. McKenna , et al. (386 additional authors not shown)

Abstract: Using a sample of 122 million Upsilon(3S) events recorded with the BaBar detector at the PEP-II asymmetric-energy e+e- collider at SLAC, we search for the $h_b(1P)$ spin-singlet partner of the P-wave chi_{bJ}(1P) states in the sequential decay Upsilon(3S) --> pi0 h_b(1P), h_b(1P) --> gamma eta_b(1S). We observe an excess of events above background in the distribution of the recoil mass against the… ▽ More Using a sample of 122 million Upsilon(3S) events recorded with the BaBar detector at the PEP-II asymmetric-energy e+e- collider at SLAC, we search for the $h_b(1P)$ spin-singlet partner of the P-wave chi_{bJ}(1P) states in the sequential decay Upsilon(3S) --> pi0 h_b(1P), h_b(1P) --> gamma eta_b(1S). We observe an excess of events above background in the distribution of the recoil mass against the pi0 at mass 9902 +/- 4(stat.) +/- 2(syst.) MeV/c^2. The width of the observed signal is consistent with experimental resolution, and its significance is 3.1sigma, including systematic uncertainties. We obtain the value (4.3 +/- 1.1(stat.) +/- 0.9(syst.)) x 10^{-4} for the product branching fraction BF(Upsilon(3S)-->pi0 h_b) x BF(h_b-->gamma eta_b). △ Less

Submitted 17 October, 2011; v1 submitted 22 February, 2011; originally announced February 2011.

Comments: 8 pages, 4 postscript figures, submitted to Phys. Rev. D (Rapid Communications)

Report number: SLAC-PUB-14378

Entropic long range order in a 3D spin glass model

Authors: Maria Chiara Angelini, Federico Ricci-Tersenghi

Abstract: We uncover a new kind of entropic long range order in finite dimensional spin glasses. We study the link-diluted version of the Edwards-Anderson spin glass model with bimodal couplings (J=+/-1) on a 3D lattice. By using exact reduction algorithms, we prove that there exists a region of the phase diagram (at zero temperature and link density low enough), where spins are long range correlated, even… ▽ More We uncover a new kind of entropic long range order in finite dimensional spin glasses. We study the link-diluted version of the Edwards-Anderson spin glass model with bimodal couplings (J=+/-1) on a 3D lattice. By using exact reduction algorithms, we prove that there exists a region of the phase diagram (at zero temperature and link density low enough), where spins are long range correlated, even if the ground states energy stiffness is null. In other words, in this region twisting the boundary conditions cost no energy, but spins are long range correlated by means of pure entropic effects. △ Less

Submitted 4 January, 2011; v1 submitted 12 July, 2010; originally announced July 2010.

Comments: 15 pages, 6 figures. v3: added a phase diagram for ferromagnetically biased couplings

Journal ref: J. Stat. Mech. P02002 (2011)