NWChem/doc/user/vscf.tex
Huub Van Dam 97f303e6f8 HvD: In the CVS era the $Id: $ tags in the source code files would
automatically be expanded to include useful information about the
checkin (including the file's revision number). With the switch over
to SVN this was lost because SVN only does this expansion if you 
explicitly ask for it (for every single file). 

I have added a script to the contrib directory that sets the appropriate
property to get SVN to do this expansion. This script will make it easy
to do this every time new source files are added. It is called
svn_expand_Id, the script contains some comments that explain the issue
and how it addresses this.  

This checkin sets this property for a subset of the relevant files
(trying to commit all files at once failed with svn crashing). 
In future the script will only affect those files for which the property
was not set before.
2010-10-29 18:04:21 +00:00

68 lines
2.3 KiB
TeX

%
% $Id$
%
\label{sec:vscf}
The VSCF module can be used to calculate the anharmonic contributions to the
vibrational modes of the molecule of interest. Energies are calculated on a
one-dimensional grid along each normal mode, on a two-dimensional grid along
each pair of normal modes, and optionally on a three-dimensional grid along
each triplet of normal modes. These energies are then used to calculate the
vibrational nuclear wavefunction at an SCF- (VSCF) and MP2-like (cc-VSCF) level
of theory.
VSCF can be used at all levels of theory, SCF and correlated methods, and DFT.
For correlated methods, only the SCF level dipole is evaluated and used to
calculate the IR intensity values.
The VSCF module is started when the task directive
\verb+TASK <theory> vscf+ is defined in the user input file. The input
format has the form:
\begin{verbatim}
VSCF
[coupling <string couplelevel default "pair">]
[ngrid <integer default 16 >]
[iexcite <integer default 1 >]
[vcfct <real default 1.0>]
END
\end{verbatim}
The order of coupling of the harmonic normal modes included in the calculation
is controlled by the specifying:
\begin{verbatim}
coupling <string couplelevel default "pair">
\end{verbatim}
For \verb+coupling=diagonal+ a one-dimensional grid along each normal mode is computed.
For \verb+coupling=pair+ a two-dimensional grid along each pair of normal modes is computed.
For \verb+coupling=triplet+ a three-dimensional grid along each triplet of normal modes is computed.
The number of grid points along each normal mode, or pair of modes can be defined
by specifying:
\begin{verbatim}
ngrid <integer default 16>
\end{verbatim}
This VSCF module by default calculates the ground state (v=0), but can also calculate excited
states (such as v=1). The number of excited states calculated is defined by specifying:
\begin{verbatim}
iexcite <integer default 1>
\end{verbatim}
With \verb+iexcite=1+ the fundamental frequencies are calculated.
With \verb+iexcite=2+ the first overtones are calculated.
With \verb+iexcite=3+ the second overtones are calculated.
In certain cases the pair coupling potentials can become larger than those for a single
normal mode. In this case the pair potentials need to be scaled down. The scaling factor
used can be defined by specifying:
\begin{verbatim}
vcfct <real default 1.0>
\end{verbatim}