| ASFIT |
ASymmetric
top FITting program
|
This
program embodies most features deemed necessary for this type of
calculation:
- Watson's reduced asymmetric top Hamiltonian
up to decadic terms, with Wang factorization
- Reductions A and S,
representations Ir and IIIr (note
the idiosyncrasy that the latter is not IIIl so that
signs of the off-diagonal quartic, sextic.... constants will be
reversed)
- J up to 250, N.lines up to 10000,
the latter can easily be increased by changing NLINES in PARAMETER statements
- Lines can have differing measurement errors
for unequal weighting and fitting of blended lines is supported
- Lines can carry annotations, both at the end
of each transition, as well as between transition declarations, in the
form of full comments, or just spacing lines
- Automatic generation of entries for the
various Q-type bands and the high-J, R-type
bands (see J.Mol.Spectrosc. 177,240(1996) and 178,125(1996))
- Publication type formatting of resulting
constants and cleaned up output of the frequency table for electronic
deposition purposes
- Output of the correlation matrix, stick
diagram of correlation coefficient distribution, contributions to
frequencies, and summary of worst lines in the fit
- Output of the hii
coefficients for identifying the most influential transitions (those
with hii>0.5 are suspicious, those
with hii>0.9 will definitely determine
the fit), and of DFBETAS diagnostics for sensitivity of fitted
constants to transitions in the data set
- Output of a constants file for use by ASROT and of .lin and .par files for use by SPFIT
- On-line modification of data: line
exclusion, insertion, deletion, as well as on-line decisions concerning
the parameters of fit
- checks for common daset errors such as
duplicated lines, incorrect blend specification and split blends
| Some background
references: |
|
| |
|
| The reduced rotational
Hamiltonian: |
J.K.G.Watson in "Vibration
Spectra and Structure" (J.Durig, Ed.), Vol.6, p.1, Elsevier, Amsterdam,
1977. |
| An introduction to the
least-squares fitting of spectroscopic data (seminal tutorial): |
D.L.Albritton, A.L.Schmeltekopf,
R.N.Zare, in "Molecuar Spectroscopy: Modern Research" (K.Narahari Rao,
Ed.), Vol.II, p.1, Academic Press, New York, 1976 |
| Introduction to the hitherto rarely
used, but very useful additional statistical tests hii,
DFBETAS
etc.: |
J.Demaison, J.Cosleou,
R.Bocquet, A.G.Lesarri, J.Mol.Spectrosc. 167,
400-418
(1994) |
| The trick for
computationally cheap derivatives by application of the Hellman-Feynman
Theorem: |
W.F.Rowe, E.B.Wilson, J.Mol.Spectrosc.
56, 163-165 (1975) |
| Calculation of the higher
order matrix elements and plot of distribution of correlation
coefficients: |
Z.Kisiel, L.Pszczolkowski, J.Mol.Spectrosc.
158, 318-327 (1993) |
| |
|
 |
| |
| |
|
| ASFIT.FOR |
The
listing - it is recommended that extension .ASF be used for the
data files, .RES for the results files, and .CON for the constants files for use by
ASROT |
| ASFIT.EXE |
Executable for those
who have problems compiling for Windows - should work for W95/Pentium I
upwards. |
| |
|
| |
Examples |
| CB.ASF |
Specimen
data set for chlorobenzene.
NOTE: In some Windows configurations
your browser might want to open the .ASF files directly with MediaPlayer or
another multimedia player, since this is one of the default extensions
for movie formats. The result is of course an error message - in such
case use the SAVE TARGET AS or the SAVE LINK AS option of the browser
(obtained with right mouse click on the link), and only then inspect
the locally saved file with a suitable editor.
|
| CB1.ASF |
Second data set for
chlorobenzene - this one goes up to J=249 for testing PENTIUM
III's etc. Note the use of several types of clarifying annotations,
which are carried over into the output files. |
| BENCH.TXT |
List of
benchmarks for this program + chlorobenzene data set recorded on
various computers: this is really a trivia item reflecting the rapid
development of computers |
| HSSH.ASF |
Some of the first
data which showed that S-reduction is really necessary (read
this note if your browser wants to open this as a multimedia file) |
| T2O.ASF |
Tests various
high-order c.d. terms, in this case reduction-A, representation
Ir (read this note if your browser wants to open
this as a multimedia file) |
| T2O.RES |
Results file for
T2O.ASF |
| T2O.CON |
Constants file
produced by ASFIT from data T2O.ASF for use by ASROT |
| |
Associated programs |
| LINASF.FOR |
Converts a .LIN
file of SPFIT into a block of frequencies in the .ASF standard. For a complete .ASF file it is still necessary to manually
add the header with values of constants. |
| LINASF.FOR |
Win32 executable |
| SORTAS.FOR |
Sorts transitions in an .ASF file according to many different
criteria. Do not use on data files containing blends, since the
constituent lines may not stay together. |
| SORTAS.EXE |
Win32 executable |
|
|
 |
Back to the
table of programs
| ASROT |
ASymmetric
ROTor predictive program
|
Calculation of frequencies and/or energy levels for an asymmetric rotor
by using Watson's reduced Hamiltonian in representation Ir,
and
reductions A and S.
- The current limits are J<300,
and 91204 lines in the output file. The line limit is partly due to the
now historic use of temporary storage of calculated lines prior to
sorting in frequency
- Calculated lines carry intensity, line
strength, lower level energy, and transition type
- All data necessary for the calculation can
be input interactively. The input list is quite long so that more often
a pipeline of the type ASROT<DATA.IN is used, where the file DATA.IN contains the necessary input
- ASROT can use constants files produced by ASFIT
- Most usefully in our practice - results
files from ASROT can be displayed graphically with ASCP
 |
| |
| ASROT.FOR |
The
listing - it is recommended that extension .ASR be used for the
results files |
| ASROT.EXE |
Executable for those who have
problems compiling for Windows - should work for W95/Pentium I upwards. |
| TRIAL.IN |
Trial input file, which can be
reedited and used for quick predictions in pipeline mode run with the
command of the type ASROT<TRIAL.IN |
| TRIAL.ASR |
The results file obtained from TRIAL.IN - this can
be viewed with ASCP |
|
|
Associated programs |
ASRBAND.FOR
ASRBAND.EXE
|
Converter of an abbreviated .ASR file into a .DAT file for use by
gle, such that
it results in an overbar with markers annotating positions of lines in
a band. Caters for reversing bands, up/down markers, longer selected
markers, and choice of quantum number determining marker length.
Documentation is at the top of the program
listing. A similar program, CATBAND, is available for use with .CAT output
from SPCAT.
|
| ASRGLE.FOR |
Converter of an .ASR file into a
data file for plotting a stick diagram of the calculated spectrum with gle. This method has now been superseded by direct gle output from ASCP, but
there are still some options, such as line blending unique to this
program. |
|
|
 |
Back
to the table of programs
| ASCP |
Graphical previewer
for spectral predictions generated by ASROT and by Pickett's SPCAT
|
Please note that the currently supported
version is ASCP_L.
Even
though that program is part of the AABS package, it is also
recommended for standalone operation. The notes below are intended to
serve as additional documentation, but the executables now have legacy
status.
We
find this program invaluable, both in MMW and FTMW work. It started
life as ASPEC, which printed the stick spectrum on a band of paper. The
resulting colour graphics version was called ASC, and the name ASCP came
into use for the version, which could display processed output from
Pickett's SPCAT.
The
current version is the most Windows compliant version so far, and has
been tested from W95 up to XP. The default program size is 800x560
points, but it is possible to change this with a suitable configuration
file. The principal features are:
- Display of up to 100000 lines from up to 20
different data sets
- The display can be either a stick diagram or
simulated spectrum or sticks + simulated spectrum. Spectral simulation
can be for zeroth, first or second derivative of Gaussian or Lorentzian
lineshape function - example, and also inclusive of Doppler doubling for jet spectra
- Input of data in .ASR format of ASROT, which can come from ASROT itself, or from SPCAT, as converted from .OUT files with the
program PISORT
- Data files in .CAT format of SPCAT
can be read directly
- Batch input mode for hands-off input and
merging of many data sets
- Filtering of input data sets according to
frequency and intensity, and possibility of outputting the filtered data
- Different data sets can be displayed in
selectable linestyles and overall relative intensity
- Complete horizontal and vertical
scrollability
- Movable cursor displays parameters of lines
currently under the cursor (all six SPFIT/SPCAT quantum numbers per
energy level, transition frequency, intensity, lower level energy)
- Subsets of transitions can be highlighted on
the basis of several different definitions of commonality: the same
value of a given common quantum number, the same data set, or the same
transition type
- Two colour coding schemes for lines are
available. Colour can be according to transition type: various shades
of red for aP, aR-
and aQ-types, greens for b-types,
blues
for c-types, and purple for unrecognized transition
type. Alternatively, colouring can be according to the number of data
set
- Several selectable overall colour schemes
for the display
- Intensities of overlapped lines can either
be summed or not
- Help screen, summary screen of input,
scrollable listing screen for the ASCII data
- Output of current screen as files for the gle program which gives a route to PostScript printouts
- The program has a half-height option so that
ASCP and SVIEW
can be used to display predictions and spectra in the form of
synchronised strips for easier assignment work - example
- ranges and names of recorded spectra can be
indicated on the plot if the file LIST generated by SLIST is in the same directory as the other input files
Summaries and brief discussions of ASCP are given in J.Chem.Phys. 109,10263-10272(1998)
and
in J.Mol.Spectrosc. 199,5-12(2000),
which can serve as citations of the use of this program. Various other
examples of its use are referred to at the top of its listing.
 |
| |
| |
Current version |
| ASCP.FOR |
The listing - this is
the currently supported version of the program, tailored for
compilation with CVF6. The program uses Microsoft graphics and the issues
associated with this are discussed separately. |
| ASCP.EXE |
The CVF6
executable of the above: 100000 lines, 20 data files. Needs about
6.5Mbytes of memory, but runs happily enough on a P120/32Mb/W95 system.
Key properties of the graphics, i.e. window size and the display font
are now read from the file ASCP.CFG (below). The colour depth of the display should be at least
8 bits (256 colours).
Running conditions and installation
are identical to those discussed in connection with SVIEW
Drag and drop operation is possible
on all three input file types, i.e. .ASR, .CAT and .INP. The file can
be dragged onto the ASCP icon and, once file type option is set as required, the
dragged file will be read in and processed appropriately.
Irrespective of the operating
system the program assumes that the default directory for file names
specified in the batch input .INP file is the directory containing that file. The ECHO.ASC file, if
specified, will also end up in this directory.
|
| ASCP.CFG |
Configuration file
for ASCP, which has to reside
in the directory C:\ROT. |
| ASCP_KEYS.PDF |
Summary of command keys for ASCP.
The pairs of cursor keys for scaling and scrolling operations have been
selected on the basis of position on the QWERTY keyboard. The intended
ergonomy of use is perturbed on some other (e.g. German, French)
keyboards, in which case a possible remedy is to switch to the EN
keyboard while using this program.
A byproduct of rapid use of these control keys is
that the user can inadvertently press the Ctrl
Q key
combination,
which stops the program and it no longer responds. If this
happens press Ctrl
S.
|
| |
ASCP Picture Gallery |
| |
Snaps for current
graphics: |
| H2O2.GIF
H2O2.PS |
Display of the .CAT file for H2O2 from the JPL
catalog with a plot of both a stick
diagram and of a simulated spectrum. The first file is an ASCP screenshot and the
second contains PostScript output generated for this display with
option G and the gle package. |
| HCCCN.GIF |
Simultaneous display
of all catalogue files for HCCCN from the Cologne Database. The cursor
displays parameters of one of the J''=50 l-doublets for the
doubly degenerate v6=1 satellite, and all lines of this
satellite are highlighted. Colour coding of datasets is used. |
| HCCCN.INP |
The batch input file
for the above. |
| C2H3NC.GIF |
Display of the
catalogue file for vinyl isocyanide from the JPL
catalog (the file c053007.cat). The
J=41<-40 aR-type
band
is highlighted, while the cursor is positioned on a line of the K-1=5<-4 aQ
band |
| |
Snaps of screens
for older graphics versions: |
| MECN.GIF |
The J=17<-16 rotational
transition of methyl cyanide at 313 GHz and room temperature. The
ground state is highlighted, while the cursor is on the upper Kl-1=0
component
of the v8=1 doubly degenerate satellite. Note that
there are five lines under the cursor, due to the hyperfine structure
which is unresolved on this scale. |
| CB.GIF |
The J=3<-2 transition in the
FTMW spectrum of parent chlorobenzene at 1K with visible chlorine
hyperfine structure. The two K-1=2 transitions
have been highlighted, and the cursor is on the strongest K-1=0
hyperfine
component. |
| R22.GIF |
MMW spectrum of the
HCFC-22 molecule (CHF2Cl), with several isotopomers and
vibrational satellites read in. The marking option is used to select
the K=33 Q-branch - see if you can spot it in the
spectrum in Fig.1 of J.Mol.Spectrosc. 184,
150 (1997). |
| |
Legacy versions |
| ASCP_old.FOR |
Listing of the last
version prior to the use of CVF6 - this is no longer supported since
there are too many differences to keep a common version for all
compilers. This version also has some limits on the size of the quantum
numbers (basically two digits).
This listing and the executables
below are really only kept for use in emergencies.
|
| ASCP_M5.EXE |
Standard compilation
of the above (MSF5, 12000 lines, VGA graphics) - this will run on all versions
of DOS/WIN which allow full screen MS-DOS mode. The program requires
about 530 k of low DOS memory, if there are problems check this with
the MEM/C
command in the DOS window (second number from the bottom) |
| ASCP_PS1.EXE |
MSPS1 compilation of ASCP_old -
more lines (100000) and VGA graphics, requires the DOSXMSF.EXE 32-bit
extender. Run this in full screen mode in an MS-DOS window, in which
the detection of Windows from DOS programs has been blocked |
| ASC.FOR |
The much older
version, which displays six, and not twelve quantum numbers but is not
limited in their values |
| ASC.EXE |
MSF5 executable for the above,
requires the ANSI.SYS utility
of MS-DOS to be loaded. |
| |
ASCP input/output
examples |
| SAMPLE.INP |
Sample input control
file for reediting. Notes:
- the filenames are to be
specified without any path, ASCP assumes that the data files are in the same directory as the
.INP file
- .ASR and .CAT input files can be mixed
- do not change any of the
descriptive comments at the beginnings of lines as these are sometimes
used by the program when sorting out input problems
|
| KATALOG.PS |
Sample PostScript
output via gle, which demonstrates how ASCP can be used to browse through a suitable catalog of spectral
species. In this case the relative intensities for the various
molecules have been adjusted for comparable intensity. Relative to
standard output only the top comment line has been added, and legends
have been changed from file names to physically meaningful descriptors. |
|
|
 |
Back
to the table of programs
| PICKETT |
Programs augmenting
H.M.Pickett's SPFIT/SPCAT
|
|
Crib sheet to
SPFIT/SPCAT
|
These are several postprocessing programs to deal with output and data
associated with SPFIT/SPCAT. They should compile
with any sensible FORTRAN compiler, but PISORT really requires
a 32-bit compiler.
PIFORM
and PISORT are tested regularly against recent versions of SPFIT/SPCAT,
but they do get caught out now and again by some format changes. If
this is found to be the case please send me a note including the
relevant .OUT or .FIT files.
The Crib-sheet is a quick-reference guide to SPFIT/SPCAT arising from
experience gained on using these programs in the IFPAN laboratory.
 |
| |
| |
PIFORM |
| PIFORM.FOR |
- To reformat the .FIT output into a
form closer to the requirements for publishing:
- the values of constants and
errors are merged into the standard x.xxxx(xx) type notation
- it is ensured that a decimal
point is preceded by a leading zero, where necessary
- readable correlation matrix is
produced
- standard errors are determined
for constants
- NEW: it is possible to put descriptive annotations in the .LIN file, that will
be echoed to the output to produce a more informative printout. An
annotation can be placed at the end of a transition declaration line,
in the form of the character !, or this character followed by descriptive text. These
annotations will be transferred to the PIFORM output, in such
a way that a lone ! will produce an empty separating line, and !
with text will echo the text between separating lines.
As an example this input (.FIT from current SPFIT) or this input (.FIT from one of
the last FORTRAN versions of SPFIT), both lead to this output.
Warning: you should not have any BAD LINES
in the fit if
the scheme for transferring annotations from the .LIN file is to work
properly.
|
| PIFORM.EXE |
Executable for any Windows
compiled with CVF6
|
| |
PISLIN |
| PISLIN.FOR |
- This program can perform two
different types of operations on .LIN files used by SPFIT:
- checking for line duplication
and for components of blended lines existing in different parts of the
dataset, as well as some automatic corrections
- sorting of the .LIN file according
to any quantum number, change in quantum number, frequency and
magnitude of error.
Sorting operations can be cumulated
as many times as necessary by means of subsorting within the criterion
of an immediately preceding sort.
One possible task is to sort lines
into blocks for a given vibrational label, such that in each block
lines are sorted into bocks of P-, Q-, and R-type
transitions,
each such sub-block is sorted in Ka,
and each Ka sub-block is then sorted in J.
The
program only rearranges the ordering of .LIN file lines, so
formatting and annotations in these lines are preserved.
|
| PISLIN.EXE |
Executable for
any Windows compiled with CVF6
|
| |
PICHAM |
| PICHAM.FOR |
- This is a reformatting program
designed to make it easier to check the actual Hamiltonian matrix
elements set up for a given parameter in the .PAR file.
- PICHAM works on SPCAT output and the procedure is as follows:
- set up a test parameter file in
which values of all parameters are zeroed by setting them to a small
value such as 1.000000000000000E-025. The value of the parameter of interest should then be set
to unity. All parameters should be blocked from being modified by the
fit by setting the parameter uncertainty ERPAR for each
parameter to some even smaller value, e.g. 1.00000000E-037.
- run SPFIT to generate a .VAR file for SPCAT.
Alternatively the editing described above for the .PAR file could have
been done directly on the .VAR file
- run SPCAT by setting sorting EGYFLG in the .INT file to 5 to
obtain a Hamiltonian dump with no diagonalisation
- run PICHAM on the
resulting .EGY file. Output will be written to the corresponding .HAM file, which
contains an echo of the constants used for the calculation, followed by
Hamiltonian matrices in square form, block by block. Full state
labelling of up to six labels is suported.
The example files HTEST.VAR, HTEST.EGY, HTEST.HAM show you how you can check that in the prolate asymmetric
rotor Hamiltonian the A rotational constant has matrix
elements equal to K2. The E+, E-, O+, O-
factorisation is also nicely visible.
|
| PICHAM.EXE |
Executable for
any Windows compiled with CVF6
|
|
PISORT |
| PISORT.FOR |
- To extract and to sort
calculated transitions from the .OUT output from SPCAT. The lines are
sorted in frequency and output can be to either:
- frequency sorted set of
unmodified .OUT lines
- data file compatible with the .ASR standard as
used by ASROT. This can then be viewed graphically with ASCP (although note that ASCP can now read .CAT files directly)
Intensities can be
recalculated for a different specified temperature, but it is assumed
that the SPCAT output is for 300K. Some rudimentary filtering of the
data into subsets of interest is possible and can be extended.
|
| PISORT.EXE |
CVF6
executable for any Windows to process files with up to 100000 lines,
graphics is not used but a 32bit compiler is required because of
internal storage requirements.
|
| |
Some other programs |
| LINASF.FOR LINASF.EXE |
Converter of line information from
the standard in .LIN files to
that of .ASF files used by the ASFIT
program. All quantum numbers other than the three asymmetric top
quantum numbers per level are discarded.
|
| CATBAND.FOR CATBAND.EXE |
This is a program
similar to ASRBAND
and serves to place line sequence/band markers on diagrams of spectra
drawn with gle. Documentation is at the top of the
program listing.
|
| PISFIT.FOR |
To extract from the .FIT output those
lines which have obs-calc deviation in excess of a specified criterion
- useful for large data sets. A similar feature is currently also
available in PIFORM
|
| |
CRIB sheet |
| CRIB.HTM |
This is a Crib sheet for
SPFIT/SPCAT,
which may be used in addition to HMP's documentation
contained in http://spec.jpl.nasa.gov/ftp/pub/calpgm/spinv.pdf
Please note that another
clarification effort is available on the Cologne group website at http://www.ph1.uni-koeln.de/cdms/pickett/index.html
|
|
|
 |
Back
to the table of programs
| CONVC |
CONVersions
of
asymmetric rotor Constants
|
| PLANM |
PLANar
Moments etc..
|
CONVC
carries out the following conversions of asymmetric rotor constants:
- Watson Ir,S to
Watson Ir,A
- Watson Ir,A to
Watson Ir,S
- vanEijck,Typke to Watson A
- Kirchhoff to Watson A
- Watson Ir,A to
Watson IIIr,A
Some
of these conversions work up to sextics, and some convert errors. This
program has never really been finished and has been modified as
necessary.
PLANM is
just a handy calculator for some quantities derivable from rotational
constants. Input is from the keyboard, output is to screen and is
echoed to file PLANM.OUT.
 |
| |
| CONVC.FOR |
The
listing. The data is to be in the .CON format as used by ASROT and
ASFIT. Main output is to screen and is also appended to
the file CONVC.OUT. Additional output in the .CON format is appended to CONVC.CON. |
| CONVC.EXE |
Windows executable |
| TRI.CON |
Sample constants file for use with
CONVC |
| CONVC.OUT |
Output for Watson's A-S
conversion of the above, compare with refitted values in Table 2 of J.Mol.Spectrosc.
178,125(1996) |
| |
|
| PLANM.FOR |
The listing |
| PLANM.EXE |
Windows executable |
| PLANM.OUT |
Sample output |
|
|
 |
Back
to the table of programs
| PLAN |
Quadratic, quartic and
sextic PLANarity relations for an asymmetric
top
|
This
program explores the validity of various order planarity relations in
Watson's reduced Hamiltonian (reduction-A, representation Ir),
using
the formulae from Watson, J.Mol.Spectrosc. 65,
123
(1977).
Input requires sextic-level constants to be in a file in the standard
as written by ASFIT. PLAN will read errors in constants and propagate them
appropriately into errors in the planarity defects Di, Dq, Ds.
The output is in the form as used in a comparison of such moments for
several molecules, made in Table 3 of Kisiel, Pszczolkowski, Z.Naturforsch.
50a, 347 (1994). Output is written
to screen and is also appended to file PLAN.OUT.
 |
| |
| PLAN.FOR |
The
listing. |
| PLAN.EXE |
An executable for any DOS/Windows. |
| TRI.CON |
Constants file for ground state
trichloroethylene from J.Mol.Spectrosc. 178,125(1996). |
| PLAN.OUT |
The output for the above. Note
that the sextic defect is not determined with significance even though
the sextics are quite precise. |
|
|
 |
Back
to the table of programs
|