Water‐mediated intermolecular interactions in 1,2‐O‐cyclohexylidene‐myo‐inositol: a quantitative analysis

The syntheses of new myo‐inositol derivatives have received much attention due to their important biological activities. 1,2‐O‐Cyclohexylidene‐myo‐inositol is an important intermediate formed during the syntheses of certain myo‐inositol derivatives. We report herein the crystal structure of 1,2‐O‐cyclohexylidene‐myo‐inositol dihydrate, C₁₂H₂₀O₆·2H₂O, which is an intermediate formed during the syntheses of myo‐inositol phosphate derivatives, to demonstrate the participation of water molecules and hydroxy groups in the formation of several intermolecular O—H…O interactions, and to determine a low‐energy conformation. The title myo‐inositol derivative crystallizes with two water molecules in the asymmetric unit in the space group C 2/c , with Z = 8. The water molecules facilitate the formation of an extensive O—H…O hydrogen‐bonding network that assists in the formation of a dense crystal packing. Furthermore, geometrical optimization and frequency analysis was carried out using density functional theory (DFT) calculations with B3LYP hybrid functionals and 6‐31G(d), 6‐31G(d,p) and 6‐311G(d,p) basis sets. The theoretical and experimental structures were found to be very similar, with only slight deviations. The intermolecular interactions were quantitatively analysed using Hirshfeld surface analysis and 2D (two‐dimensional) fingerplot plots, and the total lattice energy was calculated.     

Accurate Structure and Bonding Description of the Transition Metal‐Disulfur Monoxide Complexes [(PMe3)2M(S2O)] (M = Ni,Pd, Pt): Grimme Dispersion Corrected DFT Study

Geometry and bonding energy analysis of M–S₂O bonds in the metal‐disulfur monoxide complexes [(PMe₃)₂M(S₂O)] of nickel, palladium, and platinum were investigated at DFT, DFT‐D3, and DFT‐D3(BJ) methods using three different functionals (BP86, PBE, and TPSS). The TPSS/DFT‐D3(BJ) yields better geometry, while the BP86 geometry is least accurate for studied complexes. The geometry of platinum complex optimized at TPSS/DFT‐D3(BJ) level is in excellent agreement with the available experimental values. The M–S bonds are shorter than the M–S(O) bonds. The Mayer bond orders suggest the presence of M–S and M–S(O) single bonds. Both the M–S and M–S(O) bond lengths vary with the density functionals as TPSS‐D3(BJ) < TPSS < PBE < BP86. The Hirshfeld charge distribution indicates that the overall charge flows from metal fragment to [S₂O]. The Ni–S₂O bond has greater degree of covalent character than the ionic. The contribution of dispersion interactions is large in computing accurate bond dissociation energies between the interacting fragments. The BDEs are largest for the functional TPSS and smallest for the functional BP86. The DFT‐D3 dispersion corrections to the BDEs between the metal fragments [(PMe₃)₂M] and ligand fragment [(S₂O)] for the TPSS functional are in the range 7.1–7.3 kcal · mol^–1, which are smaller than the corresponding DFT‐D3(BJ) dispersion corrections (9.4–10.6 kcal · mol^–1).     

Assessment of dispersion‐improved exchange‐correlation functionals for the simulation of CO2 binding by alcoholamines

In this study, 12 bound complexes were selected to construct a database for testing 15 dispersion‐improved exchange‐correlation (XC) functionals, including hybrid generalized gradient approximation (GGA), modified using the Grimme's pairwise strategy, and double hybrid XC functionals, for specifically characterizing the CO₂ binding by alcoholamines. Bound complexes were selected based on the characteristics of their hydrogen bonds, dispersion, and electrostatic (particularly between the positive charge of CO₂ and the lone pair of N of alcoholamines) interactions. The extrapolated binding energy from the aug‐cc‐pVTZ (ATZ) to aug‐cc‐pVQZ (AQZ) basis set at the CCSD(T)/CBS(MP2+DZ) level was used as the reference for the XC functional comparison. M06‐2X produced the optimal agreement if the optimized geometries at MP2/ATZ level were chosen for all the test bound complexes. However, M06‐L, ωB97X, and ωB97, and were preferred if the corresponding density functional theory (DFT) optimized geometries were adapted for the benchmark. Simple bimolecular reaction between CO₂ and monoethanolamine simulated using polarizable continuum solvation model confirmed that ωB97, ωB97X, and ωB97XD qualitatively reproduced the energetics of MP2 level. The inconsistent performance of the tested XC functionals, observed when using MP2 or DFT optimized geometries, raised concerns regarding using the single‐point ab initio correction combined with DFT optimized geometry, particularly for determining the nucleophilic attack by alcoholamines to CO₂. © 2014 Wiley Periodicals, Inc.     

Density Functional Study on Adsorption of NO on AuSe (010) Surface

NO molecule adsorption on (010) surface of gold selenide (AuSe) has been studied with a periodic slab model by means of the GGA‐PW91 exchange‐correlation functional within the framework of density functional theory (DFT). Four different on‐top adsorption sites Au(1), Au(2), Se(1) and Se(2) were considered for α‐AuSe and three on‐top adsorption sites Au(1), Au(2) and Se(1) for β‐AuSe. N‐end and O‐end adsorptions of NO were investigated for the above sites. The results show that N‐end adsorptions are preferred for α‐ and β‐AuSe and O‐end adsorptions are not feasible and thought as physisorption with the weak adsorption energies from 6.0 to 10.8 kJ/mol. For the N‐end adsorptions on α‐ and β‐AuSe (010) surfaces, Au(2) sites are most favorable with the adsorption energies 89.0 and 78.0 kJ/mol for α‐ and β‐AuSe, respectively. However, the adsorptions at Au1 sites are very weak with the adsorption energies of 27.8 and 7.5 kJ/mol, respectively. In case of the adsorption of N‐down orientations of NO at Se sites for α‐ and β‐AuSe (010) surfaces, the adsorption activities of Se(1) and Se(2) sites on the α‐AuSe (010) surface and Se(1) site on the β‐AuSe (010) surface are almost the same with the adsorption energies 51.2, 52.7 and 49.2 kJ/mol. The geometric optimizations for adsorption configurations were calculated along with accounting for stretching frequency and density of states in our work.     

Generation of ground-state structures and electronic properties of ternary AlxTiyNiz clusters (x + y + z = 6) with a two-stage density functional theory global search approach

The structural and electronic properties of ternary Al_xTi_yNi_z clusters, where x, y, and z are integers and x + y + z = 6 , are investigated. Both Slater, Vosko, Wilks, and Nusair and B3LYP exchange-correlation (XC) functionals are employed in a two-stage density functional theory (DFT) calculations to generate these clusters. In the first stage, a minimum energy cluster structure is generated by an unbiased global search algorithm coupled with a DFT code using a light XC functional and small basis sets. In the second stage, the obtained cluster structure is further optimized by another round of global minimization search coupled with a DFT calculator using a heavier XC functional and more costly basis set. Electronic properties of the structures are illustrated in the form of a ternary diagram. Our DFT calculations find that the thermodynamic stability of the clusters increases with the increment in the number of constituent nickel atoms. These results provide a new insight to the structure, stability, chemical order, and electronic properties for the ternary alloy nanoclusters.     

U+与CO2气相反应的密度泛函理论研究

运用密度泛函理论(DFT)中的B3LYP方法,U原子用含相对论有效原子实势(ECP)校正的基组(SDD),C、O原子采用6-311+G(d)基组,对气相中U+和CO2的反应进行了理论研究.通过研究二重和四重自旋态的反应势能面(PESs),优化得到了两条反应路径的反应物、中间体、过渡态和产物的结构.用"两态反应"(TSR)分析反应机理,结果表明体系的优先选择路径为高自旋态进入和低自旋态离开反应,发生在四重态和二重态的自旋多重度的改变使得整个反应系统能以一个低能反应途径进行.     

Computational study toward understanding the photodissociation mechanism of sarin

The potential energy surface for sarin (C₄H₁₀FO₂P) dissociation into (CH₃)₂CHO+PO(F)(CH₃), CH₃+(CH₃)₂CHOPO(F), and F+(CH₃)₂CHOPO(CH₃) in the T₁ and S₁ states were investigated at the complete‐active‐space self‐consistent field (CASSCF) with the 6‐31G** and aug‐cc‐PVDZ basis sets. The different reaction pathways are characterized on the basis of the computed potential energy surface and surface crossing point, the time‐dependent density functional theory (TD‐DFT) was used to calculate the vertical energies based on the CAS(8,7)/6‐31G** optimized excited structures, which may provide some new insights into the mechanism of the ultraviolet photo‐degradation of sarin molecules. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Which density functional is close to CCSD accuracy to describe geometry and interaction energy of small noncovalent dimers? A benchmark study using Gaussian09

A benchmark study on all possible density functional theory (DFT) methods in Gaussian09 is done to locate functionals that agree well with CCSD/aug‐cc‐pVTZ geometry and Ave‐CCSD(T)/(Q‐T) interaction energy (E_int) for small non‐covalently interacting molecular dimers in “dispersion‐dominated” (class 1), “dipole‐induced dipole” (class 2), and “dipole‐dipole” (class 3) classes. A DFT method is recommended acceptable if the geometry showed close agreement to CCSD result (RMSD < 0.045) and E_int was within 80–120% accuracy. Among 382 tested functionals, 1–46% gave good geometry, 13–44% gave good E_int, while 1–33% satisfied geometry and energy criteria. Further screening to locate the best performing functionals for all the three classes was made by counting the acceptable values of energy and geometry given by each functionals. The meta‐generalized gradient approximation (GGA) functional M06L was the best performer with total 14 hits; seven acceptable energies and seven acceptable geometries. This was the only functional “recommended” for at least two dimers in each class. The functionals M05, B2PLYPD, B971, mPW2PLYPD, PBEB95, and CAM‐B3LYP gave 11 hits while PBEhB95, PW91B95, Wb97x, BRxVP86, BRxP86, HSE2PBE, HSEh1PBE, PBE1PBE, PBEh1PBE, and PW91TPSS gave 10 hits. Among these, M05, B971, mPW2PLYPD, Wb97x, and PW91TPSS were among the “recommended” list of at least one dimer from each class. Long‐range correction (LC) of Hirao and coworkers to exchange‐correlation functionals showed massive improvement in geometry and E_int. The best performing LC‐functionals were LC‐G96KCIS and LC‐PKZBPKZB. Our results predict that M06L is the most trustworthy DFT method in Gaussian09 to study small non‐covalently interacting systems. © 2013 Wiley Periodicals, Inc.     

Study of Cl−(H2O)n (n = 1–4) using basin‐hopping method coupled with density functional theory

Cl^?(H₂O)_n (n = 1–4) clusters were investigated using a basin‐hopping (BH) algorithm coupled with density functional theory (DFT). Structures, energetics, thermodynamics, vertical detachment energies, and vibrational frequencies were obtained from high‐level ab initio calculations. Through comparisons with previous theoretical and experimental data, it was demonstrated that the combination of the BH method and DFT could accurately predict the global and local minima of Cl^?(H₂O)_n (n = 1–4). Additionally, to optimize larger Cl^?(H₂O)_n (n > 4) clusters, several popular density functionals as well as DF‐LMP2 (Schütz et al., J. Chem. Phys. 2004, 121, 737) (second‐order Møller‐Plesset perturbation theory using local and density fitting approximations) were tested with appropriate basis sets through comparisons with MP2 optimized results. DF‐LMP2 will be used in future studies because its overall performance in describing the relative binding energies and the geometrical parameters of Cl^?(H₂O)_n (n = 1–4) was outstanding in this study. © 2013 Wiley Periodicals, Inc.     

Simple modification of the Lee?Yang?Parr correlation functional to satisfy exact nonuniform scaling requirements

In density functional theory (DFT) the exchange-correlation energy has to be approximated. One of the most widely used approximations to the correlation component, the Lee(SINGLEBOND)Yang(SINGLEBOND)Parr (LYP) functional does not obey nonuniform scaling requirements. We propose a modifying function that makes this functional satisfy nonuniform scaling constraints. As a result, two modified LYP functionals are suggested. One of them, optimized in a preliminary way, exhibits some properties better than the original LYP functional. This implies that nonuniform scaling requirements convey very important information and all approximate functionals in DFT should be made to satisfy these requirements. Our modifying factor could be used to improve the nonuniform scaling properties not only of E_c ^LYP [n] but also of most functionals which do not obey these conditions. © 1996 John Wiley & Sons, Inc.     

DFT and MO calculations of atomic and molecular chemisorption energies on surface cluster models

Summary Density functional theory (DFT) (including gradient corrections) and MCPF calculations have been performed for atomic (H, C, N, O) and molecular CH _x (x = 1–3) chemisorption on cluster models of different sites of the Cu(100) surface. The DFT and MCPF results are in good agreement once the important effects of core-valence correlation have been accounted for in the MCPF calculations by including contributions from a core polarization potential (CPP); in the DFT approach the core-valence correlation is obtained directly from the total density using the functional. Very large effects on the four-fold hollow site binding energy from core-valence correlation are found for C, N and CH. Several different DFT functionals were employed and compared in the calculations.     

Comparison of density functionals for the study of the high spin low spin gap in Fe(III) spin crossover complexes

A detailed investigation of the accuracy of different quantum mechanical methods for the study of iron(III) spin crossover complexes is presented. The energy spin state gap between the high and low spin states; ΔE (HS‐LS) of nine iron(III) quinolylsalicylaldiminate complexes were calculated with nine different DFT functionals, then compared. DFT functionals: B3LYP, B3LYP‐D3, B3LYP*, BH&HLYP, BP86, OLYP, OPBE, M06L, and TPSSh were tested with six basis sets: 3‐21G*, dgdzvp, 6‐31G**, cc‐pVDZ, Def2TZVP, and cc‐pVTZ. The cations from the X‐ray crystal structures of [Fe(qsal‐OMe)₂]Cl·MeCN·H₂O, [Fe(qsal‐OMe)₂]Cl·2MeOH·0.5H₂O, [Fe(qsal‐OMe)₂]BF₄·MeOH, [Fe(qsal‐OMe)₂]NCS·CH₂Cl₂, [Fe(qsal‐F)₂]NCS, [Fe(qsal‐Cl)₂]NCS·MeOH, [Fe(qsal‐Br)₂]NCS·MeOH, [Fe(qsal‐I)₂]OTf·MeOH, and [Fe(qsal)₂]NCS?CH₂Cl₂ were used as starting structures. The results show that B3LYP, B3LYP‐D3, OLYP, and OPBE with a 6‐31G**, Def2TZVP, and cc‐pVTZ basis set give reasonable results of ΔE (HS‐LS) compared with the experimental data. The enthalpy of [Fe(qsal‐I)₂]⁺ calculated with an OLYP functional and cc‐pVTZ basis set (1.48 kcal/mol) most closely matches the experimental data (1.34 kcal/mol). B3LYP* yields an enthalpy of 5.92 kcal/mol suggesting it may be unsuitable for these Fe(III) complexes, mirroring recent results by Kepp (Inorg . Chem ., 2016, 55 , 2717–2727).     

Estimation of isotropic nuclear magnetic shieldings in the CCSD(T) and MP2 complete basis set limit using affordable correlation calculations

A linear correlation between isotropic nuclear magnetic shielding constants for seven model molecules (CH₂O, H₂O, HF, F₂, HCN, SiH₄ and H₂S) calculated with 37 methods (34 density functionals, RHF, MP2 and CCSD(T)), with affordable pcS‐2 basis set and corresponding complete basis set results, estimated from calculations with the family of polarization‐consistent pcS‐n basis sets is reported. This dependence was also supported by inspection of profiles of deviation between CBS estimated nuclear shieldings and shieldings obtained with the significantly smaller basis sets pcS‐2 and aug‐cc‐pVTZ‐J for the selected set of 37 calculation methods. It was possible to formulate a practical approach of estimating the values of isotropic nuclear magnetic shielding constants at the CCSD(T)/CBS and MP2/CBS levels from affordable CCSD(T)/pcS‐2, MP2/pcS‐2 and DFT/CBS calculations with pcS‐n basis sets. The proposed method leads to a fairly accurate estimation of nuclear magnetic shieldings and considerable saving of computational efforts. Copyright © 2013 John Wiley & Sons, Ltd.     

Synthesis,experimental and in silico studies of N‐fluorenylmethoxycarbonyl‐O‐tert‐butyl‐N‐methyltyrosine,coupled with CSD data: a survey of interactions in the crystal structures of Fmoc–amino acids

Recently, fluorenylmethoxycarbonyl (Fmoc) amino acids (e.g. Fmoc–tyrosine or Fmoc–phenylalanine) have attracted growing interest in biomedical research and industry, with special emphasis directed towards the design and development of novel effective hydrogelators, biomaterials or therapeutics. With this in mind, a systematic knowledge of the structural and supramolecular features in recognition of those properties is essential. This work is the first comprehensive summary of noncovalent interactions combined with a library of supramolecular synthon patterns in all crystal structures of amino acids with the Fmoc moiety reported so far. Moreover, a new Fmoc‐protected amino acid, namely, 2‐{[(9H‐fluoren‐9‐ylmethoxy)carbonyl](methyl)amino}‐3‐{4‐[(2‐hydroxypropan‐2‐yl)oxy]phenyl}propanoic acid or N‐fluorenylmethoxycarbonyl‐O‐tert‐butyl‐N‐methyltyrosine, Fmoc‐N‐Me‐Tyr(t‐Bu)‐OH, C₂₉H₃₁NO₅, was successfully synthesized and the structure of its unsolvated form was determined by single‐crystal X‐ray diffraction. The structural, conformational and energy landscape was investigated in detail by combined experimental and in silico approaches, and further compared to N‐Fmoc‐phenylalanine [Draper et al. (2015). CrystEngComm, 42 , 8047–8057]. Geometries were optimized by the density functional theory (DFT) method either in vacuo or in solutio. The polarizable conductor calculation model was exploited for the evaluation of the hydration effect. Hirshfeld surface analysis revealed that H…H, C…H/H…C and O…H/H…O interactions constitute the major contributions to the total Hirshfeld surface area in all the investigated systems. The molecular electrostatic potentials mapped over the surfaces identified the electrostatic complementarities in the crystal packing. The prediction of weak hydrogen‐bonded patterns via Full Interaction Maps was computed. Supramolecular motifs formed via C—H…O, C—H…π, (fluorenyl)C—H…Cl(I), C—Br…π(fluorenyl) and C—I…π(fluorenyl) interactions are observed. Basic synthons, in combination with the Long‐Range Synthon Aufbau Modules, further supported by energy‐framework calculations, are discussed. Furthermore, the relevance of Fmoc‐based supramolecular hydrogen‐bonding patterns in biocomplexes are emphasized, for the first time.     

Potential energy surface of the reaction of imidazole with peroxynitrite: Density functional theory study

This article presents a theoretical investigation of the reaction mechanism of imidazole nitration by peroxynitrite using density functional theory calculations. Understanding this reaction mechanism will help in elucidating the mechanism of guanine nitration by peroxynitrite, which is one of the assumed chemical pathways for damaging DNA in cells. This work focuses on the analysis of the potential energy surface (PES) for this reaction in the gas phase. Calculations were carried out using Hartree–Fock (HF) and density functional theory (DFT) Hamiltonians with double‐zeta basis sets ranging from 6‐31G(d) to 6‐31++G(d,p), and the triple‐zeta basis set 6‐311G(d). The computational results reveal that the reaction of imidazole with peroxynitrite in gas phase produces the following species: (i) hydroxide ion and 2‐nitroimidazole, (ii) hydrogen superoxide ion and 2‐nitrosoimidazole, and (iii) water and 2‐nitroimidazolide. The rate‐determining step is the formation of a short‐lived intermediate in which the imidazole C₂ carbon is covalently bonded to peroxynitrite nitrogen. Three short‐lived intermediates were found in the reaction path. These intermediates are involved in a proton‐hopping transport from C₂ carbon to the terminal oxygen of the ? O? O moiety of peroxynitrite via the nitroso (ON? ) oxygen. Both HF and DFT calculations (using the Becke3–Lee–Yang–Parr functional) lead to similar reaction paths for proton transport, but the landscape details of the PES for HF and DFT calculations differ. This investigation shows that the reaction of imidazole with peroxynitrite produces essentially the same types of products (nitro‐ and nitroso‐) as observed experimentally in the reaction of guanine with peroxynitrite, which makes the former reaction a good model to study by computation the essential characteristics of the latter reaction. Nevertheless, the computationally determined activation energy for imidazole nitration by peroxynitrite in the gas phase is 84.1 kcal/mol (calculated at the B3LYP/6‐31++G(d,p) level), too large for an enzymatic reaction. Exploratory calculations on imidazole nitration in solution, and on the reaction of 9‐methylguanine with peroxynitrite in the gas phase and solution, show that solvation increases the activation energy for both imidazole and guanine, and that the modest decrease (15 kcal mol^?1) in the activation energy, due to the adjacent six member ring of guanine, is counterbalanced by solvation. These results lead to the speculation that proton tunneling may be at the origin of experimentally observed high reaction rate of guanine nitration by peroxynitrite in solution. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005