The Verlet Algorithm for Molecular Dynamics Simulations
|
|
- Morris Townsend
- 5 years ago
- Views:
Transcription
1 Cemistry Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical equations of motion are! F i = m i! a i, i = 1,2,,N. (1)! In tis equation, F i denotes te force on te it particle, m i is te mass of te particle, and a! i is te acceleration. Te arrows over te force and acceleration denote te x, y, and z components of te vectors. In te rest of tis andout, te vector notation will be dropped for convenience. Taylor Series Approximation Te solution of Newton's equations using te Verlet algoritm is based on a Taylor series expansion. Expanding te position of te it particle r i at time t+ and t in Taylor series yields, r i (t + ) = r i (t) + r i ' (t) + 2 2! r i'' (t) + 2 3! r i''' (t) + O( 4 ) (2) r i (t ) = r i (t) r i ' (t) + 2 Adding tese two equations leads to te result 2! r i '' (t) 3 3! r i ''' (t) + O( 4 ). r i (t + ) + r i (t ) = 2 r i (t) + 2 r i '' (t) + O( 4 ). (3) Note tat r i ' (t) is equivalent to te velocity of te particle, v i (t), and r i '' (t) is te acceleration of te particle, a i (t). Making tese substitutions in Eq. (3), we ave r i (t + ) + r i (t ) = 2 r i (t) + 2 a i (t) + O( 4 ). (4) Solving for r i (t + ) and using Newton's equations for te acceleration ( a i = F i /m i ) leads to r i (t + ) = 2 r i (t) r i (t ) + 2 m i F i (t) + O( 4 ). (5) Discretization of Time Te time is a continuous variable, but to implement Eq. (5) on a computer requires te solution to be generated for discrete time steps suc tat t n +1 t n =, (6) were is te step size for te time variable.
2 2 Discretizing time in Eq. (5) leads to te Verlet algoritm (1967), r i (t n +1 ) = 2 r i (t n ) r i (t n 1 ) + 2 m i F i (t n ) + O( 4 ). (7) Note tat te Verlet Algoritm relies on two previous time steps, t n and t n 1, to advance te solution forward in time, and it is a fourt-order metod. Tere are sometimes problems associated wit te Verlet algoritm due to potential loss of accuracy wen Eq. (7) is implemented. Te final term in te equation tends to be small compared to te oters, and tis term may as a result get lost in te round-off error. Notice tat te velocity does not explicitly appear in te Verlet algoritm. Tis may be a problem if properties of te system tat depend on velocity, suc as kinetic energy, are desired. Also, te velocity at time t n cannot be calculated until te position at time t n +1 is obtained. As a result te velocity is usually calculated from te positions using an approximate formula based on te derivative of te position, v i (t n ) = r i (t n +1 ) r i (t n 1 ) 2. (8) Also notice tat te Verlet algoritm is not self-starting, since it requires te positions at two previous times ( t n and t n 1 ) to get te position at te next time, t n+1. Tis is usually overcome by first completing one step of a less accurate metod tat only requires information at one previous time step, suc as te Euler metod, and ten switcing to te Verlet metod for subsequent steps. Te Verlet Leapfrog Algoritm In an attempt to correct some of te problems associated wit te original Verlet algoritm, te Verlet Leapfrog algoritm was developed. Using an approximation for te derivative, te velocity at te midpoint between times t n and t n+1 can be defined as, v i (t n+1/ 2 ) = r i (t n+1 ) r i (t n ), (9) Eq. (8) can be solved for r i (t n +1 ) to yield r i (t n +1 ) = r i (t n ) + v i (t n +1/ 2 ). (10) Similarly, defining te velocity at te midpoint between times t n 1 and t n gives v i (t n 1/ 2 ) = r i (t n ) r i (t n 1 ). (11) Te acceleration at time t n can be defined using anoter approximate formula for te derivative, a i (t n ) = v i (t n+1/ 2 ) v i (t n 1/ 2 ). (12)
3 3 Using Newton's equations for te acceleration, a i (t n ) = F i (t n ) /m i, and equating te acceleration to Eq. (12) yields v i (t n+1/ 2 ) v i (t n 1/ 2 ) = F i (t n ) m i, (13) wic can be solved for v i (t n+1/ 2 ), v i (t n+1/ 2 ) = v i (t n 1/ 2 ) + m i F i (t n ). (14) Eqs. (10) and (14) constitute te equations of te Verlet Leapfrog algoritm. Te Leapfrog algoritm as some advantages over te original Verlet algoritm, Eq. (7). First, te loss of accuracy due to round-off error in te original Verlet algoritm because of vastly different magnitudes of terms in te equation is corrected. Second, te velocities are included explicitly in te metod, unlike te original Verlet algoritm. Tere are still a few problems associated wit te Leapfrog algoritm, owever. First, it is still not self-starting. Tis can be overcome as in te original Verlet metod by completing one step of te Euler metod first, and ten switcing to te Leapfrog algoritm for subsequent steps. Toug te Leapfrog algoritm produces bot positions and velocities, te velocities are obtained at 1/2-integer time steps ( t n +1/ 2, t n + 3/ 2, t n+ 5/ 2,...), wile te positions are obtained at integer time steps ( t n, t n +1, t n + 2,...). If te velocities at integer time steps are required, an auxiliary equation may be utilized wic simply averages te 1/2- integer velocities, v i (t n ) = 1 [ 2 v i (t n 1/ 2 ) + v i (t n+1/ 2 )]. (15) Example: Verlet Leapfrog Algoritm for Bond Stretcing Consider te simple one-dimensional armonic potential as a representation of te motion of te vibrations of a bond in a molecular system, V( x) = 1 2 k x 2. (16) Te coordinate x represents te bond displacement, x = r r eq. Te Verlet leapfrog algoritm for tis system is given by x( t n +1 ) = x( t n ) + v x ( t n +1/ 2 ) v x ( t n +1/ 2 ) = v x ( t n 1/ 2 ) + µ F x( t n ). (17) Te force F x is obtained from te force field, F x = dv dx F x = kx. (18)
4 4 Substituting, te Verlet Leapfrog equations become x( t n +1 ) = x( t n ) + v x ( t n +1/ 2 ) v x ( t n +1/ 2 ) = v x ( t n 1/ 2 ) k m x ( t n). (19) Example: Carbon-Carbon Bond Stretcing Using te Verlet Leapfrog Algoritm Modeling Carbon-Carbon Bonds To apply te Verlet Leapfrog algoritm (Eq. 19) to model te motion of a carbon-carbon bond, ten te mass m would be replaced by te reduced mass µ of te bond. For a bond A-B, te reduced mass µ is defined as µ = m A m B m A + m B. (20) were m A and m B are te masses of atoms A and B. So, for a carbon-carbon bond, te reduced mass is 6.0 g/mol (or about atomic units). Te force constants for carbon-carbon bonds from te on-line CRC (Table 9-75) are listed in Table 1 in units of kcal mol 1 Å 2, N/cm, N/m (or J/m 2 ), and atomic units (artree/bor 2 ). Table 1. Force constants for typical carbon-carbon bonds. C C C Bond Force Constant (kcal mol 1 Å 2 ) Force Constant (N/cm) Force Constant (N/m) Force Constant (a.u.) C (from C 2 H 6 ) C (from CCl 2 CH 2 ) C (from C 2 ) Results: Variation of Force Constants Calculations were performed using te tree force constants given in Table 1 (k=0.289, 0.541, and a.u.), representing carbon-carbon single, double, and triple bonds, respectively. Te Verlet Leapfrog Algoritm was implemented using a value of 5 a.u. for and 500 total time steps wit an initial position of 0 a.u. and an initial velocity of a.u. [Note tat an atomic unit of time is s, or fs. An atomic unit of distance (one bor) equals Å. An atomic unit of velocity equals m/s.] Te position and velocity values are plotted as functions of time to demonstrate te dynamics of a carbon-carbon bond modeled by eac of te force constants in Figures 1-3. Note tat in eac case te position oscillates symmetrically about x=0, a caracteristic of te armonic oscillator. In addition, we see tat te velocity is out of pase wit te position (and it also oscillates symmetrically about v=0). Te figures also sow tat as te force constant increases, te frequency of oscillation of te C-C bond increases, as would be expected as a result of te relationsip between te armonic frequency and force constant. Te period of oscillation terefore decreases as te force constant increases, from ps to ps to ps for force constants of 0.289, 0.541, and a.u., respectively.
5 5 (a) (b) Figure 1. (a) Bond displacement and (b) velocity as functions of time for a model carbon-carbon single bond wit k=0.289 a.u. (a) (b) Figure 2. (a) Bond displacement and (b) velocity as functions of time for a model carbon-carbon double bond wit k=0.541 a.u. (a) (b) Figure 3. (a) Bond displacement and (b) velocity as functions of time for a model carbon-carbon triple bond wit k=0.781 a.u.
6 6 Testing te Accuracy of te Verlet Algoritm Te exact solution for te position at eac time step for te armonic oscillator may be compared to te numerical solution from te Verlet Algoritm. Comparisons were made wit te following time steps and number of points: (a.u.): # steps: Te initial conditions used for eac run were x 0 = 0.0 a.u. and v 0 = a.u. Te absolute error between te Verlet and exact results for te bond displacement coordinate x were plotted as functions of time. Results for eac of te force constants are sown in Figures 4-6. Figure 4. Absolute error in position as a function of time for a model carbon-carbon double bond wit k=0.289 a.u. Figure 5. Absolute error in position as a function of time for a model carbon-carbon double bond wit k=0.541 a.u.
7 7 Figure 6. Absolute error in position as a function of time for a model carbon-carbon double bond wit k=0.789 a.u. Te maximum absolute error in te position sown in eac of Figures 4-6 clearly increases over time. It is also cyclic, since we can see tat te largest error occurs at te maximum or minimum extension of te bond. For eac time step, te maximum absolute error over one oscillation cycle can be compared. Results are given in Table 2. Table 2. Maximum absolute error for different time steps and force constants. Maximum Absolute Error (a.u.) (a.u.) k = a.u. k = a.u. k = a.u Te results sow tat te absolute error decreases as te time step decreases. Similar trends are observed for eac of te tree force constants. Notice, owever, tat te decrease in te error is not linear. A grap sowing te maximum absolute error plotted against te time step is sown in Figure 7.
8 8 Figure 7. Maximum absolute error in position as a function of step size in time for model carboncarbon bonds. Te results sown in Figure 7 appear rougly quadratic in form. Looking at te data, we see tat wen te step size decreases from 10 to 5 a.u. (a factor of 2), te error decreases by a factor of 4. Similarly, wen te step size decreases from 10 to 2 a.u. (a factor of 5), te error decreases by a factor of 25. Te relationsip between step size and error can be explored furter. Assume tat te error as te following form, Error = C M, (21) were C is a constant and M is an integer. Taking te log of bot sides yields, log( Error) = logc + M log( ). (22) Tis equation indicates tat a plot of log(error) on te y-axis and log() on te x-axis sould be linear wit a slope of M. If te error scales quadratically wit te step size, ten te slope sould be equal to 2. Te log plot is sown in Figure 8.
9 9 Figure 8. Log-log plot of maximum error in bond displacement and step size for Verlet Leapfrog algoritm. As illustrated in Figure 8, te slope for all te k values does indeed equal 2. Tis verifies tat wen te step size is decreases by some factor S, ten te error decreases by S 2.
Numerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationPrecalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!
Precalculus Test 2 Practice Questions Page Note: You can expect oter types of questions on te test tan te ones presented ere! Questions Example. Find te vertex of te quadratic f(x) = 4x 2 x. Example 2.
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationThe Derivative as a Function
Section 2.2 Te Derivative as a Function 200 Kiryl Tsiscanka Te Derivative as a Function DEFINITION: Te derivative of a function f at a number a, denoted by f (a), is if tis limit exists. f (a) f(a + )
More informationf a h f a h h lim lim
Te Derivative Te derivative of a function f at a (denoted f a) is f a if tis it exists. An alternative way of defining f a is f a x a fa fa fx fa x a Note tat te tangent line to te grap of f at te point
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationHigher Derivatives. Differentiable Functions
Calculus 1 Lia Vas Higer Derivatives. Differentiable Functions Te second derivative. Te derivative itself can be considered as a function. Te instantaneous rate of cange of tis function is te second derivative.
More informationThe total error in numerical differentiation
AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationDerivatives of Exponentials
mat 0 more on derivatives: day 0 Derivatives of Eponentials Recall tat DEFINITION... An eponential function as te form f () =a, were te base is a real number a > 0. Te domain of an eponential function
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More informationMath 1241 Calculus Test 1
February 4, 2004 Name Te first nine problems count 6 points eac and te final seven count as marked. Tere are 120 points available on tis test. Multiple coice section. Circle te correct coice(s). You do
More information5.1 We will begin this section with the definition of a rational expression. We
Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go
More informationSection 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is
Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationExam 1 Review Solutions
Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More information1. Which one of the following expressions is not equal to all the others? 1 C. 1 D. 25x. 2. Simplify this expression as much as possible.
004 Algebra Pretest answers and scoring Part A. Multiple coice questions. Directions: Circle te letter ( A, B, C, D, or E ) net to te correct answer. points eac, no partial credit. Wic one of te following
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More information(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?
Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationCalculus I Practice Exam 1A
Calculus I Practice Exam A Calculus I Practice Exam A Tis practice exam empasizes conceptual connections and understanding to a greater degree tan te exams tat are usually administered in introductory
More informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More information3.1 Extreme Values of a Function
.1 Etreme Values of a Function Section.1 Notes Page 1 One application of te derivative is finding minimum and maimum values off a grap. In precalculus we were only able to do tis wit quadratics by find
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Matematics National Institute of Tecnology Durgapur Durgapur-7109 email: anita.buie@gmail.com 1 . Capter 8 Numerical Solution of Ordinary
More informationMath 34A Practice Final Solutions Fall 2007
Mat 34A Practice Final Solutions Fall 007 Problem Find te derivatives of te following functions:. f(x) = 3x + e 3x. f(x) = x + x 3. f(x) = (x + a) 4. Is te function 3t 4t t 3 increasing or decreasing wen
More informationDifferential Calculus (The basics) Prepared by Mr. C. Hull
Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit
More informationInvestigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001
Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x
More informationContinuity. Example 1
Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More information(4.2) -Richardson Extrapolation
(.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as
More informationChapters 19 & 20 Heat and the First Law of Thermodynamics
Capters 19 & 20 Heat and te First Law of Termodynamics Te Zerot Law of Termodynamics Te First Law of Termodynamics Termal Processes Te Second Law of Termodynamics Heat Engines and te Carnot Cycle Refrigerators,
More informationLecture 21. Numerical differentiation. f ( x+h) f ( x) h h
Lecture Numerical differentiation Introduction We can analytically calculate te derivative of any elementary function, so tere migt seem to be no motivation for calculating derivatives numerically. However
More informationSection 3.1: Derivatives of Polynomials and Exponential Functions
Section 3.1: Derivatives of Polynomials and Exponential Functions In previous sections we developed te concept of te derivative and derivative function. Te only issue wit our definition owever is tat it
More informationSolution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4.
December 09, 20 Calculus PracticeTest s Name: (4 points) Find te absolute extrema of f(x) = x 3 0 on te interval [0, 4] Te derivative of f(x) is f (x) = 3x 2, wic is zero only at x = 0 Tus we only need
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationLesson 6: The Derivative
Lesson 6: Te Derivative Def. A difference quotient for a function as te form f(x + ) f(x) (x + ) x f(x + x) f(x) (x + x) x f(a + ) f(a) (a + ) a Notice tat a difference quotient always as te form of cange
More informationNUMERICAL DIFFERENTIATION
NUMERICAL IFFERENTIATION FIRST ERIVATIVES Te simplest difference formulas are based on using a straigt line to interpolate te given data; tey use two data pints to estimate te derivative. We assume tat
More informationPractice Problem Solutions: Exam 1
Practice Problem Solutions: Exam 1 1. (a) Algebraic Solution: Te largest term in te numerator is 3x 2, wile te largest term in te denominator is 5x 2 3x 2 + 5. Tus lim x 5x 2 2x 3x 2 x 5x 2 = 3 5 Numerical
More informationTeaching Differentiation: A Rare Case for the Problem of the Slope of the Tangent Line
Teacing Differentiation: A Rare Case for te Problem of te Slope of te Tangent Line arxiv:1805.00343v1 [mat.ho] 29 Apr 2018 Roman Kvasov Department of Matematics University of Puerto Rico at Aguadilla Aguadilla,
More informationExcerpt from "Calculus" 2013 AoPS Inc.
Excerpt from "Calculus" 03 AoPS Inc. Te term related rates refers to two quantities tat are dependent on eac oter and tat are canging over time. We can use te dependent relationsip between te quantities
More informationAN IMPROVED WEIGHTED TOTAL HARMONIC DISTORTION INDEX FOR INDUCTION MOTOR DRIVES
AN IMPROVED WEIGHTED TOTA HARMONIC DISTORTION INDEX FOR INDUCTION MOTOR DRIVES Tomas A. IPO University of Wisconsin, 45 Engineering Drive, Madison WI, USA P: -(608)-6-087, Fax: -(608)-6-5559, lipo@engr.wisc.edu
More informationChapter Seven The Quantum Mechanical Simple Harmonic Oscillator
Capter Seven Te Quantum Mecanical Simple Harmonic Oscillator Introduction Te potential energy function for a classical, simple armonic oscillator is given by ZÐBÑ œ 5B were 5 is te spring constant. Suc
More informationQuantum Numbers and Rules
OpenStax-CNX module: m42614 1 Quantum Numbers and Rules OpenStax College Tis work is produced by OpenStax-CNX and licensed under te Creative Commons Attribution License 3.0 Abstract Dene quantum number.
More information2.3 Product and Quotient Rules
.3. PRODUCT AND QUOTIENT RULES 75.3 Product and Quotient Rules.3.1 Product rule Suppose tat f and g are two di erentiable functions. Ten ( g (x)) 0 = f 0 (x) g (x) + g 0 (x) See.3.5 on page 77 for a proof.
More informationSolutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014
Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More informationAMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.
Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More information1 1. Rationalize the denominator and fully simplify the radical expression 3 3. Solution: = 1 = 3 3 = 2
MTH - Spring 04 Exam Review (Solutions) Exam : February 5t 6:00-7:0 Tis exam review contains questions similar to tose you sould expect to see on Exam. Te questions included in tis review, owever, are
More informationDefinition of the Derivative
Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of
More informationA.P. CALCULUS (AB) Outline Chapter 3 (Derivatives)
A.P. CALCULUS (AB) Outline Capter 3 (Derivatives) NAME Date Previously in Capter 2 we determined te slope of a tangent line to a curve at a point as te limit of te slopes of secant lines using tat point
More informationUniversity Mathematics 2
University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at
More informationTHE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225
THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationLab 6 Derivatives and Mutant Bacteria
Lab 6 Derivatives and Mutant Bacteria Date: September 27, 20 Assignment Due Date: October 4, 20 Goal: In tis lab you will furter explore te concept of a derivative using R. You will use your knowledge
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More information4.2 - Richardson Extrapolation
. - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence
More informationWe name Functions f (x) or g(x) etc.
Section 2 1B: Function Notation Bot of te equations y 2x +1 and y 3x 2 are functions. It is common to ave two or more functions in terms of x in te same problem. If I ask you wat is te value for y if x
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationMath 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0
3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,
More informationWork and Energy. Introduction. Work. PHY energy - J. Hedberg
Work and Energy PHY 207 - energy - J. Hedberg - 2017 1. Introduction 2. Work 3. Kinetic Energy 4. Potential Energy 5. Conservation of Mecanical Energy 6. Ex: Te Loop te Loop 7. Conservative and Non-conservative
More informationFinding and Using Derivative The shortcuts
Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More information3.4 Worksheet: Proof of the Chain Rule NAME
Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are
More informationSome Review Problems for First Midterm Mathematics 1300, Calculus 1
Some Review Problems for First Midterm Matematics 00, Calculus. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd,
More informationChapter 4 Derivatives [ ] = ( ) ( )= + ( ) + + = ()= + ()+ Exercise 4.1. Review of Prerequisite Skills. 1. f. 6. d. 4. b. lim. x x. = lim = c.
Capter Derivatives Review of Prerequisite Skills. f. p p p 7 9 p p p Eercise.. i. ( a ) a ( b) a [ ] b a b ab b a. d. f. 9. c. + + ( ) ( + ) + ( + ) ( + ) ( + ) + + + + ( ) ( + ) + + ( ) ( ) ( + ) + 7
More information1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).
. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use
More informationPre-lab Quiz/PHYS 224 Earth s Magnetic Field. Your name Lab section
Pre-lab Quiz/PHYS 4 Eart s Magnetic Field Your name Lab section 1. Wat do you investigate in tis lab?. For a pair of Helmoltz coils described in tis manual and sown in Figure, r=.15 m, N=13, I =.4 A, wat
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationA Reconsideration of Matter Waves
A Reconsideration of Matter Waves by Roger Ellman Abstract Matter waves were discovered in te early 20t century from teir wavelengt, predicted by DeBroglie, Planck's constant divided by te particle's momentum,
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationGradient Descent etc.
1 Gradient Descent etc EE 13: Networked estimation and control Prof Kan) I DERIVATIVE Consider f : R R x fx) Te derivative is defined as d fx) = lim dx fx + ) fx) Te cain rule states tat if d d f gx) )
More informationExercises for numerical differentiation. Øyvind Ryan
Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can
More informationHOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS
HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS Po-Ceng Cang National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan 36 Tel: 886 3
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare ttp://ocw.mit.edu 5.74 Introductory Quantum Mecanics II Spring 9 For information about citing tese materials or our Terms of Use, visit: ttp://ocw.mit.edu/terms. Andrei Tokmakoff, MIT
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationSection 2.1 The Definition of the Derivative. We are interested in finding the slope of the tangent line at a specific point.
Popper 6: Review of skills: Find tis difference quotient. f ( x ) f ( x) if f ( x) x Answer coices given in audio on te video. Section.1 Te Definition of te Derivative We are interested in finding te slope
More information1watt=1W=1kg m 2 /s 3
Appendix A Matematics Appendix A.1 Units To measure a pysical quantity, you need a standard. Eac pysical quantity as certain units. A unit is just a standard we use to compare, e.g. a ruler. In tis laboratory
More informationlecture 35: Linear Multistep Mehods: Truncation Error
88 lecture 5: Linear Multistep Meods: Truncation Error 5.5 Linear ultistep etods One-step etods construct an approxiate solution x k+ x(t k+ ) using only one previous approxiation, x k. Tis approac enoys
More informationSFU UBC UNBC Uvic Calculus Challenge Examination June 5, 2008, 12:00 15:00
SFU UBC UNBC Uvic Calculus Callenge Eamination June 5, 008, :00 5:00 Host: SIMON FRASER UNIVERSITY First Name: Last Name: Scool: Student signature INSTRUCTIONS Sow all your work Full marks are given only
More informationPreface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.
Preface Here are my online notes for my course tat I teac ere at Lamar University. Despite te fact tat tese are my class notes, tey sould be accessible to anyone wanting to learn or needing a refreser
More informationThe Priestley-Chao Estimator
Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More information2.1 THE DEFINITION OF DERIVATIVE
2.1 Te Derivative Contemporary Calculus 2.1 THE DEFINITION OF DERIVATIVE 1 Te grapical idea of a slope of a tangent line is very useful, but for some uses we need a more algebraic definition of te derivative
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More informationRunge-Kutta methods. With orders of Taylor methods yet without derivatives of f (t, y(t))
Runge-Kutta metods Wit orders of Taylor metods yet witout derivatives of f (t, y(t)) First order Taylor expansion in two variables Teorem: Suppose tat f (t, y) and all its partial derivatives are continuous
More informationLogarithmic functions
Roberto s Notes on Differential Calculus Capter 5: Derivatives of transcendental functions Section Derivatives of Logaritmic functions Wat ou need to know alread: Definition of derivative and all basic
More informationTest 2 Review. 1. Find the determinant of the matrix below using (a) cofactor expansion and (b) row reduction. A = 3 2 =
Test Review Find te determinant of te matrix below using (a cofactor expansion and (b row reduction Answer: (a det + = (b Observe R R R R R R R R R Ten det B = (((det Hence det Use Cramer s rule to solve:
More information