Numerically constructed internal-coordinate Hamiltonian with Eckart embedding and its application for the inversion tunneling of ammonia

It is shown that the use of an Eckart-frame embedding with a kinetic energy operator expressed in curvilinear internal coordinates becomes feasible and straightforward to implement for arbitrary molecular compositions and internal coordinates if the operator is defined numerically over a (discrete variable representation) grid. The algorithm proposed utilizes the transformation method of Dymarsky and Kudin to maintain the rotational Eckart condition. In order to demonstrate the applicability and flexibility of our approach the non-rigid ammonia molecule is considered and the corresponding rotational–vibrational energy levels and wave functions are computed using kinetic energy operators with three different embeddings. Two of them fulfill the Eckart conditions corresponding to a trigonal pyramidal ( C 3v ) and a trigonal planar ( D 3h ) reference structure and the third one is a non-Eckart frame. The computed energy levels are, of course, identical, and the structure of the three different wave-function representations are analyzed in terms of the rigid rotor functions for a symmetric top. The possible advantages of one frame representation over another are discussed concerning the interpretation of the rovibrational states in terms of the traditional rigid rotor labels.