Seminar in Particles Simulations. Particle Mesh Ewald. Theodora Konstantinidou
|
|
- Isabel Dorthy Patterson
- 5 years ago
- Views:
Transcription
1 Seminar in Particles Simulations Particle Mesh Ewald Theodora Konstantinidou
2 This paper is part of the course work: seminar 2 Contents 1 Introduction Role of Molecular Dynamics - MD Molecular potential 3 3 Modelling the physical system Limitations Necessity of the calculation of the electrostatic potential Periodic boundary conditions The Ewald method for the calculation of the potential experienced by an ion Reciprocal space Reciprocal vectors Fourier Analysis Application to Ewald method Particle Mesh Ewald method Application of the method Results of the parallelization of the method in CBM d Orleans 19 8 Advantages of the PME method 20
3 1 Introduction 3 Molecular dynamics (MD) according to the definition given by [1] is a computer simulation technique where the time evolution of a set of interacting atoms is simulated by integrating their equations of motion. In molecular dynamics we follow the laws of classical mechanics, and most notably Newton s law: F i = m i a i, where m i is the atom mass, a i its acceleration, and F i the force acting upon it, due to the interactions with other atoms. 1.1 Role of Molecular Dynamics - MD Molecular dynamics could be applied to a vast area of scientific research. The following are certain areas where MD simulations have been applied with good results. Scientists are optimistic for further and greater contribution of MD in the future, especially with regards to the amazing development of hardware. [1] Liquids Defects Fracture Surfaces Friction Clusters Biomolecules Electronic properties and dynamics 2 Molecular potential In MD simulations one is interested in simulating the trajectories of interacting atoms, thus obtaining a good overview of the phenomena taking place. Most thermodynamic processes or chemical reactions are associated with specific motions of molecules. In most cases the discretization of Newton s laws gives a good approximation of the real motion. For the calculation of trajectories it is important to compute the potential (energy in a sense) of each atom as it influences considerably each atom s motion. A typical empirical force field formulation
4 of the potential energy for N atoms at positions r 1,...,r N, can be written as: 4 a i E = bonds 2 (l i l i0 ) n n angles 4ε i j [ ( σ i j b i 2 (θ i θ i0 ) 2 + ) 12 ( σ i j ) 6 + i=1 j 1 r i j r i j }{{} Lennard Jonespotential torsions V n (1 + cos(nω γ)) (1) 2 ] + 1 n n q i q j 2 i=1 j 1 r i j }{{} Coulombpotential Figure 1 illustrates the different potentials between bonds in atoms. Figure 1: Potentials appearing in a molecule [5] We have to mention here, that in MD simulations the input is not taken from real experiments, but it is the computer which generates it. In the equation 1 for example the values of the parameters a i,b i,v n should be derived from experiments whenever a different molecule is to be simulated. The characteristics of bonds among atoms however, are in most cases the same no matter to which molecule they belong. Introducing empirical values consequently, does not entail the danger of simulating a different system. Molecular dynamics is not the only approach to examine such systems. Quantum mechanics is a more rigorous one. The drawback in quantum mechanics is that it is the most computationally expensive method. In practice scientists simulate small molecules and they make the assumption that the same behavior holds for larger ones. It is a fact however, that such an assumption holds for sure for well known types of molecules, whereas there is no guarantee
5 for the rest. The molecule in figure 2 illustrates the possible movements of bonds. In most biomolecular simulations we are interested in the behavior of very big molecules as 5 Figure 2: representation of a molecule and the acting potentials [3] shown in figure 3, comprised of thousands of atoms.
6 6 Figure 3: The bovine pancreatic trypsin inhibitor protein[4] 3 Modelling the physical system As mentioned above, the physical system is completely determined by its potential, constructed from the relative positions of the atoms with respect to each other, rather than from the absolute positions. Forces are obtained as the gradients of the potential with respect to atomic displacements: F i = ri E(r i,...,r N ) Energies can be calculated using either molecular mechanics or quantum mechanics methods. It is clear that Newton s laws cannot describe situations including drastic changes in electronic structure such as bond making/breaking (see section Limitations ). For such systems a quantum mechanics non-equilibrium simulation is used. Knowledge of the atomic forces and masses can then be used to solve for the positions of each atom along a series of extremely small time steps (on the order of femtoseconds = seconds). The resulting series of snapshots of structural changes over time is called a trajectory. The use of this method to compute trajectories can be more easily seen when Newton s equation is expressed in the form, m i d 2 r i dt 2 = F i = de dr i. (2) Equation 2 gives the accelerations for each atom. In practice (for the systems we simulate) we do not know the analytical solution and hence we can only we compute numerically the velocities and trajectories from equations 3 and 4 respectively
7 7 dv i dt = a i, (3) dr i dt = v i. (4) The leapfrog method is a common numerical approach to calculating trajectories based on Newton s equation. Characteristic feature in this method is the use of two different grids, one for the velocity and one for the position. The steps can be summarized as follows: de dr i = F i = m i a i (t), v i (t + t/2) = v i (t t/2) + a i (t) t, r i (t + t) = r i (t) + v i (t + t/2) t. 3.1 Limitations The de Broglie thermal wavelength is used to define the limits between MD and Quantum Mechanics simulations [1]. 2πh Λ = 2 Mk B T Λ: Broglie thermal wavelength Λ << α. where M is the atomic mass and T the temperature h Plank s constant, k B Boltzman s constant and a the mean nearest neighbor separation. When the fraction Λ α is at maximum 0.1 one can apply Newtons laws. However, when the distance between atoms is very small (a) this fraction becomes larger which means that the system is in a non equilibrium state. As a result an MD simulation is not recommended. The above test refers to the energy of the system. One should in addition control the stiffness of the system. When simulating a stiff system we confront the problem of simulating the rapidly vanishing modes as well as the normal phenomena which finally prevail. One should then focus on the prevailing frequencies and not end the simulation before the relaxation of the quantities we are interested in. The size of the system can cause inaccuracies, in terms of space correlation functions. It can happen that the size of the simulation cell will be such, that the correlations will be totally ignored or misinterpreted. This problem can be partially reduced by a method
8 known as finite size scaling. This consist of computing a physical property A using several boxes with different sizes L, and then fitting the results on a relation 8 A(L) = A 0 + c L n L:length of each side of simulation cell. 4 Necessity of the calculation of the electrostatic potential The accurate calculation of macromolecular structures and dynamics remains a challenge [13]. Basically the inaccuracies are due to the simplifications we do on the conformation of the simulated system. In addition to that problem molecular dynamics simulations using explicit solvent molecules (that means that the solvent, most often water, is treated as discrete molecules whereas implicit solvents are treated as continuum) often exhibit unrealistic behavior on the time scale they take place. Recent work has shown that simulations using current force fields without truncation of Coulombic interactions do not exhibit this unrealistic behavior. In the article of [16] it is mentioned that molecular dynamics simulations involving explicit solvent molecules have usually been performed under one of the following boundary conditions on the Coulombic interactions: Periodic boundary conditions using a finite cutoff, similar to the approach of Verlet for Lennard Jones system (At this approach the only interactions taken into consideration are within the area of a circle whose radius Rcutoff is defined empiricaly) Periodic boundary conditions using Ewald summation Periodic boundary conditions together with a reaction field Nonperiodic boundary conditions together with some treatment of the system environment interface.
9 4.1 Periodic boundary conditions 9 The boundary conditions of the simulated system cause intrinsic difficulties. If for example the system simply terminates at the boundary then the surface effects occur, which causes problems when simulating a region from the interior of a system. An solution to this problem is to use periodic boundary conditions (PBC). When using PBC, the particles are enclosed in a simulation cell which is replicated to all three dimensions until infinity. That means that the whole space is comprised of copies of our simulation cell. This excludes any interaction with the surroundings. In other words, if a particle is located at position r in the box, we assume that this particle really represents an infinite set of particles located at: r + la + mb + nc, l,m,n (, ) integers, a,b,c : vectors along cell axis. Advantages: Avoid surface effects all physical properties remain invariant of simulation cell translation It is important to understand that this periodicity is in space not in time. Fig. visualization of this periodicity. 4 gives a Figure 4: Periodic boundary conditions. The central box is outlined by a thicker line. [6]
10 10 5 The Ewald method for the calculation of the potential experienced by an ion The electrostatic potential caused by charged bodies (i.e. atoms, ions) is given according to Coulobs law by formula 5 (we make the assumption that interactions are always between two charges). The problem with this formula is that it is only conditionally convergent. That means that if one would change the order of summation the series might not converge. E(r 1,r 2,...,r N ) = 1 2 n i j q i q j r i r j + n (5) For solving this problem, Ewald proposed a conversion (based on Adel s theorem) of the non converged series using Gaussian distributions as convergence factors. Practically, he added and then subtracted charge in the form of Gaussian distribution and thus obtained a convergent series. Then he split the resulting integral (summation) into two; One entailing only Gaussian distributions and the second one entailing both Gaussian distributions and delta functions of the charge. The first summation is taking place in the reciprocal space (see 5.1) where the Poisson s equation 2 φ a = 4πρ can very easily be solved. In the following we give the fundamental concepts of Fourier transformation. It is obvious that calculating the second derivative of exponentials provides a significant gain of computation time. The second summation is solved in normal space since the the Fourier expansion of a periodic array of delta functions do not decay for large reciprocal vectors. 5.1 Reciprocal space For a better understanding of Ewald s method we give here an overview of reciprocal space [7]. Consider the non-orthogonal lattice illustrated in figure 5. If a function f is periodic on the lattice, it means that f (r) and f at all of the points Q share the same value. The function f (r) depends only on r s location relative to the cell. The function f thus needs to throw away all the parts of r that have to do with translation, and just focus on the remaining bit of r. Thus we seek to write r in terms of our basis vectors a and b. It should be clear that r can always be written as a linear combination of a and b (i.e.,thata and b span the space), but how do we find the particular linear combination with r = αa + βb? The job of finding this linear combination is made easier by constructing the vectors a and b. a is constructed so that it is perpendicular to b and is scaled so that when dotted with
11 11 Figure 5: one of the benefits of reciprocal space is that lattice does not need to be orthogonal as shown in the picture[7] a it yields 1. The fact that, in dimensioned variables, a has the units of reciprocal length, contributes to the name reciprocal space.) Now if we dot both side of the above equation with a, the term drops out (since it is perpendicular to b) and since a a =1 we have: α = ra β = rb Now f needs to ignore the whole number parts of α andβ, as the fraction parts state where r is in the base rhombus. One way of ignoring the whole number parts is via periodic functions like the trigonometric functions sine and cosine. For example, cos[2π(nα + mβ)] = cos[(n2πa + m2πb )r] = cos[(na + mb )r] = cos[gr] 5.2 Reciprocal vectors The basic reciprocal vectors can easily be derived from the normal ones by the following formulas: a = b c ab c b = c a ab c c = a b ab c ba = 2πδ ab ca = 2πδ ca bc = 2πδ bc G = υ 1 a + υ 2 b + υ 3 c G reciprocal vector
12 12 n(r + T) = n G exp(igr)exp(igt) G since exp(igt) = 1, exp(igt) = exp[i(υ 1 a + υ 2 b + υ 3 c )(u 1 a + u 2 b + u 3 c)] n(r + T) = n(r) = exp[i2π(υ 1 u 1 + υ 2 u 2 + υ 3 u 3 )] (u i,υ j integers) Periodic functions on the real lattice can be represented on the reciprocal space by their frequencies. E.g a wave of type ψ(r) = e ikr shown in figure 6 in real space and in reciprocal or k-space Figure 6: illustration of the simplification the k-space offers [10] 5.3 Fourier Analysis In Fourier analysis the fundamental principle is that every periodic function can be represented as a linear combination of the cosine and sine functions [9]. When a function is periodic with period T the following is true: In one dimension the Fourier expansion is n(r + T) = n(r)fort R d. n(x + α) = n 0 + [C p cos(2πpx/α) + S p sin(2πpx/α)] = n 0 + [C p cos(2πpx/α + 2πp) + S p sin(2πpx/α + 2πp)] p>0 = n 0 + [C p cos(2πpx/α) + S p sin(2πpx/α)] = n(x). p>0
13 We can rewrite the above with the use of exponential function as n(x) = n p exp(i2πp/α). p 13 In three dimension we use the concept of reciprocal vectors G ana the Fourier expansion becomes n(r) = n G exp(igr) G: reciprocal vector of the 3D lattice. G The inverse Fourier Transformation is n G = Vc 1 dv n(r) exp( igr). cell 5.4 Application to Ewald method As already mentioned, Ewald just split the electrostatic potential [9] caused by the charge of atoms (or ions), into two potentials the real and the reciprocal in order to preserve convergence. The following computations reffer to the potential of a single atom at a certain position. φ = φ 1 + φ 2 φ 1 is the reciprocal space potential and is computed as the difference of the sum of the potential caused by the charge of all atoms φ a minus the potential of the charge of the atom itself φ b (see fig. 7) φ 1 = φ a φ b Figure 7 displays the splitting. In the figure, lattice 1 is the real space potential, and lattice 2 the reciprocal. Figure 7: Ewald summation [8]
14 The figure illustrates the splitting that Ewald used to achieve convergence. The parameter of the chosen distribution will not affect the result since their sum is zero. A good choice of parameters can enhance the speed of convergence. We express first the potential φ a and the charge density which causes it as a sum of Fourier series: φ a = c G exp(igr) G ρ = ρ G exp(igr). G Poisson s equation, 6, relates the potential to the charge distribution = G 14 2 φ a = 4πρ (6) G 2 c G exp(igr) = 4π ρ G exp(igr). G Details about the derivation of φ 1 (i) can be found in [9], we obtain φ 1 (i) = 4π S(G)G 2 exp( G 2 /4η) 2q i (η/π) 1/2 G The index i signifies a position in the lattice. The potential of the real space is the sum of three potentials: the point potential, the potential of the Gaussian distribution of the charge density that is in the sphere of radius r l over the l-th atom position and the potential of the Gaussian distribution of the charge density that is out the sphere of radius r l. The last term is a correction term known as background charge density. q l [ 1 r l 1 r l rl 0 ρ(r)dr r l Summing this over all atoms of the unit cell yields φ 2 = l q l r l F(η 1/2 r l ), f ρ(r) ] r dr where F(x) = (2/π 1/2 ) exp( s 2 )ds. x Finally we obtain for the potential of the atom at position i φ(i) = 4π S(G)G 2 exp( G 2 /4η) 2q i (η/π) 1/2 G }{{} φ 1 + q l F(η 1/2 r l ). l r l }{{} φ 2
15 6 Particle Mesh Ewald method 15 The Particle Mesh Ewald method is based on the Partical Mesh algorithm. The basic idea is, to solve Poisson s equation in PBCs taking advantage of fast Fourier transform (FFT) for calculating discrete Fourier transforms. The potential caused by Coulombic interactions is treated as one of a short range describing the influence of the neighboring atoms and the long range potential (in Fourier space ) as a situation in the whole domain. The four principal steps for the particle mesh calculation are: 1. Assign charge to the mesh 2. Solve the field equation on the mesh 3. Calculate the mesh defined force field 4. Interpolate to find forces on the particles In Particle Mesh Ewald method the steps are similar, the difference is on how the charge density is approximated. Here, the complex exponentials (because of PBC) of the charge are just locally interpollated on the grid points whereas in Particle Mesh method the Gaussian distributions are approximated. 1. Assign charge to the mesh using a weighting function 2. Grid transformation - Fourier space representation of the discrete charge distribution 3. Calculate the Potential in Fourier space 4. Transformation to real space 5. Interpolate the potential to the atom and differentiate for calculation of force (Smooth PME method) Figure 8 illustrates the method applied for the computation of the real space term of the potential. As one can observe only the atoms in the area defined by the radius Rcutoff are assumed to interact with the atom whose potential we wish to compute. Figure 9 illustrates graphically the above mentioned steps.
16 16 Figure 8: The real space term of the potential is calculated with R cuto f f methods [15] Figure 9: The steps a,b,c,d illustrate the time sequence of the method.[15] 6.1 Application of the method Suppose there are N point charges q 1,q 2,...,q N at positions r 1,r 2,...r N within the unit cell U satisfying q 1 + q q N = 0. The vectors a α,α = 1,2,3 which need not be orthogonal define the edges of the unit cell. The conjugate reciprocal vectors a α are defined by the relation a α a β = δ αβ. The point charge q i at position r i has fractional coordinates s αi,α = 1,2,3 defined by s αi = a αr i. The charges interact according to Coulob s law with periodic boundary conditions. Thus a point charge q i at position r i interacts with other charges q j at positions r j as well as with all of their periodic images at positions r j +n 1 α 1 +n 2 α 2 +n 3 α 3 for all integers n 1,n 2,n 3. It also interacts with its own periodic images at r i + n 1 α 1 + n 2 α 2 + n 3 α 3 non zero. The
17 electrostatic energy of the unit cell can be written : E(r 1,r 2,...,r N ) = 1 2 n i j q i q j r i r j + n We define the reciprocal lattice vectors m by m = m 1 a 1 + m 2 a 2 + m 3 a 3 and the structure factor S(m) of the unit cell given as 17 N N S(m) = q j exp(2πm r j ) = q j exp[2πi(m 1 s 1 j + m 2 s 2 j + m 3 s 3 j ]. (7) j=1 j=1 Equation 9 gives the part of the energy computed in normal (or direct) space, equation 10 is the part of the energy computed in reciprocal space (the result will be transformed back with the inverse Fourier transformation) and the energy computed by equation 11 represents a correction term. The correction term is subtracted from the total because it represents the interactions between atoms of a molecule that are very close to each other. Since this energy was added to the first two parts it should then be subtracted. E dir = 1 2 n E rec = 1 2 m N i, j er f c(β r j r i + n ) r j r i + n exp( π 2 m 2 /β 2 ) m 2 S(m)S( m) (9) 0E cor = 1 er f c(β r j r i + n ) 2 β N r (i, j) M j r i + n π i=1q 2 i (10) In order to compute the force acting on a atom at position i one has to add the above given energies and then differetiate with respect to r i. Since we have periodic boundary conditions (similar to crystal stucture), the electrostatic structure factors are approximated by interpolation of complex exponentials. As shown in figure 9, we start with a random distribution of charges, which we then interpolate at the grid points. On this gridded charge distribution we calculate the reciprocal energy. This energy is differentiated in order to get the forces. In order to compute the force acting on an atom at a position r i one has interpolate back to the atoms position. Given K 1,K 2,K 3 are positive integers and a point r in the unit cell we define the scaled fractional coordinates by u α = K α a α r for α = 1,2,3 and 0 u α < K α ( ) ( 2πim1 u 1 2πim2 u 2 exp(2πimr) = exp exp K 1 K 2 ) exp ( ) 2πim3 u 3 Let [u] denote integer part of u then we use the following equation for the interpolation of the exponentials ( ) ( ) ( ) 2πimα u α 2πimα 2πimα exp (1 (u α [u α ]))exp [u α ] +(u α [u α ])exp ([u α ]+1) K α K α K α K 3 (8) (11)
18 W 2 (u) = 1 u f or 1 u < 1 W 2 (u) = 0 f or u > 1 ( ) ( ) 2πimα u α 2πimα exp K α W 2 (u a k) exp k k= K α 18 For higher order interpolation the terms W i, i = 1,2,3 are read as W 2p (u) = j=p 1 j= p, j p (u + j k) j=p 1 j= p, j p ( j k) k u k + 1 k = p, p + 1,..., p 1 The index p defines the order of the interpolation. The linear case is for p = 1. S(m) S(m) = = K 1 1 k 1 =0 K 2 1 k 2 =0 K 3 1 k 3 =0 ( ) ( ) 2πimα u α 2πimα exp K α W 2p (u a k) exp k k= K α N q i i=1 k 1 = k 2 = k 3 = ( ) ( ) 2πim1 2πim2 exp k 1 exp k 2 K 1 K 2 [ ( m1 k 1 Q(k 1,k 2,k 3 )exp 2πi where F(Q)Fourier Transformation of Q Q(k 1,k 2,k 3 ) = N i=1 W 2p (u 1i k 1 ) W 2p (u 2i k 2 ) W 2p (u 3i k 3 ) ( ) 2πim3 exp k 3 K 3 K 1 + m 2k 2 K 2 + m 3k 1 K 3 )] = F(Q)(m 1,m 2,m 3 ) n 1,n 2,n 3 q i W 2p (u 1i k 1 n 1 K 1 ) W 2p (u 2i k 2 n 2 K 2 ) W 2p (u 3i k 3 n 3 K 3 ) E rec Ẽ rec = 1 2πV exp( π 2 m 2 /β 2 ) m 0 m 2 F(Q)(m 1,m 2,m 3 )F(Q)( m 1, m 2, m 3 ) (12) The symbol in equation 12 is justified by the interpolation. The resulting reciprocal energy will be transformed back to cartesian space.
19 19 7 Results of the parallelization of the method in CBM d Orleans The following results come from the university of Orlean (Centre Biophysique moleculaire) [15]. The simulation is of water molecules whose behavior is considered like the one of atoms. One can observe that the more the processors the less accurate the method. The ideal case for this method is for small simulations where the number of processors does not go over 8. #5000 molecules of water (15000 atoms) #Length of simulation box =10 #Number of grid points for each grid direction=25 #Interpolation order =4 #Rcutoff = 1.25 P=4 P:number of processors Energy of direct space = Energy of reciprocal space= Energy of correction term= P=8 P:number of processors Energy of direct space = Energy of reciprocal space= Energy of correction term= P=16 P:number of processors Energy of direct space = Energy of reciprocal space= Energy of correction term= P=32 P:number of processors Energy of direct space = Energy of reciprocal space= Energy of correction term= P=64 P:number of processors Energy of direct space = Energy of reciprocal space= Energy of correction term= time= time= time=time= time= time= time= time= time= time= time= time= time= time= time= time=
20 8 Advantages of the PME method 20 High accuracy: The method can provide high accuracy results ( relative force) with a relatively small augmentation to computtation time [12]. Ease of implementation: The PME method can be efficiently implemented into conventional MD Continuity: Both potentials computed with the PME method are continuous functions of positions which alleviates the problem of integrating discontinuous functions. Efficiency: The PME method is fast. For large macromolecular systems, the PME requires only about a 40% overhead over conventional truncated list-based methods to obtain relative force accuracies of References [1] Furio Ercolessi A molecular dynamics primer [2] Celeste Sagui and Thomas A.Darden Molecular Dynamics simulations of biomolecules: Long Range Electrostatic Effects, Annu.Rev. Biomol.Struct. 28: (1999) [3] molecular_mechanics_document.html [4] [5] [6] [7] T.Kirkman [8] [9] Charles Kittel, Einfuerung in Festkoerperphysic (1991) [10] official/m.mehring/condmat/reciprocal.htm [11] R. Hockney, and J. Eastwood Computer simulation using particles McGraw-Hill, New York. (1981) [12] T.Darden, D.York and L Pederson Particle mesh Ewald: An Nlog(N) for Ewald sums in large systems J.Chem.Phys.98 (1993) page 10089
21 [13] Ulrich Essman, Lalith Perera, and Max L.Berkowitz- Department of Chemistry, University of North Carolina A smoothe particle mesh Ewald method J.Chem.Phys. (1995) page 8577 [14] Brock A.Luty, Ilario G.Tironi and Wilfred F. van Gunsteren Lattice-sum methods for calculating electrostatic interactions in molecular simulations J.Chem.Phys.103(1995) page 3014 [15] /ihperf98/ihperf.htm. [16] M. Allen and D. Tildesley Computer Simulation of Liquids Oxford Science, London, (1990) 21
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 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 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 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 informationTutorial on the smooth particle-mesh Ewald algorithm
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
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 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 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 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 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 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 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 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 informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationUsing Molecular Dynamics to Compute Properties CHEM 430
Using Molecular Dynamics to Compute Properties CHEM 43 Heat Capacity and Energy Fluctuations Running an MD Simulation Equilibration Phase Before data-collection and results can be analyzed the system
More informationChapter 1. Crystal structure. 1.1 Crystal lattices
Chapter 1 Crystal structure 1.1 Crystal lattices We will concentrate as stated in the introduction, on perfect crystals, i.e. on arrays of atoms, where a given arrangement is repeated forming a periodic
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 informationCastep & Real Space Integrals
Castep & Real Space Integrals MJR March 3, 213 1 Real and Reciprocal Representations The wavefunction in Castep is expressed in a plane wave basis. That the system is periodic means that the permitted
More informationMolecular dynamics simulation. CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror
Molecular dynamics simulation CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror 1 Outline Molecular dynamics (MD): The basic idea Equations of motion Key properties of MD simulations Sample applications
More informationfree space (vacuum) permittivity [ F/m]
Electrostatic Fields Electrostatic fields are static (time-invariant) electric fields produced by static (stationary) charge distributions. The mathematical definition of the electrostatic field is derived
More information4/18/2011. Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania
1. Introduction Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania Bibliography Computer experiments Ensemble averages and time averages Molecular
More informationarxiv: v1 [cond-mat.stat-mech] 6 Jan 2014
arxiv:1401.1181v1 [cond-mat.stat-mech] 6 Jan 2014 Determination of Forces from a Potential in Molecular Dynamics (note) Bernard Monasse Mines-ParisTech, Cemef bernard.monasse@mines-paristech.fr January
More informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
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 informationClassical Electrostatics for Biomolecular Simulations
pubs.acs.org/cr Classical Electrostatics for Biomolecular Simulations G. Andreś Cisneros, Mikko Karttunen, Pengyu Ren, and Celeste Sagui*, Department of Chemistry, Wayne State University, Detroit, Michigan
More informationStructural Bioinformatics (C3210) Molecular Mechanics
Structural Bioinformatics (C3210) Molecular Mechanics How to Calculate Energies Calculation of molecular energies is of key importance in protein folding, molecular modelling etc. There are two main computational
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 informationIntroduction to first-principles modelling and CASTEP
to first-principles modelling and Phil Hasnip to + Atomistic Simulations If we know what the bonding in a material is beforehand, then we can often find good expressions for the forces between atoms, e.g.
More informationCrystals, X-rays and Proteins
Crystals, X-rays and Proteins Comprehensive Protein Crystallography Dennis Sherwood MA (Hons), MPhil, PhD Jon Cooper BA (Hons), PhD OXFORD UNIVERSITY PRESS Contents List of symbols xiv PART I FUNDAMENTALS
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 informationLecture 11: Potential Energy Functions
Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions
More informationCollective Effects. Equilibrium and Nonequilibrium Physics
Collective Effects in Equilibrium and Nonequilibrium Physics: Lecture 3, 3 March 2006 Collective Effects in Equilibrium and Nonequilibrium Physics Website: http://cncs.bnu.edu.cn/mccross/course/ Caltech
More informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
More informationGromacs Workshop Spring CSC
Gromacs Workshop Spring 2007 @ CSC Erik Lindahl Center for Biomembrane Research Stockholm University, Sweden David van der Spoel Dept. Cell & Molecular Biology Uppsala University, Sweden Berk Hess Max-Planck-Institut
More informationChapter 21. Electric Fields
Chapter 21 Electric Fields The Origin of Electricity The electrical nature of matter is inherent in the atoms of all substances. An atom consists of a small relatively massive nucleus that contains particles
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
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 informationMetropolis, 2D Ising model
Metropolis, 2D Ising model You can get visual understanding from the java applets available, like: http://physics.ucsc.edu/~peter/ising/ising.html Average value of spin is magnetization. Abs of this as
More informationMOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE
MOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE Tatsuto Kimura* and Shigeo Maruyama** *Department of Mechanical Engineering, The University of Tokyo 7-- Hongo,
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationA Nobel Prize for Molecular Dynamics and QM/MM What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential
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 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 information1. Hydrogen atom in a box
1. Hydrogen atom in a box Recall H atom problem, V(r) = -1/r e r exact answer solved by expanding in Gaussian basis set, had to solve secular matrix involving matrix elements of basis functions place atom
More informationGibbs Ensemble Computer Simulations of Vapor Liquid Equilibrium of Hydrazine
CHME 498 Undergraduate Research Final Report New Mexico State University Gibbs Ensemble Computer Simulations of Vapor Liquid Equilibrium of Hydrazine Summer 2017 Gustav Barraza Faculty Adviser: Dr. Martha
More informationA MOLECULAR DYNAMICS SIMULATION OF A BUBBLE NUCLEATION ON SOLID SURFACE
A MOLECULAR DYNAMICS SIMULATION OF A BUBBLE NUCLEATION ON SOLID SURFACE Shigeo Maruyama and Tatsuto Kimura Department of Mechanical Engineering The University of Tokyo 7-- Hongo, Bunkyo-ku, Tokyo -866,
More informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationPractical in Numerical Astronomy, SS 2012 LECTURE 9
Practical in Numerical Astronomy, SS 01 Elliptic partial differential equations. Poisson solvers. LECTURE 9 1. Gravity force and the equations of hydrodynamics. Poisson equation versus Poisson integral.
More informationSection 10 Metals: Electron Dynamics and Fermi Surfaces
Electron dynamics Section 10 Metals: Electron Dynamics and Fermi Surfaces The next important subject we address is electron dynamics in metals. Our consideration will be based on a semiclassical model.
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
More informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
More informationReciprocal Lattice. Let A(r) be some function featuring discrete translation symmetry:
Reciprocal Lattice From Quantum Mechanics we know that symmetries have kinematic implications. Each symmetry predetermines its specific quantum numbers and leads to certain constraints (conservation laws/selection
More informationCoupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation
Coupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation Bo Li Math Dept & CTBP, UCSD Li-Tien Cheng (Math, UCSD) Zhongming Wang (Math & Biochem, UCSD) Yang Xie (MAE,
More informationMolecular Dynamics Simulation of a Nanoconfined Water Film
Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.
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 informationPhonons and lattice dynamics
Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed
More informationSpectral Processing. Misha Kazhdan
Spectral Processing Misha Kazhdan [Taubin, 1995] A Signal Processing Approach to Fair Surface Design [Desbrun, et al., 1999] Implicit Fairing of Arbitrary Meshes [Vallet and Levy, 2008] Spectral Geometry
More informationd 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.
4 Legendre Functions In order to investigate the solutions of Legendre s differential equation d ( µ ) dθ ] ] + l(l + ) m dµ dµ µ Θ = 0. (4.) consider first the case of m = 0 where there is no azimuthal
More informationSolutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and
More informationLecture 11: Long-wavelength expansion in the Neel state Energetic terms
Lecture 11: Long-wavelength expansion in the Neel state Energetic terms In the last class we derived the low energy effective Hamiltonian for a Mott insulator. This derivation is an example of the kind
More informationIntroduction to Molecular Dynamics
Introduction to Molecular Dynamics Dr. Kasra Momeni www.knanosys.com Overview of the MD Classical Dynamics Outline Basics and Terminology Pairwise interacting objects Interatomic potentials (short-range
More informationTennessee s State Mathematics Standards Precalculus
Tennessee s State Mathematics Standards Precalculus Domain Cluster Standard Number Expressions (N-NE) Represent, interpret, compare, and simplify number expressions 1. Use the laws of exponents and logarithms
More informationHarbor Creek School District
Unit 1 Days 1-9 Evaluate one-sided two-sided limits, given the graph of a function. Limits, Evaluate limits using tables calculators. Continuity Evaluate limits using direct substitution. Differentiability
More informationDIFFERENTIATION RULES
3 DIFFERENTIATION RULES DIFFERENTIATION RULES Before starting this section, you might need to review the trigonometric functions. DIFFERENTIATION RULES In particular, it is important to remember that,
More information5 Topological defects and textures in ordered media
5 Topological defects and textures in ordered media In this chapter we consider how to classify topological defects and textures in ordered media. We give here only a very short account of the method following
More informationSeparation of Variables in Polar and Spherical Coordinates
Separation of Variables in Polar and Spherical Coordinates Polar Coordinates Suppose we are given the potential on the inside surface of an infinitely long cylindrical cavity, and we want to find the potential
More information6. X-ray Crystallography and Fourier Series
6. X-ray Crystallography and Fourier Series Most of the information that we have on protein structure comes from x-ray crystallography. The basic steps in finding a protein structure using this method
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture Molecular Dynamics Simulation David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu MD of hard disks intuitive Review and Preview collision
More informationCan a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water?
Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? Ruhong Zhou 1 and Bruce J. Berne 2 1 IBM Thomas J. Watson Research Center; and 2 Department of Chemistry,
More informationHigh-Resolution. Transmission. Electron Microscopy
Part 4 High-Resolution Transmission Electron Microscopy 186 Significance high-resolution transmission electron microscopy (HRTEM): resolve object details smaller than 1nm (10 9 m) image the interior of
More informationPhysics 9 Spring 2012 Midterm 1 Solutions
Physics 9 Spring 22 NAME: TA: Physics 9 Spring 22 Midterm s For the midterm, you may use one sheet of notes with whatever you want to put on it, front and back. Please sit every other seat, and please
More informationClass 29: Reciprocal Space 3: Ewald sphere, Simple Cubic, FCC and BCC in Reciprocal Space
Class 29: Reciprocal Space 3: Ewald sphere, Simple Cubic, FCC and BCC in Reciprocal Space We have seen that diffraction occurs when, in reciprocal space, Let us now plot this information. Let us designate
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 informationContinuum Limit and Fourier Series
Chapter 6 Continuum Limit and Fourier Series Continuous is in the eye of the beholder Most systems that we think of as continuous are actually made up of discrete pieces In this chapter, we show that a
More informationMolecular Dynamics. What to choose in an integrator The Verlet algorithm Boundary Conditions in Space and time Reading Assignment: F&S Chapter 4
Molecular Dynamics What to choose in an integrator The Verlet algorithm Boundary Conditions in Space and time Reading Assignment: F&S Chapter 4 MSE485/PHY466/CSE485 1 The Molecular Dynamics (MD) method
More information1. (3) Write Gauss Law in differential form. Explain the physical meaning.
Electrodynamics I Midterm Exam - Part A - Closed Book KSU 204/0/23 Name Electro Dynamic Instructions: Use SI units. Where appropriate, define all variables or symbols you use, in words. Try to tell about
More informationCS 273 Prof. Serafim Batzoglou Prof. Jean-Claude Latombe Spring Lecture 12 : Energy maintenance (1) Lecturer: Prof. J.C.
CS 273 Prof. Serafim Batzoglou Prof. Jean-Claude Latombe Spring 2006 Lecture 12 : Energy maintenance (1) Lecturer: Prof. J.C. Latombe Scribe: Neda Nategh How do you update the energy function during the
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 informationdisordered, ordered and coherent with the substrate, and ordered but incoherent with the substrate.
5. Nomenclature of overlayer structures Thus far, we have been discussing an ideal surface, which is in effect the structure of the topmost substrate layer. The surface (selvedge) layers of the solid however
More informationHandout 7 Reciprocal Space
Handout 7 Reciprocal Space Useful concepts for the analysis of diffraction data http://homepages.utoledo.edu/clind/ Concepts versus reality Reflection from lattice planes is just a concept that helps us
More informationHyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea
The Scaled Effective Solvent Method for Predicting the Equilibrium Ensemble of Structures with Analysis of Thermodynamic Properties of Amorphous Polyethylene Glycol-Water Mixtures Hyeyoung Shin a, Tod
More informationDielectrics. Lecture 20: Electromagnetic Theory. Professor D. K. Ghosh, Physics Department, I.I.T., Bombay
What are dielectrics? Dielectrics Lecture 20: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay So far we have been discussing electrostatics in either vacuum or in a conductor.
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 informationKeywords: 2015, IJARCSSE All Rights Reserved Page 315
Volume 5, Issue 3, March 2015 ISSN: 2277 128X International Journal of Advanced Research in Computer Science and Software Engineering Research Paper Available online at: www.ijarcsse.com Special Issue
More informationWest Windsor-Plainsboro Regional School District Algebra and Trigonometry Grades 11-12
West Windsor-Plainsboro Regional School District Algebra and Trigonometry Grades 11-12 Unit 1: Linear Relationships & Functions Content Area: Mathematics Course & Grade Level: Algebra & Trigonometry, 11
More informationBiomolecules are dynamic no single structure is a perfect model
Molecular Dynamics Simulations of Biomolecules References: A. R. Leach Molecular Modeling Principles and Applications Prentice Hall, 2001. M. P. Allen and D. J. Tildesley "Computer Simulation of Liquids",
More information2.2 The derivative as a Function
2.2 The derivative as a Function Recall: The derivative of a function f at a fixed number a: f a f a+h f(a) = lim h 0 h Definition (Derivative of f) For any number x, the derivative of f is f x f x+h f(x)
More informationPURE MATHEMATICS Unit 1
PURE MATHEMATICS Unit 1 FOR CAPE EXAMINATIONS DIPCHAND BAHALL CAPE is a registered trade mark of the Caribbean Examinations Council (CXC). Pure Mathematics for CAPE Examinations Unit 1 is an independent
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 information2. Acoustic Wave Equation
2. Acoustic Wave Equation 2.1. Acoustic Wave Equation Acoustic wave equation is a second order partial differential equation. The acoustic wave equation for one dimensional plane wave in Cartesian coordinates
More informationProblem Set Number 01, MIT (Winter-Spring 2018)
Problem Set Number 01, 18.377 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Thursday, March 8, 2018. Turn it in (by 3PM) at the Math.
More information1 + lim. n n+1. f(x) = x + 1, x 1. and we check that f is increasing, instead. Using the quotient rule, we easily find that. 1 (x + 1) 1 x (x + 1) 2 =
Chapter 5 Sequences and series 5. Sequences Definition 5. (Sequence). A sequence is a function which is defined on the set N of natural numbers. Since such a function is uniquely determined by its values
More informationSolvation and Macromolecular Structure. The structure and dynamics of biological macromolecules are strongly influenced by water:
Overview Solvation and Macromolecular Structure The structure and dynamics of biological macromolecules are strongly influenced by water: Electrostatic effects: charges are screened by water molecules
More informationIntroduction)! Electrostatics is the study of stationary electric charges and fields (as opposed to moving charges and currents)
Higher'Physics'1B Electricity) Electrostatics)) Introduction) Electrostatics is the study of stationary electric charges and fields (as opposed to moving charges and currents) Properties)of)Electric)Charges)
More informationStructuralBiology. October 18, 2018
StructuralBiology October 18, 2018 1 Lecture 18: Computational Structural Biology CBIO (CSCI) 4835/6835: Introduction to Computational Biology 1.1 Overview and Objectives Even though we re officially moving
More informationCourse 2BA1: Hilary Term 2006 Section 7: Trigonometric and Exponential Functions
Course 2BA1: Hilary Term 2006 Section 7: Trigonometric and Exponential Functions David R. Wilkins Copyright c David R. Wilkins 2005 Contents 7 Trigonometric and Exponential Functions 1 7.1 Basic Trigonometric
More informationAdvantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations
The Hilltop Review Volume 7 Issue 1 Winter 2014 Article 10 December 2014 Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations Tai-Hsien Wu Western Michigan University
More informationA Fast N-Body Solver for the Poisson(-Boltzmann) Equation
A Fast N-Body Solver for the Poisson(-Boltzmann) Equation Robert D. Skeel Departments of Computer Science (and Mathematics) Purdue University http://bionum.cs.purdue.edu/2008december.pdf 1 Thesis Poisson(-Boltzmann)
More informationMultiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015
Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction
More information