Monthly Archives: October 2020

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.

Mario Molina, Nobel Laureate. Rest In Peace


Prof. Mario Molina was awarded the Nobel Prize in Chemistry in 1995, the same year I started my chemistry education at the chemistry school from the National Autonomous University of Mexico, UNAM, the same school from where he got his undergraduate diploma. To be a chemistry student in the late nineties in Mexico had Prof. Molina as a sort of mythical reference, something to aspire to, a role model, the sort of representation the Latinx and other underrepresented communities still require and seldom get.

I saw him several times at UNAM, where he’d pack any auditorium almost once a year to talk about various research topics, but I remember distinctly the first time I sort of interacted with him. It was 1997 and I attended my first congress, the 5th North America Chemistry Congress. Minutes before the official inauguration which he was supposed to preside, I caught a glimpse of him in the hallways near the main conference room. Being only 19 years old, I thought it’d be a good idea to chase him, ask for his autograph and a picture. He was kind enough not to brush me off and took just a minute to shake my hand, sign my book of abstracts, and get his picture taken with me. But cameras back then relied on the user to place a roll of film correctly. I did not; so the picture, although it happened, it doesn’t exist. Because of this and other anecdotes, that congress cemented my love for chemistry. I never asked for a second picture in the few subsequent occasions I had the pleasure to hear him talk.

Prof. Molina was an advocate of green and sustainable sources of energies. His work predicted the existence of a hole in the ozone layer and his struggle brought change into the banning of CFCs and other substances which interfere with the replenishment of ozone in the sub-stratosphere. Today, his legacy remains but also do his pending battles in the quest for new policies that favor the use of green alternative forms of energy. May he rest in peace and may we continue his example.

%d bloggers like this: