Monthly Archives: November 2015

XIV Mexican Reunion on Theoretical Physical Chemistry


Each year the Mexican community who works in the realm of computational and theoretical chemistry gathers to share the most recent work done around our country. This year, I tried to live Tweet the event and although I failed miserably in doing so -as well as in convincing others to join me- I’m trying to put together the things that caught my attention. I also tried to Storify it but I cannot embed the result here in WordPress.

María Eugenia “Maru” Sandoval just came back from a short stay in Spain where she worked with Prof. David Casanova on Singlet Fission theory applied to her work on photosynthesis. Here work was presented as a poster although we would have preferred she gave a talk.

  

Also, Guillermo “Memo” Caballero presented his recent developments in reaction mechanisms.   

Below there is a list of Tweets from the conference. If you are interested in any of these items please contact me for further information, or just, you know, google the people mentioned in each Tweet, it shouldn’t be too hard.

//platform.twitter.com/widgets.js

//platform.twitter.com/widgets.js

Of course when you have a large meeting with so many people working with and on Density Functional Theory (DFT) you know that Perdew’s suggested ‘Jacob’s ladder’ of functional quality for chemical accuracy is bound to show up a few times.
I actually wrote a post that gravitates around this issue but using figurative painting as an analogue. You may find it here

//platform.twitter.com/widgets.js

//platform.twitter.com/widgets.js

//platform.twitter.com/widgets.js

//platform.twitter.com/widgets.js

//platform.twitter.com/widgets.js

That is Pt “double bond” Sn. By no means I’m equating platinum to tin. No sir. Mulliken’s population analysis should not even been brought up anymore, should it?

//platform.twitter.com/widgets.js

//platform.twitter.com/widgets.js

If there is water and ice on Mars then most definitely should be clathrates. (Please excuse the misspelling throughout, though.)

//platform.twitter.com/widgets.js

The rest are the previous announcements which were aimed to generate the momentum for the live tweeting thing.

//platform.twitter.com/widgets.js

//platform.twitter.com/widgets.js
“Not!”

I hope I can make this a thing next year during the 15th edition of RMFQT. I had the honor to be the first speaker and I will upload my presentation soon.

Advertisements

A personal artistic impression of the CompChem Landscape 


In a nutshell, computational chemistry models are about depicting, reproducing and predicting the electronic-based molecular reality. I had this conversation with my students last week and at some point I drew a parallel between them and art in terms of how such reality is approached.

Semi empirical methods
Prehistoric wall paintings depict a coarse aspect of reality without any detail but nevertheless we can draw some conclusions from the images. In the most sophisticated of these images, the cave paintings in Altamira, we can discern a bison, or could it be a bull? but definitely not a giraffe nor a whale, most in the same way Hückel´s method provides an ad hoc picture of π electron density without any regard of the σ portion of the electron density or the conformational possibilities (s-cis and s-trans 1,3-butadiene have the same Hückel description).

More sophisticated semi-empirical Hamiltonians like PM3 or PM6 have better parametrizations and hence yield better results. We are still replacing a lot of information for experimental or adjusted parameters but we still cannot truly adopt it as truthful. Take this pre-medieval painting of one of the first Kings of England, Aelred the Unready. It is, by today standards, a good children´s drawing and not a royal portrait, we now see more detail and can discern many more features yielding a better description of a human figure than those found in Altamira or Egypt.

Ethelred_the_Unready

Æthelred the Unready King of England ca. 1000 BC

 

Hartree-Fock
HF is the simplest of ab initio methods, meaning that no experimental results or adjustable parameters are introduced. Even more so, from the HF equations for a multi-electron system that complies with Pauli’s exclusion principle the exchange operator arises as a new quantum feature of matter with no classical analogue. Still, there are some shortcomings. Correlation energy is disregarded and most results vary according to the basis set employed. Take the impressionist movement, specially in France: In Monet´s Lady with Umbrella we have a more complicated composition, we observe many more features and although we have a better description of color composition some details, like her face, remain obscure. The impressionists are characterized by their broad strokes, the thicker the strokes the harder it is to observe details similar to what happens in HF when we change from a small to a large basis set, respectively.

Claude_Monet_023

Woman with a parasol – Her face or the identity of the flowers at her feet are indistinguishable yet we might be safe to say its springtime.

CI (Configurations Interaction)
Extension of HF to a multi-reference method yields better results. In CI we take the original guess wavefunction -as expressed through a Slater Determinant- and extend it with one or many more wavefunctions; thus a linear combination of Slater Determinants gives rise to a broader description of the ground state because other electronic configurations are involved to include more details like the ionic and covalent pictures (configurations). The more terms we include the more real the results feel. If we take classical figurative paintings we have a similar result; most of these paintings are constituted of many elements and the more realistically each element is captured the more real the whole composition looks even if some are just merely indicated.

Flaming_June,_by_Frederic_Lord_Leighton_(1830-1896)

Flaming June by Lord Leighton – Extreme details on the fabrics and the sea in the background makes us oblivious to the less detailed foot

CCSD(T) full-CI, CASPT2

In Edwards Much’s the scream, we might think we have lost some information again and went back to impressionism but we know this is actually an expressionist painting; we can now not only observe details of the figurative portion of the image but Munch has captured his subject´s fear in the form of distorsions on the subjective reality. In this way, CCSD(T), full-CI and CASPT2 methods provide a description of the ground as well as the excited states which -in experimental reality- are only accessed through a perturbation of the elecron density by electromagnetic radiation. Something resembling radiation has perturbated the subject in The Scream rendering him frightened and wondering how to return to his ground state or if such thing will be even possible.

The_Scream

The Scream by Edward Munch – what sort of perturbation got this guy’s fears out?

 

Density Functional Methods

At least due to its widespread use, DFT has risen as the preferred method. One of the reasons behind its success is the reduced computing time when compared to previous ab initio methods. So DFT is pretty much like photography, in which reality is captured in full but only apparently after selecting a given lens, an exposition, a filter, shutter speed and the occasional Photoshop for correcting issues such as aliasing. In photography, as in DFT, all details concerning the procedure or method for capturing an uncanny reproduction of reality must be stated in every case for reproduction purposes.

Now, in the end it all comes down to Magritte’s Pipe. Ceci n’est pas une pipe -or, ‘this is not a pipe’- reminds us that painting as with modeling we don’t get reality but rather a depiction of it. In this famous painting we look at an image that in our heads resembles that of a pipe but we cannot grab it, fill it with tobacco and smoke it.

MagrittePipe

This isn’t even Magritte’s painting! Let alone a pipe

The image above is a digital file, which translated becomes a scaled reproduction of an image painted by Magritte in which we see the 2D projection of the image of an object that reminds us of a pipe. In fact, the real name of this work is The Treachery of Images, definitely quite an epistemology problem on perception and knowledge but before I get too metaphysical I should finish this post.

Can you find where cubism or surrealism should be placed? with MPn methods, perhaps?

Partial Optimizations with Gaussian09


Sometimes you just need to optimize some fragment or moiety of your molecule for a number of reasons -whether because of its size, your current interest, or to skew the progress of a previous optimization- or maybe you want just some kind of atoms to have their positions optimized. I usually optimize hydrogen atoms when working with crystallographic files but that for some reason I want to preserve the rest of the molecule as refined, in order to keep it under a crystalline field of sorts.
Asking Gaussian to optimize some of the atoms in your molecule requires you to make a list albeit the logic behind it is not quite straightforward to me. This list is invoked by the ReadOptimize keyword in the route section and it includes all atoms by default, you can then further tell G09 which atoms are to be included or excluded from the optimization.
So, for example you want to optimize all atoms EXCEPT hydrogens, then your input should bear the ReadOptimize keyword in the route section and then, at the end of the molecule specification, the following line:

atoms notatoms=H

If you wish to selectively add some atoms to the list while excluding others, here’s an example:

atoms=C H S notatoms=5-8

This list adds, and therefore optimizes, all carbon, hydrogen and sulfur atoms, except atoms 5, 6, 7 and 8, should they be any of the previous elements in the C H S list.
The way I selectively optimize hydrogen atoms is by erasing all atoms from the list -using the noatoms instruction- and then selecting which are to be included in the list -with atoms=H-, but I haven’t tried it with only selecting hydrogen atoms from the start, as in atoms=H

noatoms atoms=H

I probably get very confused because I learned to do this with the now obsolete ReadFreeze keyword; now it sometimes may seem to me like I’m using double negatives or something – please do not optimize all atoms except if they are hydrogen atoms. You can include numbers, ranks or symbols in this list as a final line of your input file.

Common errors (by common I mean I’ve got them):

Lets look at the end of an input I just was working with:

>  AtmSel:  Line=”P  0″
>  Maximum list size exceeded in AddBin.
>  Error termination via Lnk1e in…

AtmSel is the routine which reads the atoms list and I was using a pseudopotential on phosphorous atoms, I placed the atoms list at the end of the file but it should be placed right after the coordinates and the connectivity matrix, should there be one, and thus before any external basis set or pseudopotential or any other specification to be read by Gaussian.

As a sort of test you can use the instruction:

%kjob l103
%chk=myfile.chk
...

at the Link0 section (where your checkpoint is defined). This will kill the job after the link 103 is finished, thus you will only get a list of what parameters were frozen and which were active. Then, if things look ok, you can run the job without the %kjob l103 instruction and get it done.

As usual I hope this helps. Thanks for reading except to those who didn’t read it except for the parts they did read.

%d bloggers like this: