**Take a sneak peek at the new NIST.gov and let us know what you think!***(Please note: some content may not be complete on the beta site.)*.

NIST Home > MML > Chemical Sciences Division > Chemical Informatics Research Group > Fourier Grid Hamiltonian 1D Program Interface

The program interface is divided up into a series of frames and buttons: **Data Frame****Range Frame****Mass Frame****Parameter Frame****Start Button****Thermo Button****Status Frame****Expectation Values Frame****Fitting Frame****Eigenvalues Frame****Plot Window**
The data frame provides the user with a choice of predetermined potential energy functions or reads the potential from a data file or the windows clipboard. The choices are: 1.a. Morse Potential
V(r) = D D E(v) = w e (v+1/2) - w v is the vibrational quantum number, w w h is Planck's constant, c is the speed of light and m is the mass of the particle in the well. For w w w w For the same units as above this can be written as: w See "Molecular Spectra and Molecular Structure 1. Spectra of Diatomic Molecules" G. Herzberg, 1950, D. Van Nostrand Co., New York, for further details. The constants D An approximation to the vibrational levels of HF can be described using Morse parameters: r A range of 0.5 to 2.5 Å and 60 points provides a good display. The first few energy levels should be:
This is actually used to describe any potential of the form: V(r) = C The coefficients C This can be used to describe potentials with a double minimum such as the inversion motion in ammonia (NH C A range of -1.0 to +1.0 Å and 60 points provides a good display. The first few energy levels should be:
This can be used to describe potentials of the form: V(x) = 1/2(V For the cos(nx) potential the range is set automatically to be 0 to 2p The cos(nx) potential is useful for describing internal rotations. An approximation of the Ethane (CH V mass (moment of inertia) = 1.577 amu Å The range is set automatically. 60 points provide a good display. The first few energy levels should be:
One of the nice features of the FGH method is that it has periodic boundary conditions. This makes it ideal for internal rotations with finite barriers, as the wavefunctions naturally go over to free rotation wavefunctions above the barrier.
This allows the program to get data from an external text file. The x and y values can be delimited with spaces, commas or tabs. An example text file looks like:
The first line which is the number of points. Subsequent lines are x y pairs in units of Å for x and cm The program does not understand scientific notation. The program reports the range and the number of points in the file in the Parameters frame. The program performs a cubic spline through the points to generate the points used by the FGH algorithm. The number of points and range can be adjusted in the Range frame. Setting the range outside of the range of points in the data file can lead to random potentials as the spline fit may not extrapolate well.
This allows the program to get data from the Windows clipboard. The data on the clipboard should be x y pairs, similar to the format for the Read File selection except the first line should not be the number of points.
This potential does not work well with the FGH method as the grid points used to define the basis set do not usually line up with the walls of the square well potential. This is in effect an error in the width of the well, which leads to a corresponding error in the energy levels. This frame allows the user to set the range of x used for the calculation. The units are Å for most potentials and radians for torsional potentials (in which case the mass should be entered as a moment of inertia in units of amu Å The Range frame also includes a checkbox to force the range to be from 0 to 2p(n-1)/n, where n is the number of grid points. This is useful for torsional potentials read in from a file or the clipboard. This frame allows the user to set the mass. The units are amu, or if an angle rather than a distance is used for the x coordinate then the units are amu Å The parameter frame allows the user to enter the parameters for the various predertermined potentials. The parameter frame also displays the number of points for potentials read from a file or the clipboard, as well as the minimum x and maximum x. This button starts the calculation. 60 points on a 400 MHz pc takes about 3 seconds. The progress of the calculation is reported in the Status Frame. When the calculation is complete the plot window is opened and displays the potential, the eigenvalues, and the eigenvectors. The eigenvalues are also listed in the eigenvalues frame. The Thermo button is used after the eigenvalues have been computed. The Thermo button opens the thermodynamics window and uses the eigenvalues to calculate the contribution to entropy, heat capacity and integrated heat capacity(0K to chosen temperature). The thermo window allows the user to choose a temperature, the number of eigenvalues to use, and set a symmetry number. As well as entropy, heat capacity and integrated heat capacity, the window also displays the partition function and the contribution to the partition function of the highest energy eigenvalue used. The status frame displays the progress of the calculation. Generating the Hamiltonian is generally the fastest. Due to the use of delta functions as a basis set, there are only sums, no integrals to be evaluated. Once the hamiltonian is generated and the matrix set up, the matrix is diagonalized in two steps - first tridiagonalized, and then the program solves for the eigenvalues and eigenvectors using the tridiagonal matrix. This last step is the most time consuming. This frame lists the expectations values <x> and <x This frame calculates a polynomial fit to the eigenvalues. It uses the form E(v) = T The eigenvalues frame displays v, the quantum number, the eigenvalues, and their differences. If the Fitting frame has been used the results of the fit are also displayed. The The plot Window graphically displays the potential, the eigenvalues (as horizontal lines), and the eigenvectors as colored curves. The eigenvectors are not displayed by default. To view them check the |
## Related Links:Contact
100 Bureau Drive, M/S 8320 |