Author Archives: joaquinbarroso
As far as population analysis methods goes, the Quantum Theory of Atoms in Molecules (QTAIM) a.k.a Atoms in Molecules (AIM) has become a popular option for defining atomic properties in molecular systems, however, its calculation is a bit tricky and maybe not as straightforward as Mulliken’s or NBO.
Personally I find AIM a philosophical question since, after the introduction of the molecule concept by Stanislao Cannizzaro in 1860 (although previously developed by Amadeo Avogadro who was dead at the time of the Karlsruhe congress), the questions of whether or not an atom retains its identity when bound to others? where does an atom end and the next begins? What are the connections between atoms in a molecule? are truly interesting and far deeper than we usually consider because it takes a big mental leap to think about how matter is organized to give rise to substances. Particularly I’m very interested with the concept of a Molecular Graph which in turn is concerned with the way we “draw lines” to form conceptual molecules. Perhaps in a different post we can go into the detail of the method, which is based in the Laplacian operator of the electron density, but today, I just want to collect the basic steps in getting the most basic AIM answers for any given molecule. Recently, my good friend Pezhman Zarabadi-Poor and I have used rather extensively the following procedure. We hope to have a couple of manuscripts published later on. Therefore, I’ve asked Pezhman to write a sort of guest post on how to run AIMALL, which is our selected program for the integration algorithm.
The first thing we need is a WFN or WFX file, which contains the wavefunction in a Fortran unformatted file on which the Laplacian integration is to be performed. This is achieved in Gaussian09 by incluiding the keyword output=wfn or output=wfx in the route section and adding a name for this file at the bottom line of the input file, e.g.
(NOTE: WFX is an eXtended version of WFN; particularly necessary when using pseudopotentials or ECP’s)
Analyzing this file requires the use of a third party software such as AIMALL suite of programs, of which the standard version is free of charge upon registration to their website.
OpenAIMStudio (the accompanying graphical interface) and select the AIMQB program from the run menu as shown in figure 1.
Select your WFN/WFX file on which the calculation is to be run. (Figure 2)
You can control several options for the integration of the Laplacian of the electron density as well as other features. If your molecules are simple enough, you may go through with a successful and meaningful calculation using the default settings. After the calculation is finished, several result files are obtained. We’ll work in this tutorial only with *.mpgviz (which contains information about the molecular graph, MG) and *.sum (which contains all of needed numerical data).
Visualization of the MG yields different kinds of critical points, such as: 1) Nuclear Attractor Critical Points (NACP); 2) Bond Critical Points (BCP); 3) Ring CP’s (RCP); and 4) Cage CP’s (CCP).
Of the above, BCP are the ones that indicate the presence of a chemical bond between two atoms, although this conclusion is not without controversy as pointed out by Foroutan-Njead in his paper: C. Foroutan-Nejad, S. Shahbazian and R. Marek, Chemistry – A European Journal, 2014, 20, 10140-10152. However, at a first approximation, BCP’s can help us to explore chemical interactions.
Now, let’s go back to visualizing those MGs (in our examples we’ve used methane and ethylene and acetylene). We open the corresponding *.mpgviz file in AIMStudio and export the image from the file menu and using the save as picture option (figure 3).
The labeled atoms are NACP’s while the green dots correspond to BCP’s. Multiplicity of a bond cannot be discerned within the MG; in order to find out whether a bond is a single, double or triple bond we have to look into the *.sum file, in which we’ll take a look at the bond orders between pairs of atoms in the section labeled “Diatomic Electron Pair Contributions and Delocalization Data” (Figure 4).
Delocalization indexes, DI’s, show the approximate number of electrons shared between two atoms. From the above examples we get the following DI(C,C) values: 1.93 for C2H4 and 2.87 for C2H2; on the other hand, DI(C,H) values are 0.98 for CH4, 0.97 in C2H4 and 0.96 in C2H2. These are our usual bond orders.
This is the first part of a crash tutorial on AIM, in my opinion this is the very basics anyone needs to get started with this interesting and widespread method. Thanks to all who asked about QTAIM, now you have your long answer.
Thanks a lot to my good friend Dr Pezhman Zarabadi-Poor for providing this contribution to the blog, we hope you all find it helpful. Please share and comment.
I had a blast last week at WATOC2014 in Santiago de Chile! It was a wonderful opportunity to find old friends, meet new ones and listen to some exciting research done around the world, as well as some of the classics such as Pekka Pyykkö, who was awarded the Schrödinger medal. I decided to share my talk on SlideShare.com but also here because I found at WATOC that many many people seem to like this little space of mine! I was shocked, flattered but mostly happy to know that this little blog of mine is well regarded.
So, without further ado, here is my presentation at WATOC2014, please read the captions on each image for context. Feel free to make any comments, sharing or liking. Thanks for clicking!
If you made it this far, let me tell you that this is also available at Slideshare.com :)
Thanks for reading, commenting and sharing!
I inexcusably forgot to write about my visit to WATOC2014 but as they say it’s better late than never. I’ve been here in Santiago de Chile for a day and a half now attending mostly to session 5 “applications to compelling problems”. WATOC (World Association of Theoretical and Computational Chemists) is probably the largest conference in our speciality and some of the big names are here which makes me very nervous to present my talk!
My talk will be very similar to the one presented in the last ACS meeting in San Francisco, in silico design of monomolecular drug delivery agents based on calix[n]arene macrocycles.
I wont try to report from WATOC, I’m not any good at that but will for sure be available if anyone wants to contact me.
This is my first time giving a talk at the ACS meeting so I’m a tad nervous but also excited about the possibility to listen to so many great researchers and hopefully interact with some of them, although such a massive event is hardly the best place to build a network.
My talk will deal with our recently published work on calixarene based drug carriers and the new results we are currently getting, and will take place on Tuesday 12th at 16:15 at the session regarding the use of small macrocycles in drug discovery. Feedback is very important to us so we look forward to some heavy criticism of our results.
If you are there, download the free app which will show you around the different presentations and stands.
Well, I only contributed with the theoretical section by doing electronic structure calculations, so it isn’t really a paper we can ascribe to this particular lab, however it is really nice to see my name in JACS along such a prominent researcher as Prof. Chad Mirkin from Northwestern University, in a work closely related to my area of research interest as macrocyclic recognition agents.
In this manuscript, a calixarene is allosterically opened and closed reversibly by coordinating different kinds of ligands to a platinum center linked to the macrocycle. (This approach has been referred to as the weak link approach.) I recently visited Northwestern and had a great time with José Mendez-Arroyo, the first author, who showed me around and opened the possibility for further work between our research groups.
Closed, semi-open and fully open conformations; selectivity is modulated through cavity size. (Ligands: Green = Chloride; Blue = Cyanide)
Here at UNAM we calculated the interaction energies for the two guests that were successfully inserted into the cavity: N-methyl-pyridinium (Eint = 57.4 kcal/mol) and Pyridine-N-oxide (Eint = +200.0 kcal/mol). Below you can see the electrostatic potential mapped onto the electron density isosurface for one of the adducts. Relative orientation of the hosts within the cavity follows the expected (anti-) alignment of mutual dipole moments. At this level of theory, we could easily be inclined to assert that the most stable interaction is indeed the one from the semi-open compound and that this in turn is due to the fact that host and guest are packed closer together but there is also an orbital issue: Pyridine Oxide is a better electron acceptor than N-Me-pyridinium and when we take a closer look to the (Natural Bonding) orbitals interacting it becomes evident that a closer location does not necessarily yields a stronger interaction when the electron accepting power of the ligand is weaker (which is, in my opinion, both logic and at the same time a bit counterintuitive, yet fascinating, nonetheless).
All calculations were performed at the B97D/LANL2DZ level of theory with the use of Gaussian09 and NBO3.1 as provided within the former. Computing time at UNAM’s supercomputer known as ‘Miztli‘ is fully acknowledged.
The full citation follows:
A Multi-State, Allosterically-Regulated Molecular Receptor With Switchable Selectivity
Jose Mendez-Arroyo †, Joaquín Barroso-Flores §,Alejo M. Lifschitz †, Amy A. Sarjeant †, Charlotte L. Stern †, and Chad A. Mirkin *†
Thanks to José Mendez-Arroyo for contacting me and giving me the opportunity to collaborate with his research; I’m sure this is the first of many joint projects that will mutually benefit our groups.
Just as last year, the “Dolphin Summer Internship Program” (Programa Delfín) has started and this time it coincided with #RealTimeChem week. Four students from various cities (and accents) around Mexico have come to our lab in Toluca in order to spend about 7 weeks of research in the field of molecular modeling and within our research of molecular recognition in biochemistry. Karen, Cynthia, Jesús and Marco have started their training today as they arrived to CCIQS so we went over the (very) basics of quantum chemistry, the (very) basics of Linux and the basics of Gaussian09. (I should really think about developing some web tutorials or something because this impromptu training is very exhausting!)
Their academic backgrounds are mostly centered around pharmaceutics and biochemistry although their ages range from the second to the fourth year of college education. Computational chemistry is pretty unknown to all of them; I’ll do my best to change that, while at the same time I make them aware of its power as a research tool and as a research field in itself.
Here is to a very productive summer! I hope we manage to get enough data for a paper and, more importantly, that they all get a good experience out of their time here, make new friends and learn something new that enriches their skills in this increasingly competitive world.
A couple of weeks ago I was invited to give a talk to a small university in southern Mexico called ‘Universidad de la Cañada‘ in the state of Oaxaca, one of the most underprivileged states in our nation. This institution is a rather small one but the work they are doing over there with as little resources as they have is truly remarkable . UNCA offers degrees in pharmacy, pharmacology, food sciences, clinical chemistry and other topics that aim to supply the needed human resources for the various industries that are settled in the region. There is a true feeling of togetherness at UNCA since they have little pieces of equipment yet they are all fully shared among researchers regardless of who received the finance to acquire them. Last year, two of their students came for a two months stay, after which, Alberto and Eduardo got their names on a publication of our research group. It was nice to see them again and even nicer to learn they are about to finish their studies and that they will come back again to our lab in late July.
Every year at UNCA there is a Pharmacology Day on which the students show the results to their research projects during a poster session and listen to lectures by guest speakers from various universities around Mexico. Most of their projects were aimed to the isolation of natural products from local resources and their usage in several kinds of consumer products. UNCA is in a very small town, village I might say, surrounded by mountains and vegetation; the view was spectacular as you may see from the pictures below. Thank you very much to my good friend Dr. Carmen Hernández-Galindo for inviting me to participate and share our work with their students, I hope we may go back again and keep a fruitful exchange between our groups.
During this talk, I took the opportunity to talk about the aforementioned paper in the context of molecular recognition and their in silico design but I think I should have talked more about the computational strategies that are most employed in the pharmaceutical industry. Never mind. I hope I get the opportunity to right this wrong. Still it was nice to give Alberto and Eduardo the opportunity to brag a little about being published authors.
Kudos to Rola Aburto, Dr. Margarita Bernabé, Dr. Rocío Rosas, and all the academic staff at UNCA for their invaluable dedication to teaching science against all odds, I can testify, through the hard work of their students, hat their effort is paying off.
So the World Cup is once again on top of us. I’m not a Football (Soccer) enthusiast but I’ve got to admit that the expectation of such a large and widely covered event is pretty contagious. This year, however, I’m very excited about the inaugural kick-off ceremony because a paraplegic teen will be the one to set the ball in motion, thanks to the use of an exoskeleton developed by the illustrious Brazilian researcher, Dr. Miguel Nicolelis, this patient will not only walk again but also perform a feat of equilibrium: kicking a football. More impressive than the exoskeleton itself is the brain-computer-machine interface since the patient will control the entire process by himself. Miguel Nicolelis is widely known and highly regarded in the scientific community; I’m not sure if he is that famous outside academia, but if he isn’t, he should be. The natural question about Dr. Nicolelis is what is he? Is he a robotics engineer? a neurologist? a programmer? a physician? The answer could be no other than ‘all of the above‘.
And even more impressive than all that, if that’s even possible, is the fact that this huge achievement of technology is presented at one of the most viewed sporting events on the planet. Brazilian organizers could have selected many things to kick-off this event: From Adriana Lima to Pelé; from a Samba line to aboriginal Amazonian people, but instead they chose to go with a scientific and technological breakthrough achieved by one of their own. I wonder if this is a way to tell the world they are interested in investing in science and technology as a way to pave the way of their economical and social development. Brazil is currently regarded as a fast growing nation economically although the social disparity seems to be still quite large. The message I’m getting, at least in principle, is that Brazil is a modern nation with high regard for scientific development on which they will rely their future.
Kudos to the Brazilian organizers who thought of placing this large scientific breakthrough in a sporting event, proving that this world should become boundless and the way to do it is through science.
On Friday May 30th, my good friend Dr. Josefina Aldeco, my wife and I, visited a children’s home in Querétaro (central Mexico) and brought them a few cool chemistry experiments for a short show. This event was promoted by a non-profit organization called “Anímate a estudiar” (Dare to study), namely by Mrs. Paulina Milanés who is always looking for ways to encourage kids from poor backgrounds to pursue their goals through study; among other things, they provide backpacks with school supplies to orphan kids like the girls we visited.
As a way to inspire them, we handed each girl a balloon drawn in the shape of a brain and asked them to inflate them daily by reading; by doing their homework; by asking questions all the time; by working hard in pursuit of a brighter future for which their brains are the most powerful muscles.
Many reactions took place that Friday; not only inside the flasks and beakers before our little audience but also in their faces and their engagement with us. Little by little these girls got out of their shells and became more excited, up to the point of performing their own chemical reaction themselves by polymerizing some glue with borax in hot water. This was for sure the first time they got in contact with chemistry but the true goal was to set up a spark in their minds that one day may turn into a life opportunity. We are aware that one small chemistry show can’t really have that effect, but if many more scientists reach out to these kids there is a bigger chance of creating a ripple effect that convince disenfranchised children that studying is the way to take the wheel of their own future.
Science is about development; its about spreading knowledge and the love for knowledge. Although we most times sit high on our ivory towers it is paramount to remember that there is also a social component to the scientific activity. Kids are eager to learn, but most school systems do their very best to limit their curiosity and ambition. We hope these girls find in studying a way to a better, happier and safer future. Mexico has a large economic disparity; climbing the social ladder is very hard and even more so for women which makes these girls a very vulnerable social group in the next generation.
It only takes one day. One day and some potassium iodide; some mentos on a diet-coke (sorry, Gina, for the squirt!); some cobalt chloride on paper; some balloons some glue and some borax in hot water. But above all it takes a big commitment.
I hope you readers, computational and experimental chemists alike, take some time out of your busy schedules and share your passion for science with kids, specially those with the lowest opportunities of getting in touch with real scientists. You can also contribute to this noble effort by making a small pay-pal donation to www.animateaestudiar.org or to any other similar organization in your local community.
It only takes one day.
P.S. Thanks to Josefina from Universidad Autónoma de Querétaro for providing material and reagents. Please go and check out her blogs (here and here) and encourage her to write more often! (Did I mention she published in Science a few years ago?)
A couple of weeks ago I posted a solution for a common error regarding .fchk files that will display the error below when opened with GaussView5.0. As I expected, this error has to do with the use of diffuse functions in the basis set and is related to a change of format between Gaussian versions.
CConnectionGFCHK::Parse_GFCHK() Missing or bad data: Alpha Orbital Energies Line Number 1234
Although the method described in the previous post works just fine, the following update is a better approach. Due to a change of spelling between G03 and G09 (which has been corrected for G09 but not available for GV versions prior to 5.0.9) one must change “independent” for “independant“
To make the change directly from the terminal the following command is needed:
sed -i 's/independent/independant/g' file.fchk
Alternatively you can redirect the output to a new file
sed -e 's/independent/independant/g' file.fchk > newfile.fchk
if you want to keep the old version and work with a new one.
Of course this edition can be performed manually with any text editor available (for example if you work in Windows) but solutions from the terminal always seem easier and a lot more fun to me.
Thanks to Dr. Fernando Cortés for sharing his insight into this issue.