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