Showing 1–37 of 37 results for author: Buttà, P

Asymptotics of the $φ^4_1$ measure in the sharp interface limit

Authors: Lorenzo Bertini, Paolo Buttà, Giacomo Di Gesù

Abstract: We consider the $φ^4_1$ measure in an interval of length $\ell$, defined by a symmetric double-well potential $W$ and inverse temperature $β$. Our results concern its asymptotic behavior in the joint limit $β, \ell \to \infty$, both in the subcritical regime $\ell \ll \rme^{βC_W}$ and in the supercritical regime $\ell \gg \rme^{βC_W}$, where $C_W$ denotes the surface tension. In the former case, i… ▽ More We consider the $φ^4_1$ measure in an interval of length $\ell$, defined by a symmetric double-well potential $W$ and inverse temperature $β$. Our results concern its asymptotic behavior in the joint limit $β, \ell \to \infty$, both in the subcritical regime $\ell \ll \rme^{βC_W}$ and in the supercritical regime $\ell \gg \rme^{βC_W}$, where $C_W$ denotes the surface tension. In the former case, in which the measure concentrates on the pure phases, we prove the corresponding large deviation principle. The associated rate function is the Modica-Mortola functional modified to take into account the entropy of the locations of the interfaces. Further, we provide the sharp asymptotics of the probability of having a given number of transitions between the two pure phases. In the supercritical regime, the measure does not longer concentrate and we show that the interfaces are asymptotically distributed according to a Poisson point process. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 38 pages

Report number: Roma01.Math.MP MSC Class: 82B24; 81Q20; 60F10; 60J60

Long time evolution of concentrated vortex rings with large radius

Authors: Paolo Buttà, Guido Cavallaro, Carlo Marchioro

Abstract: We study the time evolution of an incompressible fluid with axial symmetry without swirl when the vorticity is sharply concentrated on $N$ annuli of radii of the order of $r_0$ and thickness $\varepsilon$. We prove that when $r_0= |\log \varepsilon|^α$, $α>1$, the vorticity field of the fluid converges for $\varepsilon \to 0$ to the point vortex model, in an interval of time which diverges as… ▽ More We study the time evolution of an incompressible fluid with axial symmetry without swirl when the vorticity is sharply concentrated on $N$ annuli of radii of the order of $r_0$ and thickness $\varepsilon$. We prove that when $r_0= |\log \varepsilon|^α$, $α>1$, the vorticity field of the fluid converges for $\varepsilon \to 0$ to the point vortex model, in an interval of time which diverges as $\log|\log\varepsilon|$. This generalizes previous result by Cavallaro and Marchioro in [J. Math. Phys. 62, 053102, (2021)], that assumed $α>2$ and in which the convergence was proved for short times only. △ Less

Submitted 10 May, 2024; originally announced May 2024.

Comments: 24 pages

Report number: Roma01.Math.MP MSC Class: 76B47; 37N10

Leapfrogging vortex rings as scaling limit of Euler Equations

Authors: Paolo Buttà, Guido Cavallaro, Carlo Marchioro

Abstract: We consider an incompressible fluid with axial symmetry without swirl, assuming initial data such that the initial vorticity is very concentrated inside $N$ small disjoint rings of thickness $\varepsilon$, each one of vorticity mass and main radius of order $|\log\varepsilon|$. When $\varepsilon \to 0$, we show that, at least for small but positive times, the motion of the rings converges to a dyn… ▽ More We consider an incompressible fluid with axial symmetry without swirl, assuming initial data such that the initial vorticity is very concentrated inside $N$ small disjoint rings of thickness $\varepsilon$, each one of vorticity mass and main radius of order $|\log\varepsilon|$. When $\varepsilon \to 0$, we show that, at least for small but positive times, the motion of the rings converges to a dynamical system firstly introduced in [NoDEA Nonlinear Diff. Eq. Appl. 6 (1999), 473-499]. In the special case of two vortex rings with large enough main radius, the result is improved reaching longer times, in such a way to cover the case of several overtakings between the rings, thus providing a mathematical rigorous derivation of the leapfrogging phenomenon. △ Less

Submitted 17 August, 2024; v1 submitted 1 October, 2023; originally announced October 2023.

Comments: 33 pages, 1 figure; typos corrected, references updated

Report number: Roma01.Math.MP MSC Class: 76B47; 37N10

From particle systems to the BGK equation

Authors: Paolo Buttà, Mario Pulvirenti, Sergio Simonella

Abstract: In [Phys. Rev. 94 (1954), 511-525], P.L. Bhatnagar, E.P. Gross and M. Krook introduced a kinetic equation (the BGK equation), effective in physical situations where the Knudsen number is small compared to the scales where Boltzmann's equation can be applied, but not enough for using hydrodynamic equations. In this paper, we consider the stochastic particle system (inhomogeneous Kac model) underlyi… ▽ More In [Phys. Rev. 94 (1954), 511-525], P.L. Bhatnagar, E.P. Gross and M. Krook introduced a kinetic equation (the BGK equation), effective in physical situations where the Knudsen number is small compared to the scales where Boltzmann's equation can be applied, but not enough for using hydrodynamic equations. In this paper, we consider the stochastic particle system (inhomogeneous Kac model) underlying Bird's direct simulation Monte Carlo method (DSMC), with tuning of the scaled variables yielding kinetic and/or hydrodynamic descriptions. Although the BGK equation cannot be obtained from pure scaling, it does follow from a simple modification of the dynamics. This is proposed as a mathematical interpretation of some arguments in [Phys. Rev. 94 (1954), 511-525], complementing previous results in [Arch. Ration. Mech. Anal. 240 (2021), 785-808] and [Kinet. Relat. Models 16 (2023), 269-293]. △ Less

Submitted 27 January, 2023; originally announced January 2023.

Comments: 10 pages

Report number: Roma01.Math.MP MSC Class: 82C40 (Primary) 60J75; 82C22 (Secondary)

Journal ref: Ensaios Matemáticos 38 (2023), pp. 127-140

Vanishing viscosity limit for concentrated vortex rings

Authors: Paolo Buttà, Guido Cavallaro, Carlo Marchioro

Abstract: We study the time evolution of a viscous incompressible fluid with axial symmetry without swirl, when the initial vorticity is very concentrated in $N$ disjoint rings. We show that in a suitable joint limit, in which both the thickness of the rings and the viscosity tend to zero, the vorticity remains concentrated in $N$ disjointed rings, each one of them performing a simple translation along the… ▽ More We study the time evolution of a viscous incompressible fluid with axial symmetry without swirl, when the initial vorticity is very concentrated in $N$ disjoint rings. We show that in a suitable joint limit, in which both the thickness of the rings and the viscosity tend to zero, the vorticity remains concentrated in $N$ disjointed rings, each one of them performing a simple translation along the symmetry axis with constant speed. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: 29 pages, LaTex. This article is an extension of the results in arXiv:1904.04785 and arXiv:2102.07807 and therefore draws from these articles

Report number: Roma01.Math.MP MSC Class: 76D17; 37N10

Journal ref: J. Math. Phys. 63 (2022), 123103

Non-mean-field Vicsek-type models for collective behaviour

Authors: P. Buttà, B. Goddard, T. M. Hodgson, M. Ottobre, K. J. Painter

Abstract: We consider interacting particle dynamics with Vicsek type interactions, and their macroscopic PDE limit, in the non-mean-field regime; that is, we consider the case in which each particle/agent in the system interacts only with a prescribed subset of the particles in the system (for example, those within a certain distance). In this non-mean-field regime the influence between agents (i.e. the int… ▽ More We consider interacting particle dynamics with Vicsek type interactions, and their macroscopic PDE limit, in the non-mean-field regime; that is, we consider the case in which each particle/agent in the system interacts only with a prescribed subset of the particles in the system (for example, those within a certain distance). In this non-mean-field regime the influence between agents (i.e. the interaction term) can be normalised either by the total number of agents in the system (\textit{global scaling}) or by the number of agents with which the particle is effectively interacting (\textit{local scaling}). We compare the behaviour of the globally scaled and the locally scaled systems in many respects, considering for each scaling both the PDE and the corresponding particle model. In particular we observe that both the locally and globally scaled particle system exhibit pattern formation (i.e. formation of travelling-wave-like solutions) within certain parameter regimes, and generally display similar dynamics. The same is not true of the corresponding PDE models. Indeed, while both PDE models have multiple stationary states, for the globally scaled PDE such (space-homogeneous) equilibria are unstable for certain parameter regimes, with the instability leading to travelling wave solutions, while they are always stable for the locally scaled one, which never produces travelling waves. This observation is based on a careful numerical study of the model, supported by further analysis. △ Less

Submitted 17 May, 2022; originally announced June 2022.

Comments: 44 pages

A stochastic particle system approximating the BGK equation

Authors: Paolo Buttà, Mario Pulvirenti

Abstract: We consider a stochastic $N$-particle system on a torus in which each particle moving freely can instantaneously thermalize according to the particle configuration at that instant. Following [2], we show that the propagation of chaos does hold and that the one-particle distribution converges to the solution of the BGK equation. The improvement with respect to [2] consists in the fact that here, as… ▽ More We consider a stochastic $N$-particle system on a torus in which each particle moving freely can instantaneously thermalize according to the particle configuration at that instant. Following [2], we show that the propagation of chaos does hold and that the one-particle distribution converges to the solution of the BGK equation. The improvement with respect to [2] consists in the fact that here, as suggested by physical considerations, the thermalizing transition is driven only by the restriction of the particle configuration in a small neighborhood of the jumping particle. In other words, the Maxwellian distribution of the outgoing particle is computed via the empirical hydrodynamical fields associated to the fraction of particles sufficiently close to the test particle and not, as in [2], via the whole particle configuration. [2]: arXiv:2002.10535 (Journal reference: Arch. Ration. Mech. Anal. Vol. 240 (2021), pp. 785-808) △ Less

Submitted 5 January, 2022; originally announced January 2022.

Comments: 25 pages. This article is an extension of the results in arXiv: 2002.10535 and draws some insights from it

Report number: Roma01.Math.MP MSC Class: 82C40 (Primary) 60K35; 82C22 (Secondary)

Journal ref: Kinet. Relat. Models 16 (2023), pp. 269-293

Global time evolution of concentrated vortex rings

Authors: Paolo Buttà, Guido Cavallaro, Carlo Marchioro

Abstract: We study the time evolution of an incompressible fluid with axial symmetry without swirl, assuming initial data such that the initial vorticity is very concentrated inside $N$ small disjoint rings of thickness $\varepsilon$ and vorticity mass of the order of $|\log\varepsilon|^{ -1}$. When $\varepsilon \to 0$ we show that the motion of each vortex ring converges to a simple translation with consta… ▽ More We study the time evolution of an incompressible fluid with axial symmetry without swirl, assuming initial data such that the initial vorticity is very concentrated inside $N$ small disjoint rings of thickness $\varepsilon$ and vorticity mass of the order of $|\log\varepsilon|^{ -1}$. When $\varepsilon \to 0$ we show that the motion of each vortex ring converges to a simple translation with constant speed (depending on the single ring) along the symmetry axis. We obtain a sharp localization of the vorticity support at time $t$ in the radial direction, whereas we state only a concentration property in the axial direction. This is obtained for arbitrary (but fixed) intervals of time. This study is the completion of a previous paper arXiv:1904.04785 [math-ph], where a sharp localization of the vorticity support was obtained both along the radial and axial directions, but the convergence for $\varepsilon\to 0$ worked only for short times. △ Less

Submitted 5 February, 2022; v1 submitted 15 February, 2021; originally announced February 2021.

Comments: 23 pages, LaTex; small changes in abstract and introduction, typos corrected, added an appendix with the proof of Lemma 4.1, and, with respect to version 2, Lemma 3.2 corrected (which implies very slight changes in other parts of the manuscript). This article is an extension of the results arXiv:1904.04785 and therefore draws heavily from arXiv:1904.04785

Report number: Roma01.Math.MP MSC Class: 76B47; 37N10

Journal ref: Z. Angew. Math. Phys. 73 (2022), Article 70

Cucker-Smale type dynamics of infinitely many individuals with repulsive forces

Authors: Paolo Buttà, Carlo Marchioro

Abstract: We study the existence and uniqueness of the time evolution of a system of infinitely many individuals, moving in a tunnel and subjected to a Cucker-Smale type alignment dynamics with compactly supported communication kernels and to short-range repulsive interactions to avoid collisions. We study the existence and uniqueness of the time evolution of a system of infinitely many individuals, moving in a tunnel and subjected to a Cucker-Smale type alignment dynamics with compactly supported communication kernels and to short-range repulsive interactions to avoid collisions. △ Less

Submitted 13 March, 2020; originally announced March 2020.

Comments: 13 pages

Report number: Roma01.Math.MP MSC Class: 82C22 (Primary) 92D25; 37N25 (Secondary)

Journal ref: J. Stat. Phys. 181 (2020), pp. 2094-2108

Stability of the stationary solutions of the Allen-Cahn equation with non-constant stiffness

Authors: Paolo Buttà, Emilio N. M. Cirillo, Giulio Sciarra

Abstract: We study the solutions of a generalized Allen-Cahn equation deduced from a Landau energy functional, endowed with a non-constant higher order stiffness. We assume the stiffness to be a positive function of the field and we discuss the stability of the stationary solutions proving both linear and local non-linear stability. We study the solutions of a generalized Allen-Cahn equation deduced from a Landau energy functional, endowed with a non-constant higher order stiffness. We assume the stiffness to be a positive function of the field and we discuss the stability of the stationary solutions proving both linear and local non-linear stability. △ Less

Submitted 31 July, 2020; v1 submitted 1 March, 2020; originally announced March 2020.

Report number: Roma01.Math.MP

Journal ref: Wave Motion 98, 102641, 2020

Particle approximation of the BGK equation

Authors: Paolo Buttà, Maxime Hauray, Mario Pulvirenti

Abstract: In this paper we prove the convergence of a suitable particle system towards the BGK model. More precisely, we consider an interacting stochastic particle system in which each particle can instantaneously thermalize locally. We show that, under a suitable scaling limit, propagation of chaos does hold and that the one-particle distribution function converges to the solution of the BGK equation. In this paper we prove the convergence of a suitable particle system towards the BGK model. More precisely, we consider an interacting stochastic particle system in which each particle can instantaneously thermalize locally. We show that, under a suitable scaling limit, propagation of chaos does hold and that the one-particle distribution function converges to the solution of the BGK equation. △ Less

Submitted 23 November, 2020; v1 submitted 24 February, 2020; originally announced February 2020.

Comments: 21 pages

Report number: Roma01.Math.MP MSC Class: 82C40 (Primary) 60K35; 82C22 (Secondary)

Journal ref: Arch. Ration. Mech. Anal. Vol. 240 (2021), pp. 785-808

Some aspects of the inertial spin model for flocks and related kinetic equations

Authors: Dario Benedetto, Paolo Buttà, Emanuele Caglioti

Abstract: In this paper we study the macroscopic behavior of the inertial spin (IS) model. This model has been recently proposed to describe the collective dynamics of flocks of birds, and its main feature is the presence of an auxiliary dynamical variable, a sort of internal spin, which conveys the interaction among the birds with the effect of better describing the turning of flocks. After discussing the… ▽ More In this paper we study the macroscopic behavior of the inertial spin (IS) model. This model has been recently proposed to describe the collective dynamics of flocks of birds, and its main feature is the presence of an auxiliary dynamical variable, a sort of internal spin, which conveys the interaction among the birds with the effect of better describing the turning of flocks. After discussing the geometrical and mechanical properties of the IS model, we show that, in the case of constant interaction among the birds, its mean-field limit is described by a non-linear Fokker-Planck equation, whose equilibria are fully characterized. Finally, in the case of non-constant interactions, we derive the kinetic equation for the mean-field limit of the model in the absence of thermal noise, and explore its macroscopic behavior by analyzing the mono-kinetic solutions. △ Less

Submitted 6 November, 2019; originally announced November 2019.

Comments: 29 pages

Report number: Roma01.Math.MP MSC Class: 60K35; 82B40; 92D50

Journal ref: Math. Models Methods Appl. Sci. 30 (2020), pp. 1987-2022

Time evolution of concentrated vortex rings

Authors: Paolo Buttà, Carlo Marchioro

Abstract: We study the time evolution of an incompressible fluid with axisymmetry without swirl when the vorticity is sharply concentrated. In particular, we consider $N$ disjoint vortex rings of size $\varepsilon$ and intensity of the order of $|\log\varepsilon|^{-1}$. We show that in the limit $\varepsilon\to 0$, when the density of vorticity becomes very large, the movement of each vortex ring converges… ▽ More We study the time evolution of an incompressible fluid with axisymmetry without swirl when the vorticity is sharply concentrated. In particular, we consider $N$ disjoint vortex rings of size $\varepsilon$ and intensity of the order of $|\log\varepsilon|^{-1}$. We show that in the limit $\varepsilon\to 0$, when the density of vorticity becomes very large, the movement of each vortex ring converges to a simple translation, at least for a small but positive time. △ Less

Submitted 8 June, 2021; v1 submitted 9 April, 2019; originally announced April 2019.

Comments: 24 pages. This updated version provides a new Appendix B, containing the corrected proof of Lemma 3.1. For the sake of clarity, this proof has already been included in arXiv:2102.07807 (where the results of this article have been extended)

Report number: Roma01.Math.MP MSC Class: 76B47; 37N10

Journal ref: J. Math. Fluid Mech. 22 (2020), Article 9

A new and stable estimation method of country economic fitness and product complexity

Authors: Vito D. P. Servedio, Paolo Buttà, Dario Mazzilli, Andrea Tacchella, Luciano Pietronero

Abstract: We present a new metric estimating fitness of countries and complexity of products by exploiting a non-linear non-homogeneous map applied to the publicly available information on the goods exported by a country. The non homogeneous terms guarantee both convergence and stability. After a suitable rescaling of the relevant quantities, the non homogeneous terms are eventually set to zero so that this… ▽ More We present a new metric estimating fitness of countries and complexity of products by exploiting a non-linear non-homogeneous map applied to the publicly available information on the goods exported by a country. The non homogeneous terms guarantee both convergence and stability. After a suitable rescaling of the relevant quantities, the non homogeneous terms are eventually set to zero so that this new metric is parameter free. This new map almost reproduces the results of the original homogeneous metrics already defined in literature and allows for an approximate analytic solution in case of actual binarized matrices based on the Revealed Comparative Advantage (RCA) indicator. This solution is connected with a new quantity describing the neighborhood of nodes in bipartite graphs, representing in this work the relations between countries and exported products. Moreover, we define the new indicator of country net-efficiency quantifying how a country efficiently invests in capabilities able to generate innovative complex high quality products. Eventually, we demonstrate analytically the local convergence of the algorithm involved. △ Less

Submitted 12 October, 2018; v1 submitted 23 July, 2018; originally announced July 2018.

Comments: 12 pages, 8 figures

Journal ref: Entropy 2018, 20(10), 783

A non-linear kinetic model of self-propelled particles with multiple equilibria

Authors: Paolo Buttà, Franco Flandoli, Michela Ottobre, Boguslaw Zegarlinski

Abstract: We introduce and analyse a continuum model for an interacting particle system of Vicsek type. The model is given by a non-linear kinetic partial differential equation (PDE) describing the time-evolution of the density $f_t$, in the single particle phase-space, of a collection of interacting particles confined to move on the one-dimensional torus. The corresponding stochastic differential equation… ▽ More We introduce and analyse a continuum model for an interacting particle system of Vicsek type. The model is given by a non-linear kinetic partial differential equation (PDE) describing the time-evolution of the density $f_t$, in the single particle phase-space, of a collection of interacting particles confined to move on the one-dimensional torus. The corresponding stochastic differential equation for the position and velocity of the particles is a conditional McKean-Vlasov type of evolution (conditional in the sense that the process depends on its own law through its own conditional expectation). In this paper, we study existence and uniqueness of the solution of the PDE in consideration. Challenges arise from the fact that the PDE is neither elliptic (the linear part is only {\em hypoelliptic}) nor in gradient form. Moreover, for some specific choices of the interaction function and for the simplified case in which the density profile does not depend on the spatial variable, we show that the model exhibits multiple stationary states (corresponding to the particles forming a coordinated clockwise/anticlockwise rotational motion) and we study convergence to such states as well. Finally, we prove mean-field convergence of an appropriate $N$-particles system to the solution of our PDE: more precisely, we show that the empirical measures of such a particle system converge weakly, as $N \rightarrow \infty$, to the solution of the PDE. △ Less

Submitted 16 January, 2019; v1 submitted 4 April, 2018; originally announced April 2018.

Comments: 37 pages

Report number: Roma01.Math.MP MSC Class: 35A01; 35B40; 82C22; 82C40; 92D25

Journal ref: Kinetic and Relat. Models Vol. 12 (2019), pp. 791-827

On large deviations of interface motions for statistical mechanics models

Authors: Lorenzo Bertini, Paolo Buttà, Adriano Pisante

Abstract: We discuss the sharp interface limit of the action functional associated to either the Glauber dynamics for Ising systems with Kac potentials or the Glauber+Kawasaki process. The corresponding limiting functionals, for which we provide explicit formulae of the mobility and transport coefficients, describe the large deviations asymptotics with respect to the mean curvature flow. We discuss the sharp interface limit of the action functional associated to either the Glauber dynamics for Ising systems with Kac potentials or the Glauber+Kawasaki process. The corresponding limiting functionals, for which we provide explicit formulae of the mobility and transport coefficients, describe the large deviations asymptotics with respect to the mean curvature flow. △ Less

Submitted 12 February, 2018; originally announced February 2018.

Comments: 32 pages

Report number: Roma01.Math.MP MSC Class: 82C24; 60F10; 53C44

Journal ref: Ann. Henri Poincaré 20 (2019), pp. 1785-1821

Long time evolution of concentrated Euler flows with planar symmetry

Authors: Paolo Buttà, Carlo Marchioro

Abstract: We study the time evolution of an incompressible Euler fluid with planar symmetry when the vorticity is initially concentrated in small disks. We discuss how long this concentration persists, showing that in some cases this happens for quite long times. Moreover, we analyze a toy model that shows a similar behavior and gives some hints on the original problem. We study the time evolution of an incompressible Euler fluid with planar symmetry when the vorticity is initially concentrated in small disks. We discuss how long this concentration persists, showing that in some cases this happens for quite long times. Moreover, we analyze a toy model that shows a similar behavior and gives some hints on the original problem. △ Less

Submitted 21 January, 2018; v1 submitted 15 November, 2016; originally announced November 2016.

Comments: 26 pages

Report number: Roma01.Math.MP MSC Class: 76B47; 37N10; 70K65

Journal ref: SIAM J. Math. Anal. 50 (2018), pp. 735-760

Stochastic Allen-Cahn approximation of the mean curvature flow: large deviations upper bound

Authors: Lorenzo Bertini, Paolo Buttà, Adriano Pisante

Abstract: Consider the Allen-Cahn equation on the $d$-dimensional torus, $d=2,3$, in the sharp interface limit. As it is well known, the limiting dynamics is described by the motion by mean curvature of the interface between the two stable phases. Here, we analyze a stochastic perturbation of the Allen-Cahn equation and describe its large deviation asymptotics in a joint sharp interface and small noise limi… ▽ More Consider the Allen-Cahn equation on the $d$-dimensional torus, $d=2,3$, in the sharp interface limit. As it is well known, the limiting dynamics is described by the motion by mean curvature of the interface between the two stable phases. Here, we analyze a stochastic perturbation of the Allen-Cahn equation and describe its large deviation asymptotics in a joint sharp interface and small noise limit. Relying on previous results on the variational convergence of the action functional, we prove the large deviation upper bound. The corresponding rate function is finite only when there exists a time evolving interface of codimension one between the two stable phases. The zero level set of this rate function is given by the evolution by mean curvature in the sense of Brakke. Finally, the rate function can be written in terms of the sum of two non-negative quantities: the first measures how much the velocity of the interface deviates from its mean curvature, while the second is due to the possible occurrence of nucleation events. △ Less

Submitted 16 February, 2017; v1 submitted 7 April, 2016; originally announced April 2016.

Comments: 41 pages

Report number: Roma01.Math.MP MSC Class: 35R60; 60F10; 82C24; 53C44

Journal ref: Arch. Ration. Mech. Anal. Vol. 224 (2017), pp. 659-707

Dynamics of infinite classical anharmonic crystals

Authors: Paolo Buttà, Carlo Marchioro

Abstract: We consider an unbounded lattice and at each point of this lattice an anharmonic oscillator, that interacts with its first neighborhoods via a pair potential $V$ and is subjected to a restoring force of potential $U$. We assume that $U$ and $V$ are even nonnegative polynomials of degree $2σ_1$ and $2σ_2$. We study the time evolution of this system, with a control of the growth in time of the local… ▽ More We consider an unbounded lattice and at each point of this lattice an anharmonic oscillator, that interacts with its first neighborhoods via a pair potential $V$ and is subjected to a restoring force of potential $U$. We assume that $U$ and $V$ are even nonnegative polynomials of degree $2σ_1$ and $2σ_2$. We study the time evolution of this system, with a control of the growth in time of the local energy, and we give a nontrivial bound on the velocity of propagation of a perturbation. This is an extension to the case $σ_1 < 2σ_2-1$ of some already known results obtained for $σ_1 \geq 2σ_2-1$. △ Less

Submitted 3 February, 2016; originally announced February 2016.

Comments: 13 pages

Report number: Roma01.Math.MP

Journal ref: J. Stat. Phys. 164 (2016), pp. 680-692

Stochastic Allen-Cahn equation with mobility

Authors: Lorenzo Bertini, Paolo Buttà, Adriano Pisante

Abstract: We introduce a class of stochastic Allen-Cahn equations with a mobility coefficient and colored noise. For initial data with finite free energy, we analyze the corresponding Cauchy problem on the $d$-dimensional torus in the time interval $[0,T]$. Assuming that $d\le 3$ and that the potential has quartic growth, we prove existence and uniqueness of the solution as a process $u$ in $L^2$ with conti… ▽ More We introduce a class of stochastic Allen-Cahn equations with a mobility coefficient and colored noise. For initial data with finite free energy, we analyze the corresponding Cauchy problem on the $d$-dimensional torus in the time interval $[0,T]$. Assuming that $d\le 3$ and that the potential has quartic growth, we prove existence and uniqueness of the solution as a process $u$ in $L^2$ with continuous paths, satisfying almost surely the regularity properties $u\in C([0,T]; H^1)$ and $u\in L^2([0,T];H^2)$. △ Less

Submitted 28 July, 2017; v1 submitted 29 December, 2015; originally announced December 2015.

Comments: 34 pages

Report number: Roma01.Math.AP MSC Class: 35R60; 60H15

Journal ref: Nonlinear Differ. Equ. Appl. 24: 54 (2017)

Differential equations for real-structured (and unstructured) defectivity measures

Authors: Paolo Buttà, Nicola Guglielmi, Manuela Manetta, Silvia Noschese

Abstract: Let $A$ be either a complex or real matrix with all distinct eigenvalues. We propose a new method for the computation of both the unstructured and the real-structured (if the matrix is real) distance $w_{\mathbb K}(A)$ (where ${\mathbb K}=\mathbb C$ if general complex matrices are considered and ${\mathbb K} ={\mathbb R}$ if only real matrices are allowed) of the matrix $A$ from the set of defecti… ▽ More Let $A$ be either a complex or real matrix with all distinct eigenvalues. We propose a new method for the computation of both the unstructured and the real-structured (if the matrix is real) distance $w_{\mathbb K}(A)$ (where ${\mathbb K}=\mathbb C$ if general complex matrices are considered and ${\mathbb K} ={\mathbb R}$ if only real matrices are allowed) of the matrix $A$ from the set of defective matrices, that is the set of those matrices with at least a multiple eigenvalue with algebraic multiplicity larger than its geometric multiplicity. For $0 < \varepsilon \le w_{\mathbb K}(A)$, this problem is closely related to the computation of the most ill-conditioned $\varepsilon$-pseudoeigenvalues of $A$, that is points in the $\varepsilon$-pseudospectrum of $A$ characterized by the highest condition number. The method we propose couples a system of differential equations on a low rank (possibly structured) manifold which computes the $\varepsilon$-pseudoeigenvalue of $A$ which is closest to coalesce, with a fast Newton-like iteration aiming to determine the minimal value $\varepsilon$ such that such an $\varepsilon$-pseudoeigenvalue becomes defective. The method has a local behaviour; this means that in general we find upper bounds for $w_{\mathbb K}(A)$. However, they usually provide good approximations, in those (simple) cases where we can check this. The methodology can be extended to a structured matrix where it is required that the distance is computed within some manifold defining the structure of the matrix. In this paper we extensively examine the case of real matrices but we also consider pattern structures. As far as we know there do not exist methods in the literature able to compute such distance. △ Less

Submitted 14 April, 2014; originally announced April 2014.

Comments: 30 pages

Report number: Roma01.Math.NA MSC Class: 15A18; 65K05

Journal ref: SIAM J. Matrix Anal. Appl. 36 (2015), 523-548

Front fluctuations for the stochastic Cahn-Hilliard equation

Authors: L. Bertini, S. Brassesco, P. Buttà

Abstract: We consider the Cahn-Hilliard equation in one space dimension, perturbed by the derivative of a space and time white noise of intensity $ε^{\frac 12}$, and we investigate the effect of the noise, as $ε\to 0$, on the solutions when the initial condition is a front that separates the two stable phases. We prove that, given $γ< \frac 23$, with probability going to one as $ε\to 0$, the solution remain… ▽ More We consider the Cahn-Hilliard equation in one space dimension, perturbed by the derivative of a space and time white noise of intensity $ε^{\frac 12}$, and we investigate the effect of the noise, as $ε\to 0$, on the solutions when the initial condition is a front that separates the two stable phases. We prove that, given $γ< \frac 23$, with probability going to one as $ε\to 0$, the solution remains close to a front for times of the order of $ε^{-γ}$, and we study the fluctuations of the front in this time scaling. They are given by a one dimensional continuous process, self similar of order $\frac 14$ and non Markovian, related to a fractional Brownian motion and for which a couple of representations are given. △ Less

Submitted 7 March, 2014; originally announced March 2014.

Comments: 33 pages

Report number: Roma01.Math.MP MSC Class: 60H15; 35R60; 82C24

Journal ref: Braz. J. Probab. Stat. 29 (2015), 336-371

Drift of phase fluctuations in the ABC model

Authors: L. Bertini, P. Buttà

Abstract: In a recent work, Bodineau and Derrida analyzed the phase fluctuations in the ABC model. In particular, they computed the asymptotic variance and, on the basis of numerical simulations, they conjectured the presence of a drift, which they guessed to be an antisymmetric function of the three densities. By assuming the validity of the fluctuating hydrodynamic approximation, we prove the presence of… ▽ More In a recent work, Bodineau and Derrida analyzed the phase fluctuations in the ABC model. In particular, they computed the asymptotic variance and, on the basis of numerical simulations, they conjectured the presence of a drift, which they guessed to be an antisymmetric function of the three densities. By assuming the validity of the fluctuating hydrodynamic approximation, we prove the presence of such a drift, providing an analytical expression for it. This expression is then shown to be an antisymmetric function of the three densities. The antisymmetry of the drift can also be inferred from a symmetry property of the underlying microscopic dynamics. △ Less

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

Report number: Roma01.Math.MP

Journal ref: J. Stat. Phys. 152 (2013), pp. 15-36

Froth-like minimizers of a non local free energy functional with competing interactions

Authors: Paolo Buttà, Raffaele Esposito, Alessandro Giuliani, Rossana Marra

Abstract: We investigate the ground and low energy states of a one dimensional non local free energy functional describing at a mean field level a spin system with both ferromagnetic and antiferromagnetic interactions. In particular, the antiferromagnetic interaction is assumed to have a range much larger than the ferromagnetic one. The competition between these two effects is expected to lead to the sponta… ▽ More We investigate the ground and low energy states of a one dimensional non local free energy functional describing at a mean field level a spin system with both ferromagnetic and antiferromagnetic interactions. In particular, the antiferromagnetic interaction is assumed to have a range much larger than the ferromagnetic one. The competition between these two effects is expected to lead to the spontaneous emergence of a regular alternation of long intervals on which the spin profile is magnetized either up or down, with an oscillation scale intermediate between the range of the ferromagnetic and that of the antiferromagnetic interaction. In this sense, the optimal or quasi-optimal profiles are "froth-like": if seen on the scale of the antiferromagnetic potential they look neutral, but if seen at the microscope they actually consist of big bubbles of two different phases alternating among each other. In this paper we prove the validity of this picture, we compute the oscillation scale of the quasi-optimal profiles and we quantify their distance in norm from a reference periodic profile. The proof consists of two main steps: we first coarse grain the system on a scale intermediate between the range of the ferromagnetic potential and the expected optimal oscillation scale; in this way we reduce the original functional to an effective "sharp interface" one. Next, we study the latter by reflection positivity methods, which require as a key ingredient the exact locality of the short range term. Our proof has the conceptual interest of combining coarse graining with reflection positivity methods, an idea that is presumably useful in much more general contexts than the one studied here. △ Less

Submitted 18 June, 2012; originally announced June 2012.

Comments: 38 pages, 2 figures

Report number: Roma01.Math.MP

Journal ref: Commun. Math. Phys. 322 (2013), 593-632

Structured conditioning of Hamiltonian eigenvalue problems

Authors: Paolo Buttà, Silvia Noschese

Abstract: We discuss the effect of structure-preserving perturbations on complex or real Hamiltonian eigenproblems and characterize the structured worst-case effect perturbations. We derive significant expressions for both the structured condition numbers and the worst-case effect Hamiltonian perturbations. It is shown that, for purely imaginary eigenvalues, the usual unstructured perturbation analysis is s… ▽ More We discuss the effect of structure-preserving perturbations on complex or real Hamiltonian eigenproblems and characterize the structured worst-case effect perturbations. We derive significant expressions for both the structured condition numbers and the worst-case effect Hamiltonian perturbations. It is shown that, for purely imaginary eigenvalues, the usual unstructured perturbation analysis is sufficient. △ Less

Submitted 21 May, 2012; originally announced May 2012.

Comments: 12 pages, 2 figures

Report number: Roma01.Math.NA MSC Class: 65F15; 65F35

Journal ref: J. Comput. Appl. Math. 272 (2014), 304-312

Computing the structured pseudospectrum of a Toeplitz matrix and its extreme points

Authors: Paolo Buttà, Nicola Guglielmi, Silvia Noschese

Abstract: The computation of the structured pseudospectral abscissa and radius (with respect to the Frobenius norm) of a Toeplitz matrix is discussed and two algorithms based on a low rank property to construct extremal perturbations are presented. The algorithms are inspired by those considered in [SIAM J. Matrix Anal. Appl., 32 (2011), pp. 1166-1192] for the unstructured case, but their extension to struc… ▽ More The computation of the structured pseudospectral abscissa and radius (with respect to the Frobenius norm) of a Toeplitz matrix is discussed and two algorithms based on a low rank property to construct extremal perturbations are presented. The algorithms are inspired by those considered in [SIAM J. Matrix Anal. Appl., 32 (2011), pp. 1166-1192] for the unstructured case, but their extension to structured pseudospectra and analysis presents several difficulties. Natural generalizations of the algorithms, allowing to draw significant sections of the structured pseudospectra in proximity of extremal points are also discussed. Since no algorithms are available in the literature to draw such structured pseudospectra, the approach we present seems promising to extend existing software tools (Eigtool, Seigtool) to structured pseudospectra representation for Toeplitz matrices. We discuss local convergence properties of the algorithms and show some applications to a few illustrative examples. △ Less

Submitted 1 February, 2012; originally announced February 2012.

Comments: 21 pages, 11 figures

Report number: Roma01.Math.NA MSC Class: 65F15; 65L07

Journal ref: SIAM J. Matrix Anal. Appl. 33 (2012), 1300-1319

Boundary effects in the gradient theory of phase transitions

Authors: L. Bertini, P. Buttà, A. Garroni

Abstract: We consider the van der Waals' free energy functional, with a scaling small parameter epsilon, in the plane domain given by the first quadrant, and inhomogeneous Dirichlet boundary conditions. The boundary data are chosen in such a way that the interface between the pure phases tends to be horizontal and is pinned at some point on the y-axis which approaches zero as epsilon converges to zero. We s… ▽ More We consider the van der Waals' free energy functional, with a scaling small parameter epsilon, in the plane domain given by the first quadrant, and inhomogeneous Dirichlet boundary conditions. The boundary data are chosen in such a way that the interface between the pure phases tends to be horizontal and is pinned at some point on the y-axis which approaches zero as epsilon converges to zero. We show that there exists a critical scaling for the pinning point, such that, as the small parameter epsilon tends to zero, the competing effects of repulsion from the boundary and penalization of gradients play a role in determining the optimal shape of the (properly rescaled) interface. This result is achieved by means of an asymptotic development of the free energy functional. As a consequence, such analysis is not restricted to minimizers but also encodes the asymptotic probability of fluctuations. △ Less

Submitted 3 January, 2011; originally announced January 2011.

Comments: 20 pages, 1 figure

Report number: Roma01.Math.MP MSC Class: 82B26; 49J45 (Primary); 35J20; 82B24 (Secondary)

Journal ref: SIAM J. Math. Anal. 44 (2012), 926-945

Navier-Stokes equations on the flat cylinder with vorticity production on the boundary

Authors: Carlo Boldrighini, Paolo Buttà

Abstract: We study the two-dimensional Navier-Stokes system on a flat cylinder with the usual Dirichlet boundary conditions for the velocity field u. We formulate the problem as an infinite system of ODE's for the natural Fourier components of the vorticity, and the boundary conditions are taken into account by adding a vorticity production at the boundary. We prove equivalence to the original Navier-Stok… ▽ More We study the two-dimensional Navier-Stokes system on a flat cylinder with the usual Dirichlet boundary conditions for the velocity field u. We formulate the problem as an infinite system of ODE's for the natural Fourier components of the vorticity, and the boundary conditions are taken into account by adding a vorticity production at the boundary. We prove equivalence to the original Navier-Stokes system and show that the decay of the Fourier modes is exponential for any positive time in the periodic direction, but it is only power-like in the other direction. △ Less

Submitted 22 February, 2010; originally announced February 2010.

Comments: 25 pages

Report number: Roma01.Math.MP MSC Class: 35Q30; 76D05; 76N10

Journal ref: Nonlinearity 24 (2011), 2639-2662

Speedy motions of a body immersed in an infinitely extended medium

Authors: P. Buttà, G. Ferrari, C. Marchioro

Abstract: We study the motion of a classical point body of mass M, moving under the action of a constant force of intensity E and immersed in a Vlasov fluid of free particles, interacting with the body via a bounded short range potential Psi. We prove that if its initial velocity is large enough then the body escapes to infinity increasing its speed without any bound "runaway effect". Moreover, the body a… ▽ More We study the motion of a classical point body of mass M, moving under the action of a constant force of intensity E and immersed in a Vlasov fluid of free particles, interacting with the body via a bounded short range potential Psi. We prove that if its initial velocity is large enough then the body escapes to infinity increasing its speed without any bound "runaway effect". Moreover, the body asymptotically reaches a uniformly accelerated motion with acceleration E/M. We then discuss at a heuristic level the case in which Psi(r) diverges at short distances like g r^{-a}, g,a>0, by showing that the runaway effect still occurs if a<2. △ Less

Submitted 12 October, 2009; originally announced October 2009.

Comments: 15 pages

Report number: Roma01.Math.MP MSC Class: 70F40; 70B05; 70F45

Journal ref: J. Stat. Phys. 140 (2010), 1182-1194

Directed current in quasi-adiabatically ac-driven nonlinear systems

Authors: Paolo Buttà, Piero Negrini

Abstract: We rigorously prove the existence of directed transport for a certain class of ac-driven nonlinear one dimensional systems, namely the generation of transport with a preferred direction in the absence of a net driving force. We rigorously prove the existence of directed transport for a certain class of ac-driven nonlinear one dimensional systems, namely the generation of transport with a preferred direction in the absence of a net driving force. △ Less

Submitted 28 October, 2008; v1 submitted 24 October, 2008; originally announced October 2008.

Comments: 13 pages (one misprint has been corrected)

Report number: Roma01.Math.MP MSC Class: 82C70; 70K05

Journal ref: Math. Models Methods Appl. Sci. 20 (2010), 223-236

Resonances and O-curves in Hamiltonian systems

Authors: Paolo Buttà, Piero Negrini

Abstract: We investigate the problem of the existence of trajectories asymptotic to elliptic equilibria of Hamiltonian systems in the presence of resonances. We investigate the problem of the existence of trajectories asymptotic to elliptic equilibria of Hamiltonian systems in the presence of resonances. △ Less

Submitted 4 September, 2007; originally announced September 2007.

Comments: 12 pages

Report number: Roma01.Math.MP MSC Class: 70H14; 37J40

Journal ref: Regul. Chaotic Dyn. 12 (2007), 521-530

Molecular Dynamics Simulation of Vascular Network Formation

Authors: Paolo Buttà, Fiammetta Cerreti, Vito D. P. Servedio, Livio Triolo

Abstract: Endothelial cells are responsible for the formation of the capillary blood vessel network. We describe a system of endothelial cells by means of two-dimensional molecular dynamics simulations of point-like particles. Cells' motion is governed by the gradient of the concentration of a chemical substance that they produce (chemotaxis). The typical time of degradation of the chemical substance intr… ▽ More Endothelial cells are responsible for the formation of the capillary blood vessel network. We describe a system of endothelial cells by means of two-dimensional molecular dynamics simulations of point-like particles. Cells' motion is governed by the gradient of the concentration of a chemical substance that they produce (chemotaxis). The typical time of degradation of the chemical substance introduces a characteristic length in the system. We show that point-like model cells form network resembling structures tuned by this characteristic length, before collapsing altogether. Successively, we improve the non-realistic point-like model cells by introducing an isotropic strong repulsive force between them and a velocity dependent force mimicking the observed peculiarity of endothelial cells to preserve the direction of their motion (persistence). This more realistic model does not show a clear network formation. We ascribe this partial fault in reproducing the experiments to the static geometry of our model cells that, in reality, change their shapes by elongating toward neighboring cells. △ Less

Submitted 29 April, 2009; v1 submitted 7 July, 2007; originally announced July 2007.

Comments: 10 pages, 3 figures, 2 of which composite with 8 pictures each. Accepted on J.Stat.Mech. (2009). Appeared at the poster session of StatPhys23, Genoa, Italy, July 13 (2007)

Report number: Roma01.Math.MP

Journal ref: Journal of Statistical Mechanics: Theory and Experiments (2009) P05013, 10 pages

Dobrushin states in the φ^4_1 model

Authors: L. Bertini, S. Brassesco, P. Buttà

Abstract: We consider the van der Waals free energy functional in a bounded interval with inhomogeneous Dirichlet boundary conditions imposing the two stable phases at the endpoints. We compute the asymptotic free energy cost, as the length of the interval diverges, of shifting the interface from the midpoint. We then discuss the effect of thermal fluctuations by analyzing the φ^4_1-measure with Dobrushin… ▽ More We consider the van der Waals free energy functional in a bounded interval with inhomogeneous Dirichlet boundary conditions imposing the two stable phases at the endpoints. We compute the asymptotic free energy cost, as the length of the interval diverges, of shifting the interface from the midpoint. We then discuss the effect of thermal fluctuations by analyzing the φ^4_1-measure with Dobrushin boundary conditions. In particular, we obtain a nontrivial limit in a suitable scaling in which the length of the interval diverges and the temperature vanishes. The limiting state is not translation invariant and describes a localized interface. This result can be seen as the probabilistic counterpart of the variational convergence of the associated excess free energy. △ Less

Submitted 27 November, 2006; originally announced November 2006.

Comments: 34 pages

Report number: Roma01.Math.MP MSC Class: 82B24; 60K35; 81Q20

Journal ref: Arch. Ration. Mech. Anal. 190 (2008), 477-516

Soft and hard wall in a stochastic reaction diffusion equation

Authors: L. Bertini, S. Brassesco, P. Buttà

Abstract: We consider a stochastically perturbed reaction diffusion equation in a bounded interval, with boundary conditions imposing the two stable phases at the endpoints. We investigate the asymptotic behavior of the front separating the two stable phases, as the intensity of the noise vanishes and the size of the interval diverges. In particular, we prove that, in a suitable scaling limit, the front e… ▽ More We consider a stochastically perturbed reaction diffusion equation in a bounded interval, with boundary conditions imposing the two stable phases at the endpoints. We investigate the asymptotic behavior of the front separating the two stable phases, as the intensity of the noise vanishes and the size of the interval diverges. In particular, we prove that, in a suitable scaling limit, the front evolves according to a one-dimensional diffusion process with a non-linear drift accounting for a "soft" repulsion from the boundary. We finally show how a "hard" repulsion can be obtained by an extra diffusive scaling. △ Less

Submitted 27 November, 2006; originally announced November 2006.

Comments: 33 pages

Report number: Roma01.Math.MP MSC Class: 82C26; 60H15; 35K57

Journal ref: Arch. Ration. Mech. Anal. 190 (2008), 307-345

On the propagation of a perturbation in an anharmonic system

Authors: P. Buttà, E. Caglioti, S. Di Ruzza, C. Marchioro

Abstract: We give a not trivial upper bound on the velocity of disturbances in an infinitely extended anharmonic system at thermal equilibrium. The proof is achieved by combining a control on the non equilibrium dynamics with an explicit use of the state invariance with respect to the time evolution. We give a not trivial upper bound on the velocity of disturbances in an infinitely extended anharmonic system at thermal equilibrium. The proof is achieved by combining a control on the non equilibrium dynamics with an explicit use of the state invariance with respect to the time evolution. △ Less

Submitted 20 December, 2006; v1 submitted 12 July, 2006; originally announced July 2006.

Comments: 14 pages

Report number: Roma01.Math.MP MSC Class: 82C05; 82D20

Journal ref: J. Stat. Phys. 127 (2007), 313-325

Vertex Intrinsic Fitness: How to Produce Arbitrary Scale-Free Networks

Authors: Vito D. P. Servedio, Guido Caldarelli, Paolo Butta`

Abstract: We study a recent model of random networks based on the presence of an intrinsic character of the vertices called fitness. The vertices fitnesses are drawn from a given probability distribution density. The edges between pair of vertices are drawn according to a linking probability function depending on the fitnesses of the two vertices involved. We study here different choices for the probabili… ▽ More We study a recent model of random networks based on the presence of an intrinsic character of the vertices called fitness. The vertices fitnesses are drawn from a given probability distribution density. The edges between pair of vertices are drawn according to a linking probability function depending on the fitnesses of the two vertices involved. We study here different choices for the probability distribution densities and the linking functions. We find that, irrespective of the particular choices, the generation of scale-free networks is straightforward. We then derive the general conditions under which scale-free behavior appears. This model could then represent a possible explanation for the ubiquity and robustness of such structures. △ Less

Submitted 29 September, 2003; originally announced September 2003.

Comments: 4 pages, 3 figures, RevTex

Report number: Roma01.Math.MP

Journal ref: Physical Review E 70 (2004), 056126, 4 pages

Hydrodynamic Limit of Brownian Particles Interacting with Short and Long Range Forces

Authors: Paolo Butta`, Joel L. Lebowitz

Abstract: We investigate the time evolution of a model system of interacting particles, moving in a $d$-dimensional torus. The microscopic dynamics are first order in time with velocities set equal to the negative gradient of a potential energy term $Ψ$ plus independent Brownian motions: $Ψ$ is the sum of pair potentials, $V(r)+γ^d J(γr)$, the second term has the form of a Kac potential with inverse range… ▽ More We investigate the time evolution of a model system of interacting particles, moving in a $d$-dimensional torus. The microscopic dynamics are first order in time with velocities set equal to the negative gradient of a potential energy term $Ψ$ plus independent Brownian motions: $Ψ$ is the sum of pair potentials, $V(r)+γ^d J(γr)$, the second term has the form of a Kac potential with inverse range $γ$. Using diffusive hydrodynamical scaling (spatial scale $γ^{-1}$, temporal scale $γ^{-2}$) we obtain, in the limit $γ\downarrow 0$, a diffusive type integro-differential equation describing the time evolution of the macroscopic density profile. △ Less

Submitted 24 September, 1998; originally announced September 1998.

Comments: 37 pages, in TeX (compile twice), to appear on J. Stat. Phys., e-mail addresses: pbutta@math.rutgers.edu, lebowitz@math.rutgers.edu

Report number: Roma01.Math.MP

Journal ref: J. Statist. Phys. 94 (1999), pp. 653-694