CE 530 Molecular Simulation

Size: px
Start display at page:

Download "CE 530 Molecular Simulation"

Transcription

1 CE 530 Molecular Simulation Lecture 7 Beyond Atoms: Simulating Molecules David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu

2 Review Fundamentals units, properties, statistical mechanics Monte Carlo and molecular dynamics as applied to atomic systems simulating in various ensembles biasing methods for MC Molecular models for realistic (multiatomic) systems inter- and intra- atomic potentials electrostatics Now examine differences between simulations of monatomic and multiatomic molecules

3 Truncating the Potential Many molecular models employ point charges for electrostatic interactions - Potential-truncation schemes must be careful not to split the charges For a 9Å truncation distance, using water-like charges, the interaction energy for a molecule with bare charge is (huge) Always use cutoff based on molecule separation, not atom for large molecules, OK to split molecule but do not split subgroups

4 4 Volume-Scaling Moves Scaling atom displacements leads to large strain on intramolecular bonds Instead perform volume scaling moves using molecule centers-ofmass (or something similar) Let R i be COM of molecule i R i = atoms on i () i irj q j (i) be position of atom j w.r.t. R i () i () i j j For size scaling s, L new = sl old Acceptance based on change in m q = r R () i () i j (new) = s i j r R q i L Number of molecules β ( U PV) N lnv(not atoms) m

5 Rigid vs. Nonrigid Molecules MC and MD can be performed on molecules as already described MD moves advance atom positions based on current forces MC moves translate atoms and accepts based on energy change both are done considering inter- and intra-molecular forces limiting distribution has same form π / N N N β pi mi βu( r ) = dp dr e e N = Number of atoms 3N h 5 if this is all that is done, there is nothing more to say Often it is much more efficient to use a rigid-bond model MD integration then doesn t have to deal with fast intramolecular dynamics, so a larger time step can be used MC can sample configurations more efficiently using rigid-body moves (even if model does not have rigid bonds) but much care is needed to do this properly

6 Molecule Coordinate Frame 6 Molecule-frame coordinates are defined w.r.t. molecule COM with molecule in a reference orientation Simulation-frame coordinate is determined by molecule COM R and orientation ω q j R i For rigid molecules, the molecule-frame coordinates never change

7 Orientation. 7 Orientation described in terms of rotation of molecule frame y Direction cosines can be used to describe rotation Relation between same point in two frames Rotation matrix α = e e α = e e x x x y β = e e β = e e y x y y x α α x = y β β y α α A = β β = r Ar Kinematics of rigid-molecule rotation described in terms of rotation of the molecule coordinate frame (i.e. the direction cosines) r' never changes in a rigid molecule y β y y α x x α y x x x = α x α y x

8 8 Orientation. We also need to invert the relation get the simulation-frame coordinate from the molecule frame value T r = A r = A r Direction cosines are not independent in D, all can be described by just one parameter use rotation angle θ cosθ sinθ A = sinθ cosθ cosθ sinθ = sinθ cosθ inverse can be viewed as replacing θ with -θ A y y θ x x

9 9 x Euler Angles The picture in 3D is similar: (x, y, z) à (x', y', z') Nine direction cosines Three independent coordinates specify orientation Euler angles are the conventional choice ω = φθψ z three rotations give the simulation-frame orientation ζ φ ξ y ζ θ ξ η η z ψ ξ y x η z y x

10 0 3D Rotation Matrix Rotation matrix expressed in terms of Euler angles = r Ar cosψ cosφ cosθsinφsinψ cosψ sinφ cosθ cosφsinψ sinψ sinθ A = sinψ cosφ cosθsinφcosψ sinψ sinφ cosθ cosφcosψ cosψ sinθ sinθsinφ sinθ cosφ cosθ To get space-fixed coordinate, multiply molecule-fixed vector by A - again, A - = A T T r= A r

11 Transforming Coordinates. Consider a simple diatomic positions of two atoms described by x, y, z, x, y, z L z y can instead describe by molecule COM and stretch/orientation coordinates x XYZLθφ,,,,, in molecule frame, each atom position is given by r = Le z r = Le z

12 Transforming Coordinates. To get space-fixed coordinates, use rotation matrix T r = R A r R ( XYZ,, ) T r = Le L r = R A r z r = Le z z y matrix T A cosφ sinφ sinθsinφ = cosθsinφ cosθ cosφ sinθ cosφ 0 sinθ cosθ x the result is x = X Lsinθsinφ x = X Lsinθsinφ y = Y Lsinθcosφ y = Y Lsinθcosφ z = Z Lcosθ z = Z Lcosθ

13 Transforming Coordinates 3. 3 The ensemble distribution for the transformed coordinates is obtained via the Jacobian x y z x y z π = Qh Qh 3N 3N the elements of J are the derivatives X Y Z L φ θ For this transformation But we also need to transform the momenta i N N β p / m βu( r ) dp dr e i i dp dq e e = J N N β p / m βu( q ) J i e N = 8L sinθ N J αβ = r α q β 0 0 sinθ sinφ Lsinθ cosφ Lcosθ sinφ 0 0 sinθ cosφ Lsinθ sinφ Lcosθ cosφ 0 0 cosθ 0 Lsinθ 0 0 sinθ sinφ Lsinθ cosφ Lcosθ sinφ 0 0 sinθ cosφ Lsinθ sinφ Lcosθ sinφ 0 0 cosθ 0 Lsinθ x = X Lsinθ sinφ y = Y Lsinθ cosφ z = Z Lcosθ x = X Lsinθ sinφ y = Y Lsinθ cosφ z = Z Lcosθ

14 Transforming Coordinates 4. 4 Begin with the Lagrangian in the original coordinate system L = K U = (!x!y!z!x!y!z ) U (x, y,z,x, y,z ) =!r i U (r) transform to new coordinates assume m=!x = x X X! x Y Y! x Z Z! x L L! x θ θ! x φ φ!!y = etc. in general!r i =!r i = α β r i q α!q α α =!q g i!q!r i =!q G!q i r i r i!q q α q α!q β β r = x y z x y z r = q = X Y Z L θ φ

15 5 Transforming Coordinates 5. Derive momenta L =!q G!q U (q) p α L!q α = The Hamiltonian is β G aβ!q β p = G!q!q = G p!q G!q = p G p Uses G and thus G - are symmetric H = K( pq, ) U( q) = p G p U( q) The Jacobian for the momentum transformation is the reciprocal of the Jacobian for the coordinate transformation π = dp N Qh 3N q dq N J J e β p q G p q e βu (qn ) we don t have to worry about the Jacobian with the full transform

16 6 Integrating Over Momenta If we integrate out the momentum coordinates, the Jacobian again arises π ( pq, ) = π ( q) = = Qh Qh Qh 3N 3N 3N For the diatomic, this term is N N N β p G p βu( q ) dp dq e N N βu( q ) N dq e N dq e N βu( q ) / e dp e β p G p / G = J =cl sinθ c = constant G In MC simulation, the terms must be included in the construction of the transition-probability matrix

17 7 Averages with Constraints. Very stiff coordinates are sometimes treated as rigidly constrained e.g., the bond length L in the diatomic may be held at a constant value MC and MD have different ways to enforce this constraint Regardless of simulation technique, the constrained-system ensemble average may differ from the unconstrained value even when compared to the limit of an infinitely stiff bond! Why the difference? A rigid constraint implies no kinetic energy in vibration Examine the Lagrangian s = soft coordinate L = r i r i q q α q α q β U (q S ; L) β β Lα L =!qs G S!q S U (q S ; L)

18 8 Averages with Constraints. The Jacobian for the coordinate transform is the same as for the unconstrained average But the momentum Jacobian no longer has the term for the constrained coordinate Thus, in general, the distribution of unconstrained coordinates differs N l N l s qs = 3N l s s Qh The difference is N l βu( qs ) N l s = 3N l s s Qh Qh 3N l βu( q ) N s Gs ps βu qs π ( p, ; L) dp dq e e π ( q ; L) dq e dp e = dq N l s e β p N N s π ( qs) Gs = π ( q ; L) G s G / s ( ) β p G p

19 9 Averages with Constraints 3. To get correct (unconstrained-system) averages from a simulation using constraints, averages should be multiplied by this factor M unconstained = M G s G constrained Evaluating this quantity could be tedious but there is a simplification the ratio of determinants (of N-by-N and (N-l)-by-(N-l) matrices) can be given in terms of the determinant of an l-by-l matrix G s G = H a β i ri ri for the diatomic with L constrained, H = H αβ = σ σ

20 0 Rotational Dynamics For completely rigid molecules, only translation and rotation are performed Translational dynamics uses methods described previously, but now applied to the COM Rotational dynamics must consider angular velocities and accelerations Can treat via rotation of the molecule-frame coordinates in the spaced-fixed frame!e s =!ω e s Angular velocity changes in angular velocity are given via torque on molecule

21 Quaternions Rate of change of the Euler angles looks like this A problem arises when θ is near 0 no physical significance, but very inconvenient to integration of equations of motion Quaternions can be used to circumvent the problem describe orientation with 4 (non-independent) variables rotation matrix, equations of motion simply expressed in terms of these quantities note: s sinφ cosθ s cosφ cosθ s!φ = ω x ω sinθ y ω sinθ z 0 3 = q q q q q 0 cos θ φ ψ = cos sin θ φ ψ = cos 3 sin θ φ ψ = sin φ ψ 4 cos θ = sin q q q

22 Monte Carlo Rotations MC simulations of molecules include rotation moves must do this to sample orientations of rigid molecules not strictly necessary for non-rigid molecules, but very helpful very easy to do this incorrectly Trial rotation of a linear molecule Let present orientation be given by vector u Generate a unit vector v with random orientation Let new trial orientation be given by new old u = u γ v where γ is a fixed scale factor that sets the size of the perturbation Nonlinear molecule same procedure, but do perturbation on the 4-dimensional vector of quaternions

23 3 Random Vector on a Sphere Acceptance-rejection method of von Neumann Iterate (a) Generate 3 uniform random variates, r, r, r 3 on (0,) (b) Calculate z i = -r i, i=,3, so that the vector z is distributed uniformly in a cube of side, centered on the origin (c) Form the sum z = z z z 3 (d) If z <, take the random vector as (z /z, z /z, z 3 /z) and quit (e) Otherwise, reject the vector and return to (a) Alternative algorithms are possible

Mechanics Physics 151

Mechanics Physics 151 Mechanics Phsics 151 Lecture 8 Rigid Bod Motion (Chapter 4) What We Did Last Time! Discussed scattering problem! Foundation for all experimental phsics! Defined and calculated cross sections! Differential

More information

Rigid Body Rotation. Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li. Department of Applied Mathematics and Statistics Stony Brook University (SUNY)

Rigid Body Rotation. Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li. Department of Applied Mathematics and Statistics Stony Brook University (SUNY) Rigid Body Rotation Speaker: Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University (SUNY) Content Introduction Angular Velocity Angular Momentum

More information

Lecture 37: Principal Axes, Translations, and Eulerian Angles

Lecture 37: Principal Axes, Translations, and Eulerian Angles Lecture 37: Principal Axes, Translations, and Eulerian Angles When Can We Find Principal Axes? We can always write down the cubic equation that one must solve to determine the principal moments But if

More information

Lecture 38: Equations of Rigid-Body Motion

Lecture 38: Equations of Rigid-Body Motion Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can

More information

Lecture 38: Equations of Rigid-Body Motion

Lecture 38: Equations of Rigid-Body Motion Lecture 38: Equations of Rigid-Body Motion It s going to be easiest to find the equations of motion for the object in the body frame i.e., the frame where the axes are principal axes In general, we can

More information

CE 530 Molecular Simulation

CE 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 information

EE Homework 3 Due Date: 03 / 30 / Spring 2015

EE Homework 3 Due Date: 03 / 30 / Spring 2015 EE 476 - Homework 3 Due Date: 03 / 30 / 2015 Spring 2015 Exercise 1 (10 points). Consider the problem of two pulleys and a mass discussed in class. We solved a version of the problem where the mass was

More information

7 To solve numerically the equation of motion, we use the velocity Verlet or leap frog algorithm. _ V i n = F i n m i (F.5) For time step, we approxim

7 To solve numerically the equation of motion, we use the velocity Verlet or leap frog algorithm. _ V i n = F i n m i (F.5) For time step, we approxim 69 Appendix F Molecular Dynamics F. Introduction In this chapter, we deal with the theories and techniques used in molecular dynamics simulation. The fundamental dynamics equations of any system is the

More information

arxiv: v1 [math.ds] 18 Nov 2008

arxiv: v1 [math.ds] 18 Nov 2008 arxiv:0811.2889v1 [math.ds] 18 Nov 2008 Abstract Quaternions And Dynamics Basile Graf basile.graf@epfl.ch February, 2007 We give a simple and self contained introduction to quaternions and their practical

More information

Lecture 13: Forces in the Lagrangian Approach

Lecture 13: Forces in the Lagrangian Approach Lecture 3: Forces in the Lagrangian Approach In regular Cartesian coordinates, the Lagrangian for a single particle is: 3 L = T U = m x ( ) l U xi l= Given this, we can readily interpret the physical significance

More information

M2A2 Problem Sheet 3 - Hamiltonian Mechanics

M2A2 Problem Sheet 3 - Hamiltonian Mechanics MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,

More information

ROBOTICS 01PEEQW. Basilio Bona DAUIN Politecnico di Torino

ROBOTICS 01PEEQW. Basilio Bona DAUIN Politecnico di Torino ROBOTICS 01PEEQW Basilio Bona DAUIN Politecnico di Torino Kinematic Functions Kinematic functions Kinematics deals with the study of four functions(called kinematic functions or KFs) that mathematically

More information

CE 530 Molecular Simulation

CE 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 information

Differential Kinematics

Differential Kinematics Differential Kinematics Relations between motion (velocity) in joint space and motion (linear/angular velocity) in task space (e.g., Cartesian space) Instantaneous velocity mappings can be obtained through

More information

Lesson Rigid Body Dynamics

Lesson Rigid Body Dynamics Lesson 8 Rigid Body Dynamics Lesson 8 Outline Problem definition and motivations Dynamics of rigid bodies The equation of unconstrained motion (ODE) User and time control Demos / tools / libs Rigid Body

More information

Rotations about the coordinate axes are easy to define and work with. Rotation about the x axis by angle θ is. R x (θ) = 0 cos θ sin θ 0 sin θ cos θ

Rotations about the coordinate axes are easy to define and work with. Rotation about the x axis by angle θ is. R x (θ) = 0 cos θ sin θ 0 sin θ cos θ Euler Angle Formulas David Eberly Magic Software, Inc. http://www.magic-software.com Created: December 1, 1999 1 Introduction Rotations about the coordinate axes are easy to define and work with. Rotation

More information

From quantum to classical statistical mechanics. Polyatomic ideal gas.

From quantum to classical statistical mechanics. Polyatomic ideal gas. From quantum to classical statistical mechanics. Polyatomic ideal gas. Peter Košovan peter.kosovan@natur.cuni.cz Dept. of Physical and Macromolecular Chemistry Lecture 5, Statistical Thermodynamics, MC260P105,

More information

Mechanics Physics 151

Mechanics Physics 151 Mechanics Physics 151 Lecture 11 Rigid Body Motion (Chapter 5) Administravia Please fill out the midterm evaluation form Critical feedback for me to evaluate how well (or badly) I m teaching to adjust

More information

Kinematics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Kinematics Semester 1, / 15

Kinematics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Kinematics Semester 1, / 15 Kinematics Basilio Bona DAUIN Politecnico di Torino Semester 1, 2016-17 B. Bona (DAUIN) Kinematics Semester 1, 2016-17 1 / 15 Introduction The kinematic quantities used to represent a body frame are: position

More information

Lecture Notes Multibody Dynamics B, wb1413

Lecture Notes Multibody Dynamics B, wb1413 Lecture Notes Multibody Dynamics B, wb1413 A. L. Schwab & Guido M.J. Delhaes Laboratory for Engineering Mechanics Mechanical Engineering Delft University of Technolgy The Netherlands June 9, 29 Contents

More information

Engineering Physics 1 Dr. B. K. Patra Department of Physics Indian Institute of Technology-Roorkee

Engineering Physics 1 Dr. B. K. Patra Department of Physics Indian Institute of Technology-Roorkee Engineering Physics 1 Dr. B. K. Patra Department of Physics Indian Institute of Technology-Roorkee Module-05 Lecture-04 Maxwellian Distribution Law of Velocity Part 02 So, we have already told to experiment

More information

Kinematics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Kinematics Semester 1, / 15

Kinematics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Kinematics Semester 1, / 15 Kinematics Basilio Bona DAUIN Politecnico di Torino Semester 1, 2014-15 B. Bona (DAUIN) Kinematics Semester 1, 2014-15 1 / 15 Introduction The kinematic quantities used are: position r, linear velocity

More information

A Quantum Mechanical Model for the Vibration and Rotation of Molecules. Rigid Rotor

A Quantum Mechanical Model for the Vibration and Rotation of Molecules. Rigid Rotor A Quantum Mechanical Model for the Vibration and Rotation of Molecules Harmonic Oscillator Rigid Rotor Degrees of Freedom Translation: quantum mechanical model is particle in box or free particle. A molecule

More information

Kinematics. Basilio Bona. October DAUIN - Politecnico di Torino. Basilio Bona (DAUIN - Politecnico di Torino) Kinematics October / 15

Kinematics. Basilio Bona. October DAUIN - Politecnico di Torino. Basilio Bona (DAUIN - Politecnico di Torino) Kinematics October / 15 Kinematics Basilio Bona DAUIN - Politecnico di Torino October 2013 Basilio Bona (DAUIN - Politecnico di Torino) Kinematics October 2013 1 / 15 Introduction The kinematic quantities used are: position r,

More information

SOLUTIONS, PROBLEM SET 11

SOLUTIONS, PROBLEM SET 11 SOLUTIONS, PROBLEM SET 11 1 In this problem we investigate the Lagrangian formulation of dynamics in a rotating frame. Consider a frame of reference which we will consider to be inertial. Suppose that

More information

Physics 200 Lecture 4. Integration. Lecture 4. Physics 200 Laboratory

Physics 200 Lecture 4. Integration. Lecture 4. Physics 200 Laboratory Physics 2 Lecture 4 Integration Lecture 4 Physics 2 Laboratory Monday, February 21st, 211 Integration is the flip-side of differentiation in fact, it is often possible to write a differential equation

More information

Classical Mechanics. Luis Anchordoqui

Classical Mechanics. Luis Anchordoqui 1 Rigid Body Motion Inertia Tensor Rotational Kinetic Energy Principal Axes of Rotation Steiner s Theorem Euler s Equations for a Rigid Body Eulerian Angles Review of Fundamental Equations 2 Rigid body

More information

Lecture 9 Hamilton s Equa4ons

Lecture 9 Hamilton s Equa4ons Lecture 9 Hamilton s Equa4ons Informal deriva5on conjugate momentum cyclic coordinates Applica5ons/examples 1 deriva4on Start with the Euler Lagrange equa5ons d dt L q i L = λ q i j C j 1 + Q i Define

More information

Classical Mechanics Review (Louisiana State University Qualifier Exam)

Classical Mechanics Review (Louisiana State University Qualifier Exam) Review Louisiana State University Qualifier Exam Jeff Kissel October 22, 2006 A particle of mass m. at rest initially, slides without friction on a wedge of angle θ and and mass M that can move without

More information

Numerical Methods for Inverse Kinematics

Numerical Methods for Inverse Kinematics Numerical Methods for Inverse Kinematics Niels Joubert, UC Berkeley, CS184 2008-11-25 Inverse Kinematics is used to pose models by specifying endpoints of segments rather than individual joint angles.

More information

Three-Dimensional Biomechanical Analysis of Human Movement

Three-Dimensional Biomechanical Analysis of Human Movement Three-Dimensional Biomechanical Analysis of Human Movement Anthropometric Measurements Motion Data Acquisition Force Platform Body Mass & Height Biomechanical Model Moments of Inertia and Locations of

More information

Physics 2A Chapter 10 - Rotational Motion Fall 2018

Physics 2A Chapter 10 - Rotational Motion Fall 2018 Physics A Chapter 10 - Rotational Motion Fall 018 These notes are five pages. A quick summary: The concepts of rotational motion are a direct mirror image of the same concepts in linear motion. Follow

More information

Homogeneous Transformations

Homogeneous Transformations Purpose: Homogeneous Transformations The purpose of this chapter is to introduce you to the Homogeneous Transformation. This simple 4 x 4 transformation is used in the geometry engines of CAD systems and

More information

Estimating Accuracy in Classical Molecular Simulation

Estimating Accuracy in Classical Molecular Simulation Estimating Accuracy in Classical Molecular Simulation University of Illinois Urbana-Champaign Department of Computer Science Institute for Mathematics and its Applications July 2007 Acknowledgements Ben

More information

DYNAMICS OF MACHINES

DYNAMICS OF MACHINES DYNAMICS OF MACHINES Modelling Simulation Visualization Verification Ilmar Ferreira Santos, Professor, Dr.-Ing., Dr.Techn., Livre-Docente Department of Mechanical Engineering Technical University of Denmark

More information

A set of N particles forms a rigid body if the distance between any 2 particles is fixed:

A set of N particles forms a rigid body if the distance between any 2 particles is fixed: Chapter Rigid Body Dynamics.1 Coordinates of a Rigid Body A set of N particles forms a rigid body if the distance between any particles is fixed: r ij r i r j = c ij = constant. (.1) Given these constraints,

More information

Fundamental principles

Fundamental principles Dynamics and control of mechanical systems Date Day 1 (03/05) - 05/05 Day (07/05) Day 3 (09/05) Day 4 (11/05) Day 5 (14/05) Day 6 (16/05) Content Review of the basics of mechanics. Kinematics of rigid

More information

NONLINEAR MECHANICAL SYSTEMS (MECHANISMS)

NONLINEAR MECHANICAL SYSTEMS (MECHANISMS) NONLINEAR MECHANICAL SYSTEMS (MECHANISMS) The analogy between dynamic behavior in different energy domains can be useful. Closer inspection reveals that the analogy is not complete. One key distinction

More information

Pendulum with a vibrating base

Pendulum with a vibrating base Pendulum with a vibrating base Chennai, March 14, 2006 1 1 Fast perturbations- rapidly oscillating perturbations We consider perturbations of a Hamiltonian perturbed by rapid oscillations. Later we apply

More information

Lecture Note 1: Background

Lecture Note 1: Background ECE5463: Introduction to Robotics Lecture Note 1: Background Prof. Wei Zhang Department of Electrical and Computer Engineering Ohio State University Columbus, Ohio, USA Spring 2018 Lecture 1 (ECE5463 Sp18)

More information

Lecture 6. Angular Momentum and Roataion Invariance of L. x = xcosφ + ysinφ ' y = ycosφ xsinφ for φ << 1, sinφ φ, cosφ 1 x ' x + yφ y ' = y xφ,

Lecture 6. Angular Momentum and Roataion Invariance of L. x = xcosφ + ysinφ ' y = ycosφ xsinφ for φ << 1, sinφ φ, cosφ 1 x ' x + yφ y ' = y xφ, Lecture 6 Conservation of Angular Momentum (talk about LL p. 9.3) In the previous lecture we saw that time invariance of the L leads to Energy Conservation, and translation invariance of L leads to Momentum

More information

Dynamics 12e. Copyright 2010 Pearson Education South Asia Pte Ltd. Chapter 20 3D Kinematics of a Rigid Body

Dynamics 12e. Copyright 2010 Pearson Education South Asia Pte Ltd. Chapter 20 3D Kinematics of a Rigid Body Engineering Mechanics: Dynamics 12e Chapter 20 3D Kinematics of a Rigid Body Chapter Objectives Kinematics of a body subjected to rotation about a fixed axis and general plane motion. Relative-motion analysis

More information

Rotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep.

Rotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. Rotational Motion Chapter 4 P. J. Grandinetti Chem. 4300 Sep. 1, 2017 P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. 1, 2017 1 / 76 Angular Momentum The angular momentum of a particle with respect

More information

Texture, Microstructure & Anisotropy A.D. (Tony) Rollett

Texture, Microstructure & Anisotropy A.D. (Tony) Rollett 1 Carnegie Mellon MRSEC 27-750 Texture, Microstructure & Anisotropy A.D. (Tony) Rollett Last revised: 5 th Sep. 2011 2 Show how to convert from a description of a crystal orientation based on Miller indices

More information

In most robotic applications the goal is to find a multi-body dynamics description formulated

In most robotic applications the goal is to find a multi-body dynamics description formulated Chapter 3 Dynamics Mathematical models of a robot s dynamics provide a description of why things move when forces are generated in and applied on the system. They play an important role for both simulation

More information

Trajectory-tracking control of a planar 3-RRR parallel manipulator

Trajectory-tracking control of a planar 3-RRR parallel manipulator Trajectory-tracking control of a planar 3-RRR parallel manipulator Chaman Nasa and Sandipan Bandyopadhyay Department of Engineering Design Indian Institute of Technology Madras Chennai, India Abstract

More information

. D CR Nomenclature D 1

. D CR Nomenclature D 1 . D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the

More information

In this section of notes, we look at the calculation of forces and torques for a manipulator in two settings:

In this section of notes, we look at the calculation of forces and torques for a manipulator in two settings: Introduction Up to this point we have considered only the kinematics of a manipulator. That is, only the specification of motion without regard to the forces and torques required to cause motion In this

More information

UAV Coordinate Frames and Rigid Body Dynamics

UAV Coordinate Frames and Rigid Body Dynamics Brigham Young University BYU ScholarsArchive All Faculty Publications 24-- UAV oordinate Frames and Rigid Body Dynamics Randal Beard beard@byu.edu Follow this and additional works at: https://scholarsarchive.byu.edu/facpub

More information

Variation Principle in Mechanics

Variation Principle in Mechanics Section 2 Variation Principle in Mechanics Hamilton s Principle: Every mechanical system is characterized by a Lagrangian, L(q i, q i, t) or L(q, q, t) in brief, and the motion of he system is such that

More information

Using Molecular Dynamics to Compute Properties CHEM 430

Using 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 information

Lecture 11 Overview of Flight Dynamics I. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

Lecture 11 Overview of Flight Dynamics I. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Lecture 11 Overview of Flight Dynamics I Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Point Mass Dynamics Dr. Radhakant Padhi Asst. Professor

More information

Part 8: Rigid Body Dynamics

Part 8: Rigid Body Dynamics Document that contains homework problems. Comment out the solutions when printing off for students. Part 8: Rigid Body Dynamics Problem 1. Inertia review Find the moment of inertia for a thin uniform rod

More information

Ridig Body Motion Homogeneous Transformations

Ridig Body Motion Homogeneous Transformations Ridig Body Motion Homogeneous Transformations Claudio Melchiorri Dipartimento di Elettronica, Informatica e Sistemistica (DEIS) Università di Bologna email: claudio.melchiorri@unibo.it C. Melchiorri (DEIS)

More information

Physics 5153 Classical Mechanics. Canonical Transformations-1

Physics 5153 Classical Mechanics. Canonical Transformations-1 1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant

More information

Multibody simulation

Multibody simulation Multibody simulation Dynamics of a multibody system (Newton-Euler formulation) Dimitar Dimitrov Örebro University June 8, 2012 Main points covered Newton-Euler formulation forward dynamics inverse dynamics

More information

CE 530 Molecular Simulation

CE 530 Molecular Simulation CE 530 Molecular Simulation Lecture 20 Phase Equilibria David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Thermodynamic Phase Equilibria Certain thermodynamic states

More information

Lecture 7 Rolling Constraints

Lecture 7 Rolling Constraints Lecture 7 Rolling Constraints The most common, and most important nonholonomic constraints They cannot be wri5en in terms of the variables alone you must include some deriva9ves The resul9ng differen9al

More information

Physics 339 Euler Angles & Free Precession November 2017 Hamilton s Revenge

Physics 339 Euler Angles & Free Precession November 2017 Hamilton s Revenge Physics 339 Euler Angles & Free Precession November 217 Hamilton s Revenge As described in the textbook, Euler Angles are a way to specify the configuration of a 3d object. Starting from a fixed configuration

More information

Mathematical Modeling and response analysis of mechanical systems are the subjects of this chapter.

Mathematical Modeling and response analysis of mechanical systems are the subjects of this chapter. Chapter 3 Mechanical Systems A. Bazoune 3.1 INRODUCION Mathematical Modeling and response analysis of mechanical systems are the subjects of this chapter. 3. MECHANICAL ELEMENS Any mechanical system consists

More information

Chapters 10 & 11: Rotational Dynamics Thursday March 8 th

Chapters 10 & 11: Rotational Dynamics Thursday March 8 th Chapters 10 & 11: Rotational Dynamics Thursday March 8 th Review of rotational kinematics equations Review and more on rotational inertia Rolling motion as rotation and translation Rotational kinetic energy

More information

Aircraft Equations of Motion: Translation and Rotation Robert Stengel, Aircraft Flight Dynamics, MAE 331, 2018

Aircraft Equations of Motion: Translation and Rotation Robert Stengel, Aircraft Flight Dynamics, MAE 331, 2018 Aircraft Equations of Motion: Translation and Rotation Robert Stengel, Aircraft Flight Dynamics, MAE 331, 2018 Learning Objectives What use are the equations of motion? How is the angular orientation of

More information

Conservation of Linear Momentum : If a force F is acting on particle of mass m, then according to Newton s second law of motion, we have F = dp /dt =

Conservation of Linear Momentum : If a force F is acting on particle of mass m, then according to Newton s second law of motion, we have F = dp /dt = Conservation of Linear Momentum : If a force F is acting on particle of mass m, then according to Newton s second law of motion, we have F = dp /dt = d (mv) /dt where p =mv is linear momentum of particle

More information

Dynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics

Dynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics for resubmission Dynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics Aram Davtyan, 1 Gregory A. Voth, 1 2, a) and Hans

More information

QCD at hadron colliders Lecture 2: Showers, Jets and fixed-order predictions

QCD at hadron colliders Lecture 2: Showers, Jets and fixed-order predictions QCD at hadron colliders Lecture 2: Showers, Jets and fixed-order predictions Gavin Salam CERN, Princeton & LPTHE/CNRS (Paris) Maria Laach Herbtschule für Hochenenergiephysik September 20, Germany QCD lecture

More information

Chemistry 795T. NC State University. Lecture 4. Vibrational and Rotational Spectroscopy

Chemistry 795T. NC State University. Lecture 4. Vibrational and Rotational Spectroscopy Chemistry 795T Lecture 4 Vibrational and Rotational Spectroscopy NC State University The Dipole Moment Expansion The permanent dipole moment of a molecule oscillates about an equilibrium value as the molecule

More information

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Rigid body physics Particle system Most simple instance of a physics system Each object (body) is a particle Each particle

More information

CHAPTER 8 The Quantum Theory of Motion

CHAPTER 8 The Quantum Theory of Motion I. Translational motion. CHAPTER 8 The Quantum Theory of Motion A. Single particle in free space, 1-D. 1. Schrodinger eqn H ψ = Eψ! 2 2m d 2 dx 2 ψ = Eψ ; no boundary conditions 2. General solution: ψ

More information

Physical Dynamics (PHY-304)

Physical Dynamics (PHY-304) Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.

More information

Massachusetts Institute of Technology Department of Physics. Final Examination December 17, 2004

Massachusetts Institute of Technology Department of Physics. Final Examination December 17, 2004 Massachusetts Institute of Technology Department of Physics Course: 8.09 Classical Mechanics Term: Fall 004 Final Examination December 17, 004 Instructions Do not start until you are told to do so. Solve

More information

Advanced Dynamics. - Lecture 4 Lagrange Equations. Paolo Tiso Spring Semester 2017 ETH Zürich

Advanced Dynamics. - Lecture 4 Lagrange Equations. Paolo Tiso Spring Semester 2017 ETH Zürich Advanced Dynamics - Lecture 4 Lagrange Equations Paolo Tiso Spring Semester 2017 ETH Zürich LECTURE OBJECTIVES 1. Derive the Lagrange equations of a system of particles; 2. Show that the equation of motion

More information

Physics 235 Chapter 7. Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics

Physics 235 Chapter 7. Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics Many interesting physics systems describe systems of particles on which many forces are acting. Some of these forces are immediately

More information

Quaternion based Extended Kalman Filter

Quaternion based Extended Kalman Filter Quaternion based Extended Kalman Filter, Sergio Montenegro About this lecture General introduction to rotations and quaternions. Introduction to Kalman Filter for Attitude Estimation How to implement and

More information

8 Velocity Kinematics

8 Velocity Kinematics 8 Velocity Kinematics Velocity analysis of a robot is divided into forward and inverse velocity kinematics. Having the time rate of joint variables and determination of the Cartesian velocity of end-effector

More information

Phys 7221 Homework # 8

Phys 7221 Homework # 8 Phys 71 Homework # 8 Gabriela González November 15, 6 Derivation 5-6: Torque free symmetric top In a torque free, symmetric top, with I x = I y = I, the angular velocity vector ω in body coordinates with

More information

Quaternions and their applications

Quaternions and their applications Quaternions and their applications Ramachandran Subramanian February 20, 2014 Jack B. Kuipers. Quaternions and Rotation Sequences. university press, New Jersey, 1998 Princeton Outline 1 Background 2 Introduction

More information

CE 530 Molecular Simulation

CE 530 Molecular Simulation CE 530 Molecular Simulation Lecture 0 Simple Biasing Methods David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo simulation Markov chain to generate

More information

Optimal Control, Guidance and Estimation. Lecture 17. Overview of Flight Dynamics III. Prof. Radhakant Padhi. Prof.

Optimal Control, Guidance and Estimation. Lecture 17. Overview of Flight Dynamics III. Prof. Radhakant Padhi. Prof. Optimal Control, Guidance and Estimation Lecture 17 Overview of Flight Dynamics III Prof. adhakant Padhi Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Six DOF Model Prof. adhakant

More information

1.1.1 Bell Inequality - Spin correlation

1.1.1 Bell Inequality - Spin correlation January 8, 015 Lecture IV 1.1.1 Bell Inequality - Spin correlation Consider the final spin singlet state of the decay η 0 µ + µ We suppose that the η 0 decays and the muon and µ + travel in opposite directions,

More information

Curves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,

Curves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations, Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal

More information

Modern Siege Weapons: Mechanics of the Trebuchet

Modern Siege Weapons: Mechanics of the Trebuchet 1/35 Modern Siege Weapons: Mechanics of the Trebuchet Shawn Rutan and Becky Wieczorek The Project The purpose of this project is to study three simplifications of a trebuchet. We will analyze the motion

More information

d 3 r d 3 vf( r, v) = N (2) = CV C = n where n N/V is the total number of molecules per unit volume. Hence e βmv2 /2 d 3 rd 3 v (5)

d 3 r d 3 vf( r, v) = N (2) = CV C = n where n N/V is the total number of molecules per unit volume. Hence e βmv2 /2 d 3 rd 3 v (5) LECTURE 12 Maxwell Velocity Distribution Suppose we have a dilute gas of molecules, each with mass m. If the gas is dilute enough, we can ignore the interactions between the molecules and the energy will

More information

Lecture 8: Rolling Constraints II

Lecture 8: Rolling Constraints II Lecture 8: Rolling Constraints II Generaliza3ons and review Nonhorizontal surfaces Coin rolling on a slope Small sphere rolling on a larger sphere s surface Hoop rolling inside a hoop 1 What can we say

More information

DIFFERENTIAL KINEMATICS. Geometric Jacobian. Analytical Jacobian. Kinematic singularities. Kinematic redundancy. Inverse differential kinematics

DIFFERENTIAL KINEMATICS. Geometric Jacobian. Analytical Jacobian. Kinematic singularities. Kinematic redundancy. Inverse differential kinematics DIFFERENTIAL KINEMATICS relationship between joint velocities and end-effector velocities Geometric Jacobian Analytical Jacobian Kinematic singularities Kinematic redundancy Inverse differential kinematics

More information

Dynamic Model of a Badminton Stroke

Dynamic Model of a Badminton Stroke ISEA 28 CONFERENCE Dynamic Model of a Badminton Stroke M. Kwan* and J. Rasmussen Department of Mechanical Engineering, Aalborg University, 922 Aalborg East, Denmark Phone: +45 994 9317 / Fax: +45 9815

More information

Robotics & Automation. Lecture 25. Dynamics of Constrained Systems, Dynamic Control. John T. Wen. April 26, 2007

Robotics & Automation. Lecture 25. Dynamics of Constrained Systems, Dynamic Control. John T. Wen. April 26, 2007 Robotics & Automation Lecture 25 Dynamics of Constrained Systems, Dynamic Control John T. Wen April 26, 2007 Last Time Order N Forward Dynamics (3-sweep algorithm) Factorization perspective: causal-anticausal

More information

Nonholonomic Constraints Examples

Nonholonomic Constraints Examples Nonholonomic Constraints Examples Basilio Bona DAUIN Politecnico di Torino July 2009 B. Bona (DAUIN) Examples July 2009 1 / 34 Example 1 Given q T = [ x y ] T check that the constraint φ(q) = (2x + siny

More information

Physics Fall Mechanics, Thermodynamics, Waves, Fluids. Lecture 20: Rotational Motion. Slide 20-1

Physics Fall Mechanics, Thermodynamics, Waves, Fluids. Lecture 20: Rotational Motion. Slide 20-1 Physics 1501 Fall 2008 Mechanics, Thermodynamics, Waves, Fluids Lecture 20: Rotational Motion Slide 20-1 Recap: center of mass, linear momentum A composite system behaves as though its mass is concentrated

More information

Lecture D16-2D Rigid Body Kinematics

Lecture D16-2D Rigid Body Kinematics J. Peraire 16.07 Dynamics Fall 2004 Version 1.2 Lecture D16-2D Rigid Body Kinematics In this lecture, we will start from the general relative motion concepts introduced in lectures D11 and D12, and then

More information

Introduction to Electromagnetism Prof. Manoj K. Harbola Department of Physics Indian Institute of Technology, Kanpur

Introduction to Electromagnetism Prof. Manoj K. Harbola Department of Physics Indian Institute of Technology, Kanpur Introduction to Electromagnetism Prof. Manoj K. Harbola Department of Physics Indian Institute of Technology, Kanpur Lecture - 12 Line surface area and volume elements in Spherical Polar Coordinates In

More information

2 Canonical quantization

2 Canonical quantization Phys540.nb 7 Canonical quantization.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system?.1.1.lagrangian Lagrangian mechanics is a reformulation of classical mechanics.

More information

Spacecraft and Aircraft Dynamics

Spacecraft and Aircraft Dynamics Spacecraft and Aircraft Dynamics Matthew M. Peet Illinois Institute of Technology Lecture 4: Contributions to Longitudinal Stability Aircraft Dynamics Lecture 4 In this lecture, we will discuss Airfoils:

More information

Introduction and Vectors Lecture 1

Introduction and Vectors Lecture 1 1 Introduction Introduction and Vectors Lecture 1 This is a course on classical Electromagnetism. It is the foundation for more advanced courses in modern physics. All physics of the modern era, from quantum

More information

Lecture Note 12: Dynamics of Open Chains: Lagrangian Formulation

Lecture Note 12: Dynamics of Open Chains: Lagrangian Formulation ECE5463: Introduction to Robotics Lecture Note 12: Dynamics of Open Chains: Lagrangian Formulation Prof. Wei Zhang Department of Electrical and Computer Engineering Ohio State University Columbus, Ohio,

More information

Robotics. Dynamics. Marc Toussaint U Stuttgart

Robotics. Dynamics. Marc Toussaint U Stuttgart Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, Euler-Lagrange equation, Newton-Euler recursion, general robot dynamics, joint space control, reference trajectory

More information

Physical Dynamics (SPA5304) Lecture Plan 2018

Physical Dynamics (SPA5304) Lecture Plan 2018 Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle

More information

Case Study: The Pelican Prototype Robot

Case Study: The Pelican Prototype Robot 5 Case Study: The Pelican Prototype Robot The purpose of this chapter is twofold: first, to present in detail the model of the experimental robot arm of the Robotics lab. from the CICESE Research Center,

More information

Canonical transformations (Lecture 4)

Canonical transformations (Lecture 4) Canonical transformations (Lecture 4) January 26, 2016 61/441 Lecture outline We will introduce and discuss canonical transformations that conserve the Hamiltonian structure of equations of motion. Poisson

More information

Manipulator Dynamics 2. Instructor: Jacob Rosen Advanced Robotic - MAE 263D - Department of Mechanical & Aerospace Engineering - UCLA

Manipulator Dynamics 2. Instructor: Jacob Rosen Advanced Robotic - MAE 263D - Department of Mechanical & Aerospace Engineering - UCLA Manipulator Dynamics 2 Forward Dynamics Problem Given: Joint torques and links geometry, mass, inertia, friction Compute: Angular acceleration of the links (solve differential equations) Solution Dynamic

More information

Fun With Carbon Monoxide. p. 1/2

Fun With Carbon Monoxide. p. 1/2 Fun With Carbon Monoxide p. 1/2 p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results p. 1/2 Fun With Carbon Monoxide E = 0.25 ± 0.05 ev Electron beam results C V (J/K-mole) 35 30 25

More information