CE 530 Molecular Simulation


 Alicia Hart
 11 months ago
 Views:
Transcription
1 CE 530 Molecular Simulation Lecture 7 Beyond Atoms: Simulating Molecules David A. Kofke Department of Chemical Engineering SUNY Buffalo
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  Potentialtruncation schemes must be careful not to split the charges For a 9Å truncation distance, using waterlike 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 VolumeScaling Moves Scaling atom displacements leads to large strain on intramolecular bonds Instead perform volume scaling moves using molecule centersofmass (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 intramolecular 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 rigidbond 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 rigidbody moves (even if model does not have rigid bonds) but much care is needed to do this properly
6 Molecule Coordinate Frame 6 Moleculeframe coordinates are defined w.r.t. molecule COM with molecule in a reference orientation Simulationframe coordinate is determined by molecule COM R and orientation ω q j R i For rigid molecules, the moleculeframe 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 rigidmolecule 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 simulationframe 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 simulationframe 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 spacefixed coordinate, multiply moleculefixed 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 spacefixed 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 transitionprobability 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 constrainedsystem 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 (unconstrainedsystem) 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 NbyN and (Nl)by(Nl) matrices) can be given in terms of the determinant of an lbyl 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 moleculeframe coordinates in the spacedfixed 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 (nonindependent) 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 nonrigid 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 4dimensional vector of quaternions
23 3 Random Vector on a Sphere Acceptancerejection 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
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 informationRidig 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 informationLecture 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 informationDynamic 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 informationRobotics. Dynamics. Marc Toussaint U Stuttgart
Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, EulerLagrange equation, NewtonEuler recursion, general robot dynamics, joint space control, reference trajectory
More informationOptimal Control, Guidance and Estimation. Lecture 16. Overview of Flight Dynamics II. Prof. Radhakant Padhi. Prof. Radhakant Padhi
Optimal Control, Guidance and Estimation Lecture 16 Overview of Flight Dynamics II Prof. Radhakant Padhi Dept. of erospace Engineering Indian Institute of Science  Bangalore Point Mass Dynamics Prof.
More informationA Brief Introduction to Statistical Mechanics
A Brief Introduction to Statistical Mechanics E. J. Maginn, J. K. Shah Department of Chemical and Biomolecular Engineering University of Notre Dame Notre Dame, IN 46556 USA Monte Carlo Workshop Universidade
More informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
More informationDynamics. 1 Copyright c 2015 Roderic Grupen
Dynamics The branch of physics that treats the action of force on bodies in motion or at rest; kinetics, kinematics, and statics, collectively. Websters dictionary Outline Conservation of Momentum Inertia
More informationSPACECRAFT ATTITUDE PROPAGATIONWITH DIFFERENT REPRESENTATIONS
INPE93PRE/673 SPACECRAFT ATTITUDE PROPAGATIONWITH DIFFERENT REPRESENTATIONS Domingos Sávio dos Santos Rodrigues Maria Cecília Zanardi ADVANCES IN SPACE DYNAMICS : CELESTIAL MECHANICS AND ASTRONAUTICS,
More informationDynamics of spinning particles in Schwarzschild spacetime
Dynamics of spinning particles in Schwarzschild spacetime, Volker Perlick, Claus Lämmerzahl Center of Space Technology and Microgravity University of Bremen, Germany 08.05.2014 RTG Workshop, Bielefeld
More informationGeneralization of the Euler Angles
The Journal of the Astronautical Sciences, Vol. 51, No. 2, AprilJune 2003, pp. 123123 Generalization of the Euler Angles Malcolm D. Shuster 1 and F. Landis Markley 2 Abstract It is shown that the Euler
More informationeigenvalues eigenfunctions
BornOppenheimer Approximation Atoms and molecules consist of heavy nuclei and light electrons. Consider (for simplicity) a diatomic molecule (e.g. HCl). Clamp/freeze the nuclei in space, a distance r
More informationTransformations. Chapter D Transformations Translation
Chapter 4 Transformations Transformations between arbitrary vector spaces, especially linear transformations, are usually studied in a linear algebra class. Here, we focus our attention to transformation
More informationCBE 6333, R. Levicky 1. Orthogonal Curvilinear Coordinates
CBE 6333, R. Levicky 1 Orthogonal Curvilinear Coordinates Introduction. Rectangular Cartesian coordinates are convenient when solving problems in which the geometry of a problem is well described by the
More informationBrief Review of Statistical Mechanics
Brief Review of Statistical Mechanics Introduction Statistical mechanics: a branch of physics which studies macroscopic systems from a microscopic or molecular point of view (McQuarrie,1976) Also see (Hill,1986;
More informationPenning Traps. Contents. Plasma Physics Penning Traps AJW August 16, Introduction. Clasical picture. Radiation Damping.
Penning Traps Contents Introduction Clasical picture Radiation Damping Number density B and E fields used to increase time that an electron remains within a discharge: Penning, 936. Can now trap a particle
More informationGraphical Display of a Physics Simulation
Graphical Display of a Physics Simulation Steven Oetjen TJHSST Computer Systems Lab, 20062007 March 29, 2007 Abstract A physics simulation, in order to adequately demonstrate physical laws and predict
More informationChapter 1. Principles of Motion in Invariantive Mechanics
Chapter 1 Principles of Motion in Invariantive Mechanics 1.1. The EulerLagrange and Hamilton s equations obtained by means of exterior forms Let L = L(q 1,q 2,...,q n, q 1, q 2,..., q n,t) L(q, q,t) (1.1)
More informationIntroduction to Vibration. Mike Brennan UNESP, Ilha Solteira São Paulo Brazil
Introduction to Vibration Mike Brennan UNESP, Ilha Solteira São Paulo Brazil Vibration Most vibrations are undesirable, but there are many instances where vibrations are useful Ultrasonic (very high
More informationNAWCWPNS TM 8128 CONTENTS. Introduction TwoDimensinal Motion ThreeDimensional Motion Nonrotating Spherical Earth...
CONTENTS Introduction... 3 TwoDimensinal Motion... 3 ThreeDimensional Motion... 5 Nonrotating Spherical Earth...10 Rotating Spherical Earth...12 WGS84...14 Conclusion...14 Appendixes: A. Kalman Filter...15
More informationTime dilation and length contraction without relativity: The Bohr atom and the semiclassical hydrogen molecule ion
Time dilation and length contraction without relativity: The Bohr atom and the semiclassical hydrogen molecule ion Allan Walstad * Abstract Recently it was shown that classical "relativistic" particle
More informationPath Integral Quantization of the Electromagnetic Field Coupled to A Spinor
EJTP 6, No. 22 (2009) 189 196 Electronic Journal of Theoretical Physics Path Integral Quantization of the Electromagnetic Field Coupled to A Spinor Walaa. I. Eshraim and Nasser. I. Farahat Department of
More informationCase Study: Protein Folding using Homotopy Methods
Case Study: Protein Folding using Homotopy Methods james robert white Since the late 1990s, the world of genomics has exploded. Beginning with the bacterium Haemophilus influenzae in 1995, researchers
More informationThis module requires you to read a textbook such as Fowles and Cassiday on material relevant to the following topics.
Module M2 Lagrangian Mechanics and Oscillations Prerequisite: Module C1 This module requires you to read a textbook such as Fowles and Cassiday on material relevant to the following topics. Topics: Hamilton
More informationAccelerator Physics Homework #7 P470 (Problems: 14)
Accelerator Physics Homework #7 P470 (Problems: 4) This exercise derives the linear transfer matrix for a skew quadrupole, where the magnetic field is B z = B 0 a z, B x = B 0 a x, B s = 0; with B 0 a
More informationRigid Body Motion. Greg Hager Simon Leonard
Rigid ody Motion Greg Hager Simon Leonard Overview Different spaces used in robotics and why we need to get from one space to the other Focus on Cartesian space Transformation between two Cartesian coordinate
More informationChapter 3 Kinematics in Two Dimensions; Vectors
Chapter 3 Kinematics in Two Dimensions; Vectors Vectors and Scalars Units of Chapter 3 Addition of Vectors Graphical Methods Subtraction of Vectors, and Multiplication of a Vector by a Scalar Adding Vectors
More informationSpan and Linear Independence
Span and Linear Independence It is common to confuse span and linear independence, because although they are different concepts, they are related. To see their relationship, let s revisit the previous
More informationProtein Dynamics, Allostery and Function
Protein Dynamics, Allostery and Function Lecture 3. Protein Dynamics Xiaolin Cheng UT/ORNL Center for Molecular Biophysics SJTU Summer School 2017 1 Obtaining Dynamic Information Experimental Approaches
More informationProject: Vibration of Diatomic Molecules
Project: Vibration of Diatomic Molecules Objective: Find the vibration energies and the corresponding vibrational wavefunctions of diatomic molecules H 2 and I 2 using the Morse potential. equired Numerical
More informationPhysics 4A Solutions to Chapter 11 Homework
Physics 4A Solutions to Chapter 11 Homework Chapter 11 Questions:, 8, 10 Exercises & Problems: 1, 14, 4, 7, 37, 53, 66, 81, 83 Answers to Questions: Q 11 (a) 5 and 6 (b) 1 and 4 tie, then the rest tie
More informationFigure 5.16 Compound pendulum: (a) At rest in equilibrium, (b) General position with coordinate θ, Freebody
Lecture 27. THE COMPOUND PENDULUM Figure 5.16 Compound pendulum: (a) At rest in equilibrium, (b) General position with coordinate θ, Freebody diagram The term compound is used to distinguish the present
More informationLecture 21. MORE PLANAR KINEMATIC EXAMPLES
Lecture 21. MORE PLANAR KINEMATIC EXAMPLES 4.5c Another SliderCrank Mechanism Figure 4.24 Alternative slidercrank mechanism. Engineeringanalysis task: For = ω = constant, determine φ and S and their
More informationMathematics for 3D Graphics
math 1 Topics Mathematics for 3D Graphics math 1 Points, Vectors, Vertices, Coordinates Dot Products, Cross Products Lines, Planes, Intercepts References Many texts cover the linear algebra used for 3D
More informationDispersion relation for transverse waves in a linear chain of particles
Dispersion relation for transverse waves in a linear chain of particles V. I. Repchenkov* It is difficult to overestimate the importance that have for the development of science the simplest physical and
More informationPractice Test 3. Multiple Choice Identify the choice that best completes the statement or answers the question.
Practice Test 3 Multiple Choice Identify the choice that best completes the statement or answers the question. 1. A wheel rotates about a fixed axis with an initial angular velocity of 20 rad/s. During
More informationExamination paper for TMA4195 Mathematical Modeling
Department of Mathematical Sciences Examination paper for TMA4195 Mathematical Modeling Academic contact during examination: Elena Celledoni Phone: 48238584, 73593541 Examination date: 11th of December
More informationChapter 5. . Dynamics. 5.1 Introduction
Chapter 5. Dynamics 5.1 Introduction The study of manipulator dynamics is essential for both the analysis of performance and the design of robot control. A manipulator is a multilink, highly nonlinear
More informationMath 234 Exam 3 Review Sheet
Math 234 Exam 3 Review Sheet Jim Brunner LIST OF TOPIS TO KNOW Vector Fields lairaut s Theorem & onservative Vector Fields url Divergence Area & Volume Integrals Using oordinate Transforms hanging the
More informationEJL R Sβ. sum. General objective: Reduce the complexity of the analysis by exploiting symmetry. Garth J. Simpson
e sum = EJL R Sβ General objective: Reduce the complexity of the analysis by exploiting symmetry. Specific Objectives: 1. The molecular symmetry matrix S. How to populate it.. Relationships between the
More informationBIODYNAMICS: A LAGRANGIAN APPROACH
Source: STANDARD HANDBOOK OF BIOMEDICAL ENGINEERING AND DESIGN CHAPTER 7 BIODYNAMICS: A LAGRANGIAN APPROACH Donald R. Peterson Biodynamics Laboratory at the Ergonomic Technology Center, University of Connecticut
More informationReconstructing inclusions from Electrostatic Data
Reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 778433368 iharris@math.tamu.edu Joint work with: W. Rundell Purdue
More informationLAWS OF GYROSCOPES / CARDANIC GYROSCOPE
LAWS OF GYROSCOPES / CARDANC GYROSCOPE PRNCPLE f the axis of rotation of the forcefree gyroscope is displaced slightly, a nutation is produced. The relationship between precession frequency or nutation
More information8.821 String Theory Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 8.82 String Theory Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 8.82 F2008 Lecture 24 Blackhole Thermodynamics
More informationSlide 1 / 37. Rotational Motion
Slide 1 / 37 Rotational Motion Slide 2 / 37 Angular Quantities An angle θ can be given by: where r is the radius and l is the arc length. This gives θ in radians. There are 360 in a circle or 2π radians.
More informationGeneral Physics I. Lecture 8: Rotation of a Rigid Object About a Fixed Axis. Prof. WAN, Xin ( 万歆 )
General Physics I Lecture 8: Rotation of a Rigid Object About a Fixed Axis Prof. WAN, Xin ( 万歆 ) xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ New Territory Object In the past, point particle (no rotation,
More informationLecture 3. Electric Field Flux, Gauss Law
Lecture 3. Electric Field Flux, Gauss Law Attention: the list of unregistered iclickers will be posted on our Web page after this lecture. From the concept of electric field flux to the calculation of
More informationChapter 8. Rigid transformations
Chapter 8. Rigid transformations We are about to start drawing figures in 3D. There are no builtin routines for this purpose in PostScript, and we shall have to start more or less from scratch in extending
More informationLecture 10 Diatomic Vibration Spectra Harmonic Model
Chemistry II: Introduction to Molecular Spectroscopy Prof. Mangala Sunder Department of Chemistry and Biochemistry Indian Institute of Technology, Madras Lecture 10 Diatomic Vibration Spectra Harmonic
More informationGame Physics: Basic Concepts
CMSC 498M: Chapter 8a Game Physics Reading: Game Physics, by David H Eberly, 2004 Physics for Game Developers, by David M Bourg, 2002 Overview: Basic physical quantities: Mass, center of mass, moment of
More informationReview of last lecture I
Review of last lecture I t f + v x f = Q γ (f,f), (1) where f (x,v,t) is a nonegative function that represents the density of particles in position x at time t with velocity v. In (1) Q γ is the socalled
More informationThis semester. Books
Models mostly proteins from detailed to more abstract models Some simulation methods This semester Books None necessary for my group and Prof Rarey Molecular Modelling: Principles and Applications Leach,
More informationA MULTIBODY ALGORITHM FOR WAVE ENERGY CONVERTERS EMPLOYING NONLINEAR JOINT REPRESENTATION
Proceedings of the ASME 2014 33rd International Conference on Ocean, Offshore and Arctic Engineering OMAE2014 June 813, 2014, San Francisco, California, USA OMAE201423864 A MULTIBODY ALGORITHM FOR WAVE
More informationParticle Systems. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2017
Particle Systems CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2017 Particle Systems Particle systems have been used extensively in computer animation and special effects since their
More informationPhysics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text:
Physics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text: 1.3 1.6 Constraints Often times we consider dynamical systems which are defined using some kind of restrictions
More informationLarge deflection analysis of planar solids based on the Finite Particle Method
yuying@uiuc.edu 10 th US National Congress on Computational Mechanics Large deflection analysis of planar solids based on the Finite Particle Method 1, 2 Presenter: Ying Yu Advisors: Prof. Glaucio H. Paulino
More informationLectures on basic plasma physics: Hamiltonian mechanics of charged particle motion
Lectures on basic plasma physics: Hamiltonian mechanics of charged particle motion Department of applied physics, Aalto University March 8, 2016 Hamiltonian versus Newtonian mechanics Newtonian mechanics:
More informationMatrix Elements and Shower CKKW and MLM
Typeset with TeXmacs Matrix Elements and Shower CKKW and MLM P. Nason INFN, Sez. of Milano Bicocca Cern, 1182009 1 Historical approach: CKKW Catani, Krauss, Kuhn, Webber (2001), (in e + e annihilation).
More informationElastic Scattering. R = m 1r 1 + m 2 r 2 m 1 + m 2. is the center of mass which is known to move with a constant velocity (see previous lectures):
Elastic Scattering In this section we will consider a problem of scattering of two particles obeying Newtonian mechanics. The problem of scattering can be viewed as a truncated version of dynamic problem
More informationRoots and Coefficients of a Quadratic Equation Summary
Roots and Coefficients of a Quadratic Equation Summary For a quadratic equation with roots α and β: Sum of roots = α + β = and Product of roots = αβ = Symmetrical functions of α and β include: x = and
More informationLecture 35: FRI 17 APR Electrical Oscillations, LC Circuits, Alternating Current I
Physics 3 Jonathan Dowling Lecture 35: FRI 7 APR Electrical Oscillations, LC Circuits, Alternating Current I Nikolai Tesla What are we going to learn? A road map Electric charge è Electric force on other
More informationExperiment 08: Physical Pendulum. 8.01t Nov 10, 2004
Experiment 08: Physical Pendulum 8.01t Nov 10, 2004 Goals Investigate the oscillation of a real (physical) pendulum and compare to an ideal (point mass) pendulum. Angular frequency calculation: Practice
More informationThe Notsosimple Pendulum: Balancing a Pencil on its Point Peter Lynch, UCD, Dublin, May 2014
The Notsosimple Pendulum: Balancing a Pencil on its Point Peter Lynch, UCD, Dublin, May 204 arxiv:406.25v [nlin.si] 4 Jun 204 ABSTRACT. Does quantum mechanics matter at everyday scales? We generally
More informationEnergy Bounds in Linear Elasticity
Chapter 1 Energy Bounds in Linear Elasticity Focus of this Recitation: Thefocusofthisrecitationistoputtoworktheenergyboundsoflinearelasticity for the derivation of upper bound and lower bound energy solutions
More informationMechanics Physics 151
Mechanics Physics 151 Fall 003 Masahiro Morii Teaching Staff! Lecturer: Masahiro Morii! Tuesday/Thursday 11:30 1:00. Jefferson 56! Section leaders: Srinivas Paruchuri and AbdolReza Mansouri! Two or three
More informationNotes on multivariable calculus
Notes on multivariable calculus Jonathan Wise February 2, 2010 1 Review of trigonometry Trigonometry is essentially the study of the relationship between polar coordinates and Cartesian coordinates in
More informationSection 5.4 The Other Trigonometric Functions
Section 5.4 The Other Trigonometric Functions Section 5.4 The Other Trigonometric Functions In the previous section, we defined the e and coe functions as ratios of the sides of a right triangle in a circle.
More informationOvercoming the Limitations of Conservation of Linear Momentum by Including External Impulse
Overcoming the Limitations of Conservation of Linear Momentum by Including External Impulse March 12, 2008 The University of Tulsa MEPS Lunch Meeting 1 / 40 Motivation The University of Tulsa MEPS Lunch
More informationExam 3 Practice Solutions
Exam 3 Practice Solutions Multiple Choice 1. A thin hoop, a solid disk, and a solid sphere, each with the same mass and radius, are at rest at the top of an inclined plane. If all three are released at
More informationDynamic Response of Structures With Frequency Dependent Damping
Dynamic Response of Structures With Frequency Dependent Damping Blanca Pascual & S Adhikari School of Engineering, Swansea University, Swansea, UK Email: S.Adhikari@swansea.ac.uk URL: http://engweb.swan.ac.uk/
More informationQuiz Discussion. IE417: Nonlinear Programming: Lecture 12. Motivation. Why do we care? Jeff Linderoth. 16th March 2006
Quiz Discussion IE417: Nonlinear Programming: Lecture 12 Jeff Linderoth Department of Industrial and Systems Engineering Lehigh University 16th March 2006 Motivation Why do we care? We are interested in
More informationGEOMETRIC QUANTIZATION
GEOMETRIC QUANTIZATION 1. The basic idea The setting of the Hamiltonian version of classical (Newtonian) mechanics is the phase space (position and momentum), which is a symplectic manifold. The typical
More informationNormal Force. W = mg cos(θ) Normal force F N = mg cos(θ) F N
Normal Force W = mg cos(θ) Normal force F N = mg cos(θ) Note there is no weight force parallel/down the include. The car is not pressing on anything causing a force in that direction. If there were a person
More informationΑΜ Α and ΒΜ Β angular momentum basis states to form coupled ΑΒCΜ C basis states RWF Lecture #4. The WignerEckart Theorem
MIT Department of Chemistry 5.74, Spring 2004: Introductory Quantum Mechanics II Instructor: Prof. Robert Field 5.74 RWF Lecture #4 4 The WignerEcart Theorem It is always possible to evaluate the angular
More informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 3 The BornOppenheimer approximation C.K. Skylaris Learning outcomes Separate molecular Hamiltonians to electronic and nuclear parts according to the BornOppenheimer
More informationBACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigidbody spring deformation expressed in terms of field variables assembly of rigidbodies connected by spring Distinct Element Method (DEM) simple
More informationECE 145A/218A Power Amplifier Design Lectures. Power Amplifier Design 1
Power Amplifiers; Part 1 Class A Device Limitations Large signal output match Define efficiency, poweradded efficiency Class A operating conditions Thermal resistance We have studied the design of smallsignal
More informationMODEL ANSWERS TO HWK #1
MODEL ANSWERS TO HWK # Part B (a The four vertices are (,,, (,,, (,, and (,, The distance between the first two vertices is, since two coordinates differ by There are six edges, corresponding to the choice
More informationExperimental Tests and Alternative Theories of Gravity
Experimental Tests and Alternative Theories of Gravity Gonzalo J. Olmo Alba gonzalo.olmo@uv.es University of Valencia (Spain) & UWMilwaukee Experimental Tests and Alternative Theories of Gravity p. 1/2
More informationThe dynamics of a Mobile Inverted Pendulum (MIP)
The dynamics of a Mobile Inverted Pendulum (MIP) 1 Introduction Saam Ostovari, Nick Morozovsky, Thomas Bewley UCSD Coordinated Robotics Lab In this document, a Mobile Inverted Pendulum (MIP) is a robotic
More informationTechnische Universität Berlin
Technische Universität Berlin Institut für Mathematik M7  A Skateboard(v1.) Andreas Steinbrecher Preprint 216/8 PreprintReihe des Instituts für Mathematik Technische Universität Berlin http://www.math.tuberlin.de/preprints
More informationOn the Nature of Random System Matrices in Structural Dynamics
On the Nature of Random System Matrices in Structural Dynamics S. ADHIKARI AND R. S. LANGLEY Cambridge University Engineering Department Cambridge, U.K. Nature of Random System Matrices p.1/20 Outline
More informationMechatronics 1: ME 392Q6 & 348C 31Aug07 M.D. Bryant. Analogous Systems. e(t) Se: e. ef = p/i. q = p /I, p = " q C " R p I + e(t)
V +  K R +   k b V R V L L J + V C M B Analogous Systems i = q. + ω = θ. C . λ/l = q v = x F T. Se: e e(t) e = p/i R: R 1 I: I e C = q/c C = dq/dt e I = dp/dt Identical dierential equations & bond
More informationVideo 1.1 Vijay Kumar and Ani Hsieh
Video 1.1 Vijay Kumar and Ani Hsieh 1 Robotics: Dynamics and Control Vijay Kumar and Ani Hsieh University of Pennsylvania 2 Why? Robots live in a physical world The physical world is governed by the laws
More informationarxiv: v1 [quantph] 4 Nov 2008
International Journal of Modern Physics E c World Scientific Publishing Company arxiv:0811.0496v1 [quantph] 4 Nov 2008 Extended HamiltonLagrange formalism and its application to Feynman s path integral
More informationSolutions to chapter 4 problems
Chapter 9 Solutions to chapter 4 problems Solution to Exercise 47 For example, the x component of the angular momentum is defined as ˆL x ŷˆp z ẑ ˆp y The position and momentum observables are Hermitian;
More informationDetermining the Normal Modes of Vibration
Determining the ormal Modes of Vibration Introduction at the end of last lecture you determined the symmetry and activity of the vibrational modes of ammonia Γ vib 3 ) = A 1 IR, pol) + EIR,depol) the vibrational
More informationIntroduction to Mechanical Vibration
2103433 Introduction to Mechanical Vibration Nopdanai Ajavakom (NAV) 1 Course Topics Introduction to Vibration What is vibration? Basic concepts of vibration Modeling Linearization SingleDegreeofFreedom
More informationSheet 07.5: Matrices II: Inverse, Basis Transformation
Fakultät für Physik R: Rechenmethoden für Physiker, WiSe 015/16 Dozent: Jan von Delft Übungen: Benedikt Bruognolo, Dennis Schimmel, Frauke Schwarz, Lukas Weidinger http://homepages.physik.unimuenchen.de/~vondelft/lehre/15r/
More informationIntroductory Physics. Week 2015/06/26
2015/06/26 Part I Multiparticle system multiparticle system So far, we have been limiting ourselfs to the problems of the motion of single particle (except a brief introduction of coupled oscillation).
More information16.1. PROBLEM SET I 197
6.. PROBLEM SET I 97 Answers: Problem set I. a In one dimension, the current operator is specified by ĵ = m ψ ˆpψ + ψˆpψ. Applied to the left hand side of the system outside the region of the potential,
More information7.3 Heat capacities: extensive state variables (Hiroshi Matsuoka)
7.3 Heat capacities: extensive state variables (Hiroshi Matsuoka) 1 Specific heats and molar heat capacities Heat capacity for 1 g of substance is called specific heat and is useful for practical applications.
More informationChem 467 Supplement to Lecture 19 Hydrogen Atom, Atomic Orbitals
Chem 467 Supplement to Lecture 19 Hydrogen Atom, Atomic Orbitals PreQuantum Atomic Structure The existence of atoms and molecules had long been theorized, but never rigorously proven until the late 19
More informationFormation flying in elliptic orbits with the J 2 perturbation
Research in Astron. Astrophys. 2012 Vol. 12 No. 11, 1563 1575 http://www.raajournal.org http://www.iop.org/journals/raa Research in Astronomy and Astrophysics Formation flying in elliptic orbits with
More information1 Multiplicity of the ideal gas
Reading assignment. Schroeder, section.6. 1 Multiplicity of the ideal gas Our evaluation of the numbers of microstates corresponding to each macrostate of the twostate paramagnet and the Einstein model
More informationAdvanced Robotic Manipulation
Advanced Robotic Manipulation Handout CS37A (Spring 017 Solution Set # Problem 1  Redundant robot control The goal of this problem is to familiarize you with the control of a robot that is redundant with
More informationLittle Orthogonality Theorem (LOT)
Little Orthogonality Theorem (LOT) Take diagonal elements of D matrices in RG * D R D R i j G ij mi N * D R D R N i j G G ij ij RG mi mi ( ) By definition, D j j j R TrD R ( R). Sum GOT over β: * * ( )
More informationLecture 4: Affine Transformations. for Satan himself is transformed into an angel of light. 2 Corinthians 11:14
Lecture 4: Affine Transformations for Satan himself is transformed into an angel of light. 2 Corinthians 11:14 1. Transformations Transformations are the lifeblood of geometry. Euclidean geometry is based
More informationAP PHYSICS SUMMER WORK
Name AP PHYSICS SUMMER WORK For both AP Physics I and II AP Physics 1 Curriculum  Kinematics  Linear Dynamics (Newton s Laws, momentum, oscillations)  Rotational Dynamics  Energy, Work, and Power 
More information