On the Importance of Electron Correlation Effects for the Intramolecular Stacking Geometry of a Bis-Thiophene Derivative

Expanding the Range of Force Fields Available for ONIOM Calculations: The SICTWO Interface

The ONIOM scheme is one of the most popular QM/MM approaches, but its extended application has been so far hindered by the limited availability of force fields in most practical implementations. This paper describes a simple software code to overcome this limitation, and its application to three representative chemical problems. The "Shell Interface for Combining Tinker With ONIOM" (SICTWO) program gives access to all force fields available in the Tinker molecular mechanics program from a Gaussian09 or Gaussian16 calculation. The first application presented is the geometry optimization of dithienobicyclo-[4.4.1]-undeca-3,8-diene-11-one ethylene glycol ketal. A variety of force fields were tested for the MM part, and the molecular structure using the ONIOM(B3LYP:OPLS-AA) method shows good agreement with the experimental results. The second problem is related with enantioselectivity of the vanadium-catalyzed asymmetric oxidation of 1,2-bis( tert-butyl)-disulfide. ONIOM(B3LYP:MM3) results, obtained with SICTWO, are consistent with those of a previous IMOMM(B3LYP:MM3) study. In the third application, we have combined AMOEBA09 polarizable force field with ONIOM to study a water molecule binding on water ice. Calculated ONIOM(ωB97X-D:AMOEBA09) binding energies are consistent with the ωB97X-D results. These studies show SICTWO as an easy-to-use tool that can further expand the use of the multiscale ONIOM method within computational chemistry and computational biology.

Computational studies on the interactions of glycine amino acid with graphene, h-BN and h-SiC monolayers

In the present work, the adsorption of glycine amino acid and its zwitterionic form onto three different hexagonal sheets, namely graphene, boron-nitride (h-BN) and silicon carbide (h-SiC), has been investigated within the framework of density functional theory (DFT) calculations.

A charge density study of π-delocalization and intermolecular interactions

The location of bond critical points (red dots) and its associated bond path (black line) provide the evidence on the existence of the weak intermolecular interactions of the π–π interactions between triazole rings of atrz molecules in crystal with the close ring distance of 3.17 Å.

Comparison of some dispersion-corrected and traditional functionals with CCSD(T) and MP2 ab initio methods: dispersion, induction, and basis set superposition error

We compare dispersion and induction interactions for noble gas dimers and for Ne, methane, and 2-butyne with HF and LiF using a variety of functionals (including some specifically parameterized to evaluate dispersion interactions) with ab initio methods including CCSD(T) and MP2. We see that inductive interactions tend to enhance dispersion and may be accompanied by charge-transfer. We show that the functionals do not generally follow the expected trends in interaction energies, basis set superposition errors (BSSE), and interaction distances as a function of basis set size. The functionals parameterized to treat dispersion interactions often overestimate these interactions, sometimes by quite a lot, when compared to higher level calculations. Which functionals work best depends upon the examples chosen. The B3LYP and X3LYP functionals, which do not describe pure dispersion interactions, appear to describe dispersion mixed with induction about as accurately as those parametrized to treat dispersion. We observed significant differences in high-level wavefunction calculations in a basis set larger than those used to generate the structures in many of the databases. We discuss the implications for highly parameterized functionals based on these databases, as well as the use of simple potential energy for fitting the parameters rather than experimentally determinable thermodynamic state functions that involve consideration of vibrational states.

Importance of London dispersion effects for the packing of molecular crystals: a case study for intramolecular stacking in a bis-thiophene derivative

The influence of London dispersion effects in density functional theory for the computation of molecular crystals is investigated. The structure of dithienobicyclo[4.4.1]-undeca-3,8-dien-11-one ethylene glycol ketal (Resvan) that has been studied previously under isolated molecule conditions is considered as an example. The intramolecular stacking of the two thiophene rings is strongly influenced by non-covalent interactions. Previous investigations revealed significant deviations between computed molecular and experimentally observed crystal sulfur-sulfur distances. The influence of the crystal environment on this distance is estimated by calculating potential energy curves for the "gas phase" and in the crystal using the periodic plane-wave code VASP. With three common dispersion corrected (DFT-D3) density functionals (PBE, revPBE and PBEsol), an average shortening of the sulfur-sulfur distance by 0.10 A ("packing" effect) is computed which leads to a corrected "experimental" gas-phase value of R(S-S) = 4.19 A. This value is in good agreement with the results of reliable quantum chemical methods and serves as a challenging benchmark for electronic structure methods. Overall the crystal environment leads to a widening of the molecular structure due to intermolecular dispersion interaction. This is opposite to uncorrected DFT that suggest a compression due to Pauli exchange repulsion. In the crystal the intra- and inter-molecular dispersion effects are large but of opposite sign and their effect almost cancels out so that in this example dispersion-uncorrected DFT yields reasonable solid state structures. This is different for the computation of the lattice energy for which reasonable values are obtained only at the DFT-D3 level.

Atomic Cholesky decompositions: a route to unbiased auxiliary basis sets for density fitting approximation with tunable accuracy and efficiency

Cholesky decomposition of the atomic two-electron integral matrix has recently been proposed as a procedure for automated generation of auxiliary basis sets for the density fitting approximation [F. Aquilante et al., J. Chem. Phys. 127, 114107 (2007)]. In order to increase computational performance while maintaining accuracy, we propose here to reduce the number of primitive Gaussian functions of the contracted auxiliary basis functions by means of a second Cholesky decomposition. Test calculations show that this procedure is most beneficial in conjunction with highly contracted atomic orbital basis sets such as atomic natural orbitals, and that the error resulting from the second decomposition is negligible. We also demonstrate theoretically as well as computationally that the locality of the fitting coefficients can be controlled by means of the decomposition threshold even with the long-ranged Coulomb metric. Cholesky decomposition-based auxiliary basis sets are thus ideally suited for local density fitting approximations.

Microhydration of guanine...cytosine base pairs, a theoretical Study on the role of water in stability, structure and tautomeric equilibrium

The potential energy surfaces of guanine...cytosine complexes and microhydrated guanine...cytosine (one and two water molecules) were investigated by the molecular dynamics/quenching method (MD/Q), using the empirical potential Parm94 force field, implemented in the Amber program package. The calculations were conducted for all the possible combinations of the four most stable tautomers of guanine and three of cytosine (covering the canonical forms in both cases). The obtained structures were sorted by their structural motifs into three main groups: planar hydrogen-bonded; stacked; and T-shaped structures. The most stable structures found at the empirical potential energy surfaces were fully reoptimised at the second-order Møller-Plesset perturbation theory as well as using the density functional method with an empirical dispersion term (DFT-D). A combination of the canonical form of guanine and cytosine and canonical cytosine with a guanine tautomer where the hydrogen is switched from position N9 to N7 are energetically preferred in microsolvated systems as well as those without the presence of a solvent. The rising number of water molecules leads to smaller differences between the stability of the various combinations of the tautomers of bases in the base pairs. For some of the tautomer combinations (mainly the enol-enol combination), two water molecules are sufficient for the preference of stacked structures over the H-bonded ones. The interaction energies and geometries obtained by the second-order Møller-Plesset perturbation theory method and the much less computationally demanding DFT-D method are comparable, except for stacked complexes, where the interaction energies are overestimated on average by 3 kcal mol(-1) at the MP2 level.