QM/MM approaches in ab initio molecular dynamics
|
|
- Christine Ramsey
- 5 years ago
- Views:
Transcription
1 CP2K UK Workshop 2014! August, Imperial College, London /MM approaches in ab initio molecular dynamics Marcella Iannuzzi! University of Zurich!! Teodoro Laino! IBM Zurich Research Laboratory
2 MM:overview GEEP Application: covd Conclusion MM Schemes Outline Overview of the /MM methodology! Available /MM Electrostatic Schemes! GEEP: CP2K /MM driver! Charged Oxygen Vacancies in SiO2
3 Nobel Prize in Chemistry 2013 Martin Karplus, Harvard U., Cambridge, MA, USA Micheal Levitt, Stanford U., Stanford, CA, USA Arieh Warshel, U. Southern Ca., Los Angeles, CA, USA Development of Multiscale Models of Complex Chemical Systems
4 /MM Combine and MM Combine Quantum Mechanics and Molecular Mechanics /MM full atomistic by classical FF Separate system in and MM part, single out electro Vtreated (R) =byv (R) + VMM (R) + Vint (R) Define potential energy surface
5 Combine and MM : modelling of electronic rearrangements! MM: efficient inclusion of wider environment! Choice of method (semi empirical, DFT, QC)! Choice of the force field! Partitioning and treatment of the boundary
6 MM:overview GEEP Application: covd Conclusion MM Schemes 738,000 atoms - 50 nanometers Ligand binding affinity in docking! Free energy simulations! Complex biomolecular structures P.D. Blood and G.A. Voth, PNAS, 103, 2006, pp
7 MM:overview GEEP Application: covd Conclusion MM Schemes Environment effects on reaction energetics
8 MM:overview GEEP Application: covd Conclusion MM Schemes MM: overview 0.11 million atoms! 5 regions: effects of O implantation into Si! adaptive regions simox technology Yoshio Tanaka (AIST) and Aiichiro Nakano (USC)
9 MM Environment U(R N )= X i U 1 (R i )+ X i X X X U 3 (R i, R j, R k )+... j>i U 2 (R i, R j )+X i j>i k>j U(R N, n p )
10 MM Environment U(R N )= X i U 1 (R i )+ X i X X X U 3 (R i, R j, R k )+... j>i U 2 (R i, R j )+X i j>i k>j X X X X "!! #! U(R N, n p ) U(R N ) = + + X j>i X i2bonds X s2torsion k (b) i 2 (l i l i,0 ) 2 + X j2angles V s 2 (1 + cos(n s! s)) 4" ij " ij R ij 12 ij R ij 6 # k (a) j 2 ( j j,0 ) 2 + q iq j 4" 0 R ij! :[(k (b),l 0 ) bon ;(k (a), 0 ) ang ;(V s, s) tor ;(", ) pair ; q at ]
11 Topology RESI ALA 0.00 GROUP ATOM N NH1-0.47! ATOM HN H 0.31! HN-N ATOM CA CT1 0.07! HB1 ATOM HA HB 0.09! / GROUP! HA-CA--CB-HB2 ATOM CB CT3-0.27! \ ATOM HB1 HA 0.09! HB3 ATOM HB2 HA 0.09! O=C ATOM HB3 HA 0.09! GROUP! ATOM C C 0.51 ATOM O O BOND CB CA N HN N CA BOND C CA C +N CA HA CB HB1 CB HB2 CB HB3 DOUBLE O C IMPR N -C CA HN C CA +N O DONOR HN N ACCEPTOR O C IC -C CA *N HN IC -C N CA C IC N CA C +N IC +N CA *C O IC CA C +N +CA IC N C *CA CB IC N C *CA HA IC C CA CB HB IC HB1 CA *CB HB IC HB1 CA *CB HB
12 MM CP2K input &FORCE_EVAL METHOD FIST!!! &MM &FORCEFIELD PARM_FILE_NAME acn.pot PARMTYPE CHM &CHARGE ATOM CT CHARGE &END CHARGE &CHARGE ATOM YC CHARGE &END CHARGE &CHARGE ATOM YN CHARGE &END CHARGE &CHARGE ATOM HC CHARGE &END CHARGE &END FORCEFIELD!!! &POISSON &EWALD EWALD_TYPE SPME ALPHA.44 GMAX 32 O_SPLINE 6 &END EWALD &END POISSON &END MM &SUBSYS &CELL ABC &END CELL &TOPOLOGY CONNECTIVITY PSF CONN_FILE_NAME acn_topology.psf COORD_FILE_NAME acn_topology.pdb COORDINATE pdb &END TOPOLOGY &END SUBSYS STRESS_TENSOR ANALYTICAL &END FORCE_EVAL
13 Subtractive /MM E total = E MM,tot + E () E MM() MM FF also for active region! density not polarised
14 Subtractive /MM E total = E MM,tot + E () E MM() MM FF also for active region! &MULTIPLE_FORCE_EVALS FORCE_EVAL_ORDER &END MULTIPLE_FORCE_EVALS! &FORCE_EVAL METHOD MIXED &MIXED MIXING_TYPE GENMIX &GENERIC # X: Energy force_eval 2 # Y: Energy force_eval 3 # Z: Energy force_eval 4 MIXING_FUNCTION X+Y-Z VARIABLES X Y Z &END GENERIC &END MIXED &SUBSYS &TOPOLOGY CONNECTIVITY PSF CONN_FILE_NAME topo.psf COORD_FILE_NAME totsys.xyz &END TOPOLOGY &CELL ABC &END CELL &END SUBSYS &END FORCE_EVAL density not polarised! &FORCE_EVAL METHOD FIST &MM &END MM &SUBSYS &TOPOLOGY CONNECTIVITY PSF CONN_FILE_NAME topo.psf COORD_FILE_NAME totsys.xyz &END TOPOLOGY &CELL ABC &END CELL &END SUBSYS &END FORCE_EVAL! &FORCE_EVAL METHOD QS &DFT &END DFT!! &SUBSYS &TOPOLOGY COORD_FILE_NAME qmsys.xyz &END TOPOLOGY &CELL ABC &END CELL &END SUBSYS &END FORCE_EVAL! &FORCE_EVAL METHOD FIST &MM &END MM &SUBSYS &TOPOLOGY CONNECTIVITY PSF CONN_FILE_NAME qmtopo.psf COORD_FILE_NAME qmsys.xyz &END TOPOLOGY &CELL ABC &END CELL &END SUBSYS &END FORCE_EVAL
15 MM:overview GEEP Application: covd Conclusion MM Schemes Additive /MM E total = E MM,tot + E () + E /MM E MM() = E el MM() + Evdw MM() + Eb MM() Electrostatic coupling is the most involved term! Mechanical embedding possible! Linked atom scheme! vdw might need ad hoc parameterisation
16 MM:overview GEEP Application: covd Conclusion MM Schemes Outline Overview of the /MM methodology! Available /MM Electrostatic Schemes! GEEP: CP2K /MM driver! Charged Oxygen Vacancies in SiO2
17 MM:overview GEEP Application: covd Conclusion MM Schemes Available Electrostatic Schemes on the same cell on which is defined box Cost N MM * P 1 MM box
18 MM:overview GEEP Application: covd Conclusion MM Schemes Available Electrostatic Schemes Spherical Cutoff box Cost N c MM * P1 MM box
19 MM:overview GEEP Application: covd Conclusion MM Schemes Available Electrostatic Schemes Multi-pole! Expansion box MM box A. Laio, J. VandeVondele, U. Rothlisberger, J. Chem. Phys., 116, 2002, pp. 6941
20 MM:overview GEEP Application: covd Conclusion MM Schemes Outline Overview of the /MM methodology! Available /MM Electrostatic Schemes! GEEP: CP2K /MM driver! Charged Oxygen Vacancies in SiO2
21 MM:overview GEEP Application: covd Conclusion MM Schemes /MM
22 MM:overview GEEP Application: covd Conclusion MM Schemes /MM E TOT (R, R MM )=E (R )+E MM (R MM )+E /MM (R, R MM ) E /MM (R, R MM )= X MM q MM Z n(r) r R MM dr + X,MM u vdw (R, R MM )
23 MM:overview GEEP Application: covd Conclusion MM Schemes /MM Gaussians Plane Waves
24 Gaussian charge distribution n(r, R MM )= rc,mm p 3 e ( r R MM /r c,mm ) 2 v MM (r, R MM )= Erf r RMM r c,mm r R MM prevent spill out problem! accelerate calculations of electrostatics
25 MM:overview GEEP Application: covd Conclusion MM Schemes GEEP Sum of functions with different cutoffs, derived from the new Gaussian expansion of the electrostatic potential Energy [a.u.] ṽ MM (g) = 4 g 2 exp g 2 r 2 c,mm 4! Distance
26 MM:overview GEEP Application: covd Conclusion MM Schemes GEEP Sum of functions with different cutoffs, derived from the new Gaussian expansion of the electrostatic potential ṽ MM (g) = 4 g 2 exp g 2 r 2 c,mm 4! G vectors Distance T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 1, 2005, pp
27 MM:overview GEEP Application: covd Conclusion MM Schemes Multigrid Framework N i+1 =8N i T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 1, 2005, pp
28 MM:overview GEEP Application: covd Conclusion MM Schemes Multigrid Framework interpolate restrict interpolate restrict interpolate restrict Cubic Splines T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 1, 2005, pp
29 Collocation in the Box Z E /MM (R, R MM )= n(r, R )V /MM (r, R MM )dr potential on the finest grid V /MM i (r, R MM )= X MM v i MM(r, R MM ) box optimal! grid levels!! 60-80% of time
30 box
31 box
32 compact Gaussian! functions
33 compact Gaussian! functions
34 compact Gaussian! functions
35 compact Gaussian! functions
36
37
38
39
40
41
42
43
44 Scaling ~ Nc3
45 V /MM (r, R MM )= Xfine i=coarse fine Y 1 k=i I k k real space interpolation from coarsest to finest! 1V /MM i (r, R MM )
46 V /MM (r, R MM )= Xfine i=coarse fine Y 1 k=i I k k real space interpolation from coarsest to finest! 1V /MM i (r, R MM )
47 V /MM (r, R MM )= Xfine i=coarse fine Y 1 k=i I k k 1V /MM i (r, R MM )
48 V /MM (r, R MM )= Xfine i=coarse fine Y 1 k=i I k k 1V /MM i (r, R MM )
49 V /MM (r, R MM )= Xfine i=coarse fine Y 1 k=i I k k 1V /MM i (r, R MM )
50 Electrostatic Potential interpolation! 20-40% of time
51 &MM! &CELL ABC &END CELL USE_GEEP_LIB 9 ECOUPL GAUSS! &MM_KIND H RADIUS 0.44 &END MM_KIND &MM_KIND O RADIUS 0.78 &END MM_KIND! &_KIND H MM_INDEX 8 9 &END _KIND &_KIND O MM_INDEX 7 &END _KIND! &END MM! &MM.. &END MM! &DFT. &END DFT! &SUBSYS &CELL ABC &END CELL! &TOPOLOGY COORD_FILE_NAME sys.pdb COORDINATE pdb &END TOPOLOGY &END SUBSYS
52 MM:overview GEEP Application: covd Conclusion MM Schemes Extension to PBC How to handle the electrostatic potential in presence of periodic boundary conditions (PBC)? Ewald Summation scheme:
53 MM:overview GEEP Application: covd Conclusion MM Schemes Extension to PBC How to handle the electrostatic potential in presence of periodic boundary conditions (PBC)? Ewald Summation scheme: Reciprocal space Real space
54 MM:overview GEEP Application: covd Conclusion MM Schemes /MM fully periodic MM MM MM MM MM MM MM MM MM T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 2 (5), 2006, pp
55 Total ES Energy n(r) =n (r)+n MM (r)
56 Total ES Energy n(r) =n (r)+n MM (r) ±n B background charge
57 Total ES Energy n(r) =n (r)+n MM (r) ±n B background charge E TOT = 1 2 Z Z drdr 0 n(r)n(r0 ) r r 0 Z Z E MM = 1 2 E = 1 Z Z 2 E /MM = Z Z drdr 0 (nmm (r)+n B,MM )(n MM (r 0 )+n B,MM ) r r 0 drdr 0 (n (r)+n B, )(n (r 0 )+n B, ) r r 0 drdr 0 (n (r)+n B, )(n MM (r 0 )+n B,MM ) r r 0
58 MM:overview GEEP Application: covd Conclusion MM Schemes MM/MM fully periodic MM MM MM MM MM MM MM MM MM T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 2 (5), 2006, pp
59 MM:overview GEEP Application: covd Conclusion MM Schemes /MM fully periodic MM MM MM MM MM MM MM MM T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 2 (5), 2006, pp
60 MM:overview GEEP Application: covd Conclusion MM Schemes GEEP with PBC T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 2 (5), 2006, pp
61 MM:overview GEEP Application: covd Conclusion MM Schemes GEEP with PBC smooth! coarsest grid T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 2 (5), 2006, pp
62 MM:overview GEEP Application: covd Conclusion MM Schemes /MM real space term V /MM rs (r, R MM )= X LappleL cut MM 2 X 4 X Ng r RMM + L 2 3 A g exp 5 G 2 g T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 2 (5), 2006, pp
63 MM:overview GEEP Application: covd Conclusion MM Schemes /MM reciprocal space term T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 2 (5), 2006, pp
64 MM:overview GEEP Application: covd Conclusion MM Schemes /MM reciprocal space term low cutoff function! only few k vectors! needed T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 2 (5), 2006, pp
65 !!! &MM &CELL ABC &END CELL ECOUPL GAUSS USE_GEEP_LIB 6 &MM_KIND NA RADIUS &END MM_KIND &MM_KIND CL RADIUS &END MM_KIND &PERIODIC GMAX 0.5 &MULTIPOLE EWALD_PRECISION RCUT 8.0 NGRIDS ANALYTICAL_GTERM &END MULTIPOLE &END PERIODIC! &END MM
66 MM:overview GEEP Application: covd Conclusion MM Schemes GEEP Summary GEEP to speed up the evaluation of a function on a grid! The speed up factor is ~ (Nf/Nc) 3 =2 3(Ngrid-1)! Usually 3-4 grid levels are used corresponding to a speed up of ~ 10 2 times faster than the simple collocation algorithm (Interpolations and Restrictions account for a negligible amount of time)! Since the residual function is different from zero only for few k vectors, the sum in reciprocal space is restrained to few points.! Small computational overhead between the fully periodic and non-periodic
67 MM:overview GEEP Application: covd Conclusion MM Schemes Sources of Errors Cutoff of grid level appropriate to the cutoff of the mapped Gaussian (~ points per linear direction)! Error in Cubic Spline interpolation! Cutoff of the coarse grid level comparable to the cutoff of the long range function.
68 MM:overview GEEP Application: covd Conclusion MM Schemes fully periodic T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 2 (5), 2006, pp
69 MM:overview GEEP Application: covd Conclusion MM Schemes fully periodic T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 2 (5), 2006, pp
70 MM:overview GEEP Application: covd Conclusion MM Schemes De-coupling and re-coupling T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 2 (5), 2006, pp
71 MM:overview GEEP Application: covd Conclusion MM Schemes De-coupling and re-coupling T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 2 (5), 2006, pp
72 MM:overview GEEP Application: covd Conclusion MM Schemes Bloechl Scheme!!!!!! Density fitting in g-space of the total density! ˆn(r, R )= X q g (r, R ) Reproduce the correct Long-Range electrostatics! Q l = W = Z Z drr l Y l (n(r, R ) ˆn(r, R )) drr 2 (n(r, R ) ˆn(r, R )) minimise! Decoupling and Recoupling using these charges P. E. Bloechl, J. Chem. Phys., 103 (17), 1995, pp T. Laino, F. Mohamed, A. Laio and M. Parrinello, J. Chem. Th. Comp., 2 (5), 2006, pp
73 MM:overview GEEP Application: covd Conclusion MM Schemes Charged OV Migration of charged oxygen vacancy defects in silica dimer! deloc. el. E δ T. Laino, D. Donadio, I-Feng W. Kuo, Phys. Rev. B, 2007
74 MM:overview GEEP Application: covd Conclusion MM Schemes Charged OV Migration of charged oxygen vacancy defects in silica dimer! deloc. el. E δ E 1 T. Laino, D. Donadio, I-Feng W. Kuo, Phys. Rev. B, 2007
75 MM:overview GEEP Application: covd Conclusion MM Schemes Charged OV Migration of charged oxygen vacancy defects in silica dimer! deloc. el. E δ E 1 E 1 * T. Laino, D. Donadio, I-Feng W. Kuo, Phys. Rev. B, 2007
76 MM:overview GEEP Application: covd Conclusion MM Schemes
77 MM:overview MM Schemes GEEP Application: covd Conclusion
78 MM:overview GEEP Application: covd Conclusion MM Schemes NEB: Minimum Energy Path T. Laino, D. Donadio, I-Feng W. Kuo, Phys. Rev. B, 2007
79 MM:overview GEEP Application: covd Conclusion MM Schemes NEB: Minimum Energy Path T. Laino, D. Donadio, I-Feng W. Kuo, Phys. Rev. B, 2007
80 Image Charge & MM molecule + EAM metal nitrobenzene/au(111) Siepmann Sprik., JCP (1995) 102 Golze Iannuzzi Passerone Hutter, JCTC (2013) : DFT optimized structures of the five investigated molecules adsorbed on a m surface: (a) nitrobenzene, (b) thymine, (c) guanine, (d) benzene on Au(11 (e) water on Pt(111). Color code: orange: C, white: H, blue: N, red : O, y Au, grey: Pt. 52, for large N, computing the single Gaussian potentials v on the full g
81 Image Charge & MM molecule + EAM metal nitrobenzene/au(111) IC (r) = X CImet exp Imet VH (r) + VIC (r) = Z r RImet 2 (r0 ) + IC (r0 ) 0 dr = V0 0 r r IC induce polarization, solved selfconsistently Siepmann Sprik., JCP (1995) 102 Golze Iannuzzi Passerone Hutter, JCTC (2013) : DFT optimized structures of the five investigated molecules adsorbed on a m surface: (a) nitrobenzene, (b) thymine, (c) guanine, (d) benzene on Au(11 (e) water on Pt(111). Color code: orange: C, white: H, blue: N, red : O, y Au, grey: Pt. 52, for large N, computing the single Gaussian potentials v on the full g
82 Image Charge & MM molecule + EAM metal nitrobenzene/au(111) IC (r) = X CImet exp Imet VH (r) + VIC (r) = Z r RImet 2 (r0 ) + IC (r0 ) 0 dr = V0 0 r r IC induce polarization, solved selfconsistently Siepmann Sprik., JCP (1995) 102 Golze Iannuzzi Passerone Hutter, JCTC (2013) : DFT optimized structures of the five investigated molecules adsorbed on a m surface: (a) nitrobenzene, (b) thymine, (c) guanine, (d) benzene on Au(11 IC-MM IC+V:H O, y VH on Pt(111). Color code: orange: C, white: H, blue: N,Vred (e)dft water Au, grey: Pt. 52, for large N, computing the single Gaussian potentials v on the full g
83 IC distribution Z Z (VH (r) (VH (r) + VIC (r) V0 ) gi (r) + X J CJ V0 ) gi (r) = Z Z gj (r0 )gi (r) drdr0 0 r r linear set of eq. (CG iterative scheme) arameters. (a) Dependence of total energy and x-gradients (of an etal atom) on parameter. The x-gradients refer to the image charge n, see Eq. 9. (b) Execution time of a single point calculation of thymine performed on Cray XK6 system. Note that a double logarithmic scale. ations and performance ariable in the IC-/MM scheme is the width of the Gaussian charge 53
84 H2O cluster on Pt(111) H 2 O, Pt EAM, H 2 O-Pt Siepmann-Sprik + IC 1H 2 O 2H 2 O 12H 2 O kj/mol E int E ads E int E ads E H bond E int E ads E H bond /MM IC-/MM full DFT
85 1 2 e 27 of 36 Journal of Chemical Theory and Computation honeycomb arrangement 10 70% on-top site occupied density Figure 11: 29 Structural and electronic properties obtained for a water film on Pt(111). (a) Typical snapshot of30 the water/pt(111) simulation cell. (b) Snapshot of the water adlayer. (c) Plane-averaged water density rw along the surface normal. The arrows indicate the division of the water film in 31 4 different layers. (d) Plane-averaged particle density rp for O and H atoms. (e) Plane-integrated 32 difference Dr 1D between water and Pt(111). electronic density elec the water molecules in these regions are expected to be massively overstructured. Further away 35 Liquid Water at Pt(111) from the surface, the density oscillations decay rapidly. For distances larger than 10 Å, the density 55
86 H-bond distribution 56
CP2K: Selected Developments
CP2K: Selected Developments Jürg Hutter Department of Chemistry University of Zurich Outline Introduction History and Performance Current and Future Developments Post-Hartree-Fock Methods GW Methods RPA
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationAdvanced Molecular Molecular Dynamics
Advanced Molecular Molecular Dynamics Technical details May 12, 2014 Integration of harmonic oscillator r m period = 2 k k and the temperature T determine the sampling of x (here T is related with v 0
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz January 17, 2011 CP-PAW/COSMO Interface Mixed quantum/ classical molecular dynamics (QM/MM) Quantum mechanics is computationally expensive
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 informationAb initio molecular dynamics
Ab initio molecular dynamics Kari Laasonen, Physical Chemistry, Aalto University, Espoo, Finland (Atte Sillanpää, Jaakko Saukkoriipi, Giorgio Lanzani, University of Oulu) Computational chemistry is a field
More informationCCSD(T) benchmarks of non-equilibrium water clusters: the importance of monomer deformation
CCSD(T) benchmarks of non-equilibrium water clusters: the importance of monomer deformation Biswajit Santra 1, Angelos Michaelides 1,2, and Matthias Scheffler 1 1 Fritz-Haber-Institut der MPG, Berlin,
More informationLarge Scale Electronic Structure Calculations
Large Scale Electronic Structure Calculations Jürg Hutter University of Zurich 8. September, 2008 / Speedup08 CP2K Program System GNU General Public License Community Developers Platform on "Berlios" (cp2k.berlios.de)
More informationDehydrogenation of cyclohexaphenylene on Cu(111): possible reaction pathways from ab initio simulations
2nd CP2K Tutorial, ETHZ 7-11 February 2011 Materials Science & Technology Nanostructures and adsorption on metallic surfaces: Dehydrogenation of cyclohexaphenylene on Cu(111): possible reaction pathways
More informationSemi Empirical Force Fields and Their Limitations. Potential Energy Surface (PES)
Semi Empirical Force Fields and Their Limitations Ioan Kosztin Beckman Institute University of Illinois at Urbana-Champaign Potential Energy Surface (PES) Schrödinger equation: H T Ψ( r, = E Ψ( r, H =
More informationDay 1 : Introduction to AIMD and PIMD
Day 1 : Introduction to AIMD and PIMD Aug 2018 In today s exercise, we perform ab-initio molecular dynamics (AIMD) and path integral molecular dynamics (PIMD) using CP2K[?]. We will use the Zundel s cation
More informationPotentials, periodicity
Potentials, periodicity Lecture 2 1/23/18 1 Survey responses 2 Topic requests DFT (10), Molecular dynamics (7), Monte Carlo (5) Machine Learning (4), High-throughput, Databases (4) NEB, phonons, Non-equilibrium
More informationQM/MM for BIO applications
QM/MM for BIO applications February 9, 2011 I-Feng William Kuo kuo2@llnl.gov, P. O. Box 808, Livermore, CA 94551 This work performed under the auspices of the U.S. Department of Energy by under Contract
More informationDFT EXERCISES. FELIPE CERVANTES SODI January 2006
DFT EXERCISES FELIPE CERVANTES SODI January 2006 http://www.csanyi.net/wiki/space/dftexercises Dr. Gábor Csányi 1 Hydrogen atom Place a single H atom in the middle of a largish unit cell (start with a
More informationDihedral Angles. Homayoun Valafar. Department of Computer Science and Engineering, USC 02/03/10 CSCE 769
Dihedral Angles Homayoun Valafar Department of Computer Science and Engineering, USC The precise definition of a dihedral or torsion angle can be found in spatial geometry Angle between to planes Dihedral
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 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 informationNon-bonded interactions
speeding up the number-crunching continued Marcus Elstner and Tomáš Kubař December 3, 2013 why care? number of individual pair-wise interactions bonded interactions proportional to N: O(N) non-bonded interactions
More informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
More informationNon-bonded interactions
speeding up the number-crunching Marcus Elstner and Tomáš Kubař May 8, 2015 why care? key to understand biomolecular structure and function binding of a ligand efficiency of a reaction color of a chromophore
More informationA reciprocal space based method for treating long range interactions in ab initio and force-field-based calculations in clusters
A reciprocal space based method for treating long range interactions in ab initio and force-field-based calculations in clusters Glenn J. Martyna and Mark E. Tuckerman Citation: The Journal of Chemical
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 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 informationElectronic structure simulations of water solid interfaces
Electronic structure simulations of water solid interfaces Angelos Michaelides London Centre for Nanotechnology & Department of Chemistry, University College London www.chem.ucl.ac.uk/ice Main co-workers:
More informationFast Space Charge Calculations with a Multigrid Poisson Solver & Applications
Fast Space Charge Calculations with a Multigrid Poisson Solver & Applications Gisela Pöplau Ursula van Rienen Rostock University DESY, Hamburg, April 26, 2005 Bas van der Geer Marieke de Loos Pulsar Physics
More informationCharge Analysis: Atoms in Molecules
Daubechies Wavelets in Electronic Structure Calculation: BigDFT Code Tutorial CECAM - GRENOBLE : Atoms in Molecules Ali Sadeghi Basel University 21 November 2011 An output of electronic structure calculations
More informationComputational Molecular Biophysics. Computational Biophysics, GRS Jülich SS 2013
Computational Molecular Biophysics Computational Biophysics, GRS Jülich SS 2013 Computational Considerations Number of terms for different energy contributions / E el ~N ~N ~N ~N 2 Computational cost for
More informationNew Algorithms for Conventional and. Jing Kong ( 孔静 ) Q-Chem Inc. Pittsburgh, PA
New Algorithms for Conventional and Nondynamic DFT Jing Kong ( 孔静 ) Q-Chem Inc. Pittsburgh, PA XC Numerical Integration mrxc: multiresolution exchage-correlation Chunming Chang Nick Russ Phys. Rev. A.,
More informationLecture 12: Solvation Models: Molecular Mechanics Modeling of Hydration Effects
Statistical Thermodynamics Lecture 12: Solvation Models: Molecular Mechanics Modeling of Hydration Effects Dr. Ronald M. Levy ronlevy@temple.edu Bare Molecular Mechanics Atomistic Force Fields: torsion
More informationOpenDiscovery: Automated Docking of Ligands to Proteins and Molecular Simulation
OpenDiscovery: Automated Docking of Ligands to Proteins and Molecular Simulation Gareth Price Computational MiniProject OpenDiscovery Aims + Achievements Produce a high-throughput protocol to screen a
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationCP2K: Past, Present, Future. Jürg Hutter Department of Chemistry, University of Zurich
CP2K: Past, Present, Future Jürg Hutter Department of Chemistry, University of Zurich Outline Past History of CP2K Development of features Present Quickstep DFT code Post-HF methods (RPA, MP2) Libraries
More informationAmino Acids and Proteins at ZnO-water Interfaces in Molecular Dynamics Simulations: Electronic Supplementary Information
Amino Acids and Proteins at ZnO-water Interfaces in Molecular Dynamics Simulations: Electronic Supplementary Information Grzegorz Nawrocki and Marek Cieplak Institute of Physics, Polish Academy of Sciences,
More informationSupporting Information
alladium monophosphine d(h 3 ): does it really exist in solution? ietro Vidossich*, Gregori Ujaque,* and Agustí Lledós* Supporting Information Contents: 1. Schemes S1 and S2 2. Computational details 3.
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 informationOutline. Introduction: graphene. Adsorption on graphene: - Chemisorption - Physisorption. Summary
Outline Introduction: graphene Adsorption on graphene: - Chemisorption - Physisorption Summary 1 Electronic band structure: Electronic properties K Γ M v F = 10 6 ms -1 = c/300 massless Dirac particles!
More informationSupporting Online Materials: Nature of proton. transport in a water-filled carbon nanotube and in. liquid water
Supporting Online Materials: Nature of proton transport in a water-filled carbon nanotube and in liquid water Ji Chen, Xin-Zheng Li,,, Qianfan Zhang, Angelos Michaelides, and Enge Wang, ICQM and School
More informationThe Molecular Dynamics Method
H-bond energy (kcal/mol) - 4.0 The Molecular Dynamics Method Fibronectin III_1, a mechanical protein that glues cells together in wound healing and in preventing tumor metastasis 0 ATPase, a molecular
More informationUsing Molecular Dynamics to Compute Properties CHEM 430
Using Molecular Dynamics to Compute Properties CHEM 43 Heat Capacity and Energy Fluctuations Running an MD Simulation Equilibration Phase Before data-collection and results can be analyzed the system
More informationOptimizing GROMACS for parallel performance
Optimizing GROMACS for parallel performance Outline 1. Why optimize? Performance status quo 2. GROMACS as a black box. (PME) 3. How does GROMACS spend its time? (MPE) 4. What you can do What I want to
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 informationCOSMO-RS Theory. The Basics
Theory The Basics From µ to properties Property µ 1 µ 2 activity coefficient vapor pressure Infinite dilution Gas phase Pure compound Pure bulk compound Partition coefficient Phase 1 Phase 2 Liquid-liquid
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 informationSupporting Information
Supporting Information Ionic Liquid Designed for PEDOT:PSS Conductivity Enhancement Ambroise de Izarra, a,b,1 Seongjin Park, a,1 Jinhee Lee, a Yves Lansac, b,c, * Yun Hee Jang a, * a Department of Energy
More informationFragmentation methods
Fragmentation methods Scaling of QM Methods HF, DFT scale as N 4 MP2 scales as N 5 CC methods scale as N 7 What if we could freeze the value of N regardless of the size of the system? Then each method
More informationAlgebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes
Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes Elena Virnik, TU BERLIN Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov
More informationSupporting Information: Surface Polarons Reducing Overpotentials in. the Oxygen Evolution Reaction
Supporting Information: Surface Polarons Reducing Overpotentials in the Oxygen Evolution Reaction Patrick Gono Julia Wiktor Francesco Ambrosio and Alfredo Pasquarello Chaire de Simulation à l Echelle Atomique
More informationMaterial Surfaces, Grain Boundaries and Interfaces: Structure-Property Relationship Predictions
Material Surfaces, Grain Boundaries and Interfaces: Structure-Property Relationship Predictions Susan B. Sinnott Department of Materials Science and Engineering Penn State University September 16, 2016
More information3D Space Charge Routines: The Software Package MOEVE and FFT Compared
3D Space Charge Routines: The Software Package MOEVE and FFT Compared Gisela Pöplau DESY, Hamburg, December 4, 2007 Overview Algorithms for 3D space charge calculations Properties of FFT and iterative
More informationDensity Functional Theory
Density Functional Theory Iain Bethune EPCC ibethune@epcc.ed.ac.uk Overview Background Classical Atomistic Simulation Essential Quantum Mechanics DFT: Approximations and Theory DFT: Implementation using
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 informationNeighbor Tables Long-Range Potentials
Neighbor Tables Long-Range Potentials Today we learn how we can handle long range potentials. Neighbor tables Long-range potential Ewald sums MSE485/PHY466/CSE485 1 Periodic distances Minimum Image Convention:
More informationA Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane
A Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane Supplementary Information: Simulation Procedure and Physical Property Analysis Simulation Procedure The molecular
More informationvan der Waals corrections to approximative DFT functionals
van der Waals corrections to approximative DFT functionals Ari P Seitsonen Physikalisch Chemisches Institut, Universität Zürich February 8, 2011 van der Waals interactions van der Waals interactions important
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 informationCrystal Structure Prediction using CRYSTALG program
Crystal Structure Prediction using CRYSTALG program Yelena Arnautova Baker Laboratory of Chemistry and Chemical Biology, Cornell University Problem of crystal structure prediction: - theoretical importance
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 informationMolecular Dynamics Simulations
Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline Long-range Interactions Ewald Sum Fast Multipole Method Spherically Truncated Coulombic Potential Speeding up Calculations SPaSM
More informationAb Initio Molecular Dynamics
Ab Initio Molecular Dynamics Molecular Dynamics Classical equations of motion for the ions: F=ma NVE and other ensembles (NVT, NPT) Crucial for systems including liquids 'Ab Initio' or 'First Principles'
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 informationPaper No. 1: ORGANIC CHEMISTRY- I (Nature of Bonding and Stereochemistry)
Subject Chemistry Paper No and Title Paper 1: ORGANIC - I (Nature of Bonding Module No and Title Module Tag CHE_P1_M10 TABLE OF CONTENTS 1. Learning Outcomes 2. Introduction 3. Non-Covalent Interactions
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 informationElectronic State Calculation for Hydrogen Atom under Inhomogeneous Field
26 Special Issue Multiscale Simulations for Materials Research Report Electronic State Calculation for Hydrogen Atom under Inhomogeneous Field Shunsuke Yamakawa and Shi-aki Hyodo We discuss herein the
More informationSupporting information for Polymer interactions with Reduced Graphene Oxide: Van der Waals binding energies of Benzene on defected Graphene
Supporting information for Polymer interactions with Reduced Graphene Oxide: Van der Waals binding energies of Benzene on defected Graphene Mohamed Hassan, Michael Walter *,,, and Michael Moseler, Freiburg
More informationDispersion Interactions from the Exchange-Hole Dipole Moment
Dispersion Interactions from the Exchange-Hole Dipole Moment Erin R. Johnson and Alberto Otero-de-la-Roza Chemistry and Chemical Biology, University of California, Merced E. R. Johnson (UC Merced) Dispersion
More informationTools for QM studies of large systems
Tools for QM studies of large systems Automated, hessian-free saddle point search & characterization QM/MM implementation for zeolites Shaama Mallikarjun Sharada Advisors: Prof. Alexis T Bell, Prof. Martin
More informationFigure 1. Molecules geometries of 5021 and Each neutral group in CHARMM topology was grouped in dash circle.
Project I Chemistry 8021, Spring 2005/2/23 This document was turned in by a student as a homework paper. 1. Methods First, the cartesian coordinates of 5021 and 8021 molecules (Fig. 1) are generated, in
More informationElectronic properties of CdSe nanocrystals in the absence and presence of a dielectric medium
JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 11 15 MARCH 1999 Electronic properties of CdSe nanocrystals in the absence and presence of a dielectric medium Eran Rabani, Balázs Hetényi, and B. J. Berne
More informationAccurate van der Waals interactions from ground state electron density
Accurate van der Waals interactions from ground state electron density Alexandre Tkatchenko Theory Department, Fritz Haber Institut der MPG Berlin, Germany tkatchen@fhi berlin.mpg.de Haber Institute EXCITCM09,
More informationVan der Waals DFT Study of the Energetics of Alkali Metal Intercalation in Graphite
1 Van der Waals DFT Study of the Energetics of Alkali Metal Intercalation in Graphite Zhaohui Wang, Sverre M. Selbacħ and Tor Grande Department of Material Science and Engineering Trondheim, Nov. 29, 2013
More informationMolecular dynamics simulation of Aquaporin-1. 4 nm
Molecular dynamics simulation of Aquaporin-1 4 nm Molecular Dynamics Simulations Schrödinger equation i~@ t (r, R) =H (r, R) Born-Oppenheimer approximation H e e(r; R) =E e (R) e(r; R) Nucleic motion described
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 informationSoutheast University, Nanjing, China 2 Department of Applied Physics, Aalto University,
Supplementary Information to Solubility of Boron, Carbon and Nitrogen in Transition Metals: Getting Insight into Trends from First-Principles Calculations Xiaohui Hu, 1,2 Torbjörn Björkman 2,3, Harri Lipsanen
More informationNew Perspective on structure and bonding in water using XAS and XRS
New Perspective on structure and bonding in water using XAS and XRS Anders Nilsson Stanford Synchrotron Radiation Laboratory (SSRL) and Stockholm University, Sweden R. Ludwig Angew. Chem. 40, 1808 (2001)
More informationJet-FTIR-spectroscopy of pyrrole clusters
, Corey A. Rice and Martin A. Suhm Institute of Physical Chemistry, University Göttingen Tammannstr. 6, 37077 Göttingen, Germany 61 st Ohio State University International Symposium on Molecular Spectroscopy
More informationSTRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY
STRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY S. Hara, T. Kumagai, S. Izumi and S. Sakai Department of mechanical engineering, University of
More informationPotential Energy (hyper)surface
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 informationReport on TS-vdW Method, and Code Development, and. Results
Michael Pawley July, 21, 2010 Report on TS-vdW Method, and Code Development, and Results Introduction: Non-covalent Interactions have long been difficult to account for using Density Functional Theory(DFT),
More informationClassical potentials for metals
Classical potentials for metals About 80 % of all elements are metals. The crystal structures of the elements are distributed as follows: FCC 15 HCP 26 BCC 16 All other metals 13 So if we can describe
More informationWhy thermodynamics for materials?
Why thermodynamics for materials? Example p 2mkT T For = 300 K, = 1 atm ~ 10 8 site -1 s -1 p p Requires 10-12 atm to keep a clean surface clean; surface can also lose atoms Example Thermodynamic potentials
More informationPrediction of spectroscopic parameters for bio-organic and bio-inorganic intermediates in complex systems
Prediction of spectroscopic parameters for bio-organic and bio-inorganic intermediates in complex systems Erik Donovan Hedegård Department of Physics, Chemistry and Pharmacy University of Southern Denmark
More informationElectronic Supporting Information Topological design of porous organic molecules
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2017 Electronic Supporting Information Topological design of porous organic molecules Valentina Santolini,
More informationThe Projector Augmented Wave method
The Projector Augmented Wave method Advantages of PAW. The theory. Approximations. Convergence. 1 The PAW method is... What is PAW? A technique for doing DFT calculations efficiently and accurately. An
More informationApplications of Molecular Dynamics
June 4, 0 Molecular Modeling and Simulation Applications of Molecular Dynamics Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo Tohru
More informationIs there a future for quantum chemistry on supercomputers? Jürg Hutter Physical-Chemistry Institute, University of Zurich
Is there a future for quantum chemistry on supercomputers? Jürg Hutter Physical-Chemistry Institute, University of Zurich Chemistry Chemistry is the science of atomic matter, especially its chemical reactions,
More informationFrom Dynamics to Thermodynamics using Molecular Simulation
From Dynamics to Thermodynamics using Molecular Simulation David van der Spoel Computational Chemistry Physical models to describe molecules Software to evaluate models and do predictions - GROMACS Model
More informationLong-range interactions: P 3 M, MMMxD, ELC, MEMD and ICC
Long-range interactions: P 3 M, MMMxD, ELC, MEMD and ICC Axel Arnold Institute for Computational Physics Universität Stuttgart October 10, 2012 Long-range interactions gravity Coulomb interaction dipolar
More informationAccelerated Quantum Molecular Dynamics
Accelerated Quantum Molecular Dynamics Enrique Martinez, Christian Negre, Marc J. Cawkwell, Danny Perez, Arthur F. Voter and Anders M. N. Niklasson Outline Quantum MD Current approaches Challenges Extended
More informationModeling Materials. Continuum, Atomistic and Multiscale Techniques. gg CAMBRIDGE ^0 TADMOR ELLAD B. HHHHM. University of Minnesota, USA
HHHHM Modeling Materials Continuum, Atomistic and Multiscale Techniques ELLAD B. TADMOR University of Minnesota, USA RONALD E. MILLER Carleton University, Canada gg CAMBRIDGE ^0 UNIVERSITY PRESS Preface
More informationExploring the Free Energy Surface of Short Peptides by Using Metadynamics
John von Neumann Institute for Computing Exploring the Free Energy Surface of Short Peptides by Using Metadynamics C. Camilloni, A. De Simone published in From Computational Biophysics to Systems Biology
More informationLarge-scale real-space electronic structure calculations
Large-scale real-space electronic structure calculations YIP: Quasi-continuum reduction of field theories: A route to seamlessly bridge quantum and atomistic length-scales with continuum Grant no: FA9550-13-1-0113
More informationMany-body dispersion interactions between semiconducting wires
Many-body dispersion interactions between semiconducting wires Alston J. Misquitta TCM Group, Cavendish Laboratory University of Cambridge 26 th July, 2010 Collaborators James Spencer (Imperial College,
More informationTheoretical comparative study on hydrogen storage of BC 3 and carbon nanotubes
J. At. Mol. Sci. doi: 10.4208/jams.121011.011412a Vol. 3, No. 4, pp. 367-374 November 2012 Theoretical comparative study on hydrogen storage of BC 3 and carbon nanotubes Xiu-Ying Liu a,, Li-Ying Zhang
More informationAb-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University
Ab-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University Hands-on Tutorial Workshop on Ab-Initio Molecular Simulations, Fritz- Haber-Institut, Berlin, July 12-21,
More informationThe Force Field Toolkit (fftk)
Parameterizing Small Molecules Using: The Force Field Toolkit (fftk) Christopher G. Mayne, Emad Tajkhorshid Beckman Institute for Advanced Science and Technology University of Illinois, Urbana-Champaign
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 informationSupporting Information for The nature of the asymmetry in the hydrogen-bond networks of hexagonal ice and liquid water
Supporting Information for The nature of the asymmetry in the hydrogen-bond networks of hexagonal ice and liquid water Thomas D. Kühne and Rustam Z. Khaliullin, Institute of Physical Chemistry and Center
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 informationGromacs Workshop Spring CSC
Gromacs Workshop Spring 2007 @ CSC Erik Lindahl Center for Biomembrane Research Stockholm University, Sweden David van der Spoel Dept. Cell & Molecular Biology Uppsala University, Sweden Berk Hess Max-Planck-Institut
More informationCurrent address: Department of Chemistry, Hong Kong Baptist University, Kowloon Tong, Hong Kong,
Hydrolysis of Cisplatin - A Metadynamics Study Supporting Information Justin Kai-Chi Lau a and Bernd Ensing* b Department of Chemistry and Applied Bioscience, ETH Zurich, USI Campus, Computational Science,
More information