Nov Julien Michel
|
|
- Laurence Greene
- 5 years ago
- Views:
Transcription
1 Int roduct ion t o Free Ene rgy Calculat ions ov Julien Michel
2 Predict ing Binding Free Energies A A Recept or B B Recept or Can we predict if A or B will bind? Can we predict the stronger binder? Can we do this reliably and quickly?
3 W hat is in t he t oolbox? Rigorous Free Energy M e t hods Accuracy Em pirical Scoring Funct ions Approxim at e Free Energy M e t hods seconds / com pound CPU t im e days / com pound
4 The im port ance of t he ent ropy Better Lennard Jones interactions with the protein But configurat ionaly rest rict ed Experim ental Binding Free Energy difference : kcal.m ol -1 Free Energy calculations capture aspects of ligand binding that are overlooked (or crudely considered) by em pirical scoring funct ions : Ent ropic effect s Solvat ion
5 The concept of Phase space 1 dim ensional harm onic oscillator U ( x) U ( x,y) y x x U(x) = k 1 (x x o ) 2 2 dim ensional harm onic oscillator x y U(x,y) = k 1 (x x o ) 2 + k 2 (y y 0 ) 2 x Serious sim ulations typically requires thousands of degrees of freedom and hence the phase space is in thousands of dim ensions Mom otina I love you!
6 St at ist ical Physics in one equat ion P ens = P( r exp ) exp ( βu ( r )) ( ) βu ( r ) dr dr P P ens = ( r ) π ( r ) dr Boltzm ann dist ribut ion Ensem ble Average of P Value of P for one configuration r P = P (Post ulat e) macroscopic ens Point in Phase space
7 um erical int egrat ion : Quadrat ure ( ) Two atom s on a 7*7 checkerboard : 49*49 = 2401 different states contributes to eq (2) Lot s of states! In eq (2), the integral is in 3 dim ensions where is the num ber of atom s sim ulated. And there are thousands of atom s to sim ulate Direct Integration is not feasible Im portant rem ark : Many, m any of the possible states involve overlaps of atom s and hence are highly unlikely
8 um erical e int e grat ion : M ont e Carlo Generate random ly the position of the atom s on the checkerboard, calculate the property of interest I s (X i ) Repeat tim es The quantity I est is an estim ate of the integrant I I est = 1 i =1 ( X The essential advantage over quadrature techniques is that the m ethod can be applied to system s of high dim ensionality without having to build a m esh that covers space I S i )
9 W hat is t he volum e of a sphere? k 2 V π = V R k Γ k um ber of dim ensions k k V/V R 1.00E E E E E E E-69 As k increases, the vast m ajority of the points in the k-dim ension space lies outside of the sphere Random selection of points doesn t work Strong analogy with system s in the condensed phase. There are few low energy states that contributes m eaningfully to the integral and m any high energy states (e.g atom ic overlaps) that do not contribute.
10 Im port ance Sam pling Instead of drawing random points from an uniform distribution, draw points from a distribution π. The Monte Carlo integration equation becom es I est = 1 I ( X i i = 1 from π π ( X i ) ) π is selected so that points are in the region of space which contributes the m ost to the integrand (e.g, in the sphere) The bias on the selection of X i is rem oved when the contribution to the integrand is estim ated.
11 Im port ance Sam pling : Exam ple f ( x) = 3x I = f ( x) π π ( x) k k ( x) dx π π 0 1 = = 1 4x Estim ate of I after draw ing 100 sam ples w ith tw o different im portance sam pler Y Funct ion Average St d. Dev 0.5 π π X f(x) =3x**2 importance distribution uniform distribution
12 Im port ance Sam pling in St at ist ical Physics P ens = P( r exp ( βu ( r )) ( βu ( r )) dr ) exp If P(r ) does not dom inate the product in the num erator, then an ideal im portance sam pling function to estim ate equation (2) is : dr π ( r ) i = exp exp ( βu ( r )) ( βu ( r ) ) dr i Problem : Im possible to draw sam ples from π(r i ) without knowing the denom inator, which we can t (as it involves solving directly a very difficult int egral)
13 M arkov Chains ( in a nut shell) A Markov Chain is a set of probabilistic rules which governs transitions between states and is often represented as a transition m atrix Π p11 p12 p13 Probability of m oving from state 1 Π = p21 p31 p22 p32 p23 p33 to 3 Assum ing Π obeys a num ber of m athem atical properties, then the following is true π (1) ( n) = lim n p Π p (1) represent an arbitrary initial distribution (e.g, a random s starting point) and Π (n) represents n applications of the transition m atrix Π This equations m eans : Repeated applications of the transition m atrix Πconverges any initial arbitrary distribution p (1) towards a unique lim iting distribution π If the transition m atrix is properly constructed, we can draw sam ples from p without having to specify π a priori ( e.g, not have to solve the bad integral )
14 M arkov Chains M ont e Carlo : M et ropolis M ont e Carlo ( ) 1. Start in state I 2. Attem pt a m ove to state j with probabiliy pij 3. Accept this m ove with probability α ij π j exp( βu j ) / Z = = = exp( β ( U π exp( βu ) / Z i i 6. If the m ove is accepted, set i = j, otherwise i = i 7. Accum ulate any property of interest A(i) 8. Return to 1 or term inate after iterations j U i )) Perhaps the m ost im portant 20th century contribution to num erical science
15 ( U ( r )) π ( r dr A = exp + β ) Absolut e free ene rgy calculat ion : P P ens = ( r ) π ( r ) dr ( ) U ( r ) π ( r dr A = exp + β ) A = exp( + βu ( r ) The first term in the integral is significant com pared to the second term. Regions of space that contributes to the integral are not always those that have a significant value of π (r) Very slow convergence of t he free energy Anyway What is the free energy of a cactus?
16 * probability would be m ore accurate Relat ive free energy calculat ion : Zwanzig equat ion (1954) 1 AA B = ln exp β β ( ( U ) B ( r ) U A( r ) A If A and B are sim ilar, then m ost configurations of the system A will have a sim ilar energy to the system B and the energy difference is 0 Much better convergence because only a sm all num ber of configurations have to be generated (those that differs in energy* between A and B) Energy Ua Ub (Ub-Ua) Phase Space
17 The coupling param et er λ : A and B have to be very very sim ilar for the previous equation to yield converged results in a reasonable tim e Solution, m ulti stage the calculation with as m any interm ediates potential as it takes A B A A = A = 1 k= 0 1 β ln exp ( β ( U ( r ) U ( r ))) λ k + 1 λ k λ k The interm ediates potentials U k can be anything you want. Typically taken by linear com bination of the potential for system A and B U k λ = (1 λk ) U A + λ k U B
18 Therm odynam ic Int egrat ion : A = λ 1 A B = λ = 0 U λ λ dλ Another way to get the free energy difference between A and B The quantity in brackets is the ensem ble average of the gradient of the energy with respect to the coupling param eter Calculated by im plem enting λ derivatives in the sim ulation code Or by finite difference approxim ation U U λ λ If the change A to B is big, the gradients varies a lot and m any points are necessary to obtain a good estim ate of A by num erical integration
19 Sum m ary G? 0.0 λ 1.0
20 Sum m ary H λ G = 1 = λ= 0 λ λ dλ λ = 0. 0 λ = 0. 5 λ = 1. 0 G
21 Relat ive Binding Free e Energy Calculat ions by a The rm odynam ic cycle A in water A in protein G 3 G 1 G 2 G 4 B in water B in protein G = G 4 - G 3 = G 2 - G 1
22 Problem s and pit falls : Forcefield The potentials U A and U B are approxim ate and this affect the quality of the results (ligand polarisation in water and binding site?) Convergence Long sim ulation tim es required to get a good estim ate of the free energy change (= m any hours of CPUs) Many interm ediate potentials have to be run to get a god estim ate of the free energy change (= m any CPUs) It is im possible to know if your results have converged! Som e calculations give the sam e answer if you start from a different point in phase space and run with a different random num ber. Other never gives the sam e answer This restrict the applicability of the m ethod to system that are very sim ilar.
23 Diagnosing convergence: Block averaging Cut the raw data into chunks and analyse independently each chunk to get a free energy. Measure the variance of the result ing dist ribut ion of free energies. Independent runs Run the com plete sim ulation m any tim es, from different starting points Therm odynam ic cycle closure A G A B C B G A C G C B G A B - G C B G A C = 0 Deviations are a m easure of lack of convergence one of the m ethods ensure your results are trully converged
24 W hy is it so difficult? U Phase space In docking, we only care about the single lowest energy point in phase space and the search algorithm s can cut any corner to get there But in statistical physics, you need to know the probability of each low energy point in phase Very difficult to explore quickly phase space and sim ultaneously obtain a good estim ate of the probability of the low energy states
25 Current scope of free e energy calculat ions: Biot in/ Sprept avidin COX-2. COX-1 Trypsin HIV-1 prot ease P38 ( ) We need Better algorithm s to run m ore quickly Typical error in a successful free energy sim ulation = 1 kcal.m ol -1 (e.g., a 5-6 fold difference in IC 50 at room tem perature) Better forcefields to get m ore accurate results Better sam pling m ethods to get m ore precise results Better ways to sim ulate different ligands to deal with a m ore diverse set of com pounds
26 Quest ions.
Introd uction to Num erical Analysis for Eng ineers
Introd uction to Num erical Analysis for Eng ineers System s of Linear Equations Cram er s Rule Gaussian Elim ination Num erical im plem entation Num erical stability: Partial Pivoting, Equilibration,
More informationProcedures for Computing Classification Consistency and Accuracy Indices with Multiple Categories
esearcli R eport Semes 2 0 0 0 Procedures for Computing Classification Consistency and Accuracy Indices with Multiple Categories Won-Chan Lee Bradley A. Hanson Robert L. Brennan ACT O c t o b e r For additional
More informationProtein- Ligand Interactions
Protein- Ligand Interactions Outline Introduction Interm olecular bonding forces Electrostatic or ionic bonds Hydrogen bonds Van der Waals interactions Dipole- dipole and ion- dipole interactions Repulsive
More informationOPTIMISATION PROCESSES IN TIDAL ANALYSIS
OPTIMISATION PROCESSES IN TIDAL ANALYSIS by M. T. M u r r a y Liverpool Tid al In stitu te Introduction An analysis of tidal observations is norm ally carried out, not as an end in itself, but in order
More informationFew thoughts on PFA, from the calorim etric point of view
1 Few thoughts on PFA, from the calorim etric point of view discussing pad siz e, M oliere radius, distance from the interaction point. H enri Videau 2 H ave you had anytim e a close view at show ers?
More informationС-4. Simulation of Smoke Particles Coagulation in the Exhaust System of Piston Engine
Kaliningrad 2012 Simulation of Smoke Particles Coagulation in the Exhaust System of Piston Engine С-4 Sergey M. Frolov1, Konstantin A. Avdeev1, Vladislav S. Ivanov1, Branislav Basara2, Peter Priesching2,
More informationAPPLICATION OF AUTOMATION IN THE STUDY AND PREDICTION OF TIDES AT THE FRENCH NAVAL HYDROGRAPHIC SERVICE
APPLICATION OF AUTOMATION IN THE STUDY AND PREDICTION OF TIDES AT THE FRENCH NAVAL HYDROGRAPHIC SERVICE by L. R o u m é g o u x Ingénieur H ydrographe en Chef, C hief of the T idal Section The considerable
More informationGENETIC ALGORlIHMS WITH AN APPLICATION TO NONLINEAR TRANSPORTATION PROBLEMS
101 GENETIC ALGORlIHMS WITH AN APPLICATION TO NONLINEAR TRANSPORTATION PROBLEMS MATTHEW HOBBS ISOR, VICTORIA UNIVERSITY WELLINGTON Abstract In this paper I will describe the theory and implementation of
More informationTIDAL PREDICTION WITH A SMALL PERSONAL COMPUTER
International Hydrographic Review, Monaco, LXII (2), July 1985 TIDAL PREDICTION WITH A SMALL PERSONAL COMPUTER by A.S. F R A N C O 1*1 ABSTRACT This paper show s how a personal com puter with 16Kb only
More informationLarge chunks. voids. Use of Shale in Highway Embankments
Use of Shale in Highway Embankments C. W. Lovell R esearch Engineer Joint Highway R esearch Project School of Civil E ngineering P urdue University IN T R O D U C T IO N Poor perform ance of m idw estern
More informationPresented by Arkajit Dey, Matthew Low, Efrem Rensi, Eric Prawira Tan, Jason Thorsen, Michael Vartanian, Weitao Wu.
Presented by Arkajit Dey, Matthew Low, Efrem Rensi, Eric Prawira Tan, Jason Thorsen, Michael Vartanian, Weitao Wu. I ntroduction Transient Chaos Sim ulation and Anim ation Return Map I Return Map I I Modified
More informationChapter 5 Workshop on Fitting of Linear Data
Chapter 5 Workshop on Fitting of Linear Data (Contributed by E.L. Allen, SJSU) 5.0 Learning Objectives After successfully com pleting this laboratory workshop, including the assigned reading, the lab repot
More informationEC 219 SA M PLING AND INFERENCE. One and One HalfH ours (1 1 2 H ours) Answerallparts ofquestion 1,and ONE other question.
ECONOM ICS EC 219 SA M PLING AND INFERENCE One and One HalfH ours (1 1 2 H ours) Answerallparts ofquestion 1,and ONE other question. M athem aticaland Statisticaltables are provided. An approved calculatormay
More informationReal Gas Equation of State for Methane
Kaliningrad 2012 1. Commission Regulation (EC) No 692/2008 of 18 July 2008. http://ec.europa.eu/enterprise/sectors/automotive/environment/eurovi/index_ en. htm, p. 130. 2. Payne J. F. B., Skyrme G. Int.
More informationAuthor's personal copy
Physics Letters A 373 (2009) 2717 2721 Contents lists available at S c ienc ed irec t Physics Letters A w w w.elsevier.c om /loc ate/p la A search for the sim plest chaotic partial differential equation
More informationMAST ACADEMY OUTREACH. WOW (Weather on Wheels)
MAST ACADEMY OUTREACH ELEMENTARY PROGRAM Adventures Aboard WOW (Weather on Wheels) Post-Site Packet MAST Academy Maritime and Science Technology High School Miami-Dade County Public Schools Miami, Florida
More informationThe M echanism of Factor VIII Inactivation by H um an Antibodies
ANNALS O F CLINICAL A N D LABORATORY SC IE N C E, Vol. 15, No. 1 Copyright 1985, Institute for Clinical Science, Inc. The M echanism of Factor VIII Inactivation by H um an Antibodies II. The Effect of
More informationA Study of Drink Driving and Recidivism in the State of Victoria Australia, during the Fiscal Years 1992/ /96 (Inclusive)
A Study of Drink Driving and Recidivism in the State of Victia Australia, during the Fiscal Years 1992/93-1995/96 (Inclusive) M ichael M o l o n e y, D iane P a l a ia V ictia Police Fce, T raffic A lcohol
More information3B/16. L a force portante des fondations en coins
3B16 T h e U ltim a te B e a r in g C a p a c ity o f W e d g e -sh a p e d F o u n d a t io n s L a force portante des fondations en coins by P ro fe sso r G. G. M e y e r h o f, D. S c., P h. D., F.
More informationC o r p o r a t e l i f e i n A n c i e n t I n d i a e x p r e s s e d i t s e l f
C H A P T E R I G E N E S I S A N D GROWTH OF G U IL D S C o r p o r a t e l i f e i n A n c i e n t I n d i a e x p r e s s e d i t s e l f i n a v a r i e t y o f f o r m s - s o c i a l, r e l i g i
More information(2009) Journal of Rem ote Sensing (, 2006) 2. 1 (, 1999), : ( : 2007CB714402) ;
100724619 (2009) 0220183207 Journal of Rem ote Sensing 1, 2, 3, 3, 3, 1 1., 100101; 2., 100049; 3., 100080 :,,, 3, ( ),1%,, :,,, : TP79 : A 1 20,,,,;, (, 1999),,,,, (, 2006),,,,, 2 2. 1 : 2007209217; :
More informationNuclear Science Research Group
UR-NCUR 05-13 THE UNIVERSITY OF ROCHESTER Nuclear Science Research Group ROCHESTER, NEW YORK 1467-016, USA Surface Matters - Prospects for Mechanistic Studies at Intermediate Energies W. Udo Schröder Nancy
More informationArchit ect ure s w it h
Inst ruct ion Select ion for Com pilers t hat Target Archit ect ure s w it h Echo Inst ruct ions Philip Brisk Ani Nahapetian Majid philip@cs.ucla.edu ani@cs.ucla.edu Sarrafzadeh m ajid@cs.ucla.edu Em bedded
More informationThe Ability C ongress held at the Shoreham Hotel Decem ber 29 to 31, was a reco rd breaker for winter C ongresses.
The Ability C ongress held at the Shoreham Hotel Decem ber 29 to 31, was a reco rd breaker for winter C ongresses. Attended by m ore than 3 00 people, all seem ed delighted, with the lectu res and sem
More informationStraight Lines. Distance Formula. Gradients. positive direction. Equation of a Straight Line. Medians. hsn.uk.net
Distance Formula Straight Lines Distance ( ) + ( ) between points (, ) and (, ) Gradients m between the points (, ) and (, ) where Positive gradients, negative gradients, zero gradients, undefined gradients
More informationOn the convergence of solutions of the non-linear differential equation
MEMOIRS O F T H E COLLEGE O F SCIENCE, UNIVERSITY OF KYOTO, SERIES A Vol. XXVIII, Mathematics No. 2, 1953. On the convergence of solutions of the non-linear differential equation By Taro YOSHIZAWA (Received
More informationAlles Taylor & Duke, LLC Bob Wright, PE RECORD DRAWINGS. CPOW Mini-Ed Conf er ence Mar ch 27, 2015
RECORD DRAWINGS CPOW Mini-Ed Conf er ence Mar ch 27, 2015 NOMENCLATURE: Record Draw ings?????? What Hap p ened t o As- Built s?? PURPOSE: Fur n ish a Reco r d o f Co m p o n en t s Allo w Locat io n o
More informationImpact of Drink-drive Enforcement and Public Education Programs in Victoria, Australia
Impact of Drink-drive Enforcement and Public Education Programs in Victoria, Australia D avid H ealy, T ransport A ccident C om m ission B A C K G R O U N D In D ecem ber 1989, the T ransport A ccident
More informationThe Ind ian Mynah b ird is no t fro m Vanuat u. It w as b ro ug ht here fro m overseas and is now causing lo t s o f p ro b lem s.
The Ind ian Mynah b ird is no t fro m Vanuat u. It w as b ro ug ht here fro m overseas and is now causing lo t s o f p ro b lem s. Mynah b ird s p ush out nat ive b ird s, com p et ing for food and p laces
More informationBOGOLIUBOV LABORATORY OF THEORETICAL PHYSICS. Theoretical Physics at BLTP. Precision tests of the Standard Model at LHC
BOGOLIUBOV LABORATORY OF THEORETICAL PHYSICS Theoretical Physics at BLTP. Precision tests of the Standard Model at LHC Dubna, Russia, July 1, 2008 Nikola i N ik ola e vich Bog oliubov ( 1 9 0 9 1 9 9 2
More informationThe Effectiveness of the «Checkpoint Tennessee» Program
The Effectiveness of the «Checkpoint Tennessee» Program John H. Lacey*, R alph K. Jones*and Jam es C. Fell** *M id-a m erica R esearch Institute, Shepherdstow n, W est V irginia U SA **N ational H ighw
More informationLight-absorbing capacity of phytoplankton in the Gulf of Gdańsk in May, 1987
Light-absorbing capacity of phytoplankton in the Gulf of Gdańsk in May, 1987 O C E A N O L O G IA, 28, 1990 P L ISSN 0078-3234 P h y to p lan k to n C hlorophyll L ight-absorbing capacity Pigm ent index
More informationChapter - 3. ANN Approach for Efficient Computation of Logarithm and Antilogarithm of Decimal Numbers
Chapter - 3 ANN Approach for Efficient Computation of Logarithm and Antilogarithm of Decimal Numbers Chapter - 3 ANN Approach for Efficient Computation of Logarithm and Antilogarithm of Decimal Numbers
More informationOn the M in imum Spann ing Tree Determ ined by n Poin ts in the Un it Square
Vol. 14 No. CH INESE QUARTERLY JOURNAL OF M ATHEM ATICS June 1999 On the M in imum Spann ing Tree Determ ined by n Poin ts in the Un it Square Ye J ichang ( ) Xu Yinfeng ( ) Xu Chengxian ( ) (X iπan J
More informationPDF hosted at the Radboud Repository of the Radboud University Nijmegen
PDF hosted at the Radboud Repository of the Radboud University Nijmegen The following full text is a publisher's version. For additional information about this publication click this link. http://hdl.handle.net/2066/32618
More informationF O R M T H R E E K enya C ertificate of Secondary E ducation
N a m e : A d m. N o...... D a t e : C la ss:.... 565/1 FO R M 3 B U S I N E S S S T U D I E S P A P E R 1 T I M E : 2 H O U R S T R I A L 6 2 0 1 8 FO R M 3 B U S I N E S S S T U D I E S P A P E R 1 T
More informationMark Schem e (Results) March GCSE Mathem atics (Linear) 1MA0 Higher (Calculator) Paper 2H
Mark Schem e (Results) March 2013 GCSE Mathem atics (Linear) 1MA0 Higher (Calculator) Paper 2H Edexcel and BTEC Qualifications Edexcel and BTEC qualifications com e from Pearson, the world s leading learning
More informationInterpretation of Laboratory Results Using M ultidimensional Scaling and Principal C om ponent Analysis*
ANNALS OF CLINICAL AND LABORATORY SCIENCE, Vol. 17, No. 6 Copyright 1987, Institute for Clinical Science, Inc. Interpretation of Laboratory Results Using M ultidimensional Scaling and Principal C om ponent
More informationChapter 5. Self-condensation of acetophenone
Chapter 5 Self-condensation of acetophenone 5.1. Introduction D ypnone is a useful intermediate for the production o f a large range o f compounds. It has been used as a softening agent, plastisizer and
More information(IGBP) km 2,? PRO GR ESS IN GEO GRA PH Y. V o l. 20, N o. 4 D ec., 2001 : (2001) m m 3,
20 4 2001 12 PRO GR ESS IN GEO GRA PH Y V o l. 20, N o. 4 D ec., 2001 : 100726301 (2001) 0420355209, (, 100101) :, 501 mm 502 mm, 765 10 8 m 3 638 10 8 m 3, 127 10 8 m 3,,, ;, : ; ; : P42616; P343 : A
More informationM A T R IX M U L T IP L IC A T IO N A S A T E C H N IQ U E O F P O P U L A T IO N AN A L Y S IS
M A T R IX M U L T IP L IC A T IO N A S A T E C H N IQ U E O F P O P U L A T IO N AN A L Y S IS N ATH A N KEYFITZ This paper takes advantage o f the fact that the changes in numbers o f a human (o r any
More informationOscar Cubo M edina fi.upm.es)
. euristics. reedy. LocalSearch. illimbing. Meta-heuristics. Tabu Search. Sim ulated nnealing Oscar ubo M edina (ocubo@ fi.upm.es) ased on Neighbourhood Search euristic derives from the verb heuriskein
More informationEffect of Methods of Platelet Resuspension on Stored Platelets
ANNALS O F CLINICAL AND LABORATORY S C IEN C E, Vol. 14, No. 5 Copyright 1984, Institute for Clinical Science, Inc. Effect of Methods of Platelet Resuspension on Stored Platelets THOMAS KIRALY, M.A., S.B.B.
More informationHYBRIDISATION OF A WAVE PROPAGATION MODEL (SWASH) AND A MESHFREE PARTICLE METHOD (SPH) FOR REAL APPLICATIONS
e 3rd IAHR Europe Congress, Book of Proceedings, 2014, Porto - Portugal. ISBN 978-989-96479-2-3 HYBRIDISATION OF A WAVE PROPAGATION MODEL (SWASH) AND A MESHFREE PARTICLE METHOD (SPH) FOR REAL APPLICATIONS
More information(Received March 30, 1953)
MEMOIRS OF THE COLLEGE OF SCIENCE, UNIVERSITY OF KYOTO, SERIES, A V ol. XXVIII, Mathematics No. 1, 1953. On the Evaluation of the Derivatives of Solutions of y"=-f (x, y, y'). By Taro YOSHIZAWA (Received
More informationTo link to this article:
This article was downloaded by: [Duke University Libraries] On: 21 May 2015, At: 09:26 Publisher: Routledge Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer
More informationDesigning the Human Machine Interface of Innovative Emergency Handling Systems in Cars
Designing the Human Machine Interface of Innovative Emergency Handling Systems in Cars M anfred D A N G E L M A IE R * and P etra M A N G O L ** *Fraunhofer IAO, N obelstr. 12, D -70569 Stuttgart, **IA
More informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
More informationIntroduction to Machine Learning CMU-10701
Introduction to Machine Learning CMU-10701 Markov Chain Monte Carlo Methods Barnabás Póczos & Aarti Singh Contents Markov Chain Monte Carlo Methods Goal & Motivation Sampling Rejection Importance Markov
More informationThe S tructure of L iquid H ydrogen Peroxide
The S tructure of L iquid H ydrogen Peroxide By J. T. R a n d a ll, M.Sc. Communication from the Staff of the Research Laboratories of The General Electric Company Limited, Wembley, England ( Communicated
More informationECE 302 Division 1 MWF 10:30-11:20 (Prof. Pollak) Final Exam Solutions, 5/3/2004. Please read the instructions carefully before proceeding.
NAME: ECE 302 Division MWF 0:30-:20 (Prof. Pollak) Final Exam Solutions, 5/3/2004. Please read the instructions carefully before proceeding. If you are not in Prof. Pollak s section, you may not take this
More informationJia, Yudan (2011) Numerical modelling of shaft lining stability. PhD thesis, University of Nottingham.
Jia, Yudan (2011) Numerical modelling of shaft lining stability. PhD thesis, University of Nottingham. Access from the University of Nottingham repository: http://eprints.nottingham.ac.uk/11924/1/phd -_Y._Jia.pdf
More informationINERTIAL SURVEY SYSTEM (ISS)
International Hydrographie Review, M onaco, LV III (2), July 1981 INERTIAL SURVEY SYSTEM (ISS) by Lt. Cdr. Lew is A. L A PIN E, NOAA P aper presented at the NOS H ydrographie Survey C onference, Jan u
More informationBoris Backovic, B.Eng.
SOURCE-CHANNEL CODEC FOR A WCDMA BASED MULTIMEDIA SYSTEM by Boris Backovic, B.Eng. Toronto, September 2005 A Project Report presented to Ryerson University in partial fulfillment of the requirements for
More informationApplying Metrics to Rule-Based Systems
Applying Metrics to Rule-Based Systems A Thesis by: Supervisor: Paul D oyle, B.Sc. M r. Renaat Verbruggen Subm itted to D ublin C ity U niversity C om puter A pplications for the degree o f M a s te r
More informationLectures 5 & 6: Hypothesis Testing
Lectures 5 & 6: Hypothesis Testing in which you learn to apply the concept of statistical significance to OLS estimates, learn the concept of t values, how to use them in regression work and come across
More informationCHEM-UA 652: Thermodynamics and Kinetics
1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 4 I. THE ISOTHERMAL-ISOBARIC ENSEMBLE The isothermal-isobaric ensemble is the closest mimic to the conditions under which most experiments are
More informationUNIT 1: CHEMISTRY AND SOCIETY
Honors Chemistry Chemistry offers a curriculum that emphasizes students' understanding of fundamental chemistry concepts while helping them acquire tools to be conversant in a society highly influenced
More informationSome Filtering Techniques for Digital Image Processing
Some Filtering Techniques for Digital Image Processing by Chi-M ing Leung l.l.sc. McMaster University, 1975 M. A. University of British Columbia, 1971 B.Sc. T h e Chinese U niversity of Hcng Kong, 1969
More informationStrategies for E stim ating B ehaviou ral F requency in Survey Interview s
Ó M E M O R Y, 19 9 8, 6 (4 ), 33 9 ±3 66 Strategies for E stim ating B ehaviou ral F requency in Survey Interview s F re derick G. C onrad U S B ureau of Labor Statistics N orm an R. B row n U niversity
More informationPrim ary I onization Track (Gases)
Prim ary I onization Track (Gases) incoming particle ionization track Minimum- ionizing particles ( Sauli. IEEE+NSS ) ion/e - pairs Helium Argon Xenon DME GAS (STP) CH 4 de/dx (kev /cm ) n (ion pairs/
More informationJP3.7 SHORT-RANGE ENSEMBLE PRECIPITATION FORECASTS FOR NWS ADVANCED HYDROLOGIC PREDICTION SERVICES (AHPS): PARAMETER ESTIMATION ISSUES
JP3.7 SHORT-RANGE ENSEMBLE PRECIPITATION FORECASTS FOR NWS ADVANCED HYDROLOGIC PREDICTION SERVICES (AHPS): PARAMETER ESTIMATION ISSUES John Schaake*, Mary Mullusky, Edwin Welles and Limin Wu Hydrology
More informationM ultiple O rg anisations of E ven ts in M em ory
Ó M E M O R Y, 1 99 7, 5 (5 ), 5 69 ±59 9 M ultiple O rg anisations of E ven ts in M em ory Juliana S. L a nc aste r M orris Brow n C ollege, Atlanta, U SA L a w rence W. B arsa lou E m ory U niversity,
More informationHYDROGRAPHIC SOUNDINGS IN SURF AREAS
HYDROGRAPHIC SOUNDINGS IN SURF AREAS b y Captain H enri D u f f o (Ret.), French N a v y In the course of a recent survey on the W est Coast of A frica w h ich the Laboratoire Central d H y draulique de
More informationO p tim ization o f P iezo electric A ctu a to r C onfigu ration on a F lexib le F in for V ib ration C ontrol U sin g G en etic A lgorith m s
O p tim ization o f P iezo electric A ctu a to r C onfigu ration on a F lexib le F in for V ib ration C ontrol U sin g G en etic A lgorith m s by A ndrew R ader, B.E ng. A Thesis subm itted to the Faculty
More informationMonte Carlo integration (naive Monte Carlo)
Monte Carlo integration (naive Monte Carlo) Example: Calculate π by MC integration [xpi] 1/21 INTEGER n total # of points INTEGER i INTEGER nu # of points in a circle REAL x,y coordinates of a point in
More informationA GRAPHICAL METHOD OF ADJUSTING A DOUBLY-BRACED QUADRILATERAL
A GRAPHICAL METHOD OF ADJUSTING A DOUBLY-BRACED QUADRILATERAL OF MEASURED LENGTHS WHEN THEIR WEIGHTS ARE ALSO TAKEN INTO CONSIDERATION b y J.C. B h a t t a c h a r j i, M.A., F.R.A.S., F.I.G.U. Survey
More informationThe effect of pollutants of crude oil origin on the diffusive reflectance of the ocean*
The effect of pollutants of crude oil origin on the diffusive reflectance of the ocean* O C E A N O L O G IA, 23, 1986 PL ISSN 0078-3234 R adiance Oil pollution Diffusive reflectance K r z y s z t o f
More informationPDF hosted at the Radboud Repository of the Radboud University Nijmegen
PDF hosted at the Radboud Repository of the Radboud University Nijmegen The following full text is an author's version which may differ from the publisher's version. For additional information about this
More informationThe use of small airplanes to gather swan data in Alaska
Wildfowl (1973) 24:15-20 The use of small airplanes to gather swan data in Alaska JAMES G. KING A laska is th e centre of abundance for N orth A m erica s nesting sw ans. T he W histling Swan Cygnus columbianus
More informationMark Schem e (Results) Novem ber GCSE Mathem atics (Linear) 1MA0 Higher (Non-Calculator) Paper 1H
Mark Schem e (Results) Novem ber 2012 GCSE Mathem atics (Linear) 1MA0 Higher (Non-Calculator) Paper 1H Edexcel and BTEC Qualifications Edexcel and BTEC qualifications com e from Pearson, the world s leading
More information2 Semester Final Exam Study Guide
Vista Murrieta High School Advanced Placem ent Biology G. Nicholas nd 2 Semester Final Exam Study Guide The following questions will potentially be on the final exam. Answer all questions in an outline
More informationAnalysis and Experimental Study of Vibration System Characteristics of Ultrasonic Compound Electrical Machining
UDC 539.4 Analysis and Experimental Study of Vibration System Characteristics of Ultrasonic Compound Electrical Machining Z. Q. D eng, Y. W. Z h u,1 F. W ang, X. G u, an d D. Y ang Collage of Mechanical
More informationTHERMAL MASS FLOW MEASUREMENT
THERMAL MASS FLOW MEASUREMENT by Gabriel J. Duffy, B.Eng. A thesis submitted to Dublin City University for the degree of Master of Engineering S c h o o l o f E l e c t r o n ic E n g in e e r in g Dublin
More informationFr anchi s ee appl i cat i on for m
Other Fr anchi s ee appl i cat i on for m Kindly fill in all the applicable information in the spaces provided and submit to us before the stipulated deadline. The information you provide will be held
More informationBayesian Methods with Monte Carlo Markov Chains II
Bayesian Methods with Monte Carlo Markov Chains II Henry Horng-Shing Lu Institute of Statistics National Chiao Tung University hslu@stat.nctu.edu.tw http://tigpbp.iis.sinica.edu.tw/courses.htm 1 Part 3
More informationSoil disturbance during continuous flight auger piling in sand
Soil disturbance during continuous flight auger piling in sand Perturbation du sol au cours de forage à la tarière en milieu sableux M.J.Kenny & K. Z. Andrawes - Department of Civil Engineering, University
More informationLSU Historical Dissertations and Theses
Louisiana State University LSU Digital Commons LSU Historical Dissertations and Theses Graduate School 1976 Infestation of Root Nodules of Soybean by Larvae of the Bean Leaf Beetle, Cerotoma Trifurcata
More informationT h e C S E T I P r o j e c t
T h e P r o j e c t T H E P R O J E C T T A B L E O F C O N T E N T S A r t i c l e P a g e C o m p r e h e n s i v e A s s es s m e n t o f t h e U F O / E T I P h e n o m e n o n M a y 1 9 9 1 1 E T
More informationSodium-Initiated Polymerization of Alpha- Methylstyrene in the Vicinity of Its Reported Ceiling Temperature
Western Michigan University ScholarWorks at WMU Dissertations Graduate College 8-1976 Sodium-Initiated Polymerization of Alpha- Methylstyrene in the Vicinity of Its Reported Ceiling Temperature Shuenn-long
More informationPractical Baraminology
CEN Tech. J., vol. 6(2), 1992, pp. 122 137 Practical Baraminology DR KURT P. WISE ABSTRACT Baram inology is the most efficient biosystematics theory and m ethodology available to the young-earth creationist.
More informationContr, Mo..?? M arin e B iologies! L ab o rato ry H singer D anm ark
O PH ELIA, 12: 117-127 (December 1973) Insfifuü! voer Z e e ^ r c ^ i p e lijk onderzoek kclh Lï i'.:... ; í ;:3ardi î i.n ses. us o... i í.!. í.n 69 841 B re d e n e - B elgium - 7 el. 59 / 8 37 15 SEASONAL
More informationT he poisoning of a palladium catalyst by carbon m onoxide
T he poisoning of a palladium catalyst by carbon m onoxide By M. G. T. Burrows and W. H. Stookmayer Leoline Jenkins Laboratories, Jesus College, Oxford (i Communicated by D. L. Chapman, F.R.8. Received
More informationPDF hosted at the Radboud Repository of the Radboud University Nijmegen
PDF hosted at the Radboud Repository of the Radboud University Nijmegen The following full text is a publisher's version. For additional information about this publication click this link. http://hdl.handle.net/2066/35623
More informationNORTHW ESTERN UNIVERSITY. A Case S tudy of Three Universities A DISSERTATION SUBMITTED TO THE GRADUATE SCHOOL
NORTHW ESTERN UNIVERSITY C onstructing T eaching Practices A round Novel Technologies: A Case S tudy of Three Universities A DISSERTATION SUBMITTED TO THE GRADUATE SCHOOL IN PARTIAL FULFILLM ENT O F THE
More informationM a rtin H. B r e e n, M.S., Q u i T. D a n g, M.S., J o se p h T. J a in g, B.S., G reta N. B o y d,
The effect of a one for the road drink of hard liquor, beer or wine on peak breath alcohol concentration in a social drinking environment with food consumption M a rtin H. B r e e n, M.S., Q u i T. D a
More informationMonte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)
More informationGeometric Predicates P r og r a m s need t o t es t r ela t ive p os it ions of p oint s b a s ed on t heir coor d ina t es. S im p le exa m p les ( i
Automatic Generation of SS tag ed Geometric PP red icates Aleksandar Nanevski, G u y B lello c h and R o b ert H arp er PSCICO project h ttp: / / w w w. cs. cm u. ed u / ~ ps ci co Geometric Predicates
More informationBasics of Statistical Mechanics
Basics of Statistical Mechanics Review of ensembles Microcanonical, canonical, Maxwell-Boltzmann Constant pressure, temperature, volume, Thermodynamic limit Ergodicity (see online notes also) Reading assignment:
More informationToday: Fundamentals of Monte Carlo
Today: Fundamentals of Monte Carlo What is Monte Carlo? Named at Los Alamos in 940 s after the casino. Any method which uses (pseudo)random numbers as an essential part of the algorithm. Stochastic - not
More informationH am iltonian Therm ostatting Techniques for. M olecular D ynam ics Sim ulation.
H am iltonian Therm ostatting Techniques for M olecular D ynam ics Sim ulation. Thesis subm itted for the degree of Doctor of Philosophy at the University of Leicester by C hristopher R ichard Sweet BSc
More informationBasic Semantics of Prolog. CAS 706 Program m ing Language. Xinjun Wu mcmaster.ca
Basic Semantics of Prolog CAS 706 Program m ing Language Xinjun Wu 0326680 wux8@ mcmaster.ca Basic Sem antics ofprolog Logic Program m ing Language 1 Preface Logical programming or constraint programming
More informationML estimation: Random-intercepts logistic model. and z
ML estimation: Random-intercepts logistic model log p ij 1 p = x ijβ + υ i with υ i N(0, συ) 2 ij Standardizing the random effect, θ i = υ i /σ υ, yields log p ij 1 p = x ij β + σ υθ i with θ i N(0, 1)
More informationBenchm ar k One St ud y Guid e: Science Benchm ar k Wed. Oct. 2nd
Benchm ar k One St ud y Guid e: Science Benchm ar k Wed. Oct. 2nd Char act er ist ics of Science: 1. What is a gr ad uat ed cylind er used t o m easur e? Volum e of a liq uid (Cap acit y) 2. Exp lain how
More informationfor Global Optimization with a Square-Root Cooling Schedule Faming Liang Simulated Stochastic Approximation Annealing for Global Optim
Simulated Stochastic Approximation Annealing for Global Optimization with a Square-Root Cooling Schedule Abstract Simulated annealing has been widely used in the solution of optimization problems. As known
More informationMATH 3510: PROBABILITY AND STATS July 1, 2011 FINAL EXAM
MATH 3510: PROBABILITY AND STATS July 1, 2011 FINAL EXAM YOUR NAME: KEY: Answers in blue Show all your work. Answers out of the blue and without any supporting work may receive no credit even if they are
More informationMultiple time step Monte Carlo
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 18 8 NOVEMBER 2002 Multiple time step Monte Carlo Balázs Hetényi a) Department of Chemistry, Princeton University, Princeton, NJ 08544 and Department of Chemistry
More informationEssentials on the Analysis of Randomized Algorithms
Essentials on the Analysis of Randomized Algorithms Dimitris Diochnos Feb 0, 2009 Abstract These notes were written with Monte Carlo algorithms primarily in mind. Topics covered are basic (discrete) random
More informationIntroduction to Design of Experiments
Introduction to Design of Experiments Jean-Marc Vincent and Arnaud Legrand Laboratory ID-IMAG MESCAL Project Universities of Grenoble {Jean-Marc.Vincent,Arnaud.Legrand}@imag.fr November 20, 2011 J.-M.
More information7.1 Coupling from the Past
Georgia Tech Fall 2006 Markov Chain Monte Carlo Methods Lecture 7: September 12, 2006 Coupling from the Past Eric Vigoda 7.1 Coupling from the Past 7.1.1 Introduction We saw in the last lecture how Markov
More information