Free Energy Simulation Methods
|
|
- Jewel Haynes
- 6 years ago
- Views:
Transcription
1 Free Energy Simulation Methods
2
3 Free energy simulation methods Many methods have been developed to compute (relative) free energies on the basis of statistical mechanics Free energy perturbation Thermodynamic integration Umbrella sampling Slow growth Fast growth L. Landau R. W. Zwanzig J. G. Kirkwood
4 Recap of statistical mechanics: Partition function, free energy, enthalpy, & entropy Partition function of canonical (NVT) ensemble Z = exp( E / kt) dpdx E = E + E kinetic potential In classical mechanics, potential energy is independent of kinetic energy
5 Free energy, enthalpy, & entropy Free energy (Helmholtz) Internal energy A = kt ln exp( E / kt ) dx <>: Ensemble average Entropy S = ( U A) T
6 Difficulty for absolute free energy simulation Absolute free energy requires converged integration on 3N dimensions Z = exp( E / kt) dx 3N k B T No experimental absolute free energies available Relative free energies are what really matter
7 How computational chemists compute free energy? 1. Minimize geometry 2. Frequency calculation 3N-6 internal vibrations 3. Ideal gas model for translation 4. Ideal rotor model for rotation Severe problems for macromolecules Harmonic approximation Multiple conformation states Solvated states are difficult to model
8 Free energy perturbation Free energy difference between two states State 1; energy E1 State 2; energy E2
9 Free energy perturbation Free energy difference between two states ΔA1 2 Zwanzig, R. W., J. Chem. Phys. 1954, 22:
10 Free energy perturbation by Functional Calculus A( E( X)) = kt ln exp( E/ kt) dx δa( E( X)) = AE ( ) δedx E Homework: Derive Free Energy Perturbation equation using functional calculus.
11 Free energy perturbation Symmetry of FEP equations Δ A = A A ( E ) = kt ln exp( E / kt ) ( E ) = kt ln exp( E / kt ) = ( A A ) = ΔA Forward and backward simulations to reduce the error [ ] Δ A = ΔA ΔA /2 kt = ln exp( ( E E )/ kt) + ln exp( ( E E )/ kt)
12 Free energy perturbation Accuracy of FEP ρ(δe 1 2 )
13 Free energy perturbation Accuracy of FEP Δ A = kt ln exp( βδe ) ρ ( ΔE ) dδe Δ A = ktln exp( βδe ) ρ ( ΔE ) dδe
14 Free energy perturbation Technical complications Phase-space overlapping between two states Energy
15 Free energy perturbation Implementation λ 9 ΔA = ΔA 1 2 λ λ i i i+ 1 λ 1 λ 2 λ 3 λ 0 In general, free energy differences between discrete points can be determined and summed to give the total free energy difference
16 Approximated free energy perturbation Δ A = kt ln exp( βδe ) ρ ( ΔE ) dδe ( βδe ) i ln ρ1( 1 2) 1 2 i= 0 i! Δ A = kt ΔE dδe ( βδe ) 2 = kt ln 1 + ( βδ E ) ρ ( ΔE ) dδe = kt + Δ βδe1 2 βδe1 2 ln 1 β E i = kt Ci ( β ) / i! i= 1
17 Approximated free energy perturbation i Δ A1 2 = kt Ci ( β ) / i! i= 1 C = ΔE C = ΔE ΔE C = ΔE ΔE = Δ 1 2 Δ C E E C C 5 =... When distribution of ΔE is gaussian ( ) Δ A = ΔE βσ ΔE
18 Thermodynamic integration A virtual path linking two states 2 1 Kirkwood, J. G., J. Chem. Phys. 1935, 3:
19
20 Thermodynamic integration Sometimes we have a real reaction path Δ A = 0 λ ' = λ ' 0 λ ' 0 A( λ) dλ λ E( λ) dλ λ λ the area under the curve corresponds to the free energy change This is often used in simulating free energy profile and activation free energies for enzyme-catalyzed reactions.
21 Free Energy Perturbation vs Thermodynamic Integration Δ A = ktln exp( βδe ) Δ A = E( λ) λ λ dλ 1. Requirement for gradient 2. Error analysis 3. Easiness for implementation 4. Can be combined together
22 Umbrella sampling for modeling reaction processes ClCH 3 + OH- HOCH 3 + Cl-
23 Umbrella sampling E(R) E (R) P(R) Torrie, G. M., and Valleau, J. P., J. Comput. Phys. 1977, 23:
24 Umbrella sampling Modified energy function E '( R) = E( R) + V( R) The potential of mean force (PMF) of R is defined as [ ] A( R) = ktln( exp E( R)/ kt δ ( r R) dx) = kt ln ρ( R) + C PMF of biased energy function is [ ] A'( R) = ktln( exp E'( R) / kt δ ( r R) dx) = AR ( ) ktln = AR ( ) + V( R) + C [ ER ( )/ kt] δ ( r R) dx exp[ ] = kt ln( exp V( R)/ kt δ ( r R) dx) [ ] exp V( R)/ kt δ ( r R) dx
25 Umbrella sampling For two neighboring windows A'( R) = A( R) + V ( R) + C A'( R) = A( R) + V ( R) + C 2 2 2
26
27
28 A' ( R) = A( R) + V ( R) + C 1 1 1
29
30 Weighted histogram analysis method ρ Minimizing overall sampling uncertainties from linear combination of samples of different windows S β Vi( R) fi b ( u) ( R) = e ρi ( R) ρ ( R) = c ( R) ρ ( R) ci ( R) ( u) ( ) i 0 i= 1 i i S i= 1 = 1 2 ( σ ρ0 ( R) ) c ( R) = i 0 S S 2 2 ( ) i i i i= 1 i= 1 u c ( R) σ ρ ( R) μ c ( R) 1 2 ( σ ρ0 ( R) ) ci ( R) 2 ( u) ( ) ( ) 0 = = 2ci R σ ρi R μ ( i ) 2 ( u) ( ) ( ) 1 ci R = σ ρ R μ ( ) ( R) ( ) ( R) ( R) ( R) ( R) 2 ( b) i ( b) i i NiΔ ( R) ( ) ( ) 2 ( u) 2β Vi R fi 2 ( b) σ ρi R = e σ ρi R g σ ρ = ρ = g N ( ) i β Vi R fi ρ0 R e iδ
31 Umbrella sampling Overlapping between different sampling windows
32 Slow growth Replace discrete sampling configuration points to continuous quasi-equilibrium process ΔA = - WORK (done by a reversible / quasi-static isothermal process) Example: gas in a piston B Δ A = PdV = dw A B A Let s drive the system from one state to another in a very slow (quasistatic) way ΔA 1 2 N E( λi ) δλ i λ i= 0 i
33 Slow growth Noise and hysteresis in slow growth simulations ΔA + ΔA
34 Fast growth and Jarzynski s equality ΔA limw For a quasi-static slow process 0 1 = τ 0 N For a non-quasi-static process ΔA W Thermodynamic second law W For a non-quasi-static process =ΔA
35 Fast growth method J. Am. Chem. Soc., 2005, 127:
36 Thermodynamic cycle Thermodynamic cycle can be constructed to avoid the direct simulation of computationally difficult processes. Example: which of the two ligands, L1 and L2, binds stronger to the enzyme E? Approach 1: direct simulation of the two binding processes L(aq) 1 + E(aq) L-E(aq) 1 ΔG1 L(aq) + E(aq) L-E(aq) ΔG ΔΔ G = ΔG ΔG 1 2 Approach 2: simulating other two processes as Mutate free L1 into free L2 in solution Mutate L1 into L2 within the enzyme s active site L (aq) L (aq) ΔG L-E(aq) L-E(aq) ΔG
37 Thermodynamic cycle L 1 (aq) + E(aq) ΔG 1 L 1 -E(aq) ΔG 3 ΔG 4 L 2 (aq) + E(aq) ΔG 2 L 2 -E(aq) ΔG 1 ΔG + ΔG 2 4 ΔG ΔΔG = ΔG 1 = ΔG 4 2 ΔG ΔG 4 ΔG 3 3 ΔG = 3 0 Thermodynamic cycle can be constructed to avoid the direct simulation of computationally difficult processes.
38 QM/MM Methods for Complex Reaction Processes
39 Quantum mechanics In computational chemistry, the most immediate task is to compute energy and gradient for any given molecular conformational states! Computational cost Hartree-Fock: N 4 MP2: N 5 MP4: N 6 Coupled Cluster: N 7 DFT: N 3 ~N 4 Challenge developing cheap methods for the simulation of macromolecules!
40 General MM force field forms Five basic terms in all force fields V= bond stretching + angle bending + torsions + + bonds angles torsions atom pairs atom pairs electrostatic interaction van der Waals interaction Covalent interactions Non-bonded interactions
41 General MM force field forms Non-bonded interactions (electrostatic and van der Waals interactions) are calculated only between atoms not involved in direct bonding (1-2 interaction), or angle connections (1-3 interaction)
42 Example case: CO vibration ΔE 1412*(b-1.13) 2 Think about perturbation method in quantum chemistry!
43 Example case: H 2 O vibration ΔE *(θ-105) 2
44 Electrostatic potential based atomic charges Molecular electrostatic potential (ESP) is well defined: ϕ () r = ϕ () r + ϕ () r QM nuc ele nuclei Z A ρ( r ') = dr ' r r r r' A A
45 Electrostatic potential based atomic charges Fitting molecular ESP for point charges ϕ () r = MM atom A qa r r A min w( r) ϕmm ( r) ϕqm ( r) dr 2 w(r): weighting functions
46 Development of MM force fields Two fundamental assumptions: Born-Oppenheimer approximation The fast-moving electrons are adjusting instantaneously to the positions of the (relatively) slow-moving nuclei. The contributions from the electronic degrees of freedom are thus implicitly determined by the positions of nuclei. Transferable functional groups in chemistry Chemically identical function groups (for example OH, -C=O, -NH 2, ) can be described by the same set of interaction terms. Only limited number of functional groups need to be considered.
47 Development of MM force fields Unique features of MM force fields Electron-less Transferability Concept of atom type Requires extensive/tedious parameterization process Cannot describe chemical reactions in general!!!
48 Why QM/MM? QM MM Accuracy Computational cost Long-MD Polarization Bond formation / breaking Long-range vdw High High No Included Yes No Medium to low low Yes Ignored/ approximated No Yes
49 General QM/MM scheme QM MM E = E + E + E QM QM / MM MM 1. QM is used to describe the site where reactions occur, including those atoms make important and direct interactions to atoms undergoing valence change in the reactions process. 2. MM is used to describe the rest of the system. Presumably atoms in these regions contribute to the reaction moieties through a static and classical electrostatic fashion.
50 A simple approach: ONIOM method Our own n-layered Integrated molecular Orbital + molecular mechanics Method + = E QM (1) E (1+ 2) MM = E (1) + E () 2 + E (1/ 2) MM MM MM Mechanical embedding model E / (1+ 2) QM MM = E (1) + E (2) + E QM MM MM (1/ 2) = E (1) + E (1+ 2) E (1) QM MM MM
51 General QM/MM scheme MM QM H + = H + H + H QM MM QM QM/ MM MM = H + H + H QM QM/ MM, ele QM / MM, nuc l QM/ MM, vdw QM/ MM,cov a MM + H + + H + H E = E + E + E QM QM / MM MM
52 Electronic structure at QM/MM boundary Dangling QM/MM bond
53 Electronic structure at QM/MM boundary Linked hydrogen MM QM Pros: easy to implement Cons: different ensembles due to the additional hydrogen atom
54 Electronic structure at QM/MM boundary Frozen local orbital / Local SCF A strictly localized bond orbital (SLBO) for the bond. The SLBO is assumed to be transferable. SLBO is excluded from the SCF optimization and does not mix with other orbitals. To compensate for the additional electron introduced with the doubly occupied SLBO, an extra charge of 1e is placed on M1, which interacts with all other MM charges. More charges can be placed on MM atoms around to improve.
55 Electronic structure at QM/MM boundary Pseudo-bond approach A monovalent, fluorine-like boundary atom with seven valence electrons, Z = 7, was used to replace the boundary atom. An angular-momentum dependent effective core potential (ECP) was used to mimic the original C-C bonds. ECP is therefore specific to the type of the bond on the boundary and is obtained by fitting. But fitted parameters show weak dependence on the QM method.
56 Electrostatic embedding model MM QM HQM+ MM = H + H + H QM QM/ MM MM = H + H + H QM QM/ MM, ele QM / MM, nucl + H + H + H ψ H ψ = ψ H QM/ MM, vdw QM/ MM,cova MM QM+ MM QM + H + E + E QM/ MM, ele QM/ MM, nucl QM/ MM, vd W + E + E QM/ MM,cova MM When MM atoms are represented as point charges H + H = H + QM QM / MM, ele QM ψ q i i MM r r i
57 QM/MM electrostatic interactions Point charge Point MM charge has singular point at distance 0 Over-polarization of QM electrons Gaussian charge or damping of the QM/MM electrostatic interactions will improve
58 QM/MM elec. interactions: lack of Pauli interactions Pauli repulsion ensures two nuclei won t collapse into each other Lack of Pauli repulsion leads to two effects: incorrect energetics incorrect electron densities
59 QM/MM electrostatic interaction Long range electrostatic interactions Do QM images see each other? Limits the size of the system and the QM basis sets Contribution of image charge distributions to the central QM region Ewald type of method Multiple grids
60 QM/MM vdw Simple Lennard-Jones 12-6 form is often used E C C = r r 12 6 QM / MM, vdw 12 6 But dispersions of QM atoms depend on the QM electron distribution which further depends on the surrounding MM point charges T.J. Giese & D. M. York, JCP, 127:194101
61 Hierarchy of Hamiltonians Ab initio quantum chemistry Divide-and-conquer methods Fragment based approaches Effective fragment potential Molecular fractionation with conjugate caps Density fragment interaction QM/MM methods MM methods Coarse-grained methods Mesoscopic and continuum methods
62 A matrix interaction scheme
63 QM/MM: direct-sampling based free energy calculation QM SCF calculation at every MD/MC step E( rqm, rmm ) E( rqm, rmm ) E( rqm, rmm ),, and r r ψ H + H QM QM QM The forces on QM atoms r The forces on MM atoms ψ + = ψ r MM, i MM, i MM ( H H /, ) QM / MM, ele QM QM MM ele QM ψ HQM + HQM/ MM, ele ψ ψ qi/ r ri ψ = r r ψ
64 QM/MM: direct-sampling based free energy calculation Reaction coordinate R is known, computes P( R), probability distribution of R or E( rqm, rmm ), acting force on R R Umbrella sampling Thermodynamic integration Too costly for ab initio QM in general! Choice of R strongly affects the convergence
65 Example: electron-transfer reaction [ Fe( H O) ] [ Fe( H O) ] Δ A = E( η) η η dη Janak theorem Removing electron Adding electron E( η) η E( η) η = ε = ε HOMO LUMO
66 QM/MM: free-energy perturbation based free energy calculation Discrete representation of a reaction path by the QM conformations {r QM } i, i=1,,n Free energy difference between two adjacent QM conformations is computed by free energy perturbation 1 Δ A = ln exp( β E { } E { } ) β r r d r ( ) { } i + i 1 QM i+ 1 QM i MM i Zhang, Liu, & Yang, J. Chem. Phys., 2000
67 QM/MM: two computational challenges QM calculation 50 ~ 100 atoms, 500 ~ 1000 basis functions, 10 ~ 20 energy+gradient/day MM conformation Ensemble of reaction paths Convergence of MM sampling
68 QM-PMF: beyond single potential energy surface PMF of QM coordinates Integration yields free energy PMF Gradient
69 QM/MM: zero-order ESP charge approximation Decompose the QM/MM total electrostatic energy QM/MM electrostatic energy E QM internal energy ( r, r ) ESP QM / MM QM MM = q Q ( r, r ) j r - r j MMi QM QM, i MM, j i QM MM ESP-fitted QM charges E ( r, r ) = Ψ H Ψ E ( r, r ) ESP 1 QM MM eff QM / MM QM MM E 0 0 ( rqm, rmm ) E1 ( rqm, rmm ) + qq 0 0 j i( rqm, rmm) r - r j MMi QM QM, i MM, j + E ( r, r ) + E ( r, r ) + E ( r ) QM / MM, vdw QM MM QM / MM,cov QM MM MM MM
70 QM/MM: first-order polarizable QM ESP charge approximation First-order polarization of ESP charge Q ( r, v ) = Q, + χ ij vmm ( rqm j ) vmm ( r QM j ) + κ r r i QM MM QM i,, ij QM, j QM, j j QM j QM χ ij Q VMM ( QM, j ) r QM, i = N κ ij Q QM, i = rqm, j N E ( r, v ) = Q ( r, v ) v r ( ) ESP QM / MM QM MM i QM MM MM QM, i i QM = Q + χ v ( r ) v ( r ) v r i QM ( ) 0 0 QM, i ij MM QM, j MM QM, j MM QM, i j QM ( ) ( ) ( ) E ( r, v ) = E ( r, v ) χ v r v r v r QM MM 1 QM MM ij MM QM, i MM QM, j MM QM, j i QM j QM vmm ( QM, i ) v (, ) MM QM i χ ij vmm ( QM, j ) vmm ( QM, j ) 2 r r r r i QM j QM Lu & Yang, J. Chem. Phys., 2004
71 QM/MM: first-order polarizable QM ESP charge approximation First-order polarization of ESP charge E ( r, r ) = E ( r, v ) + E ( r, v ) ESP QM MM 1 QM MM QM / MM QM MM + E ( r, r ) + E ( r, r ) + E ( r ) QM / MM, vdw QM MM QM / MM,cov QM MM MM MM = E ( r, v ) + Q v ( r ) QM MM i MM QM, i i QM vmm ( QM, i ) v (, ) MM QM i χ ij vmm ( QM, j ) vmm ( QM, j ) 2 r r r r + E i, j QM QM / MM, vdw ( r, r ) + E ( r, r ) + E ( r ) QM MM QM / MM,cov QM MM MM MM E ( r, r ) E ( r, v ) v ( r ) = QM MM 1 QM MM 0 MM QM, i QQM, i rqm, i rqm, i rqm, i j QM κ 0 ij vmm ( rqm, j ) vm M( rqm, j) vmm( rqm, i) + v ( r ) v ( r ) χ v ( r ) v ( r ) r MM QM, j MM QM, j ji MM QM, i MM QM, i j QM ( EQM / MM, vdw( rqm, rm M ) + EQM / MM,cov( rqm, rmm )) r QM, i QM, i Lu & Yang, J. Chem. Phys., 2004
72 QM/MM-MFEP: straightforward optimization 1) Set initial QM structure r (0) QM, set counter n=0; 2) Increase counter n=n+1; a) Carry out MD simulation of the MM environment with QM geometry fixed at r (n-1) QM n { r τ τ = N} ( ) ( ), 1,..., MD sampling based on E ( r ) MM ref MM b) Carry out one step of QM optimization based on PMF and PMF gradient ( n-1) A( r ) ( n) ( n-1) QM QM one step in the QM optimization based on A( rqm ) and rqm r c) Update QM geometry 3) Go to step (2) until converged Hu, Lu, & Yang, J. Chem. Theory, & Comput., 2007
73 QM/MM-MFEP: path optimization NEB, Ayala-Schlegel MEP, QSM General procedures: 1. Optimize the structures of the reactant and product states, separately 2. General initial guess of the reaction path, by coordinate driving or interpolation 3. Fix the structures of the RS and PS, carry out MFEP optimization 4. Optimize the TS structure Hu, Lu, & Yang, J. Chem. Theory, & Comput., 2007
74 QM/MM Free energy perturbation Δ E = E ( r, r ) + E ( r, r ) + E ( r, r ) i j QM QM, j MM, i QM / MM, ele QM, j MM, i QM / MM, nuc QM, j MM, i - E ( r, r )- E ( r, r )- E ( r, r ) QM QM, i MM, i QM / MM, ele QM, i MM, i QM / MM, nuc QM, i MM, i + E ( r, r )- E ( r, r ) E QM / MM, vdw QM, j MM, i QM / MM, vdw QM, i MM, i ( r, r )- E ( r, r ) + QM / MM,cova QM, j MM, i QM/ MM,cov a QM, i MM, i
75 QM/MM-MFEP: bottleneck Time of MD simulations 50,000 atoms, ~ 150 ps /day (1 fs stepsize) 130,000 atoms, ~ 50 ps/day (1 fs stepsize) Time for QM calculations 800 basis functions, 10~20 energy+gradient/day Hundreds of QM optimization steps are needed Half day for one MD + QM calculation Too many repetitive MD simulations whose samplings are not efficiently utilized
76 QM/MM-FE: efficient iterative, sequential optimizations 1) Set initial QM structure r (0) QM, set counter n=0; 2) Increase counter n=n+1; a) Carry out MM optimization with QM geometry fixed at r QM (n-1) b) Carry out QM optimization with the MM conformation fixed at r QM (n-1) c) Update QM geometry 3) Go to step (2) until converged n ( ) r = arg min E r, r ( n) ( 1) MM QM MM r MM (n) ( ) r = arg min E r, r ( n) QM r QM MM QM In certain cases, sequential optimization may have unique advantages. In the QM/MM-FE method, sequential algorithm significantly reduces the number of QM calculations. Zhang, Liu, & Yang, J. Chem. Phys., 2000
77 Ab Initio QM/MM Minimum Free-Energy Path MD t
78 Iterative, sequential optimization QM MM
79 Iterative, sequential QM/MM-MFEP: algorithm 1) Set initial QM structure r (0,0) QM, set counter n=0; 2) Increase counter n=n+1; a) Carry out MD simulation of the MM environment with QM geometry fixed at r (n-1,0) QM n { r τ τ = N} ( ) ( ), 1,..., MD sampling based on E ( r ) MM ref MM b) Carry out QM optimization with the MM ensemble fixed at c) Update QM geometry 3) Go to step (2) until converged { r ( ( ) MM τ } { β ( τ ) ( τ ) } N ( ni, ) 1 1 ( ni, ) ( n,0) A ( rqm ) = Aref ln exp E QM, MM ( ) Eref MM ( ) β N r r r τ = 1 N E A ( r ) = rqm ( n) ( rqm, rmm ( τ )) ( ni, ) ( n,0) exp { β E (, ( τ) ) E ( ( τ) ) r r r r } QM, i N ( ni, ) ( n,0) exp { β E ( rqm, rmm ( τ) ) Eref ( rmm ( τ) ) } ( ni, ) QM MM ref MM QM τ = 1 τ = 1
80 Ab initio QM/MM minimum free energy path Advantages Complicated solution and enzyme reactions become gas-phase-like Removes the path-dependence of initial conformations Can be applied to solution reactions Adequate statistical sampling Reaction path optimization without explicitly defining the reaction coordinate Hu & Yang, Annu. Rev. Phys. Chem., 2008
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationStructural Bioinformatics (C3210) Molecular Mechanics
Structural Bioinformatics (C3210) Molecular Mechanics How to Calculate Energies Calculation of molecular energies is of key importance in protein folding, molecular modelling etc. There are two main computational
More informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More informationFree energy simulations
Free energy simulations Marcus Elstner and Tomáš Kubař January 14, 2013 Motivation a physical quantity that is of most interest in chemistry? free energies Helmholtz F or Gibbs G holy grail of computational
More informationCE 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 informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 09 Quantum Mechanics/Molecular Mechanics (QM/MM) Techniques Fundamentals of Sustainable Technology These notes created by David Keffer, University of Tennessee, Knoxville,
More informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationLecture 11: Potential Energy Functions
Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions
More informationFree energy calculations
Free energy calculations Berk Hess May 5, 2017 Why do free energy calculations? The free energy G gives the population of states: ( ) P 1 G = exp, G = G 2 G 1 P 2 k B T Since we mostly simulate in the
More informationAdvanced in silico drug design
Advanced in silico drug design RNDr. Martin Lepšík, Ph.D. Lecture: Advanced scoring Palacky University, Olomouc 2016 1 Outline 1. Scoring Definition, Types 2. Physics-based Scoring: Master Equation Terms
More informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationBiomolecular modeling. Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010)
Biomolecular modeling Marcus Elstner and Tomáš Kubař Theoretical Chemistry, TU Braunschweig (Dated: December 10, 2010) IX. FREE ENERGY SIMULATIONS When searching for a physical quantity that is of most
More informationAN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University
AN INTRODUCTION TO QUANTUM CHEMISTRY Mark S. Gordon Iowa State University 1 OUTLINE Theoretical Background in Quantum Chemistry Overview of GAMESS Program Applications 2 QUANTUM CHEMISTRY In principle,
More informationLecture 19: Free Energies in Modern Computational Statistical Thermodynamics: FEP and Related Methods
Statistical Thermodynamics Lecture 19: Free Energies in Modern Computational Statistical Thermodynamics: FEP and Related Methods Dr. Ronald M. Levy ronlevy@temple.edu Free energy calculations Free energy
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationFree energy calculations using molecular dynamics simulations. Anna Johansson
Free energy calculations using molecular dynamics simulations Anna Johansson 2007-03-13 Outline Introduction to concepts Why is free energy important? Calculating free energy using MD Thermodynamical Integration
More informationChemistry 334 Part 2: Computational Quantum Chemistry
Chemistry 334 Part 2: Computational Quantum Chemistry 1. Definition Louis Scudiero, Ben Shepler and Kirk Peterson Washington State University January 2006 Computational chemistry is an area of theoretical
More informationFree energy calculations and the potential of mean force
Free energy calculations and the potential of mean force IMA Workshop on Classical and Quantum Approaches in Molecular Modeling Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationClose agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations
Close agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations Alexandre V. Morozov, Tanja Kortemme, Kiril Tsemekhman, David Baker
More informationMolecular Simulation I
Molecular Simulation I Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Jeffry D. Madura Department of Chemistry & Biochemistry Center for Computational Sciences
More informationChemistry 4560/5560 Molecular Modeling Fall 2014
Final Exam Name:. User s guide: 1. Read questions carefully and make sure you understand them before answering (if not, ask). 2. Answer only the question that is asked, not a different question. 3. Unless
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More information3rd Advanced in silico Drug Design KFC/ADD Molecular mechanics intro Karel Berka, Ph.D. Martin Lepšík, Ph.D. Pavel Polishchuk, Ph.D.
3rd Advanced in silico Drug Design KFC/ADD Molecular mechanics intro Karel Berka, Ph.D. Martin Lepšík, Ph.D. Pavel Polishchuk, Ph.D. Thierry Langer, Ph.D. Jana Vrbková, Ph.D. UP Olomouc, 23.1.-26.1. 2018
More informationComputational Chemistry. An Introduction to Molecular Dynamic Simulations
Computational Chemistry An Introduction to Molecular Dynamic Simulations Computational chemistry simulates chemical structures and reactions numerically, based in full or in part on the fundamental laws
More information510 Subject Index. Hamiltonian 33, 86, 88, 89 Hamilton operator 34, 164, 166
Subject Index Ab-initio calculation 24, 122, 161. 165 Acentric factor 279, 338 Activity absolute 258, 295 coefficient 7 definition 7 Atom 23 Atomic units 93 Avogadro number 5, 92 Axilrod-Teller-forces
More informationFree energy calculations
Free energy calculations Jochen Hub & David van der Spoel Overview Free energies and Probabilities Thermodynamic cycles (Free energy perturbation (FEP)) Thermodynamic integration (TI) (Jarzynski equality
More informationICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below
ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below Introduction In statistical physics Monte Carlo methods are considered to have started in the Manhattan project (1940
More informationMolecular mechanics. classical description of molecules. Marcus Elstner and Tomáš Kubař. April 29, 2016
classical description of molecules April 29, 2016 Chemical bond Conceptual and chemical basis quantum effect solution of the SR numerically expensive (only small molecules can be treated) approximations
More informationAtomic and molecular interaction forces in biology
Atomic and molecular interaction forces in biology 1 Outline Types of interactions relevant to biology Van der Waals interactions H-bond interactions Some properties of water Hydrophobic effect 2 Types
More information8.333: Statistical Mechanics I Problem Set # 5 Due: 11/22/13 Interacting particles & Quantum ensembles
8.333: Statistical Mechanics I Problem Set # 5 Due: 11/22/13 Interacting particles & Quantum ensembles 1. Surfactant condensation: N surfactant molecules are added to the surface of water over an area
More informationMolecular simulation and structure prediction using CHARMM and the MMTSB Tool Set Free Energy Methods
Molecular simulation and structure prediction using CHARMM and the MMTSB Tool Set Free Energy Methods Charles L. Brooks III MMTSB/CTBP 2006 Summer Workshop CHARMM Simulations The flow of data and information
More informationQ-Chem 5: Facilitating Worldwide Scientific Breakthroughs
Q-Chem 5: Facilitating Worldwide Scientific Breakthroughs Founded in 1993, Q-Chem strives to bring its customers state-ofthe-art methods and algorithms for performing quantum chemistry calculations. Cutting-edge
More information4 th Advanced in silico Drug Design KFC/ADD Molecular Modelling Intro. Karel Berka, Ph.D.
4 th Advanced in silico Drug Design KFC/ADD Molecular Modelling Intro Karel Berka, Ph.D. UP Olomouc, 21.1.-25.1. 2019 Motto A theory is something nobody believes, except the person who made it An experiment
More informationThis is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry.
1 Computational Chemistry (Quantum Chemistry) Primer This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry. TABLE OF CONTENTS Methods...1 Basis
More informationBiomolecular modeling I
2016, December 6 Biomolecular structure Structural elements of life Biomolecules proteins, nucleic acids, lipids, carbohydrates... Biomolecular structure Biomolecules biomolecular complexes aggregates...
More informationMolecular Dynamics. A very brief introduction
Molecular Dynamics A very brief introduction Sander Pronk Dept. of Theoretical Physics KTH Royal Institute of Technology & Science For Life Laboratory Stockholm, Sweden Why computer simulations? Two primary
More informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationStatistical Mechanics for Proteins
The Partition Function From Q all relevant thermodynamic properties can be obtained by differentiation of the free energy F: = kt q p E q pd d h T V Q ), ( exp 1! 1 ),, ( 3 3 3 ),, ( ln ),, ( T V Q kt
More informationWhy Is Molecular Interaction Important in Our Life
Why Is Molecular Interaction Important in ur Life QuLiS and Graduate School of Science iroshima University http://www.nabit.hiroshima-u.ac.jp/iwatasue/indexe.htm Suehiro Iwata Sept. 29, 2007 Department
More informationDevelopment and Application of Combined Quantum Mechanical and Molecular Mechanical Methods
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Student Research Projects, Dissertations, and Theses - Chemistry Department Chemistry, Department of Fall 12-2014 Development
More informationQuantum mechanics can be used to calculate any property of a molecule. The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is,
Chapter : Molecules Quantum mechanics can be used to calculate any property of a molecule The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is, E = Ψ H Ψ Ψ Ψ 1) At first this seems like
More informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationk θ (θ θ 0 ) 2 angles r i j r i j
1 Force fields 1.1 Introduction The term force field is slightly misleading, since it refers to the parameters of the potential used to calculate the forces (via gradient) in molecular dynamics simulations.
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
More informationComputational Modeling of Protein-Ligand Interactions
Computational Modeling of Protein-Ligand Interactions Steven R. Gwaltney Department of Chemistry Mississippi State University Mississippi State, MS 39762 Auguste Comte, 1830 Every attempt to refer chemical
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Length Scales (c) ICAMS: http://www.icams.de/cms/upload/01_home/01_research_at_icams/length_scales_1024x780.png Goals for today: Background
More informationHyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea
The Scaled Effective Solvent Method for Predicting the Equilibrium Ensemble of Structures with Analysis of Thermodynamic Properties of Amorphous Polyethylene Glycol-Water Mixtures Hyeyoung Shin a, Tod
More informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
More information74 these states cannot be reliably obtained from experiments. In addition, the barriers between the local minima can also not be obtained reliably fro
73 Chapter 5 Development of Adiabatic Force Field for Polyvinyl Chloride (PVC) and Chlorinated PVC (CPVC) 5.1 Introduction Chlorinated polyvinyl chloride has become an important specialty polymer due to
More informationONETEP PB/SA: Application to G-Quadruplex DNA Stability. Danny Cole
ONETEP PB/SA: Application to G-Quadruplex DNA Stability Danny Cole Introduction Historical overview of structure and free energy calculation of complex molecules using molecular mechanics and continuum
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 informationCoupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation
Coupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation Bo Li Math Dept & CTBP, UCSD Li-Tien Cheng (Math, UCSD) Zhongming Wang (Math & Biochem, UCSD) Yang Xie (MAE,
More informationScuola di Chimica Computazionale
Societa Chimica Italiana Gruppo Interdivisionale di Chimica Computazionale Scuola di Chimica Computazionale Introduzione, per Esercizi, all Uso del Calcolatore in Chimica Organica e Biologica Modellistica
More informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
More informationComputing free energy: Thermodynamic perturbation and beyond
Computing free energy: Thermodynamic perturbation and beyond Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum ls
More informationStructure of Cement Phases from ab initio Modeling Crystalline C-S-HC
Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model
More informationThe Molecular Dynamics Method
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation
More informationFree-energy calculations
Measuring free-energy differences using computer simulations Theoretical and Computational Biophysics Group University of Illinois, Urbana Champaign and Équipe de dynamique des assemblages membranaires,
More informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry room B430, Chemicum 4th floor vesa.hanninen@helsinki.fi September 3, 2013 Introduction and theoretical backround September
More informationElectronic structure calculations: fundamentals George C. Schatz Northwestern University
Electronic structure calculations: fundamentals George C. Schatz Northwestern University Electronic Structure (often called Quantum Chemistry) calculations use quantum mechanics to determine the wavefunctions
More informationTransition states and reaction paths
Transition states and reaction paths Lab 4 Theoretical background Transition state A transition structure is the molecular configuration that separates reactants and products. In a system with a single
More informationReactive potentials and applications
1.021, 3.021, 10.333, 22.00 Introduction to Modeling and Simulation Spring 2011 Part I Continuum and particle methods Reactive potentials and applications Lecture 8 Markus J. Buehler Laboratory for Atomistic
More informationSubject of the Lecture:
Subject of the Lecture: Conceptual basis for the development of force fields. Implementation/validation Water - a worked example Extensions - combining molecular mechanics and quantum mechanics (QM/MM)
More informationSession 7 Overview: Part A I. Prediction of Vibrational Frequencies (IR) Part B III. Prediction of Electronic Transitions (UV-Vis) IV.
Session 7 Overview: Part A I. Prediction of Vibrational Frequencies (IR) II. Thermochemistry Part B III. Prediction of Electronic Transitions (UV-Vis) IV. NMR Predictions 1 I. Prediction of Vibrational
More informationGases and the Virial Expansion
Gases and the irial Expansion February 7, 3 First task is to examine what ensemble theory tells us about simple systems via the thermodynamic connection Calculate thermodynamic quantities: average energy,
More informationIntermolecular Forces in Density Functional Theory
Intermolecular Forces in Density Functional Theory Problems of DFT Peter Pulay at WATOC2005: There are 3 problems with DFT 1. Accuracy does not converge 2. Spin states of open shell systems often incorrect
More informationFormation Mechanism and Binding Energy for Icosahedral Central Structure of He + 13 Cluster
Commun. Theor. Phys. Beijing, China) 42 2004) pp. 763 767 c International Academic Publishers Vol. 42, No. 5, November 5, 2004 Formation Mechanism and Binding Energy for Icosahedral Central Structure of
More informationDENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY
DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY A TUTORIAL FOR PHYSICAL SCIENTISTS WHO MAY OR MAY NOT HATE EQUATIONS AND PROOFS REFERENCES
More informationElectron States of Diatomic Molecules
IISER Pune March 2018 Hamiltonian for a Diatomic Molecule The hamiltonian for a diatomic molecule can be considered to be made up of three terms Ĥ = ˆT N + ˆT el + ˆV where ˆT N is the kinetic energy operator
More informationLecture 9. Hartree Fock Method and Koopman s Theorem
Lecture 9 Hartree Fock Method and Koopman s Theorem Ψ(N) is approximated as a single slater determinant Φ of N orthogonal One electron spin-orbitals. One electron orbital φ i = φ i (r) χ i (σ) χ i (σ)
More informationMolecular Dynamics Simulation of a Nanoconfined Water Film
Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.
More informationGherman Group Meeting. Thermodynamics and Kinetics and Applications. June 25, 2009
Gherman Group Meeting Thermodynamics and Kinetics and Applications June 25, 2009 Outline Calculating H f, S, G f Components which contribute to H f, S, G f Calculating ΔH, ΔS, ΔG Calculating rate constants
More informationMolecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia
Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties
More informationSame idea for polyatomics, keep track of identical atom e.g. NH 3 consider only valence electrons F(2s,2p) H(1s)
XIII 63 Polyatomic bonding -09 -mod, Notes (13) Engel 16-17 Balance: nuclear repulsion, positive e-n attraction, neg. united atom AO ε i applies to all bonding, just more nuclei repulsion biggest at low
More informationAb initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier
Ab initio calculations for potential energy surfaces D. Talbi GRAAL- Montpellier A theoretical study of a reaction is a two step process I-Electronic calculations : techniques of quantum chemistry potential
More informationQuantum Mechanics - Molecular Mechanics (QM/MM) CHEM 430
Quantum Mechanics - Molecular Mechanics (QM/MM), CHEM 430 Quantum Mechanics In theory, a very accurate treatment of the system Largely ab initio, i.e. parameter-free Very expensive typically scales as
More informationMolecular Mechanics. I. Quantum mechanical treatment of molecular systems
Molecular Mechanics I. Quantum mechanical treatment of molecular systems The first principle approach for describing the properties of molecules, including proteins, involves quantum mechanics. For example,
More informationLecture 5: More about one- Final words about the Hartree-Fock theory. First step above it by the Møller-Plesset perturbation theory.
Lecture 5: More about one- determinant wave functions Final words about the Hartree-Fock theory. First step above it by the Møller-Plesset perturbation theory. Items from Lecture 4 Could the Koopmans theorem
More informationCalculations of band structures
Chemistry and Physics at Albany Planning for the Future Calculations of band structures using wave-function based correlation methods Elke Pahl Centre of Theoretical Chemistry and Physics Institute of
More informationMolecular Aggregation
Molecular Aggregation Structure Analysis and Molecular Simulation of Crystals and Liquids ANGELO GAVEZZOTTI University of Milano OXFORD UNIVERSITY PRESS Contents PART I FUNDAMENTALS 1 The molecule: structure,
More informationIntroduction to Hartree-Fock Molecular Orbital Theory
Introduction to Hartree-Fock Molecular Orbital Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Origins of Mathematical Modeling in Chemistry Plato (ca. 428-347
More information2m 2 Ze2. , where δ. ) 2 l,n is the quantum defect (of order one but larger
PHYS 402, Atomic and Molecular Physics Spring 2017, final exam, solutions 1. Hydrogenic atom energies: Consider a hydrogenic atom or ion with nuclear charge Z and the usual quantum states φ nlm. (a) (2
More informationBiomolecules are dynamic no single structure is a perfect model
Molecular Dynamics Simulations of Biomolecules References: A. R. Leach Molecular Modeling Principles and Applications Prentice Hall, 2001. M. P. Allen and D. J. Tildesley "Computer Simulation of Liquids",
More informationPart III: Theoretical Surface Science Adsorption at Surfaces
Technische Universität München Part III: Theoretical Surface Science Adsorption at Surfaces Karsten Reuter Lecture course: Solid State Theory Adsorption at surfaces (T,p) Phase II Phase I Corrosion Growth
More informationMachine learning the Born-Oppenheimer potential energy surface: from molecules to materials. Gábor Csányi Engineering Laboratory
Machine learning the Born-Oppenheimer potential energy surface: from molecules to materials Gábor Csányi Engineering Laboratory Interatomic potentials for molecular dynamics Transferability biomolecular
More informationJack Simons Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationThe Potential Energy Surface
The Potential Energy Surface In this section we will explore the information that can be obtained by solving the Schrödinger equation for a molecule, or series of molecules. Of course, the accuracy of
More informationHW 1 CHEM 362. Available: Jan. 16, 2008 Due: Jan. 25, 2008
HW 1 CHEM 362 Available: Jan. 16, 2008 Due: Jan. 25, 2008 1. Write an equation that can be used to define the mean S-F bond energy in SF 6. How is this value likely to be related in magnitude to the energy
More informationIntroduction to DFTB. Marcus Elstner. July 28, 2006
Introduction to DFTB Marcus Elstner July 28, 2006 I. Non-selfconsistent solution of the KS equations DFT can treat up to 100 atoms in routine applications, sometimes even more and about several ps in MD
More informationComputational Chemistry - MD Simulations
Computational Chemistry - MD Simulations P. Ojeda-May pedro.ojeda-may@umu.se Department of Chemistry/HPC2N, Umeå University, 901 87, Sweden. May 2, 2017 Table of contents 1 Basics on MD simulations Accelerated
More informationSpecific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate
Specific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate A. Arslanargin, A. Powers, S. Rick, T. Pollard, T. Beck Univ Cincinnati Chemistry Support: NSF, OSC TSRC 2016 November 2, 2016 A. Arslanargin,
More informationInteratomic Potentials. The electronic-structure problem
Interatomic Potentials Before we can start a simulation, we need the model! Interaction between atoms and molecules is determined by quantum mechanics: Schrödinger Equation + Born-Oppenheimer approximation
More informationA new extension of QM/MM methods: the adaptive buffered-force QM/MM method
A new extension of QM/MM methods: the adaptive buffered-force QM/MM method Letif Mones Engineering Department, University of Cambridge lam81@cam.ac.uk Overview Basic concept of the QM/MM methods MM vs.
More informationQUANTUM CHEMISTRY FOR TRANSITION METALS
QUANTUM CHEMISTRY FOR TRANSITION METALS Outline I Introduction II Correlation Static correlation effects MC methods DFT III Relativity Generalities From 4 to 1 components Effective core potential Outline
More informationBiomolecular modeling I
2015, December 15 Biomolecular simulation Elementary body atom Each atom x, y, z coordinates A protein is a set of coordinates. (Gromacs, A. P. Heiner) Usually one molecule/complex of interest (e.g. protein,
More informationIntroduction to multiconfigurational quantum chemistry. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Notations
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
More information