Computing resonance energies, widths, and wave functions using a Lanczos method in real arithmetic

We introduce new ideas for calculating resonance energies and widths. It is shown that a non-Hermitian-Lanczos approach can be used to compute eigenvalues of H+W, where H is the Hamiltonian and W is a complex absorbing potential (CAP), without evaluating complex matrix-vector products. This is done by exploiting the link between a CAP-modified Hamiltonian matrix and a real but nonsymmetric matrix U suggested by Mandelshtam and Neumaier [J. Theor. Comput. Chem. 1, 1 (2002)] and using a coupled-two-term Lanczos procedure. We use approximate resonance eigenvectors obtained from the non-Hermitian-Lanczos algorithm and a very good CAP to obtain very accurate energies and widths without solving eigenvalue problems for many values of the CAP strength parameter and searching for cusps. The method is applied to the resonances of HCO. We compare properties of the method with those of established approaches.     

Towards fast computations of correlated vibrational wave functions: vibrational coupled cluster response excitation energies at the two-mode coupling level

An efficient implementation of vibrational coupled cluster theory with two-mode excitations and a two-mode Hamiltonian is described. The algorithm is shown to scale cubically with respect to the number of modes which is identical to the scaling of the corresponding vibrational configuration interaction algorithm. This is achieved through the use of special intermediates. The same algorithm can also be used in vibrational M?ller-Plesset calculations. To improve performance, screening techniques have been implemented as well. Test calculations on polyaromatic hydrocarbons with up to 264 coupled modes and model systems with up to 1140 modes are used to illustrate the various features of the algorithm.     

Vibrational absorption spectra from vibrational coupled cluster damped linear response functions calculated using an asymmetric Lanczos algorithm

We report the theory and implementation of vibrational coupled cluster (VCC) damped response functions. From the imaginary part of the damped VCC response function the absorption as function of frequency can be obtained, requiring formally the solution of the now complex VCC response equations. The absorption spectrum can in this formulation be seen as a matrix function of the characteristic VCC Jacobian response matrix. The asymmetric matrix version of the Lanczos method is used to generate a tridiagonal representation of the VCC response Jacobian. Solving the complex response equations in the relevant Lanczos space provides a method for calculating the VCC damped response functions and thereby subsequently the absorption spectra. The convergence behaviour of the algorithm is discussed theoretically and tested for different levels of completeness of the VCC expansion. Comparison is made with results from the recently reported [P. Seidler, M. B. Hansen, W. Gyo?rffy, D. Toffoli, and O. Christiansen, J. Chem. Phys. 132, 164105 (2010)] vibrational configuration interaction damped response function calculated using a symmetric Lanczos algorithm. Calculations of IR spectra of oxazole, cyclopropene, and uracil illustrate the usefulness of the new VCC based method.     

Accurate estimation of correlation energies using locally dense basis sets

Locally dense basis sets using the mixed 6-311G(d, p)/3-21G basis can be used to reproduce total energies and correlation energies after empirical adjustment to 2–4 kcal/mol for a variety of small and medium size molecules containing hydrog en, carbon, and oxygen. Post-Hartree\–Fock methods can be calculated faster by this method by factors of 2–3, in general, and higher in the presence of high molecular symmetry; density functional approaches take longer and are impractical in th e locally dense basis set approach. It is shown that the correlation energy in two of the better characterized density functional approaches is generally significantly larger than that of the post-Hartree–Fock treatments studied here and appears to be insensitive to the basis set employed. © 1996 by John Wiley & Sons, Inc.     

Accurate calculations of intermolecular interaction energies using explicitly correlated wave functions

Explicitly correlated second-order M?ller-Plesset (MP2-F12) calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented and compared with standard MP2 results. The MP2 complete basis set limits are estimated using basis set extrapolation and augmented quadruple-zeta and quintuple-zeta basis sets. Already with augmented double-zeta basis sets the MP2-F12 interaction energies are found to be closer to the complete basis set limits than standard MP2 calculations with augmented quintuple-zeta basis sets. Various possible approximations in the MP2-F12 method are systematically tested. Best results are obtained with localized orbitals and the diagonal MP2-F12/C(D) ansatz. Hybrid approximations, in which some contributions of the auxiliary basis set are neglected and which considerably reduce the computational cost, have a negligible effect on the interaction energies. Also the orbital-invariant fixed-amplitude approximation of Ten-no leads to only slightly less accurate results. Preliminary results for the neon and benzene dimers, obtained with the recently proposed CCSD(T)-F12a approximation, indicate that the CCSD(T) basis set limits can also be very closely approached using augmented triple-zeta basis sets.     

Calculation of vibrational infrared intensities and Raman activities using explicit anharmonic wave functions

Methods for automatic computation of IR intensities and Raman activities are described using vibrational self-consistent field (VSCF) and vibrational configuration interaction (VCI) wave functions. Inclusion of effects due to anharmonicity in the potential energy and property surfaces are found to improve the results substantially as compared to experimental data. Sample calculations employing water and formaldehyde are presented, allowing for comparison between different vibrational methods. The convergence with respect to excitation level in VCI and the extent of mode coupling in the potential and property expansions is investigated. In addition, different electronic methods used for generating the potential and property surfaces, namely CCSD, CCSD(T), DFT/B3LYP, and DFT/CAM-B3LYP have been compared. Details of the potential and property surfaces may have significant effects on the IR and Raman intensities.     

Gaussian integral transforms of unconventional functions: contracted basis functions appropriate to molecules and metallic clusters

We report the development of new gaussian integral transforms one-electron functions. A straightforward recipe to generate gaussian contractions of these functions is introduced. The applications call into question the disseminated belief of the superiority of Slater type orbitals (STO) in molecular calculations. A first successful application of the Bessel K _v (qr) function to a molecular system is presented. One of the integral transforms that yields the 0s function is used to generate contractions that have a better performance than the STO ones in cases that diffuse functions are needed. Applications are presented to some Li _n clusters and a new conformation of the Li₃ cluster is then proposed.     

Improving the energies of approximate wave functions using the concepts of density functional theory

A simple method, derived from DFT, for improving the energy of a trial density is modified for the case of atoms. It is assumed that errors in the interelectronic repulsion are the only significant ones. The errors in the energies of single‐ζ Slater orbitals for the atoms from He to Ne are reduced an average factor of 21. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Calculating energy levels of isomerizing tetra-atomic molecules. II. The vibrational states of acetylene and vinylidene

A general, full-dimensional computational method for the accurate calculation of rotationally and vibrationally excited states of tetra-atomic molecules is further developed. The resulting computer program may be run in serial and parallel modes and is particularly appropriate for molecules executing wide-amplitude motions and isomerizations. An application to the isomerizing acetylene/vinylidene system is presented. Large-scale calculations using a coordinate system based on orthogonal satellite vectors have been performed in six dimensions and vibrational term values and wave functions for acetylene and vinylidene states up to approximately 23 000 cm(-1) above the potential minimum have been determined. This has permitted the characterization of acetylene and vinylidene states at and above the isomerization barrier. These calculations employ more extensive vibrational basis sets and hence consider a much higher density of states than in any variational calculations reported hitherto for this system. Comparison of the calculated density of states with that determined empirically suggests that our calculations are the most realistic achieved for this system to date. Indeed more states have been converged than in any previous study of this system. Calculations on lower lying excited states of acetylene based on HC-CH diatom-diatom coordinates give nearly identical results to those based on orthogonal satellite vectors. Comparisons are also made with calculations based on HH-CC diatom-diatom coordinates.     

Vibrational structure theory: new vibrational wave function methods for calculation of anharmonic vibrational energies and vibrational contributions to molecular properties

A number of recently developed theoretical methods for the calculation of vibrational energies and wave functions are reviewed. Methods for constructing the appropriate quantum mechanical Hamilton operator are briefly described before reviewing a particular branch of theoretical methods for solving the nuclear Schr?dinger equation. The main focus is on wave function methods using the vibrational self-consistent field (VSCF) as starting point, and includes vibrational configuration interaction (VCI), vibrational M?ller-Plesset (VMP) theory, and vibrational coupled cluster (VCC) theory. The convergence of the different methods towards the full vibrational configuration interaction (FVCI) result is discussed. Finally, newly developed vibrational response methods for calculation of vibrational contributions to properties, energies, and transition probabilities are discussed.     

On expectation value calculations of one-electron properties using the coupled cluster wave functions

The ability of various approximate coupled cluster (CC) methods to provide accurate first-order one-electron properties calculated as expectation values is theoretically analysed and computationally examined for BH and CO. For actual calculations the infinite number of terms of the expectation value expansion (O=¦exp (T ⁺)O exp (T)¦_c) was truncated so that T ₁ T ₂, T ₃, and (1/2) T ₂T₂ clusters were retained on both sides of O. The role of individual clusters is carefully discussed. Inclusion of T ₁, is unavoidable, but if triples are essential in the energy evaluation, they may play an even more important role in the property expansion, as shown in the case of CO. It is shown that the CC wave function, which is exact to second order, effectively satisfies the Hellmann-Feynman theorem.     

Two-electron reduced density matrices from the anti-Hermitian contracted Schrodinger equation: enhanced energies and properties with larger basis sets

Two-electron reduced density matrices (2-RDMs) have recently been directly determined from the solution of the anti-Hermitian contracted Schrodinger equation (ACSE) to obtain 95%-100% of the ground-state correlation energy of atoms and molecules, which significantly improves upon the accuracy of the contracted Schrodinger equation (CSE) [D. A. Mazziotti, Phys. Rev. Lett. 97, 143002 (2006)]. Two subsets of the CSE, the ACSE and the contraction of the CSE onto the one-particle space, known as the 1,3-CSE, have two important properties: (i) dependence upon only the 3-RDM and (ii) inclusion of all second-order terms when the 3-RDM is reconstructed as only a first-order functional of the 2-RDM. The error in the 1,3-CSE has an important role as a stopping criterion in solving the ACSE for the 2-RDM. Using a computationally more efficient implementation of the ACSE, the author treats a variety of molecules, including H2O, NH3, HCN, and HO3-, in larger basis sets such as correlation-consistent polarized double- and triple-zeta. The ground-state energy of neon is also calculated in a polarized quadruple-zeta basis set with extrapolation to the complete basis-set limit, and the equilibrium bond length and harmonic frequency of N2 are computed with comparison to experimental values. The author observes that increasing the basis set enhances the ability of the ACSE to capture correlation effects in ground-state energies and properties. In the triple-zeta basis set, for example, the ACSE yields energies and properties that are closer in accuracy to coupled cluster with single, double, and triple excitations than to coupled cluster with single and double excitations. In all basis sets, the computed 2-RDMs very closely satisfy known N-representability conditions.     

Hartree–Fock wave functions with a modified GTO basis for atoms

We have solved the atomic Hartree–Fock equations by using the algebraic approach, expanding the single-particle radial wave function in terms of a modified Gaussian type orbitals (GTOs) basis. Several atomic properties such as Kato's cusp condition for the electron density or the correct asymptotic behavior of the electron momentum density distribution are accurately verified. Additionally the energy of the atomic ground state can be obtained by using a smaller number of basis functions than in standard GTO expansions. This study has been performed for several atoms of the first three rows. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 65 : 59–64, 1997     

Theory of the expansion of wave functions in a gaussian basis

The convergence properties of the expansions of (a) the function 1/r and (b) the function exp(-αr) in an even-tempered basis of Gaussians are studied analytically. The starting points are the Gaussian integral representations of 1/r and exp(?αr). One arrives at an expansion in a finite number of Gaussians in three steps: (1) a restriction of the integration domain, (2) a variable transformation, and (3) discretization of the integral. The cutoff error goes in both cases essentially as exp(?ah), and the discretization error, as exp(?b/h). The minimum overall error is reached for the β-parameter of an even-tempered basis β ∽ exp(c/√n), where n is the dimension of the basis, and the error itself decreases as ? ∽ exp(?d√n). Different optimum basis parameters are obtained depending on which quantity one wants to minimize, e.g., the error of the energy expectation value, the distance in Hilbert space, the variance of the energy, or the density at the nucleus. © 1994 John Wiley & Sons, Inc.     

Relativistic effects determined using the Douglas-Kroll contracted basis sets and correlation consistent basis sets with small-core relativistic pseudopotentials

The coupled cluster approximation with single, double, and quasiperturbative triple excitations [CCSD(T)] was used in combination with the Douglas-Kroll contracted correlation consistent basis sets [cc-pVnZ-DK, where n = D(2), T(3), Q(4), and 5] and small-core relativistic pseudopotentials (PP) with correlation consistent polarized valence basis sets (cc-pVnZ-PP and aug-cc-pVnZ-PP) to investigate the impact of scalar relativistic corrections on energetic and structural properties of small molecules containing third-row (Ga-Kr) atoms. These molecules were taken from the Gaussian-2 extended test set for third-row atoms. Atomization energies, ionization energies, electron affinities, and proton affinities for molecules in the test set were determined and compared with nonrelativistic results which were obtained in a recent study in which the standard and augmented correlation consistent basis sets were used in combination with CCSD(T). Several schemes were used to extrapolate the energies to the complete basis set limit.     

The Cotton-Mouton effect of neon and argon: a benchmark study using highly correlated coupled cluster wave functions

The Cotton-Mouton effect (magnetic field induced linear birefringence) has been studied for neon and argon using state-of-the-art coupled cluster techniques. The coupled cluster singles, doubles and triples (CCSDT) approach has been used to obtain static benchmark results and the CC3 model with an approximate treatment of triple excitations to obtain frequency-dependent results. In the case of neon the effect of excitations beyond triples has also been estimated via coupled cluster calculations including quadruple excitations (CCSDTQ), pentuple excitations (CCSDTQP), etc. up to the full configuration-interaction level. The results obtained for the anisotropy of the hypermagnetizability Deltaeta(omega), the molecular property that determines the magnetic field induced birefringence of spherically symmetric systems, are Deltaeta=2.89 a.u. for neon and Deltaeta=24.7 a.u. for argon, with a negligible effect of frequency dispersion. For neon we could estimate an absolute error on Deltaeta of 0.1 a.u. The accuracy of these results surpasses that of recently reported experimental data.     

Long range interaction energies using Gaussian basis sets and a one center method

A one center method, based on the work of Karplus and Kolker, is discussed and used to calculate the induction energy, through O(R^?8), for the H(ls) – H⁺ interaction employing two types of Gaussian basis sets constructed from functions of the form {r^je^?αr2}. The effective hydrogen atom excitation energies and transition multipole moment matrix elements generated in these calculations are used to calculate the dispersion energy for the H(ls) – H(ls) interaction, through O(R^?10), and the R^?9 triple dipole energy corresponding to the interaction of three H(ls) atoms. The results indicate that Gaussian functions can form good basis sets for obtaining long range forces for a variety of multipole interaction energies.     

Energy levels and wave functions of weakly-bound 4Hex 20NeyH (x+y=2) systems using Pekeris coordinates and a symmetry-adapted Lanczos approach

Energy levels and wave functions of floppy triatomic rare gas hydrides are calculated using a Pekeris coordinate system and the importance of various triangular configurations is assessed through the calculation of reduced distribution functions and relative weights. The calculations are performed using a symmetry-adapted Lanczos recursion within the discrete variable representation. For the 4He2H- anion, the present results are compared with those obtained from calculations based on other methods, and the accuracy of the present method is discussed. Calculations are also performed for the case of 4He2H and 20Ne2H, as well as for the mixed 4He 20NeH neutrals. Our results show that no bound states are found for 4He2H while only one bound state is found for both the 20Ne2H and 4He 20NeH complexes. Interestingly, a very important and common property of these systems is that there is a significant contribution from linear configurations to their bound states.     

Continuum wave functions by least-squares scheme in a B-spline basis: Multicenter and multielectron formulations

Application of a purely algebraic least-squares approach in a spline basis to the determination of continuum orbitals in the single-electron molecular case illustrates the excellent numerical properties of the method. Initial work based on a single-center expansion is extended to an LCAO formulation within a single-particle Hamiltonian and to the full multielectron problem in the one-center scheme. Preliminary results on H and He photoionization are reported. © 1994 John Wiley & Sons, Inc.