Ab initio potential energy surface and vibration–rotation energy levels of germanium dicarbide,GeC2

The accurate ground‐state potential energy surface of germanium dicarbide, GeC₂, has been determined from ab initio calculations using the coupled‐cluster approach. The core–electron correlation, higher‐order valence‐electron correlation, and scalar relativistic effects were taken into account. The potential energy surface of GeC₂ was shown to be extraordinarily flat near the T‐shaped equilibrium configuration. The potential energy barrier to the linear CCGe configuration was predicted to be 1218 cm⁻¹. The vibration–rotation energy levels of some GeC₂ isotopologues were calculated using a variational method. The vibrational bending mode ν₃ was found to be highly anharmonic, with the fundamental wavenumber being only 58 cm⁻¹. Vibrational progressions due to this mode were predicted for the , , and states of GeC₂. © 2018 Wiley Periodicals, Inc.     

Low‐energy structures of benzene clusters with a novel accurate potential surface

The benzene‐benzene (Bz‐Bz) interaction is present in several chemical systems and it is known to be crucial in understanding the specificity of important biological phenomena. In this work, we propose a novel Bz‐Bz analytical potential energy surface which is fine‐tuned on accurate ab initio calculations in order to improve its reliability. Once the Bz‐Bz interaction is modeled, an analytical function for the energy of the clusters may be obtained by summing up over all pair potentials. We apply an evolutionary algorithm (EA) to discover the lowest‐energy structures of clusters (for ), and the results are compared with previous global optimization studies where different potential functions were employed. Besides the global minimum, the EA also gives the structures of other low‐lying isomers ranked by the corresponding energy. Additional ab initio calculations are carried out for the low‐lying isomers of and clusters, and the global minimum is confirmed as the most stable structure for both sizes. Finally, a detailed analysis of the low‐energy isomers of the n = 13 and 19 magic‐number clusters is performed. The two lowest‐energy isomers show S₆ and C₃ symmetry, respectively, which is compatible with the experimental results available in the literature. The structures reported here are all non‐symmetric, showing two central Bz molecules surrounded by 12 nearest‐neighbor monomers in the case of the five lowest‐energy structures. © 2015 Wiley Periodicals, Inc.     

Revealing the physical nature and the strength of charge‐inverted hydrogen bonds by SAPT(DFT), MP2, SCS‐MP2, MP2C,and CCSD(T) methods

The physical nature of charge‐inverted hydrogen bonds in H₃XH YH₃ (X = Si, Ge; Y = Al, Ga) dimer systems is studied by means of the SAPT(DFT)‐based decomposition of interaction energies and supermolecular interaction energies based on MP2, SCS‐MP2, MP2C, and CCSD(T) methods utilizing dimer‐centered aug‐cc‐pCVnZ (n = D, T, Q) basis sets as well as an extrapolation to the complete basis set limit. It is revealed that charge‐inverted hydrogen bonds are inductive in nature, although dispersion is also important. Computed interaction energies form the following relation: . It is confirmed that the aug‐cc‐pCVDZ basis set performs poorly and that very accurate values of interaction and dispersion energies require basis sets of at least quadrupole‐ζ quality. Considerably large binding energies suggest potential usefulness of charge‐inverted hydrogen bonds as an important structural motif in molecular binding. Terminology applying to σ‐ and π‐hole interactions as well as to triel and tetrel bonds is discussed. According to this new terminology the charge‐inverted hydrogen bond would become the first described case of a hydride‐triel bond. © 2017 Wiley Periodicals, Inc.     

Assessment of hydrophobic interactions and their contributions through the analysis of the methane dimer

Hydrophobic Interactions (HIs) are important in many phenomena of molecular recognition in chemistry and biology. Still, the relevance of HIs is sometimes difficult to evaluate particularly in large systems and intramolecular interactions. We put forward a method to estimate the magnitude and the different contributions of a given HI of the C···C, H? C···H, and H···H type through (i) the analysis of the electron density in the intermolecular region for eleven relative orientations of the methane dimer and (ii) the subsequent decomposition of the corresponding interaction energy in physically significant contributions using Symmetry Adapted Perturbation Theory (SAPT). Strong correlations were found between the topological properties of calculated at intermolecular bond critical points and plus its different contributions with the C···C distance of the considered orientations of (CH₄)₂. These correlations were used to construct Mollier‐like diagrams of and its components as a function of the separation between two carbons and the orientation of the groups bonded to these atoms. The ethane dimer and tert‐butylcyclohexane are used as representative examples of this new approach. Overall, we anticipate that this new method might prove useful in the study of both intramolecular and intermolecular HIs particularly of those within large systems wherein SAPT or electronic structure calculations are computationally expensive or even prohibitive. © 2014 Wiley Periodicals, Inc.     

Extension of accompanying coordinate expansion and recurrence relation method for general‐contraction basis sets

An algorithm of the accompanying coordinate expansion and recurrence relation (ACE‐RR), which is used for the rapid evaluation of the electron repulsion integral (ERI), has been extended to the general‐contraction (GC) scheme. The present algorithm, denoted by GC‐ACE‐RR, is designed for molecular calculations including heavy elements, whose orbitals consist of many primitive functions with and without higher angular momentum such as d‐ and f‐orbitals. The performance of GC‐ACE‐RR was assessed for ‐, ‐, ‐, and ‐type ERIs in terms of contraction length and the number of GC orbitals. The present algorithm was found to reduce the central processing unit time compared with the ACE‐RR algorithm, especially for higher angular momentum and highly contracted orbitals. Compared with HONDOPLUS and GAMESS program packages, GC‐ACE‐RR computations for ERIs of three‐dimensional gold clusters Au_n (n = 1, 2, …, 10, 15, 20, and 25) are more than 10 times faster. © 2014 Wiley Periodicals, Inc.     

A multipolar approach to the interatomic covalent interaction energy

Interatomic exchange‐correlation energies correspond to the covalent energetic contributions to an interatomic interaction in real space theories of the chemical bond, but their widespread use is severely limited due to their computationally intensive character. In the same way as the multipolar (mp ) expansion is customary used in biomolecular modeling to approximate the classical Coulomb interaction between two charge densities and , we examine in this work the mp approach to approximate the interatomic exchange‐correlation (xc) energies of the Interacting Quantum Atoms method. We show that the full xc mp series is quickly divergent for directly bonded atoms (1–2 pairs) albeit it works reasonably well most times for 1– n (n > 2) interactions. As with conventional perturbation theory, we show numerically that the xc series is asymptotically convergent and that, a truncated xc mp approximation retaining terms up to usually gives relatively accurate results, sometimes even for directly bonded atoms. Our findings are supported by extensive numerical analyses on a variety of systems that range from several standard hydrogen bonded dimers to typically covalent or aromatic molecules. The exact algebraic relationship between the monopole‐monopole xc mp term and the inter‐atomic bond order, as measured by the delocalization index of the quantum theory of atoms in molecules, is also established. © 2017 Wiley Periodicals, Inc.     

Quantum Monte Carlo with very large multideterminant wavefunctions

An algorithm to compute efficiently the first two derivatives of (very) large multideterminant wavefunctions for quantum Monte Carlo calculations is presented. The calculation of determinants and their derivatives is performed using the Sherman–Morrison formula for updating the inverse Slater matrix. An improved implementation based on the reduction of the number of column substitutions and on a very efficient implementation of the calculation of the scalar products involved is presented. It is emphasized that multideterminant expansions contain in general a large number of identical spin‐specific determinants: for typical configuration interaction‐type wavefunctions the number of unique spin‐specific determinants ( ) with a non‐negligible weight in the expansion is of order . We show that a careful implementation of the calculation of the N_det ‐dependent contributions can make this step negligible enough so that in practice the algorithm scales as the total number of unique spin‐specific determinants, , over a wide range of total number of determinants (here, N_det up to about one million), thus greatly reducing the total computational cost. Finally, a new truncation scheme for the multideterminant expansion is proposed so that larger expansions can be considered without increasing the computational time. The algorithm is illustrated with all‐electron fixed‐node diffusion Monte Carlo calculations of the total energy of the chlorine atom. Calculations using a trial wavefunction including about 750,000 determinants with a computational increase of ～400 compared to a single‐determinant calculation are shown to be feasible. © 2016 Wiley Periodicals, Inc.     

Quantum chemical study of the autoxidation of ascorbate

Reactions involved in the autoxidation of ascorbate have been investigated with quantum chemical first‐principles and ab initio methods. Reaction energies and Gibbs energies of the reactions were calculated at the density functional theory level applying the gradient‐corrected BP86 and the hybrid B3LYP functionals together with def2‐TZVP basis sets. Results of single‐point CC2, CCSD, and CCSD(T) calculations were used for calibration of the density functional theory data and show excellent agreement with the B3LYP values. Based on the Gibbs energy ascorbic acid AscH₂ is found to be the energetically lowest species in aqueous solution, whereas the monoanion ascorbate AscH is the most abundant one near pH = 7. Asc was found to be the preferred reducing agent for autoxidation and oxidation processes. The results also support a metal‐catalyzed synthesis of the reactive oxygen species H₂O₂ according to a redox cycling mechanism proposed in literature. © 2016 Wiley Periodicals, Inc.     

Aqueous acidities of primary benzenesulfonamides: Quantum chemical predictions based on density functional theory and SMD

Aqueous of selected primary benzenesulfonamides are predicted in a systematic manner using density functional theory methods and the SMD solvent model together with direct and proton exchange thermodynamic cycles. Some test calculations were also performed using high‐level composite CBS‐QB3 approach. The direct scheme generally does not yield a satisfactory agreement between calculated and measured acidities due to a severe overestimation of the Gibbs free energy changes of the gas‐phase deprotonation reaction by the used exchange‐correlation functionals. The relative values calculated using proton exchange method compare to experimental data very well in both qualitative and quantitative terms, with a mean absolute error of about 0.4 units. To achieve this accuracy, we find it mandatory to perform geometry optimization of the neutral and anionic species in the gas and solution phases separately, because different conformations are stabilized in these two cases. We have attempted to evaluate the effect of the conformer‐averaged free energies in the predictions, and the general conclusion is that this procedure is highly too costly as compared with the very small improvement we have gained. © 2015 Wiley Periodicals, Inc.     

Structure,stability, and nature of bonding in carbon monoxide bound complexes (E = group 14 element; X = H,F, Cl,Br, I)

A density functional theory study is performed to predict the structures and stability of carbon monoxide (CO) bound (E = C, Si, Ge, Sn, Pb; X = H, F, Cl, Br, I) complexes. The possibility of bonding through both C‐ and O‐sides of CO is considered. Thermochemical analysis reveals that all the dissociation processes producing CO and are endothermic in nature whereas most of the dissociation reactions are endergonic in nature at room temperature. The nature of bonding in E? C/O bonds is analyzed via Wiberg bond index, natural population analysis, electron density, and energy decomposition analyses in conjunction with natural orbitals for chemical valence scheme. In comparison to C? O stretching frequency ( ) in free CO, while a red shift is noted in O‐side binding, the C‐side binding results in a blue shift in . The relative change in values in CO bound complexes on changing either E or X is scrutinized and possible explanation is provided in terms of polarization in the σ‐ and π‐orbitals and the relative strength of C→E or O→E σ‐donation and E→C or E→O π‐back‐donation. © 2016 Wiley Periodicals, Inc.     

Polygonal current models for polycyclic aromatic hydrocarbons and graphene sheets of various shapes

Assuming that graphene is an “infinite alternant” polycyclic aromatic hydrocarbon resulting from tessellation of a surface by only six‐membered carbon rings, planar fragments of various size and shape (hexagon, triangle, rectangle, and rhombus) have been considered to investigate their response to a magnetic field applied perpendicularly. Allowing for simple polygonal current models, the diatropicity of a series of polycyclic textures has been reliably determined by comparing quantitative indicators, the π‐electron contribution to I^B, the magnetic field‐induced current susceptibility of the peripheral circuit, to and to , respectively the out‐of‐plane components of the magnetizability tensor and of the magnetic shielding tensor at the center of mass. Extended numerical tests and the analysis based on the polygonal model demonstrate that (i) and yield inadequate and sometimes erroneous measures of diatropicity, as they are heavily flawed by spurious geometrical factors, (ii) I^B values computed by simple polygonal models are valid quantitative indicators of aromaticity on the magnetic criterion, preferable to others presently available, whenever current susceptibility cannot be calculated ab initio as a flux integral, (iii) the hexagonal shape is the most effective to maximize the strength of π‐electron currents over the molecular perimeter, (iv) the edge current strength of triangular and rhombic graphene fragments is usually much smaller than that of hexagonal ones, (v) doping by boron and nitrogen nuclei can regulate and even inhibit peripheral ring currents, (vi) only for very large rectangular fragments can substantial current strengths be expected. © 2017 Wiley Periodicals, Inc.     

Update on phosphate and charged post‐translationally modified amino acid parameters in the GROMOS force field

In this study, we propose newly derived parameters for phosphate ions in the context of the GROMOS force field parameter sets. The non‐bonded parameters used up to now lead to a hydration free energy, which renders the dihydrogen phosphate ion too hydrophobic when compared to experimentally derived values, making a reparametrization of the phosphate moiety necessary. Phosphate species are of great importance in biomolecular simulations not only because of their crucial role in the backbone of nucleic acids but also as they represent one of the most important types of post‐translational modifications to protein side‐chains and are an integral part in many lipids. Our re‐parametrization of the free dihydrogen phosphate (H PO ) and three derivatives (methyl phosphate, dimethyl phosphate, and phenyl phosphate) leads, in conjunction with the previously updated charged side‐chains in the GROMOS parameter set 54A8, to new nucleic acid backbone parameters and a 54A8 version of the widely used GROMOS protein post‐translational modification parameter set. © 2017 Wiley Periodicals, Inc.     

Quantum supercharger library: Hyper‐parallel integral derivatives algorithms for ab initio QM/MM dynamics

This article describes an extension of the quantum supercharger library (QSL) to perform quantum mechanical (QM) gradient and optimization calculations as well as hybrid QM and molecular mechanical (QM/MM) molecular dynamics simulations. The integral derivatives are, after the two‐electron integrals, the most computationally expensive part of the aforementioned calculations/simulations. Algorithms are presented for accelerating the one‐ and two‐electron integral derivatives on a graphical processing unit (GPU). It is shown that a Hartree–Fock ab initio gradient calculation is up to 9.3X faster on a single GPU compared with a single central processing unit running an optimized serial version of GAMESS‐UK, which uses the efficient Schlegel method for ‐ and ‐orbitals. Benchmark QM and QM/MM molecular dynamics simulations are performed on cellobiose in vacuo and in a 39 Å water sphere (45 QM atoms and 24843 point charges, respectively) using the 6‐31G basis set. The QSL can perform 9.7 ps/day of ab initio QM dynamics and 6.4 ps/day of QM/MM dynamics on a single GPU in full double precision. © 2015 Wiley Periodicals, Inc.     

Structure and electrochemical properties for complexes of nitrocompounds with inorganic ions: A theoretical approach

Reduction and oxidation (redox) reactions are widely used for removal of nitrocompounds from contaminated soil and water. Structures and redox properties for complexes of nitrocompounds, such as 2,4,6‐trinitrotoluene (TNT), 2,4‐dinitrotoluene (DNT), 2,4‐dinitroanisole (DNAN), and 5‐nitro‐2,4‐dihydro‐3H?1,2,4‐triazol‐3‐one (NTO), with common inorganic ions (Na⁺, Cl^?, ) were investigated at the SMD(Pauling)/PCM(Pauling)/MPWB1K/TZVP level of theory. Atoms in molecules (AIM) theory was applied to analyze the topological properties of the bond critical points involved in the interactions between the nitrocompounds and the ions. Topological analyses show that intermolecular interactions of the types O(N)…Na⁺, C‐H…Cl^?( ), and C…Cl^?( ) may be discussed as noncovalent closed‐shell interactions, while N‐H···Cl^?( ) hydrogen bonds are partially covalent in nature. Complexation causes significant decrease of redox activity of the nitrocompounds. Analysis of the reduction potentials of the complexes obtained through application of the Pourbaix diagram of an iron/water system revealed that sodium complexes of NTO might be reduced by metallic iron. © 2016 Wiley Periodicals, Inc.     

Force constant decomposition for penta‐coordinated XH3Cl2‐ (X = C,Si, Ge) structures

Based on the energy decomposition analysis of an interacting system, we propose a method for force constant decomposition analysis with respect to the specific normal coordinate. Using the presented method, we examined the penta‐coordinated system (X = C, Si, Ge), which possesses a three‐center four‐electron bond. The origin of the difference in the stability of the penta‐coordinated D_3h structures was clearly shown to be the effect of electron delocalization–polarization term. © 2018 Wiley Periodicals, Inc.     

The solvent effect on two competing reaction mechanisms involving hypervalent iodine reagents (λ3‐iodanes): Facing the limit of the stationary quantum chemical approach

Trifluoromethylation of acetonitrile with 3,3‐dimethyl‐1‐(trifluoromethyl)?1λ³,2‐ benziodoxol is assumed to occur via reductive elimination (RE) of the electrophilic CF₃‐ligand and MeCN bound to the hypervalent iodine. Computations in gas phase showed that the reaction might also occur via an S_N2 mechanism. There is a substantial solvent effect present for both reaction mechanisms, and their energies of activation are very sensitive toward the solvent model used (implicit, microsolvation, and cluster‐continuum). With polarizable continuum model‐based methods, the S_N2 mechanism becomes less favorable. Applying the cluster‐continuum model, using a shell of solvent molecules derived from ab initio molecular dynamics (AIMD) simulations, the gap between the two activation barriers ( ) is lowered to a few kcal mol^?1 and also shows that the activation entropies ( ) and volumes ( ) for the two mechanisms differ substantially. A quantitative assessment of will therefore only be possible using AIMD. A natural bond orbital‐analysis gives further insight into the activation of the CF₃‐reagent by protonation. © 2014 Wiley Periodicals, Inc.     

Why it is sometimes difficult to determine the accurate position of a hydrogen atom by the semiexperimental method: Structure of molecules containing the OH or the CH3 group

The semiexperimental (SE) technique, whereby equilibrium rotational constants are derived from experimental ground‐state rotational constants and corrections based on an ab initio cubic force field, has the reputation to be one of the most accurate methods to determine an equilibrium structure ( ). However, in some cases, it cannot determine accurately the position of the hydrogen. To investigate the origins of this difficulty, the SE structures of several molecules containing either the OH or the CH₃ group are determined and compared to their best ab initio counterparts. It appears that an important factor is the accuracy of the geometry used to calculate the force field, in particular when the least‐squares system is not well conditioned. In this case, the mixed regression method is often an easy way to circumvent this difficulty. © 2014 Wiley Periodicals, Inc.     

GAMPMS: Genetic algorithm managed peptide mutant screening

The prominence of endogenous peptide ligands targeted to receptors makes peptides with the desired binding activity good molecular scaffolds for drug development. Minor modifications to a peptide's primary sequence can significantly alter its binding properties with a receptor, and screening collections of peptide mutants is a useful technique for probing the receptor–ligand binding domain. Unfortunately, the combinatorial growth of such collections can limit the number of mutations which can be explored using structure‐based molecular docking techniques. Genetic algorithm managed peptide mutant screening (GAMPMS) uses a genetic algorithm to conduct a heuristic search of the peptide's mutation space for peptides with optimal binding activity, significantly reducing the computational requirements of the virtual screening. The GAMPMS procedure was implemented and used to explore the binding domain of the nicotinic acetylcholine receptor (nAChR) ‐isoform with a library of 64,000 α‐conotoxin (α‐CTx) MII peptide mutants. To assess GAMPMS's performance, it was compared with a virtual screening procedure that used AutoDock to predict the binding affinity of each of the α‐CTx MII peptide mutants with the ‐nAChR. The GAMPMS implementation performed AutoDock simulations for as few as 1140 of the 64,000 α‐CTx MII peptide mutants and could consistently identify a set of 10 peptides with an aggregated binding energy that was at least 98% of the aggregated binding energy of the 10 top peptides from the exhaustive AutoDock screening. © 2015 Wiley Periodicals, Inc.     

Four faces of the interaction between ions and aromatic rings

Non‐covalent interactions between ions and aromatic rings play an important role in the stabilization of macromolecular complexes; of particular interest are peptides and proteins containing aromatic side chains (Phe, Trp, and Tyr) interacting with negatively (Asp and Glu) and positively (Arg and Lys) charged amino acid residues. The structures of the ion–aromatic‐ring complexes are the result of an interaction between the large quadrupole moment of the ring and the charge of the ion. Four attractive interaction types are proposed to be distinguished based on the position of the ion with respect to the plane of the ring: perpendicular cation–π (CP_⊥), co‐planar cation–π (CP_∥), perpendicular anion–π (AP_⊥), and co‐planar anion–π (AP_∥). To understand more than the basic features of these four interaction types, a systematic, high‐level quantum chemical study is performed, using the X^– + C₆H₆, M⁺ + C₆H₆, X^– + C₆F₆, and M⁺ + C₆F₆ model systems with X⁻ = H⁻, F⁻, Cl⁻, HCOO⁻, CH₃COO⁻ and M⁺ = H⁺, Li⁺, Na⁺, , CH₃ , whereby C₆H₆ and C₆F₆ represent an electron‐rich and an electron‐deficient π system, respectively. Benchmark‐quality interaction energies with small uncertainties, obtained via the so‐called focal‐point analysis (FPA) technique, are reported for the four interaction types. The computations reveal that the interactions lead to significant stabilization, and that the interaction energy order, given in kcal mol⁻¹ in parentheses, is CP_⊥ (23–37) > AP_⊥ (14–21) > CP_∥ (9–22) > AP_∥ (6–16). A natural bond orbital analysis performed leads to a deeper qualitative understanding of the four interaction types. To facilitate the future quantum chemical characterization of ion–aromatic‐ring interactions in large biomolecules, the performance of three density functional theory methods, B3LYP, BHandHLYP, and M06‐2X, is tested against the FPA benchmarks, with the result that the M06‐2X functional performs best. © 2017 Wiley Periodicals, Inc.