On the Binding Strength Sequence for Nucleic Acid Bases and C(60) with Density Functional and Dispersion-corrected Density Functional Theories: Whether C(60) could protect nucleic acid bases from radiation-induced damage?

The major objective of this paper is to address a controversial binding sequence between nucleic acid bases (NABs) and C(60) by investigating adsorptions of NABs and their cations on C(60) fullerene with a variety of density functional theories including two novel hybrid meta-GGA functionals, M05-2x and M06-2x, as well as a dispersion-corrected density functional, PBE-D. The M05-2x/6-311++G** provides the same binding sequence as previously reported, guanine(G) > cytosine(C) > adenine (A) > thymine (T); however, M06-2x switches the binding strengths of A and C, and PBE-D eventually results in the following sequence, G>A>T>C, which is the same as the widely accepted hierarchy for the stacking of NABs on other carbon nanomaterials such as single-walled carbon nanotube and graphite. The results indicate that the questionable relative binding strength is due to insufficient electron correlation treatment with the M05-2x or even the M06-2x method. The binding energy of G@C(60) obtained with the M06-2x/6-311++G(d,p) and the PBE-D/cc-pVDZ is -7.10 and -8.07 kcal/mol, respectively, and the latter is only slightly weaker than that predicted by the MP2/6-31G(d,p) (-8.10kca/mol). Thus, the PDE-D performs better than the M06-2x for the observed NAB@C(60) π-stacked complexes. To discuss whether C(60) could prevent NABs from radiation-induced damage, ionization potentials of NABs and C(60), and frontier molecular orbitals of the complexes NABs@C(60) and (NABs@C(60))(+) are also extensively investigated. These results revealed that when an electron escapes from the complexes, a hole was preferentially created in C(60) for T and C complexes, while for G and A the hole delocalizes over the entire complex, rather than a localization on the C(60) moiety. The interesting finding might open a new strategy for protecting DNA from radiation-induced damage and offer a new idea for designing C(60)-based antiradiation drugs.     

Theoretical Evidence for the Stronger Ability of Thymine to Disperse SWCNT than Cytosine and Adenine: self-stacking of DNA bases vs their cross-stacking with SWCNT

Self-stacking of four DNA bases, adenine (A), cytosine (C), guanine (G) and thymine (T), and their cross-stacking with (5,5) as well as (10,0) single walled carbon nanotubes (SWCNTs) were extensively investigated with a novel hybrid DFT method, MPWB1K/cc-pVDZ. The binding energies were further corrected with MP2/6-311++G(d,p) method in both gas phase and aqueous solution, where the solvent effects were included with conductor-like polarized continuum model (CPCM) model and UAHF radii. The strongest self-stacking of G and A takes displaced anti-parallel configuration, but un-displaced or "eclipsed" anti-parallel configuration is the most stable for C and T. In gas phase the self-stacking of nucleobases decreases in the sequence G>A>C>T, while because of quite different solvent effects their self-stacking in aqueous solution exhibits a distinct sequence A>G>T>C. For a given base, cross-stacking is stronger than self-stacking in both gas phase and aqueous solution. Binding energy for cross-stacking in gas phase varies as G>A>T>C for both (10,0) and (5,5) SWCNTs, and the binding of four nucleobases to (10,0) is slightly stronger than to (5,5) SWCNT by a range of 0.1-0.5 kcal/mol. The cross-stacking in aqueous solution varies differently from that gas phase: A>G>T>C for (10,0) SWCNT and G>A>T>C for (5,5) SWCNT. It is suggested that the ability of nucleobases to disperse SWCNT depends on relative strength [Formula: see text] of self-stacking and cross-stacking with SWCNT in aqueous solution. Of the four investigated nucleobases thymine (T) exhibits the highest [Formula: see text] which can well explain the experimental finding that T more efficiently functionalizes SWCNT than C and A.     

Pairwise additivity of energy components in protein-ligand binding: the HIV II protease-Indinavir case

An energy expansion (binding energy decomposition into n-body interaction terms for n ≥ 2) to express the receptor-ligand binding energy for the fragmented HIV II protease-Indinavir system is described to address the role of cooperativity in ligand binding. The outcome of this energy expansion is compared to the total receptor-ligand binding energy at the Hartree-Fock, density functional theory, and semiempirical levels of theory. We find that the sum of the pairwise interaction energies approximates the total binding energy to ～82% for HF and to >95% for both the M06-L density functional and PM6-DH2 semiempirical method. The contribution of the three-body interactions amounts to 18.7%, 3.8%, and 1.4% for HF, M06-L, and PM6-DH2, respectively. We find that the expansion can be safely truncated after n=3. That is, the contribution of the interactions involving more than three parties to the total binding energy of Indinavir to the HIV II protease receptor is negligible. Overall, we find that the two-body terms represent a good approximation to the total binding energy of the system, which points to pairwise additivity in the present case. This basic principle of pairwise additivity is utilized in fragment-based drug design approaches and our results support its continued use. The present results can also aid in the validation of non-bonded terms contained within common force fields and in the correction of systematic errors in physics-based score functions.     

Density-functional approaches to noncovalent interactions: a comparison of dispersion corrections (DFT-D), exchange-hole dipole moment (XDM) theory, and specialized functionals

A systematic study of techniques for treating noncovalent interactions within the computationally efficient density functional theory (DFT) framework is presented through comparison to benchmark-quality evaluations of binding strength compiled for molecular complexes of diverse size and nature. In particular, the efficacy of functionals deliberately crafted to encompass long-range forces, a posteriori DFT+dispersion corrections (DFT-D2 and DFT-D3), and exchange-hole dipole moment (XDM) theory is assessed against a large collection (469 energy points) of reference interaction energies at the CCSD(T) level of theory extrapolated to the estimated complete basis set limit. The established S22 [revised in J. Chem. Phys. 132, 144104 (2010)] and JSCH test sets of minimum-energy structures, as well as collections of dispersion-bound (NBC10) and hydrogen-bonded (HBC6) dissociation curves and a pairwise decomposition of a protein-ligand reaction site (HSG), comprise the chemical systems for this work. From evaluations of accuracy, consistency, and efficiency for PBE-D, BP86-D, B97-D, PBE0-D, B3LYP-D, B970-D, M05-2X, M06-2X, ωB97X-D, B2PLYP-D, XYG3, and B3LYP-XDM methodologies, it is concluded that distinct, often contrasting, groups of these elicit the best performance within the accessible double-ζ or robust triple-ζ basis set regimes and among hydrogen-bonded or dispersion-dominated complexes. For overall results, M05-2X, B97-D3, and B970-D2 yield superior values in conjunction with aug-cc-pVDZ, for a mean absolute deviation of 0.41 - 0.49 kcal/mol, and B3LYP-D3, B97-D3, ωB97X-D, and B2PLYP-D3 dominate with aug-cc-pVTZ, affording, together with XYG3/6-311+G(3df,2p), a mean absolute deviation of 0.33 - 0.38 kcal/mol.     

Blind prediction of SAMPL4 cucurbit[7]uril binding affinities with the mining minima method

Accurate methods for predicting protein–ligand binding affinities are of central interest to computer-aided drug design for hit identification and lead optimization. Here, we used the mining minima (M2) method to predict cucurbit[7]uril binding affinities from the SAMPL4 blind prediction challenge. We tested two different energy models, an empirical classical force field, CHARMm with VCharge charges, and the Poisson–Boltzmann surface area solvation model; and a semiempirical quantum mechanical (QM) Hamiltonian, PM6-DH+, coupled with the COSMO solvation model and a surface area term for nonpolar solvation free energy. Binding affinities based on the classical force field correlated strongly with the experiments with a correlation coefficient (R²) of 0.74. On the other hand, binding affinities based on the QM energy model correlated poorly with experiments (R² = 0.24), due largely to two major outliers. As we used extensive conformational search methods, these results point to possible inaccuracies in the PM6-DH+ energy model or the COSMO solvation model. Furthermore, the different binding free energy components, solute energy, solvation free energy, and configurational entropy showed significant deviations between the classical M2 and quantum M2 calculations. Comparison of different classical M2 free energy components to experiments show that the change in the total energy, i.e. the solute energy plus the solvation free energy, is the key driving force for binding, with a reasonable correlation to experiment (R² = 0.56); however, accounting for configurational entropy further improves the correlation.     

Transferable scoring function based on semiempirical quantum mechanical PM6-DH2 method: CDK2 with 15 structurally diverse inhibitors

A semiempirical quantum mechanical PM6-DH2 method accurately covering the dispersion interaction and H-bonding was used to score fifteen structurally diverse CDK2 inhibitors. The geometries of all the complexes were taken from the X-ray structures and were reoptimised by the PM6-DH2 method in continuum water. The total scoring function was constructed as an estimate of the binding free energy, i.e., as a sum of the interaction enthalpy, interaction entropy and the corrections for the inhibitor desolvation and deformation energies. The applied scoring function contains a clear thermodynamical terms and does not involve any adjustable empirical parameter. The best correlations with the experimental inhibition constants (ln K _i) were found for bare interaction enthalpy (r ² = 0.87) and interaction enthalpy corrected for ligand desolvation and deformation energies (r ² = 0.77); when the entropic term was considered, however, the correlation becomes worse but still acceptable (r ² = 0.52). The resulting correlation based on the PM6-DH2 scoring function is better than previously published function based on various docking/scoring, SAR studies or advanced QM/MM approach, however, the robustness is limited by number of available experimental data used in the correlation. Since a very similar correlation between the experimental and theoretical results was found also for a different system of the HIV-1 protease, the suggested scoring function based on the PM6-DH2 method seems to be applicable in drug design, even if diverse protein–ligand complexes have to be ranked.     

Electronic and vibrational linear and nonlinear polarizabilities of Li@C60 and [Li@C60]+

Electronic and vibrational nuclear relaxation (NR) contributions to the dipole (hyper)polarizabilities of the endohedral fullerene Li@C(60) and its monovalent cation [Li@C(60)](+) are calculated at the (U)B3LYP level. Many results are new, while others differ significantly from those reported previously using more approximate methods. The properties are compared with those of the corresponding hypothetical noninteracting systems with a valence electron transferred from Li to the cage. Whereas the NR contribution to the static linear polarizabilities is small in comparison with the corresponding electronic property, the opposite is true for the static hyperpolarizabilities. A relatively small, but non-negligible, NR contribution to the dc-Pockels effect is obtained in the infinite frequency approximation.     

Binding free energies in the SAMPL6 octa-acid host–guest challenge calculated with MM and QM methods

We have estimated free energies for the binding of eight carboxylate ligands to two variants of the octa-acid deep-cavity host in the SAMPL6 blind-test challenge (with or without endo methyl groups on the four upper-rim benzoate groups, OAM and OAH, respectively). We employed free-energy perturbation (FEP) for relative binding energies at the molecular mechanics (MM) and the combined quantum mechanical (QM) and MM (QM/MM) levels, the latter obtained with the reference-potential approach with QM/MM sampling for the MM → QM/MM FEP. The semiempirical QM method PM6-DH+ was employed for the ligand in the latter calculations. Moreover, binding free energies were also estimated from QM/MM optimised structures, combined with COSMO-RS estimates of the solvation energy and thermostatistical corrections from MM frequencies. They were performed at the PM6-DH+ level of theory with the full host and guest molecule in the QM system (and also four water molecules in the geometry optimisations) for 10–20 snapshots from molecular dynamics simulations of the complex. Finally, the structure with the lowest free energy was recalculated using the dispersion-corrected density-functional theory method TPSS-D3, for both the structure and the energy. The two FEP approaches gave similar results (PM6-DH+/MM slightly better for OAM), which were among the five submissions with the best performance in the challenge and gave the best results without any fit to data from the SAMPL5 challenge, with mean absolute deviations (MAD) of 2.4–5.2 kJ/mol and a correlation coefficient (R²) of 0.77–0.93. This is the first time QM/MM approaches give binding free energies that are competitive to those obtained with MM for the octa-acid host. The QM/MM-optimised structures gave somewhat worse performance (MAD?=?3–8 kJ/mol and R²?=?0.1–0.9), but the results were improved compared to previous studies of this system with similar methods.     

Energy relationship and the polarization of C_(60) cage in the endohedral complexes (X@C_(60))

In this paper, we carry out the calculation on the system (X@C60)(X=Li, Na, K, Kb, Cs; F, Cl, Br, I), where the position of X changes along 5 typical symmetry directions. For the calculation of quantum chemistry we use EHMO/ASED method, for the calculation of molecular mechanics we use Buckingham potential (exp-6-1) function, and for the calculation of thermo-chemical cycle we use individually isolating the processes such as the structure variation, charge transfer and charge distribution, and their interactions etc. The calculation results show that (1) In the region of radius r≈0.2 nm of the Ceo cage, the potential field is nearly spherical; (2) Except for Li and Na, the systems are the most stable with minimum energies at the center of C60 cage. For Li and Na, the systems are the most stable with minimum energies at r≈0.16 nm and r≈0.13 nm, respectively. In view of the interactive region of chemical bonds, the interactions between X and the C60 cage do not belong to the classical chemical bonds; (     

The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals

We present two new hybrid meta exchange- correlation functionals, called M06 and M06-2X. The M06 functional is parametrized including both transition metals and nonmetals, whereas the M06-2X functional is a high-nonlocality functional with double the amount of nonlocal exchange (2X), and it is parametrized only for nonmetals.The functionals, along with the previously published M06-L local functional and the M06-HF full-Hartree–Fock functionals, constitute the M06 suite of complementary functionals. We assess these four functionals by comparing their performance to that of 12 other functionals and Hartree–Fock theory for 403 energetic data in 29 diverse databases, including ten databases for thermochemistry, four databases for kinetics, eight databases for noncovalent interactions, three databases for transition metal bonding, one database for metal atom excitation energies, and three databases for molecular excitation energies. We also illustrate the performance of these 17 methods for three databases containing 40 bond lengths and for databases containing 38 vibrational frequencies and 15 vibrational zero point energies. We recommend the M06-2X functional for applications involving main-group thermochemistry, kinetics, noncovalent interactions, and electronic excitation energies to valence and Rydberg states. We recommend the M06 functional for application in organometallic and inorganometallic chemistry and for noncovalent interactions. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Contribution to the Mark S. Gordon 65th Birthday Festschrift Issue.     

Benchmark Study of the Performance of Density Functional Theory for Bond Activations with (Ni,Pd)-Based Transition-Metal Catalysts

The performance of 23 density functionals, including one LDA, four GGAs, three meta-GGAs, three hybrid GGAs, eight hybrid meta-GGAs, and ten double-hybrid functionals, was investigated for the computation of activation energies of various covalent main-group single bonds by four catalysts: Pd, PdCl⁻, PdCl₂, and Ni (all in the singlet state). A reactant complex, the barrier, and reaction energy were considered, leading to 164 energy data points for statistical analysis. Extended Gaussian AO basis sets were used in all calculations. The best functional for the complete benchmark set relative to estimated CCSD(T)/CBS reference data is PBE0-D3, with an MAD value of 1.1 kcal mol⁻¹ followed by PW6B95-D3, the double hybrid PWPB95-D3, and B3LYP-D3 (1.9 kcal mol⁻¹ each). The other tested hybrid meta-GGAs perform less well (M06-HF: 7.0 kcal mol⁻¹; M06-2X: 6.3 kcal mol⁻¹; M06: 4.9 kcal mol⁻¹) for the investigated reactions. In the Ni case, some double hybrids show larger errors due to partial breakdown of the perturbative treatment for the correlation energy in cases with difficult electronic structures (partial multi-reference character). Only double hybrids either with very low amounts of perturbative correlation (e.g., PBE0-DH) or that use the opposite-spin correlation component only (e.g., PWPB95) seem to be more robust. We also investigated the effect of the D3 dispersion correction. While the barriers are not affected by this correction, significant and mostly positive results were observed for reaction energies. Furthermore, six very recently proposed double-hybrid functionals were analyzed regarding the influence of the amount of Fock exchange as well as the type of perturbative correlation treatment. According to these results, double hybrids with <50–60 % of exact exchange and ∼30 % perturbative correlation perform best.     

Computational studies on ethylene addition to nickel bis(dithiolene)

The density functionals B3LYP, B3PW91, BMK, HSE06, LC-ωPBE, M05, M06, O3LYP, TPSS, ω-B97X, and ω-B97XD are used to optimize key transition states and intermediates for ethylene addition to Ni(edt)(2) (edt = S(2)C(2)H(2)). The efficacy of the basis sets 6-31G**, 6-31++G**, cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ is also examined. The geometric parameters optimized with different basis sets and density functionals are similar and agree well with experimental values. The ω-B97XD functional gives relative energies closest to those from CCSD, while M06 and HSE06 yield results close to those from CCSD(T). CASSCF and CASSCF-PT2 calculation results are also given. Variation of the relative energies from different density functionals appears to arise, in part, from the multireference character of this system, as confirmed by the T1 diagnostic and CASSCF calculations.     

Accurate binding energies of hydrogen,halogen, and dihydrogen bonded complexes and cation enhanced binding strengths

Interaction energy (E_int) values of a variety of hydrogen, halogen, and dihydrogen bonded complexes in the weak, medium, and strong regimes have been computed using W1BD, MP2, M06L density functional theory, and hybrid methods MP4//MP2, MP4//M06L, and CCSD(T)//MP2. W1BD level E_int and CCSD(T) results reported in the literature show very good agreement (mean absolute deviation = 0.19 kcal/mol). MP2 underestimates E_int while M06L shows accurate behavior for all except halogen and charge‐assisted hydrogen bonds. MP4//MP2, MP4//M06L, and CCSD(T)//MP2 yield E_int very close to those obtained from W1BD. The high accuracy energy data at MP4/MP2 is used to study the effect of a cation (Li⁺, NH₄⁺) on the E_int. The cation enhances electron donation from the donor to noncovalent bonding region leading to substantial enhancement in E_int (～141–566% for Li⁺ and ～105–539% for NH₄⁺) and promotes a noncovalent bond in the weak regime to medium regime and that in the medium regime to strong regime. © 2014 Wiley Periodicals, Inc.     

Description of noncovalent interactions involving π-system with high precision: An assessment of RPA,MP2, and DFT-D methods

Efficient approaches with high precision are essential for understanding the formation and stability of noncovalent interaction complexes. Here, 21 noncovalent interaction complexes involving π-system are selected and grouped in three subsets according to ETS–NOCV method: dispersion-dominated, electrostatic-dominated, and mixed. We mainly focus on examining the performance of random-phase approximation (RPA) on these π systems. The tested RPA-based method includes standard RPA and its variants including the related single excitations (SEs), renormalized single excitations (rSEs), second-order screened exchange (SOSEX), and the renormalized second-order perturbation theory (rPT2). The routine second-order Møller–Plesset perturbation theory (MP2) and three popular DFT-D functionals (M06-2X-D3, ωB97XD, and PBE-D3(BJ)) are also assessed for comparison. In this work, besides the calculation of interaction energies at Dunning-type aug-cc-pVDZ and aug-cc-pVTZ basis set, we also present a larger database of interaction energies calculated using MP2 and RPA methods with Dunning-type aug-cc-pVQZ basis set. An accurate CCSD(T)/CBS scheme is used to provide benchmark database. In addition to the high-level results, we also provide potential energy surfaces (PES) of different interaction type. Among all the tested methods, MP2 has a satisfactory performance on electrostatic-dominated and mixed-type systems, except for dispersion-dominated systems. DFT-D functionals, especially ωB97XD functional, has a balanced performance across all the tested systems. Importantly, for RPA-based methods, the calculation accuracy can be dramatically improved by taking into account SE or exchange effects, especially in the mixed complexes. We conclude that rPT2 among all the test RPA-based methods gives an overall satisfactory performance across different interaction types. © 2019 Wiley Periodicals, Inc.     

Benchmark database on isolated small peptides containing an aromatic side chain: comparison between wave function and density functional theory methods and empirical force field

A detailed quantum chemical study on five peptides (WG, WGG, FGG, GGF and GFA) containing the residues phenylalanyl (F), glycyl (G), tryptophyl (W) and alanyl (A) -- where F and W are of aromatic character -- is presented. When investigating isolated small peptides, the dispersion interaction is the dominant attractive force in the peptide backbone-aromatic side chain intramolecular interaction. Consequently, an accurate theoretical study of these systems requires the use of a methodology covering properly the London dispersion forces. For this reason we have assessed the performance of the MP2, SCS-MP2, MP3, TPSS-D, PBE-D, M06-2X, BH&H, TPSS, B3LYP, tight-binding DFT-D methods and ff99 empirical force field compared to CCSD(T)/complete basis set (CBS) limit benchmark data. All the DFT techniques with a '-D' symbol have been augmented by empirical dispersion energy while the M06-2X functional was parameterized to cover the London dispersion energy. For the systems here studied we have concluded that the use of the ff99 force field is not recommended mainly due to problems concerning the assignment of reliable atomic charges. Tight-binding DFT-D is efficient as a screening tool providing reliable geometries. Among the DFT functionals, the M06-2X and TPSS-D show the best performance what is explained by the fact that both procedures cover the dispersion energy. The B3LYP and TPSS functionals-not covering this energy-fail systematically. Both, electronic energies and geometries obtained by means of the wave-function theory methods compare satisfactorily with the CCSD(T)/CBS benchmark data.     

Synthesis,crystal structure,and BSA interaction with a new Co(II) complex based on 5-(trifluoromethyl)pyridine-2-carboxylic acid

Abstract

A new complex, Co(Htpc)₂(H₂O)₂ (1) (Htpc = 5-(trifluoromethyl)pyridine-2-carboxylic acid), has been synthesized and characterized by X-ray single-crystal diffraction, elemental analysis, infrared spectral analysis, and thermogravimetric analysis. Meanwhile, the optimized geometric structure of the ligand was determined using the M06-2X functional of density functional theory (DFT) with the 6-311?+?G(d, p) basis set. The gap energies ΔE between the frontier molecular orbitals were computed in different solvent media (water, methanol and ethanol) using the time dependent density functional theory (TD-DFT)/M06-2X by applying the Polarizable Continuum Model (PCM). The coordination sphere around Co(II) is distorted octahedral with two chelating tpc^- ligands and two coordinated water molecules. Bovine serum albumin (BSA) binding properties of the ligand, CoCl₂·6H₂O and 1 were investigated by fluorescence and UV–Vis absorption spectroscopy, revealing 1 exhibits higher binding affinity with BSA than free ligand and CoCl₂·6H₂O. ΔG, ΔH and ΔS at 298 and 308?K manifested that van der Waals interactions and hydrogen bonds were the main forces in the binding process.     

Theoretical study of hydrogen bonding in homodimers and heterodimers of amide, boronic acid, and carboxylic acid, free and in encapsulation complexes

The homodimers and the heterodimers of two amides, two boronic acids, and two carboxylic acids have been calculated in the gas phase and in N,N-dimethylformamide (DMF) and CCl(4) solvents using the DFT (M06-2X and M06-L) and the MP2 methods in conjunction with the 6-31G(d,p) and 6-311+G(d,p) basis sets. Furthermore, their pairwise coencapsulation was studied to examine its effect on the calculated properties of the hydrogen bonds at the ONIOM[M06-2X/6-31G(d,p);PM6], ONIOM[MP2/6-31G(d,p); PM6], and M06-2X/6-31G(d,p) levels of theory. The present work is directed toward the theoretical rationalization and interpretation of recent experimental results on hydrogen bonding in encaptulation complexes [D. Ajami et al. J. Am. Chem. Soc. 2011, 133, 9689-9691]. The calculated dimerization energy (ΔE) values range from 0.74 to 0.35 eV for the different dimers in the gas phase, with the ordering carboxylic homodimers > amide-carboxylic dimers > amide homodimers > boronic-carboxylic dimers > amide-boronic dimers > boronic homodimers. In solvents, generally smaller ΔE values are calculated with only small variations in the ordering. In the capsule, the ΔE values range between 0.67 and 0.33 eV with practically the same ordering as in the gas phase. The calculated % distributions of the encapsulated dimers, taking into account statistical factors, are in agreement with the experimental distribution, where the occurrence of boronic homodimer dominates, even though it is calculated to have the smallest ΔE.     

Theoretical study of M(PH3)2 complexes of C60, corannulene (C20H10), and sumanene (C21H12) (M = Pd or Pt). Unexpectedly large binding energy of M(PH3)2(C60)

DFT and MP2 to MP4(SDQ) methods were applied to M(PH3)2(C60), Pt(PH3)2(C20H10), and Pt(PH3)2(C21H12) (M = Pd or Pt, C20H10 = corannulene, and C21H12 = sumanene). The binding energy considerably fluctuates around MP2 and MP3 levels but much less upon going from MP3 to MP4(SDQ) in Pt(PH3)2(C2H4), Pt(PH3)2(C20H10), and Pt(PH3)2(C21H12). Also, the MP4(SDQ) method presents a binding energy similar to that of the CCSD(T) method in Pt(PH3)2(C2H4). Thus, it is likely that the MP4(SDQ) method is useful to evaluate binding energies of these complexes. The binding energies of Pt(PH3)2(C20H10) and Pt(PH3)2(C21H12) are evaluated to be 24.9 and 26.1 kcal/mol, respectively, by the MP4(SDQ) method and only +5.8 and -2.6 kcal/mol, respectively, by the DFT(B3LYP) method. These MP4(SDQ)-calculated binding energies of Pt(PH3)2(C20H10) and Pt(PH3)2(C21H12) are similar to that of Pt(PH3)2(C2H4), which strongly suggests that these complexes can be successfully synthesized. The binding energy of Pt(PH3)2(C60) is evaluated to be 44.8 and 45.5 kcal/mol with the ONIOM(MP4(SDQ):UFF) and ONIOM(MP4(SDQ):B3LYP) methods, respectively, and that of the Pd analogue is evaluated to be 39.9 kcal/mol with the ONIOM(MP4(SDQ):UFF) method, whereas the DFT(B3LYP), DFT(BVP86), and DFT(BPW91) methods provide much smaller binding energies. It is noted that these binding energies are much larger than those of the ethylene, corannulene, and sumanene analogues. This difference is reasonably interpreted in terms that the LUMO of C60 is at much lower energy than those of ethylene, corannulene, and sumanene. We investigated also how to separate the high level and the low level regions in the ONIOM calculation of M(PH3)2(C60) and proposed here the reasonable way to evaluate the binding energy of transition-metal complexes of C60.     

Lowest-energy structures of (C60)nX (X=Li+,Na+,K+,Cl-) and (C60)nYCl (Y=Li,Na,K) clusters for n</=13

Basin-hopping global optimization is used to find likely candidates for the lowest minima on the potential energy surface of (C(60))(n)X (X=Li(+),Na(+),K(+),Cl(-)) and (C(60))(n)YCl (Y=Li,Na,K) clusters with n相似文献