CCL Home Page
Up Directory CCL readme
These files are perl programs which extract coordinates and vibrational modes
from the output file of the appropriate program. Depending on the location
of the perl compiler executable you might need to change the first line.
Also remember to set execution bits using chmod.

The usage is:

G90modes.perl myfile.OUT   (or corresponding for DGAUSS, DMOL, ACES2)

where myfile.OUT is an output file from a Gaussian 90 program. As a result,
the 3N-6 files will be created with the coordinates and modes in the XYZ format
(1st line: no. of atoms, 2nd line: title, here a value of frequency, 
3rd and next:    At-Symbol X, Y, Z, X-mode, Y-mode, Z-mode).
The files have the root myfile and extension nux, i.e.,  myfile.nu1,
myfile.nu2, myfile.nu3... 

The files can be then read into XMol program for viewing the normal modes.
(choose vectors from the menu). To get info about xmol read
files /pub/chemistry/92/06/03 and /pub/chemistry/92/06/12.

Jan Labanowski
jkl@ccl.net

--------- announcement------
Dear Netters,

For some time now I am a "perl" enthusiast (I hope it will pass...)
Perl is a programming language for text processing. It is really a nice tool
to convert file formats, etc. It is a programming language which people
who know C and UNIX can learn in a day. The Perl is described in a man
page on UNIX system (provided that you have it installed). However, man page
is not the best way to learn perl; buy a book instead:
  Larry Wall and Randal L. Schwartz --  "Programin perl", available from
  O'Reilly and Associates, 707-829-0515, e-mail: nuts@ora.com  or
  uunet!ora!nuts (yes, they have e-mail address --- they are UNIX nuts).
The perl is a public domain software. It is made for UNIX, though probably
some ports for VAX exist. It can be taken from anonymous ftp:
archive.cis.ohio-state.edu in /pub/perl/4.0/kits@33  (33 is a patch level,
the larger the better). You get all the kits, uncompress them and do:
  sh perl.kitxx
to extract files. Installation is very easy on popular UNIX machines.

To prove that perl is usefull, I have set up a directory 
   /pub/chemistry/perl
in our archives on www.ccl.net where I placed thre perl programs in:
   /pub/chemistry/perl/vibmodes
Namely: DGAUSSmodes.perl, DMOLmodes.perl, G90modes.perl (you guess what they
do --- extract vibrational modes from corresponding OUTPUT files).
These programs will extract coordinates and vibrational modes in an XYZ format
suitable for viewing with an xmol program (for xmol availability see posts:
92/06/03, 92/06/12) --- a nice molecular viewer for X-windows.
There is also a readme file in /pub/chemistry/perl/vibmodes which tells you
how to use these programs. If you are using OSCPOST e-mail interface,
request files: ./perl/vibmodes/DGAUSSmodes.perl, ./perl/vibmodes/DMOLmodes.perl
./perl/vibmodes/G90modes.perl, ./perl/vibmodes/readme by sending message:
  send the_file from chemistry
to OSCPOST@ccl.net or OSCPOST@OHSTPY.BITNET
If you have ftp, you know what to do.

If anybody has similar perl programs, please contribute them to the archives,
so we do not repeat this boring work of choping and reformating output files.

Jan Labanowski
Ohio Supercomputer Center
jkl@ccl.net
------------------------
From lim@rani.chem.yale.edu  Wed Sep 30 20:53:50 1992
Received: from RANI.CHEM.YALE.EDU  for jkl
        by oscsunb.ccl.net (5.65c+KVa/920330.1102) id AA12438; Thu, 1 Oct 1992 \
04:53:51 GMT
Received: by rani.chem.yale.edu; Thu, 1 Oct 92 04:53:51 GMT
From: Dongchul Lim 
Message-Id: <9210010453.AA01582@rani.chem.yale.edu>
Subject: Re: perl scripts for modes extraction
To: chemistry@ccl.net (Computational Chemistry)
Date: Thu, 1 Oct 92 0:53:50 EDT
In-Reply-To: <199209302156.AA00691@krakow.ccl.net>; from "jkl@ccl.net" at Sep 3\
0, 92 5:56 pm
X-Mailer: ELM [version 2.3 PL11]


jkl@ccl.net writes:
>
> Dear Netters,
>
> There is also a readme file in /pub/chemistry/perl/vibmodes which tells you
> how to use these programs. If you are using OSCPOST e-mail interface,
> request files: ./perl/vibmodes/DGAUSSmodes.perl, ./perl/vibmodes/DMOLmodes.pe\
rl
> ./perl/vibmodes/G90modes.perl, ./perl/vibmodes/readme by sending message:
>
[... a lot of stuff deleted ...]

> #!/usr/local/bin/perl
>
> # Script to extract nomrmal modes from the G90 output file.
> # The OUT file may be the result of single calculations with
> # link line (like example 6.15 in G90 manual) or a result of
> # a separate restart with FREQ keyword from checkpoint file.
> # The first task is to find out if OUT contains two runs
> # (geometry optimization followed by frequences) or just one.
> # Then if 2 runs, skip the portion of file up to and including
> # first @ (it ends the result summary portion for geometry optimization).
> # Then search for "Z-Matrix orientation:" string, skip
> # 4 lines and collect cartesian coordinates in the array
> # until you reach line of "--------------". Count the atoms.
> # Then search for the "Frequencies ---" string, which
> # starts first portion of the table of frequences and
> # normal modes. It looks like:
[... a lot of stuff deleted ...]

> Jan Labanowski
> Ohio Supercomputer Center
> jkl@ccl.net




        The above comments were taken from 'G90modes.perl' which can be
found in 'perl' directory on www.ccl.net. I didn't try this perl
script (we haven't installed perl on our system yet.), but it seems that
the script extracts the cartesian coordinates in zmatrix orientation
rather than standard orientation.
Quoting from Gaussian 90 manual (page 172),

        When evaluation of the second derivatives is complete, link 716
        analyzes the vibrational modes. Its output includes:
        * the dipole moment and the polarizability and hyperpolarizability
          tensors, in the STANDARD ORIENTATION and atomic units.
        * the vibrational frequencies and the associated normal modes.
          these are also in the STANDARD rather than Z-matrix orientation.
                                ^^^^^^^^
(beat me if I'm wrong)
Remedy is simple. Just replace the line containing
        last if /Z-Matrix orientation:/
with
        last if /Standard orientation:/

Then this perl script will be a very useful chemistry tool, for extracting
coordinates, frequencies, normal modes etc. for further application.

-DCL

Department of Chemistry
Yale University


------------------------------------------------------
From jkl@ccl.net  Thu Oct  1 06:21:47 1992
Received: from krakow.ccl.net  for lim@rani.chem.yale.edu
        by oscsunb.ccl.net (5.65c+KVa/920330.1102) id AA18682; Thu, 1 Oct 1992 \
14:21:49 GMT
Return-Path: 
Received: by krakow.ccl.net (5.65c+KVa/901029.12)
        id AA00891; Thu, 1 Oct 1992 10:21:47 -0400
Date: Thu, 1 Oct 1992 10:21:47 -0400
From: jkl@ccl.net
Message-Id: <199210011421.AA00891@krakow.ccl.net>
To: lim@rani.chem.yale.edu
Subject: Re: perl scripts for modes extraction
In-Reply-To: Mail from 'lim@rani.chem.yale.edu'
      dated: Thu, 1 Oct 92 0:53:50 EDT
Cc: chemistry@ccl.net, jkl@ccl.net


Dongchul,

Thank you for your comments on perl modes extraction from G90. I wanted to note
here that "Standard orientation" appears in the G90 output file only if you use
symmetry. If you do not use symmetry, the "Standard orientation" never appears
in the output file and Z-matrix orientation is "Standard orientation".

I will modify the scripts to check if symmetry is used, and if so, the search
will be after "Standard orientation" rather than "Z-matrix orientation".

The script simply reflects my preference not to use symmetry even if it is
there (apparently).

Jan
Ohio Supercomputer Center
jkl@ccl.net

Modified: Mon Nov 16 17:00:00 1992 GMT
Page accessed 10363 times since Sat Apr 17 21:23:00 1999 GMT