| Q0. | How can I get the most recent Updates and Bug Fixes? |
| Q1. | Do earlier NBO 5.0 or NBO 5.G programs work with Gaussian 09? |
| No. However, NBO 5.9 can be "enabled" for either Gaussian 09 or Gaussian 03 and appears to support all non-PBC Gaussian wavefunction options in either version. The legacy NBO 5.0 program must be used for G98 and preceding Gaussian versions. The NBO 5.9 upgrade option allows current NBO 5.0/NBO 5.G license-holders to receive the necessary code for G09/NBO5 capability at reduced cost. | |
| Q2. | For some reason the NRT (or STERIC, CMO, FIXDM, NCS, NJC, etc.) keyword isn't accepted, even though I have the latest version of Gaussian. What's wrong? |
| Gaussian program packages such as Gaussian 09 are distributed with NBO 3.1 (the antiquated 1980s-vintage program that is no longer supported by the NBO development team). You must upgrade to a more recent NBO 5 version (curently, NBO 5.9) in order to have access to these newer NBO-based analysis options (but see Q3). | |
| Q3. | How can I use NBO 5.9 with G09W or other binary PC-Windows programs? |
|
The Windows-based binary program G09W contains
embedded L607 "NBO 3.1"
that cannot be replaced by NBO 5.9. However, the ARCHIVE files
produced by NBO 3.1 are still readable by all post-1980s versions
of the NBO program,
including GENNBO 5.0W. A practical strategy is therefore to
use G09W/NBO3 to generate ARCHIVE (.47) files, which are subsequently used as
input to GENNBO 5.0W for NRT, STERIC, CMO, PLOT,
and other newer options. [Note, however, that the intrinsically
interactive options (NCS, NJC, NEDA, $DEL) cannot be
performed by stand-alone GENNBO 5.0W.] Details of FILE47
format and general GENNBO program usage are given in the NBO 5.0 manual.
The specific steps are as follows: First, prepare the NBO archive (.47) file from your Gaussian run. To do so, just include the "archive" keyword and chosen filename "file=myjob" in the $NBO keylist of your Gaussian input file (with the usual POP=NBOREAD on the route card). You should find the "myjob.47" file produced by this run. Second, modify the myjob.47 file by adding desired keyword options to the $NBO...$END keylist in the second line of this file. Then, with myjob.47 in the same directory, either (a) run GENNBO 5.0W and follow the prompts, or (b) run your OS-specific GENNBO (created with ENABLE.F and system compilers) from the OS prompt, with myjob.47 as standard input (viz., GENNBO < myjob.47 > myjob.out). |
|
| Q4. | With Gaussian I tried to evaluate NPA charges at a correlated MP2 level, but the results seemed to be identical to uncorrelated HF results. What's wrong? |
| You must include "DENSITY=CURRENT" (or "DENSITY=MP2") in the Gaussian route card to analyze higher-level corrections to the HF density. | |
| Q5. | My Gaussian job didn't produce any NBO output but instead gave the message, "NBO is unable to handle linearly dependant basis sets." What does this mean, and what can I do? |
| The Gaussian program checks for numerical instabilities due to near-linear dependence of basis functions (chiefly due to inclusion of diffuse functions) and reduces the dimension of the basis set if necessary. When this occurs, Gaussian prints the above error message and bypasses entry to NBO. One possible solution is to run the same job with G94 or earlier G9X versions. (This usually works, because NBO has its own checks for numerical instabilities.) Another solution is to include IOP(3/32=2) to bypass the Gaussian linear dependence test. However, if linear dependence is truly severe, the only alternative may be to remove "+" or other basis functions until the Gaussian linear dependency fix is not triggered. | |
| Q6. | When I perform $DEL deletions with Gaussian, the program says that the SCF is "not converged." Does this indicate an error? |
| No. The $DEL procedure uses the converged Fock (or Kohn-Sham) operator for single-pass energy evaluation of a deleted density that differs from the "converged" (full) density. Hence, this message can be safely ignored. | |
| Q7. | I compared NPA from a job run on Jaguar and Gaussian and found that many NAO populations differ in the two programs. Who's right? |
| Both are. If internal input coordinates are employed, different program systems often use a different choice of cartesian coordinate system. Under these conditions, the meaning, e.g., of "x" is different (relative to molecular orientation) in the two calculations, and the populations of "p(x)" NAO will differ accordingly. Nevertheless, the two programs should give all the same NAO populations if the coordinate axes are chosen consistently. | |
| Q8. | My MP2 calculation under Gaussian led to many warning messages about "non-physical" occupancies and eventually stopped. What's wrong? |
| MP2 and related methods lead to a perturbative approximation to density that is not physically consistent with an N-particle wavefunction (i.e., not N-representable) order by order. These perturbative inconsistencies are insignificant when the perturbative series is well behaved, but they sometimes become serious (forcing abandonment of NBO analysis) for spin-contaminated open-shell systems and other ill-converged cases. In such cases the NBO error messages warn that the MP2 density is unreliable. The FIXDM keyword can suppress some of these error messages, but cannot "correct" for a physically non-convergent perturbation expansion. | |
| Q9. | My NRT job goes into an infinite hang that requires killing the job. How can I recognize that the job is hung, and what can I do to prevent it? |
| The hang probably results from a case where apparent hypervalency was detected and the program attempted to restart NRT with the NRTFDM option using the full (rather than valence-only) density matrix, but without sufficient memory to accomodate necessary reference structures. (Look for the "apparent hypervalency" message near the bottom of the .LOG file to see if this type of hang is likely.) The solution is to provide LOTS of memory for possible hypervalent cases, and/or to include the NRTFDM keyword so that the full density matrix is used from the outset. In this case, inadequate memory will lead to an early abort, rather than an infinite hang. | |
| Q10. | I tried to go to a higher-level MP2 or CASSCF treatment, but suddenly there were no NBO orbital energies and no table of 2nd-order perturbative energies? What's wrong? |
| NBO evaluates "orbital energies" and 2nd-order stabilization energies only when there is a well-defined 1-electron effective Hamiltonian operator (e.g., Fock or Kohn-Sham operator). Such an operator is unavailable for correlated descriptions, except those of DFT type. | |
| Q11. | I often see a warning message about "population inversion" after NPA analysis. Should I be worried? |
| Probably not. The "inversion" (occupancy ordering inconsistent with energy ordering) probably occurs for Rydberg-type orbitals of very low occupancy or near-degeneracy in occupancy or energy, so that no significant physical effects are indicated. The contrary case of an "inversion" involving high-occupancy, non-degenerate orbitals indicates an excited state. | |
| Q12. | How can I get orbital diagrams of NBOs or other natural localized orbitals? |
| If the ESS has integrated NBO and a read-write file, the checkpointing (or SPARTAN) options may allow you to use standard MO-plotting methods to plot the localized orbitals written (over the MOs) into this file. The NBOView utility program (see order info) provides a more general way to obtain 1-d amplitude profiles, 2-d contour plots, or rendered 3-d photograph-like images for any desired localized orbitals, using the input files created by the PLOT keyword. NBOView orbital imagery is illustrated in the homepage logo and elsewhere throughout this website. | |
| Q13. | I can't compile gennbo.f with the g77 compiler on my linux system. What's wrong? |
|
Use the compiler directive g77 -Wno-globals -fno-globals gennbo.f -o gennbo to bypass checks for strict consistency between defined vs. called subroutine argument lists. |
|
| Q14. | My nbo_59.src file contains many "^@" control characters that were not on the original CD. What happened? |
| The .src textfile may have been transferred between different operating systems using a Windows-type ftp utility with "T/B" (rather than text "T") as the default file-type setting. | |
| Q15. | I am working with lanthanides and/or actinides. Does the NPA partitioning scheme recognize the 5d (for lanthanides) and 6d (for actinides) subshell as a "valence" subshell, and will this affect my NPA charges? |
| All recent NBO 5 versions (including NBO 5.9) recognize these subshells as valence orbitals, consistent with the observed partial occupancy of these subshells in certain ground-state atoms of the "f-block" lanthanides and actinides [J. Comput. Chem. 28: 198-203, 2007]. (The NPA partitioning for lanthanides/actinides is now consistent with long-established treatment of main and transition blocks and known ground-state atomic configurations. Older NBO versions conformed to the Madelung Rule, which forces some lanthanides and actinides to be analyzed in terms of excited-state atomic configurations.) This change can affect the NPA charges [J. Chem. Phys. 121, 2563-2570, 2004]. For the code-change to make older NBO versions consistent with current NBO 5.9 in this respect, please see the update. | |
| Q16. | How can I use NBO 5.9 with NWChem? |
| The interface code to NWChem 3.2 no longer works with the current NWChem versions. However, the nwnbo.F code produced from NBO 5.9 can be patched with nwnbo.patch (patch -p0 <nwnbo.patch) and compiled with NWChem 5.1.1 (June 2008 release). For users with NWChem 5.1, the NBO_Fix.tar tar-file contains the nwnbo.patch and all needed files to enable NBO to be compiled into NWChem 5.1. The tar-file includes a README document outlining how to perform the patch operation. [Thanks to Dr. Bert de Jong (PNNL, bert.dejong@pnl.gov) for these updates and utilities.] | |
| Q17. | The NBOView for MacOSX installation instructions failed to work on my newer model Macintosh. What's wrong? |
|
The libraries used by NBOView are not consistently available on newer
Macintosh models. You can use Xcode 3.1.2 (downloaded from Apple
Developers Web Site) and the gnu compiler gcc4.4 for Intel macs (downloaded
from http://hpc.sourceforge.net/) to compile and link the NBOView program
with the following commands (using the Terminal application):
(1) Compile the draw.c program with the gnu compiler using the following gcc command: gcc -c -D_XWINDOWS -D_MACOSX draw.c (2) Now compile the nboview.f fortran code and link it to the draw.o object file and X-windows library with the gnu compiler using the following gfortran command: gfortran nboview.f draw.o -L/usr/X11R6/lib -lX11 -o nboview (3) Finally, remove the draw.o object file: rm draw.o You should now have an executable nboview program that can be invoked by typing ./nboview from the Terminal application. |
|
| Q18. | My NICS-type NMR calculation using G03/NBO5G with the NCS keyword failed with the error message "Subroutine NAOANL could not find a s-type valence orbital on atom gh ..." What's wrong? |
| The problem stems from a Gaussian change in handling Bq ghost atoms between revisions B and C of G03. Early NBO 5.G distributions give such error when linked to G03 Rev. C (or later). The problem was corrected for later G03/NBO5.G combinations, but apparently recurs in current G09/NBO5.9. Until further notice, NICS users should revert to older G0X/NBO5 code combinations that reproduce trusted literature values. | |
| Q19. | My attempt to compile GENNBO for Mac OSX platorm failed with g77. With gfortran, compilation apparently succeeded (no error message), but the resulting gennbo.exe program terminated abruptly. What's wrong? |
|
For obscure reasons, 64-bit Mac OSX systems must be selected with
the "32-bit integer" option in the enable.f menu procedure for use with
available open-source compilers. The resulting gennbo.f code can then
be compiled with gfortran (not g77, which is apparently at the end of
its development cycle), using, e.g., gfortran gennbo.f -o gennbo_32 (Ignore the many warnings about HOLLERITH to INTEGER conversion.) This produces a gennbo_32.exe that should run successfully as a 32-bit application. |
|
| Q20. | My DFT deletions (POP=DEL) jobs with G09 no longer agree with corresponding G03 jobs. What's wrong? |
|
An integration option was modified in G09 that affects
all $DEL evaluations with DFT methods. All G09 users should re-run
such jobs with the IOp work-around illustrated below:
#B3LYP/6-31+G* Pop=NBODel SCF=NoVarAcc IOp(5/48=100000) The error affects current (12/15/2011) revision C and former revisions A, B of G09, but is to be remedied in forthcoming revisions. |
|
| Q21. | My attempt to compile gennbo.f on Ubuntu 11.10 (32 bit x86) with gfortran gave an executable that did not work with ORCA. I then installed and tried the g77 compiler (as in Q13), but this failed with "cannot find" error messages for crtl.o, crti.o ("No such file or directory") and -lgcc_s. What's wrong? |
|
The g77 compiler is not finding needed shared libraries. Two steps are needed to fix this: Step 1 (in terminal): LIBRARY_PATH=/usr/lib/i386-linux-gnu:/usr/lib/gcc/i486-linux-gnu:$LIBRARY_PATH export LIBRARY_PATH Step 2 (in terminal, in i486-linux-gnu directory): The broken symlink to the missing libgcc_s.so library (/usr/lib/gcc/i486-linux-gnu/libgcc_s.so) can be repaired with the command: sudo ln -s /lib/i386-linux-gnu/libgcc_s.so.1 libgcc_s.so Then the g77 compile command of Q13 should work correctly. |