Excluded Volume Effects in Gene Stretching. Pui-Man Lam Physics Department, Southern University Baton Rouge, Louisiana
|
|
- Preston Norton
- 5 years ago
- Views:
Transcription
1 Excluded Volume Effects in Gene Stretcing Pui-Man Lam Pysics Department, Soutern University Baton Rouge, Louisiana 7083 Abstract We investigate te effects excluded volume on te stretcing of a single DNA in solution. We find tat for small force F, te extension is not linear in F but proportion to F γ, wit γ(-ν)/ν, were ν is te well-known universal correlation lengt exponent. A freely joint cain model wit te segment lengt cosen to reproduce te small extension beavior gives excellent fit to te experimental data of λ-page DNA over te wole experimental range. We sow tat excluded volume effects are stronger in two dimensions and also derive results in two dimensions wic are different from te tree dimensional results. Tis suggests experiments to be performed in tese lower dimensions. I. INTRODUCTION Many forms of polymers or macromolecules suc as proteins, lipids, fatty acids and DNA, are vital for life. Te DNA is a very large molecule. For instance te DNA in one uman cell can encode approximately 00,000 genes. It is partitioned into 46 cromosomes, 3 from eac parent and is normally twisted and folded into te nucleus of a cell a few micrometers in widt. In solution te uman DNA consists of two strands of polymeried nucleotides, twisted into a rigt and elix. New tecnologies makes it possible to study single molecules of DNA. Suc studies can give insigts for
2 development of laboratory tecniques for analying, fractionating and sequencing DNA. One suc study investigates ow a single DNA stretces wen it is pulled wit a force at its two ends, using micromanipulation tools suc as optical tweeers or microneedles to grasp te molecule by its two ends. Stretcing te DNA can provide a muc clearer picture under te microscope. As te average end-to-end separation increases, a force arises tat opposes te stretcing. Tis force is entropic in nature since te number of configurations of te polymer decreases as increases, leading to a decrease in te entropy. Suc a force versus extension curve F( ) for a force applied in te -direction ad been directly measured for a 97-kb λ-page DNA dimer []. Te experimental data can be fitted remarkably well using a worm-like cain (WLC) polymer model []. However, in te WLC model, excluded volume effect was not taken into account. For a polymer in solution, excluded volume effect is always present, except peraps at te tricritical θ-point. Terefore tere is in general no reason for tis effect to be negligible for DNA in solution. In tis paper we investigate te effect of excluded volume in DNA stretcing. We find tat for a week force F applied at te ends of te molecule, te extension in te direction of te force is not linear in F, but rater proportional to F γ, wit γ(-ν)/ν, were ν is te well-known universal correlation lengt exponent. Using te result ν0.6, [3] tis gives γ/3. We find tat te experimental data for λ-page DNA at small F is actually in agreement wit tis beavior wit γ Including te excluded volume effect in a freely joint cain (FJC) model by coosing te segment lengt to reproduce te small extension beavior gives excellent fit to te full range of experimental data.
3 In two dimensions, te exponent ν is exactly known to be ¾[3]. Tis will give an extension proportional to F /3 in a week force F. In te FJC model, te functional form of te force versus extension curve F( ) is different from tat of te tree dimensional case. Again te excluded volume effects can be included by coosing te segment lengt in te FJC to reproduce te small extension beavior. Tis suggests experiments to be performed in two-dimensional geometry to test te stronger effect of excluded volume in lower dimensions. In section II we will review te models in wic excluded volume effect is neglected, including te Gaussian cain and te FJC model. In section III we will discuss DNA model wit te excluded volume effect and use it to fit te experimental data. In section IV we will discuss te two dimensional case. Section V is te conclusion. II. MODELS WITHOUT EXCLUDED VOLUME EFFECT A. GAUSSIAN CHAIN In te Gaussian cain model, a long DNA molecule at equilibrium in solution as random walk statistics. Te probability distribution p() of its end-to-end distance is approximately Gaussian [4]: p ) π 3 / 3 3 exp, < R > < R > ( () were <R> is te average end-to-end distance. For fixed, te entropy of te DNA is Sk B log[p()], were k B is te Boltmann s constant. For a fixed end-to-end distance, te force in te -direction is given by 3
4 F S T 3 k / < R > () For a polymer of wit a contour lengt L, te persistence lengt P can be defined in te Gaussian cain as P<R> /(L) [4]. Tis persistence lengt is independent of te contour lengt since for a Gaussian cain, <R> L. In terms of te persistence lengt P, equation () can be written as LP 3k T B F. (3) Te extension is linear in te force, so te DNA beaves as a Hookean spring wit ero natural lengt and a temperature dependent effective spring constant. B. FREELY JOINT CHAIN In te freely joint model [5], te cain of contour lengt L consists of (L/b) freely joint segments eac of lengt b. Wen tis cain is stretced at its two ends wit a force F, eac segment is pulled at bot ends wit tis force F and tus tends to align it in te direction of F, wile termal fluctuations tend to orient it in random directions. For a force in te direction and denoting te angle between F and a particular segment as θ, te partition function of te cain can be written as Z π bf cos( θ ) exp sin( ) d k π θ θ 0 L / b 4πk bf sin bf k L / b (4) Te average extension of te cain in te direction is given by k B T log Z F bf k Lcot k bf bf LΛ, k (5) 4
5 5 were Λ(x)cot(x)-x - is te Langevin function. For small x, Λ(x) x/3. Terefore for small F, LbF/(3k B T). Comparison wit equation (3) gives bp. Substituting into equation (5) ten gives II. MODELS WITH EXCLUDED VOLUME EFFECT For a polymer wit excluded volume effect, te probability distribution of te end-to- end distance is given by [6] were A is a normaliation constant. Te exponent is related to te universal exponent for te average end-to-end distance by /(-ν) [7]. Te exponent α is strictly ero at te upper critical dimension d c 4, wen ν/ and te probability distribution reduces to te form given in equation (). We will sow in te following tat α sould be ero in all dimensions. Te entropy is given by Te force exerted at te ends of te cain is given by Λ T k FP L B ) (6, 3 exp ) ( > < > < α R R A p (7) > < > < + α R R A k p T k S B B 3 log log ) ( log (8) > < > < 3 α R R T k S T F B ) (9
6 Unless α0, te first term would imply te unpysical result tat te force F becomes infinite as te extension goes to ero. Putting α0 in equation (7), te component of te force in te direction is given by F S T 3 k ) / ( ) ( < R > (0) Te average force F can be obtained by taking x y : F < R > < R > 3 / k ( ) / 3 k (3 ) () Te extension of te cain in te direction due to a force F in te direction is given by F < R > / 3 k () Terefore wit excluded volume effect, te extension is not linear in F for. In Figure we sow te double logaritmic plot of te experimental data for te extension versus force of a 97-kb λ-page DNA dimer []. Te slope of te curve at small extension is about Tis is in good agreement wit te exponent /(-)/3 obtained using te known value of ν0.6 [3]. Tis sows tat excluded volume effect is important for DNA. In te freely joint cain model we again assume a cain consisting of L/b freely joint segments eac of lengt b. Ten te average extension in te direction of te force is 6
7 again given by equation (5). Te effects of excluded volume can be included in te FJC model by coosing te segment lengt b to reproduce te small extension beavior. Comparing te small extension result of tat equation LbF/(3k B T) wit equation (), we obtain for te segment lengt b 3k < > F R / FL 3 k (3) Substituting tis into equation (5) we obtain te expression for te extension in te direction of te force for a cain wit excluded volume effect F < R > LΛ / L 3 k 3 (4) In contrast to te case of a Gaussian cain, te persistence lengt defined by P<R> /(L) is not a constant but would depend on te contour lengt L. Terefore we use ere te average end-to-end distance <R> as a parameter, rater tan te persistence lengt. We ave fitted te experimental data []to equation (4), wit /(-ν).5, but varying te parameters L and <R> to minimie te least square deviation. Te best fit is obtained wit <R>3.9µm and L3.8µm. In Figure we sow te result of our fit wit te experimental data of te extension versus force. Te dased lines are te results of te freely joint cain wit no excluded volume effect. Te value P3nm corresponds to best least square fit wit te data and P50nm corresponds to a fit to te small extension data. We see tat te fit to te experimental data is excellent for te model wit excluded volume effect. A least square fit varying all tree parameters <R>, L and yields.6. Tis sows tat te form /(-ν) is in fact very good. 7
8 We ave also fitted te more recent data of Wang et al [8 ] to equation (4). Teir experimental data were obtained using a double stranded DNA of muc sorter sie (~µm). Wang et al. found tat teir data could be fitted well wit te teoretical forms of Marko and Siggia [9], (wit persistence lengt P43.5nm, contour lengt L.3µm) and Odjik [0] (wit persistence lengt P4.nm, contour lengt L.3µm). Bot forms take into account also elasticity teory. However, tey found tat teir experimental data were not well fit by te freely joint cain expression over any force region. Tis finding in fact supports conclusions of previous work wit double-stranded DNA molecules. Te freely joint cain model, owever can be successfully applied to describe singly stranded DNA. Since our teory wit excluded volume is based on te freely joint cain model we do not expect it to fit well te data of a double stranded DNA. In spite of tis we sow in Figure 3 our best fit to te data of Wang et al. [8]. Te dased curve is obtained using te freely joint cain model witout excluded volume, wit persistence lengt P4nm and contour lengt L.3µm. Te solid curve is te best fit using our expression (4) wit <R>0.394µm, L.3µm and.5. Even toug bot te freely joint cain and our equation (4) do not fit te data well we find tat our expression seem to fit better over most of te data range. We also find tat te largest slope in te data is about 0.5. Terefore even toug te force-extension relation is nonlinear, te exponent is muc smaller tan te value /3 predicted by our teory. Tis is again due to te double stranded nature of te DNA. It sould be pointed out tat our expression (4) applies only wen te cain is not yet close to fully extended. Wen te cain is close to fully extended, te freely joint cain expression (6) sould be better, because ten te excluded volume effect must be 8
9 negligible. One could, in principle tink of a way of bridging te two expressions (6) and (4). However, it is not known at wat extension, or force, is te excluded volume effect in te cain negligible. Only wen te cain is fully extended it is clear tat te excluded volume effect is negligible. But te fully extended configuration is only one among a large number of configurations wose number increases exponentially wit te lengt of te cain. For a cain tat is close to but not yet fully extended te critical force for te crossover can only be determined by comparing te fit of te experimental data to (6) and (4) to see at wat extension is te expression (6) a better fit tan (4). Looking at Figure, one can see tat te fits to te data from (6) and (4) are almost indistinguisable. Terefore we prefer to leave (4) as it is, keeping tis condition in mind. III. DNA IN TWO DIMENSIONS A. WITHOUT EXCLUDED VOLUME EFFECT Te partition function for a freely joint cain given in equation (4) is valid only for tree dimensions. In two dimensions it as to be replaced by Z L / b π Fbcosθ Fb dθ exp I 0 k π k 0 L / b (5) were I 0 is te modified Bessel function of order ero. Te average extension in te direction of te force is given by k Fb Fb T log Z L I I (6) F k k B 0 9
10 were I is te modified Bessel function of order one. From te beaviors of I 0 and I for small arguments one obtains LFb, F 0 k T B (7) Comparing tis to te Gaussian cain result for te two dimensional case, PLF /(k B T), we find bp. Substituting tis into equation (6), we find Fb L I k 0 Fb I k (8) for te extension of a two-dimensional DNA in direction of te force F. B. WITH EXCLUDED VOLUME EFFECT To take into account excluded volume effect we use equation (0) wit x + and find F k < R > ( ) k < R > (9) From tis we find ( ) / F < R > k (0) In a freely joint cain model wit (L/b) segments, eac of lengt b, te extension in te direction of te force is given by (6), wit te small F beavior given by equation (7). Again te excluded volume effect can be included by coosing te segment lengt in te FJC model to reproduce te small extension beavior. From equations (7) and 0
11 (0) we find b ( ) / k T F < R > B () LF k T B Substituting tis value of b into equation (6) we ave ( ) / F < R > LI L k ( ) / F < R > I0 L k () In two dimensions, te exponent for te end-to-end distance is known exactly to be ν3/4 [3]. Tis gives /(-ν)4, a muc larger deviation from te value in te case of no excluded volume effect. Only two parameters L and <R> in equation () are left to fit experimental data. However since no experimental data is known for DNA in two dimensions, we will use te same values of L and <R> obtained in tree dimensions. Te results are sown in Figure 4 for bot wit and witout excluded volume effect. We suggest experiments to be performed in two dimensions to test te stronger effects of excluded volume in lower dimensions. IV. CONCLUSION We investigate te effects excluded volume on te stretcing of a single DNA in solution. We find tat for small force F, te extension is not linear in F but proportion to F γ, wit γ(-ν)/ν, were ν is te well-known universal correlation lengt exponent. A freely joint cain model wit te segment lengt cosen to reproduce te small extension beavior gives excellent fit to te experimental data of λ-page DNA over te wole
12 experimental range. Te parameters used in te fitting are te contour lengt L and te average end-to-end distance <R>. Te best fit is obtained wit <R>3.9µm and L3.8µm. Te fit to te experimental data is just as good as te worm like cain model [], in wic only te elastic energy is taken into account by treating te cain as a uniform elastic rod. It is possible tat bot effects are present in DNA. We sow tat excluded volume effects are stronger in two dimensions and also derive results in two dimensions wic are different from te tree dimensional results. Tis suggests experiments to be performed in tese lower dimensions Acknowledgement: Tis work as been supported by te Department of Energy grant DE-FG0-97ER5343 REFERENCES [] S. B. Smit, L. Fini, C. Bustemante, Science 58, (99) [] C. Bustamante, J.F. Marko, E.D. Siggia and S. Smit, Science 65, 599 (994) [3] H.E. Stanley, Introduction to Pase Transition and Critical Penomena, Oxford University Press (997) [4] R. H. Austin, J.P. Brody, E.C. Cox, T. Duke and W. Volkmut, Pysics Today, February, 3 (997) [5] A. Yu. Grosberg and A.R. Koklov, Statistical Pysics of Macromolecules (AIP Press, New York, 994); J.M. Scurr, S.B. Smit, Biopolymers 9, 6 (990) [6] C. Domb, J. Gillis and G. Wilmer, Proc. Pys. Soc. (London) 85, 65 (965) [7] M.E. Fiser, J. Cem. Pys. 44, 66 (966) [8] M. Wang et al. Biopys. J. 7, 335 (997) [9] J.F. Marko and E.D. Siggia, Macromolecules 8, 8759 (995) [0] T. Odjik, Macromolecules 8, 706 (995)
13 FIGURE CAPTIONS Figure : Double logaritmic plot of extension versus force for λ-page DNA. Te slope of te curve at small extension is about 0.64 Figure : Extension versus force. Te solid dots are experimental data of λ-page DNA. Te solid line is for freely-joint cain wit excluded volume effects, using te best least square fit values for <R> and L and fixed value of.5. Te dased curves are results of freely joint cain, for two values of te persistence lengt. Figure 3: Extension versus force. Te solid dots are experimental data of ref. [8]. Te solid line is for free-joint cain wit excluded volume effects, using te best least square fit values of <R> and L and fixed value of.5. Te dased curve is te result of te freely joint cain. Figure 4: Extension versus force for a two dimensional polymer in te freely joint cain model. Te dased curve is for te freely joint cain model in two dimensions witout excluded volume, wit persistence lengt P50 nm. Te dotted curve is te same ting for tree dimensions, for comparison. Te solid curve is for te freely joint cain model in two dimensions wit excluded volume, using parameter values of <R> and L obtained in tree dimensions. 3
14 EXTENSION( µ m) 0 slope FORCE(pN) Figure 4
15 35 30 P50(nm) EXTENSION( µ m) P3(nm) <R>3.9µm L3.8µm FORCE(pN) Figure 5
16 .5 extension( µm) P4(nm) slope0.5 <R>0.394(µm) L.3(µm) force(pn) Figure 3 6
17 35 30 Gaussian EXTENSION( µ m) dimensions Gaussian 3 dimensions excluded volume 5 dimensions FORCE(pN) Figure 4 7
A = 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 informationWhy gravity is not an entropic force
Wy gravity is not an entropic force San Gao Unit for History and Pilosopy of Science & Centre for Time, SOPHI, University of Sydney Email: sgao7319@uni.sydney.edu.au Te remarkable connections between gravity
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 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 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 informationQuantum Theory of the Atomic Nucleus
G. Gamow, ZP, 51, 204 1928 Quantum Teory of te tomic Nucleus G. Gamow (Received 1928) It as often been suggested tat non Coulomb attractive forces play a very important role inside atomic nuclei. We can
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
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 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 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 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 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 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 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 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 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 informationChapter 2 Ising Model for Ferromagnetism
Capter Ising Model for Ferromagnetism Abstract Tis capter presents te Ising model for ferromagnetism, wic is a standard simple model of a pase transition. Using te approximation of mean-field teory, te
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 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 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 informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More informationProblem Set 4: Whither, thou turbid wave SOLUTIONS
PH 253 / LeClair Spring 2013 Problem Set 4: Witer, tou turbid wave SOLUTIONS Question zero is probably were te name of te problem set came from: Witer, tou turbid wave? It is from a Longfellow poem, Te
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationPolynomials 3: Powers of x 0 + h
near small binomial Capter 17 Polynomials 3: Powers of + Wile it is easy to compute wit powers of a counting-numerator, it is a lot more difficult to compute wit powers of a decimal-numerator. EXAMPLE
More informationSection 15.6 Directional Derivatives and the Gradient Vector
Section 15.6 Directional Derivatives and te Gradient Vector Finding rates of cange in different directions Recall tat wen we first started considering derivatives of functions of more tan one variable,
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 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 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 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 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 informationCubic Functions: Local Analysis
Cubic function cubing coefficient Capter 13 Cubic Functions: Local Analysis Input-Output Pairs, 378 Normalized Input-Output Rule, 380 Local I-O Rule Near, 382 Local Grap Near, 384 Types of Local Graps
More information. If lim. x 2 x 1. f(x+h) f(x)
Review of Differential Calculus Wen te value of one variable y is uniquely determined by te value of anoter variable x, ten te relationsip between x and y is described by a function f tat assigns a value
More informationModels for Single Molecule Experiments
Models for Single Molecule Experiments Binding Unbinding Folding - Unfolding -> Kramers Bell theory 1 Polymer elasticity: - Hooke s law (?) - FJC (Freely joined chain) - WLC (Worm-like chain) 3 Molecular
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 informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 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
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More informationNotes on wavefunctions II: momentum wavefunctions
Notes on wavefunctions II: momentum wavefunctions and uncertainty Te state of a particle at any time is described by a wavefunction ψ(x). Tese wavefunction must cange wit time, since we know tat particles
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 informationMath 262 Exam 1 - Practice Problems. 1. Find the area between the given curves:
Mat 6 Exam - Practice Problems. Find te area between te given curves: (a) = x + and = x First notice tat tese curves intersect wen x + = x, or wen x x+ =. Tat is, wen (x )(x ) =, or wen x = and x =. Next,
More informationADCP MEASUREMENTS OF VERTICAL FLOW STRUCTURE AND COEFFICIENTS OF FLOAT IN FLOOD FLOWS
ADCP MEASUREMENTS OF VERTICAL FLOW STRUCTURE AND COEFFICIENTS OF FLOAT IN FLOOD FLOWS Yasuo NIHEI (1) and Takeiro SAKAI (2) (1) Department of Civil Engineering, Tokyo University of Science, 2641 Yamazaki,
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
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 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 informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
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 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 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 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 information1 2 x Solution. The function f x is only defined when x 0, so we will assume that x 0 for the remainder of the solution. f x. f x h f x.
Problem. Let f x x. Using te definition of te derivative prove tat f x x Solution. Te function f x is only defined wen x 0, so we will assume tat x 0 for te remainder of te solution. By te definition of
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 informationHT TURBULENT NATURAL CONVECTION IN A DIFFERENTIALLY HEATED VERTICAL CHANNEL. Proceedings of 2008 ASME Summer Heat Transfer Conference HT2008
Proceedings of 2008 ASME Summer Heat Transfer Conference HT2008 August 10-14, 2008, Jacksonville, Florida USA Proceedings of HT2008 2008 ASME Summer Heat Transfer Conference August 10-14, 2008, Jacksonville,
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 informationSin, Cos and All That
Sin, Cos and All Tat James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 9, 2017 Outline Sin, Cos and all tat! A New Power Rule Derivatives
More informationAverage Rate of Change
Te Derivative Tis can be tougt of as an attempt to draw a parallel (pysically and metaporically) between a line and a curve, applying te concept of slope to someting tat isn't actually straigt. Te slope
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 informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More information1 + t5 dt with respect to x. du = 2. dg du = f(u). du dx. dg dx = dg. du du. dg du. dx = 4x3. - page 1 -
Eercise. Find te derivative of g( 3 + t5 dt wit respect to. Solution: Te integrand is f(t + t 5. By FTC, f( + 5. Eercise. Find te derivative of e t2 dt wit respect to. Solution: Te integrand is f(t e t2.
More informationMAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016
MAT244 - Ordinary Di erential Equations - Summer 206 Assignment 2 Due: July 20, 206 Full Name: Student #: Last First Indicate wic Tutorial Section you attend by filling in te appropriate circle: Tut 0
More informationProperties of the Spin-flip Amplitude of Hadron Elastic Scattering and Possible Polarization Effects at RHIC
Properties of te Spin-flip Amplitude of Hadron Elastic Scattering and Possible Polarization Effects at RHIC arxiv:ep-p/0210418v1 30 Oct 2002 O. V. Selyugin 1 Joint Institute for Nuclear Researc, Dubna,
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.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 informationMATH 155A FALL 13 PRACTICE MIDTERM 1 SOLUTIONS. needs to be non-zero, thus x 1. Also 1 +
MATH 55A FALL 3 PRACTICE MIDTERM SOLUTIONS Question Find te domain of te following functions (a) f(x) = x3 5 x +x 6 (b) g(x) = x+ + x+ (c) f(x) = 5 x + x 0 (a) We need x + x 6 = (x + 3)(x ) 0 Hence Dom(f)
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 informationWYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More information1. Questions (a) through (e) refer to the graph of the function f given below. (A) 0 (B) 1 (C) 2 (D) 4 (E) does not exist
Mat 1120 Calculus Test 2. October 18, 2001 Your name Te multiple coice problems count 4 points eac. In te multiple coice section, circle te correct coice (or coices). You must sow your work on te oter
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
More informationQuaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 01 Aug 08, 2016.
Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals Gary D. Simpson gsim1887@aol.com rev 1 Aug 8, 216 Summary Definitions are presented for "quaternion functions" of a quaternion. Polynomial
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 information2.3 Algebraic approach to limits
CHAPTER 2. LIMITS 32 2.3 Algebraic approac to its Now we start to learn ow to find its algebraically. Tis starts wit te simplest possible its, and ten builds tese up to more complicated examples. Fact.
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 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 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 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 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 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 informationSolution for the Homework 4
Solution for te Homework 4 Problem 6.5: In tis section we computed te single-particle translational partition function, tr, by summing over all definite-energy wavefunctions. An alternative approac, owever,
More informationRules of Differentiation
LECTURE 2 Rules of Differentiation At te en of Capter 2, we finally arrive at te following efinition of te erivative of a function f f x + f x x := x 0 oing so only after an extene iscussion as wat te
More informationCalculus I Homework: The Derivative as a Function Page 1
Calculus I Homework: Te Derivative as a Function Page 1 Example (2.9.16) Make a careful sketc of te grap of f(x) = sin x and below it sketc te grap of f (x). Try to guess te formula of f (x) from its grap.
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
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 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 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 information3. Using your answers to the two previous questions, evaluate the Mratio
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0219 2.002 MECHANICS AND MATERIALS II HOMEWORK NO. 4 Distributed: Friday, April 2, 2004 Due: Friday,
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
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 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 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 informationContinuous Stochastic Processes
Continuous Stocastic Processes Te term stocastic is often applied to penomena tat vary in time, wile te word random is reserved for penomena tat vary in space. Apart from tis distinction, te modelling
More informationACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019
ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS MATH00030 SEMESTER 208/209 DR. ANTHONY BROWN 6. Differential Calculus 6.. Differentiation from First Principles. In tis capter, we will introduce
More informationhal , version 1-29 Jan 2008
ENERGY OF INTERACTION BETWEEN SOLID SURFACES AND LIQUIDS Henri GOUIN L. M. M. T. Box 3, University of Aix-Marseille Avenue Escadrille Normandie-Niemen, 3397 Marseille Cedex 0 France E-mail: enri.gouin@univ-cezanne.fr
More informationDerivatives of trigonometric functions
Derivatives of trigonometric functions 2 October 207 Introuction Toay we will ten iscuss te erivates of te si stanar trigonometric functions. Of tese, te most important are sine an cosine; te erivatives
More informationIEOR 165 Lecture 10 Distribution Estimation
IEOR 165 Lecture 10 Distribution Estimation 1 Motivating Problem Consider a situation were we ave iid data x i from some unknown distribution. One problem of interest is estimating te distribution tat
More informationTime (hours) Morphine sulfate (mg)
Mat Xa Fall 2002 Review Notes Limits and Definition of Derivative Important Information: 1 According to te most recent information from te Registrar, te Xa final exam will be eld from 9:15 am to 12:15
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 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 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 informationCHAPTER (A) When x = 2, y = 6, so f( 2) = 6. (B) When y = 4, x can equal 6, 2, or 4.
SECTION 3-1 101 CHAPTER 3 Section 3-1 1. No. A correspondence between two sets is a function only if eactly one element of te second set corresponds to eac element of te first set. 3. Te domain of a function
More informationTechnology-Independent Design of Neurocomputers: The Universal Field Computer 1
Tecnology-Independent Design of Neurocomputers: Te Universal Field Computer 1 Abstract Bruce J. MacLennan Computer Science Department Naval Postgraduate Scool Monterey, CA 9393 We argue tat AI is moving
More information7.1 Using Antiderivatives to find Area
7.1 Using Antiderivatives to find Area Introduction finding te area under te grap of a nonnegative, continuous function f In tis section a formula is obtained for finding te area of te region bounded between
More informationBlueprint Algebra I Test
Blueprint Algebra I Test Spring 2003 2003 by te Commonwealt of Virginia Department of Education, James Monroe Building, 101 N. 14t Street, Ricmond, Virginia, 23219. All rigts reserved. Except as permitted
More information