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.     

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.     

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.     

Spin–orbit effects on magnetically induced current densities in the () clusters

This study reports the spin–orbit effects on the aromaticity of the , , , , , and anionic clusters via the magnetically induced current‐density method. All‐electron density functional theory (DFT) calculations were carried out using the four‐component Dirac‐Coulomb (DC) hamiltonian, including scalar and spin–orbit relativistic effects. The magnetic index of aromaticity was calculated by numerical integration over the current flow between two atoms in the pentagonal ring. These values were compared to the spin‐free values (spin–orbit coupling switched off), in order to assess the spin–orbit effect on aromaticity. It was found that in the heavy anions, and , there is a significant influence of the spin–orbit coupling. © 2018 Wiley Periodicals, Inc.     

Calculation of photoionization differential cross sections using complex Gauss‐type orbitals

Accurate theoretical calculation of photoelectron angular distributions for general molecules is becoming an important tool to image various chemical reactions in real time. We show in this article that not only photoionization total cross sections but also photoelectron angular distributions can be accurately calculated using complex Gauss‐type orbital (cGTO) basis functions. Our method can be easily combined with existing quantum chemistry techniques including electron correlation effects, and applied to various molecules. The so‐called two‐potential formula is applied to represent the transition dipole moment from an initial bound state to a final continuum state in the molecular coordinate frame. The two required continuum functions, the zeroth‐order final continuum state and the first‐order wave function induced by the photon field, have been variationally obtained using the complex basis function method with a mixture of appropriate cGTOs and conventional real Gauss‐type orbitals (GTOs) to represent the continuum orbitals as well as the remaining bound orbitals. The complex orbital exponents of the cGTOs are optimized by fitting to the outgoing Coulomb functions. The efficiency of the current method is demonstrated through the calculations of the asymmetry parameters and molecular‐frame photoelectron angular distributions of and . In the calculations of , the static exchange and random phase approximations are employed, and the dependence of the results on the basis functions is discussed. © 2017 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.     

Fast and accurate determination of the relative binding affinities of small compounds to HIV‐1 protease using non‐equilibrium work

The fast pulling ligand (FPL) out of binding cavity using non‐equilibrium molecular dynamics (MD) simulations was demonstrated to be a rapid, accurate and low CPU demand method for the determination of the relative binding affinities of a large number of HIV‐1 protease (PR) inhibitors. In this approach, the ligand is pulled out of the binding cavity of the protein using external harmonic forces, and the work of pulling force corresponds to the relative binding affinity of HIV‐1 PR inhibitor. The correlation coefficient between the pulling work and the experimental binding free energy of shows that FPL results are in good agreement with experiment. It is thus easier to rank the binding affinities of HIV‐1 PR inhibitors, that have similar binding affinities because the mean error bar of pulling work amounts to . The nature of binding is discovered using the FPL approach. © 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.     

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.     

Molecules under external electric field: On the changes in the electronic structure and validity limits of the theoretical predictions

A detailed analysis of the electronic structure of the ground and first excited spin state of three diatomic molecules ( and ) under static applied electric field is performed at CCSD(T), DFT, MRCI and MRCI(Q) levels of theory. Our findings have revealed that by boosting the applied field one induces changes in the occupation numbers of molecular orbitals, giving rise to changes in the equilibrium geometry and in the HOMO–LUMO energy gap. Specifically, singlet to triplet spin transition can be induced by increasing the applied electric field beyond a critical value. Accordingly, affecting the accuracy of the widely used expression of energy expanded in Taylor series with respect to the applied electric field. © 2018 Wiley Periodicals, Inc.     

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.     

Exploring possible reaction pathways for the o‐atom transfer reactions to unsaturated substrates catalyzed by a [Ni‐NO2] ↔ [Ni‐NO] redox couple using DFT methods

The (nitro)(N‐methyldithiocarbamato)(trimethylphospane)nickel(II), [Ni(NO₂)(S₂CNHMe)(PMe₃)] complex catalyses efficiently the O‐atom transfer reactions to CO and acetylene. Energetically feasible sequence of elementary steps involved in the catalytic cycle of the air oxidation of CO and acetylene are proposed promoted by the Ni(NO₂)(S₂CNHMe)(PMe₃)] ↔ Ni(NO₂)(S₂CNHMe)(PMe₃) redox couple using DFT methods both in vacuum and dichloromethane solutions. The catalytic air oxidation of HC≡CH involves formation of a five‐member metallacycle intermediate, via a [3 + 2] cyclo‐addition reaction of HC≡CH to the Ni‐N = O moiety of the Ni(NO₂)(S₂CNHMe)(PMe₃)] complex, followed by a β H‐atom migration toward the C_α carbon atom of the coordinated acetylene and release of the oxidation product (ketene). The geometric and energetic reaction profile for the reversible [Ni( ‐NO₂)(S₂CNHMe)(PMe₃)] [Ni( ‐ONO)(S₂CNHMe)(PMe₃)] linkage isomerization has also been modeled by DFT calculations. © 2017 Wiley Periodicals, Inc.     

Development of spin‐orbit coupling for stochastic configuration interaction techniques

To perform spin‐orbit coupling calculations on atoms and molecules, good zeroth‐order wavefunctions are necessary. Here, we present the software development of the Monte Carlo Configuration Interaction (MCCI) method, to enable calculation of such properties, where MCCI iteratively constructs a multireference wavefunction using a stochastic procedure. In this initial work, we aim to establish the efficacy of this technique in predicting the splitting of otherwise degenerate energy levels on a range of atoms and small diatomic molecules. It is hoped that this work will subsequently act as a gateway toward using this method to investigate singlet‐triplet interactions in larger multireference molecules. We show that MCCI can generate very good results using highly compact wavefunctions compared to other techniques, with no prior knowledge of important orbitals. Higher‐order relativistic effects are neglected and spin‐orbit coupling effects are incorporated using first‐order degenerate perturbation theory with the Breit‐Pauli Hamiltonian and effective nuclear charges in the one‐electron operator. Results are obtained and presented for B, C, O, F, Si, S, and Cl atoms and OH, CN, NO, and C₂ diatomic radicals including spin‐orbit coupling constants and the relative splitting of the lowest energy degenerate state for each species. Convergence of MCCI to the full configuration interaction result is demonstrated on the multireference problem of stretched OH. We also present results from the singlet‐triplet interaction between the and both the and states of the O₂ molecule. © 2017 Wiley Periodicals, Inc.     

Nuclear quantum effects induce metallization of dense solid molecular hydrogen

We present an accurate computational study of the electronic structure and lattice dynamics of solid molecular hydrogen at high pressure. The band‐gap energies of the , Pc, and structures at pressures of 250, 300, and 350 GPa are calculated using the diffusion quantum Monte Carlo (DMC) method. The atomic configurations are obtained from ab initio path‐integral molecular dynamics (PIMD) simulations at 300 K and 300 GPa to investigate the impact of zero‐point energy and temperature‐induced motion of the protons including anharmonic effects. We find that finite temperature and nuclear quantum effects reduce the band‐gaps substantially, leading to metallization of the and Pc phases via band overlap; the effect on the band‐gap of the structure is less pronounced. Our combined DMC‐PIMD simulations predict that there are no excitonic or quasiparticle energy gaps for the and Pc phases at 300 GPa and 300 K. Our results also indicate a strong correlation between the band‐gap energy and vibron modes. This strong coupling induces a band‐gap reduction of more than 2.46 eV in high‐pressure solid molecular hydrogen. Comparing our DMC‐PIMD with experimental results available, we conclude that none of the structures proposed is a good candidate for phases III and IV of solid hydrogen. © 2017 Wiley Periodicals, Inc.     

Fast search algorithms for computational protein design

One of the main challenges in computational protein design (CPD) is the huge size of the protein sequence and conformational space that has to be computationally explored. Recently, we showed that state‐of‐the‐art combinatorial optimization technologies based on Cost Function Network (CFN) processing allow speeding up provable rigid backbone protein design methods by several orders of magnitudes. Building up on this, we improved and injected CFN technology into the well‐established CPD package Osprey to allow all Osprey CPD algorithms to benefit from associated speedups. Because Osprey fundamentally relies on the ability of to produce conformations in increasing order of energy, we defined new strategies combining CFN lower bounds, with new side‐chain positioning‐based branching scheme. Beyond the speedups obtained in the new ‐CFN combination, this novel branching scheme enables a much faster enumeration of suboptimal sequences, far beyond what is reachable without it. Together with the immediate and important speedups provided by CFN technology, these developments directly benefit to all the algorithms that previously relied on the DEE/ combination inside Osprey* and make it possible to solve larger CPD problems with provable algorithms. © 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.     

FEREBUS: Highly parallelized engine for kriging training

FFLUX is a novel force field based on quantum topological atoms, combining multipolar electrostatics with IQA intraatomic and interatomic energy terms. The program FEREBUS calculates the hyperparameters of models produced by the machine learning method kriging. Calculation of kriging hyperparameters ( θ and p ) requires the optimization of the concentrated log‐likelihood . FEREBUS uses Particle Swarm Optimization (PSO) and Differential Evolution (DE) algorithms to find the maximum of . PSO and DE are two heuristic algorithms that each use a set of particles or vectors to explore the space in which is defined, searching for the maximum. The log‐likelihood is a computationally expensive function, which needs to be calculated several times during each optimization iteration. The cost scales quickly with the problem dimension and speed becomes critical in model generation. We present the strategy used to parallelize FEREBUS, and the optimization of through PSO and DE. The code is parallelized in two ways. MPI parallelization distributes the particles or vectors among the different processes, whereas the OpenMP implementation takes care of the calculation of , which involves the calculation and inversion of a particular matrix, whose size increases quickly with the dimension of the problem. The run time shows a speed‐up of 61 times going from single core to 90 cores with a saving, in one case, of ～98% of the single core time. In fact, the parallelization scheme presented reduces computational time from 2871 s for a single core calculation, to 41 s for 90 cores calculation. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

An anisotropic coarse‐grained model based on Gay–Berne and electric multipole potentials and its application to simulate a DMPC bilayer in an implicit solvent model

In this work, we aim at optimizing the performance of the anisotropic GBEMP model, which adopts a framework by combining a Gay–Berne (GB) anisotropic potential with an electric multipole (EMP) potential, in simulating a DMPC lipid bilayer in an implicit solvent model. First, the Gay–Berne parameters were initially obtained by fitting to atomistic profiles of van der Waals interactions between homodimers of molecular fragments while EMP parameters was directly derived from the expansion of point multipoles at predefined EMP sites. Second, the GB and EMP parameters for DMPC molecule were carefully optimized to be comparable to AMBER atomistic model in the calculations of the dipole moments of DMPC monomers adopting different conformations as well as the nonbonded interactions between two DMPC molecules adopting different conformations and separated at various distances. Finally, the GB parameters for DMPC were slightly adjusted in simulating a 72 DMPC bilayer system so that our GBEMP model would be able to reproduce a few important structural properties, namely, thickness (), area per lipid ( ) and volume per lipid ( ). Meanwhile, the atomistic and experimental results for electron density profiles and order parameters were reproduced reasonably well by the GBEMP model, demonstrating the promising feature of GBEMP model in modeling lipid systems. Finally, we have shown that current GBEMP model is more efficient by a factor of about 25 than AMBER atomistic point charge model. © 2015 Wiley Periodicals, Inc.     

Computational study of basis set and electron correlation effects on anapole magnetizabilities of chiral molecules

In the presence of a static, nonhomogeneous magnetic field, represented by the axial vector at the origin of the coordinate system and by the polar vector , assumed to be spatially uniform, the chiral molecules investigated in this paper carry an orbital electronic anapole, described by the polar vector . The electronic interaction energy of these molecules in nonordered media is a cross term, coupling and via , one third of the trace of the anapole magnetizability a_αβ tensor, that is, . Both and W^BC have opposite sign in the two enantiomeric forms, a fact quite remarkable from the conceptual point of view. The magnitude of predicted in the present computational investigation for five chiral molecules is very small and significantly biased by electron correlation contributions, estimated at the density functional level via three different functionals. © 2016 Wiley Periodicals, Inc.