Interaction energies for the water dimer by supermolecular methods and symmetry-adapted perturbation theory: the role of bond functions and convergence of basis subsets

 Using a systematic series of basis sets in supermolecular and symmetry-adapted intermolecular perturbation theory calculations it is examined how interaction energies of various water dimer structures change upon addition and shifting of bond functions. Their addition to augmented double- and triple-zeta basis sets brings the sum of the electron correlation contributions to the second-order interaction energy nearly to convergence, while accurate first-order electrostatic and exchange contributions require better than augmented quadruple-zeta quality. A scheme which combines the different perturbation energy contributions as computed in different basis subsets performs uniformly well for the various dimer structures. It yields a symmetry-adapted perturbation theory value of −21.08 kJ/mol for the energy of interaction of two vibrationally averaged water molecules compared to −21.29 kJ/mol when the full augmented triple-zeta basis set is used throughout. Received: 4 November 1999 / Accepted: 8 February 2000 / Published online: 12 May 2000     

Restricted Hartree-Fock and unrestricted Hartree-Fock as reference states in many-body perturbation theory: a critical comparison of the two approaches

Summary The paper deals with two topics related to the problem which reference state is better for many-body perturbation theory: restricted Hartree-Fock (RHF) or unrestricted Hartree-Fock (UHF)? The first topic concerns the potential surfaces. Several examples are presented to show shortcomings of the two approaches and a simple way is presented which seems to give a useful potential curve in the whole range of interatomic distances by a composition of RHF and UHF potential curves. The second topic concerns the many-body perturbation theory for open-shell systems in the RHF formalism. The method is critically examined and compared with the ordinary many-body perturbation theory using UHF as the reference. This examination of many-body techniques provides also some insight into the problems inherent of the SCF theory: spin contamination from higher multiplets, localization of orbitals, and self-consistency effects.     

Intruder states in multireference perturbation theory: the ground state of manganese dimer

A detailed analysis of a severe intruder state problem in the multistate multireference perturbation theory (MS-MRPT) calculations on the ground state of manganese dimer is presented. An enormous number of detected intruder states (> 5000) do not permit finding even an approximate shape of the X(1)Sigma(g) (+) potential energy curve. The intruder states are explicitly demonstrated to originate from quasidegeneracies in the zeroth-order Hamiltonian spectrum. The electronic configurations responsible for appearance of the quasidegeneracies are identified as single and double excitations from the active orbitals to the external orbitals. It is shown that the quasidegeneracy problem can be completely eliminated using shift techniques despite of its severity. The resultant curves are smooth and continuous. Unfortunately, strong dependence of the spectroscopic parameters of the X(1)Sigma(g) (+) state on the shift parameter is observed. This finding rises serious controversies regarding validity of employing shift techniques for solving the intruder state problem in MS-MRPT. Various alternative approaches of removing intruder states (e.g., modification of the basis set or changing the active space) are tested. None of these conventional techniques is able to fully avoid the quasidegeneracies. We believe that the MS-MRPT calculations on the three lowest A(g) states of manganese dimer constitute a perfect benchmark case for studying the behavior of MRPT in extreme situations.     

An application of second-order n-electron valence state perturbation theory to the calculation of excited states

n–electron valence state perturbation theory (NEVPT) is a form of multireference perturbation theory where all the zero-order wave functions are of multireference nature, being generated as eigenfunctions of a two–electron model Hamiltonian. The absence of intruder states makes NEVPT an interesting choice for the calculation of electronically excited states. Test calculations have been performed on several valence and Rydberg transitions for the formaldehyde and acetone molecules; the results are in good accordance with the best calculations and with the existing experimental data.Contribution to the Jacopo Tomasi Honorary Issue     

Information theory of D‐dimensional hydrogenic systems: Application to circular and Rydberg states

The analytic information theory of quantum systems includes the exact determination of their spatial extension or multidimensional spreading in both position and momentum spaces by means of the familiar variance and its generalization, the power and logarithmic moments, and, more appropriately, the Shannon entropy and the Fisher information. These complementary uncertainty measures have a global or local character, respectively, because they are power‐like (variance, moments), logarithmic (Shannon) and gradient (Fisher) functionals of the corresponding probability distribution. Here we explicitly discuss all these spreading measures (and their associated uncertainty relations) in both position and momentum for the main prototype in D‐dimensional physics, the hydrogenic system, directly in terms of the dimensionality and the hyperquantum numbers which characterize the involved states. Then, we analyze in detail such measures for s‐states, circular states (i.e., single‐electron states of highest angular momenta allowed within an electronic manifold characterized by a given principal hyperquantum number), and Rydberg states (i.e., states with large radial hyperquantum numbers n). © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Nature and importance of three-body interactions in the (H2O)2HCl trimer

 The nature and importance of nonadditive three-body interactions in the (H₂O)₂HCl cluster have been studied by the supermolecule coupled-cluster method and by symmetry-adapted perturbation theory (SAPT). The convergence of the SAPT expansion was tested by comparison with the results obtained from the supermolecule coupled-cluster calculations including single, double, and noniterative triple excitations [CCSD(T)]. It is shown that the SAPT results reproduce the converged CCSD(T) results within 3% at worst. The SAPT method has been used to analyze the three-body interactions for various geometries of the (H₂O)₂HCl cluster. It is shown that the induction nonadditivity is dominant, but it is partly quenched by the first-order Heitler–London-type exchange and higher-order exchange–induction/deformation terms. This implies that the classical induction term alone is not a reliable approximation to the nonadditive energy and that it will be difficult to approximate the three-body potential for (H₂O)₂HCl by a simple analytical expression. The three-body energy represents as much as 21–27% of the pair CCSD(T) intermolecular energy. Received: 15 September 1999 / Accepted: 3 February 2000 / Published online: 2 May 2000     

Cyanovinyl radical: an illustration of the poor performance of unrestricted perturbation theory and density functional theory procedures in calculating radical stabilization energies

Stabilization energies for the 1-cyanovinyl radical (CH₂=CCN) have been calculated using a variety of conventional ab initio (M?ller–Plesset, quadratic configuration interaction and coupled-cluster) and density functional theory (B-LYP, B3-LYP) procedures, as well as with a range of compound methods. Compared with a high-level benchmark value (that predicts a stabilization energy of 17.1 kJ mol⁻¹), UMP2 and UMP4 give the wrong sign and magnitude of the stabilization energy (both methods predicting desta- bilization instead of stabilization), while B-LYP and B3-LYP overestimate the degree of stabilization. The RMP2, RMP4, QCISD(T) and CCSD(T) techniques, and several, but not all, variants of G2 and CBS theories give radical stabilization energies in good agreement with the benchmark value. Received: 15 June 1998 / Accepted: 19 August 1998 /  Published online: 15 February 1999     

Generalization of coupled-cluster response theory to multireference expansion spaces: application of the coupled-cluster singles and doubles effective Hamiltonian

Employing separate cluster ansatz in time-independent and time-dependent wave-operators, coupled-cluster (CC) response theory is generalized to multireference (MR) expansion spaces. For state energies, this corresponds to the MR secular problem with an arbitrary similarity-transformed effective Hamiltonian, H˜=Ω⁻¹ HΩ. The effective Hamiltonian can be generated via size-extensive CC methods. Thus the states in MR linear response theory (MRLRT) maintain the usual CC core-extensive properties. We have used the Gelfand unitary group basis of the spin-adapted configurations to construct the matrix of H˜ in the MR excitation space. As a preliminary application, the CC singles and doubles effective Hamiltonian is applied to excitation and photoionization energies of the CH⁺ and N₂ molecules, and is compared with experimental results and results from other numerical procedures including conventional CC linear response theory (CC-LRT), MR and full configuration interaction (MRCI and FCI) methods. The numerical results indicate that MRLRT reproduces valence and external excited states quantitatively, combining the best features of CC-LRT and MRCI. Received: 2 July 1998 / Accepted: 28 August 1998 / Published online: 11 November 1998     

Theoretical study of the π→π* excited states of linear polyenes: The energy gap between 11Bu+ and 21Ag− states and their character

Multireference perturbation theory with complete active space self-consistent field (CASSCF) reference functions is applied to the study of the valence π→π* excited states of 1,3-butadiene, 1,3,5-hexatriene, 1,3,5,7-octatetraene, and 1,3,5,7,9-decapentaene. Our focus was put on determining the nature of the two lowest-lying singlet excited states, 1¹B_u⁺ and 2¹A_g⁻, and their ordering. The 1¹B_u⁺ state is a singly excited state with an ionic nature originating from the HOMO→LUMO one-electron transition while the covalent 2¹A_g⁻ state is the doubly excited state which comes mainly from the (HOMO)²→(LUMO)² transition. The active-space and basis-set effects are taken into account to estimate the excitation energies of larger polyenes. For butadiene, the 1¹B_u⁺ state is calculated to be slightly lower by 0.1 eV than the doubly excited 2¹A_g⁻ state at the ground-state equilibrium geometry. For hexatriene, our calculations predict the two states to be virtually degenerate. Octatetraene is the first polyene for which we predict that the 2¹A_g⁻ state is the lowest excited singlet state at the ground-state geometry. The present theory also indicates that the 2¹A_g⁻ state lies clearly below the 1¹B_u⁺ state in decapentaene with the energy gap of 0.4 eV. The 0–0 transition and the emission energies are also calculated using the planar C_2h relaxed excited-state geometries. The covalent 2¹A_g⁻ state is much more sensitive to the geometry variation than is the ionic 1¹B_u⁺ state, which places the 2¹A_g⁻ state significantly below the 1¹B_u⁺ state at the relaxed geometry. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 157–175, 1998     

Generalized vibrational perturbation theory for rotovibrational energies of linear,symmetric and asymmetric tops: Theory,approximations, and automated approaches to deal with medium‐to‐large molecular systems

Models going beyond the rigid‐rotor and the harmonic oscillator levels are mandatory for providing accurate theoretical predictions for several spectroscopic properties. Different strategies have been devised for this purpose. Among them, the treatment by perturbation theory of the molecular Hamiltonian after its expansion in power series of products of vibrational and rotational operators, also referred to as vibrational perturbation theory (VPT), is particularly appealing for its computational efficiency to treat medium‐to‐large systems. Moreover, generalized (GVPT) strategies combining the use of perturbative and variational formalisms can be adopted to further improve the accuracy of the results, with the first approach used for weakly coupled terms, and the second one to handle tightly coupled ones. In this context, the GVPT formulation for asymmetric, symmetric, and linear tops is revisited and fully generalized to both minima and first‐order saddle points of the molecular potential energy surface. The computational strategies and approximations that can be adopted in dealing with GVPT computations are pointed out, with a particular attention devoted to the treatment of symmetry and degeneracies. A number of tests and applications are discussed, to show the possibilities of the developments, as regards both the variety of treatable systems and eligible methods. © 2015 The Authors International Journal of Quantum Chemistry Published by Wiley Periodicals, Inc.