![]() Gaussview will automatically format the chk file (this can cause some troubles when crossing architectures or platforms, so be careful to generate the proper formatted chk file!) Once load select the orbital you need to visualize and go to the ‘Visualize’ tab and click update the orbital displayed is the Natural Bond Orbital. Click on the ‘New’ tab and load your checkpoint file. This will open a new window with four tabs at the bottom of the molecule image. ![]() Load the output file with gaussview normally and then click on the Molecular Orbital icon (or go to Edit -> MO’s). This will save the NBO coefficients into the checkpoint file. #P opt HF/6-31G(d) freq=noraman pop=(ful,nboread,savenbo) geom=connectivity When performing the NBO analysis on Gaussian (by the way, Gaussian 09 cannot be visualized with gaussview 3.x) use the savenbo option in the rout section, for example: This is the gaussview version I currently work with. Some people complain about the look and feel of the latest molekel version so they stick to this old bugged one, so that is why I’m posting this method. ![]() Sometimes the window closes because some sort of resolution problem, specially when taking snapshots (interestingly enough this happened to me when the background color was set to white) decrease the resolution of your monitor before taking the snapshot to prevent this problem. Molekel 4.3, though, is filled with bugs that will make it close unexpectedly, specially when running under Windows Vista. This should do it! As an interesting exercise try computing the same orbitals (lets say HOMO and HOMO-1) with and without loading filename.47 in order to observe the difference between the shape of the MO’s and NBO’s. Now, go to Compute -> Orbital and now select the Natural Bond Orbital you are interested in. Then go to Load -> nbo orb and load filename.47. That aside, load your output file ( filename.out or filename.log) on Molekel 4.3 as usual. If you are using Gaussian 03 or Gaussian 09 and STILL want to use Molekel 4.3 then you probably already know you have to change 03 for 98 on the header of the output file: The one you want to pay attention to is filename.47 The PLOT option will generate a series of files with numerical extensions. You may include other keywords such as BNDIDX which generates a Wiberg bond index ( order) matrix or BOAO which generates the same matrix but in the Natural Atomic Orbitals basis. NOTE: I have my own issues and concerns about the use of DFT along with NBO but maybe that will be treated in another post.Īt the end of the input file, after the infamous blank line type the following: This calculation requests a geometry optimization followed by a Natural Bond Orbital population analysis (with keywords to be read) using the Restricted Hartree-Fock method with the split valence basis set 6-31G(d,p) #P opt rhf/6-31G(d,p) pop=(full,nboread) gfoldprint geom=connectivity Make sure that the route section in your input file includes the following options: pop=nboread gfoldprint (in case you are using G03 or even G09) or gfprint (in case you are still running G98), e.g.: Most of the answers are concerned with what visualization programs to use but they seldom provide step by step instructions, furthermore most manuals are a bit cryptic about the procedure to plot this orbitals. Visualizing them, however is never a straightforward process, and quite often we see more questions on the CCL than answers trying to address the matter. So far I have only worked with Gaussian 98, Gaussian 03 and only recently with Gaussian 09 in calculating NBO’s although it is possible also to perform them with GAMESS and the standalone NBO5.0 program created by Frank Weinhold. Another advantage of the NBO analysis is that it provides a localized depiction of the electron density over a molecule, making it more related to chemists intuition. is highly sensitive to basis set effects. The Natural Bond Orbitals analysis is a powerful tool in population analysis calculations which is more robust than the traditional Mulliken approach, if for no better reason because its almost insensitive to the change of basis set while Mulliken’s P.A. This is the second post on a series which will try to address common technical questions in computational chemistry that recursively appear on the CCL.
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |