A Toy Model. Viscosity
|
|
- Barbra Atkins
- 6 years ago
- Views:
Transcription
1 A Toy Model for Sponsoring Viscosity Jean BELLISSARD Westfälische Wilhelms-Universität, Münster Department of Mathematics SFB 878, Münster, Germany Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics
2 Main References C. A. Angell Formation of Glasses from Liquids and Biopolymers, Science, 267, No (Mar. 31, 1995), pp T. Egami, Elementary Excitation and Energy Landscape in Simple Liquids, Mod. Phys. Lett. B, 28, (2014), :1-19. J. Bellissard Delone Sets and Material Science: a Program, in Mathematics of Aperiodic Order, J. Kellendonk, D. Lenz, J. Savinien, Eds., Progress in Mathematics series, 309, Birkhäuser, pp , (2015).
3 Content Warning This talk is reporting on a work in progress. 1. Motivation 2. Anankeons and Phonons 3. Constructing the Model 4. Computing Maxwell s Time 5. Conclusion
4 I - Motivation
5 Motivation C. A. Angell, Formation of Glasses from Liquids and Biopolymers, Science, 267, No (Mar. 31, 1995), pp
6 Motivation J. H. Magill, D. J. Plazek, Physical Properties of Aromatic Hydrocarbons. II. Solidification Behavior of 1,3,5-Tri-α-Naphthylbenzene, J. Chem. Phys., 45, (1967),
7 Motivation Change of slope in Arrhenius behavior
8 Motivation Comparison between the Maxwell relaxation time τ M and the relaxation time τ LC for local configurations. T. Iwashita, D. M. Nicholson, T. Egami, Phys. Rev. Lett., 110, :5, (2013).
9 Motivation Specific Heat for Fragile Glasses M. H. Cohen, D. Turnbull, J. Chem. Phys., 31, 1164 (1959).
10 Motivation Specific Heat for Fragile Glasses: contributions of phonons & anankeons
11 Motivation At temperature T > T co larger than the crossover temperature, only the anankeons contribute to the specific heat. If T g < T < T co phonons and anankeon interact. Question: How?
12 II - Anankeons and Phonons
13 Atomic Configurations The set L of position of atomic nuclei is a Delone set, namely The minimum distance between atoms is 2r 0 > 0. The maximum diameter of a hole without atoms is 2r 1 <.
14 Atomic Configurations
15 Atomic Configurations
16 Atomic Configurations
17 Atomic Configurations
18 Atomic Configurations
19 Atomic Configurations
20 Atomic Configurations
21 Voronoi Cells Let L be Delone. If x L its Voronoi cell is defined by V(x) = {y R d ; y x < y x x L, x x} V(x) is open. Its closure T(x) = V(x) is called the Voronoi tile of x
22 Voronoi Cells Let L be Delone. If x L its Voronoi cell is defined by V(x) = {y R d ; y x < y x x L, x x} V(x) is open. Its closure T(x) = V(x) is called the Voronoi tile of x
23 Voronoi Cells Let L be Delone. If x L its Voronoi cell is defined by V(x) = {y R d ; y x < y x x L, x x} V(x) is open. Its closure T(x) = V(x) is called the Voronoi tile of x Proposition: If L is Delone, the Voronoi tile of any x L is a convex polytope
24 Voronoi Cells Let L be Delone. If x L its Voronoi cell is defined by V(x) = {y R d ; y x < y x x L, x x} V(x) is open. Its closure T(x) = V(x) is called the Voronoi tile of x Proposition: If L is Delone, the Voronoi tile of any x L is a convex polytope containing the closed ball B(x; r 0 ) and contained in the ball B(x; r 1 )
25 The Delone Graph Proposition: the Voronoi tiles of a Delone set touch face-to-face
26 The Delone Graph Proposition: the Voronoi tiles of a Delone set touch face-to-face Two atoms are nearest neighbors if their Voronoi tiles touch along a face of maximal dimension.
27 The Delone Graph Proposition: the Voronoi tiles of a Delone set touch face-to-face Two atoms are nearest neighbors if their Voronoi tiles touch along a face of maximal dimension. An edge is a pair of nearest neighbors. E denotes the set of edges.
28 The Delone Graph Proposition: the Voronoi tiles of a Delone set touch face-to-face Two atoms are nearest neighbors if their Voronoi tiles touch along a face of maximal dimension. An edge is a pair of nearest neighbors. E denotes the set of edges. The family G = (L, E) is the Delone graph.
29 The Delone Graph taken from J. D. Bernal, Nature, 183, , (1959)
30 The Delone Graph Modulo graph isomorphism, the Delone graph encodes the local topology
31 Atomic Movements small, deterministic, conserves local topology large, unpredictable, quick jump, change the local topology
32 III - Constructing the Model
33 Harmonic Motion The vibration of an atom around its equilibrium position is assumed to be harmonic. To simplify further, the oscillator will be supposed to be onedimensional. The frequency ω of the harmonic oscillator is defined by the curvature of the potential energy near the equilibrium position. Let q denotes the position of the harmonic particle relative to the equilibrium position. Then X will denote the phase space vector X = ωq q = [ u v ]
34 Harmonic Motion Equation of motion m d2 q dt 2 + kq = 0, ω = k m. Equivalently dx dt = ωjx, J = [ ] Or else X(t) = e ωtj X(0) = [ cos ωt sin ωt sin ωt cos ωt ] [ u(0) v(0) ]
35 Anankeon Interaction At random times < τ n 1 < τ n < τ n+1 <, the atom jumps quickly from one potential well to another one. Each jump will be considered as instantaneous. In the new potential well, the curvature is different, hence, after the time τ n, the new frequency is ω n. After the jump, the new relative phase space position is changed by a vector ξ n, namely X(τ n + 0) = X(τ n 0) + ξ n
36 Randomness: Assumptions The random times are Poissonian, namely the variables τ n+1 τ n are i.i.d, with exponential distribution and average E{τ n+1 τ n } = τ n+1 τ n = τ LC The frequencies ω n are also random and i.i.d, such that E(ω n ) = ω, E{(ω n ω) 2 } = σ 2 The phase-space initial positions ξ n are also random and i.i.d, with Gibbs distribution Prob{ξ n Λ} = e βm ξ 2 /2 d 2 ξ 2πk B T/m, β = 1 k B T Λ
37 Correlation Function The goal is to compute the stress-stress correlation function. It is sufficient to compute C f (t) = E{ f (X(t)) f (X(0))} for any complex valued function f defined on the phase space and vanishing at infinity. The viscosity is given by the Green-Kubo formula for a suitable f η = V k B T 0 C f (t) dt
38 Correlation Function Then, the goal is to show that C f (t) t e t/τ M Hence η τ M, which allows to interpret τ M as the Maxwell relaxation time.
39 Correlation Function The dissipative evolution operator P t acting on the set of functions f is defined by P t f (x) = E{ f (X(t)) X(0) = x} Then C f (t) = R 2 f (x) P t f (x) d 2 x
40 IV - Computing Maxwell s Time
41 Laplace Transform The Laplace transform of C(t) is defined by LC(ζ) = 0 e tζ C(t) dt The function C admits the asymptotic C(t) t e t/τ M if and only if LC(ζ) is holomorphic w.r.t. ζ in the domain Rζ > 1/τ M In practice, 1/τ M is the singularity nearest to the origin in the complex plane.
42 Dual Actions Phase Space Rotation: If f is a function, and x = (u, v) a point in the phase space, then f ( e ωt J x ) = ( e ωtj f ) (x), J = v u u v Hence J is the phase-space angular momentum. Phase Space Translation: Similarly, if ξ = (ξ 1, ξ 2 ) is a phasespace vector, f (x + ξ) = ( e ξ f ) (x), ξ = ξ 1 u + ξ 2 v Hence ξ is the phase-space momentum along the vector ξ.
43 Stochastic Evolution Let X(t) the stochastic value of the phase-space position at time t, with initial condition X(0) = x. If n anankeons occurred during this time then τ n t < τ n+1, with τ 0 = 0, so that f (X(t)) = n ( e (τ j τ j 1 )ω j 1 J e ξ ) j e (t τ n)ω n J f (x) j=1 In this expression the the τ j τ j 1 s, the ω j s and the ξ j s are random, independent and identically distributed.
44 Stochastic Evolution Averaging and taking the Laplace transform leads to LP ζ f (x) = 0 e tζ P t f (x) dt, ζ C Equivalently LP ζ f (x) = E n=0 τn+1 τ n e tζ n ( e (τ j τ j 1 )ω j 1 J e ) ξ j e (t τ n)ω n J f (x) dt j=1 Remark: τ 0 = 0, thus t = t τ n +(τ n τ n 1 )+ +(τ 1 τ 0 ). Hence ω j 1 J can be replaced by ω j 1 J + ζ and the exponential pre-factor e tζ disappears.
45 Stochastic Evolution First evaluate the integral over time between τ n and τ n+1, by setting s = t τ n τn+1 τn 0 e sζ e sω nj ds = 1 e(τ n+1 τ n )(ζ+ω n J) ζ + ω n J Second, average over the τ j τ j 1 s, using the formula (here A is an operator) E τ { e (τ j τ j 1 )A } = 0 e s/τ LC sa ds τ LC = τ LC A
46 Stochastic Evolution Reminder: if two random variables X, Y are stochastically independent then E{XY} = E{X} E{Y}. The average over the τ j τ j 1 s gives LP ζ f (x) = τ LC E n=0 n ( ) τ LC (ζ + ω j 1 J) eξ j j= τ LC (ζ + ω n J) f (x) This gives a new averaging over a product of independent variables.
47 Stochastic Evolution Averaging over the ξ j s can be done using (Gibbs average is a Gaussian integral) E ξ { e ξ j } = e k BT /2m, = = 2 u + 2 v Averaging over the ω j s leads to defining the following operator 1 A(ζ) = E 1 + τ LC (ζ + ω j 1 J) Remark: A(ζ) does not depend on j since all the ω j s have same distribution.
48 Stochastic Evolution Inserting into the expression of the Laplace transform leads to Questions: LP ζ f (x) = τ LC n=0 { A(ζ) e k B T /2m } n A(ζ) f (x) = τ LC 1 1 A(ζ) e k BT /2m A(ζ) f (x) How do we evaluate this function? How can we compute the domain of analyticity in ζ?
49 Angular Momentum Trick: the operator A(ζ) is a function of the phase-space angular momentum J. The polar coordinates in the 2D-phase space are given by It follows that { u = r cos θ v = r sin θ J = θ { r 2 = u 2 + v 2 tan θ = v/u
50 Angular Momentum Consequently the eigenvalues of J are given by ıl with l = 0, ±1, ±2,, namely l Z. The eigenfunctions have the form g l (r, θ) = ĝ l (r) e ılθ Projecting a function f onto the eigenspace of eigenvalue ıl is given by 2π Π l f (r, θ) = e ılθ f (r, θ) e ılθ dθ 2π Hence the spectral decomposition gives J = ı l Π l 0 l Z
51 Angular Momentum Any function of J can be written as F(J) = F(ıl) Π l l Z It leads to A(ζ) = a l (ζ) Π l l Z with a l (ζ) a complex number given by 1 a l (ζ) = E ω 1 + τ LC (ζ + ılω j )
52 Angular Momentum New Trick: the operator commutes with J, more precisely, the polar decomposition gives = p 2 r + J2 r 2 p 2 r = 2 r + 1 r r This gives 1 LP ζ = τ LC l Z 1 a l (ζ)e k BT(p 2 r +l2 /r 2 ) a l(ζ) Π l
53 Analyticity Remark that if ω > 0 the function (1 + τ LC (ζ + ılω)) 1 admits a pole at ζ = 1 τ LC ılω It follows that a l (ζ) is analytic in R(ζ) > 1/τ LC The operator p 2 r + l2 /r 2 is positive, and its spectrum is the entire positive real line. Hence, LP ζ is analytic in the domain for which there is no l Z nor any p 0 such that a l (ζ) = e p2 1. It means the set of ζ must satisfy a l (ζ) [1, )
54 Analyticity The domain of analyticity of a l (ζ) is given by Rζ > 1/τ LC
55 Analyticity The correlation function for X(t) involves only the angular momentum l = ±1. Assume that the distribution of ω j is uniform with average ω and variance σ. Then the ω j s are uniformly distributed in the interval [ω 3 σ, ω + 3 σ] and 1 a ±1 = 2ı 3τ LC σ ln 1 + τ LC(ζ + ıω + ı 3σ) 1 + τ LC (ζ + ıω ı 3σ) If ζ + ıω is not real, then the imaginary part of the r.h.s. is non zero, so that a ±1 [1, ).
56 Analyticity If ξ = ζ + ıω is real then 1 a ±1 = 2ı 3τ LC σ ln 1 + τ LC(ζ + ıω + ı 3σ) 1 + τ LC (ζ + ıω ı 3σ) = θ 3τLCσ where tan θ = 3τLCσ 1 + τ LC ξ θ < π 2
57 Maxwell Relaxation Time After some algebra, this gives τ M = τ LC (1 τ LC 3τLC σ tan 3τ LC σ if 3τLCσ π/2 ) 1 > τlc otherwise Using a uniform distribution of oscillator frequencies is a reasonable approximation. If τ LC σ decreases to zero as T 0, then This gives a crossover temperature T co above which the anankeon dominates and τ LC = τ M for T T co If T < T co, the phonons resist and τ M /τ LC > 1.
58 Maxwell Relaxation Time The variance σ of the random frequencies depend upon the modification of the local landscape by anankeons. However, each anankeon involves several atoms, at least d + 2 atoms in dimension d. So the landscape is modified in a region that might be large compared with the mean atomic distance. For this reason, σ is expected to be proportional to a Gibbs factor, namely to follow also an Arrhenius law σ = σ e W v/k B T σ = lim T σ(t) Question: What is the meaning of W v?
59 Maxwell Relaxation Time Similarly the local configuration time τ LC is also given by a similar expression, thanks to Kramers formula, where now W is the potential energy barrier to be crossed when an anankeon occurs τ LC = τ e W/k BT τ = lim T τ(t) If the hypothesis made on σ is correct, then K(T) = τ LC (T)σ(T) = K e (W v W)/k B T K = σ τ The condition K(T) T 0 0 requires W v > W
60 Maxwell Relaxation Time Hence, if T g < T < T co τ M τ LC = 1 1 3K(T) tan ( 3K(T)) T 0 e2(w v W)/k B T K 2 Similarly the crossover temperature is reached whenever 3τ LC σ = π/2, which gives (W v W) = k BT co 10.2
61 V - Conclusion
62 To Summarize The liquid phase of fragile glasses is dominated by anankeons at least above the crossover temperature T > T co. If T g < T < T co, the phonon-anankeon interaction becomes essential to explain the difference between the Maxwell and the local configuration times. Hence the change of the Arrhenius behavior of the viscosity is explained through a dynamical effect. As the temperature decreases, phonons become more coherent and limit the dissipative effect of the anankeons. The crossover temperature is related to the difference W s W between the activation energies associated with the phonon frequency fluctuation and the anankeon potential barrier.
arxiv: v1 [cond-mat.soft] 19 Aug 2017
Anankeon theory and viscosity of liquids: a toy model arxiv:1708.06624v1 [cond-mat.soft] 19 Aug 2017 Jean V. Bellissard Georgia Institute of Technology, School of Mathematics, School of Physics, Atlanta,
More informationHappy Birthday Yosi!
Happy Birthday Yosi! MODELING LIQUID METALS Sponsoring and BULK METALLIC GLASSES This material is based upon work supported by the National Science Foundation Grant No. DMS-1160962 Jean BELLISSARD Georgia
More informationA Simple Theory of Viscosity in Liquids
A Simple Theory of Viscosity in Liquids Jean Bellissard * School of Mathematics and School of Physics, Georgia Institute of Technology, Atlanta GA, USA, and Fachbereich Mathematik und Informatik, Westfälische
More informationMODELING LIQUID METALS BULK METALLIC GLASSES. and. Sponsoring. Jean BELLISSARD
Friedrich-Schiller-Universität Jena May 19, 2015 1 MODELING LIQUID METALS Sponsoring and BULK METALLIC GLASSES This material is based upon work supported by the National Science Foundation Grant No. DMS-1160962
More informationMODELING LIQUID METALS BULK METALLIC GLASSES. and. Sponsoring. Jean BELLISSARD
MODELING LIQUID METALS Sponsoring and BULK METALLIC GLASSES This material is based upon work supported by the National Science Foundation Grant No. DMS-1160962 Jean BELLISSARD Georgia Institute of Technology,
More informationWANNIER TRANSFORM APERIODIC SOLIDS. for. Jean BELLISSARD. Collaboration:
WANNIER TRANSFORM for APERIODIC SOLIDS Jean BELLISSARD Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics e-mail: jeanbel@math.gatech.edu Collaboration: G. DE NITTIS (SISSA,
More informationRIEMANNIAN GEOMETRY COMPACT METRIC SPACES. Jean BELLISSARD 1. Collaboration:
RIEMANNIAN GEOMETRY of COMPACT METRIC SPACES Jean BELLISSARD 1 Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics Collaboration: I. PALMER (Georgia Tech, Atlanta) 1 e-mail:
More informationVARIOUS MATHEMATICAL ASPECTS TILING SPACES. Jean BELLISSARD 1 2. Collaborations: Georgia Institute of Technology & Institut Universitaire de France
GaTech January 24 2005 1 VARIOUS MATHEMATICAL ASPECTS of TILING SPACES Jean BELLISSARD 1 2 Georgia Institute of Technology & Institut Universitaire de France Collaborations: D. SPEHNER (Essen, Germany)
More informationCOHOMOLOGY. Sponsoring. Jean BELLISSARD
Sponsoring Grant no. 0901514 COHOMOLOGY Jean BELLISSARD CRC 701, Bielefeld Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics e-mail: jeanbel@math.gatech.edu Main References
More informationLAPLACIANS COMPACT METRIC SPACES. Sponsoring. Jean BELLISSARD a. Collaboration:
LAPLACIANS on Sponsoring COMPACT METRIC SPACES Jean BELLISSARD a Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics Collaboration: I. PALMER (Georgia Tech, Atlanta) a e-mail:
More informationDISSIPATIVE TRANSPORT APERIODIC SOLIDS KUBO S FORMULA. and. Jean BELLISSARD 1 2. Collaborations:
Vienna February 1st 2005 1 DISSIPATIVE TRANSPORT and KUBO S FORMULA in APERIODIC SOLIDS Jean BELLISSARD 1 2 Georgia Institute of Technology, Atlanta, & Institut Universitaire de France Collaborations:
More informationThe Topology of Tiling Spaces
East Lansing October 16, 2009 1 The Topology of Tiling Spaces Jean BELLISSARD East Lansing October 16, 2009 2 East Lansing October 16, 2009 3 East Lansing October 16, 2009 4 The Topology of Tiling Spaces
More informationTHE TOPOLOGY of TILING SPACES
Seoul National University March20, 2014 1 THE TOPOLOGY of Sponsoring TILING SPACES This material is based upon work supported by the National Science Foundation Grant No. DMS-1160962 Jean BELLISSARD Georgia
More informationAperiodic Hamiltonians
Spectrum Approximation for Aperiodic Hamiltonians Sponsoring Jean BELLISSARD Westfälische Wilhelms-Universität, Münster Department of Mathematics Georgia Institute of Technology, Atlanta School of Mathematics
More informationNON-UNIVERSALITY OF THE NAZAROV-SODIN CONSTANT
NON-UNIVERSALITY OF THE NAZAROV-SODIN CONSTANT Abstract. We prove that the Nazarov-Sodin constant, which up to a natural scaling gives the leading order growth for the expected number of nodal components
More informationPeriodic Approximants. Aperiodic Hamiltonians
Periodic Approximants to Sponsoring Aperiodic Hamiltonians Jean BELLISSARD CRC 701, Bielefeld, Germany Westfälische Wilhelms-Universität, Münster Department of Mathematics Georgia Institute of Technology,
More informationNONCOMMUTATIVE. GEOMETRY of FRACTALS
AMS Memphis Oct 18, 2015 1 NONCOMMUTATIVE Sponsoring GEOMETRY of FRACTALS Jean BELLISSARD Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics e-mail: jeanbel@math.gatech.edu
More informationTheory of Aperiodic Solids:
Theory of Aperiodic Solids: Sponsoring from 1980 to present Jean BELLISSARD jeanbel@math.gatech.edu Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics Content 1. Aperiodic
More informationCOMPUTATIONAL. NONCOMMUTATIVE GEOMETRY The work of E. Prodan. Sponsoring. Jean BELLISSARD
COMPUTATIONAL Sponsoring NONCOMMUTATIVE GEOMETRY The work of E. Prodan This material is based upon work supported by the National Science Foundation Grant No. DMS-1160962 Jean BELLISSARD Georgia Institute
More informationA mathematical model for Mott s Variable Range Hopping
Santiago November 23-28, 2009 1 A mathematical model for Mott s Variable Range Hopping Jean BELLISSARD 1 2 Georgia Institute of Technology, Atlanta, http://www.math.gatech.edu/ jeanbel/ Collaborations:
More informationNodal lines of Laplace eigenfunctions
Nodal lines of Laplace eigenfunctions Spectral Analysis in Geometry and Number Theory on the occasion of Toshikazu Sunada s 60th birthday Friday, August 10, 2007 Steve Zelditch Department of Mathematics
More informationSymmetries 2 - Rotations in Space
Symmetries 2 - Rotations in Space This symmetry is about the isotropy of space, i.e. space is the same in all orientations. Thus, if we continuously rotated an entire system in space, we expect the system
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationDELONE SETS AND MATERIAL SCIENCE: A PROGRAM
DELONE SETS AND MATERIAL SCIENCE: A PROGRAM JEAN BELLISSARD Abstract. These notes are proposing a program liable to provide physicists working in material science, especially metallic liquids and glasses,
More informationTopics for the Qualifying Examination
Topics for the Qualifying Examination Quantum Mechanics I and II 1. Quantum kinematics and dynamics 1.1 Postulates of Quantum Mechanics. 1.2 Configuration space vs. Hilbert space, wave function vs. state
More informationF r (t) = 0, (4.2) F (t) = r= r a (t)
Chapter 4 Stochastic Equations 4.1 Langevin equation To explain the idea of introduction of the Langevin equation, let us consider the concrete example taken from surface physics: motion of an atom (adatom)
More informationQUANTUM TRANSPORT. KUBO S FORMULA for. Jean BELLISSARD 1 2. Collaborations:
Oberwolfach May 4th, 2004 1 KUBO S FORMULA for QUANTUM TRANSPORT Jean BELLISSARD 1 2 Georgia Institute of Technology, Atlanta, & Institut Universitaire de France Collaborations: H. SCHULZ-BALDES (T.U.
More informationLecture 17 Brownian motion as a Markov process
Lecture 17: Brownian motion as a Markov process 1 of 14 Course: Theory of Probability II Term: Spring 2015 Instructor: Gordan Zitkovic Lecture 17 Brownian motion as a Markov process Brownian motion is
More informationNonequilibrium transitions in glassy flows. Peter Schall University of Amsterdam
Nonequilibrium transitions in glassy flows Peter Schall University of Amsterdam Liquid or Solid? Liquid or Solid? Example: Pitch Solid! 1 day 1 year Menkind 10-2 10 0 10 2 10 4 10 6 10 8 10 10 10 12 10
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationJ10M.1 - Rod on a Rail (M93M.2)
Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight
More informationA model of alignment interaction for oriented particles with phase transition
A model of alignment interaction for oriented particles with phase transition Amic Frouvelle Archimedes Center for Modeling, Analysis & Computation (ACMAC) University of Crete, Heraklion, Crete, Greece
More informationLinear Systems Theory
ME 3253 Linear Systems Theory Review Class Overview and Introduction 1. How to build dynamic system model for physical system? 2. How to analyze the dynamic system? -- Time domain -- Frequency domain (Laplace
More informationStochastic nonlinear Schrödinger equations and modulation of solitary waves
Stochastic nonlinear Schrödinger equations and modulation of solitary waves A. de Bouard CMAP, Ecole Polytechnique, France joint work with R. Fukuizumi (Sendai, Japan) Deterministic and stochastic front
More informationPeriodic Approximations for Energy Spectra in Aperiodic Media
Periodic Approximations for Energy Spectra in Aperiodic Media Sponsoring Jean BELLISSARD Westfälische Wilhelms-Universität, Münster Department of Mathematics Georgia Institute of Technology, Atlanta School
More informationBranching Processes II: Convergence of critical branching to Feller s CSB
Chapter 4 Branching Processes II: Convergence of critical branching to Feller s CSB Figure 4.1: Feller 4.1 Birth and Death Processes 4.1.1 Linear birth and death processes Branching processes can be studied
More informationDynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology.
Dynamic Response Assoc. Prof. Enver Tatlicioglu Department of Electrical & Electronics Engineering Izmir Institute of Technology Chapter 3 Assoc. Prof. Enver Tatlicioglu (EEE@IYTE) EE362 Feedback Control
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More information2012 NCTS Workshop on Dynamical Systems
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,
More informationarxiv: v7 [quant-ph] 22 Aug 2017
Quantum Mechanics with a non-zero quantum correlation time Jean-Philippe Bouchaud 1 1 Capital Fund Management, rue de l Université, 75007 Paris, France. (Dated: October 8, 018) arxiv:170.00771v7 [quant-ph]
More informationPolymer Dynamics and Rheology
Polymer Dynamics and Rheology 1 Polymer Dynamics and Rheology Brownian motion Harmonic Oscillator Damped harmonic oscillator Elastic dumbbell model Boltzmann superposition principle Rubber elasticity and
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationPeriodic Approximant. Aperiodic Hamiltonians
Sponsoring Periodic Approximant to Aperiodic Hamiltonians CRC 701, Bielefeld, Germany Jean BELLISSARD Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics e-mail: jeanbel@math.gatech.edu
More informationTILINGS APERIODIC MEDIA
Duke February 28 2003 1 TILINGS APERIODIC MEDIA and their NONCOMMUTATIVE GEOMETRY Jean BELLISSARD 1 2 Georgia Institute of Technology & Institut Universitaire de France Collaborations: D. SPEHNER (Essen,
More informationBrownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539
Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory
More informationIn-class exercises Day 1
Physics 4488/6562: Statistical Mechanics http://www.physics.cornell.edu/sethna/teaching/562/ Material for Week 11 Exercises due Mon Apr 16 Last correction at April 16, 2018, 11:19 am c 2018, James Sethna,
More informationBloch Theory for 1D-FLC Aperiodic Media
Sponsoring Bloch Theory for 1D-FLC Aperiodic Media CRC 701, Bielefeld, Germany Jean BELLISSARD Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics e-mail: jeanbel@math.gatech.edu
More information1+t 2 (l) y = 2xy 3 (m) x = 2tx + 1 (n) x = 2tx + t (o) y = 1 + y (p) y = ty (q) y =
DIFFERENTIAL EQUATIONS. Solved exercises.. Find the set of all solutions of the following first order differential equations: (a) x = t (b) y = xy (c) x = x (d) x = (e) x = t (f) x = x t (g) x = x log
More informationA model of alignment interaction for oriented particles with phase transition
A model of alignment interaction for oriented particles with phase transition Amic Frouvelle ACMAC Joint work with Jian-Guo Liu (Duke University, USA) and Pierre Degond (Institut de Mathématiques de Toulouse,
More informationSupplementary Information for Observation of dynamic atom-atom correlation in liquid helium in real space
3 4 5 6 7 8 9 0 3 4 5 6 7 8 9 0 Supplementary Information for Observation of dynamic atom-atom correlation in liquid helium in real space Supplementary Note : Total PDF The total (snap-shot) PDF is obtained
More informationLINEAR RESPONSE THEORY
MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior
More informationT. Egami. Model System of Dense Random Packing (DRP)
Introduction to Metallic Glasses: How they are different/similar to other glasses T. Egami Model System of Dense Random Packing (DRP) Hard Sphere vs. Soft Sphere Glass transition Universal behavior History:
More informationIn terms of measures: Exercise 1. Existence of a Gaussian process: Theorem 2. Remark 3.
1. GAUSSIAN PROCESSES A Gaussian process on a set T is a collection of random variables X =(X t ) t T on a common probability space such that for any n 1 and any t 1,...,t n T, the vector (X(t 1 ),...,X(t
More informationChapter Three Theoretical Description Of Stochastic Resonance 24
Table of Contents List of Abbreviations and Symbols 5 Chapter One Introduction 8 1.1 The Phenomenon of the Stochastic Resonance 8 1.2 The Purpose of the Study 10 Chapter Two The Experimental Set-up 12
More informationOn a class of stochastic differential equations in a financial network model
1 On a class of stochastic differential equations in a financial network model Tomoyuki Ichiba Department of Statistics & Applied Probability, Center for Financial Mathematics and Actuarial Research, University
More information1. Stochastic Processes and filtrations
1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S
More informationThe glass transition as a spin glass problem
The glass transition as a spin glass problem Mike Moore School of Physics and Astronomy, University of Manchester UBC 2007 Co-Authors: Joonhyun Yeo, Konkuk University Marco Tarzia, Saclay Mike Moore (Manchester)
More information11. Recursion Method: Concepts
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 16 11. Recursion Method: Concepts Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative
More informationHomework 3 Solutions Problem 1 (a) The technique is essentially that of Homework 2, problem 2. The situation is depicted in the figure:
Homework 3 Solutions Problem (a) The technique is essentially that of Homework 2, problem 2. The situation is depicted in the figure: θ photon vdt A θ d Figure : The figure shows the system at time t.
More informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
More informationto 1D Periodic Approximants Aperiodic Hamiltonians Sponsoring Jean BELLISSARD CRC 701, Bielefeld, Germany
Sponsoring Periodic Approximants to 1D Aperiodic Hamiltonians CRC 701, Bielefeld, Germany Jean BELLISSARD Georgia Institute of Technology, Atlanta School of Mathematics & School of Physics e-mail: jeanbel@math.gatech.edu
More informationLocal time path integrals and their application to Lévy random walks
Local time path integrals and their application to Lévy random walks Václav Zatloukal (www.zatlovac.eu) Faculty of Nuclear Sciences and Physical Engineering Czech Technical University in Prague talk given
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
More informationNodal Sets of High-Energy Arithmetic Random Waves
Nodal Sets of High-Energy Arithmetic Random Waves Maurizia Rossi UR Mathématiques, Université du Luxembourg Probabilistic Methods in Spectral Geometry and PDE Montréal August 22-26, 2016 M. Rossi (Université
More informationNov : Lecture 22: Differential Operators, Harmonic Oscillators
14 Nov. 3 005: Lecture : Differential Operators, Harmonic Oscillators Reading: Kreyszig Sections:.4 (pp:81 83),.5 (pp:83 89),.8 (pp:101 03) Differential Operators The idea of a function as something that
More informationQuantum fields close to black hole horizons
Quantum fields close to black hole horizons Kinematics Accelerated scalar fields and inertial forces Photons in Rindler space vs thermal photons Interactions Static interactions of scalar, electric and
More informationHandout 10: Inductance. Self-Inductance and inductors
1 Handout 10: Inductance Self-Inductance and inductors In Fig. 1, electric current is present in an isolate circuit, setting up magnetic field that causes a magnetic flux through the circuit itself. This
More informationKadath: a spectral solver and its application to black hole spacetimes
Kadath: a spectral solver and its application to black hole spacetimes Philippe Grandclément Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris F-92195 Meudon, France philippe.grandclement@obspm.fr
More informationFundamental Inequalities, Convergence and the Optional Stopping Theorem for Continuous-Time Martingales
Fundamental Inequalities, Convergence and the Optional Stopping Theorem for Continuous-Time Martingales Prakash Balachandran Department of Mathematics Duke University April 2, 2008 1 Review of Discrete-Time
More informationRadiation by a dielectric wedge
Radiation by a dielectric wedge A D Rawlins Department of Mathematical Sciences, Brunel University, United Kingdom Joe Keller,Cambridge,2-3 March, 2017. We shall consider the problem of determining the
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationWeb Course Physical Properties of Glass. Range Behavior
Web Course Physical Properties of Glass Glass Transformation- Range Behavior Richard K. Brow Missouri University of Science & Technology Department of Materials Science & Engineering Glass Transformation-1
More informationIntroduction to Process Control
Introduction to Process Control For more visit :- www.mpgirnari.in By: M. P. Girnari (SSEC, Bhavnagar) For more visit:- www.mpgirnari.in 1 Contents: Introduction Process control Dynamics Stability The
More informationThe Equipartition Theorem
Chapter 8 The Equipartition Theorem Topics Equipartition and kinetic energy. The one-dimensional harmonic oscillator. Degrees of freedom and the equipartition theorem. Rotating particles in thermal equilibrium.
More information2 Lyapunov Exponent exponent number. Lyapunov spectrum for colour conductivity of 8 WCA disks.
3 2 Lyapunov Exponent 1-1 -2-3 5 1 15 exponent number Lyapunov spectrum for colour conductivity of 8 WCA disks. 2 3 Lyapunov exponent 2 1-1 -2-3 5 1 15 exponent number Lyapunov spectrum for shear flow
More informationFunctional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals
Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico
More informationChapter 6: The Rouse Model. The Bead (friction factor) and Spring (Gaussian entropy) Molecular Model:
G. R. Strobl, Chapter 6 "The Physics of Polymers, 2'nd Ed." Springer, NY, (1997). R. B. Bird, R. C. Armstrong, O. Hassager, "Dynamics of Polymeric Liquids", Vol. 2, John Wiley and Sons (1977). M. Doi,
More informationPHYSICAL SCIENCES PART A
PHYSICAL SCIENCES PART A 1. The calculation of the probability of excitation of an atom originally in the ground state to an excited state, involves the contour integral iωt τ e dt ( t τ ) + Evaluate the
More informationSlightly off-equilibrium dynamics
Slightly off-equilibrium dynamics Giorgio Parisi Many progresses have recently done in understanding system who are slightly off-equilibrium because their approach to equilibrium is quite slow. In this
More informationWhat is Q? Interpretation 1: Suppose A 0 represents wave amplitudes, then
What is Q? Interpretation 1: Suppose A 0 represents wave amplitudes, then A = A 0 e bt = A 0 e ω 0 t /(2Q ) ln(a) ln(a) = ln(a 0 ) ω0 t 2Q intercept slope t Interpretation 2: Suppose u represents displacement,
More informationON FRACTIONAL RELAXATION
Fractals, Vol. 11, Supplementary Issue (February 2003) 251 257 c World Scientific Publishing Company ON FRACTIONAL RELAXATION R. HILFER ICA-1, Universität Stuttgart Pfaffenwaldring 27, 70569 Stuttgart,
More informationSymmetry Properties and Exact Solutions of the Fokker-Planck Equation
Nonlinear Mathematical Physics 1997, V.4, N 1, 13 136. Symmetry Properties and Exact Solutions of the Fokker-Planck Equation Valery STOHNY Kyïv Polytechnical Institute, 37 Pobedy Avenue, Kyïv, Ukraïna
More informationInformation and Credit Risk
Information and Credit Risk M. L. Bedini Université de Bretagne Occidentale, Brest - Friedrich Schiller Universität, Jena Jena, March 2011 M. L. Bedini (Université de Bretagne Occidentale, Brest Information
More informationBessel Functions Michael Taylor. Lecture Notes for Math 524
Bessel Functions Michael Taylor Lecture Notes for Math 54 Contents 1. Introduction. Conversion to first order systems 3. The Bessel functions J ν 4. The Bessel functions Y ν 5. Relations between J ν and
More informationLECTURE 8: DYNAMICAL SYSTEMS 7
15-382 COLLECTIVE INTELLIGENCE S18 LECTURE 8: DYNAMICAL SYSTEMS 7 INSTRUCTOR: GIANNI A. DI CARO GEOMETRIES IN THE PHASE SPACE Damped pendulum One cp in the region between two separatrix Separatrix Basin
More informationThursday, August 4, 2011
Chapter 16 Thursday, August 4, 2011 16.1 Springs in Motion: Hooke s Law and the Second-Order ODE We have seen alrealdy that differential equations are powerful tools for understanding mechanics and electro-magnetism.
More informationSecond Order Systems
Second Order Systems independent energy storage elements => Resonance: inertance & capacitance trade energy, kinetic to potential Example: Automobile Suspension x z vertical motions suspension spring shock
More informationA system that is both linear and time-invariant is called linear time-invariant (LTI).
The Cooper Union Department of Electrical Engineering ECE111 Signal Processing & Systems Analysis Lecture Notes: Time, Frequency & Transform Domains February 28, 2012 Signals & Systems Signals are mapped
More informationSimple Harmonic Oscillator
Classical harmonic oscillator Linear force acting on a particle (Hooke s law): F =!kx From Newton s law: F = ma = m d x dt =!kx " d x dt + # x = 0, # = k / m Position and momentum solutions oscillate in
More informationLatent voter model on random regular graphs
Latent voter model on random regular graphs Shirshendu Chatterjee Cornell University (visiting Duke U.) Work in progress with Rick Durrett April 25, 2011 Outline Definition of voter model and duality with
More informationAsymptotics of minimax stochastic programs
Asymptotics of minimax stochastic programs Alexander Shapiro Abstract. We discuss in this paper asymptotics of the sample average approximation (SAA) of the optimal value of a minimax stochastic programming
More informationThe Quantum Hall Conductance: A rigorous proof of quantization
Motivation The Quantum Hall Conductance: A rigorous proof of quantization Spyridon Michalakis Joint work with M. Hastings - Microsoft Research Station Q August 17th, 2010 Spyridon Michalakis (T-4/CNLS
More informationExam 4 SCORE. MA 114 Exam 4 Spring Section and/or TA:
Exam 4 Name: Section and/or TA: Last Four Digits of Student ID: Do not remove this answer page you will return the whole exam. You will be allowed two hours to complete this test. No books or notes may
More informationThis is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or
Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects
More informationGeometric inequalities for black holes
Geometric inequalities for black holes Sergio Dain FaMAF-Universidad Nacional de Córdoba, CONICET, Argentina. 3 August, 2012 Einstein equations (vacuum) The spacetime is a four dimensional manifold M with
More informationCourse: ESO-209 Home Work: 1 Instructor: Debasis Kundu
Home Work: 1 1. Describe the sample space when a coin is tossed (a) once, (b) three times, (c) n times, (d) an infinite number of times. 2. A coin is tossed until for the first time the same result appear
More informationScale invariant fluid dynamics for the dilute Fermi gas at unitarity
Scale invariant fluid dynamics for the dilute Fermi gas at unitarity Thomas Schaefer North Carolina State University Fluids: Gases, Liquids, Plasmas,... Hydrodynamics: Long-wavelength, low-frequency dynamics
More information