Blog Archives

Density Keyword in Excited State Calculations with Gaussian

I have written about extracting information from excited state calculations but an important consideration when analyzing the results is the proper use of the keyword density.

This keyword let’s Gaussian know which density is to be used in calculating some results. An important property to be calculated when dealing with excited states is the change in dipole moment between the ground state and any given state. The Transition Dipole Moment is an important quantity that allows us to predict whether any given electronic transition will be allowed or not. A change in the dipole moment (i.e. non-zero) of a molecule during an electronic transition helps us characterize said transition.

Say you perform a TD-DFT calculation without the density keyword, the default will provide results on the lowest excited state from all the requested states, which may or may not be the state of interest to the transition of interest; you may be interested in the dipole moment of all your excited states.

Three separate calculations would be required to calculate the change of dipole moment upon an electronic transition:

1) A regular DFT for the ground state as a reference
2) TD-DFT, to calculate the electronic transitions; request as many states as you need/want, analyze it and from there you can see which transition is the most important.
3) Request the density of the Nth state of interest to be recovered from the checkpoint file with the following route section:

# TD(Read,Root=N) LOT Density=Current Guess=Read Geom=AllCheck

replace N for the Nth state which caught your eye in step number 2) and LOT for the Level of Theory you’ve been using in the previous steps. That should give you the dipole moment for the structure of the Nth excited state and you can compare it with the one in the ground state calculated in 1). Again, if density=current is not used, only properties of N=1 will be printed.

XVI Mexican Meeting on Phys.Chem.

A yearly tradition of this Comp.Chem. lab and many others throughout our nation is to attend the Mexican Meeting on Theoretical Physical Chemistry to share news, progress and also a few drinks and laughs. This year the RMFQT was held in Puebla and although unfortunately I was not able to attend this lab was proudly represented by its current members. Gustavo Mondragón gave a talk about his progress on his photosynthesis research linking to the previous work of María Eugenia Sandoval already presented in previous editions; kudos to Gustavo for performing remarkably and thanks to all those who gave us their valuable feedback and criticism. Also, five posters were presented successfully, I can only thank the entire team for representing our laboratory in such an admirable way, and a special mention to the junior members, I hope this was the first of many scientific events they attend and may you deeply enjoy each one of them.

Among the invited speakers, the RMFQT had the honor to welcome Prof. John Perdew (yes, the P in PBE); the team took the opportunity of getting a lovely picture with him.

Here is the official presentation of the newest members of our group:

Alejandra Barrera (hyperpolarizabilty calculations on hypothetical poly-calyx[n]arenes for the search of NLO materials)


Fernando Uribe (Interaction energy calculations for non-canonical nucleotides)


Juan Guzmán (Reaction mechanisms calculations for catalyzed organic reactions)


We thank the organizing committee for giving us the opportunity to actively participate in this edition of the RMFQT, we eagerly await for next year as every year.


Grimme’s Dispersion DFT-D3 in Gaussian #CompChem

I was just asked if it is possible to perform DFT-D3 calculations in Gaussian and my first answer was to use the following  keyword:


which is available in G16 and G09 only in revision D, apparently.

There are also some overlays that can be used to invoke the use dispersion in various scenarios:

IOp(3/74=x) Exchange and Correlation Potentials





DSD-PBEP86 (double hybrid, DFT-D3).


B2PLYP-D3 (double hybrid, DFT-D3).

B97-D (DFT-D3).

IOp(3/76=x) Mixing of HF and DFT.

-33 PW6B95 and PW6B95-D3 coefficients.

IOp(3/124=x) Empirical dispersion term.




Force dispersion type 3 (Grimme DFT-D3).

Force dispersion type 4 (Grimme DFT-D3(BJ)).

Force dispersion type 5 (Grimme D3, PM7 version).


The D3 correction method of Grimme defines the van der Waals energy like:

$\displaystyle E_{\rm disp} = -\frac{1}{2} \sum_{i=1}^{N_{at}} \sum_{j=1}^{N_{at...
...{6ij}} {r_{ij,{L}}^6} +f_{d,8}(r_{ij,L})\,\frac{C_{8ij}} {r_{ij,L}^8} \right ),$

where coefficients $ C_{6ij}$ are adjusted depending on the geometry of atoms i and j. The damping D3 function for is:

$\displaystyle f_{d,n}(r_{ij}) = \frac{s_n}{1+6(r_{ij}/(s_{R,n}R_{0ij}))^{-\alpha_{n}}},$

where the values of s are adjustable parameters fit for the exchange-correlation functionals used in each calculation.

%d bloggers like this: