Gallery up to 2017

My invited seminar at Università degli Studi di Milano – Segrate (November 16th 2017)

2017-11-16-Milano_v1-1

My invited seminar at Università degli Studi di Napoli Federico II (April 20th 2017)

2017-04-20-Napoli_v1

My presentation at MMM2016, Dijon (October 11st 2016)

2016-10-MMM

My presentation at STATPHYS26, Lyon (July 21st 2016)

2016-07-STATPHYS26

My presentation at the MECO 41 Conference, Vienna (February 15th 2016)

Slides from a seminar at the Theoretical Physics department in Universitat de Barcelona (13th July 2015)

My 4-works poster at PUMPS2015 Barcelona (best poster prize)

Fe-PUMPS2015

My presentation at Viscous Liquids IV, Montpellier, (5th May 2015)

My poster at the Avalanches, Intermittency, and Nonlinear Response in Far-From-Equilibrium Solids program in KITP Santa Barabara, California (November 2014)

FeMaBa-KITP2014

Parallel kinetic Monte Carlo simulation of Coulomb glasses

coulomb_gap — Pseudo-equilibrium (a) and time-dependent (b) single-particle density of states for a system size $N = 512 \times 512$ , as function of temperature and time respectively. Data in (b) corresponds to $T = 0.04$ .

Pseudo-equilibrium (a) and time-dependent (b) single-particle density of states for a system size $N = 512 \times 512$ , as function of temperature and time respectively. Data in (b) corresponds to $T = 0.04$ .

polarization_ESlaw — (a) Inset: time evolution of the polarization starting from a random initial configuration at different temperatures for $L = 512$ . The steady-state conductivity $\sigma_0 = \lim t\to\infty (dP/dt)/(NE)$ is reached after a transient time $\tau \simeq 1/(T^2\sigma_0)$ . Main figure: same data rescaled as a function of $t/\tau$ . (b) Fit to the Efros-Shklovskii law $\sigma_0 = bT^{-\lambda} \exp[−(T_0/T)^{1/2}]$ for $\lambda = 1,2$ .

performance — Computing time [ms] per executed hop for the serial and parallel codes in three different CPU+GPU platforms. We benchmark a fair serial single CPU core implementation against our parallel GPU implementation (with the same CPU core as host), both with double precision floating point operations. (a) Temperature dependence of the average rejection time for $L = 256$ . (b) Size dependence of the average update time of the $N = L^2$ local energies. The insets show the speedup of the GPU over the CPU implementation. While the speedup starts to saturate with L in (b), it is still strongly growing with $1/T$ in (a).

We develop a parallel rejection algorithm to tackle the problem of low acceptance in Monte Carlo methods, and apply it to the simulation of the hopping conduction in Coulomb glasses using Graphics Processing Units, for which we also parallelize the update of local energies. In two dimensions, our parallel code achieves speedups of up to two orders of magnitude in computing time over an equivalent serial code. We find numerical evidence of a scaling relation for the relaxation of the conductivity at different temperatures.

View full article…

Numerical approaches on driven elastic interfaces in random media

Linking transport and geometry. (a) Snapshot of a domain wall in a two-dimensional ferromagnet. (b) Typical velocity–force characteristics. (c) Crossover lengths l_opt and l_av representing the optimal excitation and the deterministic avalanches, respectively. (d) Geometric crossover diagram.

Low temperature dynamics of the driven elastic string below the depinning threshold.

We present the main concepts behind the statistical and dynamical properties of elastic systems in disordered media, focused on the relation between the rough geometry and collective transport properties in driven steady-states. We review the numerical approaches that allow to analyze the equilibrium, creep, and depinning regimes of motion in these models.

View full article…

Nonsteady relaxation and critical exponents at the depinning transition

This image has an empty alt attribute; its file name is FeBuKo2013_Fig7-300x191.png — The function $\gamma(t) = w(t) / t v(t)$ showing finite size effects in the QEW line relaxation.

The function $\gamma(t) = w(t) / t v(t)$ showing finite size effects in the QEW line relaxation.

This image has an empty alt attribute; its file name is FeBuKo2013_Fig15-300x175.png — Heuristic connection between steady state
and the nonsteady universal relaxation at the thermodynamical critical
force $f_c$

This image has an empty alt attribute; its file name is FeBuKo2013_Fig13-300x178.png — String velocity v(t) as a function of time, and rescaling with the critical exponents $\beta, \nu, z$ .

We study the nonsteady relaxation of a driven one-dimensional elastic interface at the depinning transition by extensive numerical simulations concurrently implemented on graphics processing units. Above a first, nonuniversal microscopic time regime, we find a nontrivial long crossover towards the nonsteady macroscopic critical regime. In order to avoid fitting effective exponents with a systematic bias we implement a practical criterion of consistency and perform large-scale ( $L\sim 2^{25}$ ) simulations for the nonsteady dynamics of the continuum displacement quenched Edwards-Wilkinson equation, getting accurate and consistent depinning exponents for this class: $\beta = 0.245 \pm 0.006$ , $z = 1.433 \pm 0.007$ , $\zeta = 1.250 \pm 0.005$ , and $\nu = 1.333 \pm 0.007$ .

View full article…

$q$ -state Potts model metastability study using optimized GPU-based Monte Carlo algorithms

This image has an empty alt attribute; its file name is FeThesis2011_Fig7.16-300x121.png — Schematic behavior for the internal energy as a function of the temperature at the $q$ -state Potts model first order phase transition. The two scenarios proposed by Binder.

This image has an empty alt attribute; its file name is FeDeWoCa2012_Fig2-300x186.png — Energy and magnetization as a function of temperature for different values of $q$ .

This image has an empty alt attribute; its file name is FeThesis2011_Fig7.18-300x203.png — Log–log plot of energy differences versus temperatures $T > T_t$

We implemented a GPU-based parallel code to perform Monte Carlo simulations of the two-dimensional $q$ -state Potts model. The possibility of performing high speed simulations at large enough system sizes allowed us to tackle an old standing question posed by k. Binder in 1981, regarding the existence or not of specific heat singularities at spinodal temperatures different of the transition one.

View full article and our code…
Also find here our approach to the spinodal points in the $q$ -state Potts model.

Relaxation in the $q$ -state Potts model

We studied the nonequilibrium dynamics of the $q$ -state Potts model in the square lattice, after a quench to subcritical temperatures. By means of a continuous time Monte Carlo algorithm ͑nonconserved order parameter dynamics͒ we analyzed the long term behavior of the energy and relaxation time for a wide range of quench temperatures and system sizes. For $q > 4$ we found the existence of different dynamical regimes, according to quench temperature range.