Velocity distribution in active particles systems, Supplemental Material
|
|
- Janis Robbins
- 5 years ago
- Views:
Transcription
1 Velocity distribution in active particles systems, Supplemental Material Umberto Marini Bettolo Marconi, 1 Claudio Maggi, 2 Nicoletta Gnan, 2 Matteo Paoluzzi, 3 and Roberto Di Leonardo 2 1 Scuola di Scienze e Tecnologie, Università di Camerino, Via Madonna delle Carceri, 62032, Camerino, INFN Perugia, Italy 2 Universita Sapienza, Rome, Italy 3 Syracuse University, NY, USA (Dated: February 19, 2016 A simple derivation of the stationary velocity distribution of a collection of active particles subject to non equilibrium colored noise is provided in this note. Moreover simulation results for the 2- dimensional case are reported. I. DERIVATION OF THE STATIONARY VELOCITY DISTRIBUTION The eective dynamics for space coordinates of an assembly of active spheres 1,2 is ṙ i (t = 1 γ F i(r 1,..., r N + u i (t (1 where the velocities u i evolve according to the law: u i (t = 1 u i(t + D1/2 a η i (t (2 The force F i = i U acting on the i-th particle is conservative and associated to the potential U(r 1,..., r N, γ is the drag coecient, whereas the stochastic vectors η i (t are Gaussian and Marovian processes distributed with zero mean and moments η i (tη j (t = 2δ ij δ(t t. where d is the spatial dimensionality. The coecient D a due to the activity is related to the correlation of the Ornstein-Uhlenbec process u i (t via u i (tu j (t = d D a δ ij exp( t t /. where d is the spatial dimension. To simplify the notation we switch from r i to an array x i as done in previous publications: ẋ i (t = 1 γ F i(x 1,..., x N + u i (t (3 Dierentiate again ẍ i (t = 1 γ F i x ẋ + u i (t (4 Let us introduce the variable v i and recast eq.(5 as: ẍ i (t = 1 F i ẋ 1 γ x [ẋ i F i γ ] + D1/2 a η i (5 ẋ i = v i (6 v i = 1 F i v 1 γ x [v i F i γ ] + D1/2 a η i (7
2 2 A. One particle Before presenting the multidimensional result, we digress to illustrate the inetic method of solution in a simple one-dimensional case. We begin with a single particle in one dimension and drop the index i. We dierentiate eq. (1 with respect to time and introduce the velocity variable v = ẋ so that instead of the original system (1 and (2 we have: We obtain the Kramers equation for the phase-space distribution f(x, ẋ; t: v = 1 (1 df γ dx v + 1 γ F + D1/2 a η ẋ = v (8 f t + v f x U f γ v = 1 v ( D v + Γ(xv f (9 with Γ(x = 1 + γ 2 U(x. By multiplying and integrating over v eq. (9 and considering only time independent solutions f 0 (x, v one obtains: dvv 2 f 0 (x, v + U dvf 0 (x, v = 1 x γ Γ(x dvvf 0 (x, v (10 Such an integro-dierential equation can be solved by the following ( ansatz f 0 (x, v = Π(v x P(x where Π is the Gaussian velocity distribution : Π(v x = ( 2πD 1/2 Γ(x exp 2D Γ(xv2 whose width depends on the particle position and the average velocity v vanishes. After substituting the factorization f 0 into (9 and evaluating the velocity variance dvv 2 Π(v (x = D a 1 (11 Γ(x we arrive at the following dierential equation determining the steady state coordinate distribution P(x : d dx [Γ 1 (xp(x] + U(x P(x = 0 (12 Dγ which is identical to the dierential equation determining the stationary coordinate distribution in the unied color noise approximation (UCNA: with P(x = 1 Z 1 exp( H(x T s (13 H(x = U(x + 2γ (U(x x 2 T s ln[1 + γ d 2 U(x dx 2 ] (14 and T s = D a γ and Z 1 a normalization constant dened as: ( Z 1 = dx exp H(x. (15 T s The method can be easily extended to the multidimensional case. II. A KINETIC DERIVATION OF THE STEADY DISTRIBUTION, VELOCITY CORRELATIONS FOR MANY PARTICLE SYSTEMS We now generalize the inetic argument above illustrated to a many particle system and write the following multidimensional Kramers equation 3 describing the evolution of the phase-space distribution f N (x 1,..., x N ;, v 1..., v N ; t f N t + i v i f N x i + i F i f N γ v i = 1 i ( Da v i v i + Γ i v f N (16
3 3 with the non dimensional friction matrix Γ i dened as Γ i = δ i + γ 2 U x i x. (17 The probability conservation and the momentum balance equations are straightforwardly obtained by projection, i.e. by integrating eq.(16 over the dn dimensional velocity space after multiplying by 1 and v i, respectively: J i ({x i }, t t + P N ({x i }, t t + i J i ({x i }, t x i = 0 (18 p i ({x i }, t F i({x i }, t P N ({x i }, t = 1 Γ i ({x i }, tj ({x i }, t (19 x γ where P N ({x i }, t is the marginalized distribution giving the distribution of positions of the particles regardless their velocities: P N ({x i }, t = d N v f N ({x i }, ({v i }, t, (20 while the current J i ({x i }, t is the N-dimensional vector: J i ({x i }, t = d N v v i f N ({x i }, ({v i }, t. (21 Finally, p ij ({x i }, t is the N N dimensional generalized inetic pressure tensor p ij ({x i }, t = d N v v i v j f N ({x i }, ({v i }, t (22 Equations (18-(19 are an exact consequence of (16, but require the nowledge of p ij which can be obtained by continuing the projection procedure in velocity space to higher order in v i and closing the hierarchy by a suitable truncation ansatz, such as the introduction of phenomenological transport coecients. On the other hand, if we limit ourselves to study the stationary regime the approach is fruitful and leads to a dierent interpretation of the MUCNA equations. Let us rst, notice that the operator featuring in the right hand side of (16 posseses a null eigenvalue whose associated eigenfunction (the ground state is the multivariate Gaussian velocity distribution : Π(({v i } {x i } = ( N/2 ( det Γ exp v i Γ ij ({x}v j (23 2πD a 2D a In other words Π({v i } {x i } is the conditional probability distribution of velocities given the positions tae on the values {x i }. We now construct a time-independent trial phase-space distribution having a factorized form: ij f (K T ({x i }, ({v i } = Π(({v i } {x i } P st ({x i } (24 and corresponding to zero current J i (x 1,..., x N. For generic potentials the velocity moments of G are non constant in space, but depend on coordinates x i and dierent velocity components may be correlated and are given by the formula: v i v j = d N vv i v j Π({v i } {x i } = D a Γ 1 ij ({x i} (25 so that the "pressure " tensor reads p ij ({x i } = Da Γ 1 ij ({x i}p st ({x i }. Using eq. condition of vanishing currents, J i = 0, we have the following balance equations: D a (19 in conjunction with the p i ({x i } F i({x i } P st ({x i } = 0 (26 x D a γ Eliminating the pressure tensor using (25 we arrive at the dierential equation determining P st ({x i } : [Γ 1 i x ({x i}p st ({x i }] + 1 U({x i } P st ({x i } = 0 (27 D a γ x i
4 Finally, by multiplying by the matrix Γ equation (27 we arrive at an equation identical to the MUCNA equation (see ref 2 ( PN ( (D a + D t γ P N ln det Γ P N x n x n δ n + γ 2 U U(x1,..., x N = 0. (28 x n x x in the case D t = 0, so that we may identify the congurational part of the trial solution P st with P N derived by the MUCNA method 2 : P N (x 1,..., x N = 1 Z N { exp 1 [ U(x 1,..., x N + D a γ 2γ N ( U(x1,..., x N 2 ]} Da γ ln det Γ i x where Z N is a normalization constant. We turn now to the formula for the velocity covariance matrix for many particles. Using the Gaussianity of the distribution we can immediately write the averages of the velocity as < v i v j >= d N x v i v j P N ({x i } = D a 4 (29 d N xγ 1 ij ({x i}p N ({x i } (30 where stands for the double average over velocities and over space since the velocity distribution depends on the positions, unlie the equilibrium case. In the general case we must evaluate the matrix elements Γ ij. As shown in the appendix this can be done to rst order in /γ with the result for the velocity self-correlations: < v i αv i α >= D a d N rp st (r 1,..., r N [1 w αα (r i, r j ] (31 γ j i Finally, averaging over the whole system and introducing the positional pair correlation function g(r r one obtains 1 N N i=1 < vαv i α i >= D a (1 γ ρ drg 2 (rw αα (r (32 In analogy with granular gases we can dene the second moment of the velocity distribution to be the average inetic temperature, T of the system, via T = 1 N N i v iv i. III. APPROXIMATION FOR THE DETERMINANT AND VELOCITY CORRELATIONS The exact evaluation of the determinant Γ associated with the Hessian matrix is beyond the authors capabilities and we loo for approximations in order to evaluate the eective forces. We consider the associated determinant in the case of two spatial dimensions and vanishing external potential: [1 + γ wxx(r1, rj] wxy(r1, rj wxx(r1, r2... wxy(r1, rn γ γ γ wyx(r1, rj [1 + γ γ wyy(r1, rj] wyx(r1, r2... wyy(r1, rn γ γ wxx(r2, r1 wyx(r2, r1 [1 + γ γ γ j 2 wxx(r2, rj]... wxy(r2, rn γ γ wxy(rn, r1 γ wyy(rn, r [1 + γ j N wyy(rn, rj] It is interesting to remar that the o-diagonal elements contain only one term, while the diagonal elements and their neighbors contain N terms. Thus in the limit of N we expect that the matrix becomes eectively diagonal. However, even with such a limit, we see that to order /γ the inverse matrix, which is directly related to velocity
5 5 correlator v αi v βj (see (25, reads: [1 γ wxx(r1, rj] wxy(r1, rj wxx(r1, r2... wxy(r1, rn γ γ γ wyx(r1, rj [1 γ γ wyy(r1, rj] wyx(r1, r2... wyy(r1, rn γ γ wxx(r2, r1 wyx(r2, r1 D v [1 γ γ γ j 2 wxx(r2, rj]... wxy(r2, rn γ αiv βj wxy(rn, r1 γ γ wyy(rn, r [1 γ j N wyy(rn, rj] We, now, assume translational invariance and perform the average over the spatial distribution P N (r 1,..., r N, and introduce the pair correlation as P 2 (r 1, r 2 = ρ2 N(N 1 g(r 1 r 2. The sought average is: D < v αiv βj >= 1 ρ γ xx(rg(r ρ γ xy(rg(r ρ Nγ xx(rg(r... ρ Nγ xy(rg(r ρ γ yx(rg(r 1 ρ γ yy(rg(r ρ Nγ yx(r... ρ Nγ yy(rg(r ρ rw Nγ xx(rg(r ρ Nγ yx(rg(r 1 ρ γ xx(rg(r... ρ Nγ xy(rg(r One can see that in the limit of large N the o-diagonal terms become neglegible because they carry factors 1/N and one obtains a bloc diagonal matrix whose blocs are 2 2 matrices: ( < vxv x > < v xv y > = D ( 1 ρ γ xx(rg(r ρ γ xy(rg(r < v yv x > < v yv y > yx(rg(r 1 ρ γ yy(rg(r If the system is rotationally invariant < v x v x >=< v y v y > and < v x v y >=< v y v x >= 0. ρ γ A. Evaluation of the ensemble average of the determinant to linear order in /γ and velocity variance Notice that to linear order /γ the ensemble average of the determinant is: dr 1... dr N P N (r 1,..., r N det Γ(r 1,..., r N ([1 + γ ρ xx (rg(r][1 + γ ρ yy (rg(r] N (33 dr 1... dr N P N (r 1,..., r N det Γ(r 1,..., r N (1 + γ ρ dr[w xx (r + w yy (r]g(r N (34 To order /γ: dr 1... dr N P N (r 1,..., r N ln det Γ N ln (1+ γ ρ ( D dr[w xx (r+w yy (r]g(r N ln 1 ( D +N ln < v x v x > 1 < v y v y > (35 IV. 2D SIMULATION RESULTS We simulate a 2d system driven by GCN and composed by N = 1000 particles with periodic boundary conditions. The particles interact via the purely repulsive pair potential ϕ(r 1 r 2 = r 1 r We x = 1 and we explore several values of density ρ = N/L 2 and diusivity D. The results are reported in Fig. 1. Fig. 1(a shows the normalized velocity variance as a function of density for dierent values of D. It is seen that, as in the 1d case, the quantity ẋ 2 /(D/ is a decreasing function of ρ. However in the 1d case ẋ 2 /(D/ is also a decreasing function of D (at xed ρ, while in 2d a non-monotonic behavior in D can be observed at high ρ. Fig. 1(b shows the amplitude of the Fourier-transformed density uctuations at low q (q (20σ 1 which shows a very similar behaviour to the 1d case. Fig. 1(c and (d show two snapshots at high D and intermediate densities where evident clustering of the particles is observed but not a full phase separation. This is the case also upon further increasing to very high values as shown in Fig. 2.
6 6 Figure 1: (a Velocity variance as a function of density at dierent values of D. (b The colormap represents the values of the structure factor at low q for several values of the density and the diusivity. (c and (d Snapshot of the system at dierent values of D and ρ, an evident clustering is observed. 1 C. Maggi, U. M. B. Marconi, N. Gnan, and R. Di Leonardo, Scientic Reports ( U. M. B. Marconi and C. Maggi, Soft matter ( H. Risen, Foer-Planc Equation (Springer, 1984.
7 Figure 2: Snapshot of the system at high for two dierent values of D and ρ, an evident clustering is observed without a full phase separation. 7
Statistical Mechanics of Active Matter
Statistical Mechanics of Active Matter Umberto Marini Bettolo Marconi University of Camerino, Italy Naples, 24 May,2017 Umberto Marini Bettolo Marconi (2017) Statistical Mechanics of Active Matter 2017
More informationMultidimensional Stationary Probability Distribution for Interacting Active Particles, Supplemental Material
Multidimensional Stationary Probability istribution for Interacting Active Particles, Supplemental Material Claudio Maggi 1, a Umberto Marini Bettolo Marconi 2, Nicoletta Gnan 3, and Roberto i Leonardo
More informationSoft Matter Accepted Manuscript
Accepted Manuscript This is an Accepted Manuscript, which has been through the Royal Society of Chemistry peer review process and has been accepted for publication. Accepted Manuscripts are published online
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationDynamics that trigger/inhibit cluster formation in a one-dimensional granular gas
Physica A 342 (24) 62 68 www.elsevier.com/locate/physa Dynamics that trigger/inhibit cluster formation in a one-dimensional granular gas Jose Miguel Pasini a;1, Patricio Cordero b; ;2 a Department of Theoretical
More informationEEL 5544 Noise in Linear Systems Lecture 30. X (s) = E [ e sx] f X (x)e sx dx. Moments can be found from the Laplace transform as
L30-1 EEL 5544 Noise in Linear Systems Lecture 30 OTHER TRANSFORMS For a continuous, nonnegative RV X, the Laplace transform of X is X (s) = E [ e sx] = 0 f X (x)e sx dx. For a nonnegative RV, the Laplace
More informationMultiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions
Multiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions March 6, 2013 Contents 1 Wea second variation 2 1.1 Formulas for variation........................
More informationLesson 4: Stationary stochastic processes
Dipartimento di Ingegneria e Scienze dell Informazione e Matematica Università dell Aquila, umberto.triacca@univaq.it Stationary stochastic processes Stationarity is a rather intuitive concept, it means
More informationAsymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks
Physica A 275 (2000) 483 504 www.elsevier.com/locate/physa Asymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks V. Kumaran Department of Chemical Engineering, Indian
More information16. Working with the Langevin and Fokker-Planck equations
16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationNotes on Density Functional Theory
Notes on Density Functional Theory Rocco Martinazzo E-mail: rocco.martinazzo@unimi.it Contents 1 Introduction 1 Density Functional Theory 7 3 The Kohn-Sham approach 11 1 Introduction We consider here a
More informationU U B U P x
Smooth Quantum Hydrodynamic Model Simulation of the Resonant Tunneling Diode Carl L. Gardner and Christian Ringhofer y Department of Mathematics Arizona State University Tempe, AZ 8587-184 Abstract Smooth
More informationp D q EA Figure 3: F as a function of p d plotted for p d ' q EA, t = 0:05 and y =?0:01.
F p D q EA Figure 3: F as a function of p d plotted for p d ' q EA, t = :5 and y =?:. 5 F p A Figure : F as a function of p d plotted for < p d < 2p A, t = :5 and y =?:. F p A p max p p p D q EA Figure
More informationarxiv: v1 [cond-mat.soft] 31 Jan 2017
arxiv:7.932v [cond-mat.soft] 3 Jan 27 Effective equilibrium states in the colored-noise model for active matter I. Pairwise forces in the and Unified Colored Noise Approximations René Wittmann,, a C. Maggi,
More informationobtained is essentially depends on the degeneracy of one-electron levels (on the degree of symmetry). The lower symmetry, the larger secular equation
FRENKEL-EXCITON IN FINITE-GROUP-LATTICE Stan F.Harlapenko, Slava V. Rotkin Ioe Physico{Technical Institute, 26, Politekhnicheskaya st., 194021 St.Petersburg, Russia. E-mail: ROTKIN@THEORY.IOFFE.RSSI.RU
More informationMODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione
MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione, Univ. di Roma Tor Vergata, via di Tor Vergata 11,
More informationLinear Regression and Its Applications
Linear Regression and Its Applications Predrag Radivojac October 13, 2014 Given a data set D = {(x i, y i )} n the objective is to learn the relationship between features and the target. We usually start
More informationBROWNIAN DYNAMICS SIMULATIONS WITH HYDRODYNAMICS. Abstract
BROWNIAN DYNAMICS SIMULATIONS WITH HYDRODYNAMICS Juan J. Cerdà 1 1 Institut für Computerphysik, Pfaffenwaldring 27, Universität Stuttgart, 70569 Stuttgart, Germany. (Dated: July 21, 2009) Abstract - 1
More informationcorresponding velocity eld to consider in the Fokker-Planck equation shows singularities in the locations of the nodes of the wave function. These sin
Exact solutions of Fokker-Planck equations associated to quantum wave functions Nicola Cufaro Petroni, Salvatore De Martino and Silvio De Siena y ABSTRACT: We analyze the non stationary solutions of the
More informationUNIVERSITY OF MARYLAND Department of Physics College Park, Maryland. PHYSICS Ph.D. QUALIFYING EXAMINATION PART II
UNIVERSITY OF MARYLAND Department of Physics College Park, Maryland PHYSICS Ph.D. QUALIFYING EXAMINATION PART II January 22, 2016 9:00 a.m. 1:00 p.m. Do any four problems. Each problem is worth 25 points.
More information5 Applying the Fokker-Planck equation
5 Applying the Fokker-Planck equation We begin with one-dimensional examples, keeping g = constant. Recall: the FPE for the Langevin equation with η(t 1 )η(t ) = κδ(t 1 t ) is = f(x) + g(x)η(t) t = x [f(x)p
More informationParallel Additive Gaussian Channels
Parallel Additive Gaussian Channels Let us assume that we have N parallel one-dimensional channels disturbed by noise sources with variances σ 2,,σ 2 N. N 0,σ 2 x x N N 0,σ 2 N y y N Energy Constraint:
More informationB. Physical Observables Physical observables are represented by linear, hermitian operators that act on the vectors of the Hilbert space. If A is such
G25.2651: Statistical Mechanics Notes for Lecture 12 I. THE FUNDAMENTAL POSTULATES OF QUANTUM MECHANICS The fundamental postulates of quantum mechanics concern the following questions: 1. How is the physical
More informationChapter 6 - Random Processes
EE385 Class Notes //04 John Stensby Chapter 6 - Random Processes Recall that a random variable X is a mapping between the sample space S and the extended real line R +. That is, X : S R +. A random process
More informationRandom Matrix Eigenvalue Problems in Probabilistic Structural Mechanics
Random Matrix Eigenvalue Problems in Probabilistic Structural Mechanics S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
More informationwhere (E) is the partition function of the uniform ensemble. Recalling that we have (E) = E (E) (E) i = ij x (E) j E = ij ln (E) E = k ij ~ S E = kt i
G25.265: Statistical Mechanics Notes for Lecture 4 I. THE CLASSICAL VIRIAL THEOREM (MICROCANONICAL DERIVATION) Consider a system with Hamiltonian H(x). Let x i and x j be specic components of the phase
More informationThe Time{Dependent Born{Oppenheimer Approximation and Non{Adiabatic Transitions
The Time{Dependent Born{Oppenheimer Approximation and Non{Adiabatic Transitions George A. Hagedorn Department of Mathematics, and Center for Statistical Mechanics, Mathematical Physics, and Theoretical
More informationLecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions
Lecture Models for heavy-ion collisions (Part III: transport models SS06: Dynamical models for relativistic heavy-ion collisions Quantum mechanical description of the many-body system Dynamics of heavy-ion
More informationin order to insure that the Liouville equation for f(?; t) is still valid. These equations of motion will give rise to a distribution function f(?; t)
G25.2651: Statistical Mechanics Notes for Lecture 21 Consider Hamilton's equations in the form I. CLASSICAL LINEAR RESPONSE THEORY _q i = @H @p i _p i =? @H @q i We noted early in the course that an ensemble
More informationn i,j+1/2 q i,j * qi+1,j * S i+1/2,j
Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/TERMO-5-97 DATE: December 9,997 TITLE A comparison of complete vs. simplied viscous terms in boundary layer
More informationGaussian random variables inr n
Gaussian vectors Lecture 5 Gaussian random variables inr n One-dimensional case One-dimensional Gaussian density with mean and standard deviation (called N, ): fx x exp. Proposition If X N,, then ax b
More informationBrownian Motion: Fokker-Planck Equation
Chapter 7 Brownian Motion: Fokker-Planck Equation The Fokker-Planck equation is the equation governing the time evolution of the probability density of the Brownian particla. It is a second order differential
More informationLoop Corrections in Non-Linear Cosmological Perturbation Theory. Roman Scoccimarro
Loop Corrections in Non-Linear Cosmological Perturbation Theory Roman Scoccimarro Department of Physics and Enrico Fermi Institute, University of Chicago, Chicago, IL 6067 and NASA/Fermilab Astrophysics
More informationViscous Fluids. Amanda Meier. December 14th, 2011
Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass density, velocity and pressure. An Eulerian description of uids focuses on the transport
More informationLecture 2. Spring Quarter Statistical Optics. Lecture 2. Characteristic Functions. Transformation of RVs. Sums of RVs
s of Spring Quarter 2018 ECE244a - Spring 2018 1 Function s of The characteristic function is the Fourier transform of the pdf (note Goodman and Papen have different notation) C x(ω) = e iωx = = f x(x)e
More informationIterative procedure for multidimesional Euler equations Abstracts A numerical iterative scheme is suggested to solve the Euler equations in two and th
Iterative procedure for multidimensional Euler equations W. Dreyer, M. Kunik, K. Sabelfeld, N. Simonov, and K. Wilmanski Weierstra Institute for Applied Analysis and Stochastics Mohrenstra e 39, 07 Berlin,
More information11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.
C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a
More informationTransformation of Probability Densities
Transformation of Probability Densities This Wikibook shows how to transform the probability density of a continuous random variable in both the one-dimensional and multidimensional case. In other words,
More informationAperture Photometry in Band Limited. Data. Robert Lupton. March 4, Band limited data, that is, data that has non-zero Fourier components in
Aperture Photometry in Band Limited Data Robert Lupton March 4, 996 Band Limited Data Band limited data, that is, data that has non-zero Fourier components in only a restricted part of k-space jkj < k,
More informationDynamic System Identification using HDMR-Bayesian Technique
Dynamic System Identification using HDMR-Bayesian Technique *Shereena O A 1) and Dr. B N Rao 2) 1), 2) Department of Civil Engineering, IIT Madras, Chennai 600036, Tamil Nadu, India 1) ce14d020@smail.iitm.ac.in
More informationBifurcation analysis of incompressible ow in a driven cavity F.W. Wubs y, G. Tiesinga z and A.E.P. Veldman x Abstract Knowledge of the transition point of steady to periodic ow and the frequency occurring
More informationMultiple Random Variables
Multiple Random Variables This Version: July 30, 2015 Multiple Random Variables 2 Now we consider models with more than one r.v. These are called multivariate models For instance: height and weight An
More informationFokker-Planck Equation with Detailed Balance
Appendix E Fokker-Planck Equation with Detailed Balance A stochastic process is simply a function of two variables, one is the time, the other is a stochastic variable X, defined by specifying: a: the
More informationRandom Eigenvalue Problems Revisited
Random Eigenvalue Problems Revisited S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
More informationMachine Learning. 7. Logistic and Linear Regression
Sapienza University of Rome, Italy - Machine Learning (27/28) University of Rome La Sapienza Master in Artificial Intelligence and Robotics Machine Learning 7. Logistic and Linear Regression Luca Iocchi,
More informationThe Derivative. Appendix B. B.1 The Derivative of f. Mappings from IR to IR
Appendix B The Derivative B.1 The Derivative of f In this chapter, we give a short summary of the derivative. Specifically, we want to compare/contrast how the derivative appears for functions whose domain
More informationDynamical Systems & Scientic Computing: Homework Assignments
Fakultäten für Informatik & Mathematik Technische Universität München Dr. Tobias Neckel Summer Term Dr. Florian Rupp Exercise Sheet 3 Dynamical Systems & Scientic Computing: Homework Assignments 3. [ ]
More informationThe Gaussian distribution
The Gaussian distribution Probability density function: A continuous probability density function, px), satisfies the following properties:. The probability that x is between two points a and b b P a
More informationBayesian decision theory Introduction to Pattern Recognition. Lectures 4 and 5: Bayesian decision theory
Bayesian decision theory 8001652 Introduction to Pattern Recognition. Lectures 4 and 5: Bayesian decision theory Jussi Tohka jussi.tohka@tut.fi Institute of Signal Processing Tampere University of Technology
More informationPlan of Class 4. Radial Basis Functions with moving centers. Projection Pursuit Regression and ridge. Principal Component Analysis: basic ideas
Plan of Class 4 Radial Basis Functions with moving centers Multilayer Perceptrons Projection Pursuit Regression and ridge functions approximation Principal Component Analysis: basic ideas Radial Basis
More informationActive Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: Hydrodynamics of SP Hard Rods
Active Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: of SP Hard Rods M. Cristina Marchetti Syracuse University Baskaran & MCM, PRE 77 (2008);
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 1 Defining statistical mechanics: Statistical Mechanics provies the connection between microscopic motion of individual atoms of matter and macroscopically
More informationp(z)
Chapter Statistics. Introduction This lecture is a quick review of basic statistical concepts; probabilities, mean, variance, covariance, correlation, linear regression, probability density functions and
More informationThe Multivariate Gaussian Distribution
The Multivariate Gaussian Distribution Chuong B. Do October, 8 A vector-valued random variable X = T X X n is said to have a multivariate normal or Gaussian) distribution with mean µ R n and covariance
More informationCHAPTER V. Brownian motion. V.1 Langevin dynamics
CHAPTER V Brownian motion In this chapter, we study the very general paradigm provided by Brownian motion. Originally, this motion is that a heavy particle, called Brownian particle, immersed in a fluid
More informationOn prediction and density estimation Peter McCullagh University of Chicago December 2004
On prediction and density estimation Peter McCullagh University of Chicago December 2004 Summary Having observed the initial segment of a random sequence, subsequent values may be predicted by calculating
More informationPhysics 403. Segev BenZvi. Parameter Estimation, Correlations, and Error Bars. Department of Physics and Astronomy University of Rochester
Physics 403 Parameter Estimation, Correlations, and Error Bars Segev BenZvi Department of Physics and Astronomy University of Rochester Table of Contents 1 Review of Last Class Best Estimates and Reliability
More informationECONOMICS 207 SPRING 2006 LABORATORY EXERCISE 5 KEY. 8 = 10(5x 2) = 9(3x + 8), x 50x 20 = 27x x = 92 x = 4. 8x 2 22x + 15 = 0 (2x 3)(4x 5) = 0
ECONOMICS 07 SPRING 006 LABORATORY EXERCISE 5 KEY Problem. Solve the following equations for x. a 5x 3x + 8 = 9 0 5x 3x + 8 9 8 = 0(5x ) = 9(3x + 8), x 0 3 50x 0 = 7x + 7 3x = 9 x = 4 b 8x x + 5 = 0 8x
More informationFrom Hamiltonian particle systems to Kinetic equations
From Hamiltonian particle systems to Kinetic equations Università di Roma, La Sapienza WIAS, Berlin, February 2012 Plan of the lectures 1 Particle systems and BBKGY hierarchy (the paradigm of the Kinetic
More informationChapter 5. Sound Waves and Vortices. 5.1 Sound waves
Chapter 5 Sound Waves and Vortices In this chapter we explore a set of characteristic solutions to the uid equations with the goal of familiarizing the reader with typical behaviors in uid dynamics. Sound
More informationOrdering periodic spatial structures by non-equilibrium uctuations
Physica A 277 (2000) 327 334 www.elsevier.com/locate/physa Ordering periodic spatial structures by non-equilibrium uctuations J.M.G. Vilar a;, J.M. Rub b a Departament de F sica Fonamental, Facultat de
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More informationBMIR Lecture Series on Probability and Statistics Fall 2015 Discrete RVs
Lecture #7 BMIR Lecture Series on Probability and Statistics Fall 2015 Department of Biomedical Engineering and Environmental Sciences National Tsing Hua University 7.1 Function of Single Variable Theorem
More informationActive and Driven Soft Matter Lecture 3: Self-Propelled Hard Rods
A Tutorial: From Langevin equation to Active and Driven Soft Matter Lecture 3: Self-Propelled Hard Rods M. Cristina Marchetti Syracuse University Boulder School 2009 A Tutorial: From Langevin equation
More informationMixture Models & EM. Nicholas Ruozzi University of Texas at Dallas. based on the slides of Vibhav Gogate
Mixture Models & EM icholas Ruozzi University of Texas at Dallas based on the slides of Vibhav Gogate Previously We looed at -means and hierarchical clustering as mechanisms for unsupervised learning -means
More informationThe effect of mutual angular misalignment in the quantized sliding of solid lubricants
Facoltà di Scienze e Tecnologie Laurea Triennale in Fisica The effect of mutual angular misalignment in the quantized sliding of solid lubricants Relatore: Prof. Nicola Manini Correlatore: Prof. Rosario
More informationThus, the volume element remains the same as required. With this transformation, the amiltonian becomes = p i m i + U(r 1 ; :::; r N ) = and the canon
G5.651: Statistical Mechanics Notes for Lecture 5 From the classical virial theorem I. TEMPERATURE AND PRESSURE ESTIMATORS hx i x j i = kt ij we arrived at the equipartition theorem: * + p i = m i NkT
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationAppendix 5. Derivatives of Vectors and Vector-Valued Functions. Draft Version 24 November 2008
Appendix 5 Derivatives of Vectors and Vector-Valued Functions Draft Version 24 November 2008 Quantitative genetics often deals with vector-valued functions and here we provide a brief review of the calculus
More informationLecture 15. Theory of random processes Part III: Poisson random processes. Harrison H. Barrett University of Arizona
Lecture 15 Theory of random processes Part III: Poisson random processes Harrison H. Barrett University of Arizona 1 OUTLINE Poisson and independence Poisson and rarity; binomial selection Poisson point
More informationOne-Loop Integrals at vanishing external momenta and appli. Higgs potentials reconstructions CALC Mikhail Dolgopolov. Samara State University
One-Loop Integrals at vanishing external momenta and applications for extended Higgs potentials reconstructions CALC 2012 Samara State University 1. Basic examples and an idea for extended Higgs potentials
More informationquantum semiconductor devices. The model is valid to all orders of h and to rst order in the classical potential energy. The smooth QHD equations have
Numerical Simulation of the Smooth Quantum Hydrodynamic Model for Semiconductor Devices Carl L. Gardner and Christian Ringhofer y Department of Mathematics Arizona State University Tempe, AZ 8587-84 Abstract
More informationMixture Models & EM. Nicholas Ruozzi University of Texas at Dallas. based on the slides of Vibhav Gogate
Mixture Models & EM icholas Ruozzi University of Texas at Dallas based on the slides of Vibhav Gogate Previously We looed at -means and hierarchical clustering as mechanisms for unsupervised learning -means
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationToday: Fundamentals of Monte Carlo
Today: Fundamentals of Monte Carlo What is Monte Carlo? Named at Los Alamos in 940 s after the casino. Any method which uses (pseudo)random numbers as an essential part of the algorithm. Stochastic - not
More information5. Random Vectors. probabilities. characteristic function. cross correlation, cross covariance. Gaussian random vectors. functions of random vectors
EE401 (Semester 1) 5. Random Vectors Jitkomut Songsiri probabilities characteristic function cross correlation, cross covariance Gaussian random vectors functions of random vectors 5-1 Random vectors we
More information3. Probability and Statistics
FE661 - Statistical Methods for Financial Engineering 3. Probability and Statistics Jitkomut Songsiri definitions, probability measures conditional expectations correlation and covariance some important
More informationImplicit sampling for particle filters. Alexandre Chorin, Mathias Morzfeld, Xuemin Tu, Ethan Atkins
0/20 Implicit sampling for particle filters Alexandre Chorin, Mathias Morzfeld, Xuemin Tu, Ethan Atkins University of California at Berkeley 2/20 Example: Try to find people in a boat in the middle of
More informationSimultaneous Robust Design and Tolerancing of Compressor Blades
Simultaneous Robust Design and Tolerancing of Compressor Blades Eric Dow Aerospace Computational Design Laboratory Department of Aeronautics and Astronautics Massachusetts Institute of Technology GTL Seminar
More informationDependence. Practitioner Course: Portfolio Optimization. John Dodson. September 10, Dependence. John Dodson. Outline.
Practitioner Course: Portfolio Optimization September 10, 2008 Before we define dependence, it is useful to define Random variables X and Y are independent iff For all x, y. In particular, F (X,Y ) (x,
More informationLecture 5: Moment Generating Functions
Lecture 5: Moment Generating Functions IB Paper 7: Probability and Statistics Carl Edward Rasmussen Department of Engineering, University of Cambridge February 28th, 2018 Rasmussen (CUED) Lecture 5: Moment
More informationARTIFICIAL INTELLIGENCE LABORATORY. and CENTER FOR BIOLOGICAL INFORMATION PROCESSING. A.I. Memo No August Federico Girosi.
MASSACHUSETTS INSTITUTE OF TECHNOLOGY ARTIFICIAL INTELLIGENCE LABORATORY and CENTER FOR BIOLOGICAL INFORMATION PROCESSING WHITAKER COLLEGE A.I. Memo No. 1287 August 1991 C.B.I.P. Paper No. 66 Models of
More informationA Probability Review
A Probability Review Outline: A probability review Shorthand notation: RV stands for random variable EE 527, Detection and Estimation Theory, # 0b 1 A Probability Review Reading: Go over handouts 2 5 in
More information7 To solve numerically the equation of motion, we use the velocity Verlet or leap frog algorithm. _ V i n = F i n m i (F.5) For time step, we approxim
69 Appendix F Molecular Dynamics F. Introduction In this chapter, we deal with the theories and techniques used in molecular dynamics simulation. The fundamental dynamics equations of any system is the
More informationModelling high-frequency economic time series
Physica A 288 (2000) 444 450 www.elsevier.com/locate/physa Modelling high-frequency economic time series Lei-Han Tang a;, Zhi-Feng Huang a;b a Department of Physics and Center for Nonlinear Studies, Hong
More informationTHE USE OF ROBUST OBSERVERS IN THE SIMULATION OF GAS SUPPLY NETWORKS
UNIVERSITY OF READING DEPARTMENT OF MATHEMATICS THE USE OF ROBUST OBSERVERS IN THE SIMULATION OF GAS SUPPLY NETWORKS by SIMON M. STRINGER This thesis is submitted forthe degree of Doctorof Philosophy NOVEMBER
More informationTurbulence and Bose-Einstein condensation Far from Equilibrium
Turbulence and Bose-Einstein condensation Far from Equilibrium Dénes Sexty Uni Heidelberg Collaborators: Jürgen Berges, Sören Schlichting March, 01, Seattle Gauge Field Dynamics In and Out of Equilibrium
More informationNational Sun Yat-Sen University CSE Course: Information Theory. Maximum Entropy and Spectral Estimation
Maximum Entropy and Spectral Estimation 1 Introduction What is the distribution of velocities in the gas at a given temperature? It is the Maxwell-Boltzmann distribution. The maximum entropy distribution
More informationStochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno
Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.
More information/95 $ $.25 per page
Fields Institute Communications Volume 00, 0000 McGill/95-40 gr-qc/950063 Two-Dimensional Dilaton Black Holes Guy Michaud and Robert C. Myers Department of Physics, McGill University Montreal, Quebec,
More informationPY 351 Modern Physics - Lecture notes, 3
PY 351 Modern Physics - Lecture notes, 3 Copyright by Claudio Rebbi, Boston University, October 2016. These notes cannot be duplicated and distributed without explicit permission of the author. Time dependence
More information8.1 Concentration inequality for Gaussian random matrix (cont d)
MGMT 69: Topics in High-dimensional Data Analysis Falll 26 Lecture 8: Spectral clustering and Laplacian matrices Lecturer: Jiaming Xu Scribe: Hyun-Ju Oh and Taotao He, October 4, 26 Outline Concentration
More information04. Random Variables: Concepts
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 215 4. Random Variables: Concepts Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationPolarons. University of Ljubljana Faculty of Mathematics and Physics. Seminar - Ib, 1. year II. cycle degrees
University of Ljubljana Faculty of Mathematics and Physics Seminar - Ib, 1. year II. cycle degrees Polarons Author: Jaka Vodeb Advisor: prof. dr. Viktor Kabanov Kranj, 017 Abstract After the discovery
More informationSecond-order approximation of dynamic models without the use of tensors
Second-order approximation of dynamic models without the use of tensors Paul Klein a, a University of Western Ontario First draft: May 17, 2005 This version: January 24, 2006 Abstract Several approaches
More informationA priori tests of a stochastic mode reduction strategy
Physica D 170 (2002) 206 252 A priori tests of a stochastic mode reduction strategy A. Majda a,b, I. Timofeyev a,b, E. Vanden-Eijnden a, a Courant Institute of Mathematical Sciences, New Yor University,
More informationMultivariate Distribution Models
Multivariate Distribution Models Model Description While the probability distribution for an individual random variable is called marginal, the probability distribution for multiple random variables is
More informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
More information