computational quantum chemistry tools
TRANSCRIPT
![Page 1: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/1.jpg)
Computational Quantum Chemistry Tools
Chun-Min ChangResearch Consultant
Institute for Cyber-Enabled Research
Download this presentation: https://wiki.hpcc.msu.edu/display/TEAC/iCER+Training+Home
![Page 2: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/2.jpg)
Outline
• Overview of Quantum Mechanical simulation & tools
• Build molecules & run calcuations in interactive mode
– gaussview and gaussian tutorial for computational chemistry
• Calculate electronic structures with more resource by
job submission in batch mode
• Checkpoint run your jobs with more time in HPCC
• Summary & Questions
![Page 3: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/3.jpg)
Overview of Quantum
Mechanical Simulation & Tools
![Page 4: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/4.jpg)
Electron behaving as waves
The wave-like nature of the electron allows it to pass
through two parallel slits simultaneously, rather than just one
slit as would be the case for a classical particle.
In quantum mechanics, the wave-like property of one
particle can be described mathematically as a complex-
valued function, the wave function, commonly denoted by
the Greek letter psi (ψ). When the absolute value of this
function is squared, it gives the probability that a particle will
be observed near a location—a probability density ( )ρ x
![Page 5: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/5.jpg)
Hartree–Fock method
In atomic units, and with r denoting electronic and R denoting nucleardegrees of freedom, the electronic SchrÄodinger equation is
with a wavefunction of the general form which is known as a HartreeProduct:
The generalization to N electrons is
Slater determinant
Electrons are indistinguishable
![Page 6: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/6.jpg)
Hartree–Fock method
One-electron operator h as follows
Coulomb operator as Exchange operator
Fock operator, as
Introducing a basis set transforms
Hartree-Fock-Roothaan equations can be written in matrix form
Eigenvalues
![Page 7: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/7.jpg)
DFT is the most popular and versatile quantum mechanical
modeling method used to investigate the electronic structure
(principally the ground state) in physics and chemistry.
Despite recent improvements, there are still difficulties in using
density functional theory to properly describe intermolecular
interactions, especially van der Waals forces (dispersion); transition
states, some other strongly correlated systems; and in calculations
of the band gap and ferromagnetism in semiconductors.
Density Functional Theory
![Page 8: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/8.jpg)
In density functional theory, the total energy of a system is expressed as a functional of the charge density as
where Ts is the Kohn–Sham kinetic energy which is expressed in terms of the Kohn–Sham orbitals as
vext is the external Potential acting on the interacting system, VH is the Hartree (or Coulomb) energy,
DFT Calculation
![Page 9: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/9.jpg)
to yield the Kohn–Sham potential as
where the last term is the exchange-correlation potential:
and Exc is the exchange-correlation energy. The Kohn–Sham equations are found by varying the total energy expression with respect to a set of orbitals:
DFT Method
![Page 10: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/10.jpg)
How we treat electrons in different simulation systems?
Electron wave function described by a basis set:
*Local Atomic Orbital Basis Set (with Gaussian functions):
Non-periodic molecular systems
(periodic system? Looking for fast method DFTB, QEM)
*Plane Wave Basis Set : Periodic systems (solid state)
All-electron:
system with small atoms (need to use LAO Basis Set)
Pseudo-potentials:
system with heavy atoms. (mostly use Plane Wave Basis Set)
Electrons are divided in two groups: valence electrons and inner
core electrons.
![Page 11: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/11.jpg)
• Electrons are treated as charge density distributions (probability) within certain ranges or PBC.
• Consider electron correlation interactions.
• Diagonalization of Hamiltonian and self consistent field (SCF) iterationsto derive the stationary energy state.
• Accurate but slow, good for small molecular systems, charge transfer, bond forming and breaking.
• Electrons and nuclei are usually considered as point charges or multipole moments.
• Electron correlation interactions are replaced by force field.
• No matrix diagonalization or SCF.
• Less accurate but fast, for big molecules, large systems & statistical behaviors.
Quantum MechanicsMolecular Mechanics
Some Comparison between QM and MM Calculation
![Page 12: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/12.jpg)
Quantum Mechanical Simulation Packages
Package License† Basis Periodic‡ Mol. mech. Semi-emp. HF Post-HF DFT
CP2K GPLHybrid GTO /PW
3d Yes Yes Yes Yes Yes
CRYSTAL Commercial GTO Any Yes No Yes Yes10 Yes
GAMESS Academic GTO No Yes2 Yes Yes Yes Yes
Gaussian Commercial GTO Any Yes Yes Yes Yes Yes
MOLPRO Commercial GTO No No No Yes Yes Yes
NWChem ECL v2 GTO, PW Yes(PW) No(GTO) Yes No Yes Yes Yes
PSI GPL GTO No No No Yes Yes Yes
Q-Chem Commercial GTO No Yes Yes Yes Yes Yes
Quantum ESPRESSO
GPL PW 3d Yes No Yes No Yes
Spartan Commercial GTO No Yes Yes Yes Yes Yes
VASP Commercial PW 3d Yes No Yes Yes Yes
![Page 13: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/13.jpg)
Program Developer(s) License Platforms Info
Avogadro Avogadro project team GPL Linux, Mac OS X, Windows
3D molecule editor and visualizer
ChemDraw CambridgeSoft proprietary Mac OS X, Windows Chemical structure and reaction editor
Discovery Studio Accelrys proprietary Windows
freeware version available; includes name2structure and structure2name, InChInaming, and canonical SMILES
GaussView Gaussian Commercial Linux, Mac OS X, Windows
IQmol ANU GPL Linux, Mac OS X, Windows
MOE Chemical Computing Group Commercial Linux, Mac OS X,
Windows
SPARTAN Wavefunction, Inc. proprietaryLinux, Mac OS X, Windows
VMD UIUC proprietary Linux, Mac OS X, Windows
Molecule Editor Software
![Page 14: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/14.jpg)
Building Molecular Geometry
and Computing Electronic
Structure Interactively
![Page 15: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/15.jpg)
An electronic structure package capable of predicting many properties of atoms, molecules, and reactive systems, e.g.
•single point energy and properties (electron density, dipole moment, …)
•geometry optimization
•frequency
•reaction path following
utilizing ab initio, density functional theory, semi-empirical, molecular mechanics, and hybrid methods.
Calculations with Gaussian
![Page 16: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/16.jpg)
Calculations with Gaussian
• semi-empirical: AM1, PM3, MNDO, …
• density functional theory: B3LYP, MPW1PW91, …
• ab initio: HF, MP2, CCSD, CCSD(T), …
• hybrid: G2, G3, …
Levels of Theory Available:
• Pople-type: 3-21G, 6-31G, 6-311G(d,p), …
• Dunning: cc-pVDZ, aug-cc-pVTZ, …
• Huzinaga and Others: MIDIX, …
• User-defined or General Basis Set
Basis Sets Available:
![Page 17: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/17.jpg)
GaussView --- graphical interface for Gaussian 09 ---
• build molecules or reactive systems
• setup Gaussian 09 input files
• graphically examine results
Visualization
Electronic Structure Calculations
Molecule Building with GaussView
![Page 18: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/18.jpg)
• Instead of typing all the coordinates, theory, basis set, etc., we can use GaussView.
• The calculation is specified by pointing and clicking to build the molecule, and using pull-down menus to select the calculation type, level of theory and basis set.
• GaussView generates the Gaussian input file, and can run Gaussian without ever returning to the Unix prompt.
• GaussView can also be used to read Gaussian output files and visualize the results.
Building with GaussView
![Page 19: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/19.jpg)
• Open a terminal and Connect to HPCC:
$ ssh –X <MSU Net ID>@hpcc.msu.edu
and ssh to a development node (such as dev-gfx10).
• For a fast graphical interface, use remote desktop toconnect rdp.hpcc.msu.edu and ssh to a developmentnode.
• We use GaussView to create molecular geometry.
• In the command line, please do$ module load Gaussian/g09
$ gview
Getting Started with GaussView
![Page 20: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/20.jpg)
Builder Window
Viewer Window
Start GaussView
![Page 21: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/21.jpg)
Use Your Mouse
![Page 22: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/22.jpg)
Steps to Building
• Choose an atom or fragment.
• Choose the location of the fragment that will attach to your molecule in the builder window.
• Select/add the fragment to your molecule in the viewer window.
• Repeat until molecular building is done.
• Click the Clean button (the broom icon ). This option will use Molecular Mechanics calculations to optimize the geometry (very fast).
![Page 23: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/23.jpg)
The atom labeled “Hot” is is where the fragment will attach to the system you are building.
Molecule is put together here.
Fragments are selected here.
Elements are selected here
![Page 24: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/24.jpg)
Submitting jobs through GaussView
Save Input File
Choose Method & Basis Set
Choose Job Type
Select Calculate
Submit
![Page 25: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/25.jpg)
Run your Gaussian job
• If you haven’t saved your input file yet, Gaussviewwill Prompt you to save your input file.
• You will name the input file, and then submit the calculation.
• When Gaussian is finished running, you will receive a message in Gaussview.
• BE CAREFUL how long you run interactively.
![Page 26: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/26.jpg)
How to view your output
• Your submitted calculation will run in the background.
• When it is complete, Gaussview will inform you and ask you if you wish to view an output file.
• From the list of files, you can pick your output (.log or .chk file).
• To open an output file again, click File on the menu then click Open.
![Page 27: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/27.jpg)
What can be visualized
• Geometry• Vibrations• Orbitals• Electron density• Electrostatic potential
Main results menu
![Page 28: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/28.jpg)
Vibration Frequency & IR Spectrum
Selecting Vibrations
from the main-window
Results menu will
calculate and then
display the vibration
frequency, plus a
calculated IR spectrum.
![Page 29: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/29.jpg)
Display Molecular Orbitals
To view molecular orbitals, choose Surfaces/Contours from the main-window Resultsmenu. Begin by generating "cube" files for the HOMO and LUMO to display orbitals,electron densities, electrostatic potentials, etc.
![Page 30: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/30.jpg)
Submitting Calculations by Command Line
Interactively use command:g09 < input_file > output_file
Notes:• Input and output files can use absolute or relative pathnames.
• Gaussian temporary files can be redirected by setting the GAUSS_SCRDIR environment variable, i.e.
GAUSS_SCRDIR=/mnt/scratch/$USER/GAUSS_SCRDIR
• Schedule more time to run your jobs interactively:https://wiki.hpcc.msu.edu/display/hpccdocs/Scheduling+Interactive+Jobs
![Page 31: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/31.jpg)
Quantum Mechanical Simulations
with More Resource Usage by Job
Submission in Batch Mode
![Page 32: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/32.jpg)
More Resource Request
• In Gaussian input file, you can specify more cores and memory for your calculation.
• Create a job script with Gaussian command to execute your QM calcuations.
• Submit the calculation jobs with wall-time less than 168 hours (7 days).
%NProcShared=8%Mem=15GB%Chk=mychk#P Freq …
input continues…
The Gaussian input file request:
• 8 cores
• 15GB RAM memory
to do the calculation
![Page 33: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/33.jpg)
Job Sheduling Script
https://wiki.hpcc.msu.edu/display/hpccdocs/Scheduling+Jobs
How to write a job script to schedule a QM Simulation job?
Job Scriptsqsub options
Scheduling TipsWhen will my job start?Why was my job killed?
Queuing Policies
![Page 34: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/34.jpg)
A Simple Job Script
#!/bin/bash --login
###### Give the job a name
#PBS -N GaussianJob
###### Number of nodes and the number of processors per node
#PBS -l nodes=1:ppn=8
###### Memory needed by the job
#PBS -l mem=18gb
###### Time job will take to execute (HH:MM:SS format)
#PBS -l walltime=72:00:00
###### The location for Output file
#PBS -o JobOut
###### Make output and error files the same file
#PBS -j oe
###### Setup the clusters available for the job
#PBS -l feature='intel14|intel11|intel10|gfx10'
###### Send an email when a job is aborted, begins or ends
#PBS -m n
Resource Request for HPCC
Tip:Request more memory in job script than the request from Gaussian input file.
![Page 35: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/35.jpg)
A Simple Job Script
InputFile= <Gaussian input file name>
OutputFile= <Gaussian output file name>
cd ${PBS_O_WORKDIR}
module load Gaussian/g09
#GAUSS_SCRDIR=/mnt/scratch/$USER/GAUSS_SCRDIR
#mkdir -p ${GAUSS_SCRDIR}
g09 < ${InputFile} > ${OutputFile}
# Print out final statistics about resource uses before job
echo; echo "Job is finished. Print out resource usage:“
qstat -f ${PBS_JOBID}
echo; echo 'Exit!!!'
Commands to
be executed !
![Page 36: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/36.jpg)
Use Scratch Directory
• The I/O in home and research directory is not fast. For fast job running (especially for request of many cores), we need to use Scratch directory.
• Scratch space can do very fast I/O. Move all of your Gaussian input files and Job script there.
• You may set up Gaussian Scratch directory to a proper place by setting up the variable GAUSS_SCRDIR in your job script and making the directory: mkdir –p ${GAUSS_SCRDIR}
• After jobs are done, copy your files back to home or research directory. Files on scratch will be purged after 45 days.
![Page 37: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/37.jpg)
Checkpoint Run your Jobs with
More Time in HPCC
What if you need calculation time more than 7 days?
![Page 38: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/38.jpg)
Using Checkpoint Run
• HPCC maintains buy-in nodes, which are able to run non-priority (non buy-in) jobs that are less than four hours in walltime. If you submit a job which requests less than four hours of walltime, it is likely that the short job will start fairly quickly.
• In the input file, specify restart mode for your calculation (checkpoint run).
• Submit the calculation jobs with wall-time less than 4 hours.
• If the job does not finish before the job ends, submit another job again and restart from the place where last calculation ends.
• Keep the job submissions until the calculation finishes.
![Page 39: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/39.jpg)
Restart Previous Gaussian Job
For example, in a frequency
calculation you would include
following commands:
The following input file
will restart a previous
job:
%RWF=myrwf%NoSave%Chk=mychk#P Freq …
input continues…
%RWF=myrwf%NoSave%Chk=mychk#P Restart
• The checkpoint file can be placed in the directory of the output file, which might have only a moderate amount of free space. However, the read-write file for a job large enough to be worth restarting should be on a large, local scratch file system.
• The file (such as myrwf from %RWF) before %NoSave will be deleted if the job finishes normally
![Page 40: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/40.jpg)
Job Script for Continuous Submission
• With the Checkpoint Run, we can restart a Gaussian calculation where it is not finished last time.
• It is better if we can have a job script to keep submitting restarting jobs automatically if the calculation is not finished.
• The job script contains two major processes:
A foreground process to execute Gaussian calculation
A background process in sleep until the end of wall-time to submit another.
• If the foreground process is not finished at the end of wall-time the background process will submit another and kill the running job.
• See example: /mnt/ls15/scratch/users/changc81/QCW/example.pbs
![Page 41: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/41.jpg)
Shell Script to Create a PBS script
• You could also use a Shell Script to Create a PBS script for Gaussian checkpoint run according to your Gaussian input file.
• See example: /mnt/ls15/scratch/users/changc81/QCW/qsubg09
• To use the Shell Script, use the command with a input of the Gaussian input file name:
./qsubg09 <Gaussian input file name>
• Some possible selections:
-t Wall-Time (HH:MM:SS format)
-o Job output file name
-N Job Name
-q PBS script file name
![Page 42: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/42.jpg)
How to Specify Job Resource & Check Job Status
• Check HPCC resources of computing nodes:• https://wiki.hpcc.msu.edu/display/hpccdocs/Documentation+and+User+Manual#DocumentationandUserManual-SystemConfiguration
• Check usage of HPCC resource now by
module load powertools
and Use commands, such as … clusterstate, node_status
• According to the HPCC resource, specify your job resource.
• Check your job status use commands (load powertools first):
qs, sq, qsq
qstat -f <PBS_JOBID>, checkjob –v <PBS JOBID>
• More from powertools: powertools, bigbrother, memorymonitor, …
![Page 43: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/43.jpg)
Buy-In Opportunities
• We will maintain your computers for you• Researchers get exclusive use of their nodes within 4 hours
of submitting a job• Buy-in jobs will automatically overflow into the general
resources.• 2016 Buy-In is coming.
![Page 44: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/44.jpg)
Getting Help
• Documentation and User Manual – wiki.hpcc.msu.edu• Contact HPCC and iCER Staff for:
– Reporting System Problems– HPC Program writing/debugging Consultation– Help with HPC grant writing– System Requests– Other General Questions
• Primary form of contact - http://contact.icer.msu.edu/• HPCC Request tracking system – rt.hpcc.msu.edu• HPCC Phone – (517) 353-9309• HPCC Office – 1400 PBS• Open Office Hours – 1pm to 2pm Monday and Thursdays (BPS 1440)
![Page 45: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/45.jpg)
Q & A
![Page 46: Computational Quantum Chemistry Tools](https://reader031.vdocument.in/reader031/viewer/2022022013/6210da8c4194c723637b530e/html5/thumbnails/46.jpg)
Thanks!