Search

Many-body perturbation theory for a direct calculation of the electrostatic interaction energy is developed. Since no multipole expansion is used, the obtained electrostatic energy includes the short-range contributions resulting from the overlap (penetration) of monomers’ charge distributions. The influence of intramonomer electronic correlation is accounted for by the perturbation expansion in terms of the Mo\ller–Plesset type fluctuation potentials for the interacting molecules. Two types of expansions are introduced: one based on the standard Mo\ller–Plesset expansion of the electron density, and the second accounting for the perturbation induced modifications of the monomer’s Fock operators, i.e., for the so-called response or orbital relaxation effects. Explicit orbital expressions for the terms through the fourth order in the intramonomer fluctuation potentials are derived. In this way the leading three-particle correlation contribution to the electrostatic energy is taken into account. Numerical results of perturbative calculations through the fourth order for the He2, (H2)2, and He–H2 systems are presented and compared with the complete configuration interaction results obtained using the same basis sets. It is found that the convergence of the many-body perturbation expansion of the electrostatic energy is fast. The sum of corrections through the fourth order reproduces to within few percent the correlation part of the electrostatic energy computed with the full configuration interaction wave function. [ABSTRACT FROM AUTHOR]

The pulsed-laser photolysis/laser-induced fluorescence method is used to study the kinetics of the reaction of NH2 with H2O2 to yield a second-order rate constant of (2.42 ± 0.55) × 10−14 cm3 molecule−1 s−1 at 412 K in 10–22 mbar of Ar bath gas. There are no prior measurements for comparison. To check this value and enable reliable extrapolation to other temperatures, we also compute thermal rate constants for this process over the temperature range 298–3000 K via multi-structural canonical variational transition-state theory with small-curvature multidimensional tunneling (MS-CVT/SCT). The CVT/SCT rate constants are derived using a dual-level direct dynamics approach utilizing single-point CCSD(T)-F12b/cc-pVQZ-F12 energies—corrected for core-valence and scalar relativistic effects—and M06-2X/MG3S geometries, gradients, and Hessians—for all stationary and non-stationary points along the reaction path. The multistructural method with torsional anharmonicity, based on a coupled torsional potential, is then employed to calculate correction factors for the rate constants, accounting for the comprehensive effects of torsional anharmonicity on the kinetics of this reaction system. The final MS-CVT/SCT rate constants are found to be in good agreement with our measurements and can be expressed in modified Arrhenius form as 2.13 × 10−15 (T/298 K)4.02 exp(−513 K/T) cm3 molecule−1 s−1 over the temperature range of 298–3000 K. [ABSTRACT FROM AUTHOR]

This work examines the viability of matrix completion methods as cost-effective alternatives to full nuclear Hessians for calculating quantum and variational effects in chemical reactions. The harmonic variety-based matrix completion (HVMC) algorithm, developed in a previous study [S. J. Quiton et al., J. Chem. Phys. 153, 054122 (2020)], exploits the low-rank character of the polynomial expansion of potential energy to recover vibrational frequencies (square roots of eigenvalues of nuclear Hessians) constituting the reaction path using a small sample of its entities. These frequencies are essential for calculating rate coefficients using variational transition state theory with multidimensional tunneling (VTST-MT). HVMC performance is examined for four SN2 reactions and five hydrogen transfer reactions, with each H-transfer reaction consisting of at least one vibrational mode strongly coupled to the reaction coordinate. HVMC is robust and captures zero-point energies, vibrational free energies, zero-curvature tunneling, and adiabatic ground state and free energy barriers as well as their positions on the reaction coordinate. For medium to large reactions involving H-transfer, with the sole exception of the most complex Ir catalysis system, less than 35% of total eigenvalue information is necessary for accurate recovery of key VTST-MT observables. [ABSTRACT FROM AUTHOR]

Symmetry-adapted perturbation theory (SAPT) has become an invaluable tool for studying the fundamental nature of non-covalent interactions by directly computing the electrostatics, exchange (steric) repulsion, induction (polarization), and London dispersion contributions to the interaction energy using quantum mechanics. Further application of SAPT is primarily limited by its computational expense, where even its most affordable variant (SAPT0) scales as the fifth power of system size [ O ( N 5 ) ] due to the dispersion terms. The algorithmic scaling of SAPT0 is reduced from O ( N 5 ) → O ( N 4 ) by replacing these terms with the empirical D3 dispersion correction of Grimme and co-workers, forming a method that may be termed SAPT0-D3. Here, we optimize the damping parameters for the -D3 terms in SAPT0-D3 using a much larger training set than has previously been considered, namely, 8299 interaction energies computed at the complete-basis-set limit of coupled cluster through perturbative triples [CCSD(T)/CBS]. Perhaps surprisingly, with only three fitted parameters, SAPT0-D3 improves on the accuracy of SAPT0, reducing mean absolute errors from 0.61 to 0.49 kcal mol−1 over the full set of complexes. Additionally, SAPT0-D3 exhibits a nearly 2.5× speedup over conventional SAPT0 for systems with ∼300 atoms and is applied here to systems with up to 459 atoms. Finally, we have also implemented a functional group partitioning of the approach (F-SAPT0-D3) and applied it to determine important contacts in the binding of salbutamol to G-protein coupled β1-adrenergic receptor in both active and inactive forms. SAPT0-D3 capabilities have been added to the open-source Psi4 software. [ABSTRACT FROM AUTHOR]

We present quasiclassical trajectory simulations of O2 + O collisions under conditions representative of a crossed molecular beam experiment [Lahankar et al., J. Phys. Chem. A 120, 5348–5359 (2016)]. These calculations are compared to experimental data in order to further validate Potential Energy Surfaces (PESs) recently developed at the University of Minnesota [Z. Varga, Y. Paukku, and D. G. Truhlar, J. Chem. Phys. 147, 154312 (2017)]. Spin and spatial symmetries allow nine adiabatic PESs to represent the ground state interactions of O2 + O. We simulate trajectories adiabatically on all nine surfaces and perform analysis on data from each surface separately and the combined set of nine. It is shown that aggregated adiabatic calculations from nine surfaces agree better with the experiment for both inelastic and exchange collisions than a previous study that used a single surface and lie within the reported experimental uncertainty at almost all points. Distributions for exchange reactions using all nine PESs are interpreted using each surface's opacity function and activation energy for exchange. Rovibrationally resolved product distributions are then studied for a better understanding of energy relaxation in exchange collisions and may prove useful for further quasiclassical trajectory study and new experiments that use rovibrational spectroscopy to resolve the internal energy of the scattered products. [ABSTRACT FROM AUTHOR]

Symmetry-adapted perturbation theory (SAPT) is a method for calculations of intermolecular (noncovalent) interaction energies. The set of SAPT codes that is described here, the current version named SAPT2020, includes virtually all variants of SAPT developed so far, among them two-body SAPT based on perturbative, coupled cluster, and density functional theory descriptions of monomers, three-body SAPT, and two-body SAPT for some classes of open-shell monomers. The properties of systems governed by noncovalent interactions can be predicted only if potential energy surfaces (force fields) are available. SAPT is the preferred approach for generating such surfaces since it is seamlessly connected to the asymptotic expansion of interaction energy. SAPT2020 includes codes for automatic development of such surfaces, enabling generation of complete dimer surfaces with a rigid monomer approximation for dimers containing about one hundred atoms. These codes can also be used to obtain surfaces including internal degrees of freedom of monomers. [ABSTRACT FROM AUTHOR]

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook. [ABSTRACT FROM AUTHOR]

A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods such as the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree–Fock, and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as it is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite such as GAMESS is maintained and developed is considered. Future developments are briefly summarized. [ABSTRACT FROM AUTHOR]

An analytic full-dimensional diabatic potential energy matrix (DPEM) for the lowest three singlet states of thiophenol (C6H5SH) at geometries accessible during photodissociation is constructed using the anchor points reactive potential (APRP) scheme. The data set used for modeling is obtained from electronic structure calculations including dynamic correlation via excitations into the virtual space of a three-state multiconfiguration self-consistent field calculation. The resulting DPEM is a function of all the internal coordinates of thiophenol. The DPEM as a function of the S−H bond stretch and C−C−S−H torsion and the diabatic couplings along two in-plane bend modes and nine out-of-plane distortion modes are computed using extended multiconfigurational quasidegenerate perturbation theory followed by the fourfold way determination of diabatic molecular orbitals and model space diabatization by configurational uniformity, and this dependence of the DPEM is represented by general functional forms. Potentials along 31 tertiary internal degrees of freedom are modeled with system-dependent, primary-coordinate-dependent nonreactive molecular mechanics-type force fields that are parameterized by Cartesian Hessians calculated by generalized Kohn-Sham density functional theory. Adiabatic potential energy surfaces (PESs) and nonadiabatic couplings are obtained by a transformation of the DPEM. The topography of the APRP PESs is characterized by vertical excitation energies, equilibrium geometries, vibrational frequencies, and conical intersections, and we find good agreement with available reference data. This analytic DPEM is suitable for full-dimensional electronically nonadiabatic molecular dynamics calculations of the photodissociation of thiophenol with analytic gradients in either the adiabatic or diabatic representation. [ABSTRACT FROM AUTHOR]

The gas-phase kinetics of S(³P) atoms with H2 and D2 have been studied via the laser flash photolysis—resonance fluorescence technique. S atoms were generated by pulsed photolysis of CS2 at 193 nm and monitored by time-resolved fluorescence at 181 nm. The rate coefficients for H2 (k1) and D2 (k2), respectively, are summarized as k1(600-1110 K) = 3.0 × 109 exp (-1.317×105-2.703×107K=T/8.314 T=K) cm³ molecule-1 s-1 and k2(770-1110 K) = 2.2 × 10-14 (T/298 K)3.55 exp(-5420 K/T) cm³ molecule-1 s-1. Error limits are discussed in the text. The rate coefficients for formation of SH(SD) + H(D) on a newly developed triplet potential energy surface were characterized via ring polymer molecular dynamics and canonical variational transition-state theory. There is excellent agreement above about 1000 K between theory and experiment. At lower temperatures, the experimental rate coefficient is substantially larger than the results computed for the adiabatic reaction, suggesting a significant role for intersystem crossing to the singlet potential energy surface at lower temperatures. [ABSTRACT FROM AUTHOR]

We measured the rates of abstraction of a hydrogen atom from specific sites in propane C3H8, 2-methyl propane (i-C4H10), and butane (n-C4H10); the sites are a primary hydrogen of C3H8 and i-C4H10 and a secondary hydrogen of n-C4H10. The excellent reproducibility of conditions of a diaphragm-less shock tube enabled us to conduct comparative measurements of the evolution of H atoms in three mixtures--(i) 0.5 ppm C2H5I + Ar, (ii) 0.5 ppm C2H5I + 50-100 ppm alkane as C3H8 or i-C4H10 or n-C4H10 + Ar, and (iii) the same concentrations of alkane + Ar without C2H5I--in the temperature range 1000-1200 K and at a pressure of 2.0 bars. The net profile of rise and decay of H atoms in the C2H5I + alkane mixture was derived on subtracting the absorbance of (iii) from that of (ii). Measurements of the mixture (iii) are important because the absorption of alkanes at 121.6 nm is not negligible. In the temperature range 1000-1100 K, the rate of decomposition of C2H5I was evaluated directly on analyzing the exponential growth of H atoms in the mixture (i). The rate of decomposition of C2H5I is summarized as ln(k/s-1) = (33.12 ± 1.4) - (25.23 ± 1.5) 10³/T (T = 1000-1100 K, P = 2.0 bars); the broadening factor F(T) in the Lindemann-Hinshelwood formula was evaluated in the fall-off region. The site-specific rates of H + (C3-C4) alkanes are summarized as follows: H + C3H8 → H2 + 1-C3H7, ln(k1a) = -(21.34 ± 0.86) - (5.39 ± 0.93)103/T, H + i-C4H10 → H2 + i-C4H9, ln(k2a) = -(20.50 ± 1.36) - (6.14 ± 0.13)10³/T, H + n-C4H10 → H2 + 2-C4H9, ln(k3b) = -(21.37 ± 1.15) - (4.83 ± 1.26)10³/T. The present experimental results are compared with published results from quantum-chemical calculations of potential-energy surfaces and transition-state theory. The present experiments are consistent with those calculations for the reaction rates for the attack at the primary site for H + C3H8 and H + i-C4H10, but for the attack at the secondary site of n-C4H10, our results are substantially smaller than the computational prediction, which might indicate a hindrance by the C-H bonds of the primary sites that serves to decrease the rate of abstraction from the secondary site of n-C4H10. The influence on the total rates of reactions H + alkane and the group additivity rule are discussed. [ABSTRACT FROM AUTHOR]

An ab initio investigation has been employed to determine the structural and spectroscopic parameters, such as rotational constants, vibrational frequencies, vertical excitation energies, and the stability of the ethyl-water complex. The ethyl-water complex has a binding energy of 1.15 kcal·mol−1. The interaction takes place between the hydrogen of water and the unpaired electron of the radical. This interaction is found to produce a red shift in the OH stretching bands of water of ca. 84 cm−1, and a shift of all UV absorption bands to higher energies. [ABSTRACT FROM AUTHOR]

The adsorption of noble gases on metallic surfaces represents a paradigmatic case of van-der-Waals (vdW) interaction due to the role of screening effects on the corrugation of the interaction potential [J. L. F. Da Silva et al., Phys. Rev. Lett. 90, 066104 (2003)]. The extremely small adsorption energy of He atoms on the Mg(0001) surface (below 3 meV) and the delocalized nature and mobility of the surface electrons make the He/Mg(0001) system particularly challenging, even for state-of-the-art vdW-corrected density functional-based (vdW-DFT) approaches [M. P. de Lara-Castells et al., J. Chem. Phys. 143, 194701 (2015)]. In this work, we meet this challenge by applying two different procedures. First, the dispersion-corrected second-order Möller-Plesset perturbation theory (MP2C) approach is adopted, using bare metal clusters of increasing size. Second, the method of increments [H. Stoll, J. Chem. Phys. 97, 8449 (1992)] is applied at coupled cluster singles and doubles and perturbative triples level, using embedded cluster models of the metal surface. Both approaches provide clear evidences of the anti-corrugation of the interaction potential: the He atom prefers on-top sites, instead of the expected hollow sites. This is interpreted as a signature of the screening of the He atom by the metal for the on-top configuration. The strong screening in the metal is clearly reflected in the relative contribution of successively deeper surface layers to the main dispersion contribution. Aimed to assist future dynamical simulations, a pairwise potential model for the He/surface interaction as a sum of effective He-Mg pair potentials is also presented, as an improvement of the approximation using isolated He-Mg pairs. [ABSTRACT FROM AUTHOR]

The accuracy and transferability of the electronic structure approach combining dispersionless density functional theory (DFT) [K. Pernal et al., Phys. Rev. Lett. 103, 263201 (2009)] with the method of increments [H. Stoll, J. Chem. Phys. 97, 8449 (1992)], are validated for the interaction between the noble-gas Ne, Ar, Kr, and Xe atoms and coronene/graphene/graphite surfaces. This approach uses the method of increments for surface cluster models to extract intermonomer dispersion-like (2- and 3-body) correlation terms at coupled cluster singles and doubles and perturbative triples level, while periodic dispersionless density functionals calculations are performed to estimate the sum of Hartree-Fock and intramonomer correlation contributions. Dispersion energy contributions are also obtained using DFT-based symmetry-adapted perturbation theory [SAPT(DFT)]. An analysis of the structure of the X/surface (X = Ne, Ar, Kr, and Xe) interaction energies shows the excellent transferability properties of the leading intermonomer correlation contributions across the sequence of noble-gas atoms, which are also discussed using the Drude oscillator model. We further compare these results with van der Waals-(vdW)-corrected DFT-based approaches. As a test of accuracy, the energies of the low-lying nuclear bound states supported by the laterally averaged X/graphite potentials (X = 3He, 4He, Ne, Ar, Kr, and Xe) are calculated and compared with the best estimations from experimental measurements and an atom-bond potential model using the ab initio-assisted finetuning of semiempirical parameters. The bound-state energies determined differ by less than 6-7 meV (6%) from the atom-bond potential model. The crucial importance of including incremental 3-body dispersion-type terms is clearly demonstrated, showing that the SAPT(DFT) approach effectively account for these terms. With the deviations from the best experimental-based estimations smaller than 2.3 meV (1.9%), the accuracy of the combined DFT and post-HF incremental scheme is established for all the noble-gas atoms. With relative deviations smaller than 4% and 11%, good agreement is also achieved by applying the vdW-corrected DFT treatments PBE-D3 and vdW-DF2 for noble-gas atoms heavier than neon [ABSTRACT FROM AUTHOR]

Ground state potential energy curves for homonuclear and heteronuclear dimers consisting of noble gas atoms from He to Kr were calculated within the symmetry adapted perturbation theory based on the density functional theory (DFT-SAPT). These potentials together with spectroscopic data derived from them were compared to previous high-precision coupled cluster with singles and doubles including the connected triples theory calculations (or better if available) as well as to experimental data used as the benchmark. The impact of midbond functions on DFT-SAPT results was tested to study the convergence of the interaction energies. It was shown that, for most of the complexes, DFT-SAPT potential calculated at the complete basis set (CBS) limit is lower than the corresponding benchmark potential in the region near its minimum and hence, spectroscopic accuracy cannot be achieved. The influence of the residual term δ(HF) on the interaction energy was also studied. As a result, we have found that this term improves the agreement with the benchmark in the repulsive region for the dimers considered, but leads to even larger overestimation of potential depth De. Although the standard hybrid exchange-correlation (xc) functionals with asymptotic correction within the second order DFT-SAPT do not provide the spectroscopic accuracy at the CBS limit, it is possible to adjust empirically basis sets yielding highly accurate results. [ABSTRACT FROM AUTHOR]

Theory of one-electron transition density matrices has been formulated within the time-independent coupled cluster method for the polarization propagator [R. Moszynski, P. S. Zuchowski, and B. Jeziorski, Coll. Czech. Chem. Commun. 70, 1109 (2005)]. Working expressions have been obtained and implemented with the coupled cluster method limited to single, double, and linear triple excitations (CC3). Selected dipole and quadrupole transition probabilities of the alkali earth atoms, computed with the new transition density matrices are compared to the experimental data. Good agreement between theory and experiment is found. The results obtained with the new approach are of the same quality as the results obtained with the linear response coupled cluster theory. The one-electron density matrices for the ground state in the CC3 approximation have also been implemented. The dipole moments for a few representative diatomic molecules have been computed with several variants of the new approach, and the results are discussed to choose the approximation with the best balance between the accuracy and computational efficiency. [ABSTRACT FROM AUTHOR]

Binding energies for the complexes of the S12L database by Grimme [Chem. Eur. J. 18, 9955 (2012)] were calculated using intermolecular symmetry-adapted perturbation theory combined with a density-functional theory description of the interacting molecules. The individual interaction energy decompositions revealed no particular change in the stabilisation pattern as compared to smaller dimer systems at equilibrium structures. This demonstrates that, to some extent, the qualitative description of the interaction of small dimer systems may be extrapolated to larger systems, a method that is widely used in force-fields in which the total interaction energy is decomposed into atomatom contributions. A comparison of the binding energies with accurate experimental reference values from Grimme, the latter including thermodynamic corrections from semiempirical calculations, has shown a fairly good agreement to within the error range of the reference binding energies. [ABSTRACT FROM AUTHOR]

A systematic examination of the computational expense and accuracy of Symmetry-Adapted Perturbation Theory (SAPT) for the prediction of non-covalent interaction energies is provided with respect to both method [SAPT0, DFT-SAPT, SAPT2, SAPT2+, SAPT2+(3), and SAPT2+3; with and without CCD dispersion for the last three] and basis set [Dunning cc-pVDZ through aug-cc-pV5Z wherever computationally tractable, including truncations of diffuse basis functions]. To improve accuracy for hydrogen-bonded systems, we also include two corrections based on exchange-scaling (sSAPT0) and the supermolecular MP2 interaction energy (δMP2). When considering the best error performance relative to computational effort, we recommend as the gold, silver, and bronze standard of SAPT: SAPT2+(3)δMP2/aug-cc-pVTZ, SAPT2+/aug-cc-pVDZ, and sSAPT0/jun-ccpVDZ. Their respective mean absolute errors in interaction energy across the S22, HBC6, NBC10, and HSG databases are 0.15 (62.9), 0.30 (4.4), and 0.49 kcal mol-1 (0.03 h for adenine · thymine complex). [ABSTRACT FROM AUTHOR]

The performance of second-order symmetry-adapted perturbation theory (SAPT) calculations using Kohn-Sham (KS) orbitals is evaluated against benchmark results for intermolecular interactions. Unlike previous studies of this "SAPT(KS)" methodology, the present study uses non-empirically tuned long-range corrected (LRC) functionals for the monomers. The proper vxc(r) → 0 asymptotic limit is achieved by tuning the range separation parameter in order to satisfy the condition that the highest occupied KS energy level equals minus the molecule's ionization energy, for each monomer unit. Tests for He2, Ne2, and the S22 and S66 data sets reveal that this condition is important for accurate prediction of the non-dispersion components of the energy, although errors in SAPT(KS) dispersion energies remain unacceptably large. In conjunction with an empirical dispersion potential, however, the SAPT(KS) method affords good results for S22 and S66, and also accurately predicts the whole potential energy curve for the sandwich isomer of the benzene dimer. Tuned LRC functionals represent an attractive alternative to other asymptotic corrections that have been employed in density-functional-based SAPT calculations, and we recommend the use of tuned LRC functionals in both coupled-perturbed SAPT(DFT) calculations and dispersion-corrected SAPT(KS) calculations. [ABSTRACT FROM AUTHOR]

We report the definition and refinement of a new first principles potential for the acetylene dimer. The ab initio calculations were performed with the DFT-SAPT combination of symmetry-adapted intermolecular perturbation method and density functional theory, and fitted to a model site-site functional form. Comparison of the calculated microwave spectrum with experimental data revealed that the barriers to isomerization were too low. This potential was refined by fitting the model parameters in order to reproduce the observed transitions, an excellent agreement within ∼1 MHz being achieved. [ABSTRACT FROM AUTHOR]

Symmetry-adapted perturbation theory (SAPT) provides a means of probing the fundamental nature of intermolecular interactions. Low-orders of SAPT (here, SAPT0) are especially attractive since they provide qualitative (sometimes quantitative) results while remaining tractable for large systems. The application of density fitting and Laplace transformation techniques to SAPT0 can significantly reduce the expense associated with these computations and make even larger systems accessible. We present new factorizations of the SAPT0 equations with density-fitted two-electron integrals and the first application of Laplace transformations of energy denominators to SAPT. The improved scalability of the DF-SAPT0 implementation allows it to be applied to systems with more than 200 atoms and 2800 basis functions. The Laplace-transformed energy denominators are compared to analogous partial Cholesky decompositions of the energy denominator tensor. Application of our new DF-SAPT0 program to the intercalation of DNA by proflavine has allowed us to determine the nature of the proflavine-DNA interaction. Overall, the proflavine-DNA interaction contains important contributions from both electrostatics and dispersion. The energetics of the intercalator interaction are are dominated by the stacking interactions (two-thirds of the total), but contain important contributions from the intercalator-backbone interactions. It is hypothesized that the geometry of the complex will be determined by the interactions of the intercalator with the backbone, because by shifting toward one side of the backbone, the intercalator can form two long hydrogen-bonding type interactions. The long-range interactions between the intercalator and the next-nearest base pairs appear to be negligible, justifying the use of truncated DNA models in computational studies of intercalation interaction energies. [ABSTRACT FROM AUTHOR]

We report a fragment-based electronic structure method, intended for the study of clusters and molecular liquids, that incorporates electronic polarization (induction) in a self-consistent fashion but treats intermolecular exchange and dispersion interactions perturbatively, as post-self-consistent field corrections, using a form of pairwise symmetry-adapted perturbation theory. The computational cost of the method scales quadratically as a function of the number of fragments (monomers), but could be made to scale linearly by exploiting distance-dependent thresholds. Extensive benchmark calculations are reported using the S22 database of high-level ab initio binding energies for dimers, and we find that average errors can be reduced to <1 kcal/mol with a suitable choice of basis set. Comparison to ab initio benchmarks for water clusters as large as (H2O)20 demonstrates that the method recovers >=90% of the binding energy in these systems, at a tiny fraction of the computational cost. As such, this approach represents a promising path toward accurate, systematically improvable, and parameter-free simulation of molecular liquids. [ABSTRACT FROM AUTHOR]

An accurate description of dispersion interactions is required for reliable theoretical studies of many noncovalent complexes. This can be obtained with the wave function-based formulation of symmetry-adapted perturbation theory (SAPT) provided that the contribution of triple excitations to dispersion is included. Unfortunately, this triples dispersion correction limits the applicability of SAPT due to its O(N7) scaling. The efficiency of the evaluation of this correction can be greatly improved by removing virtual orbitals from the computation. The error incurred from truncating the virtual space is reduced if second-order Mo\ller-Plesset perturbation theory (MP2) natural orbitals are used in place of the canonical Hartree-Fock molecular orbitals that are typically used. This approximation is further improved if the triples correction to dispersion is scaled to account for the smaller virtual space. If virtual MP2 natural orbitals are removed according to their occupation numbers, in practice, roughly half of the virtual orbitals can be removed (with the aug-cc-pVDZ basis set) with negligible errors if the remaining triples dispersion contribution is scaled. This typically leads to speedups of 15-20 times for the cases considered here. By combining the truncated virtual space with the frozen core approximation, the triples correction can be evaluated approximately 50 times faster than the canonical computation. These approximations cause less than 1% error (or at most 0.02 kcal mol-1) for the cases considered. Truncation of greater fractions of the virtual space is possible for larger basis sets (leading to speedups of over 40 times before additional speedups from the frozen core approximation). [ABSTRACT FROM AUTHOR]

Symmetry-adapted perturbation theory (SAPT) offers insight into the nature of intermolecular interactions. In addition, accurate energies can be obtained from the wave function-based variant of SAPT provided that intramonomer electron correlation effects are included. We apply density-fitting (DF) approximations to the intramonomer correlation corrections in SAPT. The introduction of this approximation leads to an improvement in the computational cost of SAPT by reducing the scaling of certain SAPT terms, reducing the amount of disk I/O, and avoiding the explicit computation of certain types of MO integrals. We have implemented all the intramonomer correlation corrections to SAPT through second-order under the DF approximation. Additionally, leading third-order terms are also implemented. The accuracy of this truncation of SAPT is tested against the S22 test set of Hobza and co-workers [Phys. Chem. Chem. Phys. 8, 1985 (2006)]. When the intramonomer corrections to dispersion are included in SAPT, a mean absolute deviation of 0.3–0.4 kcal mol-1 is observed for the S22 test set when using an aug-cc-pVDZ basis. The computations on the adenine-thymine complexes in the S22 test set with an aug-cc-pVDZ basis represent the largest SAPT computations to date that include this degree of intramonomer correlation. Computations of this size can now be performed routinely with our newly developed DF-SAPT program. [ABSTRACT FROM AUTHOR]

Density fitting (DF) approximations have been used to increase the efficiency of several quantum mechanical methods. In this work, we apply DF and a related approach, Cholesky decomposition (CD), to wave function-based symmetry-adapted perturbation theory (SAPT). We also test the one-center approximation to the Cholesky decomposition. The DF and CD approximations lead to a dramatic improvement in the overall computational cost of SAPT, while introducing negligible errors. For typical target accuracies, the Cholesky basis needed is noticeably smaller than the DF basis (although the cost of constructing the Cholesky vectors is slightly greater than that of constructing the three-index DF integrals). The SAPT program developed in this work is applied to the interactions between acenes previously studied by Grimme [Angew. Chem., Int. Ed. 47, 3430 (2008)], expanding the cases studied by adding the pentacene dimer. The SAPT decomposition of the acene interactions provides a more realistic picture of the interactions than that from the energy decomposition analysis previously reported. The data suggest that parallel-displaced and T-shaped acene dimers both feature a special stabilizing π-π interaction arising from electron correlation terms which are significantly more stabilizing than expected on the basis of pairwise -C6R-6 estimates. These terms are qualitatively the same in T-shaped as in parallel-displaced geometries, although they are roughly a factor of 2 smaller in T-shaped geometries because of the larger distances between the intermolecular pairs of electrons. [ABSTRACT FROM AUTHOR]

A distributed point polarizable model (DPP2) for water, with explicit terms for charge penetration, induction, and charge transfer, is introduced. The DPP2 model accurately describes the interaction energies in small and large water clusters and also gives an average internal energy per molecule and radial distribution functions of liquid water in good agreement with experiment. A key to the success of the model is its accurate description of the individual terms in the n-body expansion of the interaction energies. [ABSTRACT FROM AUTHOR]

A density matrix-based Laplace reformulation of coupled-perturbed self-consistent field (CPSCF) theory is presented. It allows a direct, instead of iterative, solution for the integral-independent part of the density matrix-based CPSCF (D-CPSCF) equations [J. Kussmann and C. Ochsenfeld, J. Chem. Phys. 127, 054103 (2007)]. In this way, the matrix-multiplication overhead compared to molecular orbital-based solutions is reduced to a minimum, while at the same time, the linear-scaling behavior of D-CPSCF theory is preserved. The present Laplace-based equation solver is expected to be of general applicability. [ABSTRACT FROM AUTHOR]

The nature of interactions of furan with various hydrides (Y) (Y=HF,HCl,H2O,H2S,NH3,PH3) is investigated using ab initio calculations. The contribution of attractive (electrostatic, inductive, and dispersive) and repulsive (exchange) components to the interactions energy is analyzed. HF, H2O, and NH3 favor σo-type H bonding, while HCl, H2S, and PH3 favor π-type H bonding. Interaction energy decomposition reveals that σo-type complexes interactions are predominantly electrostatic in nature, while the dispersion and electrostatic interactions dominate the π-type complexes. [ABSTRACT FROM AUTHOR]

The application of the frozen-core approximation (FCA) and effective core potentials (ECPs) within symmetry-adapted perturbation theory (SAPT) has been investigated and implemented. Unlike in the case of conventional electronic-structure theories, the development of a frozen-core version of SAPT is not straightforward. In particular, the FCA realizations neglecting excitations from core orbitals and restricting all summation indices to valence orbitals only are no longer equivalent. It is shown that it is necessary in SAPT to keep some terms containing products of the valence orbitals of one monomer and the core orbitals of the other one in the exchange-energy components. When these terms are included or, equivalently, the “infinite-excitation-energy” approximation omitting only the excitations from the core orbitals is used, the accuracy of the frozen-core approximation in SAPT matches that obtained in supermolecular perturbational and coupled-cluster methods. If these terms are neglected, i.e., within the “index-range-restriction” approximation, several exchange corrections are significantly underestimated. When ECPs are used in SAPT, the accuracy of the interaction energies is as good as in conventional supermolecular methods, provided that the residual supermolecular Hartree-Fock term is included. We have found that only some types of ECPs can be reliably used for calculations of interaction energies both in SAPT and in supermolecular approaches. For systems containing heavy atoms, both FCA and the use of ECPs lead to very significant savings of computer time. [ABSTRACT FROM AUTHOR]

On the basis of density functional theory (DFT) and high level ab initio theory, we report the structures, binding energies, thermodynamic quantities, IR spectra, and electronic properties of the hydride anion hydrated by up to six water molecules. Ground state DFT molecular dynamics simulations (based on the Born-Oppenheimer potential surface) show that as the temperature increases, the surface-bound hydride anion changes to the internally bound structure. Car-Parrinello molecular dynamics simulations are also carried out for the spectral analysis of the monohydrated hydride. Excited-state ab initio molecular dynamics simulations show that the photoinduced charge-transfer-to-solvent phenomena are accompanied by the formation of the excess electron-water clusters and the detachment of the H radical from the clusters. The dynamics of the detachment process of a hydrogen radical upon the excitation is discussed. [ABSTRACT FROM AUTHOR]

In this work the authors present an experimental and theoretical study about the Q-branch lines’ broadening coefficients of N2 perturbed by H2. Experimental values for these parameters have been obtained at 440 and 580 K, and quantum calculations have been performed using a new ab initio potential energy surface, obtained by quantum chemistry methods. The results of these calculations are compared to experimental data obtained previously at 77 and 298 K [L. Gomez et al., Mol. Phys. 104, 1869 (2006)] and to the present measurements. A satisfactory agreement is obtained for the whole range of temperatures used in the experiments. [ABSTRACT FROM AUTHOR]

Equilibrium geometries, interaction energies, atomic charge, and charge transfer for the intermolecular interactions between furan and dihalogen molecules XY(X; Y=F,Cl,Br) were studied at the MP2/aug-cc-pVDZ level. Three types of geometry are observed in these interactions: the π-type geometry (I), in which the XY lies above the furan ring and almost perpendicularly to the C4–C5 bond of furan; the σ-type geometry (II), where the X atom is pointed toward the nonbonding electron pair (n pair) of oxygen atom in furan; and the χ-type geometry (III), describing a blueshift hydrogen bond formed between the hydrogen atom of furan and dihalogen molecules XY. The calculated interaction energies show that the π-type structures are more stable than the corresponding σ-type and χ-type structures. To study the nature of the intermolecular interactions, an energy decomposition analysis was carried out and the results indicate that both the π-type and σ-type interactions are dominantly inductive energy in nature, while dispersion energy governs the χ-type interactions. [ABSTRACT FROM AUTHOR]

One-electron density matrices resulting from the explicitly connected commutator expansion of the expectation value were implemented at the singles and doubles coupled cluster (CCSD) level. In the proposed approach the one-electron density matrix is obtained at a little extra cost in comparison to the calculation of the CCSD correlation energy. Therefore, in terms of the computational time the new method is significantly less demanding than the conventional linear-response CCSD theory which requires additionally an expensive calculation of the left-hand solution of the CCSD equations. The quality of the new density matrices was investigated by computing a set of one-electron properties for a series of molecules of varying sizes and comparing the results with data obtained using the full configuration interaction method or higher level coupled cluster theory. It has been found that the results obtained using the new approach are of the same quality as those predicted by the linear-response CCSD method. The novel one-electron density matrices have also been applied to study the energy of the electrostatic interaction for a number of van der Waals complexes, including the benzene and azulene dimers. [ABSTRACT FROM AUTHOR]

The interacting quantum atoms approach [IQA, as presented by Blanco et al., J. Chem. Theory Comput. 1, 1096 (2005)] is applied to standard hydrogen bonded dimers. IQA is an interpretation tool based on a real space energy decomposition scheme fully consistent with the quantum theory of atoms in molecules. It provides a partition of every physical term present in the Hamiltonian into atomic and interatomic contributions. The procedure is orbital-free and self-contained, needing neither external references nor artificial intermediate states. Binding is the result of a competition between the destabilizing deformations suffered by the interacting fragments upon interaction and the stabilizing interaction energy itself. According to IQA, there is no incompatibility between the prevalent electrostatic image of hydrogen bonded systems and that favoring important covalent contributions. Depending on how we gather the different energetic terms, we may recover electrostatic or covalent pictures from the same underlying quantum mechanical description. Our results show that the nonclassical contributions to hydrogen bonding are spatially localized, involving only the H atom and its two nearest neighbors. IQA is well suited as a comparative tool. Its thin energetic decomposition allows us to recover exactly (or to a very good approximation) the quantities of the most widely used energy decomposition schemes. Such a comparison sheds light on the virtues and faults of the different methods and on the origin of the 50 years old debate regarding the covalent/electrostatic nature of the hydrogen bond. [ABSTRACT FROM AUTHOR]

We present an extension of many-body symmetry-adapted perturbation theory (SAPT) by including all third-order polarization and exchange contributions obtained with the neglect of intramonomer correlation effects. The third-order polarization energy, which naturally decomposes into the induction, dispersion, and mixed, induction-dispersion components, is significantly quenched at short range by electron exchange effects. We propose a decomposition of the total third-order exchange energy into the exchange-induction, exchange-dispersion, and exchange-induction-dispersion contributions which provide the quenching for the corresponding individual polarization contributions. All components of the third-order energy have been expressed in terms of molecular integrals and orbital energies. The obtained formulas, valid for both dimer- and monomer-centered basis sets, have been implemented within the general closed-shell many-electron SAPT program. Test calculations for several small dimers have been performed and their results are presented. For dispersion-bound dimers, the inclusion of the third-order effects eliminates the need for a hybrid SAPT approach, involving supermolecular Hartree-Fock calculations. For dimers consisting of strongly polar monomers, the hybrid approach remains more accurate. It is shown that, due to the extent of the quenching, the third-order polarization effects should be included only together with their exchange counterparts. Furthermore, the latter have to be calculated exactly, rather than estimated by scaling the second-order values. [ABSTRACT FROM AUTHOR]

Empirically correlated density matrices of N-electron systems are investigated. Closed-form expressions are derived for the one- and two-electron reduced density matrices from a pairwise correlated wave function. Approximate expressions are then proposed which reflect dispersive interactions between closed-shell centrosymmetric subsystems. Said expressions clearly illustrate the consequences of second-order correlation effects on the reduced density matrices. Application is made to a simple example: the He2 system. Reduced density matrices are explicitly calculated, correct to second order in correlation, and compared with approximations of independent electrons and independent electron pairs. The models proposed allow for variational calculations of interaction energies and equilibrium distance as well as a clear interpretation of dispersive effects on electron distributions. Both exchange and second order correlation effects are shown to play a critical role on the quality of the results. [ABSTRACT FROM AUTHOR]

A symmetry-adapted perturbation theory based on Kohn–Sham determinants [SAPT(KS)] and utilizing asymptotically corrected exchange-correlation potentials has been applied to the He2, Ne2, (H2O)2, and (CO2)2 dimers. It is shown that SAPT(KS) is able to recover the electrostatic, first-order exchange, second-order induction, and exchange-induction energies with an accuracy approaching and occasionally surpassing that of regular SAPT at the currently programmed theory level. The use of the asymptotic corrections is critical to achieve this accuracy. The SAPT(KS) results can be obtained at a small fraction of the time needed for regular SAPT calculations. The robustness of the SAPT(KS) method with respect to the basis set size is also demonstrated. A theoretical justification for high accuracy of SAPT(KS) predictions for the electrostatic, first-order exchange, and second-order induction energies has been provided. [ABSTRACT FROM AUTHOR]

The ΔK=0 and 1 subbands of the (4,0)←(0,0) transition of (HF)2, near 14 700 cm-1, have been measured by molecular-beam intracavity laser-induced fluorescence. The hydrogen interchange tunneling is basically quenched in (4, 0) for both K=0 and 1 levels, consistent with the early suggestion from a phenomenological model [H.-C. Chang and W. Klemperer, J. Chem. Phys. 104, 7830 (1996)]. The band origin ν0=14 700.458(7) cm-1 and rotational constant (B+C)/2=0.222 78(31) cm-1 are determined for K=0 of the (4, 0) mode. From the observed ΔK=1←0 spectrum, we determined that A=24.3 cm-1, (B+C)/2=0.222 96(20) cm-1, and (B-C)=4.5(2)×10-3 cm-1. The predissociation linewidths of both K=0 and 1 levels are 470(30) MHz with no apparent rotational dependence. [ABSTRACT FROM AUTHOR]

A four-dimensional potential energy surface (PES) for the CO dimer consisting of rigid molecules has been calculated, using a scheme that combines density functional theory to describe the monomers and symmetry adapted perturbation theory for the interaction energy (DFT-SAPT). The potential is fitted in terms of analytic functions, and the fitted potential is used to compute the lowest rovibrational states of the dimer. The quality of the PES is comparable to that of a previously published surface [G. W. M. Vissers, P. E. S. Wormer, and A. van der Avoird, Phys. Chem. Chem. Phys., 5, 4767 (2003)], which was calculated with the coupled cluster single double and perturbative triples [CCSD(T)] method. It is shown that a weighted average of the DFT-SAPT and the CCSD(T) potential gives results that are in very good agreement with experimental data, for both (12CO)2 and (13CO)2. The relative weight was determined by adjusting the energy gap between the origins of the lowest two stacks of rotational levels of (12CO)2 to the measured value.© 2005 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The previously developed DFT-SAPT approach, which combines symmetry-adapted intermolecular perturbation theory (SAPT) with a density-functional theory (DFT) representation of the monomers, has been implemented by using density fitting of two-electron objects. This approach, termed DF-DFT-SAPT, scales with the fifth power of the molecular size and with the third power upon increase of the basis set size for a given dimer, thus drastically reducing the cost of the conventional DFT-SAPT method. The accuracy of the density fitting approximation has been tested for the ethyne dimer. It has been found that the errors in the interaction energies due to density fitting are below 10-3 kcal/mol with suitable auxiliary basis sets and thus one or two orders of magnitude smaller than the errors due to the use of a limited atomic orbital basis set. An investigation of three prominent structures of the benzene dimer, namely, the T shaped, parallel displaced, and sandwich geometries, employing basis sets of up to augmented quadruple-ζ quality shows that DF-DFT-SAPT outperforms second-order Møller-Plesset theory (MP2) and gives total interaction energies which are close to the best estimates infered from combining the results of MP2 and coupled-cluster theory with single, double, and perturbative triple excitations.© 2005 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The ab initio three-dimensional potential energy surface (PES) for the weak interaction of hydrogen molecule with bromine anion is presented. The surface was obtained by the supermolecular method at the coupled cluster with single and double excitations and noniterative correction to triple excitations (CCSD(T)) level of theory. Our calculations indicate the van der Waals (vdW) system for the linear orientation at R=3.37 Å with a well depth of De=660.1 cm-1. The presented PES reveals also transition state for the perpendicular orientation at R=4.22 Å with a barrier of 607.1 cm-1. The physical origin of the stability of vdW H2...Br- structure with respect to the H2...Cl- one was analyzed by the symmetry adapted perturbation theory based on the single determinant Hartree–Fock (HF) wave function. The separation of the interaction energy shows that the dispersion forces play slightly more important role in the stabilization of the vdW system with Br- than with Cl-. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Five two-dimensional potential energy surfaces for the interaction of He with cyanoacetylene (HCCCN) are presented, obtained from ab initio calculations using symmetry-adapted perturbation theory and the supermolecular method at different levels of electron correlation. HCCCN is taken to be a rigid linear molecule with the interatomic distances fixed at the experimental “r[sub 0]” geometry extracted from ground-state rotational constants. The complex was found to have a global minimum at a T-shaped configuration and a secondary minimum at the linear configuration with the He atom facing the H atom. Two saddle points were also located. There is good agreement between the positions of the stationary points on each of the five surfaces though their energies differ by up to 19%. Rovibrational bound state calculations were performed for the [sup 4]He–HCCCN and [sup 3]He–HCCCN complexes. Spectra (including intensities) and wave functions of [sup 4]He–HCCCN obtained from these calculations are presented. The effective rotational constant of HCCCN solvated in a helium droplet was estimated by minimizing the energy of He[sub n]–HCCCN for n=2–12, selecting the n=7 complex as giving the largest magnitude of interaction energy per He, and shifting the resulting ring of He atoms to the position corresponding to the average geometry of the ground state of the He–HCCCN dimer. This estimate is within 4.8% of the measured value. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The interface between water and Langmuir films of long chain aliphatic molecules is investigated using accurate intermolecular potentials. The stabilities of various ice structures which could form at the interface are examined. Antiferroelectric ice is found to be the most stable, but this stability depends crucially on the first layer of water. Ferroelectric structures are found to collapse upon relaxation. Our model was not able to differentiate between the different nucleation properties of C[sub 31]H[sub 63]OH and C[sub 30]H[sub 61]OH. A better description of the alcohol–water interaction is probably required to account for this difference. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The energy dependence of the capture cross section and the temperature dependence of the capture rate constants for inverse power attractive potentials V∝ - R[sup -n] is considered in the regime where the quantum character of the relative motion of colliding partners is important. For practically interesting cases n = 4 and n = 6, a simple formula for the cross section is suggested which interpolates between the classical and the quantum Bethe limits. We have shown that the classical approximation for the capture cross section performs well far below the simple estimations of the onset the quantum regime. This seemingly "classical" feature of the cross section and the rate constant is due to the large quantum effects of the waves in transmission through and reflection above the centrifugal potential barriers. [ABSTRACT FROM AUTHOR]

In this work we test two ab initio methodologies which allow the decomposition of the total intermolecular interaction energy into physically meaningful contributions, namely, the symmetry adapted perturbation theory (SAPT) and the use of localized orbitals within a Møller-Plesset perturbation scheme. The accuracy of the two different methods is compared to supermolecular results, within MP2 and coupled-cluster theory within single and double excitations, with perturbative estimates of the amplitudes of triple excitations [CCSD(T)]. Some relations between the different approaches are conjectured from theoretical considerations, and are confirmed by numerical results. The corresponding calculations have been performed for three model dimers: two NH[SUB3] … H[SUB2]O dimers, with NH[SUB3] acting once as a proton acceptor and once as a proton donor, and the NH[SUB4,SUP+]…H[SUB2]O considered as a prototype of the ion -- molecule interaction. We may conclude that third-order terms in SAPT help significantly to reproduce the Hartree -- Fock induction and the relaxed, total dispersion in the LMP2 decomposition. [ABSTRACT FROM AUTHOR]

We calculated the intermolecular potential for the complex methane-Ar in which the methane is deformed along the ν[sub 1] and the ν[sub 3] vibrational normal modes. We present a fit of the potential energy surface as a function of the three argon coordinates and the four normal mode coordinates q[sub ν][sub 3x], q[sub ν][sub 3y], q[sub ν][sub 3z], and q[sub ν][sub 1]. This potential is used in the following paper [J. Chem. Phys. 117, 7562 (2002)] dealing with the calculation of the lifetimes of the quasibound states of the dimer with the methane monomer in the vibrationally excited ν[sub 3] state. [ABSTRACT FROM AUTHOR]

Physisorption of water on the (111) surface of barium fluoride is investigated using accurate intermolecular potentials. A revised version of the successful ASP-W2 water potential is developed together with a new potential describing the interaction between water and the barium fluoride surface. Isolated water molecules are found to have a binding energy of -39.8 kJ mol[sup -1]. Monolayer and multilayer coverages are also investigated and compared with previous experimental and theoretical work. We find no evidence to support epitaxial growth of ferroelectric ice on this surface and suggest reasons for this. [ABSTRACT FROM AUTHOR]

Symmetry-adapted perturbation theory has been applied to compute the intermolecular potential energy surface of the He–CO[sub 2] complex. The ab initio potential has a global minimum of ε[sub m]=-50.38 cm[sup -1] at R[sub m]=5.81 bohr for the “T”-shaped geometry of the complex, and a local one of ε[sub m]=-28.94 cm[sup -1] at R[sub m]=8.03 bohr for the linear He···O==C==O structure. The computed potential energy surface has been analytically fitted and used in converged variational calculations to generate bound rovibrational states of the He–CO[sub 2] complex and the infrared spectrum corresponding to the simultaneous excitation of the ν[sub 3] vibration and internal rotation in the CO[sub 2] subunit within the complex. The complex was shown to be a semirigid asymmetric top and the rovibrational energy levels could be classified with the asymmetric top quantum numbers. The computed frequencies of the infrared transitions in the ν[sub 4] band of the spectrum are in very good agreement with the high resolution experimental data of Weida et al. [J. Chem. Phys. 101, 8351 (1994)]. The energy levels corresponding to the ν[sub 5] bending mode of the complex have been used to compute the transition frequencies in the ν[sub 5] hot band of He–CO[sub 2]. A tentative assignment of the transitions observed in the ν[sub 5] band with the quantum numbers of the asymmetric rotor is presented. As a further test of the ab initio potential we also report the pressure broadening coefficients of the R branch rotational lines of the ν[sub 3] spectrum of CO[sub 2] in a helium bath at various temperatures. Very good agreement is found with the wealth of experimental results for various rotational states of CO[sub 2] at different temperatures. Finally, we also tested the potential by computing the second virial coefficients at various temperatures. Again, the agreement between theory an... [ABSTRACT FROM AUTHOR]

Examines the applicability of various density functional theory (DFT) methods to describe the anisotropy of the intermolecular potential energy surfaces of hydrogen-bonded and van der Waals complexes. Comparison between supermolecule coupled-cluster method and perturbational symmetry-adapted perturbation theory.

Studies the nature of the pi-H interaction in ethene-H[sub 2]O, benzene-H[sub 2]O and benzene-(H[sub 2]O)[sub 2]. Minimum geometries and vibrational frequencies; Binding energy of the complexes; Results from symmetry adapted perturbation theory on the interaction of water with benzene and ethene.