Showing 1–50 of 165 results for author: Marinari, E

A Very Effective and Simple Diffusion Reconstruction for the Diluted Ising Model

Authors: Stefano Bae, Enzo Marinari, Federico Ricci-Tersenghi

Abstract: Diffusion-based generative models are machine learning models that use diffusion processes to learn the probability distribution of high-dimensional data. In recent years, they have become extremely successful in generating multimedia content. However, it is still unknown if such models can be used to generate high-quality datasets of physical models. In this work, we use a Landau-Ginzburg-like di… ▽ More Diffusion-based generative models are machine learning models that use diffusion processes to learn the probability distribution of high-dimensional data. In recent years, they have become extremely successful in generating multimedia content. However, it is still unknown if such models can be used to generate high-quality datasets of physical models. In this work, we use a Landau-Ginzburg-like diffusion model to infer the distribution of a $2D$ bond-diluted Ising model. Our approach is simple and effective, and we show that the generated samples reproduce correctly the statistical and critical properties of the physical model. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Multiscaling in the 3D critical site-diluted Ising ferromagnet

Authors: E. Marinari, V. Martin-Mayor, G. Parisi, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo

Abstract: We have studied numerically the appearance of multiscaling behavior in the three-dimensional ferromagnetic Ising site diluted model, in the form of a multifractal distribution of the decay exponents for the spatial correlation functions at the critical temperature. We have computed the exponents of the long-distance decay of higher moments of the correlation function, up to the 10th power, by stud… ▽ More We have studied numerically the appearance of multiscaling behavior in the three-dimensional ferromagnetic Ising site diluted model, in the form of a multifractal distribution of the decay exponents for the spatial correlation functions at the critical temperature. We have computed the exponents of the long-distance decay of higher moments of the correlation function, up to the 10th power, by studying three different quantities: global susceptibilities, local susceptibilities and correlation functions. We have found very clear evidences for multiscaling behavior. △ Less

Submitted 16 January, 2024; v1 submitted 13 September, 2023; originally announced September 2023.

Comments: 19 pages and 5 figures. Final version of the paper

Journal ref: J. Stat. Mech. (2024) 013301

Eigenvector Dreaming

Authors: Marco Benedetti, Louis Carillo, Enzo Marinari, Marc Mèzard

Abstract: Among the performance-enhancing procedures for Hopfield-type networks that implement associative memory, Hebbian Unlearning (or dreaming) strikes for its simplicity and its clear biological interpretation. Yet, it does not easily lend itself to a clear analytical understanding. Here we show how Hebbian Unlearning can be effectively described in terms of a simple evolution of the spectrum and the e… ▽ More Among the performance-enhancing procedures for Hopfield-type networks that implement associative memory, Hebbian Unlearning (or dreaming) strikes for its simplicity and its clear biological interpretation. Yet, it does not easily lend itself to a clear analytical understanding. Here we show how Hebbian Unlearning can be effectively described in terms of a simple evolution of the spectrum and the eigenvectors of the coupling matrix. We use these ideas to design new dreaming algorithms that are effective from a computational point of view, and are analytically far more transparent than the original scheme. △ Less

Submitted 9 August, 2023; originally announced August 2023.

Quantifying memory in spin glasses

Authors: Janus Collaboration, I. Paga, J. He, M. Baity-Jesi, E. Calore, A. Cruz, L. A. Fernandez, J. M. Gil-Narvion, I. Gonzalez-Adalid Pemartin, A. Gordillo-Guerrero, D. Iñiguez, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Moreno-Gordo, A. Muñoz Sudupe, D. Navarro, R. L. Orbach, G. Parisi, S. Perez-Gaviro, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, D. L. Schlagel, B. Seoane , et al. (2 additional authors not shown)

Abstract: Rejuvenation and memory, long considered the distinguishing features of spin glasses, have recently been proven to result from the growth of multiple length scales. This insight, enabled by simulations on the Janus~II supercomputer, has opened the door to a quantitative analysis. We combine numerical simulations with comparable experiments to introduce two coefficients that quantify memory. A thir… ▽ More Rejuvenation and memory, long considered the distinguishing features of spin glasses, have recently been proven to result from the growth of multiple length scales. This insight, enabled by simulations on the Janus~II supercomputer, has opened the door to a quantitative analysis. We combine numerical simulations with comparable experiments to introduce two coefficients that quantify memory. A third coefficient has been recently presented by Freedberg et al. We show that these coefficients are physically equivalent by studying their temperature and waiting-time dependence. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 10 pages, 8 figures

Multifractality in spin glasses

Authors: Janus Collaboration, M. Baity-Jesi, E. Calore, A. Cruz, L. A. Fernandez, J. M. Gil-Narvion, I. Gonzalez-Adalid Pemartin, A. Gordillo-Guerrero, D. Iñiguez, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Moreno-Gordo, A. Muñoz Sudupe, D. Navarro, I. Paga, G. Parisi, S. Perez-Gaviro, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, B. Seoane, A. Tarancon, D. Yllanes

Abstract: We unveil the multifractal behavior of Ising spin glasses in their low-temperature phase. Using the Janus II custom-built supercomputer, the spin-glass correlation function is studied locally. Dramatic fluctuations are found when pairs of sites at the same distance are compared. The scaling of these fluctuations, as the spin-glass coherence length grows with time, is characterized through the comp… ▽ More We unveil the multifractal behavior of Ising spin glasses in their low-temperature phase. Using the Janus II custom-built supercomputer, the spin-glass correlation function is studied locally. Dramatic fluctuations are found when pairs of sites at the same distance are compared. The scaling of these fluctuations, as the spin-glass coherence length grows with time, is characterized through the computation of the singularity spectrum and its corresponding Legendre transform. A comparatively small number of site pairs controls the average correlation that governs the response to a magnetic field. We explain how this scenario of dramatic fluctuations (at length scales smaller than the coherence length) can be reconciled with the smooth, self-averaging behavior that has long been considered to describe spin-glass dynamics. △ Less

Submitted 22 January, 2024; v1 submitted 7 June, 2023; originally announced June 2023.

Comments: 15 pages, 15 figures

Journal ref: PNAS 121, e2312880120 (2024)

Towards a robust criterion of anomalous diffusion

Authors: Vittoria Sposini, Diego Krapf, Enzo Marinari, Raimon Sunyer, Felix Ritort, Fereydoon Taheri, Christine Selhuber-Unkel, Rebecca Benelli, Matthias Weiss, Ralf Metzler, Gleb Oshanin

Abstract: Anomalous-diffusion, the departure of the spreading dynamics of diffusing particles from the traditional law of Brownian-motion, is a signature feature of a large number of complex soft-matter and biological systems. Anomalous-diffusion emerges due to a variety of physical mechanisms, e.g., trapping interactions or the viscoelasticity of the environment. However, sometimes systems dynamics are err… ▽ More Anomalous-diffusion, the departure of the spreading dynamics of diffusing particles from the traditional law of Brownian-motion, is a signature feature of a large number of complex soft-matter and biological systems. Anomalous-diffusion emerges due to a variety of physical mechanisms, e.g., trapping interactions or the viscoelasticity of the environment. However, sometimes systems dynamics are erroneously claimed to be anomalous, despite the fact that the true motion is Brownian -- or vice versa. This ambiguity in establishing whether the dynamics as normal or anomalous can have far-reaching consequences, e.g., in predictions for reaction- or relaxation-laws. Demonstrating that a system exhibits normal- or anomalous-diffusion is highly desirable for a vast host of applications. Here, we present a criterion for anomalous-diffusion based on the method of power-spectral analysis of single trajectories. The robustness of this criterion is studied for trajectories of fractional-Brownian-motion, a ubiquitous stochastic process for the description of anomalous-diffusion, in the presence of two types of measurement errors. In particular, we find that our criterion is very robust for subdiffusion. Various tests on surrogate data in absence or presence of additional positional noise demonstrate the efficacy of this method in practical contexts. Finally, we provide a proof-of-concept based on diverse experiments exhibiting both normal and anomalous-diffusion. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 13 pages, 6 figures, RevTeX

On the superposition principle and non-linear response in spin glasses

Authors: I. Paga, Q. Zhai, M. Baity-Jesi, E. Calore, A. Cruz, C. Cummings, L. A. Fernandez, J. M. Gil-Narvion, I. Gonzalez-Adalid Pemartin, A. Gordillo-Guerrero, D. Iñiguez, G. G. Kenning, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Moreno-Gordo, A. Muñoz-Sudupe, D. Navarro, R. L. Orbach, G. Parisi, S. Perez-Gaviro, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, D. L. Schlagel , et al. (3 additional authors not shown)

Abstract: The extended principle of superposition has been a touchstone of spin glass dynamics for almost thirty years. The Uppsala group has demonstrated its validity for the metallic spin glass, CuMn, for magnetic fields $H$ up to 10 Oe at the reduced temperature $T_\mathrm{r}=T/T_\mathrm{g} = 0.95$, where $T_\mathrm{g}$ is the spin glass condensation temperature. For $H > 10$ Oe, they observe a departure… ▽ More The extended principle of superposition has been a touchstone of spin glass dynamics for almost thirty years. The Uppsala group has demonstrated its validity for the metallic spin glass, CuMn, for magnetic fields $H$ up to 10 Oe at the reduced temperature $T_\mathrm{r}=T/T_\mathrm{g} = 0.95$, where $T_\mathrm{g}$ is the spin glass condensation temperature. For $H > 10$ Oe, they observe a departure from linear response which they ascribe to the development of non-linear dynamics. The thrust of this paper is to develop a microscopic origin for this behavior by focusing on the time development of the spin glass correlation length, $ξ(t,t_\mathrm{w};H)$. Here, $t$ is the time after $H$ changes, and $t_\mathrm{w}$ is the time from the quench for $T>T_\mathrm{g}$ to the working temperature $T$ until $H$ changes. We connect the growth of $ξ(t,t_\mathrm{w};H)$ to the barrier heights $Δ(t_\mathrm{w})$ that set the dynamics. The effect of $H$ on the magnitude of $Δ(t_\mathrm{w})$ is responsible for affecting differently the two dynamical protocols associated with turning $H$ off (TRM, or thermoremanent magnetization) or on (ZFC, or zero field-cooled magnetization). In this paper, we display the difference between the zero-field cooled $ξ_{\text {ZFC}}(t,t_\mathrm{w};H)$ and the thermoremanent magnetization $ξ_{\text {TRM}}(t,t_\mathrm{w};H)$ correlation lengths as $H$ increases, both experimentally and through numerical simulations, corresponding to the violation of the extended principle of superposition in line with the finding of the Uppsala Group. △ Less

Submitted 26 June, 2023; v1 submitted 21 July, 2022; originally announced July 2022.

Comments: 23 pages and 18 figures

Journal ref: Phys. Rev. B 107, 214436- Published 23 June 2023

Memory and rejuvenation in spin glasses: aging systems are ruled by more than one length scale

Authors: Janus Collaboration, M. Baity-Jesi, E. Calore, A. Cruz, L. A. Fernandez, J. M. Gil-Narvion, I. Gonzalez-Adalid Pemartin, A. Gordillo-Guerrero, D. Iñiguez, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Moreno-Gordo, A. Muñoz-Sudupe, D. Navarro, I. Paga, G. Parisi, S. Perez-Gaviro, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, B. Seoane, A. Tarancon, R. Tripiccione, D. Yllanes

Abstract: Memory and rejuvenation effects in the magnetic response of off-equilibrium spin glasses have been widely regarded as the doorway into the experimental exploration of ultrametricity and temperature chaos (maybe the most exotic features in glassy free-energy landscapes). Unfortunately, despite more than twenty years of theoretical efforts following the experimental discovery of memory and rejuvenat… ▽ More Memory and rejuvenation effects in the magnetic response of off-equilibrium spin glasses have been widely regarded as the doorway into the experimental exploration of ultrametricity and temperature chaos (maybe the most exotic features in glassy free-energy landscapes). Unfortunately, despite more than twenty years of theoretical efforts following the experimental discovery of memory and rejuvenation, these effects have thus far been impossible to simulate reliably. Yet, three recent developments convinced us to accept this challenge: first, the custom-built Janus II supercomputer makes it possible to carry out "numerical experiments" in which the very same quantities that can be measured in single crystals of CuMn are computed from the simulation, allowing for parallel analysis of the simulation/experiment data. Second, Janus II simulations have taught us how numerical and experimental length scales should be compared. Third, we have recently understood how temperature chaos materializes in aging dynamics. All three aspects have proved crucial for reliably reproducing rejuvenation and memory effects on the computer. Our analysis shows that (at least) three different length scales play a key role in aging dynamics, while essentially all theoretical analyses of the aging dynamics emphasize the presence and the crucial role of a single glassy correlation length. △ Less

Submitted 30 August, 2022; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: Main text:13 pages, 4 figures SM: 7 pages, 8 figures

Journal ref: Nature Physics 19, 978-985 (2023)

Noise-to-signal ratio of single-trajectory spectral densities in centered Gaussian processes

Authors: Alessio Squarcini, Enzo Marinari, Gleb Oshanin, Luca Peliti, Lamberto Rondoni

Abstract: We discuss the statistical properties of a single-trajectory power spectral density $S(ω,\mathcal{T})$ of an arbitrary real-valued centered Gaussian process $X(t)$, where $ω$ is the angular frequency and $\mathcal{T}$ the observation time. We derive a double-sided inequality for its noise-to-signal ratio and obtain the full probability density function of $S(ω,\mathcal{T})$. Our findings imply tha… ▽ More We discuss the statistical properties of a single-trajectory power spectral density $S(ω,\mathcal{T})$ of an arbitrary real-valued centered Gaussian process $X(t)$, where $ω$ is the angular frequency and $\mathcal{T}$ the observation time. We derive a double-sided inequality for its noise-to-signal ratio and obtain the full probability density function of $S(ω,\mathcal{T})$. Our findings imply that the fluctuations of $S(ω,\mathcal{T})$ exceed its average value $μ(ω,\mathcal{T})$. This implies that using $μ(ω,\mathcal{T})$ to describe the behavior of these processes can be problematic. We finally evaluate the typical behavior of $S(ω,\mathcal{T})$ and find that it deviates markedly from the average $μ(ω,\mathcal{T})$ in most cases. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Comments: 10 pages

Frequency-frequency correlations of single-trajectory spectral densities of Gaussian processes

Authors: Alessio Squarcini, Enzo Marinari, Gleb Oshanin, Luca Peliti, Lamberto Rondoni

Abstract: We investigate the stochastic behavior of the single-trajectory spectral density $S(ω,\mathcal{T})$ of several Gaussian stochastic processes, i.e., Brownian motion, the Ornstein-Uhlenbeck process, the Brownian gyrator model and fractional Brownian motion, as a function of the frequency $ω$ and the observation time $\mathcal{T}$. We evaluate in particular the variance and the frequency-frequency co… ▽ More We investigate the stochastic behavior of the single-trajectory spectral density $S(ω,\mathcal{T})$ of several Gaussian stochastic processes, i.e., Brownian motion, the Ornstein-Uhlenbeck process, the Brownian gyrator model and fractional Brownian motion, as a function of the frequency $ω$ and the observation time $\mathcal{T}$. We evaluate in particular the variance and the frequency-frequency correlation of $S(ω,\mathcal{T})$ for different values of $ω$. We show that these properties exhibit different behaviors for different physical cases and can therefore be used as a sensitive probe discriminating between different kinds of random motion. These results may prove quite useful in the analysis of experimental data. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Comments: 29 pages, 7 figures

Spectral fingerprints of non-equilibrium dynamics: The case of a Brownian gyrator

Authors: S. Cerasoli, S. Ciliberto, E. Marinari, G. Oshanin, L. Peliti, L. Rondoni

Abstract: The same system can exhibit a completely different dynamical behavior when it evolves in equilibrium conditions or when it is driven out-of-equilibrium by, e.g., connecting some of its components to heat baths kept at different temperatures. Here we concentrate on an analytically solvable and experimentally-relevant model of such a system -- the so-called Brownian gyrator -- a two-dimensional nano… ▽ More The same system can exhibit a completely different dynamical behavior when it evolves in equilibrium conditions or when it is driven out-of-equilibrium by, e.g., connecting some of its components to heat baths kept at different temperatures. Here we concentrate on an analytically solvable and experimentally-relevant model of such a system -- the so-called Brownian gyrator -- a two-dimensional nanomachine that performs a systematic, on average, rotation around the origin under non-equilibrium conditions, while no net rotation takes place in equilibrium. On this example, we discuss a question whether it is possible to distinguish between two types of a behavior judging not upon the statistical properties of the trajectories of components, but rather upon their respective spectral densities. The latter are widely used to characterize diverse dynamical systems and are routinely calculated from the data using standard built-in packages. From such a perspective, we inquire whether the power spectral densities possess some "fingerprint" properties specific to the behavior in non-equilibrium. We show that indeed one can conclusively distinguish between equilibrium and non-equilibrium dynamics by analyzing the cross-correlations between the spectral densities of both components in the short frequency limit, or from the spectral densities of both components evaluated at zero frequency. Our analytical predictions, corroborated by experimental and numerical results, open a new direction for the analysis of a non-equilibrium dynamics. △ Less

Submitted 6 August, 2022; v1 submitted 13 January, 2022; originally announced January 2022.

Comments: 23 pages, 6 figures

Journal ref: Phys. Rev. E 106, 014137 (2022)

Supervised perceptron learning vs unsupervised Hebbian unlearning: Approaching optimal memory retrieval in Hopfield-like networks

Authors: Marco Benedetti, Enrico Ventura, Enzo Marinari, Giancarlo Ruocco, Francesco Zamponi

Abstract: The Hebbian unlearning algorithm, i.e. an unsupervised local procedure used to improve the retrieval properties in Hopfield-like neural networks, is numerically compared to a supervised algorithm to train a linear symmetric perceptron. We analyze the stability of the stored memories: basins of attraction obtained by the Hebbian unlearning technique are found to be comparable in size to those obtai… ▽ More The Hebbian unlearning algorithm, i.e. an unsupervised local procedure used to improve the retrieval properties in Hopfield-like neural networks, is numerically compared to a supervised algorithm to train a linear symmetric perceptron. We analyze the stability of the stored memories: basins of attraction obtained by the Hebbian unlearning technique are found to be comparable in size to those obtained in the symmetric perceptron, while the two algorithms are found to converge in the same region of Gardner's space of interactions, having followed similar learning paths. A geometric interpretation of Hebbian unlearning is proposed to explain its optimal performances. Because the Hopfield model is also a prototypical model of disordered magnetic system, it might be possible to translate our results to other models of interest for memory storage in materials. △ Less

Submitted 7 March, 2022; v1 submitted 31 December, 2021; originally announced January 2022.

Comments: 11 pages, 13 figures, 43 references

Journal ref: J. Chem. Phys. 156, 104107 (2022)

How we are leading a 3-XORSAT challenge: from the energy landscape to the algorithm and its efficient implementation on GPUs

Authors: M. Bernaschi, M. Bisson, M. Fatica, E. Marinari, V. Martin-Mayor, G. Parisi, F. Ricci-Tersenghi

Abstract: A recent 3-XORSAT challenge required to minimize a very complex and rough energy function, typical of glassy models with a random first order transition and a golf course like energy landscape. We present the ideas beyond the quasi-greedy algorithm and its very efficient implementation on GPUs that are allowing us to rank first in such a competition. We suggest a better protocol to compare algorit… ▽ More A recent 3-XORSAT challenge required to minimize a very complex and rough energy function, typical of glassy models with a random first order transition and a golf course like energy landscape. We present the ideas beyond the quasi-greedy algorithm and its very efficient implementation on GPUs that are allowing us to rank first in such a competition. We suggest a better protocol to compare algorithmic performances and we also provide analytical predictions about the exponential growth of the times to find the solution in terms of free-energy barriers. △ Less

Submitted 24 February, 2021; v1 submitted 18 February, 2021; originally announced February 2021.

Comments: 7 pages, 7 figure, EPL format + SM (2 pages)

Journal ref: EPL, 133 (2021) 60005

Recognition Capabilities of a Hopfield Model with Auxiliary Hidden Neurons

Authors: Marco Benedetti, Victor Dotsenko, Giulia Fischetti, Enzo Marinari, Gleb Oshanin

Abstract: We study the recognition capabilities of the Hopfield model with auxiliary hidden layers, which emerge naturally upon a Hubbard-Stratonovich transformation. We show that the recognition capabilities of such a model at zero-temperature outperform those of the original Hopfield model, due to a substantial increase of the storage capacity and the lack of a naturally defined basin of attraction. The m… ▽ More We study the recognition capabilities of the Hopfield model with auxiliary hidden layers, which emerge naturally upon a Hubbard-Stratonovich transformation. We show that the recognition capabilities of such a model at zero-temperature outperform those of the original Hopfield model, due to a substantial increase of the storage capacity and the lack of a naturally defined basin of attraction. The modified model does not fall abruptly in a regime of complete confusion when memory load exceeds a sharp threshold. △ Less

Submitted 13 January, 2021; originally announced January 2021.

Journal ref: Phys. Rev. E 103, 060401 (2021)

Spin-glass dynamics in the presence of a magnetic field: exploration of microscopic properties

Authors: I. Paga, Q. Zhai, M. Baity-Jesi, E. Calore, A. Cruz, L. A. Fernandez, J. M. Gil-Narvion, I. Gonzalez-Adalid Pemartin, A Gordillo-Guerrero, D. Iñiguez, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Moreno-Gordo, A. Muñoz-Sudupe, D. Navarro, R. L. Orbach, G. Parisi, S. Perez-Gaviro, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, D. L. Schlagel, B. Seoane, A. Tarancon , et al. (2 additional authors not shown)

Abstract: The synergy between experiment, theory, and simulations enables a microscopic analysis of spin-glass dynamics in a magnetic field in the vicinity of and below the spin-glass transition temperature $T_\mathrm{g}$. The spin-glass correlation length, $ξ(t,t_\mathrm{w};T)$, is analysed both in experiments and in simulations in terms of the waiting time $t_\mathrm{w}$ after the spin glass has been cool… ▽ More The synergy between experiment, theory, and simulations enables a microscopic analysis of spin-glass dynamics in a magnetic field in the vicinity of and below the spin-glass transition temperature $T_\mathrm{g}$. The spin-glass correlation length, $ξ(t,t_\mathrm{w};T)$, is analysed both in experiments and in simulations in terms of the waiting time $t_\mathrm{w}$ after the spin glass has been cooled down to a stabilised measuring temperature $T<T_\mathrm{g}$ and of the time $t$ after the magnetic field is changed. This correlation length is extracted experimentally for a CuMn 6 at. % single crystal, as well as for simulations on the Janus II special-purpose supercomputer, the latter with time and length scales comparable to experiment. The non-linear magnetic susceptibility is reported from experiment and simulations, using $ξ(t,t_\mathrm{w};T)$ as the scaling variable. Previous experiments are reanalysed, and disagreements about the nature of the Zeeman energy are resolved. The growth of the spin-glass magnetisation in zero-field magnetisation experiments, $M_\mathrm{ZFC}(t,t_\mathrm{w};T)$, is measured from simulations, verifying the scaling relationships in the dynamical or non-equilibrium regime. Our preliminary search for the de Almeida-Thouless line in $D=3$ is discussed. △ Less

Submitted 10 March, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: 49 pages, figures 25

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

Temperature chaos is present in off-equilibrium spin-glass dynamics

Authors: Janus Collaboration, M. Baity-Jesi, E. Calore, A. Cruz, L. A. Fernandez, J. M. Gil-Narvion, I. Gonzalez-Adalid Pemartin, A. Gordillo-Guerrero, D. Iñiguez, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Moreno-Gordo, A. Muñoz-Sudupe, D. Navarro, I. Paga, G. Parisi, S. Perez-Gaviro, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, B. Seoane, A. Tarancon, R. Tripiccione, D. Yllanes

Abstract: We find a dynamic effect in the non-equilibrium dynamics of a spin glass that closely parallels equilibrium temperature chaos. This effect, that we name dynamic temperature chaos, is spatially heterogeneous to a large degree. The key controlling quantity is the time-growing spin-glass coherence length. Our detailed characterization of dynamic temperature chaos paves the way for the analysis of rec… ▽ More We find a dynamic effect in the non-equilibrium dynamics of a spin glass that closely parallels equilibrium temperature chaos. This effect, that we name dynamic temperature chaos, is spatially heterogeneous to a large degree. The key controlling quantity is the time-growing spin-glass coherence length. Our detailed characterization of dynamic temperature chaos paves the way for the analysis of recent and forthcoming experiments. This work has been made possible thanks to the most massive simulation to date of non-equilibrium dynamics, carried out on the Janus~II custom-built supercomputer. △ Less

Submitted 6 July, 2021; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: Version accepted for publication in Communication Physics 10 pages, 9 figures

Journal ref: Communications Physics 4, 74 (2021)

Scaling law describes the spin-glass response in theory, experiments and simulations

Authors: Q. Zhai, I. Paga, M. Baity-Jesi, E. Calore, A. Cruz, L. A. Fernandez, J. M. Gil-Narvion, I. Gonzalez-Adalid Pemartin, A. Gordillo-Guerrero, D. Iñiguez, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Moreno-Gordo, A. Muñoz-Sudupe, D. Navarro, R. L. Orbach, G. Parisi, S. Perez-Gaviro, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, D. L. Schlagel, B. Seoane, A. Tarancon , et al. (2 additional authors not shown)

Abstract: The correlation length $ξ$, a key quantity in glassy dynamics, can now be precisely measured for spin glasses both in experiments and in simulations. However, known analysis methods lead to discrepancies either for large external fields or close to the glass temperature. We solve this problem by introducing a scaling law that takes into account both the magnetic field and the time-dependent spin-g… ▽ More The correlation length $ξ$, a key quantity in glassy dynamics, can now be precisely measured for spin glasses both in experiments and in simulations. However, known analysis methods lead to discrepancies either for large external fields or close to the glass temperature. We solve this problem by introducing a scaling law that takes into account both the magnetic field and the time-dependent spin-glass correlation length. The scaling law is successfully tested against experimental measurements in a CuMn single crystal and against large-scale simulations on the Janus II dedicated computer. △ Less

Submitted 30 November, 2020; v1 submitted 7 July, 2020; originally announced July 2020.

Comments: Revised version, including supplemental material

Journal ref: Phys. Rev. Lett. 125, 237202 (2020)

Competing endogenous RNA crosstalk at system level

Authors: Mattia Miotto, Enzo Marinari, Andrea De Martino

Abstract: microRNAs (miRNAs) regulate gene expression at post-transcriptional level by repressing target RNA molecules. Competition to bind miRNAs tends in turn to correlate their targets, establishing effective RNA-RNA interactions that can influence expression levels, buffer fluctuations and promote signal propagation. Such a potential has been characterized mathematically for small motifs both at steady… ▽ More microRNAs (miRNAs) regulate gene expression at post-transcriptional level by repressing target RNA molecules. Competition to bind miRNAs tends in turn to correlate their targets, establishing effective RNA-RNA interactions that can influence expression levels, buffer fluctuations and promote signal propagation. Such a potential has been characterized mathematically for small motifs both at steady state and \red{away from stationarity}. Experimental evidence, on the other hand, suggests that competing endogenous RNA (ceRNA) crosstalk is rather weak. Extended miRNA-RNA networks could however favour the integration of many crosstalk interactions, leading to significant large-scale effects in spite of the weakness of individual links. To clarify the extent to which crosstalk is sustained by the miRNA interactome, we have studied its emergent systemic features in silico in large-scale miRNA-RNA network reconstructions. We show that, although generically weak, system-level crosstalk patterns (i) are enhanced by transcriptional heterogeneities, (ii) can achieve high-intensity even for RNAs that are not co-regulated, (iii) are robust to variability in transcription rates, and (iv) are significantly non-local, i.e. correlate weakly with miRNA-RNA interaction parameters. Furthermore, RNA levels are generically more stable when crosstalk is strongest. As some of these features appear to be encoded in the network's topology, crosstalk may functionally be favoured by natural selection. These results suggest that, besides their repressive role, miRNAs mediate a weak but resilient and context-independent network of cross-regulatory interactions that interconnect the transcriptome, stabilize expression levels and support system-level responses. △ Less

Submitted 21 October, 2019; originally announced October 2019.

Comments: 25 pages, includes Supporting Information; to appear in PLoS Comp Biol

Passive advection of fractional Brownian motion by random layered flows

Authors: Alessio Squarcini, Enzo Marinari, Gleb Oshanin

Abstract: We study statistical properties of the process $Y(t)$ of a passive advection by quenched random layered flows in situations when the inter-layer transfer is governed by a fractional Brownian motion $X(t)$ with the Hurst index $H \in (0,1)$. We show that the disorder-averaged mean-squared displacement of the passive advection grows in the large time $t$ limit in proportion to $t^{2 - H}$, which def… ▽ More We study statistical properties of the process $Y(t)$ of a passive advection by quenched random layered flows in situations when the inter-layer transfer is governed by a fractional Brownian motion $X(t)$ with the Hurst index $H \in (0,1)$. We show that the disorder-averaged mean-squared displacement of the passive advection grows in the large time $t$ limit in proportion to $t^{2 - H}$, which defines a family of anomalous super-diffusions. We evaluate the disorder-averaged Wigner-Ville spectrum of the advection process $Y(t)$ and demonstrate that it has a rather unusual power-law form $1/f^{3 - H}$ with a characteristic exponent which exceed the value $2$. Our results also suggest that sample-to-sample fluctuations of the spectrum can be very important. △ Less

Submitted 21 September, 2019; originally announced September 2019.

Comments: 18 pages, 4 figures

Dimensional crossover in the aging dynamics of spin glasses in a film geometry

Authors: L. A. Fernandez, E. Marinari, V. Martin-Mayor, I. Paga, J. J. Ruiz-Lorenzo

Abstract: Motivated by recent experiments of exceptional accuracy, we study numerically the spin-glass dynamics in a film geometry. We cover all the relevant time regimes, from picoseconds to equilibrium, at temperatures at and below the 3D critical point. The dimensional crossover from 3D to 2D dynamics, that starts when the correlation length becomes comparable to the film thickness, consists of four dyna… ▽ More Motivated by recent experiments of exceptional accuracy, we study numerically the spin-glass dynamics in a film geometry. We cover all the relevant time regimes, from picoseconds to equilibrium, at temperatures at and below the 3D critical point. The dimensional crossover from 3D to 2D dynamics, that starts when the correlation length becomes comparable to the film thickness, consists of four dynamical regimes. Our analysis, based on a Renormalization Group transformation, finds consistent the overall physical picture employed by Orbach et al. in the interpretation of their experiments. △ Less

Submitted 14 February, 2020; v1 submitted 4 June, 2019; originally announced June 2019.

Comments: 8 pages, 8 figures. Version accepted for publication in Phys. Rev. B

Journal ref: Phys. Rev. B 100, 184412 (2019)

The Essential Role of Thermodynamics in metabolic network modeling: physical insights and computational challenges

Authors: A. De Martino, D. De Martino, E. Marinari

Abstract: Quantitative studies of cell metabolism are often based on large chemical reaction network models. A steady state approach is suited to analyze phenomena on the timescale of cell growth and circumvents the problem of incomplete experimental knowledge on kinetic laws and parameters, but it shall be supported by a correct implementation of thermodynamic constraints. In this article we review the lat… ▽ More Quantitative studies of cell metabolism are often based on large chemical reaction network models. A steady state approach is suited to analyze phenomena on the timescale of cell growth and circumvents the problem of incomplete experimental knowledge on kinetic laws and parameters, but it shall be supported by a correct implementation of thermodynamic constraints. In this article we review the latter aspect highlighting its computational challenges and physical insights. The simple introduction of Gibbs inequalities avoids the presence of unfeasible loops allowing for correct timescale analysis but leads to possibly non-convex feasible flux spaces, whose exploration needs efficient algorithms. We shorty review on the implementation of thermodynamics through variational principles in constraints based models of metabolic networks. △ Less

Submitted 19 February, 2019; originally announced February 2019.

Comments: Contribution to "Chemical Kinetics: Beyond the Textbook", edited by Katja Lindenberg, Ralf Metzler and Gleb Oshanin (World Scientific, Singapore 2019)

Covariance of the running range of a Brownian trajectory

Authors: Brandon Annesi, Enzo Marinari, Gleb Oshanin

Abstract: The question how the extremal values of a stochastic process achieved on different time intervals are correlated to each other has been discussed within the last few years on examples of the running maximum of a Brownian motion, of a Brownian Bridge and of a Slepian process. Here, we focus on the two-time correlations of the running range of Brownian motion - the maximal extent of a Brownian traje… ▽ More The question how the extremal values of a stochastic process achieved on different time intervals are correlated to each other has been discussed within the last few years on examples of the running maximum of a Brownian motion, of a Brownian Bridge and of a Slepian process. Here, we focus on the two-time correlations of the running range of Brownian motion - the maximal extent of a Brownian trajectory on a finite time interval. We calculate exactly the covariance function of the running range and analyse its asymptotic behaviour. Our analysis reveals non-trivial correlations between the value of the largest descent (rise) of a BM from the top to a bottom on some time interval, and the value of this property on a larger time interval. △ Less

Submitted 19 February, 2019; originally announced February 2019.

Comments: 13 pages, 3 figures

Spectral content of a single non-Brownian trajectory

Authors: D. Krapf, N. Lukat, E. Marinari, R. Metzler, G. Oshanin, C. Selhuber-Unkel, A. Squarcini, L. Stadler, M. Weiss, X. Xu

Abstract: Time-dependent processes are often analysed using the power spectral density (PSD), calculated by taking an appropriate Fourier transform of individual trajectories and finding the associated ensemble-average. Frequently, the available experimental data sets are too small for such ensemble averages, and hence it is of a great conceptual and practical importance to understand to which extent releva… ▽ More Time-dependent processes are often analysed using the power spectral density (PSD), calculated by taking an appropriate Fourier transform of individual trajectories and finding the associated ensemble-average. Frequently, the available experimental data sets are too small for such ensemble averages, and hence it is of a great conceptual and practical importance to understand to which extent relevant information can be gained from $S(f,T)$, the PSD of a single trajectory. Here we focus on the behavior of this random, realization-dependent variable, parametrized by frequency $f$ and observation-time $T$, for a broad family of anomalous diffusions---fractional Brownian motion (fBm) with Hurst-index $H$---and derive exactly its probability density function. We show that $S(f,T)$ is proportional---up to a random numerical factor whose universal distribution we determine---to the ensemble-averaged PSD. For subdiffusion ($H<1/2$) we find that $S(f,T)\sim A/f^{2H+1}$ with random-amplitude $A$. In sharp contrast, for superdiffusion $(H>1/2)$ $S(f,T)\sim BT^{2H-1}/f^2$ with random amplitude $B$. Remarkably, for $H>1/2$ the PSD exhibits the same frequency-dependence as Brownian motion, a deceptive property that may lead to false conclusions when interpreting experimental data. Notably, for $H>1/2$ the PSD is ageing and is dependent on $T$. Our predictions for both sub- and superdiffusion are confirmed by experiments in live cells and in agarose hydrogels, and by extensive simulations. △ Less

Submitted 1 February, 2019; originally announced February 2019.

Comments: 13 pages, 5 figures, Supplemental Material can be found at https://journals.aps.org/prx/supplemental/10.1103/PhysRevX.9.011019/prx_SM_final.pdf

Journal ref: Phys. Rev. X 9, 011019 (2019)

Forgetting Memories and their Attractiveness

Authors: Enzo Marinari

Abstract: We study numerically the memory which forgets, introduced in 1986 by Parisi by bounding the synaptic strength, with a mechanism which avoid confusion, allows to remember the pattern learned more recently and has a physiologically very well defined meaning. We analyze a number of features of the learning at finite number of neurons and finite number of patterns. We discuss how the system behaves in… ▽ More We study numerically the memory which forgets, introduced in 1986 by Parisi by bounding the synaptic strength, with a mechanism which avoid confusion, allows to remember the pattern learned more recently and has a physiologically very well defined meaning. We analyze a number of features of the learning at finite number of neurons and finite number of patterns. We discuss how the system behaves in the large but finite N limit. We analyze the basin of attraction of the patterns that have been learned, and we show that it is exponentially small in the age of the pattern. This is a clearly non physiological feature of the model. △ Less

Submitted 31 May, 2018; originally announced May 2018.

The out-equilibrium 2D Ising spin glass: almost, but not quite, a free-field theory

Authors: L. A. Fernandez, E. Marinari, V. Martin-Mayor, G. Parisi, J. J. Ruiz-Lorenzo

Abstract: We consider the spatial correlation function of the two-dimensional Ising spin glass under out-equilibrium conditions. We pay special attention to the scaling limit reached upon approaching zero temperature. The field-theory of a non-interacting field makes a surprisingly good job at describing the spatial shape of the correlation function of the out-equilibrium Edwards-Anderson Ising model in two… ▽ More We consider the spatial correlation function of the two-dimensional Ising spin glass under out-equilibrium conditions. We pay special attention to the scaling limit reached upon approaching zero temperature. The field-theory of a non-interacting field makes a surprisingly good job at describing the spatial shape of the correlation function of the out-equilibrium Edwards-Anderson Ising model in two dimensions. △ Less

Submitted 24 October, 2018; v1 submitted 22 May, 2018; originally announced May 2018.

Comments: 20 pages + 5 Figures

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

An experiment-oriented analysis of 2D spin-glass dynamics: a twelve time-decades scaling study

Authors: L. A. Fernandez, E. Marinari, V. Martin-Mayor, G. Parisi, J. J. Ruiz-Lorenzo

Abstract: Recent high precision experimental results on spin-glass films ask for a detailed understanding of the domain-growth dynamics of two-dimensional spin glasses. To achieve this goal, we numerically simulate the out-equilibrium dynamics of the Ising spin glass for a time that spans close to twelve orders of magnitude (from picoseconds to order of a second), in systems large enough to avoid finite-siz… ▽ More Recent high precision experimental results on spin-glass films ask for a detailed understanding of the domain-growth dynamics of two-dimensional spin glasses. To achieve this goal, we numerically simulate the out-equilibrium dynamics of the Ising spin glass for a time that spans close to twelve orders of magnitude (from picoseconds to order of a second), in systems large enough to avoid finite-size effects. We find that the time-growth of the size of the glassy domains is excellently described by a single scaling function. A single time-scale $τ(T)$ controls the dynamics. $τ(T)$ diverges upon approaching the $T=0$ critical point. The divergence of $τ(T\to 0)$ is Arrhenius-like, with a barrier height that depends very mildly on temperature. The growth of this barrier-height is best described by critical dynamics. As a side product we obtain an impressive confirmation of universality of the equilibrium behavior of two-dimensional spin-glasses. △ Less

Submitted 6 May, 2019; v1 submitted 17 May, 2018; originally announced May 2018.

Comments: 21 pages, 9 figures. Updated references. Added DOI and Journal ref

Journal ref: L A Fernandez et al 2019 J. Phys. A: Math. Theor. 52 224002

The Mpemba effect in spin glasses is a persistent memory effect

Authors: Janus collaboration, M. Baity-Jesi, E. Calore, A. Cruz, L. A. Fernandez, J. M. Gil-Narvion, A. Gordillo-Guerrero, D. Iñiguez, A. Lasanta, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Moreno-Gordo, A. Muñoz-Sudupe, D. Navarro, G. Parisi, S. Perez-Gaviro, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, B. Seoane, A. Tarancon, R. Tripiccione, D. Yllanes

Abstract: The Mpemba effect occurs when a hot system cools faster than an initially colder one, when both are refrigerated in the same thermal reservoir. Using the custom built supercomputer Janus II, we study the Mpemba effect in spin glasses and show that it is a non-equilibrium process, governed by the coherence length ξof the system. The effect occurs when the bath temperature lies in the glassy phase,… ▽ More The Mpemba effect occurs when a hot system cools faster than an initially colder one, when both are refrigerated in the same thermal reservoir. Using the custom built supercomputer Janus II, we study the Mpemba effect in spin glasses and show that it is a non-equilibrium process, governed by the coherence length ξof the system. The effect occurs when the bath temperature lies in the glassy phase, but it is not necessary for the thermal protocol to cross the critical temperature. In fact, the Mpemba effect follows from a strong relationship between the internal energy and ξthat turns out to be a sure-tell sign of being in the glassy phase. Thus, the Mpemba effect presents itself as an intriguing new avenue for the experimental study of the coherence length in supercooled liquids and other glass formers. △ Less

Submitted 5 July, 2019; v1 submitted 20 April, 2018; originally announced April 2018.

Comments: Version accepted for publication in PNAS. 6 pages, 7 figures

Journal ref: PNAS 116 (2019), 15350-15355

Aging rate of spin glasses from simulations matches experiments

Authors: Janus Collaboration, M. Baity-Jesi, E. Calore, A. Cruz, L. A. Fernandez, J. M. Gil-Narvion, A. Gordillo-Guerrero, D. Iñiguez, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Moreno-Gordo, A. Muñoz-Sudupe, D. Navarro, G. Parisi, S. Perez-Gaviro, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, B. Seoane, A. Tarancon, R. Tripiccione, D. Yllanes

Abstract: Experiments on spin glasses can now make precise measurements of the exponent $z(T)$ governing the growth of glassy domains, while our computational capabilities allow us to make quantitative predictions for experimental scales. However, experimental and numerical values for $z(T)$ have differed. We use new simulations on the Janus II computer to resolve this discrepancy, finding a time-dependent… ▽ More Experiments on spin glasses can now make precise measurements of the exponent $z(T)$ governing the growth of glassy domains, while our computational capabilities allow us to make quantitative predictions for experimental scales. However, experimental and numerical values for $z(T)$ have differed. We use new simulations on the Janus II computer to resolve this discrepancy, finding a time-dependent $z(T, t_w)$, which leads to the experimental value through mild extrapolations. Furthermore, theoretical insight is gained by studying a crossover between the $T = T_c$ and $T = 0$ fixed points. △ Less

Submitted 18 June, 2018; v1 submitted 6 March, 2018; originally announced March 2018.

Comments: Version accepted for publication in PRL. 12 pages, 9 figures

Journal ref: Phys. Rev. Lett. 120, 267203 (2018)

Power spectral density of a single Brownian trajectory: What one can and cannot learn from it

Authors: Diego Krapf, Enzo Marinari, Ralf Metzler, Gleb Oshanin, Xinran Xu, Alessio Squarcini

Abstract: The power spectral density (PSD) of any time-dependent stochastic processes $X_t$ is a meaningful feature of its spectral content. In its text-book definition, the PSD is the Fourier transform of the covariance function of $X_t$ over an infinitely large observation time $T$, that is, it is defined as an ensemble-averaged property taken in the limit $T \to \infty$. A legitimate question is what inf… ▽ More The power spectral density (PSD) of any time-dependent stochastic processes $X_t$ is a meaningful feature of its spectral content. In its text-book definition, the PSD is the Fourier transform of the covariance function of $X_t$ over an infinitely large observation time $T$, that is, it is defined as an ensemble-averaged property taken in the limit $T \to \infty$. A legitimate question is what information on the PSD can be reliably obtained from single-trajectory experiments, if one goes beyond the standard definition and analyzes the PSD of a \textit{single} trajectory recorded for a \textit{finite} observation time $T$. In quest for this answer, for a $d$-dimensional Brownian motion we calculate the probability density function of a single-trajectory PSD for arbitrary frequency $f$, finite observation time $T$ and arbitrary number $k$ of projections of the trajectory on different axes. We show analytically that the scaling exponent for the frequency-dependence of the PSD specific to an ensemble of Brownian motion trajectories can be already obtained from a single trajectory, while the numerical amplitude in the relation between the ensemble-averaged and single-trajectory PSDs is a fluctuating property which varies from realization to realization. The distribution of this amplitude is calculated exactly and is discussed in detail. Our results are confirmed by numerical simulations and single particle tracking experiments, with remarkably good agreement. In addition we consider a truncated Wiener representation of Brownian motion, and the case of a discrete-time lattice random walk. We highlight some differences in the behavior of a single-trajectory PSD for Brownian motion and for the two latter situations. The framework developed herein will allow for meaningful physical analysis of experimental stochastic trajectories. △ Less

Submitted 9 January, 2018; originally announced January 2018.

Comments: 24 pages, 13 figures

Journal ref: New J Phys 20, 023029 (2018)

Numerical evidences of universal trap-like aging dynamics

Authors: Chiara Cammarota, Enzo Marinari

Abstract: Trap models have been initially proposed as toy models for dynamical relaxation in extremely simplified rough potential energy landscapes. Their importance has considerably grown recently thanks to the discovery that the trap like aging mechanism is directly controlling the out-of-equilibrium relaxation processes of more sophisticated spin models, that are considered as the solvable counterpart of… ▽ More Trap models have been initially proposed as toy models for dynamical relaxation in extremely simplified rough potential energy landscapes. Their importance has considerably grown recently thanks to the discovery that the trap like aging mechanism is directly controlling the out-of-equilibrium relaxation processes of more sophisticated spin models, that are considered as the solvable counterpart of real disordered systems. Establishing on a firmer ground the connection between these spin model out-of-equilibrium behavior and the trap like aging mechanism would shed new light on the properties, still largely mysterious, of the activated out-of-equilibrium dynamics of disordered systems. In this work we discuss numerical evidences of emergent trap-like aging behavior in a variety of disordered models. Our numerical results are backed by analytic derivations and heuristic discussions. Such exploration reveals some of the tricks needed to analyze the trap behavior in spite of the occurrence of secondary processes, of the existence of dynamical correlations and of finite system's size effects. △ Less

Submitted 22 October, 2017; originally announced October 2017.

Comments: 25 pages, 15 figures

Dynamic Variational Study of Chaos: Spin Glasses in Three Dimensions

Authors: A. Billoire, L. A. Fernandez, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Moreno-Gordo, G. Parisi, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo

Abstract: We have introduced a variational method to improve the computation of integrated correlation times in the Parallel Tempering Dynamics, obtaining a better estimate (a lower bound, at least) of the exponential correlation time. Using this determination of the correlation times, we revisited the problem of the characterization of the chaos in temperature in finite dimensional spin glasses by means of… ▽ More We have introduced a variational method to improve the computation of integrated correlation times in the Parallel Tempering Dynamics, obtaining a better estimate (a lower bound, at least) of the exponential correlation time. Using this determination of the correlation times, we revisited the problem of the characterization of the chaos in temperature in finite dimensional spin glasses by means of the study of correlations between different chaos indicators computed in the static and the correlation times of the Parallel Tempering dynamics. The sample-distribution of the characteristic time for the Parallel Tempering dynamics turns out to be fat-tailed and it obeys finite-size scaling. △ Less

Submitted 13 December, 2017; v1 submitted 28 September, 2017; originally announced September 2017.

Comments: New appendix added. Version to appear in JSTAT

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

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

Matching microscopic and macroscopic responses in glasses

Authors: Janus Collaboration, M. Baity-Jesi, E. Calore, A. Cruz, L. A. Fernandez, J. M. Gil-Narvion, A. Gordillo-Guerrero, D. Iñiguez, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Monforte-Garcia, A. Muñoz-Sudupe, D. Navarro, G. Parisi, S. Perez-Gaviro, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, B. Seoane, A. Tarancon, R. Tripiccione, D. Yllanes

Abstract: We first reproduce on the Janus and Janus II computers a milestone experiment that measures the spin-glass coherence length through the lowering of free-energy barriers induced by the Zeeman effect. Secondly we determine the scaling behavior that allows a quantitative analysis of a new experiment reported in the companion Letter [S. Guchhait and R. Orbach, Phys. Rev. Lett. 118, 157203 (2017)]. The… ▽ More We first reproduce on the Janus and Janus II computers a milestone experiment that measures the spin-glass coherence length through the lowering of free-energy barriers induced by the Zeeman effect. Secondly we determine the scaling behavior that allows a quantitative analysis of a new experiment reported in the companion Letter [S. Guchhait and R. Orbach, Phys. Rev. Lett. 118, 157203 (2017)]. The value of the coherence length estimated through the analysis of microscopic correlation functions turns out to be quantitatively consistent with its measurement through macroscopic response functions. Further, non-linear susceptibilities, recently measured in glass-forming liquids, scale as powers of the same microscopic length. △ Less

Submitted 25 April, 2017; originally announced April 2017.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. Lett. 118, 157202 (2017)

Numerical construction of the Aizenman-Wehr metastate

Authors: A. Billoire, L. A. Fernandez, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Moreno-Gordo, G. Parisi, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo

Abstract: Chaotic size dependence makes it extremely difficult to take the thermodynamic limit in disordered systems. Instead, the metastate, which is a distribution over thermodynamic states, might have a smooth limit. So far, studies of the metastate have been mostly mathematical. We present a numerical construction of the metastate for the d=3 Ising spin glass. We work in equilibrium, below the critical… ▽ More Chaotic size dependence makes it extremely difficult to take the thermodynamic limit in disordered systems. Instead, the metastate, which is a distribution over thermodynamic states, might have a smooth limit. So far, studies of the metastate have been mostly mathematical. We present a numerical construction of the metastate for the d=3 Ising spin glass. We work in equilibrium, below the critical temperature. Leveraging recent rigorous results, our numerical analysis gives evidence for a "dispersed" metastate, supported on many thermodynamic states. △ Less

Submitted 5 April, 2017; originally announced April 2017.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. Lett. 119, 037203 (2017)

A yield-cost tradeoff governs Escherichia coli's decision between fermentation and respiration in carbon-limited growth

Authors: Matteo Mori, Enzo Marinari, Andrea De Martino

Abstract: Many microbial systems are known to actively reshape their proteomes in response to changes in growth conditions induced e.g. by nutritional stress or antibiotics. Part of the re-allocation accounts for the fact that, as the growth rate is limited by targeting specific metabolic activities, cells simply respond by fine-tuning their proteome to invest more resources into the limiting activity (i.e.… ▽ More Many microbial systems are known to actively reshape their proteomes in response to changes in growth conditions induced e.g. by nutritional stress or antibiotics. Part of the re-allocation accounts for the fact that, as the growth rate is limited by targeting specific metabolic activities, cells simply respond by fine-tuning their proteome to invest more resources into the limiting activity (i.e. by synthesizing more proteins devoted to it). However, this is often accompanied by an overall re-organization of metabolism, aimed at improving the growth yield under limitation by re-wiring resource through different pathways. While both effects impact proteome composition, the latter underlies a more complex systemic response to stress. By focusing on E. coli's `acetate switch', we use mathematical modeling and a re-analysis of empirical data to show that the transition from a predominantly fermentative to a predominantly respirative metabolism in carbon-limited growth results from the trade-off between maximizing the growth yield and minimizing its costs in terms of required the proteome share. In particular, E. coli's metabolic phenotypes appear to be Pareto-optimal for these objective functions over a broad range of dilutions. △ Less

Submitted 2 March, 2017; originally announced March 2017.

Comments: 13 pages (includes Supporting Text)

A statics-dynamics equivalence through the fluctuation-dissipation ratio provides a window into the spin-glass phase from nonequilibrium measurements

Authors: Janus Collaboration, M. Baity-Jesi, E. Calore, A. Cruz, L. A. Fernandez, J. M. Gil-Narvion, A. Gordillo-Guerrero, D. Iñiguez, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Monforte-Garcia, A. Muñoz-Sudupe, D. Navarro, G. Parisi, S. Perez-Gaviro, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, B. Seoane, A. Tarancon, R. Tripiccione, D. Yllanes

Abstract: The unifying feature of glass formers (such as polymers, supercooled liquids, colloids, granulars, spin glasses, superconductors, ...) is a sluggish dynamics at low temperatures. Indeed, their dynamics is so slow that thermal equilibrium is never reached in macroscopic samples: in analogy with living beings, glasses are said to age. Here, we show how to relate experimentally relevant quantities wi… ▽ More The unifying feature of glass formers (such as polymers, supercooled liquids, colloids, granulars, spin glasses, superconductors, ...) is a sluggish dynamics at low temperatures. Indeed, their dynamics is so slow that thermal equilibrium is never reached in macroscopic samples: in analogy with living beings, glasses are said to age. Here, we show how to relate experimentally relevant quantities with the experimentally unreachable low-temperature equilibrium phase. We have performed a very accurate computation of the non-equilibrium fluctuation-dissipation ratio for the three-dimensional Edwards-Anderson Ising spin glass, by means of large-scale simulations on the special-purpose computers Janus and Janus II. This ratio (computed for finite times on very large, effectively infinite, systems) is compared with the equilibrium probability distribution of the spin overlap for finite sizes. The resulting quantitative statics-dynamics dictionary, based on observables that can be measured with current experimental methods, could allow the experimental exploration of important features of the spin-glass phase without uncontrollable extrapolations to infinite times or system sizes. △ Less

Submitted 24 April, 2017; v1 submitted 5 October, 2016; originally announced October 2016.

Comments: Version accepted for publication in PNAS. Reuploaded to fix a typo in the list of authors

Journal ref: Proc. Natl. Acad. Sci. USA 114 (2017), 1838-1843

Constrained Allocation Flux Balance Analysis

Authors: Matteo Mori, Terence Hwa, Olivier C. Martin, Andrea De Martino, Enzo Marinari

Abstract: New experimental results on bacterial growth inspire a novel top-down approach to study cell metabolism, combining mass balance and proteomic constraints to extend and complement Flux Balance Analysis. We introduce here Constrained Allocation Flux Balance Analysis, CAFBA, in which the biosynthetic costs associated to growth are accounted for in an effective way through a single additional genome-w… ▽ More New experimental results on bacterial growth inspire a novel top-down approach to study cell metabolism, combining mass balance and proteomic constraints to extend and complement Flux Balance Analysis. We introduce here Constrained Allocation Flux Balance Analysis, CAFBA, in which the biosynthetic costs associated to growth are accounted for in an effective way through a single additional genome-wide constraint. Its roots lie in the experimentally observed pattern of proteome allocation for metabolic functions, allowing to bridge regulation and metabolism in a transparent way under the principle of growth-rate maximization. We provide a simple method to solve CAFBA efficiently and propose an "ensemble averaging" procedure to account for unknown protein costs. Applying this approach to modeling E. coli metabolism, we find that, as the growth rate increases, CAFBA solutions cross over from respiratory, growth-yield maximizing states (preferred at slow growth) to fermentative states with carbon overflow (preferred at fast growth). In addition, CAFBA allows for quantitatively accurate predictions on the rate of acetate excretion and growth yield based on only 3 parameters determined by empirical growth laws. △ Less

Submitted 1 July, 2016; originally announced July 2016.

Comments: 21 pages, 6 figures (main) + 33 pages, various figures and tables (supporting); for the supplementary MatLab code, see http://tinyurl.com/h763es8

Journal ref: PLoS Comput Biol 12(6): e1004913 (2016)

Sample-to-sample fluctuations of power spectrum of a random motion in a periodic Sinai model

Authors: David S. Dean, Antonio Iorio, Enzo Marinari, Gleb Oshanin

Abstract: The Sinai model of a tracer diffusing in a quenched Brownian potential is a much studied problem exhibiting a logarithmically slow anomalous diffusion due to the growth of energy barriers with the system size. However, if the potential is random but periodic, the regime of anomalous diffusion crosses over to one of normal diffusion once a tracer has diffused over a few periods of the system. Here… ▽ More The Sinai model of a tracer diffusing in a quenched Brownian potential is a much studied problem exhibiting a logarithmically slow anomalous diffusion due to the growth of energy barriers with the system size. However, if the potential is random but periodic, the regime of anomalous diffusion crosses over to one of normal diffusion once a tracer has diffused over a few periods of the system. Here we consider a system in which the potential is given by a Brownian Bridge on a finite interval $(0,L)$ and then periodically repeated over the whole real line, and study the power spectrum $S(f)$ of the diffusive process $x(t)$ in such a potential. We show that for most of realizations of $x(t)$ in a given realization of the potential, the low-frequency behavior is $S(f) \sim {\cal A}/f^2$, i.e., the same as for standard Brownian motion, and the amplitude ${\cal A}$ is a disorder-dependent random variable with a finite support. Focusing on the statistical properties of this random variable, we determine the moments of ${\cal A}$ of arbitrary, negative or positive order $k$, and demonstrate that they exhibit a multi-fractal dependence on $k$, and a rather unusual dependence on the temperature and on the periodicity $L$, which are supported by atypical realizations of the periodic disorder. We finally show that the distribution of ${\cal A}$ has a log-normal left tail, and exhibits an essential singularity close to the right edge of the support, which is related to the Lifshitz singularity. Our findings are based both on analytic results and on extensive numerical simulations of the process $x(t)$. △ Less

Submitted 30 June, 2016; originally announced July 2016.

Comments: 8 pages, 5 figures

Journal ref: Phys. Rev. E 94, 032131 (2016)

Temperature chaos is a non-local effect

Authors: L. A. Fernandez, E. Marinari, V. Martin-Mayor, G. Parisi, D. Yllanes

Abstract: Temperature chaos plays a role in important effects, like for example memory and rejuvenation, in spin glasses, colloids, polymers. We numerically investigate temperature chaos in spin glasses, exploiting its recent characterization as a rare-event driven phenomenon. The peculiarities of the transformation from periodic to anti-periodic boundary conditions in spin glasses allow us to conclude that… ▽ More Temperature chaos plays a role in important effects, like for example memory and rejuvenation, in spin glasses, colloids, polymers. We numerically investigate temperature chaos in spin glasses, exploiting its recent characterization as a rare-event driven phenomenon. The peculiarities of the transformation from periodic to anti-periodic boundary conditions in spin glasses allow us to conclude that temperature chaos is non-local: no bounded region of the system causes it. We precise the statistical relationship between temperature chaos and the free-energy changes upon varying boundary conditions. △ Less

Submitted 11 October, 2016; v1 submitted 10 May, 2016; originally announced May 2016.

Comments: 15 pages, 8 figures. Version accepted for publication in JSTAT

Journal ref: J. Stat. Mech. (2016) 123301

Universal critical behavior of the 2d Ising spin glass

Authors: L. A. Fernandez, E. Marinari, V. Martin-Mayor, G. Parisi, J. J. Ruiz-Lorenzo

Abstract: We use finite size scaling to study Ising spin glasses in two spatial dimensions. The issue of universality is addressed by comparing discrete and continuous probability distributions for the quenched random couplings. The sophisticated temperature dependency of the scaling fields is identified as the major obstacle that has impeded a complete analysis. Once temperature is relinquished in favor of… ▽ More We use finite size scaling to study Ising spin glasses in two spatial dimensions. The issue of universality is addressed by comparing discrete and continuous probability distributions for the quenched random couplings. The sophisticated temperature dependency of the scaling fields is identified as the major obstacle that has impeded a complete analysis. Once temperature is relinquished in favor of the correlation length as the basic variable, we obtain a reliable estimation of the anomalous dimension and of the thermal critical exponent. Universality among binary and Gaussian couplings is confirmed to a high numerical accuracy. △ Less

Submitted 15 April, 2016; originally announced April 2016.

Comments: 11 pages,7 figures

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

Noise Processing by MicroRNA-Mediated Circuits: the Incoherent Feed-Forward Loop, Revisited

Authors: Silvia Grigolon, Francesca Di Patti, Andrea De Martino, Enzo Marinari

Abstract: The intrinsic stochasticity of gene expression is usually mitigated in higher eukaryotes by post-transcriptional regulation channels that stabilise the output layer, most notably protein levels. The discovery of small non-coding RNAs (miRNAs) in specific motifs of the genetic regulatory network has led to identifying noise buffering as the possible key function they exert in regulation. Recent in… ▽ More The intrinsic stochasticity of gene expression is usually mitigated in higher eukaryotes by post-transcriptional regulation channels that stabilise the output layer, most notably protein levels. The discovery of small non-coding RNAs (miRNAs) in specific motifs of the genetic regulatory network has led to identifying noise buffering as the possible key function they exert in regulation. Recent in vitro} and in silico studies have corroborated this hypothesis. It is however also known that miRNA-mediated noise reduction is hampered by transcriptional bursting in simple topologies. Here, using stochastic simulations validated by analytical calculations based on van Kampen's expansion, we revisit the noise-buffering capacity of the miRNA-mediated Incoherent Feed Forward Loop (IFFL), a small module that is widespread in the gene regulatory networks of higher eukaryotes, in order to account for the effects of intermittency in the transcriptional activity of the modulator gene. We show that bursting considerably alters the circuit's ability to control static protein noise. By comparing with other regulatory architectures, we find that direct transcriptional regulation significantly outperforms the IFFL in a broad range of kinetic parameters. This suggests that, under pulsatile inputs, static noise reduction may be less important than dynamical aspects of noise and information processing in characterising the performance of regulatory elements. △ Less

Submitted 11 April, 2016; originally announced April 2016.

Comments: 25 pages (Main Text and Supplementary Information), 5 figures

Journal ref: Heliyon 2 (2016) e00095

Probing the limits to microRNA-mediated control of gene expression

Authors: Araks Martirosyan, Matteo Figliuzzi, Enzo Marinari, Andrea De Martino

Abstract: According to the `ceRNA hypothesis', microRNAs (miRNAs) may act as mediators of an effective positive interaction between long coding or non-coding RNA molecules, carrying significant potential implications for a variety of biological processes. Here, inspired by recent work providing a quantitative description of small regulatory elements as information-conveying channels, we characterize the eff… ▽ More According to the `ceRNA hypothesis', microRNAs (miRNAs) may act as mediators of an effective positive interaction between long coding or non-coding RNA molecules, carrying significant potential implications for a variety of biological processes. Here, inspired by recent work providing a quantitative description of small regulatory elements as information-conveying channels, we characterize the effectiveness of miRNA-mediated regulation in terms of the optimal information flow achievable between modulator (transcription factors) and target nodes (long RNAs). Our findings show that, while a sufficiently large degree of target derepression is needed to activate miRNA-mediated transmission, (a) in case of differential mechanisms of complex processing and/or transcriptional capabilities, regulation by a post-transcriptional miRNA-channel can outperform that achieved through direct transcriptional control; moreover, (b) in the presence of large populations of weakly interacting miRNA molecules the extra noise coming from titration disappears, allowing the miRNA-channel to process information as effectively as the direct channel. These observations establish the limits of miRNA-mediated post-transcriptional cross-talk and suggest that, besides providing a degree of noise buffering, this type of control may be effectively employed in cells both as a failsafe mechanism and as a preferential fine tuner of gene expression, pointing to the specific situations in which each of these functionalities is maximized. △ Less

Submitted 26 January, 2016; originally announced January 2016.

Comments: 16 pages

Journal ref: PLoS Comput Biol 12: e1004715 (2016)

A numerical study of planar arrays of correlated spin islands

Authors: I. Maccari, A. Maiorano, E. Marinari, J. J. Ruiz-Lorenzo

Abstract: We present our analysis of a system of interacting islands of XY spins on a triangular lattice that has been introduced a few years ago by Eley et al. to account for the phenomenology in experiments on tunable arrays of proximity coupled long superconductor-normal metal-superconductor junctions. The main features of the model are the separation of a local and a global interaction energy scale and… ▽ More We present our analysis of a system of interacting islands of XY spins on a triangular lattice that has been introduced a few years ago by Eley et al. to account for the phenomenology in experiments on tunable arrays of proximity coupled long superconductor-normal metal-superconductor junctions. The main features of the model are the separation of a local and a global interaction energy scale and the mesoscopic character of the spin islands. Upon lowering the temperature the model undergoes two crossovers corresponding to an increasing phase coherence on a single island and to the onset of global coherence across the array; the latter is a thermodynamical phase transition in the Ising universality class. The dependence of the second transition on the island edge-to-edge spacing is related to the proximity-effect of the coupling constant. △ Less

Submitted 11 April, 2016; v1 submitted 15 September, 2015; originally announced September 2015.

Comments: 6 pages, 5 figures, accepted for publication in EPJ B

Journal ref: Eur. Phys. J. B (2016) 89: 127

Cross correlations of the American baby names

Authors: Paolo Barucca, Jacopo Rocchi, Enzo Marinari, Giorgio Parisi, Federico Ricci-Tersenghi

Abstract: The quantitative description of cultural evolution is a challenging task. The most difficult part of the problem is probably to find the appropriate measurable quantities that can make more quantitative such evasive concepts as, for example, dynamics of cultural movements, behavior patterns and traditions of the people. A strategy to tackle this issue is to observe particular features of human act… ▽ More The quantitative description of cultural evolution is a challenging task. The most difficult part of the problem is probably to find the appropriate measurable quantities that can make more quantitative such evasive concepts as, for example, dynamics of cultural movements, behavior patterns and traditions of the people. A strategy to tackle this issue is to observe particular features of human activities, i.e. cultural traits, such as names given to newborns. We study the names of babies born in the United States of America from 1910 to 2012. Our analysis shows that groups of different correlated states naturally emerge in different epochs, and we are able to follow and decrypt their evolution. While these groups of states are stable across many decades, a sudden reorganization occurs in the last part of the twentieth century. We think that this kind of quantitative analysis can be possibly extended to other cultural traits: although databases covering more than one century (as the one we used) are rare, the cultural evolution on shorter time scales can be studied thanks to the fact that many human activities are usually recorded in the present digital era. △ Less

Submitted 13 May, 2015; v1 submitted 10 October, 2014; originally announced October 2014.

Comments: submitted for consideration to PNAS

Journal ref: PNAS 112, 7943-7947 (2015)

Spontaneous energy-barrier formation in an entropy-driven glassy dynamics

Authors: Chiara Cammarota, Enzo Marinari

Abstract: The description of activated relaxation of glassy systems in the multidimensional configurational space is a long-standing open problem. We develop a phenomenological description of the out-of-equilibrium dynamics of a model with a rough potential energy landscape and we analyse it both numerically and analytically. The model provides an example of dynamics where typical relaxation channels go ove… ▽ More The description of activated relaxation of glassy systems in the multidimensional configurational space is a long-standing open problem. We develop a phenomenological description of the out-of-equilibrium dynamics of a model with a rough potential energy landscape and we analyse it both numerically and analytically. The model provides an example of dynamics where typical relaxation channels go over finite potential energy barriers despite the presence of less-energy-demanding escaping paths in configurational space; we expect this phenomenon to be also relevant in the thermally activated regime of realistic models of glass-formers. In this case, we found that typical dynamical paths episodically reach an high fixed threshold energy unexpectedly giving rise to a simple thermally activated aging phenomenology. In order to unveil this peculiar aging behavior we introduce a novel description of the dynamics in terms of spontaneously emerging dynamical basins. △ Less

Submitted 8 October, 2014; originally announced October 2014.

Comments: 5 pages, 3 figures

The cumulative overlap distribution function in realistic spin glasses

Authors: A. Billoire, A. Maiorano, E. Marinari, V. Martin-Mayor, D. Yllanes

Abstract: We use a sample-dependent analysis, based on medians and quantiles, to analyze the behavior of the overlap probability distribution of the Sherrington-Kirkpatrick and 3D Edwards-Anderson models of Ising spin glasses. We find that this approach is an effective tool to distinguish between RSB-like and droplet-like behavior of the spin-glass phase. Our results are in agreement with a RSB-like behavio… ▽ More We use a sample-dependent analysis, based on medians and quantiles, to analyze the behavior of the overlap probability distribution of the Sherrington-Kirkpatrick and 3D Edwards-Anderson models of Ising spin glasses. We find that this approach is an effective tool to distinguish between RSB-like and droplet-like behavior of the spin-glass phase. Our results are in agreement with a RSB-like behavior for the 3D Edwards-Anderson model. △ Less

Submitted 1 September, 2014; v1 submitted 6 June, 2014; originally announced June 2014.

Comments: Version accepted in PRB. 12 pages, 10 figures

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

The three dimensional Ising spin glass in an external magnetic field: the role of the silent majority

Authors: Janus Collaboration, M. Baity-Jesi, R. A. Banos, A. Cruz, L. A. Fernandez, J. M. Gil-Narvion, A. Gordillo-Guerrero, D. Iniguez, A. Maiorano, F. Mantovani, E. Marinari, V. Martin-Mayor, J. Monforte-Garcia, A. Munoz Sudupe, D. Navarro, G. Parisi, S. Perez-Gaviro, M. Pivanti, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, B. Seoane, A. Tarancon, R. Tripiccione, D. Yllanes

Abstract: We perform equilibrium parallel-tempering simulations of the 3D Ising Edwards-Anderson spin glass in a field. A traditional analysis shows no signs of a phase transition. Yet, we encounter dramatic fluctuations in the behaviour of the model: Averages over all the data only describe the behaviour of a small fraction of it. Therefore we develop a new approach to study the equilibrium behaviour of th… ▽ More We perform equilibrium parallel-tempering simulations of the 3D Ising Edwards-Anderson spin glass in a field. A traditional analysis shows no signs of a phase transition. Yet, we encounter dramatic fluctuations in the behaviour of the model: Averages over all the data only describe the behaviour of a small fraction of it. Therefore we develop a new approach to study the equilibrium behaviour of the system, by classifying the measurements as a function of a conditioning variate. We propose a finite-size scaling analysis based on the probability distribution function of the conditioning variate, which may accelerate the convergence to the thermodynamic limit. In this way, we find a non-trivial spectrum of behaviours, where a part of the measurements behaves as the average, while the majority of them shows signs of scale invariance. As a result, we can estimate the temperature interval where the phase transition in a field ought to lie, if it exists. Although this would-be critical regime is unreachable with present resources, the numerical challenge is finally well posed. △ Less

Submitted 2 April, 2014; v1 submitted 11 March, 2014; originally announced March 2014.

Comments: 42 pages, 19 figures. Minor changes and added figure (results unchanged)

Journal ref: J. Stat. Mech. (2014) P05014

RNA-based regulation: dynamics and response to perturbations of competing RNAs

Authors: Matteo Figliuzzi, Andrea De Martino, Enzo Marinari

Abstract: The observation that, through a titration mechanism, microRNAs (miRNAs) can act as mediators of effective interactions among their common targets (competing endogenous RNAs or ceRNAs) has brought forward the idea ('ceRNA hypothesis') that RNAs can regulate each other in extended 'cross-talk' networks. Such an ability might play a major role in post-transcriptional regulation (PTR) in shaping a cel… ▽ More The observation that, through a titration mechanism, microRNAs (miRNAs) can act as mediators of effective interactions among their common targets (competing endogenous RNAs or ceRNAs) has brought forward the idea ('ceRNA hypothesis') that RNAs can regulate each other in extended 'cross-talk' networks. Such an ability might play a major role in post-transcriptional regulation (PTR) in shaping a cell's protein repertoire. Recent work focusing on the emergent properties of the cross-talk networks has emphasized the high flexibility and selectivity that may be achieved at stationarity. On the other hand, dynamical aspects, possibly crucial on the relevant time scales, are far less clear. We have carried out a dynamical study of the ceRNA hypothesis on a model of PTR. Sensitivity analysis shows that ceRNA cross-talk is dynamically extended, i.e. it may take place on time scales shorter than those required to achieve stationairity even in cases where no cross-talk occurs in the steady state, and is possibly amplified. Besides, in case of large, transfection-like perturbations the system may develop strongly non-linear, threshold response. Finally, we show that the ceRNA effect provides a very efficient way for a cell to achieve fast positive shifts in the level of a ceRNA when necessary. These results indicate that competition for miRNAs may indeed provide an elementary mechanism to achieve system-level regulatory effects on the transcriptome over physiologically relevant time scales. △ Less

Submitted 19 December, 2013; originally announced December 2013.

Comments: Main text: 10 pages, 13 figures. Supporting Text: 3 pages, 6 figures

Energy metabolism and glutamate-glutamine cycle in the brain: a stoichiometric modeling perspective

Authors: Francesco Alessandro Massucci, Mauro DiNuzzo, Federico Giove, Bruno Maraviglia, Isaac Perez Castillo, Enzo Marinari, Andrea De Martino

Abstract: The energetics of cerebral activity critically relies on the functional and metabolic interactions between neurons and astrocytes. Important open questions include the relation between neuronal versus astrocytic energy demand, glucose uptake and intercellular lactate transfer, as well as their dependence on the level of activity. We have developed a large-scale, constraint-based network model of t… ▽ More The energetics of cerebral activity critically relies on the functional and metabolic interactions between neurons and astrocytes. Important open questions include the relation between neuronal versus astrocytic energy demand, glucose uptake and intercellular lactate transfer, as well as their dependence on the level of activity. We have developed a large-scale, constraint-based network model of the metabolic partnership between astrocytes and glutamatergic neurons that allows for a quantitative appraisal of the extent to which stoichiometry alone drives the energetics of the system. We find that the velocity of the glutamate-glutamine cycle ($V_{cyc}$) explains part of the uncoupling between glucose and oxygen utilization at increasing $V_{cyc}$ levels. Thus, we are able to characterize different activation states in terms of the tissue oxygen-glucose index (OGI). Calculations show that glucose is taken up and metabolized according to cellular energy requirements, and that partitioning of the sugar between different cell types is not significantly affected by $V_{cyc}$. Furthermore, both the direction and magnitude of the lactate shuttle between neurons and astrocytes turn out to depend on the relative cell glucose uptake while being roughly independent of $V_{cyc}$. These findings suggest that, in absence of ad hoc activity-related constraints on neuronal and astrocytic metabolism, the glutamate-glutamine cycle does not control the relative energy demand of neurons and astrocytes, and hence their glucose uptake and lactate exchange. △ Less

Submitted 24 October, 2013; originally announced October 2013.

Comments: 21 pages, incl. supporting text and tables

Journal ref: BMC Systems Biology 7:103 (2013)

Counting and correcting thermodynamically infeasible flux cycles in genome-scale metabolic networks

Authors: Daniele De Martino, Fabrizio Capuani, Matteo Mori, Andrea De Martino, Enzo Marinari

Abstract: Thermodynamics constrains the flow of matter in a reaction network to occur through routes along which the Gibbs energy decreases, implying that viable steady-state flux patterns should be void of closed reaction cycles. Identifying and removing cycles in large reaction networks can unfortunately be a highly challenging task from a computational viewpoint. We propose here a method that accomplishe… ▽ More Thermodynamics constrains the flow of matter in a reaction network to occur through routes along which the Gibbs energy decreases, implying that viable steady-state flux patterns should be void of closed reaction cycles. Identifying and removing cycles in large reaction networks can unfortunately be a highly challenging task from a computational viewpoint. We propose here a method that accomplishes it by combining a relaxation algorithm and a Monte Carlo procedure to detect loops, with ad hoc rules (discussed in detail) to eliminate them. As test cases, we tackle (a) the problem of identifying infeasible cycles in the E. coli metabolic network and (b) the problem of correcting thermodynamic infeasibilities in the Flux-Balance-Analysis solutions for 15 human cell-type specific metabolic networks. Results for (a) are compared with previous analyses of the same issue, while results for (b) are weighed against alternative methods to retrieve thermodynamically viable flux patterns based on minimizing specific global quantities. Our method on one hand outperforms previous techniques and, on the other, corrects loopy solutions to Flux Balance Analysis. As a byproduct, it also turns out to be able to reveal possible inconsistencies in model reconstructions. △ Less

Submitted 14 October, 2013; originally announced October 2013.

Comments: 10 pages; see http://chimera.roma1.infn.it/SYSBIO/ for supporting files

Journal ref: Metabolites 3:946 (2013)