Source author record

Matteo Lulli

Matteo Lulli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.dis-nn cond-mat.soft cond-mat.stat-mech physics.flu-dyn nlin.CG physics.comp-ph gr-qc hep-th quant-ph

Catalog footprint

What is connected

10works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Mesoscale Modelling of the Tolman Length in Multi-component Systems

In this paper we analyze the curvature corrections to the surface tension in the context of the Shan-Chen (SC) multi-component Lattice Boltzmann method (LBM). We demonstrate that the same techniques recently applied in the context of the Shan-Chen multi-phase model can be applied to multi-component mixtures. We implement, as a new application, the calculation of the surface of tension radius $R_s$ through the minimization of the generalized surface tension $σ[R]$. In turn we are able to estimate the Tolman length, i.e. the first order coefficient of the curvature expansion of the surface tension $σ(R)$, as well as the higher order corrections, i.e. the curvature- and the Gaussian-rigidity coefficients. The SC multi-component model allows to model both fully-symmetric as well as asymmetric interactions among the components. By performing an extensive set of simulations we present a first example of tunable Tolman length in the mesoscopic model, being zero for symmetric interactions and different from zero otherwise. This result paves the way for controlling such interface properties which are paramount in presence of thermal fluctuations. All reported results can be independently reproduced through the "idea.deploy" framework available at https://github.com/lullimat/idea.deploy.

preprint2021arXiv

A Mesoscale Perspective on the Tolman Length

We demonstrate that the multi-phase Shan-Chen lattice Boltzmann method (LBM) yields a curvature dependent surface tension $σ$ as computed from three-dimensional hydrostatic droplets/bubbles simulations. Such curvature dependence is routinely characterized, at first order, by the so-called {\it Tolman length} $δ$. LBM allows to precisely compute $σ$ at the surface of tension $R_s$ and determine the Tolman length from the coefficient of the first order correction. The corresponding values of $δ$ display universality for different equations of state, following a power-law scaling near the critical temperature. The Tolman length has been studied so far mainly via computationally demanding molecular dynamics (MD) simulations or by means of density functional theory (DFT) approaches playing a pivotal role in extending Classical Nucleation Theory. The present results open a new hydrodynamic-compliant mesoscale arena, in which the fundamental role of the Tolman length, alongside real-world applications to cavitation phenomena, can be effectively tackled. All the results can be independently reproduced through the "idea.deploy" framework.

preprint2021arXiv

TLBfind: a Thermal Lattice Boltzmann code for concentrated emulsions with FINite-size Droplets

In this paper, we present TLBfind, a GPU code for simulating the hydrodynamics of droplets along with a dynamic temperature field. TLBfind hinges on a two-dimensional multi-component lattice Boltzmann (LB) model simulating a concentrated emulsion with finite-size droplets evolving in a thermal convective state, just above the transition from conduction to convection. The droplet concentration of the emulsion system is tenable, and at the core of the code lies the possibility to measure a large number of physical observables characterising the flow and droplets. Furthermore, TLBfind includes a parallel implementation on GPU of the Delaunay triangulation useful for the detection of droplets' plastic rearrangements, and several types of boundary conditions, supporting simulations of channels with structured rough walls.

preprint2020arXiv

Fragile Glasses Associated with a Dramatic Drop of Entropy under Supercooling

We perform kinetic Monte Carlo simulations of a distinguishable-particle lattice model of structural glasses with random particle interactions. By varying the interaction distribution and the average particle hopping energy barrier, we obtain an extraordinarily wide range of kinetic fragility. A stretching exponent, characterizing structural relaxation, is found to decrease with the kinetic fragility in agreement with experiments. The most fragile glasses are those exhibiting low hopping barriers and, more importantly, dramatic drops of entropies upon cooling toward the glass transition temperatures. The entropy drops reduce possible kinetic pathways and lead to dramatic slowdowns in the dynamics. In addition, the kinetic fragility is shown to correlate with a thermodynamic fragility.

preprint2020arXiv

Generalized Holographic Principle, Gauge Invariance and the Emergence of Gravity a la Wilczek

We show that a generalized version of the holographic principle can be derived from the Hamiltonian description of information flow within a quantum system that maintains a separable state. We then show that this generalized holographic principle entails a general principle of gauge invariance. When this is realized in an ambient Lorentzian space-time, gauge invariance under the Poincare group is immediately achieved. We apply this pathway to retrieve the action of gravity. The latter is cast a la Wilczek through a similar formulation derived by MacDowell and Mansouri, which involves the representation theory of the Lie groups SO(3,2) and SO(4,1).

preprint2019arXiv

Configuration-tree Theoretical Calculation of the Mean-Squared Displacement of Particles in Glass Formers

We report an analytical evaluation of the mean-squared displacement (MSD) of the particles in glasses based on their coarse grained trajectories. The calculation is conducted by means of a local random configuration-tree theory that was recently proposed by one of us [C.-H. Lam, J. Stat. Mech. \textbf{2018}, 023301 (2018)]. Results are compared with the numerical simulations of a lattice glass model, and good quantitative agreement has been obtained over a wide range of temperatures in the entire region of time with virtually no free parameters. To the best of our knowledge, the calculation is the first in its kind.

preprint2019arXiv

Spatial heterogeneities in structural temperature cause Kovacs expansion gap paradox in aging of glasses

Volume and enthalpy relaxation of glasses after a sudden temperature change has been extensively studied since Kovacs seminal work. One observes an asymmetric approach to equilibrium upon cooling versus heating and, more counter-intuitively, the expansion gap paradox, i.e. a dependence on the initial temperature of the effective relaxation time even close to equilibrium when heating. Here we show that a distinguishable-particles lattice model can capture both the asymmetry and the expansion gap. We quantitatively characterize the energetic states of the particles configurations using a physical realization of the fictive temperature called the structural temperature, which, in the heating case, displays a strong spatial heterogeneity. The system relaxes by nucleation and expansion of warmer mobile domains having attained the final temperature, against cooler immobile domains maintained at the initial temperature. A small population of these cooler regions persists close to equilibrium, thus explaining the paradox.

preprint2016arXiv

Fluidisation and plastic activity in a model soft-glassy material flowing in micro-channels with rough walls

By means of mesoscopic numerical simulations of a model soft-glassy material, we investigate the role of boundary roughness on the flow behaviour of the material, probing the bulk/wall and global/local rheologies. We show that the roughness reduces the wall slip induced by wettability properties and acts as a source of fluidisation for the material. A direct inspection of the plastic events suggests that their rate of occurrence grows with the fluidity field, reconciling our simulations with kinetic elasto-plastic descriptions of jammed materials. Notwithstanding, we observe qualitative and quantitative differences in the scaling, depending on the distance from the rough wall and on the imposed shear. The impact of roughness on the orientational statistics is also studied.

preprint2015arXiv

Off-equilibrium finite-size method for critical behavior analyses

We present a new dynamic off-equilibrium method for the study of continuous transitions, which represents a dynamic generalization of the usual equilibrium cumulant method. Its main advantage is that critical parameters are derived from numerical data obtained much before equilibrium has been attained. Therefore, the method is particularly useful for systems with long equilibration times, like spin glasses. We apply it to the three-dimensional Ising spin-glass model, obtaining accurate estimates of the critical exponents and of the critical temperature with a limited computational effort.

preprint2014arXiv

Highly optimized simulations on single- and multi-GPU systems of 3D Ising spin glass

We present a highly optimized implementation of a Monte Carlo (MC) simulator for the three-dimensional Ising spin-glass model with bimodal disorder, i.e., the 3D Edwards-Anderson model running on CUDA enabled GPUs. Multi-GPU systems exchange data by means of the Message Passing Interface (MPI). The chosen MC dynamics is the classic Metropolis one, which is purely dissipative, since the aim was the study of the critical off-equilibrium relaxation of the system. We focused on the following issues: i) the implementation of efficient access patterns for nearest neighbours in a cubic stencil and for lagged-Fibonacci-like pseudo-Random Numbers Generators (PRNGs); ii) a novel implementation of the asynchronous multispin-coding Metropolis MC step allowing to store one spin per bit and iii) a multi-GPU version based on a combination of MPI and CUDA streams. We highlight how cubic stencils and PRNGs are two subjects of very general interest because of their widespread use in many simulation codes. Our code best performances ~3 and ~5 psFlip on a GTX Titan with our implementations of the MINSTD and MT19937 respectively.