Source author record

Alberto Maiocchi

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

13 published item(s)

preprint2018arXiv

A large probability averaging Theorem for the defocousing NLS

We consider the nonlinear Schroedinger equation on the one dimensional torus, with a defocousing polynomial nonlinearity and study the dynamics corresponding to initial data in a set of large measure with respect to the Gibbs measure. We prove that along the corresponding solutions the modulus of the Fourier coefficients is approximately constant for times of order $β^{2+ς}$, $β$ being the inverse of the temperature and $ς$ a positive number (we prove $ς= 1/10$). The proof is obtained by adapting to the context of Gibbs measure for PDEs some tools of Hamiltonian perturbation theory.

preprint2016arXiv

Infrared optical properties of $α$ quartz by molecular dynamics simulations

This paper is concerned with theoretical estimates of the refractive--index curves for quartz, obtained by the Kubo formulæ in the classical approximation, through MD simulations for the motions of the ions. Two objectives are considered. The first one is to understand the role of nonlinearities in situations where they are very large, as at the $α$--$β$ structural phase transition. We show that on the one hand they don't play an essential role in connection with the form of the spectra in the infrared. On the other hand they play an essential role in introducing a chaoticity which involves a definite normal mode. This might explain why that mode is Raman active in the $α$ phase, but not in the $β$ phase. The second objective concerns whether it is possible in a microscopic model to obtain normal mode frequencies, or peak frequencies in the optical spectra, that are in good agreement with the experimental data for quartz. Notwithstanding a lot of effort, we were unable to find results agreeing better than about 6%, as apparently also occurs in the whole available literature. We interpret this fact as indicating that some essential qualitative feature is lacking in all models which consider, as the present one, only short--range repulsive potentials and unretarded long--range electric forces.

preprint2015arXiv

A replacement of the Lorentz law for the shape of the spectral lines in the infrared region

We propose a new phenomenological law for the shape of the spectral lines in the infrared, which accounts for the exponential decay of the extinction coefficient in the high frequency region, observed in many spectra. We apply this law to the measured infrared spectra of LiF, NaCl and MgF$_2$, finding a good agreement, over a wide range of frequencies.

preprint2015arXiv

Agreement of classical Kubo theory with the infrared dispersion curves $n(ω)$ of ionic crystals

The theoretical dispersion curves $n(ω)$ (refractive index versus frequency) of ionic crystals in the infrared domain are expressed, within the Green--Kubo theory, in terms of a time correlation function involving the motion of the ions only. The aim of this paper is to investigate how well the experimental data are reproduced by a classical approximation of the theory, in which the time correlation functions are expressed in terms of the ions orbits. We report the results of molecular dynamics (MD) simulations for the ions motions of a LiF lattice of 4096 ions at room temperature. The theoretical curves thus obtained are in surprisingly good agreement with the experimental data, essentially over the whole infrared domain. This shows that at room temperature the motion of the ions develops essentially in a classical regime.

preprint2015arXiv

The effective equation method

In this chapter we present a general method of constructing the effective equation which describes the behaviour of small-amplitude solutions for a nonlinear PDE in finite volume, provided that the linear part of the equation is a hamiltonian system with a pure imaginary discrete spectrum. The effective equation is obtained by retaining only the resonant terms of the nonlinearity (which may be hamiltonian, or may be not); the assertion that it describes the limiting behaviour of small-amplitude solutions is a rigorous mathematical theorem. In particular, the method applies to the three-- and four--wave systems. We demonstrate that different possible types of energy transport are covered by this method, depending on whether the set of resonances splits into finite clusters (this happens, e.g. in case of the Charney-Hasegawa-Mima equation), or is connected (this happens, e.g. in the case of the NLS equation if the space-dimension is at least two). For equations of the first type the energy transition to high frequencies does not hold, while for equations of the second type it may take place. In the case of the NLS equation we use next some heuristic approximation from the arsenal of wave turbulence to show that under the iterated limit "the volume goes to infinity", taken after the limit "the amplitude of oscillations goes to zero", the energy spectrum of solutions for the effective equation is described by a Zakharov-type kinetic equation. Evoking the Zakharov ansatz we show that stationary in time and homogeneous in space solutions for the latter equation have a power law form. Our method applies to various weakly nonlinear wave systems, appearing in plasma, meteorology and oceanology.

preprint2015arXiv

The Fermi-Pasta-Ulam system as a model for glasses

We show that the standard Fermi--Pasta--Ulam system, with a suitable choice for the interparticle potential, constitutes a model for glasses, and indeed an extremely simple and manageable one. Indeed, it allows one to describe the landscape of the minima of the potential energy and to deal concretely with any one of them, determining the spectrum of frequencies and the normal modes. A relevant role is played by the harmonic energy $\mathcal E$ relative to a given minimum, i.e., the expansion of the Hamiltonian about the minimum up to second order. Indeed we find that there exists an energy threshold in $\mathcal E$ such that below it the harmonic energy $\mathcal E$ appears to be an approximate integral of motion for the whole observation time. Consequently, the system remains trapped near the minimum, in what may be called a vitreous or glassy state. Instead, for larger values of $\mathcal E$ the system rather quickly relaxes to a final equilibrium state. Moreover we find that the vitreous states present peculiar statistical behaviors, still involving the harmonic energy $\mathcal E$. Indeed, the vitreous states are described by a Gibbs distribution with an effective Hamiltonian close to $\mathcal E$ and with a suitable effective inverse temperature. The final equilibrium state presents instead statistical properties which are in very good agreement with the Gibbs distribution relative to the full Hamiltonian of the system.

preprint2015arXiv

Time-averaging for weakly nonlinear CGL equations with arbitrary potentials

Consider weakly nonlinear complex Ginzburg--Landau (CGL) equation of the form: $$ u_t+i(-Δu+V(x)u)=εμΔu+ε\mathcal{P}( u),\quad x\in {R^d}\,, \quad(*) $$ under the periodic boundary conditions, where $μ\geqslant0$ and $\mathcal{P}$ is a smooth function. Let $\{ζ_1(x),ζ_2(x),\dots\}$ be the $L_2$-basis formed by eigenfunctions of the operator $-Δ+V(x)$. For a complex function $u(x)$, write it as $u(x)=\sum_{k\geqslant1}v_kζ_k(x)$ and set $I_k(u)=\frac{1}{2}|v_k|^2$. Then for any solution $u(t,x)$ of the linear equation $(*)_{ε=0}$ we have $I(u(t,\cdot))=const$. In this work it is proved that if equation $(*)$ with a sufficiently smooth real potential $V(x)$ is well posed on time-intervals $t\lesssim ε^{-1}$, then for any its solution $u^ε(t,x)$, the limiting behavior of the curve $I(u^ε(t,\cdot))$ on time intervals of order $ε^{-1}$, as $ε\to0$, can be uniquely characterized by a solution of a certain well-posed effective equation: $$ u_t=εμ\triangle u+εF(u), $$ where $F(u)$ is a resonant averaging of the nonlinearity $\mathcal{P}(u)$. We also prove a similar results for the stochastically perturbed equation, when a white in time and smooth in $x$ random force of order $\sqrtε$ is added to the right-hand side of the equation. The approach of this work is rather general. In particular, it applies to equations in bounded domains in $R^d$ under Dirichlet boundary conditions.

preprint2014arXiv

Resonant averaging for small solutions of stochastic NLS equations

We consider the free linear Schrödinger equation on a torus $\mathbb T^d$, perturbed by a hamiltonian nonlinearity, driven by a random force and damped by a linear damping: $$ u_t -iΔu +iνρ|u|^{2q_*}u = - νf(-Δ) u + \sqrtν\,\frac{d}{d t}\sum_{k\in \mathbb Z^d} b_lβ^k(t)e^{ik\cdot x} \ . $$ Here $u=u(t,x),\ x\in\mathbb T^d$, $0<ν\ll 1$, $q_*\in\mathbb N$, $f$ is a positive continuous function, $ρ$ is a positive parameter and $β^k(t)$ are standard independent complex Wiener processes. We are interested in limiting, as $ν\to0$, behaviour of distributions of solutions for this equation and of its stationary measure. Writing the equation in the slow time $τ=νt$, we prove that the limiting behaviour of the both is described by the effective equation $$ u_τ+ f(-Δ) u = -iF(u)+\frac{d}{dτ}\sum b_kβ^k(τ)e^{ik\cdot x} \, $$ where the nonlinearity $F(u)$ is made out of the resonant terms of the monomial $ |u|^{2q_*}u$. We explain the relevance of this result for the problem of weak turbulence

preprint2014arXiv

The limit of small Rossby numbers for randomly forced quasi-geostrophic equation on $β$-plane

We consider the 2d quasigeostrophic equation on the $β$-plane for the stream function $ψ$, with dissipation and a random force: $$ (*)\qquad (-Δ+K)ψ_t - ρJ(ψ, Δψ) -βψ_x= \langle \text{random force}\rangle -κΔ^2ψ+Δψ, $$ where $ψ=ψ(t,x,y), \ x\in\mathbb{R}/2πL\mathbb{Z}, \ y\in \mathbb{R}/2π\mathbb{Z}$. For typical values of the horizontal period $L$ we prove that the law of the action-vector of a solution for $(*)$ (formed by the halves of the squared norms of its complex Fourier coefficients) converges, as $β\to\infty$, to the law of an action-vector for solution of an auxiliary effective equation, and the stationary distribution of the action-vector for solutions of $(*)$ converges to that of the effective equation. Moreover, this convergence is uniform in $κ\in(0,1]$. The effective equation is an infinite system of stochastic equations which splits into invariant subsystems of complex dimension $\le3$; each of these subsystems is an integrable hamiltonian system, coupled with a Langevin thermostat. Under the iterated limits $\lim_{L=ρ\to\infty} \lim_{β\to\infty}$ and $\lim_{κ\to 0} \lim_{β\to\infty}$ we get similar systems. In particular, none of the three limiting systems exhibits the energy cascade to high frequencies.

preprint2013arXiv

An averaging theorem for FPU in the thermodynamic limit

Consider an FPU chain composed of $N\gg 1$ particles, and endow the phase space with the Gibbs measure corresponding to a small temperature $β^{-1}$. Given a fixed $K<N$, we construct $K$ packets of normal modes whose energies are adiabatic invariants (i.e., are approximately constant for times of order $β^{1-a}$, $a>0$) for initial data in a set of large measure. Furthermore, the time autocorrelation function of the energy of each packet does not decay significantly for times of order $β$. The restrictions on the shape of the packets are very mild. All estimates are uniform in the number $N$ of particles and thus hold in the thermodynamic limit $N\to\infty$, $β>0$.

preprint2013arXiv

Derivation of the Kolmogorov-Zakharov equation from the resonant-averaged stochastic NLS equation

We suggest a new derivation of a kinetic equation of Kolmogorov-Zakharov (KZ) type for the spectrum of the weakly nonlinear Schrödinger equation with stochastic forcing. The kynetic equation is obtained as a result of a double limiting procedure. Firstly, we consider the equation on a finite box with periodic boundary conditions and send the size of the nonlinearity and of the forcing to zero, while the time is correspondingly rescaled; then, the size of the box is sent to infinity (with a suitable rescaling of the solution). We report here the results of the first limiting procedure, analyzed with full rigour in arXiv:1311.6793, and show how the second limit leads to a kinetic equation for the spectrum, if some further hypotheses (commonly employed in the weak turbulence theory) are accepted. Finally we show how to derive from these equations the KZ spectra.

preprint2013arXiv

Resonant averaging for weakly nonlinear stochastic Schrödinger equations

We consider the free linear Schroedinger equation on a torus $\mathbb T^d$, perturbed by a Hamiltonian nonlinearity, driven by a random force and damped by a linear damping: $$u_t -iΔu +iνρ|u|^{2q_*}u = - νf(-Δ) u + \sqrtν\,\frac{d}{d t}\sum_{k\in \mathbb Z^d} b_kβ^k(t)e^{ik\cdot x} \ . $$ Here $u=u(t,x),\ x\in\mathbb T^d$, $0<ν\ll1$, $q_*\in\mathbb N\cup\{0\}$, $f$ is a positive continuous function, $ρ$ is a positive parameter and $β^k(t)$ are standard independent complex Wiener processes. We are interested in limiting, as $ν\to0$, behaviour of solutions for this equation and of its stationary measure. Writing the equation in the slow time $τ=νt$, we prove that the limiting behaviour of the both is described by the effective equation $$ u_τ+ f(-Δ) u = -iF(u)+\frac{d}{dτ}\sum b_kβ^k(τ)e^{ik\cdot x} \, $$ where the nonlinearity $F(u)$ is made out of the resonant terms of the monomial $ |u|^{2q_*}u$. We explain the relevance of this result for the problem of weak turbulence.

preprint2012arXiv

Transition to chaos in magnetized, weakly coupled plasmas

We report the results of numerical simulations for a model of a one component plasma (a system of N point electrons with mutual Coulomb interactions) in a uniform stationary magnetic field. We take N up to 512, with periodic boundary conditions, and macroscopic parameters corresponding to the weak coupling regime, with a coupling parameter Γ=1/64. We find that a transition from order to chaos takes place when the density is increased or the field decreased so that the ratio ω_p/ω_c between plasma and cyclotron frequencies becomes of order 1 (or equivalently the ratio r_L/λ_D between Larmor radius and Debye length becomes of order 1). The result is in agreement with the theoretical prediction obtained in [1], on the basis of an old estimate of Iglesias, Lebowitz and MacGowan [2] for the intensity of the electric field acting on one electron and due to all the other ones. A comparison can be made with the threshold obtained from kinetic theory arguments, which corresponds to the condition ν_{ee}/ω_c=1, where ν_{ee} is the electron collision frequency. The latter threshold has a completely different dependence on the physical parameters and, for Γ=1/64, gives a critical value of ω_p about 80 times larger.