Saturday, March 23, 2013

DFT-D, DFT-D2 and DFT-D3: Density functional dispersion correction

Dispersion-corrected DFT, such as B3LYP-D3, are not new functionals but a mix of conventional functionals and an add-on energy term.  For example, B3LYP-D3 denotes a calculation with the usual B3LYP functional plus a D3 dispersion correction energy term.  The dispersion correction energy term is a relatively simply function of interatomic distances and contain adjustable parameters that are fitted to conformational and interaction energies computed using CCSD(T)/CBS.  The fitting is done for a given functional. DFT-D and DFT-D2 energy corrections consider all pairs of atoms while DFT-D3 also consider triplets of atoms to account for three-body effects.

Because the dispersion correction is an add-on term it does not directly alter the wavefunction or any other molecular property.  However, geometry optimizations with dispersion correction will lead to a different geometry than without because the dispersion correction contributes to the forces acting on the atoms.

Dispersion corrections can lead to significant improvements in accuracy and the computational cost associated with dispersion corrections are negligible, so if your favorite code supports dispersion corrections for your functional of choice there is little reason not to use it.

Similar (but not identical) dispersion corrections have also been developed for semi-empirical methods such as PM6, usually in combinations with analogous add-on energy terms to improve hydrogen bonding, e.g. PM6-DH+ and PM6-DH2.

Further reading:
A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu
Dispersion corrections and bio-molecular structure and reactivity

Creative Commons License
This work is licensed under a Creative Commons Attribution 3.0 

Thursday, February 28, 2013

The Molecule Calculator v.1.1


+Jimmy Charnley Kromann has put a new version of the Molecule Calculator online.  The main new features are:

1. JSmol instead of Jmol.  This means that MolCalc now works on iPhones and iPads.  The only thing that doesn't seem to work there is measuring distances and angles.  This also means that the load button had disappeared.

2. New starting points for molecule building: the default starting structure is still methane, but this can now now be changed to benzene or the water dimer.  It is also possible to search for other structures using common names or SMILES.

3. One can now compute other thermodynamic properties in addition to the heat of formation.

The guts of the calculations are still the same as for version 1.0 except that the MMFF force field is used for pre-optimization instead of the UFF force field.

As always the source code is available on Github

The development of MolCalc is supported by the University of Copenhagen through the Education at its Best initiative (Den gode uddannelse).

Creative Commons License
This work is licensed under a Creative Commons Attribution 3.0 Unported License.  

Sunday, February 3, 2013

Getting started with JSmol

Readers of this blog will know I am quite fond of Jmol. But Jmol uses Java, which has some securities issues and doesn't work on iOS.  Luckily, Jmol has been converted to JavaScript to create a program called JSmol.

JSmol accepts all the same commands as Jmol, but setting up a webpage with JSmol is a little bit different.  Here is a very simple page to help get you started and the HTML code and be found below.

The location of the jsmoltest.html is in my home directory while the jsmol code is in a folder called "jsmol" also sitting in my hope directory.  The jsmol folder is created when you unzip the JSmol zip file you can download here.

The JSmol distribution also includes a simple example of having two JSmol windows on one webpage.  You can find that file here.

Creative Commons License
This work is licensed under a Creative Commons Attribution 3.0