A general code for fitting global potential energy surfaces via CHIPR method: Direct-Fit Diatomic and tetratomic molecules

Abstract

A general program to fit global adiabatic potential energy surfaces of up to tetratomic molecules to ab initio points and available spectroscopic data for simple diatomics is reported. It is based on the Combined-Hyperbolic-Inverse-Power-Representation (CHIPR) method. The final form describes all dissociating fragments and long-range/valence interactions, while obeying the system permutational symmetry. The code yields as output a Fortran 90 subroutine that readily evaluates the potential and gradient at any arbitrary geometry. Program summaryProgram title: CHIPR-4.0CPC Library link to program files:http://dx.doi.org/10.17632/8wdv87gt5x.2Licensing provisions: GPLv3Programming language: Fortran 90Journal reference of previous version: C. M. R. Rocha, A. J. C. Varandas, Comput. Phys. Commun. 247 (2020) 106913Does the new version supersede the previous version?: YesReasons for the new version: CHIPR-4.0 is a major extension of the previous CHIPR-3.0 code [1] which, besides diatomics and triatomics, also includes the possibility of fitting global adiabatic potential energy surfaces (PESs) of A4-, AB3-, A2B2-, ABC2- and ABCD-type tetratomic molecules [3]. Additionally, a new feature is that it allows the user to fine-tune ab initio diatomic curves using available spectroscopic data.Summary of revisions:• Implementation of CHIPR’s polynomial form and permutation operators for tetratomic molecules• Implementation of a direct-fit module to fine-tune ab initio potential energy curves using spectroscopic data• Automatic global minimum search on the diatomic, triatomic or tetratomic final PESs• Automatic harmonic vibrational analysis at the optimum molecular geometriesNature of problem: This version of the CHIPR code, CHIPR-4.0, provides a set of subroutines to fit global adiabatic potential energy surfaces of up to tetratomic molecules using ab initio and (for diatomics) optional fine-tuning experimental data.Solution method: CHIPR-4.0 uses the Combined-Hyperbolic-Inverse-Power-Representation (CHIPR) [2] method to interpolate and extrapolate ab initio data points.Additional comments including restrictions and unusual features: For triatomic and tetratomic fits, the user must supply subroutines containing the analytic forms defining the sum of two- and two-plus-three-bodies, respectively. These are utilized internally to calculate (from ab initio data) the corresponding three- and four-body energies to be fitted. The current version of this code is strictly based on single-sheeted analytical PESs. High-accuracy in diatomic refinements is currently limited to the 1Σ single-curve case.