| 
 
 NWChem Version 4.0 Notes
NOTE: These release notes are not a
substitute for reading the
User Manual!  They are meant to give
a thumbnail sketch of the new capabilities and bug fixes that
are available in NWChem Version 4.0.  When there is a conflict
between the release notes and the
User Manual, the
User Manual takes precedence.
 
The changes for version 4.0.1 are listed at the bottom of this file.
 
The TITLE directive is no longer backward compatible.
The syntax is
 
title "string title"
 
Several new major capabilities are available in NWChem, including several
new modules.  The new modules are:
 
For each of these capabilities, the user is referred to the
User Manual for more information. ONIOM - ONIOM hybrid mthod of Morokuma and co-workers, COSMO - the continuum solvation 'Conductor-Like Screening' Model
    of A. Klamt and G. Schuurmann to describe dielectric screening effects in
    solvents, NBO - Natural Bond Orbital analysis is now available inside NWChem (not
    just the interface which was in the previous versions), SODFT - Spin-orbit density functional theory, and MD - While this is not a new capability, the MD module has been rewritten
    to allow for new force fields to be incorporated in a modular fashion. 
Additional major additions to the capabilities include:
 
Again, for each of these capabilities, the user is referred to the
User Manual for more information. The Douglas-Kroll(-Hess) approach is implemented. It can be used at
    all levels of theory. Gradients have been implemented allowing geometry
    optimizations and frequency calculations using the Douglas-Kroll
    Hamiltonian, Inclusion of code to evaluate relativistic integrals based on Dyall's
    modified Dirac (DmD) equation.  A one-electron approximation very similar
    to the Douglas-Kroll-Hess approach has been implemented with gradients.
    Two-electron scalar relativistic corrections have been implemented for
    energies only, Vibrational analysis now allows different masses to be used for a
    previously calculated hessian - see VIB description, Input files to standard graphics packages such as RasMol or XMol can
    be generated to animate vibrational modes - see VIB description, and Several new DFT functionals have been added:
  
   Exchange-Correlation PBE0, HCTH402 and Becke1997GGA-1 functionals 
A major bug in version 3.3.1 that has been fixed in version 4.0:
 
A bug was found in the distribution of multidimensional Global Array's.
    This bug created uneven distributions of array elements (i.e. some nodes
    got many elements and some nodes got a few elements), and this had a
    significant impact on the parallel performance through-out the code.
    Note: the final calculated results were correct, the performance was not
    as good as could be expected. 
Listed below are the other major and many minor changes for
each module with significant changes.  These descriptions are somewhat
terse and more information is generally available in the
User Manual.
 
Top level directives:
 
 The TITLE directive has different syntax.  This change is incompatible
    with previous versions of the directive so the user will need to change
    the TITLE directive in any current input files. The CONTINUE directive is no longer supported. Addition of RELATIVISTIC directive to control evaluation of relativistic
    integrals. The default memory has been increased on most workstations.BUG FIX: Task shell is now working under Linux.BUG FIX: The words start and restart in the TITLE directive no longer cause
    problems.BUG FIX: Linux versions no longer give negative time information. 
Initial orbital Guess:
 
BUG FIX: Projection from a smaller basis set to a larger basis set now
    works for open shells. 
Basis Sets:
 
Basis set library has been updated to include all new basis sets at the
EMSL Gaussian Basis Set Order Form page of David Feller Relativistic contractions of the cc-pvnz, n = 2--5, basis sets for the
    Douglas-Kroll and DmD approximations. Dyall double-zeta basis sets for the 4p, 5p and 6p elements addedBUG FIX: Some of the basis sets were corrected:
  
   Basis Set "Ag_Hay-Wadt_VDZ_(n+1)_ECP", correction coefficient in first
      s-contraction Basis Set "Ag_LANL2DZ_ECP", correction coefficient in first s-contraction Complete Basis Set Family "CRENBS_ECP", complete revision of all basis
      sets ECP "Ce_SBKJC_VDZ_ECP", correction switched integer powers for d-f
      component 
Geometry:
 
Inclusion of a finite Gaussian nuclear charge distribution.Inclusion of velocities for QMD and PSPW.More redundant internals are allowed.Allow linear angles.Mixed Cartesians in zmatrix input is more robust.BUG FIX: D5H classification is now
 correct.BUG FIX: C5H operations are recognized. 
DFT:
 
 New functionals: PBE0, HCTH402 and Becke1997GGA-1 Self-interaction corrections. Lebedev + Mura grids are now the default. XC fitting can now use sphericals.BUG FIX: Code correctly executes when the number of iterations is 0. 
SCF:
 
CHANGE: New line search for optimizing the wavefunction.  Does not
    abort when the number of steps is exhausted.  
Driver:
 
BUG FIX: Added code to avoid accidental degeneracies in the Hessian. 
Properties:
 
Note that spherical functions do not work yet with properties. The user is now allowed to choose the point of expansion for dipole,
    quadrupole, and octupole.  See CENTER directive. The center for expansion for dipole, quadrupole and octupole has changed
    from the center of mass to the center of charge.  This behavior can be
    changed by using the CENTER directive. Relativistic corrections to multipole integrals in the Dyall-modified-Dirac
    approach have been added. Performance fix: GIAO convergence has been improved.BUG FIX: Several parallel bugs that caused calculations to fail have
    been fixed.BUG FIX: GIAO code now recognizes when an unsupported wavefunction is
    trying to be used with GIAO. 
MP2:
 
BUG FIX: RIMP2 now has more rigorous error handling. 
Hessian and Frequencies:
 
 New VIB input group to allow for hessian reuse. 
 
Prepare:
 
 PREPARE module can not be used on WIN32 platforms because of problems with
    the Global Array tools. Change in parameter file format. Do not use new files with previous
    versions of NWChem. Z-matrix rules to assign coordinates of missing atoms. Charge separation data from parameter files. Improved performance for large systems. Optional hydrogen build. User defined default HIS protonation state defined in .nwchemrc. Enhanced options for solute alignment, centering, translation, repeat,
    solvent removal, pmf definitions. Automatic counter ion positioning. Option to read old format restart files. Option to write xyz and pov formats. 
MD:
 
 Redesigned MD module. Load balancing takes the cost of I/O into account. New, single trajectory format. Potential of mean force calculations. Improved communication by data packing. Calculation history on restart file. Single point energy calculation prints out solute non-bonded energies,
    analytic and numerical forces. Positional atom constraints. Profiling information. Anisotropic pressure scaling implemented. Browne-Clark leap frog integration implemented (guided MD). Option to center in selected dimension only. Option to exclude pressure scaling in z-dimension. On-the-fly analysis removed.BUG FIX: Modified sub-box redistribution to prevent 'lbbl increase' errors.BUG FIX: The .nwchemrc file now allows for extraneous white spaces.KNOWN BUG: Free energy calculations can not be performed with
    polarizable force fields. 
Analysis:
 
 ANALYSIS module can not be used on WIN32 platforms because of problems with
    the Global Array tools. Essential dynamics analysis. Structural analysis: internal degrees of freedom, RMS deviations, atom
    and charge distributions, atomic distances. Protein structural analysis: Ramachandran and back-bone hydrogen bonds. Electrostatic potential trajectory generation, including periodic
    boundary conditions. Output trajectories in pov format with rendering and color options. Option to read old format trajectory files. 
Integral changes:
 
 Douglas-Kroll(-Hess) integral routines have been added. Relativistic one- and two-electron integral routines for the DmD approach
    have been added. General contraction and multipassing have been implemented for the
    Mc-Murchie-Davidson 2e integrals. Spin-orbit integrals for SODFT enabled. Effective charge potential for COSMO added. 
Dplot:
 
 Gaussian cube file can now be output. Format for InsightII is compatible with version 97.0. Spherical functions can be plotted for MOs and density. 
Tests, manuals and misc script changes:
 
Updated User's Manual to reflect new and modified functionality in NWChem.QA tests have been updated to reflect the new functionality and the
TITLE change. Updated the example directory. INSTALL instructions have been rewritten to be more user friendly. Updates in the contrib directory to reflect more recent machines and
    capabilities, especially the Python capability. 
Installation availability:
 
The current tested platforms and O/S versions are:
 
  IBM SP with P2SC nodes, AIX 4.3.2, and PSSP 3.1 IBM SP with silver nodes (SMP nodes with two 604e processors), AIX
    4.3.2, and PSSP 3.1 IBM SP with Power 3 nodes (SMP nodes with four processors), AIX 4.3.2
    and PSSP 3.1 IBM RS6000 workstation, AIX 4.1 and 4.3 Cray T3E, 2.0.4.61 UNICOSMK SGI R8000/10000/12000 IRIX 6.5 SUN workstations Solaris 5.6 Compaq DEC alpha workstion (600 MHz EV6), Tru64 V4.0 and V5.0,
    DEC C V5.8-009, Digital Fortran V5.2 Linux with Intel x86 cpus. Since there are at least 8 popular
    distributions of the Linux operating system and numerous others in
    existence, including downloading everything and building your own Linux OS,
    it is impossible to test all possible versions of Linux with NWChem.
    NWChem Release 4.0 has been tested on, RedHat 6.0 and 6.2, Debian 2.2.1,
    Suse 6.4, TurboLinux 6.04, Mandrake based on RedHat 6.0, FreeBSD 4.1.4,
    and RedHat 6.0 for the Power PC Macintosh. These all use the EGCS
    compilers at different levels. However, the Portland Group Compiler has
    been tested in a less robust manner. Linux with Alpha cpus has been tested with RedHat 6.1 and 6.2.
    Usage of the Compaq Fortran and C compilers is necessary for compiling. HP 9000/800 workstations with  HPUX B.11.00. f90 must be used for
    compilation. Intel x86 with Windows 98 and NT has been tested with Compaq Visual Fortran
    6.0 and 6.1 with WMPI 1.3 or NT-Mpich. 
 Changes for NWChem 4.0 to 4.0.1
 
 Fix for PeIGS so cross-compiling works on SP systems. PSPW now works on HP Fix for PSPW simulations of charged system Level shift is now excluded from DFT unoccupied orbital energies Fix for reading cartesian hessian from a previous frequency calculation
  when run in parallel QM MD restarts no longer fail under MPI MD restart files are no longer overwritten DNA and RNA fragments were missing some ring closure bonds. In the prepare module, the distance restraint definition is fixed. In the MD module, some printing has been removed when performing
   distance restraints. Order of variables for the LINK card in the PDB file is now correct COSMO gradients are now always calculated numerically. UDFT calculations with COSMO behave better. NWChem now dies gracefully when COSMO is run with either spherical
    basis functions or with ECPs. Documentation has been updated based on the above changes. 
 NWChem |
Capabilities |
Platforms |
Download |
User's Manual |
Programmer's Manual |
Release Notes |
FAQ
 Known Bugs |
Support |
Tutorial |
Contributors |
Benchmarks |
Search |
Mol Sci. Soft. Group |
Citation
 Contact: NWChem Support
 Updated: March 8, 2005
 |