Editing large molecules on a seemingly simple visualizer as GaussView can be a bit daunting. I’m working on a follow up of that project we recently published in JACS but now we require to attach two macrocycles to the organometallic moiety; the only caveat is that this time we don’t have any crystallographic data with which to start. Generating a 3D model of this structure is already hard enough and even when you managed to do it there are many degrees of freedom that in some cases can lead to unrealistic geometries after optimization.
I recently came across a simple way to edit a large complicated molecule by optimizing the fragments separately and then joining them in a new molecule by using the clipboard. This rather simple method, that I for one had never exploited has just saved me a few good hours.
Copy a molecule (CTRL+C) and it will go to the clipboard as a molecular fragment for which you can define a new hot atom and thus bind it to the other fragment as you would with the regular builder. I strongly suggest to use a “New Molecule Group” instead of editing over an existing molecule. Also, if you are using the “paste” button, observe that it has three different options; you may want to use the last one “append to existing molecule” or you will have your fragments in different windows.
And remember, dihedral angles are your best friends.
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.
It’s been a long time since I last posted something and so many things have happened in our research group! I should catch up with them in short but times have just been quite hectic.
Here is a contribution from Igor Marques at the University of Aveiro in Portugal (Group Website); the original text can be found as a comment in the original NBO Visualization post but it is pretty much the same thing you can find in this post. Here is a link to Chemcraft’s website. Thanks for sharing this, Igor!
=> Examples provided by Igor Marques used Chemcraft Version 1.7, build 365 <=
In the Gaussian input, with the NBORead option included under the population keyword, we should include the PLOT option as illustrated below. The gfoldprint keyword will print the basis set to the output file in the old G03 format. Some visualization programs require a certain format of the basis set to be printed to the output file in order to plot orbitals and other surfaces like the electron density; therefore, if you want to play safe, use gfoldprint, gfprint and gfinput in the same line. gfprint will print the basis set as a list but in the new G09 format, whereas gfinput will print the basis set using Gaussian’s own input format. (The used level of theory and number of shared processors are shown as illustrations only; also the Opt keyword is not fundamental to the visualization of the NBO’s)
%chk=filename.chk %nprocshared=8 #P b3lyp/6-311++g** Opt pop=(full,nboread) gfoldprint filename 0 1 molecular coordinates $NBO BNDIDX PLOT $END
this will generate files from *.31 to *.41
For the visualization of NBOs, you’ll need FILE.31 and FILE.37. Open FILE.31 from chemcraft. It will automatically detect FILE.37 (if in the same directory).
Tools > Orbitals > Render molecular orbitals
select the NBOs of interest (whcih are in the same order of the output),
Adjust settings > OK
On the left side of the window, select the NBO of interest and then click on ‘show isosurface’. Adjust the remaining settings. To represent another orbital, click on ‘keep this surface’ and then select another orbital from the rendered set and follow the previous steps.
> It’s possible to open a formated checkpoint file, containing the NBOs, in chemcraft.
%Chk=filename.chk %nprocshared=4 #P b3lyp/6-311++g** Opt pop=(full,nboread,savenbo) gfoldprint filename 0 1 molecular coordinates $NBO BNDIDX $END
the procedure is identical, but it is only necessary to read the *fchk file and then render the desired orbitals.
However, two problems might arise:
a) Orbitals in the checkpoint are reordered, thus requiring some careful inspection of the output.
b) Sometimes, for a larger molecule, the checkpoint might not be properly saved and the Gaussian job (as previously reported – http://goo.gl/DrSgA ) will end with:
Failed in SchOr1 in NBStor.
Error termination via Lnk1e in /data/programs/g09/l607.exe at Wed Mar 6 15:27:33 2013.
As usual, thanks to all for reading/commenting/rating this and other posts in this blog!
Due to extensive popular demand, I hereby make available the necessary files to run Molekel in its old 4.3 version. The program has been compiled to work under Windows 32 bit architecture. Just extract it and place the main folder (as provided here) in any location and run the .exe file located inside. You can generate a direct access to it from your desktop and it even includes a small icon to be used for this purpose.
Also, the manual in pdf format is included; in it you can find the proper citation which must be included in any publication that makes use of Molekel. Just in case you can’t find it, here it is:
MOLEKEL 4.3, P. Flükiger, H.P. Lüthi, S. Portmann, J. Weber, Swiss Center for Scientific
Computing, Manno (Switzerland), 2000-2002.
Stefan Portmann & Hans Peter Lüthi. MOLEKEL: An Interactive Molecular Graphics Tool.
CHIMIA (2000) 54 766-770.
Now some considerations:
- This is an old program. It was generated back in the WindowsXP days, so it wouldn’t be a surprise if it doesn’t work in more recent platforms or under any other 32 bit OS.
- The manual is included. Please read it. This blog is not Molekel’s help desk; I may help but I can’t solve everything, I just don’t have the time for it.
- I didn’t participate/collaborate/helped or got involved in any way in the development of this program, i.e., don’t shoot the messenger. The Molekel homepage is: http://molekel.cscs.ch/wiki/pmwiki.php
- I strongly recommend to make use of the NEW version. A bit more obscure but also great once you figure it out, plus there is available support for it from the actual developers.
- I also strongly recommend to look over the internet for other visualization softwares. I don’t recall having reviewed any in this blog. Perhaps some other time.
- Since this is not my development I will remove it from the server upon the request of the rightful owners A.S.A.P! My guess is they wont mind all that much since its an old version and it was given away for free from their server anyway.
- I can’t think of anything else to put on this list right now but I reserve the right to come back to it and add something more. I just don’t want any trouble.
So, here it is! Right click on the link and download it; Use it to generate nice plots of your orbitals, densities, electrostatic potentials, etc. Consider this a Happy New Year’s gift!
Rate and comment this and all the other posts you find interesting in this blog. Please!
UPDATE: Thanks to Yuekui Wang for the following information.
This copy doen’t work on some WinXP machine with ATI monitor card. The original copy is still available on the cscs web site. Download link is as follows:
It works fine on many macjines, I am sure.