LE JOURNAL DE PHYSIQUE
|
|
- Stella Lang
- 6 years ago
- Views:
Transcription
1 Nous The Tome 43 No 5 MAI 1982 LE JOURNAL DE PHYSIQUE J. Physique 43 (1982) MAi 1982, Classification Physics Abstracts Approach to equilibrium in a chain of nonlinear oscillators F. Fucito (*) (+), F. Marchesoni (**), E. Marinari (*) (+), G. Parisi (***), L. Peliti (*) (+ +), S. Ruffo (**) and A. Vulpiani (*) (*) Istituto di Fisica «G. Marconi», Università di Roma, Italy. (**) Istituto di Fisica, Università di Pisa, Italy, and INFN, Pisa. (***) Istituto di Fisica, Facoltà di Ingegneria, Università di Roma, Italy and INFN, Frascati. (+) INFN, Roma. (++) GNSMCNR, Unità di Roma. (Recu le 6 aout 1981, révisé le 30 décembre, accepté le 18 janvier 1982) 2014 Résumé. considérons l approche à l équilibre thermodynamique d une chaîne d oscillateurs non linéaires couplés où l énergie est initialement fournie aux modes de plus grande longueur d onde. La considération des effets des singularités dans le plan complexe de la solution des équations du mouvement permet de prédire qu à des temps relativement courts le spectre d énergie a une dépendance exponentielle en fonction du nombre d onde k lorsque celuici est grand. Elle permet également de distinguer entre deux régions temporelles : une à temps courts où la pente de l exponentielle dépend linéairement de ln t, et une à des temps intermédiaires où elle est proportionnelle à (In t)1/2. Ces prédictions sont en accord avec les simulations numériques que nous avons faites sur ce système. Abstract We consider the approach to thermodynamical equilibrium in a chain of coupled nonlinear oscillators when energy is initially fed only to the longest wavelength modes. Consideration of the effects of singularities in the complex plane of the solutions of the equations of motion allows us to predict that at not too long times the energy spectrum has an exponential tail in the wavenumber k and to distinguish between two time regions : a short time region where the slope of the exponential shows a linear dependence on In t, and an intermediate region where it is proportional to (In t)1/2. These predictions are successfully compared with the numerical simulations we have performed on the system. 1. Introduction. approach to equilibrium of a classical hamiltonian system is still an open and interesting problem. The main rigorous results are the KAM [1] and Sinai [2] theorems. They describe the behaviour of two opposite extreme cases of a nonintegrable hamiltonian system. The Sinai theorem shows the validity of the ergodic hypothesis for a gas of hard spheres. The main consequence of the KAM theorem is the permanence of invariant tori in presence of a small perturbation of an integrable system. In this case the long time behaviour of the system will not be described by a Boltzmann distribution. The range of validity of the KAM theorem for a system with a large number of degrees of freedom is not yet understood. Numerical simulations have been extensively examined for a linear chain of coupled nonlinear oscillators, starting from the famous work of Fermi, Pasta and Ulam [3]. More recent works [46] have shown that the regions of phase space where ordered motion takes place are much more extended than on the basis of the KAM theorem. Some expected authors [6] stress the existence of a threshold in energy which separates ordered from disordered behaviour. All these works study long times in order to distinguish between ordered and stochastic asymptotic behaviour. The present work analyses short times in order to understand the way chaotic behaviour sets in. Our system is a discrete version of a classical onedimensional field obeying a nonlinear Klein Article published online by EDP Sciences and available at
2 We 708 Gordon equation. This system has already been considered in reference [6]. We consider initial conditions such that energy is concentrated in the longest wavelength modes and we study how it propagates to shorter wavelengths. In the times we consider the bulk of the energy is still contained in the initially excited modes, and only a small quantity is fed to the shorter wavelengths. In order to allow for an approximate analytic treatment we focus our attention on the process by which this small amount is distributed. We are then able to compare our predictions with the results of the numerical simulations. as follows : the model is The paper is organized introduced in section 2, where the analytic treatment at short and intermediate times is given. The results of the numerical simulations as well as their comparison with our theoretical predictions are given in section 3. In section 4 we draw some conclusions and point out some possible directions for further work. 2. Analytical treatment. study the onedimensional nonlinear KleinGordon equation : where the real, onecomponent field p(x, t) is defined on the interval L/2 x L/2 with periodic boundary conditions. The corresponding hamiltonian density is given by : where the field 7r is canonically conjugate to po The system may be either considered as a onedimensional elastic string in an external anharmonic potential, or as a onedimensional classical field theory. We study the time evolution of the field ~ as a function of the initial condition p(x, 0), It is known that if H is finite at t 0, the solution of equation (2.1) exists and is unique [7]. It is useful to define the space Fourier transform of the field p as follows : We expect the field p to reach asymptotically a thermal equilibrium distribution, given by a Boltzmann factor exp( f3h) for some value of the inverse temperature f3 determined by the initial conditions. In this case, at values of the wavenumber k so large that the mass and nonlinear terms of H are negligible, one would have : This behaviour of W corresponds to functions p(x, t) which are not differentiable with respect to x. Now it is known that since T(x, 0) is analytical as a function of x, the solution T(x, t) will remain analytical at any finite time t [7]. Equation (2.6) can only be valid for infinite time. This means that, as time goes on, singularities of p(x, t) appear in the complex x plane which creep toward the real axis and accumulate onto it at infinite times. We show below that these singularities are simple poles. If we assume this to be the case, then we may relate these singularities to the large k behaviour of W by means of the theorem of residues. Let us consider the integral defining 0 (eq. (2.3)) in the complex x plane. At positive (negative) k we close the contour in the lower (upper) complex half plane. We then obtain the following expression for 0(k, t) : where the sum runs over all poles located in the relevant half plane, Rj being their residue and xj + iy; their location. We therefore obtain the following asymptotic behaviour of W(k, t ) at large k : where ys(t) is the imaginary part of the location of the pole which lies nearest to the real axis. Our strategy is then to evaluate the most likely value of ys(t) by extending the approach of Frisch and Morf [8] to a deterministic partial differential equation. We should in principle evaluate the analytic continuation to complex x of the solution p(x, t) of equation (2.1). We prefer instead to consider the solution 9(x + iy, t ) of the complex partial differential equation : We shall consider as a typical initial condition : We are interested in the behaviour, in the region I k I > ki, of the spectrum W(k, t), defined by : where ~ is an analytical function of z x + iy. Such an equation is the complexification of equation (2.1) with the hypothesis that T(x, t ) can be analytically continued. In other words we follow the commutative diagram given below :
3 Force 709 We first consider the case of an ordinary differential equation obtained from equation (2.9) by neglecting the laplacian >> term a2paz2. Let us consider the solution 4>«fJo, t) of this equation at real times t, as a function of its (complex) initial condition (po : If we write 4> 4>R + i4>i we have the following differential equations for 4>R 4>1 : imaginary axisx we may introduce a potential according to : The time needed to reach infinity starting from the point (po i(p" 0 may be then easily computed from the conservation of «energy» in the onedimensional motion : The force field associated to this system is drawn in figure 1 (only the first quadrant is shown). Although the origin is a stable point, we find a pair of unstable points on the imaginary axis, from which a trajectory escapes to infinity. The point at infinity appears as a saddle point, and may be reached in a finite times from points on the imaginary axis. This is the origin of the singularities of 4>(cpo, t). It is then clear that these singularities are simple poles, located at that value of To from which infinity may be reached in time t. We may estimate the location of these poles by computing the time needed to reach infinity starting from a given point on the imaginary 45 axis. Since the motion is in this case onedimensional (all along the At large ~o the time T behaves as : On the contrary, no singularity will appear at any time if for the time T diverges (although slowly) as Im (po I approaches this value from above. From equation (2.14) we therefore obtain the location of the pole of 4>( CPo, t) at large CPo s, therefore at short t s : If the effect of the laplacian term in equation (2.9) may be neglected at short enough times, the above analysis may be also applied to the onedimensional chain as follows : the solution T(z, t) of equation (2. 9) will be given by 4>(cpo(z), t), where CPo(z) is the analytical continuation of the initial condition T(x, 0) to the complex x plane. In particular, a typical initial condition yields in the complex x plane We have therefore at large y I : z Fig. 1. field of equation (2.11). Only the upper right quadrant is shown. The position of the instable point on the imaginary axis (1m cp m/gl/2) is marked. According to our hypothesis, the singularity of CPo(z, t) will be located where Im CPo, as estimated
4 A We 710 from equation (2. 19), will be of the order of that given by equation (2.16). Therefore We now consider the effect of the previously neglected laplacian term. We claim that the analysis we just made retains its validity at very short times. The laplacian term will not in fact be able to prevent a 9(z, t) from escaping to infinity if its initial imaginary part is large enough. Localized o bursts» will therefore appear in (p(x, t), whose spatial coherence is only due to correlations present in the initial conditions. The laplacian term will on the contrary be essential in introducing singularities in the part of the complex x plane which corresponds to small initial values of Im cpo If the term were not present, Im T would remain bounded at all times. If however! Im cp [ were to become larger than m/gi/2, the mechanism we just described would drive it to infinity in a very short time, and a singularity will be produced. We must now estimate the time needed for Im T to become larger than m/g 1/2. The analysis becomes simpler if we consider a harmonic chain ( g 0) in the limit of infinite length (L oo). We expect this analysis to be also valid for our case if g is not too large. It is known that the onepoint probability distribution function of a classical harmonic field in onedimensional is gaussian [9]. The variance Q2 of the gaussian will be then given at small y s by : cannot therefore draw conclusions about the behaviour at very long times before thermal equilibrium is reached. We summarize the results for W(k, t) at short and intermediate times : 3. Numerical results. present in this section the results of numerical simulations of equation (2.1 ) for different initial conditions, in order to verify the validity of the analytical results we obtained in the previous section. We consider a chain of N points with periodic boundary conditions which represents a discretization of equation (2.1 ). We set L (the length of the chain) equal to 1 ; and we consider the field cp(x) to be represented by gj cpu Ax 1/2), where Ax 1/N. The periodic boundary conditions then imply 9jlN 9j The Fourier transform of the discrete field Pj is defined as follows : The probability p, for I 1m cp to be larger than m/g 1/2 will be therefore given by : The time needed for this to happen will therefore be of order 1/pc. This yields the following estimate. of YS : The main effect of the nonlinear terms in this regime will be to change the value of (12. If a2 were time independent, equation (2.21) will be essentially correct. Let us distinguish between the role of short and long wavelength modes. At the times we are interested in, most of the energy is contained in the long wavelength modes, which may be assumed to be in a kind of thermal equilibrium among themselves. Their contribution to (J2 may be then considered as essential modes will ly constant in time. The short wavelength however also contribute to (J2. As long as W(k, t) is small in the large k region, their contribution is negligible. As time goes on, however, W(k, t ) will start increasing, what will increase the value of Q2 and fasten therefore the transfer of energy to short wavelength modes. This triggers a catastrophic process which our analytical tools are unable to handle. We Fig. 2. typical spectrum, In Wn(t) vs. n for A 5, g 5, N64, t2, T 1.
5 Slope 711 The quantity corresponding to the W(k, t ) of the previous section will be : We have performed simulations for N 64 and N 128. The initial conditions were chosen to be : significant differences have been seen in runs with different At. The results we report have been obtained via computations with 15 or 16 significant digits. Since part of the arguments in the previous section are probabilistic, it is convenient to introduce the average of Wn(t) centred around the time t : all other Fourier coefficient having been set to zero. The integration of equation (2.1 ) was performed by means of the central difference algorithm [10] : where the discretized force Fj({ T }) is given by : In the limit At 0 equation (3. 3) reduces to In order to eliminate oscillations, T must be chosen to be of the order of magnitude of the period of oscillation of a mode of wavenumber kn 2 nn/ L in a harmonic chain. Since m is quite small in our case, this yields comparatively large values of T, of order 0.1 : 0.2 in the short time region and 2 : 10 in the intermediate time region. The deviations of W,,(t) from its average may be understood on the basis of the statistical properties of the field (p(x, t) at large k, which is analysed in detail in a forthcoming paper [11]. We now present the results for short times. We have verified that the spectrum Wn(t ) can be very well described by an exponential in the region of not too small n. A typical example is shown in figure 2. It is easy to obtain the value of the slope S (t ) of this exponential from the data : The initial condition cpj(o) 0 corresponds to (pj( et) w;(0). We have done simulations keeping m and varying A and g in the range 0.5 : 10. At has been taken in the range x We have checked that the energy is conserved within 0.1 %; no This slope is related to ys(t) of the previous section by : We show in figure 3 the dependence of S (t ) on In (tagi/2) with different values of t, g and A. We see Fig. 3. S (arbitrary units) vs. In (Ag1/2 t), N 64, T 0.1 : 0.2.
6 Slope Slope 712 Fig. 4. S (arbitrary units) vs. (A2 g In t) 1/2, N 64, T 5. that the data gather on an universal curve, which is essentially a straight line, giving therefore a confirmation to the analysis of short times performed in section 2. This plot may be therefore considered as a check of equation (2.20). All error bars represent essentially uncertainties in the estimation of S(t). We have verified that the value of Wn(t ) coincides with its time average at short times. This is due to the dominance of a single complex singularity. As the time increases and we enter the intermediate time region, we expect that more than one complex singularity at roughly the same distance from the real x axis will Fig. 5. z S (arbitrary units) vs. (In t) 112, N 128, A 3, g 5, T 10 except for t 8 where T 2t
7 We We 713 appear. Oscillations will therefore be present in Wn(t) as a consequence of the interference among such singularities. In the intermediate time region we expect that S (t) oc (A2 g In t) 112. We have plotted S vs. this variable in figure 4, varying either t at fixed A and g, or g at fixed t and A. We see again that the data gather on a single straight line. Therefore the functional dependence of S (t ) is as expected on the basis of equation (2.23) for not too long times. The behaviour at longer times of a solution with a given value of A and g is shown in figure 5. One sees both the (In t) li2 behaviour and a break at longer times. This is connected with the observation done at the end of the previous section. 4. Conclusions. have seen that it is possible to give definite predictions about the behaviour of a nonintegrable hamiltonian system with an infinite number of degrees of freedom in the nonasymptotic time region. This has been achieved by the use of the analytical continuation of the equations of motion in the complex domain, a method introduced in a similar context by Frisch and Morf [8]. The role of the infinite volume limit has also been clarified. The free evolution of the system is indeed such that 9 may only become arbitrarily large in the infinite volume limit. In this case no linear perturbation may be considered small if one follows the system at long enough times. It is possible that intermittance concepts are relevant in this context [8, 11]. Further work aims at clarifying this connection [11]. One of our main results is that the system approaches equilibrium with a logarithmic dependence on t, so that the nonequilibrium spectrum may persist for extremely long times, and may be mistaken for a stationary state if the observation time is not sufficiently long. It is amusing to remark that this quasiequilibrium distribution is similar to Wien s law for black body radiation, with a slowly varying «Planck s constant». Acknowledgments. thank U. Frisch for having made available to us the paper in ref. [8] before publication. We also thank L. Galgani, A. Giorgilli and M. Vitaletti for helpful conversations and constructive criticism. References [1] KOLMOGOROV, A. N., Dokl. Akad. Nauk SSSR 98 (1954) 527. ARNOL D, V. I., Russ. Math. Surv. 18 (1963) 9. MOSER, J., Nachr. Akad. Wiss. Göttingen, Math. Phys. K1., 2, 1 (1962) 15. [2] SINAI, Y. G., Russ. Math. Surv. 25 (1970) 137. [3] FERMI, E., PASTA, J. and ULAM, S., Los Alamos Scientific Laboratory Report No. La1940 (1955). [4] HÉNON, H. and HEILES, C., Astron. J. 69 (1964) 73. [5] FORD, J., Adv. Chem. Phys. 24 (1973) 155. GUSTAVSON, F. G., Astron. J. 71 (1966) 670. BENETTIN, G., GALGANI, L. and STRELCYN, J. M., Phys. Rev. A14 (1976) BENETTIN, G., LO VECCHIO, G. and TENENBAUM, A., Phys. Rev. A22 (1980) [6] BUTERA, P., GALGANI, L., GIORGILLI, A., TAGLIANI, A. and SABATA, H., Nuovo Cimento 59B (1980) 81. [7] BREZIS, H., in : M. Atteia, D. Bancel and P. Gumovski (Eds.) : Nonlinear Problems of Analysis in geometry and Mechanics (Boston : Potman Adv. Publ. Program) 1981, p. 3. [8] FRISCH, U. and MORF, R., Phys. Rev. A 23 (1981) [9] VAN HEMMEN, J. L., Phys. Reports 65 (1980) 43. [10] See for example BENETTIN, G. et al. (1980) in ref. [5]. [11] FUCITO, F., MARCHESONI, F., SPERPEGLIONE, M. and VULPIANI, A., Pise preprint (1982).
LE JOURNAL DE PHYSIQUE - LETTRES
Les The Tome 46 No 6 15 MARS 1985 LE JOURNAL DE PHYSIQUE - LETTRES J. Physique Lett. 46 (1985) L-217 - L-222 15 MARS 1985, : L-217 Classification Physics Abstracts 75.50K Random free energies in spin glasses
More informationLE JOURNAL DE PHYSIQUE - LETTRES
Nous We Tome 46? 21 1 er NOVEMBRE 1985 LE JOURNAL DE PHYSIQUE - LETTRES J. Physique Lett. 46 (1985) L-985 - L-989 1 er NOVEMBRE 1985, L-985 Classification Physics Abstracts 64.70P - 05.40-75.40 Diffusion
More informationPasta-Ulam problem. The birth of nonlinear physics and dynamical chaos. Paolo Valentini
The Fermi-Pasta Pasta-Ulam problem The birth of nonlinear physics and dynamical chaos Paolo Valentini SUMMARY The Fermi-Pasta Pasta-Ulam experiment on MANIAC at LANL in 1955 The ergodic hypothesis The
More informationGlobal theory of one-frequency Schrödinger operators
of one-frequency Schrödinger operators CNRS and IMPA August 21, 2012 Regularity and chaos In the study of classical dynamical systems, the main goal is the understanding of the long time behavior of observable
More informationReflection of light by a random layered system
Nous We J. Physique 47 (1986) 1467-1475 SEPTEMBRE 1986, 1467 Classification Physics Abstracts 41. 10H - 71. 55J Reflection of light by a random layered system E. Bouchaud and M. Daoud Laboratoire Léon
More informationQuantum Quenches in Extended Systems
Quantum Quenches in Extended Systems Spyros Sotiriadis 1 Pasquale Calabrese 2 John Cardy 1,3 1 Oxford University, Rudolf Peierls Centre for Theoretical Physics, Oxford, UK 2 Dipartimento di Fisica Enrico
More informationStrong Chaos without Buttery Eect. in Dynamical Systems with Feedback. Via P.Giuria 1 I Torino, Italy
Strong Chaos without Buttery Eect in Dynamical Systems with Feedback Guido Boetta 1, Giovanni Paladin 2, and Angelo Vulpiani 3 1 Istituto di Fisica Generale, Universita di Torino Via P.Giuria 1 I-10125
More information1 Equal-time and Time-ordered Green Functions
1 Equal-time and Time-ordered Green Functions Predictions for observables in quantum field theories are made by computing expectation values of products of field operators, which are called Green functions
More informationKhinchin s approach to statistical mechanics
Chapter 7 Khinchin s approach to statistical mechanics 7.1 Introduction In his Mathematical Foundations of Statistical Mechanics Khinchin presents an ergodic theorem which is valid also for systems that
More informationCapillary rise between closely spaced plates : effect of Van der Waals forces
Capillary rise between closely spaced plates : effect of Van der Waals forces B. Legait, P.G. De Gennes To cite this version: B. Legait, P.G. De Gennes. Capillary rise between closely spaced plates : effect
More informationOn the smoothness of the conjugacy between circle maps with a break
On the smoothness of the conjugacy between circle maps with a break Konstantin Khanin and Saša Kocić 2 Department of Mathematics, University of Toronto, Toronto, ON, Canada M5S 2E4 2 Department of Mathematics,
More informationLee Yang zeros and the Ising model on the Sierpinski gasket
J. Phys. A: Math. Gen. 32 (999) 57 527. Printed in the UK PII: S35-447(99)2539- Lee Yang zeros and the Ising model on the Sierpinski gasket Raffaella Burioni, Davide Cassi and Luca Donetti Istituto Nazionale
More informationIrregular wavy flow due to viscous stratification
Irregular wavy flow due to viscous stratification T. Shlang, G.I. Sivashinsky, A.J. Babchin, A.L. Frenkel To cite this version: T. Shlang, G.I. Sivashinsky, A.J. Babchin, A.L. Frenkel. Irregular wavy flow
More informationLE JOURNAL DE PHYSIQUE - LETTRES
05.50 On We Tome 45 No 14 15 JUILLET 1984 LE JOURNAL DE PHYSIQUE - LETTRES J. Physique Lett. 45 (1984) L-701 - L-706 15 JUILLET 1984, L-701 Classification Physics Abstracts - - 05.20 64.60 On a model of
More information2 A. Morbidelli and A. Giorgilli relevant consequence is that the speed of the Arnold diusion (if any) turns out to be of superexponentially small siz
SUPEREXPONENTIAL STABILITY OF KAM TORI ALESSANDRO MORBIDELLI CNRS, Observatoire de la C^ote d'azur BP 229, 06304 { NICE Cedex 4, France. ANTONIO GIORGILLI Dipartimento di Matematica dell' Universita di
More informationAccurate critical exponents from the ϵ-expansion
Accurate critical exponents from the ϵ-expansion J.C. Le Guillou, J. Zinn-Justin To cite this version: J.C. Le Guillou, J. Zinn-Justin. Accurate critical exponents from the ϵ-expansion. Journal de Physique
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 3, 3 March 2006 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationOn the nonrelativistic binding energy for positive ions
On the nonrelativistic binding energy for positive ions G.I. Plindov, I.K. Dmitrieva To cite this version: G.I. Plindov, I.K. Dmitrieva. On the nonrelativistic binding energy for positive ions. Journal
More informationAgeostrophic instabilities of a front in a stratified rotating fluid
8 ème Congrès Français de Mécanique Grenoble, 27-3 août 27 Ageostrophic instabilities of a front in a stratified rotating fluid J. Gula, R. Plougonven & V. Zeitlin Laboratoire de Météorologie Dynamique
More informationMath 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
More informationExercise sheet n Compute the eigenvalues and the eigenvectors of the following matrices. C =
L2 - UE MAT334 Exercise sheet n 7 Eigenvalues and eigenvectors 1. Compute the eigenvalues and the eigenvectors of the following matrices. 1 1 1 2 3 4 4 1 4 B = 1 1 1 1 1 1 1 1 1 C = Which of the previous
More informationANNALES DE L I. H. P., SECTION A
ANNALES DE L I. H. P., SECTION A G. VELO An existence theorem for a massive spin one particle in an external tensor field Annales de l I. H. P., section A, tome 22, n o 3 (1975), p. 249-255
More informationNonintegrability and the Fourier heat conduction law
Nonintegrability and the Fourier heat conduction law Giuliano Benenti Center for Nonlinear and Complex Systems, Univ. Insubria, Como, Italy INFN, Milano, Italy In collaboration with: Shunda Chen, Giulio
More informationLE JOURNAL DE PHYSIQUE
On A Tome 46 N 8 AOÛT 1985 LE JOURNAL DE PHYSIQUE J. Physique 46 (1985) 1269-1275 Ao8T 1985, 1269 Classification Physics Abstracts 05.50-75.50K Separable and non-separable spin glass models F. Benamira
More informationAdsorption of chain molecules with a polar head a scaling description
Adsorption of chain molecules with a polar head a scaling description S. Alexander To cite this version: S. Alexander. Adsorption of chain molecules with a polar head a scaling description. Journal de
More informationMicroscopic Hamiltonian dynamics perturbed by a conservative noise
Microscopic Hamiltonian dynamics perturbed by a conservative noise CNRS, Ens Lyon April 2008 Introduction Introduction Fourier s law : Consider a macroscopic system in contact with two heat baths with
More informationStochastic equations for thermodynamics
J. Chem. Soc., Faraday Trans. 93 (1997) 1751-1753 [arxiv 1503.09171] Stochastic equations for thermodynamics Roumen Tsekov Department of Physical Chemistry, University of Sofia, 1164 Sofia, ulgaria The
More informationPoisson s ratio effect of slope stability calculations
Poisson s ratio effect of slope stability calculations Murray Fredlund, & Robert Thode SoilVision Systems Ltd., Saskatoon, SK, Canada ABSTRACT This paper presents the results of a study on the effect of
More informationBasics of Statistical Mechanics
Basics of Statistical Mechanics Review of ensembles Microcanonical, canonical, Maxwell-Boltzmann Constant pressure, temperature, volume, Thermodynamic limit Ergodicity (see online notes also) Reading assignment:
More informationSome Collision solutions of the rectilinear periodically forced Kepler problem
Advanced Nonlinear Studies 1 (2001), xxx xxx Some Collision solutions of the rectilinear periodically forced Kepler problem Lei Zhao Johann Bernoulli Institute for Mathematics and Computer Science University
More informationTemperature-concentration diagram of polymer solutions
Temperatureconcentration diagram of polymer solutions M. Daoud, G. Jannink To cite this version: M. Daoud, G. Jannink. Temperatureconcentration diagram of polymer solutions. Journal de Physique, 1976,
More informationEffect of thermal noise and initial conditions in the dynamics of a mean-field ferromagnet
Effect of thermal noise and initial conditions in the dynamics of a mean-field ferromagnet O. Golinelli, B. Derrida To cite this version: O. Golinelli, B. Derrida. Effect of thermal noise and initial conditions
More informationUNDERSTANDING BOLTZMANN S ANALYSIS VIA. Contents SOLVABLE MODELS
UNDERSTANDING BOLTZMANN S ANALYSIS VIA Contents SOLVABLE MODELS 1 Kac ring model 2 1.1 Microstates............................ 3 1.2 Macrostates............................ 6 1.3 Boltzmann s entropy.......................
More informationThe fractal nature of a diffusion front and the relation to percolation
The fractal nature of a diffusion front and the relation to percolation Bernard Sapoval, Michel Rosso, JeanFrançois Gouyet To cite this version: Bernard Sapoval, Michel Rosso, JeanFrançois Gouyet. The
More informationA Reflexive toy-model for financial market
A Reflexive toy-model for financial market An alternative route towards intermittency. Luigi Palatella Dipartimento di Fisica, Università del Salento via Arnesano, I-73100 Lecce, Italy SM&FT 2011 - The
More informationGENERALIZED OPTICAL BISTABILITY AND CHAOS IN A LASER WITH A SATURABLE ABSORBER
GENERALIZED OPTICAL BISTABILITY AND CHAOS IN A LASER WITH A SATURABLE ABSORBER E. Arimondo, F. De Tomasi, B. Zambon, F. Papoff, D. Hennequin To cite this version: E. Arimondo, F. De Tomasi, B. Zambon,
More informationEffect of the different parts of the interaction potential on the line profile in neutral atom-atom collisions
Effect of the different parts of the interaction potential on the line profile in neutral atomatom collisions N.F. Allard, Y.G. Biraud To cite this version: N.F. Allard, Y.G. Biraud. Effect of the different
More informationarxiv: v1 [cond-mat.stat-mech] 7 Nov 2017
About thermometers and temperature arxiv:1711.02349v1 [cond-mat.stat-mech] 7 Nov 2017 M. Baldovin Dipartimento di Fisica, Sapienza Università di Roma, p.le A. Moro 2, 00185 Roma, Italy E-mail: marco.baldovin@roma1.infn.it
More informationTHE MANDELSTAM REPRESENTATION IN PERTURBATION THEORY
THE MANDELSTAM REPRESENTATION IN PERTURBATION THEORY P. V. Landshoff, J. C. Polkinghorne, and J. C. Taylor University of Cambridge, Cambridge, England (presented by J. C. Polkinghorne) 1. METHODS The aim
More informationSPLITTING OF SEPARATRICES FOR (FAST) QUASIPERIODIC FORCING. splitting, which now seems to be the main cause of the stochastic behavior in
SPLITTING OF SEPARATRICES FOR (FAST) QUASIPERIODIC FORCING A. DELSHAMS, V. GELFREICH, A. JORBA AND T.M. SEARA At the end of the last century, H. Poincare [7] discovered the phenomenon of separatrices splitting,
More informationClassical dynamics of a non-integrable Hamiltonian near coupling-induced resonance islands
J. Phys. A: Math. Gen. 29 (1996) 5963 5977. Printed in the UK Classical dynamics of a non-integrable Hamiltonian near coupling-induced resonance islands Marc Joyeux Laboratoire de Spectrométrie Physique
More information2 Formulation. = arg = 2 (1)
Acoustic di raction by an impedance wedge Aladin H. Kamel (alaahassan.kamel@yahoo.com) PO Box 433 Heliopolis Center 11757, Cairo, Egypt Abstract. We consider the boundary-value problem for the Helmholtz
More informationOn the large time behavior of solutions of fourth order parabolic equations and ε-entropy of their attractors
On the large time ehavior of solutions of fourth order paraolic equations and ε-entropy of their attractors M.A. Efendiev & L.A. Peletier Astract We study the large time ehavior of solutions of a class
More informationEnergy Level Sets for the Morse Potential
Energy Level Sets for the Morse Potential Fariel Shafee Department of Physics Princeton University Princeton, NJ 08540 Abstract: In continuation of our previous work investigating the possibility of the
More informationCopyright (c) 2006 Warren Weckesser
2.2. PLANAR LINEAR SYSTEMS 3 2.2. Planar Linear Systems We consider the linear system of two first order differential equations or equivalently, = ax + by (2.7) dy = cx + dy [ d x x = A x, where x =, and
More informationLandau Bogolubov Energy Spectrum of Superconductors
Landau Bogolubov Energy Spectrum of Superconductors L.N. Tsintsadze 1 and N.L. Tsintsadze 1,2 1. Department of Plasma Physics, E. Andronikashvili Institute of Physics, Tbilisi 0128, Georgia 2. Faculty
More informationCascades and statistical equilibrium in shell models of turbulence
PHYSICAL REVIEW E VOLUME 53, UMBER 5 MAY 1996 Cascades and statistical equilibrium in shell models of turbulence P. D. Ditlevsen and I. A. Mogensen The iels Bohr Institute, Department for Geophysics, University
More informationarxiv:cond-mat/ v1 29 Dec 1996
Chaotic enhancement of hydrogen atoms excitation in magnetic and microwave fields Giuliano Benenti, Giulio Casati Università di Milano, sede di Como, Via Lucini 3, 22100 Como, Italy arxiv:cond-mat/9612238v1
More informationChopping time of the FPU α-model
Chopping time of the FPU α-model A. Carati A. Ponno November 24, 2017 Abstract We study, both numerically and analytically, the time needed to observe the breaking of an FPU α-chain in two or more pieces,
More informationA Study of the Thermal Properties of a One. Dimensional Lennard-Jones System
A Study of the Thermal Properties of a One Dimensional Lennard-Jones System Abstract In this study, the behavior of a one dimensional (1D) Lennard-Jones (LJ) system is simulated. As part of this research,
More information~ij is a function of the positions of the points
64.40C Nous We J. Phys. France 49 (1988) 2019-2025 DÉCEMBRE 1988, 2019 Classification Physics Abstracts - - 02.10 75.10NR The Euclidean matching problem Marc Mézard (1) and Giorgio Parisi (2) (1) Laboratoire
More informationChaos, Quantum Mechanics and Number Theory
Chaos, Quantum Mechanics and Number Theory Peter Sarnak Mahler Lectures 2011 Hamiltonian Mechanics (x, ξ) generalized coordinates: x space coordinate, ξ phase coordinate. H(x, ξ), Hamiltonian Hamilton
More informationHamiltonian Dynamics
Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;
More informationAbstracts. Furstenberg The Dynamics of Some Arithmetically Generated Sequences
CHAOS AND DISORDER IN MATHEMATICS AND PHYSICS Monday 10:00-11:00 Okounkov Algebraic geometry of random surfaces 11:30-12:30 Furstenberg Dynamics of Arithmetically Generated Sequences 12:30-14:30 lunch
More informationarxiv: v1 [math.ap] 20 Nov 2007
Long range scattering for the Maxwell-Schrödinger system with arbitrarily large asymptotic data arxiv:0711.3100v1 [math.ap] 20 Nov 2007 J. Ginibre Laboratoire de Physique Théorique Université de Paris
More informationNumerical Analysis of 2-D Ising Model. Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011
Numerical Analysis of 2-D Ising Model By Ishita Agarwal Masters in Physics (University of Bonn) 17 th March 2011 Contents Abstract Acknowledgment Introduction Computational techniques Numerical Analysis
More informationA Two-dimensional Mapping with a Strange Attractor
Commun. math. Phys. 50, 69 77 (1976) Communications in Mathematical Physics by Springer-Verlag 1976 A Two-dimensional Mapping with a Strange Attractor M. Henon Observatoire de Nice, F-06300 Nice, France
More informationBasics of Statistical Mechanics
Basics of Statistical Mechanics Review of ensembles Microcanonical, canonical, Maxwell-Boltzmann Constant pressure, temperature, volume, Thermodynamic limit Ergodicity (see online notes also) Reading assignment:
More informationInequivalent Representations of a q-oscillator Algebra in a Quantum q-gas
CBPF-NF-028/95 Inequivalent Representations of a q-oscillator Algebra in a Quantum q-gas M.R-Monteiro a and L.M.C.S. Rodrigues b Centro Brasileiro de Pesquisas Físicas - CBPF Rua Dr. Xavier Sigaud, 50
More informationA Search for the Simplest Chaotic Partial Differential Equation
A Search for the Simplest Chaotic Partial Differential Equation C. Brummitt University of Wisconsin-Madison, Department of Physics cbrummitt@wisc.edu J. C. Sprott University of Wisconsin-Madison, Department
More informationMean field theory for Heisenberg spin glasses
Mean field theory for Heisenberg spin glasses G. Toulouse, M. Gabay To cite this version: G. Toulouse, M. Gabay. Mean field theory for Heisenberg spin glasses. Journal de Physique Lettres, 1981, 42 (5),
More informationSummary of free theory: one particle state: vacuum state is annihilated by all a s: then, one particle state has normalization:
The LSZ reduction formula based on S-5 In order to describe scattering experiments we need to construct appropriate initial and final states and calculate scattering amplitude. Summary of free theory:
More informationarxiv:hep-th/ v3 16 May 1996
BNL-63106 An Exact Solution for Quantum Tunneling in a Dissipative System arxiv:hep-th/9605081v3 16 May 1996 Li Hua Yu National Synchrotron Light Source, Brookhaven National Laboratory, N.Y.11973 Abstract
More informationIn Arnold s Mathematical Methods of Classical Mechanics (1), it
Near strongly resonant periodic orbits in a Hamiltonian system Vassili Gelfreich* Mathematics Institute, University of Warwick, Coventry CV4 7AL, United Kingdom Communicated by John N. Mather, Princeton
More informationLow-temperature state of Ising spin glass
Low-temperature state of Ising spin glass L.B. Ioffe, M.V. Feigel Man To cite this version: L.B. Ioffe, M.V. Feigel Man. Low-temperature state of Ising spin glass. Journal de Physique Lettres, 1983, 44
More informationOn the Onsager-Yang-Value of the Spontaneous Magnetization
On the Onsager-Yang-Value of the Spontaneous Magnetization G. Benettin, G. Gallavotti*, G. Jona-Lasinio, A. L. Stella Istituto di Fisica del Universita, Padova, Italy Received November 10, 1972 Abstract:
More informationBLOW-UP OF COMPLEX SOLUTIONS OF THE 3-d NAVIER-STOKES EQUATIONS AND BEHAVIOR OF RELATED REAL SOLUTIONS
BLOW-UP OF COMPLEX SOLUTIONS OF THE 3-d NAVIER-STOKES EQUATIONS AND BEHAVIOR OF RELATED REAL SOLUTIONS BLOW-UP OF COMPLEX SOLUTIONS OF THE 3-d NAVIER-STOKES EQUATIONS AND BEHAVIOR OF RELATED REAL SOLUTIONS
More informationA Cantor set of tori with monodromy near a focus focus singularity
INSTITUTE OF PHYSICS PUBLISHING Nonlinearity 17 (2004) 1 10 NONLINEARITY PII: S0951-7715(04)65776-8 A Cantor set of tori with monodromy near a focus focus singularity Bob Rink Mathematics Institute, Utrecht
More informationTime fractional Schrödinger equation
Time fractional Schrödinger equation Mark Naber a) Department of Mathematics Monroe County Community College Monroe, Michigan, 48161-9746 The Schrödinger equation is considered with the first order time
More informationgoverning the dynamics of wrinkled flame fronts
47.20 J. Physique 46 (1985) 14851494 SEPTEMBRE 1985, 1485 Classification Physics Abstracts 05.45 47.70F Application of pole decomposition to an equation governing the dynamics of wrinkled flame fronts
More informationRelaxation times in an anharmonic crystal with diluted. impurities. Lapo Casetti. Scuola Normale Superiore, Piazza dei Cavalieri 7, Pisa, Italy
Relaxation times in an anharmonic crystal with diluted impurities Lapo Casetti Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy Roberto Livi y Dipartimento di Fisica, Universita di Bologna,
More informationWhat really matters in Hilbert-space stochastic processes
What really matters in Hilbert-space stochastic processes February 6, 207 arxiv:702.04524v [quant-ph] 5 Feb 207 Giancarlo Ghirardi a, Oreste Nicrosini b and Alberto Rimini c a Abdus Salam ICTP, Strada
More informationOPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS
Stochastics and Dynamics c World Scientific Publishing Company OPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS BRIAN F. FARRELL Division of Engineering and Applied Sciences, Harvard University Pierce Hall, 29
More informationObservations on the ponderomotive force
Observations on the ponderomotive force D.A. Burton a, R.A. Cairns b, B. Ersfeld c, A. Noble c, S. Yoffe c, and D.A. Jaroszynski c a University of Lancaster, Physics Department, Lancaster LA1 4YB, UK b
More informationExact Functional Renormalization Group for Wetting Transitions in Dimensions.
EUROPHYSICS LETTERS Europhys. Lett., 11 (7), pp. 657-662 (1990) 1 April 1990 Exact Functional Renormalization Group for Wetting Transitions in 1 + 1 Dimensions. F. JULICHER, R. LIPOWSKY (*) and H. MULLER-KRUMBHAAR
More informationAPPLICATION OF A STOCHASTIC REPRESENTATION IN NUMERICAL STUDIES OF THE RELAXATION FROM A METASTABLE STATE
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 7 (1), 83-90 (2001) APPLICATION OF A STOCHASTIC REPRESENTATION IN NUMERICAL STUDIES OF THE RELAXATION FROM A METASTABLE STATE FERDINANDO DE PASQUALE 1, ANTONIO
More informationEquilibrium correlations and heat conduction in the Fermi-Pasta-Ulam chain.
Equilibrium correlations and heat conduction in the Fermi-Pasta-Ulam chain. Abhishek Dhar International centre for theoretical sciences TIFR, Bangalore www.icts.res.in Suman G. Das (Raman Research Institute,
More informationarxiv: v1 [nlin.ps] 9 May 2015
Scaling properties of generalized two-dimensional Kuramoto-Sivashinsky equations V. Juknevičius Institute of Theoretical Physics and Astronomy, Vilnius University, A. Goštauto 2, LT-008 Vilnius, Lithuania
More informationThe Fermi Pasta Ulam problem as a challenge for the foundations of physics
The Fermi Pasta Ulam problem as a challenge for the foundations of physics A. Carati L. Galgani Università di Milano, Dipartimento di Matematica, Via Saldini 50, 20133 Milano, Italy A. Giorgilli Università
More informationLecture 11 : Overview
Lecture 11 : Overview Error in Assignment 3 : In Eq. 1, Hamiltonian should be H = p2 r 2m + p2 ϕ 2mr + (p z ea z ) 2 2 2m + eφ (1) Error in lecture 10, slide 7, Eq. (21). Should be S(q, α, t) m Q = β =
More informationOPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS
Stochastics and Dynamics, Vol. 2, No. 3 (22 395 42 c World Scientific Publishing Company OPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS Stoch. Dyn. 22.2:395-42. Downloaded from www.worldscientific.com by HARVARD
More informationAnalysis of Hermite s equation governing the motion of damped pendulum with small displacement
Vol. 10(12), pp. 364-370, 30 June, 2015 DOI: 10.5897/IJPS2015.4364 Article Number: 9A709F454056 ISSN 1992-1950 Copyright 2015 Author(s) retain the copyright of this article http://www.academicjournals.org/ijps
More informationPHYSICAL REVIEW LETTERS
PHYSICAL REVIEW LETTERS VOLUME 80 1 JUNE 1998 NUMBER 22 Field-Induced Stabilization of Activation Processes N. G. Stocks* and R. Mannella Dipartimento di Fisica, Università di Pisa, and Istituto Nazionale
More informationInvaded cluster dynamics for frustrated models
PHYSICAL REVIEW E VOLUME 57, NUMBER 1 JANUARY 1998 Invaded cluster dynamics for frustrated models Giancarlo Franzese, 1, * Vittorio Cataudella, 1, * and Antonio Coniglio 1,2, * 1 INFM, Unità di Napoli,
More informationNON LINEAR PULSE EVOLUTION IN SEEDED AND CASCADED FELS
NON LINEAR PULSE EVOLUTION IN SEEDED AND CASCADED FELS L. Giannessi, S. Spampinati, ENEA C.R., Frascati, Italy P. Musumeci, INFN & Dipartimento di Fisica, Università di Roma La Sapienza, Roma, Italy Abstract
More informationarxiv:cond-mat/ v1 [cond-mat.dis-nn] 29 Mar 2006
Conductance fluctuations in the localized regime A. M. Somoza, J. Prior and M. Ortuño arxiv:cond-mat/0603796v1 [cond-mat.dis-nn] 29 Mar 2006 Departamento de Física, Universidad de Murcia, Murcia 30.071,
More informationON THE ASYMPTOTICS OF GREEN S FUNCTIONS OF ELLIPTIC OPERATORS WITH CONSTANT COEFFICIENTS. Shmuel Agmon
Séminaires & Congrès 9, 2004, p. 13 23 ON THE ASYMPTOTICS OF GREEN S FUNCTIONS OF ELLIPTIC OPERATORS WITH CONSTANT COEFFICIENTS by Shmuel Agmon Abstract. In this paper we discuss the following problem.
More informationNonlinear Single-Particle Dynamics in High Energy Accelerators
Nonlinear Single-Particle Dynamics in High Energy Accelerators Part 4: Canonical Perturbation Theory Nonlinear Single-Particle Dynamics in High Energy Accelerators There are six lectures in this course
More informationGlobal well-posedness for semi-linear Wave and Schrödinger equations. Slim Ibrahim
Global well-posedness for semi-linear Wave and Schrödinger equations Slim Ibrahim McMaster University, Hamilton ON University of Calgary, April 27th, 2006 1 1 Introduction Nonlinear Wave equation: ( 2
More informationInelastic scattering of 30 Mev polarized protons from 112Cd
Inelastic scattering of 30 Mev polarized protons from 112Cd R. De Swiniarski, G. Bagieu, DinhLien Pham, M. Massaad, J.Y. Grossiord, A. Guichard To cite this version: R. De Swiniarski, G. Bagieu, DinhLien
More informationStatistical switching kinetics in ferroelectrics
Statistical switching kinetics in ferroelectrics X.J. Lou Department of Materials Science and Engineering, National University of Singapore, 117574, Singapore Abstract: By assuming a more realistic nucleation
More informationThe Smaller (SALI) and the Generalized (GALI) Alignment Index Methods of Chaos Detection: Theory and Applications. Haris Skokos
The Smaller (SALI) and the Generalized (GALI) Alignment Index Methods of Chaos Detection: Theory and Applications Haris Skokos Max Planck Institute for the Physics of Complex Systems Dresden, Germany E-mail:
More informationStability of Feedback Control Systems: Absolute and Relative
Stability of Feedback Control Systems: Absolute and Relative Dr. Kevin Craig Greenheck Chair in Engineering Design & Professor of Mechanical Engineering Marquette University Stability: Absolute and Relative
More informationarxiv: v1 [math.ds] 19 Dec 2012
arxiv:1212.4559v1 [math.ds] 19 Dec 2012 KAM theorems and open problems for infinite dimensional Hamiltonian with short range Xiaoping YUAN December 20, 2012 Abstract. Introduce several KAM theorems for
More informationMathematical Structures of Statistical Mechanics: from equilibrium to nonequilibrium and beyond Hao Ge
Mathematical Structures of Statistical Mechanics: from equilibrium to nonequilibrium and beyond Hao Ge Beijing International Center for Mathematical Research and Biodynamic Optical Imaging Center Peking
More informationarxiv:cond-mat/ v1 4 Aug 2003
Conductivity of thermally fluctuating superconductors in two dimensions Subir Sachdev arxiv:cond-mat/0308063 v1 4 Aug 2003 Abstract Department of Physics, Yale University, P.O. Box 208120, New Haven CT
More informationarxiv: v1 [math.ds] 16 Jan 2018
arxiv:1801.05199v1 [math.ds] 16 Jan 2018 The Fermi Pasta Ulam problem and its underlying integrable dynamics: an approach through Lyapunov Exponents G. Benettin, S. Pasquali, A. Ponno Abstract FPU models,
More informationRandom walks on the Sierpinski Gasket
Random walks on the Sierpinski Gasket R. Friedberg, O. Martin To cite this version: R. Friedberg, O. Martin. Random walks on the Sierpinski Gasket. Journal de Physique, 1986, 47 (10), pp.16631669. .
More informationDETERMINATION OF MODEL VALID PREDICTION PERIOD USING THE BACKWARD FOKKER-PLANCK EQUATION
.4 DETERMINATION OF MODEL VALID PREDICTION PERIOD USING THE BACKWARD FOKKER-PLANCK EQUATION Peter C. Chu, Leonid M. Ivanov, and C.W. Fan Department of Oceanography Naval Postgraduate School Monterey, California.
More informationarxiv:cond-mat/ v1 [cond-mat.soft] 19 May 2004
arxiv:cond-mat/0405440v1 [cond-mat.soft] 19 May 2004 Condensate localization in a quasi-periodic structure Y. Eksioglu, P. Vignolo and M.P. Tosi NEST-INFM and Scuola Normale Superiore, Piazza dei Cavalieri
More information