Source author record

Andrea Carati

Andrea Carati 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
8topics
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)

preprint2020arXiv

Tsallis distributions, their relaxations and the relation $Δt \cdot ΔE \simeq h$, in the dynamical fluctuations of a classical model of a crystal

We report the results of a numerical investigation, performed in the frame of dynamical systems' theory, for a realistic model of a ionic crystal for which, due to the presence of long--range Coulomb interactions, the Gibbs distribution is not well defined. Taking initial data with a Maxwell-Boltzmann distribution for the mode-energies $E_k$, we study the dynamical fluctuations, computing the moduli of the the energy-changes $|E_k(t)-E_k(0)|$. The main result is that they follow Tsallis distributions, which relax to distributions close to Maxwell-Boltzmann ones; indications are also given that the system remains correlated. The relaxation time $τ$ depends on specific energy $\varepsilon$, and for the curve $τ$ vs, $\varepsilon$ one has two results. First, there exists an energy threshold $\varepsilon_0$, above which the curve has the form $$ τ\cdot \varepsilon \simeq h\ , $$ where, unexpectedly, Planck's constant $h$ shows up. In terms of the standard deviation $ΔE$ of a mode-energy (for which one has $ΔE=\varepsilon$), denoting by $Δt$ the relaxation time $τ$, the relation reads $Δt \cdot ΔE \simeq h$, which reminds of the Heisenberg uncertainty relation. Moreover, the threshold corresponds to zero-point energy. Indeed, the quantum value of the latter is $hν/2$ ( where $ν$ is the characterisic infrared frequency of the system), while we find $\varepsilon \simeq hν/4$, so that one only has a discrepancy of a factor 2. So it seems that lack of full chaoticity manifests itself, in Statistical Thermodynamics, through quantum-like phenomena.

preprint2019arXiv

Approach to equilibrium via Tsallis distributions in a realistic ionic--crystal model and in the FPU model

In Statistical Mechanics, Tsallis distributions were apparently conceived in connection with systems presenting long--range interactions. In fact, they were observed in numerical computations for models of such a type, as occurring in the approach to equilibrium, i.e., to a Maxwell--Boltzmann distribution. Here we exhibit two apparently new results. The first one is that Tsallis distributions occur also in an ionic--crystal model with long--range Coulomb forces, which is so realistic as to reproduce in an impressively good way the experimental infrared spectra. Thus such distributions may be expected to be actual physical features of crystals. The second result is that Tsallis distributions occur in the standard short--range FPU model too, so that the presence of long--range interactions is not a necessary condition for Tsallis distributions to occur. In fact, this is in agreement with a previous result of the first author in connection with the statistics of return times for the classical FPU model. We thus confirm the thesis advanced by Tsallis himself, that the relevant property for a dynamical system to present Tsallis distributions is that its dynamics should be not fully chaotic, a property which is known to actually pertain to long--range systems.

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.

preprint2016arXiv

Numerical study of the Transverse Diffusion coefficient for a one component model of a plasma

We report the results of MD numerical simulations for a one component model of a plasma in the weakly coupled regime, at different values of temperature $T$ and applied magnetic field $\vec B$, in which the diffusion coefficient $D_{\perp}$ transverse to the field is estimated. We find that there exists a threshold in temperature, at which an inversion occurs, namely, for $T$ above the threshold the diffusion coefficient $D_{\perp}$ starts decreasing as $T$ increases. This is at variance with the behavior predicted by the Bohm law $D_{\perp}\sim T/B$, which actually holds below the threshold. In addition we find that, for temperatures above such a threshold, another transition occurs, now with respect to the values of the magnetic field: for weak magnetic fields the diffusion coefficients scales as $1/B^2$, in agreement with the predictions of the standard kinetics theory, while it apparently saturates when the field strength is sufficiently increased.

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 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.

preprint2014arXiv

Classical microscopic theory of dispersion, emission and absorption of light in dielectrics

This paper is a continuation of a recent one in which, apparently for the first time, the existence of polaritons in ionic crystals was proven in a microscopic electrodynamic theory. This was obtained through an explicit computation of the dispersion curves. Here the main further contribution consists in studying electric susceptibility, from which the spectrum can be inferred. We show how susceptibility is obtained by the Green--Kubo methods of Hamiltonian statistical mechanics, and give for it a concrete expression in terms of time--correlation functions. As in the previous paper, here too we work in a completely classical framework, in which the electrodynamic forces acting on the charges are all taken into account, both the retarded forces and the radiation reaction ones. So, in order to apply the methods of statistical mechanics, the system has to be previously reduced to a Hamiltonian one. This is made possible in virtue of two global properties of classical electrodynamics, namely, the Wheeler--Feynman identity and the Ewald resummation properties, the proofs of which were already given for ordered system. The second contribution consists in formulating the theory in a completely general way, so that in principle it applies also to disordered systems such as glasses, or liquids or gases, provided the two general properties mentioned above continue to hold. A first step in this direction is made here by providing a completely general proof of the Wheeler--Feynman identity, which is shown to be the counterpart of a general causality property of classical electrodynamics. Finally it is shown how a line spectrum can appear at all in classical systems, as a counterpart of suitable stability properties of the motions, with a broadening due to a coexistence of chaoticity.

preprint2014arXiv

Theoretical thermodynamic analysis of a closed-cycle process for the conversion of heat into electrical energy by means of a distiller and an electrochemical cell

We analyse a device aimed at the conversion of heat into electrical energy, based on a closed cycle in which a distiller generates two solutions at different concentrations, and an electrochemical cell consumes the concentration difference, converting it into electrical current. We first study an ideal model of such a process. We show that, if the device works at a single fixed pressure (i.e. with a ``single effect''), then the efficiency of the conversion of heat into electrical power can approach the efficiency of a reversible Carnot engine operating between the boiling temperature of the concentrated solution and that of the pure solvent. When two heat reservoirs with a higher temperature difference are available, the overall efficiency can be incremented by employing an arrangement of multiple cells working at different pressures (``multiple effects''). We find that a given efficiency can be achieved with a reduced number of effects by using solutions with a high boiling point elevation.

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$.

preprint2012arXiv

A series expansion for the time autocorrelation of dynamical variables

We present here a general iterative formula which gives a (formal) series expansion for the time autocorrelation of smooth dynamical variables, for all Hamiltonian systems endowed with an invariant measure. We add some criteria, theoretical in nature, which enable one to decide whether the decay of the correlations is exponentially fast or not. One of these criteria is implemented numerically for the case of the Fermi-Pasta-Ulam system, and we find indications which might suggest a sub-exponential decay of the time autocorrelation of a relevant dynamical variable.

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.

preprint2011arXiv

Exponentially long stability times for a nonlinear lattice in the thermodynamic limit

In this paper, we construct an adiabatic invariant for a large 1--$d$ lattice of particles, which is the so called Klein Gordon lattice. The time evolution of such a quantity is bounded by a stretched exponential as the perturbation parameters tend to zero. At variance with the results available in the literature, our result holds uniformly in the thermodynamic limit. The proof consists of two steps: first, one uses techniques of Hamiltonian perturbation theory to construct a formal adiabatic invariant; second, one uses probabilistic methods to show that, with large probability, the adiabatic invariant is approximately constant. As a corollary, we can give a bound from below to the relaxation time for the considered system, through estimates on the autocorrelation of the adiabatic invariant.