Tutorial on the smooth particle-mesh Ewald algorithm
|
|
- George Andrews
- 5 years ago
- Views:
Transcription
1 Tutorial on the smooth particle-mesh Ewald algorithm These notes explain the smooth particle-mesh Ewald algorithm (SPME) in detail. The accompanying library libpme6 is a fully-featured implementation of SPME that can be used for experimentation and reference. This content can also be found in PDF format here. Introduction Motivation and overview of the SPME algorithm Consider two atoms, labeled and, contained in a parallelepiped with sides determined by the (positively oriented) vectors l, l, l 3 R 3 (see Figure ). Then the matrix L [l, l, l 3 ] R 3 3 is such Figure : Simulation box spanned by three vectors. that V det L > 0. ow suppose that atom i is at a position r i L x i where x i [0, ) 3 for i, and that r ij r i r j (see Figure ). Between these two atoms there is a force that depends on Figure : Two atoms at a distance r ij. their distance and whose corresponding potential energy, U ij U ij (r ij ), is shown in Figure. If instead of just two atoms we had a total of > atoms (see Figure ), then the total potential energy would be equal to H i j i j U ij (r ij ). () By solving ewton s second law we would be able to simulate the Force equals mass times acceleration.
2 tutorial on spme U(r) r r 6 U 0 (r) r U (r) r 6 Figure 3: Decomposition of the Lennard Jones -6 potential into repulsive (powers of ) and dispersive (powers of 6) contributions. Energy Distance (r) Figure 4: Argon gas ( 0 3 atoms). molecular system of atoms that interact according to the (pairwise) forces given by F U. In practice, whatever value of that can be tractable with current computers is still far away from the more realistic figure of A (Avogadro s number) which would capture the real physics of the system. One trick to somehow get closer to A is to use a tractable value of and impose periodic boundary conditions (PBCs). This means that we now have an infinite lattice of simulation boxes with identical copies of the particles in each box. By computing the interactions between the atoms in the original box and their copies (only the copies in the neighboring simulation boxes are typically taken into account this is called the minimum image convention ), we can make our simulations more realistic.
3 tutorial on spme 3 Figure 5: The protein DHFR solvated in water (left) and two views (center, right) without solvent. This system has 3536 atoms, so 08 and log 05. The total potential energy () then becomes H Uij (kri r j k) + i j i6 j Uij (kri r j + Lnk). n Z3 n 6 0 i j The naive approach to PBCs requires O( ) calculations and is unfeasible, hence the interest of finding efficient ways of carrying out this task. Ewald summation takes advantage of ideas from Harmonic analysis to reduce the complexity to O( 3/ ) operations, which is still impractical if we want to solve the equations of motion of a large system. The Particle-Mesh Ewald (PME) method 3 uses the Fast Fourier Transform to bring the complexity down to O( log ), and this algorithm has been a staple of Molecular Simulation since its appearance. The Smooth Particle-Mesh Ewald (SPME) improves on PME by giving a sufficiently smooth energy function whose derivative can be obtained analytically, which results in more realistic simulations due to improved energy conservation. Some of the people involved in the development of the SPME algorithm appear in Figure. P. P. Ewald. Die Berechnung optischer und elektrostatischer Gitterpotentiale. Annalen der Physik, 369(3):53 87, 9 T. Darden, D. York, and L. Pedersen. Particle mesh Ewald: An log( ) method for Ewald sums in large systems. The Journal of Chemical Physics, 98():0089, Smooth Particle-Mesh Ewald We derive the formulas for Ewald summation following 4 and then we introduce the smooth Particle-Mesh Ewald method from 5. Decomposition of the potential energy function Ewald summation can be applied whenever the pairwise potential energy U U (r ) is proportional to r n for some n (a concrete instance of U appears in Figure ). In that case, we can write U as a sum of the form U (r ) U0 (r ) + U (r ), () D. E. Williams. Accelerated convergence treatment of R n lattice sums. In Uri Shmueli, editor, International Tables for Crystallography, volume B, pages Kluwer Academic Publishers, U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, H. Lee, and L. G. Pedersen. A smooth particle mesh Ewald method. The Journal of Chemical Physics, 03(9):8577, 995 4
4 tutorial on spme 4 Figure 6: Some of the people behind SPME. Pictures reproduced from 0/0/00/wx0006/ and psc.edu/science/darden/. (a) Paul Peter Ewald (b) Tom Darden and ( ). Lee Pedersen. (c) Mike Crowley. where the short-range term, U 0, decays fast as r goes to infinity (so that it can be truncated) and the long-range term, U, is very spread out and smooth, so that its Fourier transform is concentrated (and can also be truncated). There are many ways of splitting the potential energy but 6 establishes the optimality for a class of functions involving the (incomplete) Gamma function. If we define 6 B. R. A. ijboer and F. W. De Wette. On the calculation of lattice sums. Physica, 3(-5):309 3, 957 ϕ n (r) Γ(n/, α π r ), Γ(n/) where α is a parameter to be discussed later, then we clearly have r n r n (ϕ n(r) + ( ϕ n (r))). We are going to concern ourselves with the case n 6 (see Figure ), and to define ψ(r) r 6 ϕ 6 (r). The free parameter α is determined by numerically solving the equation r 6 ( ϕ 6 (r cutoff ; α)) ε with prescribed values of r cutoff and ε (the function ϕ n must be decreasing for r > r cutoff ). This makes the short-range part of the computation involving effectively vanish for r > r cutoff.
5 tutorial on spme 5 0. r 6 ϕ 6 (r) r 6 ( ϕ 6 (r)) r 6 Figure 7: Decomposition of the function r 6 into short-range and long-range terms. Energy r Recall that in the case of > particles, we have H i j i j U ij ( r i r j ) + n Z 3 i j n 0 U ij ( r i r j + Ln ). where r ij is the distance between atoms i and j. We can split the pairwise potential U ij above as in (). Ewald summation The summation of the long-range terms (with b i > 0), U i j i j can be rearranged as U b i b j ψ( r i r j ) + n Z 3 i j where ψ(0) lim r 0 + ψ(r). n Z 3 i j n 0 b i b j ψ( r i r j + Ln ) b i b j ψ( r i r j + Ln ) (3) bi ψ(0), (4) i Lemma (Poisson s summation formula). If D(r) n Z 3 δ(r Ln), then D(h) V d Z 3 δ(h L T d). Lemma. Consider the distributions ρ(r) i b i δ(r r i ), η(r) ρ( r), and D(r) n Z 3 δ(r Ln). The following identity holds: n Z 3 i j b i b j ψ( r i r j + Ln ) ψ(r) (ρ η D) (r) dr, where r r.
6 tutorial on spme 6 D Figure 8: Example of Poisson s summation formula D Proof. First, ρ η (r) ρ(s) η(r s) ds ρ(s) ρ(s r) ds i j i j b i b j δ(s r i ) δ(r s r j ) ds b i b j δ(r (r i r j )). By the associative property of the convolution, we have so Finally, ρ η D (ρ η) D, ρ η D (r) ρ η (s) D(s r) ds i j i j ψ(r) (ρ η D) (r) dr ψ(r) b i b j n Z 3 i j n Z 3 i j δ(s (r i r j )) D(s r) ds b i b j D(s (r i r j )) b i b j δ(s (r i r j + Ln)). b i b j δ(s (r i r j + Ln)) dr n Z 3 i j b i b j ψ( r i r j + Ln ).
7 tutorial on spme 7 At this point we can go back to (4) and rewrite it as U ψ(r) (ρ η D) (r) dr bi ψ(0), (5) i Furthermore, using Plancherel s Theorem 7 and the fact that ρ η D is real-valued, we can write ψ(r) (ρ η D) (r) dr ψ(h) (ρ η D)(h) dh, (6) 7 T. H. Wolff. Lectures on harmonic analysis, volume 9 of University Lecture Series. American Mathematical Society, Providence, RI, 003 where h h. Consequently, equation (5) becomes U ψ(h) (ρ η D)(h) dh bi ψ(0), (7) i By the properties of the Fourier transform of a convolution 8 and Poisson s summation formula, (ρ η D)(h) ρ(h) ρ( h) V d Z 3 δ(h L T d) ρ(h) V d Z 3 δ(h L T d). 8 T. H. Wolff. Lectures on harmonic analysis, volume 9 of University Lecture Series. American Mathematical Society, Providence, RI, 003 From the above and (6), we have ψ (ρ η D) ψ(h) ( ρ(h) ) V δ(h L T d) dh d Z 3 V d Z 3 ψ( L T d ) ρ(l T d) V m V m 0 ψ(m) ρ(m) ψ(m) ρ(m) + V lim ψ(m) ρ(0) m 0 V m 0 ψ(m) ρ(m) + V lim ψ(m) m 0 ( ) b i, i where m L T d and m m. So, to sum up, we have arrived at the identity: U i j i j b i b j ψ( r i r j ) + n Z 3 i j n 0 b i b j ψ( r i r j + Ln ) V ψ(m) ρ(m) ( ) + V ψ(0) b i ψ(0) bi, m 0 i i
8 tutorial on spme 8 Fast Fourier Transform and Particle-Mesh method Ewald summation using optimal parameters requires O( 3/ ) operations 9 but it can be modified so that it involves only about O( log ) operations by using the Fast Fourier Transform 0. We now concern ourselves with summing (and later taking derivatives of) the expression ψ(m) ρ(m). (8) m 0 By the definition of the (inverse) Fourier Transform, we see that ρ(m) b j e πi mt r j, j where r j L x j for some x j [0, ) 3. Moreover, 9 S. W. de Leeuw, J. W. Perram, and E. R. Smith. Simulation of Electrostatic Systems in Periodic Boundary Conditions. I. Lattice Sums and Dielectric Constants. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 373(75):7 56, October U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, H. Lee, and L. G. Pedersen. A smooth particle mesh Ewald method. The Journal of Chemical Physics, 03(9):8577, 995; and R. W. Hockney and J. J. W. Eastwood. Computer Simulation Using Particles. Adam Hilger, 988 ρ(m) b j e πi mt r j j b j e πi dt L L x j j b j e πi dt x j. j We are going to approximate the complex exponential terms above by means of Cardinal B-Splines. An order p B-spline function is the p-fold convolution of the uniform density in the unit interval (see Figure ). By construction, each order p B-spline is of class C p, their I. J. Schoenberg. Cardinal Spline Interpolation. Society for Industrial and Applied Mathematics, 973 Order Order Figure 9: Basic B-splines y 0 y x x Order 3 Order 4 y 0 y x x integral equals, and they have compact support. Remark. Observe that these splines play the role of approximations to the identity when used in a meshing of the unit cell. It is in this sense that they play a role in approximating the Dirac delta functions in ρ.
9 tutorial on spme 9 ow, to discretize the simulation box we consider K, K, K 3 and let [ ] K K K. For each r Lx, we set u Kx and we have e πi mtr e πi dtx e πi dt K u πi dtk l W(u l) e l Z 3 k k k K 3 n Z 3 W(u (k + Kn)) e πi dt K (k+kn) πi dtn W(u (k + Kn)) e πi dtk k e }{{} n Z 3 n Z 3 W(u (k + Kn)) e πi dt K k, where we have used the Euclidean algorithm to write l k + Kn with k i {0,,..., K i } and n Z 3. We also define the forward (unnormalized) DFT of the 3D array Q k as K K K 3 Q d DFT[Q k ] d Q k e πi dt K k. k 0 k 0 k 3 0 And, likewise, the backward DFT is IDFT[ Q d ] k K d 0 K d 0 Using the previous results, we write ρ(m) b j e πi dt x j j k j K 3 d 3 0 Q d e πi dt K k. b j W(u (k + Kn)) n Z } 3 {{} Q k Q k e πi dt K k Q d. k e πi dt K k Remark. It is important to notice that since the splines have compact support, the infinite series above is actually a finite sum. Remark. The following diagram depicts the relationships between the different objects involved in the SPME method: IFT ρ Q IDFT ρ Q
10 tutorial on spme 0 Going back to (8), we carry out the following approximation: ψ(m) ρ(m) ψ(m) ρ( m) ρ(m) ψ( L T d ) Q d Q d m 0 m 0 d 0 ow define with ϕ( L T h ), if h 0, A d 0, if h 0. d i, if 0 d i < h i (d) K i, K i d i, if K i d i < K i. (note that K i h i < K i ). Then, K d K K d K K 3 ψ( L T d ) Q d Q d d 3 K 3 K d 0 K d 0 K 3 d 3 0 A d Q d Q d. Applying the following property of the Fourier transform (see ), f dg ĝ d f, T. H. Wolff. Lectures on harmonic analysis, volume 9 of University Lecture Series. American Mathematical Society, Providence, RI, 003 we write (using the notation of the previous section and keeping in mind that the result works in the distributional sense too) ) ψ η ρ ψ η ρ (Ĉ ρ ρ where we have used that f g f ĝ and that, for f f (x), f (x) f ( x) (again, see 3 ). Therefore, the algorithm boils down to: Q k Q d DFT(Q) d 3 T. H. Wolff. Lectures on harmonic analysis, volume 9 of University Lecture Series. American Mathematical Society, Providence, RI, 003 A d DFT(Q) d IDFT(A DFT(Q)) k (IDFT(A) Q) k (IDFT(A) Q) k Q k, Assorted results in Harmonic analysis Theorem 3 (The convolution is commutative). Theorem 4 (The Fourier transform of a convolution is the product of the Fourier transforms). Theorem 5. 4 For absolutely continuous measures µ, ν in R n. We have ˆµ dν ˆν dµ 4 T. H. Wolff. Lectures on harmonic analysis, volume 9 of University Lecture Series. American Mathematical Society, Providence, RI, 003 Theorem 6 (Plancherel s theorem).
11 tutorial on spme References [dlps80] S. W. de Leeuw, J. W. Perram, and E. R. Smith. Simulation of Electrostatic Systems in Periodic Boundary Conditions. I. Lattice Sums and Dielectric Constants. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 373(75):7 56, October 980. [DYP93] T. Darden, D. York, and L. Pedersen. Particle mesh Ewald: An log() method for Ewald sums in large systems. The Journal of Chemical Physics, 98():0089, 993. [EPB + 95] U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, H. Lee, and L. G. Pedersen. A smooth particle mesh Ewald method. The Journal of Chemical Physics, 03(9):8577, 995. [Ewa] P. P. Ewald. Die Berechnung optischer und elektrostatischer Gitterpotentiale. Annalen der Physik, 369(3):53 87, 9. [HE88] R. W. Hockney and J. J. W. Eastwood. Computer Simulation Using Particles. Adam Hilger, 988. [D57] B. R. A. ijboer and F. W. De Wette. On the calculation of lattice sums. Physica, 3(-5):309 3, 957. [Sch73] I. J. Schoenberg. Cardinal Spline Interpolation. Society for Industrial and Applied Mathematics, 973. [Wil06] D. E. Williams. Accelerated convergence treatment of R n lattice sums. In Uri Shmueli, editor, International Tables for Crystallography, volume B, pages Kluwer Academic Publishers, 006. [Wol03] T. H. Wolff. Lectures on harmonic analysis, volume 9 of University Lecture Series. American Mathematical Society, Providence, RI, 003.
Non-bonded interactions
speeding up the number-crunching continued Marcus Elstner and Tomáš Kubař December 3, 2013 why care? number of individual pair-wise interactions bonded interactions proportional to N: O(N) non-bonded interactions
More informationNon-bonded interactions
speeding up the number-crunching Marcus Elstner and Tomáš Kubař May 8, 2015 why care? key to understand biomolecular structure and function binding of a ligand efficiency of a reaction color of a chromophore
More informationAdvanced Molecular Molecular Dynamics
Advanced Molecular Molecular Dynamics Technical details May 12, 2014 Integration of harmonic oscillator r m period = 2 k k and the temperature T determine the sampling of x (here T is related with v 0
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationNeighbor Tables Long-Range Potentials
Neighbor Tables Long-Range Potentials Today we learn how we can handle long range potentials. Neighbor tables Long-range potential Ewald sums MSE485/PHY466/CSE485 1 Periodic distances Minimum Image Convention:
More informationN-body simulations. Phys 750 Lecture 10
N-body simulations Phys 750 Lecture 10 N-body simulations Forward integration in time of strongly coupled classical particles or (rigid-body composite objects) dissipation Brownian motion Very general
More informationLecture 12: Solvation Models: Molecular Mechanics Modeling of Hydration Effects
Statistical Thermodynamics Lecture 12: Solvation Models: Molecular Mechanics Modeling of Hydration Effects Dr. Ronald M. Levy ronlevy@temple.edu Bare Molecular Mechanics Atomistic Force Fields: torsion
More informationComputational Molecular Biophysics. Computational Biophysics, GRS Jülich SS 2013
Computational Molecular Biophysics Computational Biophysics, GRS Jülich SS 2013 Computational Considerations Number of terms for different energy contributions / E el ~N ~N ~N ~N 2 Computational cost for
More informationNotes on Ewald summation techniques
February 3, 011 Notes on Ewald summation techniques Adapted from similar presentation in PHY 71 he total electrostatic potential energy of interaction between point charges {q i } at the positions {r i
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationLong-range interactions: P 3 M, MMMxD, ELC, MEMD and ICC
Long-range interactions: P 3 M, MMMxD, ELC, MEMD and ICC Axel Arnold Institute for Computational Physics Universität Stuttgart October 10, 2012 Long-range interactions gravity Coulomb interaction dipolar
More informationMOLECULAR DYNAMICS SIMULATIONS OF BIOMOLECULES: Long-Range Electrostatic Effects
Annu. Rev. Biophys. Biomol. Struct. 1999. 28:155 79 Copyright c 1999 by Annual Reviews. All rights reserved MOLECULAR DYNAMICS SIMULATIONS OF BIOMOLECULES: Long-Range Electrostatic Effects Celeste Sagui
More informationWhy Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.
Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. 25, 2002 Molecular Dynamics: Introduction At physiological conditions, the
More informationSmoothed Particle Hydrodynamics (SPH) 4. May 2012
Smoothed Particle Hydrodynamics (SPH) 4. May 2012 Calculating density SPH density estimator Weighted summation over nearby particles: ρ(r) = N neigh b=1 m bw (r r b, h) W weight function with dimension
More informationOptimizing GROMACS for parallel performance
Optimizing GROMACS for parallel performance Outline 1. Why optimize? Performance status quo 2. GROMACS as a black box. (PME) 3. How does GROMACS spend its time? (MPE) 4. What you can do What I want to
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationAn FPGA Implementation of Reciprocal Sums for SPME
An FPGA Implementation of Reciprocal Sums for SPME Sam Lee and Paul Chow Edward S. Rogers Sr. Department of Electrical and Computer Engineering University of Toronto Objectives Accelerate part of Molecular
More informationIntermolecular Forces and Monte-Carlo Integration 열역학특수연구
Intermolecular Forces and Monte-Carlo Integration 열역학특수연구 2003.3.28 Source of the lecture note. J.M.Prausnitz and others, Molecular Thermodynamics of Fluid Phase Equiliria Atkins, Physical Chemistry Lecture
More informationNonlinear Stationary Subdivision
Nonlinear Stationary Subdivision Michael S. Floater SINTEF P. O. Box 4 Blindern, 034 Oslo, Norway E-mail: michael.floater@math.sintef.no Charles A. Micchelli IBM Corporation T.J. Watson Research Center
More informationarxiv: v2 [physics.comp-ph] 19 Oct 2008
P 3 M algorithm for dipolar interactions Juan J. Cerdà, 1, V. Ballenegger, 2 O. Lenz, 1, 3 and C. Holm 1, 3 1 Frankfurt Inst. for Advanced Studies, (FIAS), arxiv:0805.4783v2 [physics.comp-ph] 19 Oct 2008
More informationSeminar in Particles Simulations. Particle Mesh Ewald. Theodora Konstantinidou
Seminar in Particles Simulations Particle Mesh Ewald Theodora Konstantinidou This paper is part of the course work: seminar 2 Contents 1 Introduction 3 1.1 Role of Molecular Dynamics - MD.......................
More informationEconometrics I. Lecture 10: Nonparametric Estimation with Kernels. Paul T. Scott NYU Stern. Fall 2018
Econometrics I Lecture 10: Nonparametric Estimation with Kernels Paul T. Scott NYU Stern Fall 2018 Paul T. Scott NYU Stern Econometrics I Fall 2018 1 / 12 Nonparametric Regression: Intuition Let s get
More informationThe Sommerfeld Polynomial Method: Harmonic Oscillator Example
Chemistry 460 Fall 2017 Dr. Jean M. Standard October 2, 2017 The Sommerfeld Polynomial Method: Harmonic Oscillator Example Scaling the Harmonic Oscillator Equation Recall the basic definitions of the harmonic
More information1 Introduction. 2 Measure theoretic definitions
1 Introduction These notes aim to recall some basic definitions needed for dealing with random variables. Sections to 5 follow mostly the presentation given in chapter two of [1]. Measure theoretic definitions
More informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
More informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
More informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
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 informationCoarse-Grained Models!
Coarse-Grained Models! Large and complex molecules (e.g. long polymers) can not be simulated on the all-atom level! Requires coarse-graining of the model! Coarse-grained models are usually also particles
More informationLong range interactions in computer simulations. 1 (Part 3)
Long range interactions in computer simulations. 1 (Part 3) Laboratoire de Physique Théorique (UMR-867), Université Paris-Sud 11 et CNRS, Bâtiment 10, 91405 Orsay Cedex, France. (.@th.u-psud.fr) 1 st Workshop
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationWe have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma and (3.55) with j = 1. We can write any f W 1 as
88 CHAPTER 3. WAVELETS AND APPLICATIONS We have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma 3..7 and (3.55) with j =. We can write any f W as (3.58) f(ξ) = p(2ξ)ν(2ξ)
More informationTHE MATTILA INTEGRAL ASSOCIATED WITH SIGN INDEFINITE MEASURES. Alex Iosevich and Misha Rudnev. August 3, Introduction
THE MATTILA INTEGRAL ASSOCIATED WITH SIGN INDEFINITE MEASURES Alex Iosevich and Misha Rudnev August 3, 005 Abstract. In order to quantitatively illustrate the role of positivity in the Falconer distance
More informationA new method for the solution of scattering problems
A new method for the solution of scattering problems Thorsten Hohage, Frank Schmidt and Lin Zschiedrich Konrad-Zuse-Zentrum Berlin, hohage@zibde * after February 22: University of Göttingen Abstract We
More informationPhysics 342 Lecture 30. Solids. Lecture 30. Physics 342 Quantum Mechanics I
Physics 342 Lecture 30 Solids Lecture 30 Physics 342 Quantum Mechanics I Friday, April 18th, 2008 We can consider simple models of solids these highlight some special techniques. 30.1 An Electron in a
More informationMolecular Modeling -- Lecture 15 Surfaces and electrostatics
Molecular Modeling -- Lecture 15 Surfaces and electrostatics Molecular surfaces The Hydrophobic Effect Electrostatics Poisson-Boltzmann Equation Electrostatic maps Electrostatic surfaces in MOE 15.1 The
More informationA review: The Laplacian and the d Alembertian. j=1
Chapter One A review: The Laplacian and the d Alembertian 1.1 THE LAPLACIAN One of the main goals of this course is to understand well the solution of wave equation both in Euclidean space and on manifolds
More informationLecture Notes on ELECTROMAGNETIC FIELDS AND WAVES
Lecture Notes on ELECTROMAGNETIC FIELDS AND WAVES (227-0052-10L) Prof. Dr. Lukas Novotny ETH Zürich, Photonics Laboratory January 12, 2018 Introduction The properties of electromagnetic fields and waves
More informationAn introduction to Molecular Dynamics. EMBO, June 2016
An introduction to Molecular Dynamics EMBO, June 2016 What is MD? everything that living things do can be understood in terms of the jiggling and wiggling of atoms. The Feynman Lectures in Physics vol.
More informationAnalytic families of multilinear operators
Analytic families of multilinear operators Mieczysław Mastyło Adam Mickiewicz University in Poznań Nonlinar Functional Analysis Valencia 17-20 October 2017 Based on a joint work with Loukas Grafakos M.
More informationEfficient multiple time step method for use with Ewald and particle mesh Ewald for large biomolecular systems
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 5 1 AUGUST 2001 Efficient multiple time step method for use with Ewald and particle mesh Ewald for large biomolecular systems Ruhong Zhou IBM Thomas J. Watson
More informationPreliminaries: what you need to know
January 7, 2014 Preliminaries: what you need to know Asaf Pe er 1 Quantum field theory (QFT) is the theoretical framework that forms the basis for the modern description of sub-atomic particles and their
More informationTOOLS FROM HARMONIC ANALYSIS
TOOLS FROM HARMONIC ANALYSIS BRADLY STADIE Abstract. The Fourier transform can be thought of as a map that decomposes a function into oscillatory functions. In this paper, we will apply this decomposition
More informationEfficient Algorithms for Electrostatic Interactions Including Dielectric Contrasts
Entropy 2013, 15, 4569-4588; doi:10.3390/e15114569 Article OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Efficient Algorithms for Electrostatic Interactions Including Dielectric Contrasts
More informationHarmonic Analysis Homework 5
Harmonic Analysis Homework 5 Bruno Poggi Department of Mathematics, University of Minnesota November 4, 6 Notation Throughout, B, r is the ball of radius r with center in the understood metric space usually
More informationAll-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010
All-atom Molecular Mechanics Trent E. Balius AMS 535 / CHE 535 09/27/2010 Outline Molecular models Molecular mechanics Force Fields Potential energy function functional form parameters and parameterization
More informationA reciprocal space based method for treating long range interactions in ab initio and force-field-based calculations in clusters
A reciprocal space based method for treating long range interactions in ab initio and force-field-based calculations in clusters Glenn J. Martyna and Mark E. Tuckerman Citation: The Journal of Chemical
More informationABSTRACT SCATTERING SYSTEMS: SOME SURPRISING CONNECTIONS
ABSTRACT SCATTERING SYSTEMS: SOME SURPRISING CONNECTIONS Cora Sadosky Department of Mathematics Howard University Washington, DC, USA csadosky@howard.edu Department of Mathematics University of New Mexico
More informationTrace Class Operators and Lidskii s Theorem
Trace Class Operators and Lidskii s Theorem Tom Phelan Semester 2 2009 1 Introduction The purpose of this paper is to provide the reader with a self-contained derivation of the celebrated Lidskii Trace
More informationBIOC : Homework 1 Due 10/10
Contact information: Name: Student # BIOC530 2012: Homework 1 Due 10/10 Department Email address The following problems are based on David Baker s lectures of forces and protein folding. When numerical
More informationHOW TO LOOK AT MINKOWSKI S THEOREM
HOW TO LOOK AT MINKOWSKI S THEOREM COSMIN POHOATA Abstract. In this paper we will discuss a few ways of thinking about Minkowski s lattice point theorem. The Minkowski s Lattice Point Theorem essentially
More informationQuadratic reciprocity (after Weil) 1. Standard set-up and Poisson summation
(September 17, 010) Quadratic reciprocity (after Weil) Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ I show that over global fields (characteristic not ) the quadratic norm residue
More informationLecture 27 Frequency Response 2
Lecture 27 Frequency Response 2 Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/6/12 1 Application of Ideal Filters Suppose we can generate a square wave with a fundamental period
More informationCHAPTER 6 bis. Distributions
CHAPTER 6 bis Distributions The Dirac function has proved extremely useful and convenient to physicists, even though many a mathematician was truly horrified when the Dirac function was described to him:
More informationReproducing formulas associated with symbols
Reproducing formulas associated with symbols Filippo De Mari Ernesto De Vito Università di Genova, Italy Modern Methods of Time-Frequency Analysis II Workshop on Applied Coorbit space theory September
More informationAccelerated Stokesian Dynamics simulations
J. Fluid Mech. (2001), vol. 448, pp. 115 146. c 2001 Cambridge University Press DOI: 10.1017/S0022112001005912 Printed in the United Kingdom 115 Accelerated Stokesian Dynamics simulations By ASIMINA SIEROU
More informationFree energy simulations
Free energy simulations Marcus Elstner and Tomáš Kubař January 14, 2013 Motivation a physical quantity that is of most interest in chemistry? free energies Helmholtz F or Gibbs G holy grail of computational
More informationCS 542G: The Poisson Problem, Finite Differences
CS 542G: The Poisson Problem, Finite Differences Robert Bridson November 10, 2008 1 The Poisson Problem At the end last time, we noticed that the gravitational potential has a zero Laplacian except at
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationFundamental Solutions of Stokes and Oseen Problem in Two Spatial Dimensions
J. math. fluid mech. c 6 Birkhäuser Verlag, Basel DOI.7/s-5-9-z Journal of Mathematical Fluid Mechanics Fundamental Solutions of Stokes and Oseen Problem in Two Spatial Dimensions Ronald B. Guenther and
More informationA universal thin film model for Ginzburg-Landau energy with dipolar interaction
A universal thin film model for Ginzburg-Landau energy with dipolar interaction Cyrill B. Muratov Dedicated to V. V. Osipov on the occasion of his 77th birthday Abstract We present an analytical treatment
More informationA Banach Gelfand Triple Framework for Regularization and App
A Banach Gelfand Triple Framework for Regularization and Hans G. Feichtinger 1 hans.feichtinger@univie.ac.at December 5, 2008 1 Work partially supported by EUCETIFA and MOHAWI Hans G. Feichtinger hans.feichtinger@univie.ac.at
More informationGreen s Function for Problems with Layers
Chapter 9 Green s Function for Problems with Layers 9. Continuity conditions The medium comprises two or more straight or curved layers. The properties of the medium are constant within each layer; they
More information1. Poisson-Boltzmann 1.1. Poisson equation. We consider the Laplacian. which is given in spherical coordinates by (2)
1. Poisson-Boltzmann 1.1. Poisson equation. We consider the Laplacian operator (1) 2 = 2 x + 2 2 y + 2 2 z 2 which is given in spherical coordinates by (2) 2 = 1 ( r 2 ) + 1 r 2 r r r 2 sin θ θ and in
More informationTheoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 11
WiSe 22..23 Prof. Dr. A-S. Smith Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg Problem. Theoretische Physik 2:
More informationRANDOM PROPERTIES BENOIT PAUSADER
RANDOM PROPERTIES BENOIT PAUSADER. Quasilinear problems In general, one consider the following trichotomy for nonlinear PDEs: A semilinear problem is a problem where the highest-order terms appears linearly
More informationBochner s Theorem on the Fourier Transform on R
Bochner s heorem on the Fourier ransform on Yitao Lei October 203 Introduction ypically, the Fourier transformation sends suitable functions on to functions on. his can be defined on the space L ) + L
More information4: The Pandemic process
4: The Pandemic process David Aldous July 12, 2012 (repeat of previous slide) Background meeting model with rates ν. Model: Pandemic Initially one agent is infected. Whenever an infected agent meets another
More informationThe Gaussians Distribution
CSE 206A: Lattice Algorithms and Applications Winter 2016 The Gaussians Distribution Instructor: Daniele Micciancio UCSD CSE 1 The real fourier transform Gaussian distributions and harmonic analysis play
More informationMultiple time step Monte Carlo simulations: Application to charged systems with Ewald summation
JOURNAL OF CHEMICAL PHYSICS VOLUME 11, NUMBER 1 1 JULY 004 Multiple time step Monte Carlo simulations: Application to charged systems with Ewald summation Katarzyna Bernacki a) Department of Chemistry
More informationSimulations of bulk phases. Periodic boundaries. Cubic boxes
Simulations of bulk phases ChE210D Today's lecture: considerations for setting up and running simulations of bulk, isotropic phases (e.g., liquids and gases) Periodic boundaries Cubic boxes In simulations
More informationAppendix lecture 9: Extra terms.
Appendi lecture 9: Etra terms The Hyperolic method has een used to prove that d n = log + 2γ + O /2 n This can e used within the Convolution Method to prove Theorem There eists a constant C such that n
More informationMultiresolution analysis by infinitely differentiable compactly supported functions. N. Dyn A. Ron. September 1992 ABSTRACT
Multiresolution analysis by infinitely differentiable compactly supported functions N. Dyn A. Ron School of of Mathematical Sciences Tel-Aviv University Tel-Aviv, Israel Computer Sciences Department University
More informationAn Inverse Mass Expansion for Entanglement Entropy. Free Massive Scalar Field Theory
in Free Massive Scalar Field Theory NCSR Demokritos National Technical University of Athens based on arxiv:1711.02618 [hep-th] in collaboration with Dimitris Katsinis March 28 2018 Entanglement and Entanglement
More informationSolutions to Tutorial 11 (Week 12)
THE UIVERSITY OF SYDEY SCHOOL OF MATHEMATICS AD STATISTICS Solutions to Tutorial 11 (Week 12) MATH3969: Measure Theory and Fourier Analysis (Advanced) Semester 2, 2017 Web Page: http://sydney.edu.au/science/maths/u/ug/sm/math3969/
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 informationFOURIER INVERSION. an additive character in each of its arguments. The Fourier transform of f is
FOURIER INVERSION 1. The Fourier Transform and the Inverse Fourier Transform Consider functions f, g : R n C, and consider the bilinear, symmetric function ψ : R n R n C, ψ(, ) = ep(2πi ), an additive
More informationMolecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia
Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties
More informationMolecular Dynamics. A very brief introduction
Molecular Dynamics A very brief introduction Sander Pronk Dept. of Theoretical Physics KTH Royal Institute of Technology & Science For Life Laboratory Stockholm, Sweden Why computer simulations? Two primary
More informationSome Background Material
Chapter 1 Some Background Material In the first chapter, we present a quick review of elementary - but important - material as a way of dipping our toes in the water. This chapter also introduces important
More informationSYSM 6303: Quantitative Introduction to Risk and Uncertainty in Business Lecture 4: Fitting Data to Distributions
SYSM 6303: Quantitative Introduction to Risk and Uncertainty in Business Lecture 4: Fitting Data to Distributions M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu
More informationMolecular Dynamics Simulations
Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline Long-range Interactions Ewald Sum Fast Multipole Method Spherically Truncated Coulombic Potential Speeding up Calculations SPaSM
More informationOn some weighted fractional porous media equations
On some weighted fractional porous media equations Gabriele Grillo Politecnico di Milano September 16 th, 2015 Anacapri Joint works with M. Muratori and F. Punzo Gabriele Grillo Weighted Fractional PME
More informationPoisson summation and the discrete Fourier transform John Kerl University of Arizona, Math 523B April 14, 2006
Poisson summation and the discrete Fourier transform John Kerl University of Arizona, Math 523B April 4, 2006 In Euclidean space, given a vector,... 2 ... we can put down a coordinate frame (say an orthonormal
More informationAbel Inversion using the Maximum Entropy Method
Abel Inversion using the Maximum Entropy Method Danilo R. Neuber Wolfgang von der Linden 3rd October 2003 Inst. für Theoretische Physik, Tel.: +43/3 16/8 73-81 88; neuber@itp.tu-graz.ac.at Inst. für Theoretische
More informationChern forms and the Fredholm determinant
CHAPTER 10 Chern forms and the Fredholm determinant Lecture 10: 20 October, 2005 I showed in the lecture before last that the topological group G = G (Y ;E) for any compact manifold of positive dimension,
More informationFourier transform. Stefano Ferrari. Università degli Studi di Milano Methods for Image Processing. academic year
Fourier transform Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Methods for Image Processing academic year 27 28 Function transforms Sometimes, operating on a class of functions
More informationand the seeds of quantisation
noncommutative spectral geometry algebra doubling and the seeds of quantisation m.s.,.,, stabile, vitiello, PRD 84 (2011) 045026 mairi sakellariadou king s college london noncommutative spectral geometry
More informationMultiscale Modelling, taking into account Collisions
Multiscale Modelling, taking into account Collisions Andreas Adelmann (Paul Scherrer Institut) October 28, 2017 Multiscale Modelling, taking into account Collisions October 28, 2017 Page 1 / 22 1 Motivation
More informationFast Multipole Methods: Fundamentals & Applications. Ramani Duraiswami Nail A. Gumerov
Fast Multipole Methods: Fundamentals & Applications Ramani Duraiswami Nail A. Gumerov Week 1. Introduction. What are multipole methods and what is this course about. Problems from physics, mathematics,
More informationPhysics 342 Lecture 22. The Hydrogen Atom. Lecture 22. Physics 342 Quantum Mechanics I
Physics 342 Lecture 22 The Hydrogen Atom Lecture 22 Physics 342 Quantum Mechanics I Friday, March 28th, 2008 We now begin our discussion of the Hydrogen atom. Operationally, this is just another choice
More informationGeneral Instructions
Math 240: Real Analysis Qualifying Exam May 28, 2015 Name: Student ID#: Problems/Page Numbers Total Points Your Score Problem 1 / Page 2-3 40 Points Problem 2 / Page 4 Problem 3 / Page 5 Problem 4 / Page
More informationSubdiffusion in a nonconvex polygon
Subdiffusion in a nonconvex polygon Kim Ngan Le and William McLean The University of New South Wales Bishnu Lamichhane University of Newcastle Monash Workshop on Numerical PDEs, February 2016 Outline Time-fractional
More informationA DaFT (Mixed Radix) FFT. for DL_POLY_4
A DaFT (Mixed Radix) FFT Table of Contents for DL_POLY_4 Ian J. Bush Nag Ltd. Jordan Hill Oxford OX2 8DR 1Introduction and Background...2 2The Ewald Summation and SPME...3 3DaFT...4 4DaFT (Mixed Radix)...7
More informationTopics in Harmonic Analysis Lecture 1: The Fourier transform
Topics in Harmonic Analysis Lecture 1: The Fourier transform Po-Lam Yung The Chinese University of Hong Kong Outline Fourier series on T: L 2 theory Convolutions The Dirichlet and Fejer kernels Pointwise
More informationTopic 4 Notes Jeremy Orloff
Topic 4 Notes Jeremy Orloff 4 auchy s integral formula 4. Introduction auchy s theorem is a big theorem which we will use almost daily from here on out. Right away it will reveal a number of interesting
More informationMGA Tutorial, September 08, 2004 Construction of Wavelets. Jan-Olov Strömberg
MGA Tutorial, September 08, 2004 Construction of Wavelets Jan-Olov Strömberg Department of Mathematics Royal Institute of Technology (KTH) Stockholm, Sweden Department of Numerical Analysis and Computer
More informationEECS 117 Lecture 18: Magnetic Energy and Inductance
University of California, Berkeley EECS 117 Lecture 18 p. 1/2 EECS 117 Lecture 18: Magnetic Energy and Inductance Prof. Niknejad University of California, Berkeley University of California, Berkeley EECS
More informationFunction spaces and multiplier operators
University of Wollongong Research Online Faculty of Informatics - Papers (Archive Faculty of Engineering and Information Sciences 2001 Function spaces and multiplier operators R. Nillsen University of
More information