Showing 1–50 of 90 results for author: Leuzzi, L

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

pastamarkers: astrophysical data visualization with pasta-like markers

Authors: PASTA Collaboration, N. Borghi, E. Ceccarelli, A. Della Croce, L. Leuzzi, L. Rosignoli, A. Traina

Abstract: We aim at facilitating the visualization of astrophysical data for several tasks, such as uncovering patterns, presenting results to the community, and facilitating the understanding of complex physical relationships to the public. We present pastamarkers, a customized Python package fully compatible with matplotlib, that contains unique pasta-shaped markers meant to enhance the visualization of a… ▽ More We aim at facilitating the visualization of astrophysical data for several tasks, such as uncovering patterns, presenting results to the community, and facilitating the understanding of complex physical relationships to the public. We present pastamarkers, a customized Python package fully compatible with matplotlib, that contains unique pasta-shaped markers meant to enhance the visualization of astrophysical data. We prove that using different pasta types as markers can improve the clarity of astrophysical plots by reproducing some of the most famous plots in the literature. △ Less

Submitted 29 March, 2024; originally announced March 2024.

Comments: 3 pages, 3 figures

Logarithmic critical slowing down in complex systems: from statics to dynamics

Authors: Luca Leuzzi, Tommaso Rizzo

Abstract: We consider second-order phase transitions in which the order parameter is a replicated overlap matrix. We focus on a tricritical point that occurs in a variety of mean-field models and that, more generically, describes higher order liquid-liquid or liquid-glass transitions. We show that the static replicated theory implies slowing down with a logarithmic decay in time. The dynamical equations tur… ▽ More We consider second-order phase transitions in which the order parameter is a replicated overlap matrix. We focus on a tricritical point that occurs in a variety of mean-field models and that, more generically, describes higher order liquid-liquid or liquid-glass transitions. We show that the static replicated theory implies slowing down with a logarithmic decay in time. The dynamical equations turn out to be those predicted by schematic Mode Coupling Theory for supercooled viscous liquids at a $A_3$ singularity, where the parameter exponent is $λらむだ=1$. We obtain a quantitative expression for the parameter $μみゅー$ of the logarithmic decay in terms of cumulants of the overlap, which are physically observable in experiments or numerical simulations. △ Less

Submitted 23 May, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: 22 pages, 2 figures

Journal ref: Phys. Rev. B 109, 174211(2024)

Euclid Preparation XXXIII. Characterization of convolutional neural networks for the identification of galaxy-galaxy strong lensing events

Authors: Euclid Collaboration, L. Leuzzi, M. Meneghetti, G. Angora, R. B. Metcalf, L. Moscardini, P. Rosati, P. Bergamini, F. Calura, B. Clément, R. Gavazzi, F. Gentile, M. Lochner, C. Grillo, G. Vernardos, N. Aghanim, A. Amara, L. Amendola, S. Andreon, N. Auricchio, S. Bardelli, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia , et al. (194 additional authors not shown)

Abstract: Forthcoming imaging surveys will potentially increase the number of known galaxy-scale strong lenses by several orders of magnitude. For this to happen, images of tens of millions of galaxies will have to be inspected to identify potential candidates. In this context, deep learning techniques are particularly suitable for the finding patterns in large data sets, and convolutional neural networks (… ▽ More Forthcoming imaging surveys will potentially increase the number of known galaxy-scale strong lenses by several orders of magnitude. For this to happen, images of tens of millions of galaxies will have to be inspected to identify potential candidates. In this context, deep learning techniques are particularly suitable for the finding patterns in large data sets, and convolutional neural networks (CNNs) in particular can efficiently process large volumes of images. We assess and compare the performance of three network architectures in the classification of strong lensing systems on the basis of their morphological characteristics. We train and test our models on different subsamples of a data set of forty thousand mock images, having characteristics similar to those expected in the wide survey planned with the ESA mission \Euclid, gradually including larger fractions of faint lenses. We also evaluate the importance of adding information about the colour difference between the lens and source galaxies by repeating the same training on single-band and multi-band images. Our models find samples of clear lenses with $\gtrsim 90\%$ precision and completeness, without significant differences in the performance of the three architectures. Nevertheless, when including lenses with fainter arcs in the training set, the three models' performance deteriorates with accuracy values of $\sim 0.87$ to $\sim 0.75$ depending on the model. Our analysis confirms the potential of the application of CNNs to the identification of galaxy-scale strong lenses. We suggest that specific training with separate classes of lenses might be needed for detecting the faint lenses since the addition of the colour information does not yield a significant improvement in the current analysis, with the accuracy ranging from $\sim 0.89$ to $\sim 0.78$ for the different models. △ Less

Submitted 26 January, 2024; v1 submitted 17 July, 2023; originally announced July 2023.

Comments: 23 pages,12 figures

Journal ref: Astronomy & Astrophysics,2024, 681, A68

Intensity pseudo-localized phase in the glassy random laser

Authors: Jacopo Niedda, Luca Leuzzi, Giacomo Gradenigo

Abstract: Evidence of an emergent pseudo-localized phase characterizing the low-temperature replica symmetry breaking phase of the complex disordered models for glassy light is provided in the mode-locked random laser model. A pseudo-localized phase corresponds to a state in which the intensity of light modes is neither equipartited among all modes nor really localized on few of them. Such a hybrid phase ha… ▽ More Evidence of an emergent pseudo-localized phase characterizing the low-temperature replica symmetry breaking phase of the complex disordered models for glassy light is provided in the mode-locked random laser model. A pseudo-localized phase corresponds to a state in which the intensity of light modes is neither equipartited among all modes nor really localized on few of them. Such a hybrid phase has been recently characterized in other models, such as the Discrete Non-Linear Schrödinger equation, just as a finite size effect, while in the low temperature phase of the glassy random laser it seems to be robust in the limit of large size. △ Less

Submitted 24 February, 2023; v1 submitted 9 December, 2022; originally announced December 2022.

Comments: 21 pages, 6 figures

Hopfield-like open channel retrieval for disordered optical media

Authors: Marco Leonetti, Luca Leuzzi, Giancarlo Ruocco

Abstract: The measurement of the optical Transmission Matrix (TM) enables to access "open channels": input patterns, specific to each scattering structure, capable to deliver very high transmission. Various approaches, based either on multiple interferometric measurements or on systematic random testing of incident wavefronts, enable to estimate the inputs required to excite these open channels. Here, we pr… ▽ More The measurement of the optical Transmission Matrix (TM) enables to access "open channels": input patterns, specific to each scattering structure, capable to deliver very high transmission. Various approaches, based either on multiple interferometric measurements or on systematic random testing of incident wavefronts, enable to estimate the inputs required to excite these open channels. Here, we provide for the first time an approach enabling the complete and reference-less retrieval of the open channels. It is based on the full mapping all the pairwise interference terms resulting from all the input modes couples. We show that these interference terms are organized into a bi-dyadic coupling matrix whose eigenvalues enables to access the open channel. A disordered optical system, is thus behaving exactly like an Hopfield neural network, where a specific input vector (an eigenvalue of the neurons' coupling matrix) enables to retrieve a specific memory pattern. The proposed Hopfield like open-channel-retrieval approach, enables to reach almost 100$\%$ of the theoretically expected value of the Intensity. Moreover employing a digital micromirror device to modulate light, we demonstrate high speed laser scanning at the back of a disordered medium. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Comments: 6 pages, 5 figures

LeMoN: Lens Modelling with Neural networks -- I. Automated modelling of strong gravitational lenses with Bayesian Neural Networks

Authors: Fabrizio Gentile, Crescenzo Tortora, Giovanni Covone, Léon V. E. Koopmans, Rui Li, Laura Leuzzi, Nicola R. Napolitano

Abstract: The unprecedented number of gravitational lenses expected from new-generation facilities such as the ESA Euclid telescope and the Vera Rubin Observatory makes it crucial to rethink our classical approach to lens-modelling. In this paper, we present LeMoN (Lens Modelling with Neural networks): a new machine-learning algorithm able to analyse hundreds of thousands of gravitational lenses in a reason… ▽ More The unprecedented number of gravitational lenses expected from new-generation facilities such as the ESA Euclid telescope and the Vera Rubin Observatory makes it crucial to rethink our classical approach to lens-modelling. In this paper, we present LeMoN (Lens Modelling with Neural networks): a new machine-learning algorithm able to analyse hundreds of thousands of gravitational lenses in a reasonable amount of time. The algorithm is based on a Bayesian Neural Network: a new generation of neural networks able to associate a reliable confidence interval to each predicted parameter. We train the algorithm to predict the three main parameters of the Singular Isothermal Ellipsoid model (the Einstein radius and the two components of the ellipticity) by employing two simulated datasets built to resemble the imaging capabilities of the Hubble Space Telescope and the forthcoming Euclid satellite. In this work, we assess the accuracy of the algorithm and the reliability of the estimated uncertainties by applying the network to several simulated datasets of 10.000 images each. We obtain accuracies comparable to previous studies present in the current literature and an average modelling time of just 0.5s per lens. Finally, we apply the LeMoN algorithm to a pilot dataset of real lenses observed with HST during the SLACS program, obtaining unbiased estimates of their SIE parameters. The code is publicly available on GitHub (https://github.com/fab-gentile/LeMoN). △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: Submitted to MNRAS; The code will be available on GitHub upon acceptance of the paper

Multi-mode fiber reservoir computing overcomes shallow neural networks classifiers

Authors: Daniele Ancora, Matteo Negri, Antonio Gianfrate, Dimitris Trypogeorgos, Lorenzo Dominici, Daniele Sanvitto, Federico Ricci-Tersenghi, Luca Leuzzi

Abstract: In the field of disordered photonics, a common objective is to characterize optically opaque materials for controlling light delivery or performing imaging. Among various complex devices, multi-mode optical fibers stand out as cost-effective and easy-to-handle tools, making them attractive for several tasks. In this context, we leverage the reservoir computing paradigm to recast these fibers into… ▽ More In the field of disordered photonics, a common objective is to characterize optically opaque materials for controlling light delivery or performing imaging. Among various complex devices, multi-mode optical fibers stand out as cost-effective and easy-to-handle tools, making them attractive for several tasks. In this context, we leverage the reservoir computing paradigm to recast these fibers into random hardware projectors, transforming an input dataset into a higher dimensional speckled image set. The goal of our study is to demonstrate that using such randomized data for classification by training a single logistic regression layer improves accuracy compared to training on direct raw images. Interestingly, we found that the classification accuracy achieved using the reservoir is also higher than that obtained with the standard transmission matrix model, a widely accepted tool for describing light transmission through disordered devices. We find that the reason for such improved performance could be due to the fact that the hardware classifier operates in a flatter region of the loss landscape when trained on fiber data, which aligns with the current theory of deep neural networks. These findings strongly suggest that multi-mode fibers possess robust generalization properties, positioning them as promising tools for optically-assisted neural networks. With this study, in fact, we want to contribute to advancing the knowledge and practical utilization of these versatile instruments, which may play a significant role in shaping the future of machine learning. △ Less

Submitted 26 May, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Universality class of the mode-locked glassy random laser

Authors: Jacopo Niedda, Giacomo Gradenigo, Luca Leuzzi, Giorgio Parisi

Abstract: By means of enhanced Monte Carlo numerical simulations parallelized on GPU's we study the critical properties of the spin-glass-like model for the mode-locked glassy random laser, a $4$-spin model with complex spins with a global spherical constraint and quenched random interactions. Using two different boundary conditions for the mode frequencies we identify the critical points and the critical i… ▽ More By means of enhanced Monte Carlo numerical simulations parallelized on GPU's we study the critical properties of the spin-glass-like model for the mode-locked glassy random laser, a $4$-spin model with complex spins with a global spherical constraint and quenched random interactions. Using two different boundary conditions for the mode frequencies we identify the critical points and the critical indices of the random lasing phase transition using , with finite size scaling techniques. The outcome of the scaling analysis is that the mode-locked random laser is in a mean-field universality class, though different from the mean-field class of the Random Energy Model and the glassy random laser in the narrow band approximation, that is, the fully connected version of the present model. The low temperature (high pumping) phase is finally characterized by means of the overlap distribution and evidence for the onset of replica symmetry breaking in the lasing regime is provided. △ Less

Submitted 24 February, 2023; v1 submitted 9 October, 2022; originally announced October 2022.

Comments: 17 pages, 11 figures

Journal ref: SciPost Phys. 14, 144 (2023)

Replica symmetry breaking in random lasers: experimental measurement of the overlap distribution

Authors: Claudio Conti, Neda Ghofraniha, Luca Leuzzi, Giancarlo Ruocco

Abstract: In this chapter we report on the measurements of the overlap distribution of the replica symmetry breaking solution in complex disordered systems. After a general introduction to the problem of the experimental validation of the Parisi order parameter, we focus on the systems where the measurement has been possible for the first time: random lasers. Starting from first principles of light-matter i… ▽ More In this chapter we report on the measurements of the overlap distribution of the replica symmetry breaking solution in complex disordered systems. After a general introduction to the problem of the experimental validation of the Parisi order parameter, we focus on the systems where the measurement has been possible for the first time: random lasers. Starting from first principles of light-matter interaction we sketch the main steps leading to the construction of the statistical mechanical model for the dynamics of light modes in a random laser, a spherical multi-p-spin model with complex spins. A new overlap is introduced, the intensity fluctuation overlap, whose probability distribution, under specific assumptions, is equivalent to the Parisi overlap distribution. The experimental protocol for measuring this overlap is based on the possibility of experimentally realizing real replicas. After a description of the first experiment on the random laser made of T5CO$_x$ grains we review and discuss various experiments measuring the overlap distribution, as well as the possible connection with Levy-like distribution of the intensity of the light modes around the laser threshold, the connection with turbulence in fiber lasers and the role of spatial etherogeneities of light modes in random media. △ Less

Submitted 12 September, 2022; v1 submitted 8 September, 2022; originally announced September 2022.

Comments: Contribution to the edited volume "Spin Glass Theory and Far Beyond - Replica Symmetry Breaking after 40 years", World Scientific. 30 pages, 16 figures

A generalized Hopfield model to store and retrieve mismatched memory patterns

Authors: Luca Leuzzi, Alberto Patti, Federico Ricci-Tersenghi

Abstract: We study a class of Hopfield models where the memories are represented by a mixture of Gaussian and binary variables and the neurons are Ising spins. We study the properties of this family of models as the relative weight of the two kinds of variables in the patterns varies. We quantitatively determine how the retrieval phase squeezes towards zero as the memory patterns contain a larger fraction o… ▽ More We study a class of Hopfield models where the memories are represented by a mixture of Gaussian and binary variables and the neurons are Ising spins. We study the properties of this family of models as the relative weight of the two kinds of variables in the patterns varies. We quantitatively determine how the retrieval phase squeezes towards zero as the memory patterns contain a larger fraction of mismatched variables. As the memory is purely Gaussian retrieval is lost for any positive storage capacity. It is shown that this comes about because of the spherical symmetry of the free energy in the Gaussian case. Introducing two different memory pattern overlaps between spin configurations and each contribution to the pattern from the two kinds of variables one can observe that the Gaussian parts of the patterns act as a noise, making retrieval more difficult. The basins of attraction of the states, the accuracy of the retrieval and the storage capacity are studied by means of Monte Carlo numerical simulations. We uncover that even in the limit where the network capacity shrinks to zero, the (few) retrieval states maintain a large basin of attraction and large overlaps with the mismatched patterns. So the network can be used for retrieval, but with a very small capacity. △ Less

Submitted 20 June, 2022; v1 submitted 9 April, 2022; originally announced April 2022.

Comments: 18 pages, 27 figures

Journal ref: J. Stat. Mech. (2022) 073301

Statistical mechanics of biomolecular condensates via cavity methods

Authors: N. Lauber, O. Tichacek, R. Bose, C. Flamm, L. Leuzzi, T-Y Dora tang, K. Ruiz-Mirazo, D. De Martino

Abstract: Physical mechanisms of phase separation in living systems can play key physiological roles and have recently been the focus of intensive studies. The strongly heterogeneous and disordered nature of such phenomena in the biological domain poses difficult modeling challenges that require going beyond mean field approaches based on postulating a free energy landscape. The alternative pathway we take… ▽ More Physical mechanisms of phase separation in living systems can play key physiological roles and have recently been the focus of intensive studies. The strongly heterogeneous and disordered nature of such phenomena in the biological domain poses difficult modeling challenges that require going beyond mean field approaches based on postulating a free energy landscape. The alternative pathway we take in this work is to tackle the full statistical mechanics problem of calculating the partition function in these systems, starting from microscopic interactions, by means of cavity methods. We illustrate the procedure first on the simple binary case, and we then apply it successfully to ternary systems, in which the naive mean field approximations are proved inadequate. We then demonstrate the agreement with lattice model simulations, to finally contrast our theory also with experiments of coacervate formation by associative de-mixing of nucleotides and poly-lysine in aqueous solution. In this way, different types of evidence are provided to support cavity methods as ideal tools for quantitative modeling of biomolecular condensation, giving an optimal balance between the accurate consideration of spatial aspects of the microscopic dynamics and the fast computational results rooted in their analytical tractability. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: 8 pages, 5 figures

Worldwide bilateral geopolitical interactions network inferred from national disciplinary profiles

Authors: Maria Grazia Izzo, Cinzia Daraio, Luca Leuzzi, Giammarco Quaglia, Giancarlo Ruocco

Abstract: A disciplinary profile of a country is defined as the versor whose components are the numbers of articles produced in a given discipline divided the overall production of the country. Starting from the Essential Science Indicators (ESI) schema of classification of subject area, we obtained the yearly disciplinary profiles of a worldwide graph, where on each node sits a country, in the two time int… ▽ More A disciplinary profile of a country is defined as the versor whose components are the numbers of articles produced in a given discipline divided the overall production of the country. Starting from the Essential Science Indicators (ESI) schema of classification of subject area, we obtained the yearly disciplinary profiles of a worldwide graph, where on each node sits a country, in the two time intervals [1980-1988] and [1992-2017], the fall of the Berlin Wall being the watershed. We analyse the empirical pairwise cross-correlation matrices of the time series of disciplinary profiles. The contrast with random matrix theory proves that, beyond measurement noise, the empirical cross-correlation matrices bring genuine information. Arising from the Shannon theorem as the least-structured model consistent with the measured pairwise correlations, the stationary probability distribution of disciplinary profiles can be described by a Boltzmann distribution related to a generalized $n_d$-dimensional Heisenberg model. The set of network interactions of the Heisenberg model have been inferred and to it they have been applied two clusterization methods, hierarchical clustering and principal component analysis. On a geopolitical plane this allow to obtain a characterization of the worldwide bilateral interactions based on physical modeling. A simple geopolitical analysis reveals the consistency of the results obtained and gives a boost to deeper historical analysis. In order to obtain the optimal set of pairwise interactions we used a Pseudo-Likelihood approach. We analytically computed the Pseudo-Likelihood and its gradient. The analytical computations deserve interest in whatever inference Bayesian problem involving a $n_d$-dimensional Heisenberg model. △ Less

Submitted 19 January, 2022; v1 submitted 3 January, 2022; originally announced January 2022.

Single-step transmission matrix retrieval for fast imaging through multi-mode fibers

Authors: Daniele Ancora, Lorenzo Dominici, Antonio Gianfrate, Paolo Cazzato, Milena De Giorgi, Dario Ballarini, Daniele Sanvitto, Luca Leuzzi

Abstract: Recovering the transmission matrix of a disordered medium is a challenging problem in disordered photonics. Usually, its reconstruction relies on a complex inversion that aims at connecting a fully-controlled input to the deterministic interference of the light field scrambled by the device. At the moment, iterative phase-retrieval protocols provide the fastest reconstructing frameworks, convergin… ▽ More Recovering the transmission matrix of a disordered medium is a challenging problem in disordered photonics. Usually, its reconstruction relies on a complex inversion that aims at connecting a fully-controlled input to the deterministic interference of the light field scrambled by the device. At the moment, iterative phase-retrieval protocols provide the fastest reconstructing frameworks, converging in a few tens of iterations. Exploiting the knowledge of speckle correlations, we construct a new phase retrieval algorithm that reduces the computational cost to a single iteration. Besides being faster, our method is effective also using less measurements than state-of-the-art protocols. Thanks to reducing computation time by one order of magnitude, our result can be groundbreaking for real-time optical operations in medical imaging. △ Less

Submitted 16 December, 2021; originally announced December 2021.

Comments: 9 pages, 7 figures

Journal ref: Photon. Res. 10, 2349-2358 (2022)

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)

Optical computation of a spin glass dynamics with tunable complexity

Authors: M. Leonetti, E. Hörmann, L. Leuzzi, G. Parisi, G. Ruocco

Abstract: Spin Glasses (SG) are paradigmatic models for physical, computer science, biological and social systems. The problem of studying the dynamics for SG models is NP hard, i.e., no algorithm solves it in polynomial time. Here we implement the optical simulation of a SG, exploiting the N segments of a wavefront shaping device to play the role of the spin variables, combining the interference at downstr… ▽ More Spin Glasses (SG) are paradigmatic models for physical, computer science, biological and social systems. The problem of studying the dynamics for SG models is NP hard, i.e., no algorithm solves it in polynomial time. Here we implement the optical simulation of a SG, exploiting the N segments of a wavefront shaping device to play the role of the spin variables, combining the interference at downstream of a scattering material to implement the random couplings between the spins (the J ij matrix) and measuring the light intensity on a number P of targets to retrieve the energy of the system. By implementing a plain Metropolis algorithm, we are able to simulate the spin model dynamics, while the degree of complexity of the potential energy landscape and the region of phase diagram explored is user-defined acting on the ratio the P/N = αあるふぁ. We study experimentally, numerically and analytically this peculiar system displaying a paramagnetic, a ferromagnetic and a SG phase, and we demonstrate that the transition temperature T g to the glassy phase from the paramagnetic phase grows with αあるふぁ. With respect to standard in silico approach, in the optical SG interaction terms are realized simultaneously when the independent light rays interferes at the target screen, enabling inherently parallel measurements of the energy, rather than computations scaling with N as in purely in silico simulations. △ Less

Submitted 1 July, 2020; v1 submitted 15 June, 2020; originally announced June 2020.

MSC Class: 78-05; 68W20

Journal ref: PNAS May 25, 2021 118 (21) e2015207118

Spin glasses in a field show a phase transition varying the distance among real replicas (and how to exploit it to find the critical line in a field)

Authors: Maddalena Dilucca, Luca Leuzzi, Giorgio Parisi, Federico Ricci-Tersenghi, Juan J. Ruiz-Lorenzo

Abstract: We discuss a phase transition in spin glass models which have been rarely considered in the past, namely the phase transition that may take place when two real replicas are forced to be at a larger distance (i.e. at a smaller overlap) than the typical one. In the first part of the work, by solving analytically the Sherrington-Kirkpatrick model in a field close to its critical point, we show that e… ▽ More We discuss a phase transition in spin glass models which have been rarely considered in the past, namely the phase transition that may take place when two real replicas are forced to be at a larger distance (i.e. at a smaller overlap) than the typical one. In the first part of the work, by solving analytically the Sherrington-Kirkpatrick model in a field close to its critical point, we show that even in a paramagnetic phase the forcing of two real replicas to an overlap small enough leads the model to a phase transition where the symmetry between replicas is spontaneously broken. More importantly, this phase transition is related to the de Almeida-Thouless (dAT) critical line. In the second part of the work, we exploit the phase transition in the overlap between two real replicas to identify the critical line in a field in finite-dimensional spin glasses. This is a notoriously difficult computational problem, because of huge finite-size corrections. We introduce a new method of analysis of Monte Carlo data for disordered systems, where the overlap between two real replicas is used as a conditioning variate. We apply this analysis to equilibrium measurements collected in the paramagnetic phase in a field, $h>0$ and $T_c(h)<T<T_c(h=0)$, of the $d=1$ spin glass model with long-range interactions decaying fast enough to be outside the regime of validity of the mean-field theory. We thus provide very reliable estimates for the thermodynamic critical temperature in a field. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Journal ref: Entropy 22, 250 (2020)

Self-starting nonlinear mode-locking in random lasers

Authors: Fabrizio Antenucci, Giovanni Lerario, Blanca Silva Fernandéz, Milena De Giorgi, Dario Ballarini, Daniele Sanvitto, Luca Leuzzi

Abstract: In ultra-fast multi-mode lasers, mode-locking is implemented by means of ad hoc devices, like saturable absorbers or modulators, allowing for very short pulses. This comes about because of nonlinear interactions induced among modes at different, well equispaced, frequencies. Theory predicts that the same locking of modes would occur in random lasers but, in absence of any device, its detection is… ▽ More In ultra-fast multi-mode lasers, mode-locking is implemented by means of ad hoc devices, like saturable absorbers or modulators, allowing for very short pulses. This comes about because of nonlinear interactions induced among modes at different, well equispaced, frequencies. Theory predicts that the same locking of modes would occur in random lasers but, in absence of any device, its detection is unfeasible so far. Because of the general interest in the phenomenology and understanding of random lasers and, moreover, because it is a first example of self-starting mode-locking we devise and test a way to measure such peculiar non-linear coupling. Through a detailed analysis of multi-mode correlations we provide clear evidence for the occurrence of nonlinear mode-coupling in the cavity-less random laser made of a powder of GaAs crystals and its self-starting mode-locking nature. The behavior of multi-point correlations among intensity peaks is tested against the nonlinear frequency matching condition equivalent to the one underlying phase-locking in ordered ultrafast lasers. Non-trivially large multi-point correlations are clearly observed for spatially overlapping resonances and turn out to sensitively depend on the frequency matching being satisfied, eventually demonstrating the occurrence of non-linear mode-locked mode-coupling. △ Less

Submitted 31 March, 2019; originally announced April 2019.

Comments: 13 pages, 8 figures, 1 table

Journal ref: Phys. Rev. Lett. 126, 173901 (2021)

Transmission Matrix Inference via Pseudolikelihood Decimation

Authors: Daniele Ancora, Luca Leuzzi

Abstract: One of the biggest challenges in the field of biomedical imaging is the comprehension and the exploitation of the photon scattering through disordered media. Many studies have pursued the solution to this puzzle, achieving light-focusing control or reconstructing images in complex media. In the present work, we investigate how statistical inference helps the calculation of the transmission matrix… ▽ More One of the biggest challenges in the field of biomedical imaging is the comprehension and the exploitation of the photon scattering through disordered media. Many studies have pursued the solution to this puzzle, achieving light-focusing control or reconstructing images in complex media. In the present work, we investigate how statistical inference helps the calculation of the transmission matrix in a complex scrambling environment, enabling its usage like a normal optical element. We convert a linear input-output transmission problem into a statistical formulation based on pseudolikelihood maximization, learning the coupling matrix via random sampling of intensity realizations. Our aim is to uncover insights from the scattering problem, encouraging the development of novel imaging techniques for better medical investigations, borrowing a number of statistical tools from spin-glass theory. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Glassiness and Lack of Equipartition in Random Lasers: the common roots of Ergodicity Breaking in Disordered and Non-linear Systems

Authors: Giacomo Gradenigo, Fabrizio Antenucci, Luca Leuzzi

Abstract: We present here for the first time a unifying perspective for the lack of equipartition in non-linear ordered systems and the low temperature phase-space fragmentation in disordered systems. We demonstrate that they are just two manifestation of the same underlying phenomenon: ergodicity breaking. Inspired by recent experiments, suggesting that lasing in optically active disordered media is relate… ▽ More We present here for the first time a unifying perspective for the lack of equipartition in non-linear ordered systems and the low temperature phase-space fragmentation in disordered systems. We demonstrate that they are just two manifestation of the same underlying phenomenon: ergodicity breaking. Inspired by recent experiments, suggesting that lasing in optically active disordered media is related to an ergodicity-breaking transition, we studied numerically a statistical mechanics model for the nonlinearly coupled light modes in a disordered medium under external pumping. Their collective behavior appears to be akin to the one displayed around the ergodicity-breaking transition in glasses, as we show measuring the glass order parameter of the replica-symmetry-breaking theory. Most remarkably, we also find that at the same critical point a breakdown of energy equipartition among light modes occurs, the typical signature of ergodicity breaking in non-linear systems as the celebrated Fermi-Pasta-Ulam model. The crucial ingredient of our system which allows us to find equipartition breakdown together with replica symmetry breaking is that the amplitudes of light modes are locally unbounded, i.e., they are only subject to a global constraint. The physics of random lasers appears thus as a unique test-bed to develop under a unifying perspective the study of ergodicity breaking in statistical disordered systems and non-linear ordered ones. △ Less

Submitted 5 February, 2021; v1 submitted 31 January, 2019; originally announced February 2019.

Comments: 14 pages, 4 figures

Journal ref: Phys. Rev. Research 2, 023399 (2020)

Learning Direct and Inverse Transmission Matrices

Authors: Daniele Ancora, Luca Leuzzi

Abstract: Linear problems appear in a variety of disciplines and their application for the transmission matrix recovery is one of the most stimulating challenges in biomedical imaging. Its knowledge turns any random media into an optical tool that can focus or transmit an image through disorder. Here, converting an input-output problem into a statistical mechanical formulation, we investigate how inference… ▽ More Linear problems appear in a variety of disciplines and their application for the transmission matrix recovery is one of the most stimulating challenges in biomedical imaging. Its knowledge turns any random media into an optical tool that can focus or transmit an image through disorder. Here, converting an input-output problem into a statistical mechanical formulation, we investigate how inference protocols can learn the transmission couplings by pseudolikelihood maximization. Bridging linear regression and thermodynamics let us propose an innovative framework to pursue the solution of the scattering-riddle. △ Less

Submitted 2 February, 2019; v1 submitted 15 January, 2019; originally announced January 2019.

Improved Pseudolikelihood Regularization and Decimation methods on Non-linearly Interacting Systems with Continuous Variables

Authors: Alessia Marruzzo, Payal Tyagi, Fabrizio Antenucci, Andrea Pagnani, Luca Leuzzi

Abstract: We propose and test improvements to state-of-the-art techniques of Bayeasian statistical inference based on pseudolikelihood maximization with $\ell_1$ regularization and with decimation. In particular, we present a method to determine the best value of the regularizer parameter starting from a hypothesis testing technique. Concerning the decimation, we also analyze the worst case scenario in whic… ▽ More We propose and test improvements to state-of-the-art techniques of Bayeasian statistical inference based on pseudolikelihood maximization with $\ell_1$ regularization and with decimation. In particular, we present a method to determine the best value of the regularizer parameter starting from a hypothesis testing technique. Concerning the decimation, we also analyze the worst case scenario in which there is no sharp peak in the tilded-pseudolikelihood function, firstly defined as a criterion to stop the decimation. Techniques are applied to noisy systems with non-linear dynamics, mapped onto multi-variable interacting Hamiltonian effective models for waves and phasors. Results are analyzed varying the number of available samples and the externally tunable temperature-like parameter mimicing real data noise. Eventually the behavior of inference procedures described are tested against a wrong hypothesis: non-linearly generated data are analyzed with a pairwise interacting hypothesis. Our analysis shows that, looking at the behavior of the inverse graphical problem as data size increases, the methods exposed allow to rule out a wrong hypothesis. △ Less

Submitted 23 May, 2018; v1 submitted 2 August, 2017; originally announced August 2017.

Comments: 37 pages, 24 figures

Journal ref: SciPost Phys. 5, 002 (2018)

Phase transitions in integer linear problems

Authors: S. Colabrese, D. De Martino, L. Leuzzi, E. Marinari

Abstract: The resolution of linear system with positive integer variables is a basic yet difficult computational problem with many applications. We consider sparse uncorrelated random systems parametrised by the density $c$ and the ratio $αあるふぁ=N/M$ between number of variables $N$ and number of constraints $M$. By means of ensemble calculations we show that the space of feasible solutions endows a Van-Der-Waals… ▽ More The resolution of linear system with positive integer variables is a basic yet difficult computational problem with many applications. We consider sparse uncorrelated random systems parametrised by the density $c$ and the ratio $αあるふぁ=N/M$ between number of variables $N$ and number of constraints $M$. By means of ensemble calculations we show that the space of feasible solutions endows a Van-Der-Waals phase diagram in the plane ($c$, $αあるふぁ$). We give numerical evidence that the associated computational problems become more difficult across the critical point and in particular in the coexistence region. △ Less

Submitted 17 May, 2017; originally announced May 2017.

Comments: 15 pages, 6 figures, comments are welcome

Inverse problem for multi-body interaction of nonlinear waves

Authors: Alessia Marruzzo, Payal Tyagi, Fabrizio Antenucci, Andrea Pagnani, Luca Leuzzi

Abstract: The inverse problem is studied in multi-body systems with nonlinear dynamics representing, e.g., phase-locked wave systems, standard multimode and random lasers. Using a general model for four-body interacting complex-valued variables we test two methods based on pseudolikelihood, respectively with regularization and with decimation, to determine the coupling constants from sets of measured config… ▽ More The inverse problem is studied in multi-body systems with nonlinear dynamics representing, e.g., phase-locked wave systems, standard multimode and random lasers. Using a general model for four-body interacting complex-valued variables we test two methods based on pseudolikelihood, respectively with regularization and with decimation, to determine the coupling constants from sets of measured configurations. We test statistical inference predictions for increasing number of sampled configurations and for an externally tunable {\em temperature}-like parameter mimicing real data noise and helping minimization procedures. Analyzed models with phasors and rotors are generalizations of problems of real-valued spherical problems (e.g., density fluctuations), discrete spins (Ising and vectorial Potts) or finite number of states (standard Potts): inference methods presented here can, then, be straightforward applied to a large class of inverse problems. △ Less

Submitted 1 January, 2017; v1 submitted 28 July, 2016; originally announced July 2016.

Comments: 9 pages, 9 figures

Journal ref: Scientific Reports 7, Article number: 3463 (2017)

Regularization and decimation pseudolikelihood approaches to statistical inference in $XY$-spin models

Authors: Payal Tyagi, Alessia Marruzzo, Andrea Pagnani, Fabrizio Antenucci, Luca Leuzzi

Abstract: We implement a pseudolikelyhood approach with l2-regularization as well as the recently introduced pseudolikelihood with decimation procedure to the inverse problem in continuous spin models on arbitrary networks, with arbitrarily disordered couplings. Performances of the approaches are tested against data produced by Monte Carlo numerical simulations and compared also from previously studied full… ▽ More We implement a pseudolikelyhood approach with l2-regularization as well as the recently introduced pseudolikelihood with decimation procedure to the inverse problem in continuous spin models on arbitrary networks, with arbitrarily disordered couplings. Performances of the approaches are tested against data produced by Monte Carlo numerical simulations and compared also from previously studied fully-connected mean-field-based inference techniques. The results clearly show that the best network reconstruction is obtained through the decimation scheme, that also allows to dwell the inference down to lower temperature regimes. Possible applications to phasor models for light propagation in random media are proposed and discussed. △ Less

Submitted 23 March, 2016; v1 submitted 16 March, 2016; originally announced March 2016.

Comments: 10 pages, 12 figures

Journal ref: Phys. Rev. B 94, 024203 (2016)

Multi-body quenched disordered $XY$ and $p$-clock models on random graphs

Authors: Alessia Marruzzo, Luca Leuzzi

Abstract: The $XY$ model with four-body quenched disordered interactions and its discrete $p$-clock proxy is studied on bipartite random graphs by means of the cavity method. The phase diagrams are determined from the ordered case to the spin-glass case. Dynamic, spinodal and thermodynamic transition lines are identified by analyzing free energy, complexity and tree reconstruction functions as temperature a… ▽ More The $XY$ model with four-body quenched disordered interactions and its discrete $p$-clock proxy is studied on bipartite random graphs by means of the cavity method. The phase diagrams are determined from the ordered case to the spin-glass case. Dynamic, spinodal and thermodynamic transition lines are identified by analyzing free energy, complexity and tree reconstruction functions as temperature and disorder are changed. The study of the convergence of the $p$-clock model to the $XY$ model is performed down to temperature low enough to determine all relevant transition points for different connectivity. △ Less

Submitted 23 March, 2016; v1 submitted 7 December, 2015; originally announced December 2015.

Comments: 15 pages, 6 figures

Journal ref: Phys. Rev. B 93, 094206 (2016)

Statistical mechanics models for multimode lasers and random lasers

Authors: F. Antenucci, A. Crisanti, M. Ibáñez Berganza, A. Marruzzo, L. Leuzzi

Abstract: We review recent statistical mechanical approaches to multimode laser theory. The theory has proved very effective to describe standard lasers. We refer of the mean field theory for passive mode locking and developments based on Monte Carlo simulations and cavity method to study the role of the frequency matching condition. The status for a complete theory of multimode lasing in open and disordere… ▽ More We review recent statistical mechanical approaches to multimode laser theory. The theory has proved very effective to describe standard lasers. We refer of the mean field theory for passive mode locking and developments based on Monte Carlo simulations and cavity method to study the role of the frequency matching condition. The status for a complete theory of multimode lasing in open and disordered cavities is discussed and the derivation of the general statistical models in this framework is presented. When light is propagating in a disordered medium, the system can be analyzed via the replica method. For high degrees of disorder and nonlinearity, a glassy behavior is expected at the lasing threshold, providing a suggestive link between glasses and photonics. We describe in details the results for the general Hamiltonian model in mean field approximation and mention an available test for replica symmetry breaking from intensity spectra measurements. Finally, we summary some perspectives still opened for such approaches. △ Less

Submitted 23 September, 2015; originally announced September 2015.

Comments: 16 pages, 7 figures

MSC Class: 78-06

Journal ref: Philosophical Magazine Volume 96, Issue 7-9, 2016

Dynamical arrest with zero complexity: the unusual behavior of the spherical Blume Emery Griffiths disordered model

Authors: Corrado Rainone, Ulisse Ferrari, Matteo Paoluzzi, Luca Leuzzi

Abstract: The short- and long-time dynamics of model systems undergoing a glass transition with apparent inversion of Kauzmann and dynamical arrest glass transition lines is investigated. These models belong to the class of the spherical mean-field approximation of a spin-$1$ model with $p$-body quenched disordered interaction, with $p>2$, termed spherical Blume-Emery-Griffiths models. Depending on temperat… ▽ More The short- and long-time dynamics of model systems undergoing a glass transition with apparent inversion of Kauzmann and dynamical arrest glass transition lines is investigated. These models belong to the class of the spherical mean-field approximation of a spin-$1$ model with $p$-body quenched disordered interaction, with $p>2$, termed spherical Blume-Emery-Griffiths models. Depending on temperature and chemical potential the system is found in a paramagnetic or in a glassy phase and the transition between these phases can be of a different nature. In specific regions of the phase diagram coexistence of low density and high density paramagnets can occur, as well as the coexistence of spin-glass and paramagnetic phases. The exact static solution for the glassy phase is known to be obtained by the one-step replica symmetry breaking ansatz. Different scenarios arise for both the dynamic and the thermodynamic transitions. These include: (i) the usual random first- order transition (Kauzmann-like) for mean-field glasses preceded by a dynamic transition, (ii) a thermodynamic first-order transition with phase coexistence and latent heat and (iii) a regime of apparent inversion of static transition line and dynamic transition lines, the latter defined as a non-zero complexity line. The latter inversion, though, turns out to be preceded by a novel dynamical arrest line at higher temperature. Crossover between different regimes is analyzed by solving mode coupling theory equations throughout the space of external thermodynamic parameters and the relationship with the underlying statics is discussed. △ Less

Submitted 17 June, 2015; originally announced June 2015.

Comments: 19 pages, 8 figures

Journal ref: Phys. Rev. E 92, 062150 (2015)

Inference for interacting linear waves in ordered and random media

Authors: P. Tyagi, A. Pagnani, F. Antenucci, M. Ibáñez Berganza, L. Leuzzi

Abstract: A statistical inference method is developed and tested for pairwise interacting systems whose degrees of freedom are continuous angular variables, such as planar spins in magnetic systems or wave phases in optics and acoustics. We investigate systems with both deterministic and quenched disordered couplings on two extreme topologies: complete and sparse graphs. To match further applications in opt… ▽ More A statistical inference method is developed and tested for pairwise interacting systems whose degrees of freedom are continuous angular variables, such as planar spins in magnetic systems or wave phases in optics and acoustics. We investigate systems with both deterministic and quenched disordered couplings on two extreme topologies: complete and sparse graphs. To match further applications in optics also complex couplings and external fields are considered and general inference formulas are derived for real and imaginary parts of Hermitian coupling matrices from real and imaginary parts of complex correlation functions. The whole procedure is, eventually, tested on numerically generated correlation functions and local magnetizations by means of Monte Carlo simulations. △ Less

Submitted 31 January, 2015; originally announced February 2015.

Comments: 17 pages, 7 figures

Journal ref: J. Stat. Mech. (2015) P05031

The glassy random laser: replica symmetry breaking in the intensity fluctuations of emission spectra

Authors: Fabrizio Antenucci, Andrea Crisanti, Luca Leuzzi

Abstract: The behavior of a newly introduced overlap parameter is analyzed, measuring the correlation between intensity fluctuations of waves in random media in different physical regimes, with varying amount of disorder and non-linearity. Its relationship is established to the standard Parisi overlap order parameter in replica theory for spin-glasses. In the complex spherical spin-glass model, describing t… ▽ More The behavior of a newly introduced overlap parameter is analyzed, measuring the correlation between intensity fluctuations of waves in random media in different physical regimes, with varying amount of disorder and non-linearity. Its relationship is established to the standard Parisi overlap order parameter in replica theory for spin-glasses. In the complex spherical spin-glass model, describing the onset and behavior of random lasers, replica symmetry breaking in the intensity fluctuation overlap is shown to occur at high pumping or low temperature. This order parameter identifies the laser transition in random media and describes its glassy nature in terms of emission spectra data, the only data so far accessible in random laser measurements. The theoretical analysis is, eventually, compared to recent intensity fluctuation overlap measurements demonstrating the validity of the theory and providing a straightforward interpretation of different spectral behaviors in different random lasers. △ Less

Submitted 11 June, 2015; v1 submitted 30 December, 2014; originally announced December 2014.

Comments: 11 pages, 3 figures

Journal ref: Sci. Rep. 5, 16792 (2015)

The Complex Spherical 2+4 Spin Glass: a Model for Nonlinear Optics in Random Media

Authors: Fabrizio Antenucci, Andrea Crisanti, Luca Leuzzi

Abstract: A disordered mean field model for multimode laser in open and irregular cavities is proposed and discussed within the replica analysis. The model includes the dynamics of the mode intensity and accounts also for the possible presence of a linear coupling between the modes, due, e.g., to the leakages from an open cavity. The complete phase diagram, in terms of disorder strength, source pumping and… ▽ More A disordered mean field model for multimode laser in open and irregular cavities is proposed and discussed within the replica analysis. The model includes the dynamics of the mode intensity and accounts also for the possible presence of a linear coupling between the modes, due, e.g., to the leakages from an open cavity. The complete phase diagram, in terms of disorder strength, source pumping and non-linearity, consists of four different optical regimes: incoherent fluorescence, standard mode locking, random lasing and the novel spontaneous phase locking. A replica symmetry breaking phase transition is predicted at the random lasing threshold. For a high enough strength of non-linearity, a whole region with nonvanishing complexity anticipates the transition, and the light modes in the disordered medium display typical discontinuous glassy behavior, i.e., the photonic glass has a multitude of metastable states that corresponds to different mode-locking processes in random lasers. The lasing regime is still present for very open cavities, though the transition becomes continuous at the lasing threshold. △ Less

Submitted 19 May, 2015; v1 submitted 30 December, 2014; originally announced December 2014.

Comments: 26 pages, 13 figures

Journal ref: Phys. Rev. A 91, 053816 (2015)

Statistical physics of nonlinear wave interaction

Authors: Fabrizio Antenucci, Miguel Ibáñez Berganza, Luca Leuzzi

Abstract: The thermodynamic properties of vector (O(2) and Complex Spherical) models with four-body interactions are analyzed. When defined in dense topologies, these are effective models for the nonlinear interaction of scalar fields in the presence of a stochastic noise, as has been well established for the case of the mode locking laser formation in a closed cavity. With the help of a novel efficient Mon… ▽ More The thermodynamic properties of vector (O(2) and Complex Spherical) models with four-body interactions are analyzed. When defined in dense topologies, these are effective models for the nonlinear interaction of scalar fields in the presence of a stochastic noise, as has been well established for the case of the mode locking laser formation in a closed cavity. With the help of a novel efficient Monte Carlo algorithm we show how beyond the fully connected case novel and rich phenomenology emerges. Below a certain dilution threshold, the spherical model condensates in a non-equipartite way, while in the XY model the transition becomes continuous and the O(2) symmetry remains unbroken, we attribute this fact to the invariance under local gauge transformations. The introduction of topological inhomogeneities in the network of quadruplets induces novel features: again symmetry conservation; the vanishing of two-point correlators; and a dynamical correlation function presenting two timescales, the large one being related to the transition between different degenerated configurations, connected by nonlocal gauge transformations. We discuss possible experimental implications of these results in the context of nonlinear optics. △ Less

Submitted 30 December, 2014; originally announced December 2014.

Comments: 21 pages, 24 figures

Journal ref: Phys. Rev. B 92, 014204 (2015)

Infinite volume extrapolation in the one-dimensional bond diluted Levy spin-glass model near its lower critical dimension

Authors: L. Leuzzi, G. Parisi, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo

Abstract: We revisited, by means of numerical simulations, the one dimensional bond diluted Levy Ising spin glasses outside the limit of validity of mean field theories. In these models the probability that two spins at distance $r$ interact (via a disordered interactions, $J_{ij}=\pm 1$) decays as $r^{-ρろー}$. We have estimated, using finite size scaling techniques, the infinite volume correlation length and… ▽ More We revisited, by means of numerical simulations, the one dimensional bond diluted Levy Ising spin glasses outside the limit of validity of mean field theories. In these models the probability that two spins at distance $r$ interact (via a disordered interactions, $J_{ij}=\pm 1$) decays as $r^{-ρろー}$. We have estimated, using finite size scaling techniques, the infinite volume correlation length and spin glass susceptibility for $ρろー=5/3$ and $ρろー=9/5$. We have obtained strong evidence for divergences of the previous observables at a non zero critical temperature. We discuss the behavior of the critical exponents, especially when approaching the value $ρろー=2$, corresponding to a critical threshold beyond which the model has no phase transition. Finally, we numerically study the model right at the threshold value $ρろー=2$. △ Less

Submitted 12 December, 2014; originally announced December 2014.

Comments: 9 pages and 16 figures

Journal ref: Phys. Rev. B 91, 064202 (2015)

Nonlinear $XY$ and $p$-clock models on sparse random graphs: mode-locking transition of localized waves

Authors: Alessia Marruzzo, Luca Leuzzi

Abstract: A statistical mechanic study of the $XY$ model with nonlinear interaction is presented on bipartite sparse random graphs. The model properties are compared to those of the $p$-clock model, in which the planar continuous spins are discretized into $p$ values. We test the goodness of the discrete approximation to the XY spins to be used in numerical computations and simulations and its limits of con… ▽ More A statistical mechanic study of the $XY$ model with nonlinear interaction is presented on bipartite sparse random graphs. The model properties are compared to those of the $p$-clock model, in which the planar continuous spins are discretized into $p$ values. We test the goodness of the discrete approximation to the XY spins to be used in numerical computations and simulations and its limits of convergence in given, $p$-dependent, temperature regimes. The models are applied to describe the mode-locking transition of the phases of light-modes in lasers at the critical lasing threshold. A frequency is assigned to each variable node and function nodes implement a frequency matching condition. A non-trivial unmagnetized phase-locking occurs at the phase transition, where the frequency dependence of the phases turns out to be linear in a broad range of frequencies, as in standard mode-locking multimode laser at the optical power threshold. △ Less

Submitted 17 November, 2014; originally announced November 2014.

Comments: 13 pages, 13 figures

Journal ref: Phys. Rev. B 91, 054201 (2015)

Statistical physical theory of mode-locking laser generation with a frequency comb

Authors: Fabrizio Antenucci, M. Ibáñez Berganza, Luca Leuzzi

Abstract: A study of the Mode-locking lasing pulse formation in closed cavities is presented within a statistical mechanical framework where the onset of laser coincides with a thermodynamic phase transition driven by the optical power pumped into the system. Electromagnetic modes are represented by classical degrees of freedom of a Hamiltonian model at equilibrium in an effective ensemble corresponding to… ▽ More A study of the Mode-locking lasing pulse formation in closed cavities is presented within a statistical mechanical framework where the onset of laser coincides with a thermodynamic phase transition driven by the optical power pumped into the system. Electromagnetic modes are represented by classical degrees of freedom of a Hamiltonian model at equilibrium in an effective ensemble corresponding to the stationary laser regime. By means of optimized Monte Carlo numerical simulations, the system properties are analyzed varying mode interaction dilution, gain profile and number of modes. Novel properties of the resulting mode-locking laser phase are presented, not observable by previous mean-field approaches. For strong dilution of the nonlinear interaction network, power condensation occurs as the whole optical intensity is taken by a few electromagnetic modes, whose number does not depend on the size of the system. For all reported cases laser thresholds, intensity spectra, and ultra-fast electromagnetic pulses are computed. △ Less

Submitted 23 April, 2015; v1 submitted 22 September, 2014; originally announced September 2014.

Comments: 7 pages, 6 figures + 2 videos in attachment. To view the videos download them from the ancillary files list or download and extract the gzipped tar source file listed under "Other formats" to retrieve the supplemental file

Journal ref: Phys. Rev. A 91, 043811 (2015)

Experimental evidence of replica symmetry breaking in random lasers

Authors: N. Ghofraniha, I. Viola, F. Di Maria, G. Barbarella, G. Gigli, L. Leuzzi, C. Conti

Abstract: Spin-glass theory is one of the leading paradigms of complex physics and describes condensed matter, neural networks and biological systems, ultracold atoms, random photonics, and many other research fields. According to this theory, identical systems under identical conditions may reach different states and provide different values for observable quantities. This effect is known as Replica Symmet… ▽ More Spin-glass theory is one of the leading paradigms of complex physics and describes condensed matter, neural networks and biological systems, ultracold atoms, random photonics, and many other research fields. According to this theory, identical systems under identical conditions may reach different states and provide different values for observable quantities. This effect is known as Replica Symmetry Breaking and is revealed by the shape of the probability distribution function of an order parameter named the Parisi overlap. However, a direct experimental evidence in any field of research is still missing. Here we investigate pulse-to-pulse fluctuations in random lasers, we introduce and measure the analogue of the Parisi overlap in independent experimental realizations of the same disordered sample, and we find that the distribution function yields evidence of a transition to a glassy light phase compatible with a replica symmetry breaking. △ Less

Submitted 16 February, 2015; v1 submitted 21 July, 2014; originally announced July 2014.

Comments: 10 pages, 5 figures

Journal ref: Nature Communications 6, 6058 (2015); CORRIGENDUM: Nature Communications 6, 6300 (2015)

General phase-diagram of multimodal ordered and disordered lasers in closed and open cavities

Authors: Fabrizio Antenucci, Claudio Conti, Andrea Crisanti, Luca Leuzzi

Abstract: We present a unified approach to the theory of multimodal laser cavities including a variable amount of structural disorder. A general mean-field theory is studied for waves in media with variable non-linearity and randomness. Phase diagrams are reported in terms of optical power, degree of disorder and degree of non-linearity, tuning between closed and open cavity scenario's. In the thermodynamic… ▽ More We present a unified approach to the theory of multimodal laser cavities including a variable amount of structural disorder. A general mean-field theory is studied for waves in media with variable non-linearity and randomness. Phase diagrams are reported in terms of optical power, degree of disorder and degree of non-linearity, tuning between closed and open cavity scenario's. In the thermodynamic limit of infinitely many modes the theory predicts four distinct regimes: a continuous wave behavior for low power, a standard mode-locking laser regime for high power and weak disorder, a random laser for high pumped power and large disorder and an intermediate regime of phase locking occurring in presence of disorder below the lasing threshold. △ Less

Submitted 30 June, 2014; originally announced June 2014.

Comments: 9 pages, 3 figures

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

A simple spin model for three steps relaxation and secondary proccesses in glass formers

Authors: A. Crisanti, L. Leuzzi

Abstract: A number of general trends are known to occur in systems displaying secondary processes in glasses and glass formers. Universal features can be identified as components of large and small cooperativeness whose competition leads to excess wings or apart peaks in the susceptibility spectrum. To the aim of understanding such rich and complex phenomenology we analyze the behavior of a model combining… ▽ More A number of general trends are known to occur in systems displaying secondary processes in glasses and glass formers. Universal features can be identified as components of large and small cooperativeness whose competition leads to excess wings or apart peaks in the susceptibility spectrum. To the aim of understanding such rich and complex phenomenology we analyze the behavior of a model combining two apart glassy components with a tunable different cooperativeness. The model salient feature is, indeed, based on the competition of the energetic contribution of groups of dynamically relevant variables, e.g., density fluctuations, interacting in small and large sets. We investigate how the model is able to reproduce the secondary processes physics without further ad hoc ingredients, displaying known trends and properties under cooling or pressing. △ Less

Submitted 7 July, 2014; v1 submitted 27 May, 2014; originally announced May 2014.

Comments: 11 Pages, 11 Figures

Journal ref: Journal of Non-Crystalline Solids 407, 110-117 (2015)

Finite size correction to the spectrum of regular random graphs: an analytical solution

Authors: Fernando L. Metz, Giorgio Parisi, Luca Leuzzi

Abstract: We develop a thorough analytical study of the $O(1/N)$ correction to the spectrum of regular random graphs with $N \rightarrow \infty$ nodes. The finite size fluctuations of the resolvent are given in terms of a weighted series over the contributions coming from loops of all possible lengths, from which we obtain the isolated eigenvalue as well as an analytical expression for the $O(1/N)$ correcti… ▽ More We develop a thorough analytical study of the $O(1/N)$ correction to the spectrum of regular random graphs with $N \rightarrow \infty$ nodes. The finite size fluctuations of the resolvent are given in terms of a weighted series over the contributions coming from loops of all possible lengths, from which we obtain the isolated eigenvalue as well as an analytical expression for the $O(1/N)$ correction to the continuous part of the spectrum. The comparison between this analytical formula and direct diagonalization results exhibits an excellent agreement, confirming the correctness of our expression. △ Less

Submitted 17 November, 2014; v1 submitted 11 March, 2014; originally announced March 2014.

Comments: Extended version with an extra appendix explaining the connection with rigorous results

Journal ref: Phys. Rev. E 90, 052109 (2014)

Small clusters Renormalization Group in 2D and 3D Ising and BEG models with ferro, antiferro and quenched disordered magnetic interactions

Authors: Fabrizio Antenucci, Andrea Crisanti, Luca Leuzzi

Abstract: The Ising and BEG models critical behavior is analyzed in 2D and 3D by means of a renormalization group scheme on small clusters made of a few lattice cells. Different kinds of cells are proposed for both ordered and disordered model cases. In particular, cells preserving a possible antiferromagnetic ordering under decimation allow for the determination of the Néel critical point and its scaling i… ▽ More The Ising and BEG models critical behavior is analyzed in 2D and 3D by means of a renormalization group scheme on small clusters made of a few lattice cells. Different kinds of cells are proposed for both ordered and disordered model cases. In particular, cells preserving a possible antiferromagnetic ordering under decimation allow for the determination of the Néel critical point and its scaling indices. These also provide more reliable estimates of the Curie fixed point than those obtained using cells preserving only the ferromagnetic ordering. In all studied dimensions, the present procedure does not yield the strong disorder critical point corresponding to the transition to the spin-glass phase. This limitation is thoroughly analyzed and motivated. △ Less

Submitted 22 December, 2014; v1 submitted 24 February, 2014; originally announced February 2014.

Comments: 14 pages, 12 figures

Journal ref: Physical Review E 90 (1), 012112 (2014)

The renormalization flow of the hierarchical Anderson model at weak disorder

Authors: F. L. Metz, L. Leuzzi, G. Parisi

Abstract: We study the flow of the renormalized model parameters obtained from a sequence of simple transformations of the 1D Anderson model with long-range hierarchical hopping. Combining numerical results with a perturbative approach for the flow equations, we identify three qualitatively different regimes at weak disorder. For a sufficiently fast decay of the hopping energy, the Cauchy distribution is th… ▽ More We study the flow of the renormalized model parameters obtained from a sequence of simple transformations of the 1D Anderson model with long-range hierarchical hopping. Combining numerical results with a perturbative approach for the flow equations, we identify three qualitatively different regimes at weak disorder. For a sufficiently fast decay of the hopping energy, the Cauchy distribution is the only stable fixed-point of the flow equations, whereas for sufficiently slowly decaying hopping energy the renormalized parameters flow to a delta peak fixed-point distribution. In an intermediate range of the hopping decay, both fixed-point distributions are stable and the stationary solution is determined by the initial configuration of the random parameters. We present results for the critical decay of the hopping energy separating the different regimes. △ Less

Submitted 4 November, 2013; originally announced November 2013.

Comments: 8 pages, 6 figures

Journal ref: Phys. Rev. B 89, 064201 (2014)

Time-resolved dynamics of granular matter by random laser emission

Authors: Viola Folli, Neda Ghofraniha, Andrea Puglisi, Luca Leuzzi, Claudio Conti

Abstract: Because of the huge commercial importance of granular systems, the second-most used material in industry after water, intersecting the industry in multiple trades, like pharmacy and agriculture, fundamental research on grain-like materials has received an increasing amount of attention in the last decades. In photonics, the applications of granular materials have been only marginally investigated.… ▽ More Because of the huge commercial importance of granular systems, the second-most used material in industry after water, intersecting the industry in multiple trades, like pharmacy and agriculture, fundamental research on grain-like materials has received an increasing amount of attention in the last decades. In photonics, the applications of granular materials have been only marginally investigated. We report the first phase-diagram of a granular as obtained by laser emission. The dynamics of vertically-oscillated granular in a liquid solution in a three-dimensional container is investigated by employing its random laser emission. The granular motion is function of the frequency and amplitude of the mechanical solicitation, we show how the laser emission allows to distinguish two phases in the granular and analyze its spectral distribution. This constitutes a fundamental step in the field of granulars and gives a clear evidence of the possible control on light-matter interaction achievable in grain-like system. △ Less

Submitted 22 October, 2013; originally announced October 2013.

Comments: 16 pages, 7 figures

Journal ref: Scientific Reports 3, 2251 (2013)

Imry-Ma criterion for long-range random field Ising model: short-/long-range equivalence in a field

Authors: Luca Leuzzi, Giorgio Parisi

Abstract: The Ising model in a random field and with power-law decaying ferromagnetic bonds is studied at zero temperature. Comparing the scaling of the energy contributions of the ferromagnetic domain wall flip and of the random field a la Imry-Ma we obtain a threshold value for the power $ρろー$ of the long-range interaction, beyond which no critical behavior occurs. The critical threshold value is $ρろー_c=3/2$,… ▽ More The Ising model in a random field and with power-law decaying ferromagnetic bonds is studied at zero temperature. Comparing the scaling of the energy contributions of the ferromagnetic domain wall flip and of the random field a la Imry-Ma we obtain a threshold value for the power $ρろー$ of the long-range interaction, beyond which no critical behavior occurs. The critical threshold value is $ρろー_c=3/2$, at a difference with the zero field model in which $ρろー_c=2$. This prediction is confirmed by numerical computation of the ground states below, at, and above this threshold value. Some possible implications for the critical behavior of spin-glasses in a field are conjectured. △ Less

Submitted 25 March, 2013; originally announced March 2013.

Comments: 5 pages, 6 figures

Journal ref: Phys. Rev. B 88, 224204 (2013)

Renormalization Group on hierarchical lattices in finite dimensional disordered Ising and Blume-Emery-Griffiths Models

Authors: F. Antenucci, A. Crisanti, L. Leuzzi

Abstract: Renormalization group on hierarchical lattices is often considered a valuable tool to understand the critical behavior of more complicated statistical mechanical models. In presence of quenched disorder, however, in many model cases predictions obtained with the Migdal-Kadanoff bond removal approach fail to quantitatively and qualitatively reproduce critical properties obtained in the mean-field a… ▽ More Renormalization group on hierarchical lattices is often considered a valuable tool to understand the critical behavior of more complicated statistical mechanical models. In presence of quenched disorder, however, in many model cases predictions obtained with the Migdal-Kadanoff bond removal approach fail to quantitatively and qualitatively reproduce critical properties obtained in the mean-field approximation or by numerical simulations in finite dimensions. In order to critically review this limitation we analyze the behavior of Ising and Blume-Emery-Griffiths models on more complicated hierarchical lattices. We find that, apart from some exceptions, the different behavior appears not only limited to Midgal-Kadanoff-like cells but is associated right to the hierarchization of Bravais lattices in small cells also when in-cell loops are considered. △ Less

Submitted 10 January, 2014; v1 submitted 15 March, 2013; originally announced March 2013.

Comments: 36 pages, 18 figures, 10 tables

Journal ref: Journal of Statistical Physics 155 (5), 909-931 (2014)

Transition between localized and extended states in the hierarchical Anderson model

Authors: F. L. Metz, L. Leuzzi, G. Parisi, V. Sacksteder IV

Abstract: We present strong numerical evidence for the existence of a localization-delocalization transition in the eigenstates of the 1-D Anderson model with long-range hierarchical hopping. Hierarchical models are important because of the well-known mapping between their phases and those of models with short range hopping in higher dimensions, and also because the renormalization group can be applied exac… ▽ More We present strong numerical evidence for the existence of a localization-delocalization transition in the eigenstates of the 1-D Anderson model with long-range hierarchical hopping. Hierarchical models are important because of the well-known mapping between their phases and those of models with short range hopping in higher dimensions, and also because the renormalization group can be applied exactly without the approximations that generally are required in other models. In the hierarchical Anderson model we find a finite critical disorder strength Wc where the average inverse participation ratio goes to zero; at small disorder W < Wc the model lies in a delocalized phase. This result is based on numerical calculation of the inverse participation ratio in the infinite volume limit using an exact renormalization group approach facilitated by the model's hierarchical structure. Our results are consistent with the presence of an Anderson transition in short-range models with D > 2 dimensions, which was predicted using renormalization group arguments. Our finding should stimulate interest in the hierarchical Anderson model as a simplified and tractable model of the Anderson localization transition which occurs in finite-dimensional systems with short-range hopping. △ Less

Submitted 22 April, 2013; v1 submitted 8 March, 2013; originally announced March 2013.

Comments: 7 pages, 7 figures

Exactly solvable spin-glass models with ferromagnetic couplings: the spherical multi-$p$-spin model in a self-induced field

Authors: Andrea Crisanti, Luca Leuzzi

Abstract: We report some results on the quenched disordered Spherical multi-$p$-Spin Model in presence of ferromagnetic couplings. In particular, we present the phase diagrams of some representative cases that schematically describe, in the mean-field approximation, the behavior of most known transitions in glassy materials, including dynamic arrest in super-cooled liquids, amorphous-amorphous transitions a… ▽ More We report some results on the quenched disordered Spherical multi-$p$-Spin Model in presence of ferromagnetic couplings. In particular, we present the phase diagrams of some representative cases that schematically describe, in the mean-field approximation, the behavior of most known transitions in glassy materials, including dynamic arrest in super-cooled liquids, amorphous-amorphous transitions and spin-glass transitions. A simplified notation is introduced in order to compute systems properties in terms of an effective, self-induced, field encoding the whole ferromagnetic information. △ Less

Submitted 8 January, 2013; originally announced January 2013.

Comments: 34 pages, 7 figures

Journal ref: Nucl. Phys. B 870, 176-204 (2013)

Critical behavior of the XY model in complex topologies

Authors: Miguel Ibáñez Berganza, Luca Leuzzi

Abstract: The critical behavior of the O(2) model on dilute Levy graphs built on a 2D square lattice is analyzed. Different qualitative cases are probed, varying the exponent rho governing the dependence on the distance of the connectivity probability distribution. The mean-field regime, as well as the long-range and short-range non-mean-field regimes are investigated by means of high-performance parallel M… ▽ More The critical behavior of the O(2) model on dilute Levy graphs built on a 2D square lattice is analyzed. Different qualitative cases are probed, varying the exponent rho governing the dependence on the distance of the connectivity probability distribution. The mean-field regime, as well as the long-range and short-range non-mean-field regimes are investigated by means of high-performance parallel Monte-Carlo numerical simulations running on GPUs. The relationship between the long-range rho exponent and the effective dimension of an equivalent short-range system with the same critical behavior is investigated. Evidence is provided for the effective short-range dimension to coincide with the spectral dimension of the Levy graph for the XY model in the mean-field regime. △ Less

Submitted 15 August, 2013; v1 submitted 16 November, 2012; originally announced November 2012.

Comments: 18 pages, 22 figures. Minor corrections, improved estimations of critical quantities

Journal ref: Physical Review B 88 (14), 144104 (2013)

Shaken Granular Lasers

Authors: Viola Folli, Andrea Puglisi, Luca Leuzzi, Claudio Conti

Abstract: Granular materials have been studied for decades, also driven by industrial and technological applications. These very simple systems, composed by agglomerations of mesoscopic particles, are characterized, in specific regimes, by a large number of metastable states and an extreme sensitivity (e.g., in sound transmission) on the arrangement of grains; they are not substantially affected by thermal… ▽ More Granular materials have been studied for decades, also driven by industrial and technological applications. These very simple systems, composed by agglomerations of mesoscopic particles, are characterized, in specific regimes, by a large number of metastable states and an extreme sensitivity (e.g., in sound transmission) on the arrangement of grains; they are not substantially affected by thermal phenomena, but can be controlled by mechanical solicitations. Laser emission from shaken granular matter is so far unexplored; here we provide experimental evidence that it can be affected and controlled by the status of motion of the granular, we also find that competitive random lasers can be observed. We hence demonstrate the potentialities of gravity affected moving disordered materials for optical applications, and open the road to a variety of novel interdisciplinary investigations, involving modern statistical mechanics and disordered photonics. △ Less

Submitted 27 May, 2012; originally announced May 2012.

Comments: 4 pages, 3 figures. To be published in Physical Review Letters

Journal ref: Phys. Rev. Lett. 108, 248002 (2012)

Following states in temperature in the spherical s+p-spin glass model

Authors: YiFan Sun, Andrea Crisanti, Florent Krzakala, Luca Leuzzi, Lenka Zdeborová

Abstract: In many mean-field glassy systems, the low-temperature Gibbs measure is dominated by exponentially many metastable states. We analyze the evolution of the metastable states as temperature changes adiabatically in the solvable case of the spherical $s+p$-spin glass model, extending the work of Barrat, Franz and Parisi J. Phys. A 30, 5593 (1997). We confirm the presence of level crossings, bifurcati… ▽ More In many mean-field glassy systems, the low-temperature Gibbs measure is dominated by exponentially many metastable states. We analyze the evolution of the metastable states as temperature changes adiabatically in the solvable case of the spherical $s+p$-spin glass model, extending the work of Barrat, Franz and Parisi J. Phys. A 30, 5593 (1997). We confirm the presence of level crossings, bifurcations, and temperature chaos. For the states that are at equilibrium close to the so-called dynamical temperature $T_d$, we find, however, that the following state method (and the dynamical solution of the model as well) is intrinsically limited by the vanishing of solutions with non-zero overlap at low temperature. △ Less

Submitted 17 April, 2012; originally announced April 2012.

Journal ref: J. Stat. Mech. (2012) P07002

Critical slowing down exponents in quenched disordered spin models for structural glasses: Random Orthogonal and related models

Authors: Francesco Caltagirone, Ulisse Ferrari, Luca Leuzzi, Giorgio Parisi, Tommaso Rizzo

Abstract: An important prediction of Mode-Coupling-Theory (MCT) is the relationship between the power- law decay exponents in the βべーた regime. In the original structural glass context this relationship follows from the MCT equations that are obtained making rather uncontrolled approximations and λらむだ has to be treated like a tunable parameter. It is known that a certain class of mean-field spin-glass models is ex… ▽ More An important prediction of Mode-Coupling-Theory (MCT) is the relationship between the power- law decay exponents in the βべーた regime. In the original structural glass context this relationship follows from the MCT equations that are obtained making rather uncontrolled approximations and λらむだ has to be treated like a tunable parameter. It is known that a certain class of mean-field spin-glass models is exactly described by MCT equations. In this context, the physical meaning of the so called parameter exponent λらむだ has recently been unveiled, giving a method to compute it exactly in a static framework. In this paper we exploit this new technique to compute the critical slowing down exponents in a class of mean-field Ising spin-glass models including, as special cases, the Sherrington-Kirkpatrick model, the p-spin model and the Random Orthogonal model. △ Less

Submitted 10 March, 2012; originally announced March 2012.

Comments: 7 pages, 3 figures

Journal ref: Phys. Rev. B 86, 064204 (2012)