1 / 21

Molecular Modeling Using HPC and Gaussian: Density-Functional Theory and Noncovalent Interactions

Molecular Modeling Using HPC and Gaussian: Density-Functional Theory and Noncovalent Interactions. Gino A. DiLabio National Institute for Nanotechnology 11421 Saskatchewan Drive Edmonton, Alberta T5T 5A6 Gino.DiLabio@nrc.ca www.ualberta.ca/~gdilabio. Westgrid Seminar Series

inez
Download Presentation

Molecular Modeling Using HPC and Gaussian: Density-Functional Theory and Noncovalent Interactions

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Molecular Modeling Using HPC and Gaussian: Density-Functional Theory and Noncovalent Interactions Gino A. DiLabio National Institute for Nanotechnology 11421 Saskatchewan Drive Edmonton, Alberta T5T 5A6 Gino.DiLabio@nrc.ca www.ualberta.ca/~gdilabio Westgrid Seminar Series University of Alberta Feb. 5, 2014

  2. Geckos have evolved one of the most versatile and effective adhesives known. The mechanism of dry adhesion in the millions of setae on the toes of geckos has been the focus of scientific study for over a century. We provide the first direct experimental evidence for dry adhesion of gecko setae by noncovalent interactions (also called van der Waals forces) 12252–12256 PNAS September 17, 2002 vol. 99 no. 19 www.pnas.orgcgidoi10.1073pnas.192252799

  3. Why do we care about weak interactions? Formation of 1D and 2D organic lines on Si Nature, 2005, 435, 658-661.

  4. Why do we care about weak interactions? Heavy oil/bitumen upgrading Mackie and DiLabio, Energy and Fuels, 2010, 24, 6468-6475.

  5. Why do we care about weak interactions? Organic electronic material J. Phys. Chem. C, 2010, 114. 10952.

  6. How do we model these systems? • Large systems necessitate the use of density-functional theory (DFT) • B3LYP was one of the first hybrid DFT that was implemented in most computational chemistry programs. It has since found general acceptance and useand works well for thermochemistry. • However, most conventional DFT methods, including B3LYP,cannot accurately treat weak, non-covalently bonded systems – specifically “dispersion” interactions.

  7. How bad are conventional DFT methods for weak interactions? Chem. Phys. Lett. 2006, 419, 333-339

  8. Percent Error in Binding Energy vs. Interaction

  9. Improving the performance of DFT methods: Dispersion-corrected DFT in Gaussian • DFT-D/D3: Add-on to the DFT energy an empirical van der Waals term. Developed for use with many functionals. (Grimme, J. Chem. Phys. 2010, 132, 154104) • M06-2X: A DFT method parameterized to reproduce dispersion binding, among other properties. (Zhao and Truhlar, Theor. Chim. Acta, 2008, 120, 215.) • DCP: Dispersion correcting potentials correct the erroneous noncovalent behaviour of a small number of DFT methods. (Torres and DiLabio, J. Phys. Chem. Lett. 2012,3, 1738) – Compute Canada RAC supported work. Review: J. Phys. Org. Chem. 2009, 22, 1127-1135. See also: 2014 version of Reviews in Computational Chemistry

  10. How to incorporate dispersion corrections into your Gaussian DFT calculations Grimme “D3” approach Truhlar M06-2X approach #B3LYP Gen EmpiricalDispersion=GD3BJ SCF=(Conver=6) Water Dimer 0 1 O -0.702196054 -0.056060256 0.009942262 H -1.022193224 0.846775782 -0.011488714 H 0.257521062 0.042121496 0.005218999 O 2.220871067 0.026716792 0.000620476 H 2.597492682 -0.411663274 0.766744858 H 2.593135384 -0.449496183 -0.744782026 H O 0 6-31+G(2d,2p) **** #M062X Gen SCF=(Conver=6) Water Dimer 0 1 O -0.702196054 -0.056060256 0.009942262 H -1.022193224 0.846775782 -0.011488714 H 0.257521062 0.042121496 0.005218999 O 2.220871067 0.026716792 0.000620476 H 2.597492682 -0.411663274 0.766744858 H 2.593135384 -0.449496183 -0.744782026 H O 0 6-31+G(2d,2p) ****

  11. #B3LYP Gen Pseudo=Read SCF=(Conver=6) Water Dimer 0 1 O -0.702196054 -0.056060256 0.009942262 H -1.022193224 0.846775782 -0.011488714 H 0.257521062 0.042121496 0.005218999 O 2.220871067 0.026716792 0.000620476 H 2.597492682 -0.411663274 0.766744858 H 2.593135384 -0.449496183 -0.744782026 H O 0 6-31+G(2d,2p) **** H 0 H 1 0 P an up 3 2 0.120883601 0.000231333 2 0.044528578 -0.000070677 2 0.005658790 -0.000000451 S-P 1 2 0.174740501 -0.000049845 O 0 O 3 0 F an up 3 2 0.192168931 0.000242019 2 0.166560549 0.003000000 2 0.016734867 -0.000000131 S-F 1 2 0.039337457 0.000016041 P-F 1 2 0.123780982 -0.000129441 D-F 1 2 0.061418151 -0.003520223 How to use DCPs in Gaussian: Conventional Input DCPs Input generating utility for Gaussian input files at: www.ualberta.ca/~gdilabio

  12. How DCPs work • Based on (old) ECP technology: • e.g. Iodine = 1s22s22p6…5s24d105p5 [Potential] 5s24d105p5 • In the case of DCPs, the potentials don’t replace any core electrons but instead are constructed such that the long-range behaviour of a DFT is corrected.

  13. How DCPs are generated Compute Canada RAC work #B3LYP Gen Pseudo=Read SCF=(Conver=6) Water Dimer 0 1 O -0.702196054 -0.056060256 0.009942262 H -1.022193224 0.846775782 -0.011488714 H 0.257521062 0.042121496 0.005218999 O 2.220871067 0.026716792 0.000620476 H 2.597492682 -0.411663274 0.766744858 H 2.593135384 -0.449496183 -0.744782026 H O 0 6-31+G(2d,2p) **** H 0 H 1 0 P an up 3 2 0.120883601 0.000231333 2 0.044528578 -0.000070677 2 0.005658790 -0.000000451 S-P 1 2 0.174740501 -0.000049845 O 0 O 3 0 F an up 3 2 0.192168931 0.000242019 2 0.166560549 0.003000000 2 0.016734867 -0.000000131 S-F 1 2 0.039337457 0.000016041 P-F 1 2 0.123780982 -0.000129441 D-F 1 2 0.061418151 -0.003520223 H 0 H 1 0 P an up 3 2 0.120883601 0.000231333 2 0.044528578 -0.000070677 2 0.005658790 -0.000000451 S-P 1 2 0.174740501 -0.000049845

  14. How DCPs are generated Compute Canada RAC work • Steps: • Develop an initial guess for a DCP for an atom/DFT method/basis set combination. • Select the fitting data to which DCPs will be optimized. • Submit optimization script to Grex. The script: • Builds Gaussian input files based on the information in 1&2. • Builds queue files for Gaussian runs • Submits queue files to the queue • Monitors status of runs • Extracts energies of noncovalently bonded systems and their monomers • Adjusts values of Ci and ζi according to some optimization scheme • Select next atom and repeat from step 2.

  15. Past Successes B3LYP-DCP versus other dispersion-corrected DFT methods Performance on the S66 benchmark set: Occurrences in %error in BE

  16. Past Successes B3LYP-DCP versus other dispersion-corrected DFT methods Performance on the S66 benchmark set: Occurrences in %error in BE

  17. Past Successes B3LYP-DCP versus other dispersion-corrected DFT methods Performance on the S66 benchmark set: Occurrences in %error in BE

  18. Past Successes B3LYP-DCP versus other dispersion-corrected DFT methods Performance on the S66 benchmark set: Occurrences in %error in BE

  19. Past Successes More than just noncovalent interactions. DCPs can improve thermochemical properties: LC-ωPBE/6-31+G(2d,2p) Bond Dissociation Mean Absolute Error (kcal/mol)

  20. Future Efforts on DCPs “Low cost” Density-functional theory methods + various computational chemistry/physics codes. Thermochemistry, kinetics and noncovalent interactions. Develop a deeper understanding of how DCPs work. Challenges Optimization scripts are complicated and platform dependent. Conventional resource allocation scheme is not ideal for this work.

  21. Thanks to: Dr. Edmanuel Torres - NINT Dr. Iain D. Mackie - NINT Prof. Erin R. Johnson – UC Merced

More Related