Diffusion monte carlo study of the hydrogen molecules adsorbed on C4H3Li

Diffusion Monte Carlo (DMC) simulations were used to calculate the binding energies for hydrogen molecules adsorbed on the lithium metal–organic complex C₄H₃Li. The calculations use all‐electron DMC techniques where every electron is explicitly included in the simulation. Also we have systematically studied it using density functional theory (DFT) methods, revealing that each C₄H₃Li can hold up to four H₂ molecules and the adsorption distance is about 2.2 Å. The DMC binding energies are in the range of 0.055–0.143 eV and are compared with those obtained with DFT using various exchange‐correlation functionals, with values ranging from 0.029 to 0.504 eV. These results indicate that caution is required applying DFT methods to weakly bound systems such as hydrogen storage materials based on lithium‐doped metal–organic frameworks. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

The Nature of the Binding of Au, Ag, and Pd to Benzene, Coronene, and Graphene: From Benchmark CCSD(T) Calculations to Plane-Wave DFT Calculations

The adsorption of Ag, Au, and Pd atoms on benzene, coronene, and graphene has been studied using post Hartree-Fock wave function theory (CCSD(T), MP2) and density functional theory (M06-2X, DFT-D3, PBE, vdW-DF) methods. The CCSD(T) benchmark binding energies for benzene-M (M = Pd, Au, Ag) complexes are 19.7, 4.2, and 2.3 kcal/mol, respectively. We found that the nature of binding of the three metals is different: While silver binds predominantly through dispersion interactions, the binding of palladium has a covalent character, and the binding of gold involves a subtle combination of charge transfer and dispersion interactions as well as relativistic effects. We demonstrate that the CCSD(T) benchmark binding energies for benzene-M complexes can be reproduced in plane-wave density functional theory calculations by including a fraction of the exact exchange and a nonempirical van der Waals correction (EE+vdW). Applying the EE+vdW method, we obtained binding energies for the graphene-M (M = Pd, Au, Ag) complexes of 17.4, 5.6, and 4.3 kcal/mol, respectively. The trends in binding energies found for the benzene-M complexes correspond to those in coronene and graphene complexes. DFT methods that use empirical corrections to account for the effects of vdW interactions significantly overestimate binding energies in some of the studied systems.     

Enhancement of hydrogen physisorption on graphene and carbon nanotubes by Li doping

Density-functional calculations of the adsorption of molecular hydrogen on a planar graphene layer and on the external surface of a (4,4) carbon nanotube, undoped and doped with lithium, have been carried out. Hydrogen molecules are physisorbed on pure graphene and on the nanotube with binding energies about 80-90 meV/molecule. However, the binding energies increase to 160-180 meV/molecule for many adsorption configurations of the molecule near a Li atom in the doped systems. A charge-density analysis shows that the origin of the increase in binding energy is the electronic charge transfer from the Li atom to graphene and the nanotube. The results support and explain qualitatively the enhancement of the hydrogen storage capacity observed in some experiments of hydrogen adsorption on carbon nanotubes doped with alkali atoms.     

Assessing the accuracy of quantum Monte Carlo and density functional theory for energetics of small water clusters

We present a detailed study of the energetics of water clusters (H(2)O)(n) with n ≤ 6, comparing diffusion Monte Carlo (DMC) and approximate density functional theory (DFT) with well converged coupled-cluster benchmarks. We use the many-body decomposition of the total energy to classify the errors of DMC and DFT into 1-body, 2-body and beyond-2-body components. Using both equilibrium cluster configurations and thermal ensembles of configurations, we find DMC to be uniformly much more accurate than DFT, partly because some of the approximate functionals give poor 1-body distortion energies. Even when these are corrected, DFT remains considerably less accurate than DMC. When both 1- and 2-body errors of DFT are corrected, some functionals compete in accuracy with DMC; however, other functionals remain worse, showing that they suffer from significant beyond-2-body errors. Combining the evidence presented here with the recently demonstrated high accuracy of DMC for ice structures, we suggest how DMC can now be used to provide benchmarks for larger clusters and for bulk liquid water.     

A benchmark quantum Monte Carlo study of the ground state chromium dimer

We report calculations of the ground state energy and binding curve of the chromium dimer using the variational and diffusion quantum Monte Carlo (VMC and DMC) methods. We examined various single‐determinant and multideterminant wavefunctions multiplied by a Jastrow factor as a trial/guiding wavefunction for VMC/DMC. The molecular orbitals in the single determinants were calculated using restricted or unrestricted Hartree–Fock or density functional theory (DFT) calculations where five commonly used local (SVWN5), semilocal (PW91 and BLYP), and hybrid (B1LYP and B3LYP) functionals were examined. The multideterminant expansions were obtained from the generalized valence bond and (truncated) unrestricted configuration interaction with single and double excitations (UCISD) methods. We also examined a UCISD wavefunction in which UCISD expansions were added to the UB3LYP single‐determinant reference, and their coefficients were optimized at the VMC level. In addition to the wavefunction dependence, the effects of pseudopotentials and backflow transformation were also investigated. The UB3LYP single‐determinant and multideterminant wavefunctions were found to give the variationally best DMC energies within the framework of single‐determinant and multideterminants, respectively, though both the DMC energies were higher than twice the DMC atomic energy. Some of the VMC binding curves show a flat or quite shallow well bottom, which gets recovered deeper by DMC. All the DMC binding curves have a minimum indicating a bound state, but the unrestricted ones overestimate the equilibrium bond length. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

On the accuracy of density-functional theory exchange-correlation functionals for H bonds in small water clusters: benchmarks approaching the complete basis set limit

The ability of several density-functional theory (DFT) exchange-correlation functionals to describe hydrogen bonds in small water clusters (dimer to pentamer) in their global minimum energy structures is evaluated with reference to second order Moller-Plesset perturbation theory (MP2). Errors from basis set incompleteness have been minimized in both the MP2 reference data and the DFT calculations, thus enabling a consistent systematic evaluation of the true performance of the tested functionals. Among all the functionals considered, the hybrid X3LYP and PBE0 functionals offer the best performance and among the nonhybrid generalized gradient approximation functionals, mPWLYP and PBE1W perform best. The popular BLYP and B3LYP functionals consistently underbind and PBE and PW91 display rather variable performance with cluster size.     

Noncovalent interactions between cisplatin and graphene prototypes

Cisplatin (CP) has been widely used as an anticancer drug for more than 30 years despite severe side effects due to its low bioavailability and poor specificity. For this reason, it is paramount to study and design novel nanomaterials to be used as vectors capable to effectively deliver the drug to the biological target. The CP square‐planar geometry, together with its low water solubility, suggests that it could be possibly easily adsorbed on 2D graphene nanostructures through the interaction with the related highly conjugated π‐electron system. In this work, pyrene has been first selected as the minimum approximation to the graphene plane, which allows to properly study the noncovalent interactions determining the CP adsorption. In particular, electronic structure calculations at the MP2C and DFT‐SAPT levels of theory have allowed to obtain benchmark interaction energies for some limiting configurations of the CP–pyrene complex, as well as to assess the role of the different contributions to the total interaction: it has been found that the parallel configurations of the aggregate are mainly stabilized around the minimum region by dispersion, in a similar way as for complexes bonded through π‐π interactions. Then, the benchmark interaction energies have been used to test corresponding estimations obtained within the less expensive DFT to validate an optimal exchange‐correlation functional which includes corrections to take properly into account for the dispersion contribution. Reliable DFT interaction energies have been therefore obtained for CP adsorbed on graphene prototypes of increasing size, ranging from coronene, ovalene, and up to C₁₅₀H₃₀. Finally, DFT geometry optimizations and frequency calculations have also allowed a reliable estimation of the adsorption enthalpy of CP on graphene, which is found particularly favorable (about −20 kcal/mol at 298 K and 1 bar) being twice that estimated for the corresponding benzene adsorption. © 2017 Wiley Periodicals, Inc.     

Quantum Monte Carlo study of the reaction: Cl+CH3OH-->CH2OH+HCl

A theoretical study is reported of the Cl+CH3OH-->CH2OH+HCl reaction based on the diffusion Monte Carlo (DMC) variant of the quantum Monte Carlo method. Using a DMC trial function constructed as a product of Hartree-Fock and correlation functions, we have computed the barrier height, heat of reaction, atomization energies, and heats of formation of reagents and products. The DMC heat of reaction, atomization energies, and heats of formation are found to agree with experiment to within the error bounds of computation and experiment. M?ller-Plesset second order perturbation theory (MP2) and density functional theory, the latter in the B3LYP generalized gradient approximation, are found to overestimate the experimental heat of reaction. Intrinsic reaction coordinate calculations at the MP2 level of theory demonstrate that the reaction is predominantly direct, i.e., proceeds without formation of intermediates, which is consistent with a recent molecular beam experiment. The reaction barrier as determined from MP2 calculations is found to be 2.24 kcal/mol and by DMC it is computed to be 2.39(49) kcal/mol.     

Hydrogen bonding in the urea dimers and adenine–thymine DNA base pair: anharmonic effects in the intermolecular H-bond and intramolecular H-stretching vibrations

The equilibrium structures, binding energies, vibrational harmonic frequencies, and the anharmonic corrections for two different (cyclic and asymmetric) urea dimers and for the adenine–thymine DNA base pair system have been studied using the second-order Møller–Plesset perturbation theory (MP2) method and different density functional theory (DFT) exchange–correlation (XC) functionals (BLYP, B3LYP, PBE, HCTH407, KMLYP, and BH and HLYP) with the D95V, D95V**, and D95V++** basis sets. The widely used a posteriori Boys–Bernardi or counterpoise correction scheme for basis set superposition error (BSSE) has been included in the calculations to take into account the BSSE effects during geometry optimization (on structure), on binding energies and on the different levels of approximation used for calculating the vibrational frequencies. The results obtained with the ab initio MP2 method are compared with those calculated with different DFT XC functionals; and finally the suitability of these DFT XC functionals to describe intermolecular hydrogen bonds as well as harmonic frequencies and the anharmonic corrections is assessed and discussed.     

Halogen‐Bonding Interactions with π Systems: CCSD(T), MP2, and DFT Calculations

Halogen bonding is a noncovalent interaction between a halogen atom and a nucleophilic site. Interactions involving the π electrons of aromatic rings have received, up to now, little attention, despite the large number of systems in which they are present. We report binding energies of the interaction between either NCX or PhX (X=F, Cl, Br, I) and the aromatic benzene system as determined with the coupled cluster with perturbative triple excitations method [CCSD(T)] extrapolated at the complete basis set limit. Results are compared with those obtained by Møller–Plesset perturbation theory to second order (MP2) and density functional theory (DFT) calculations by using some of the most common functionals. Results show the important role of DFT in studying this interaction.     

Computational study of hydrogen binding by metal-organic framework-5

We report the results of quantum chemistry calculations on H(2) binding by the metal-organic framework-5 (MOF)-5. Density functional theory calculations were used to calculate the atomic positions, lattice constant, and effective atomic charges from the electrostatic potential for the MOF-5 crystal structure. Second-order M?ller-Plesset perturbation theory was used to calculate the binding energy of H(2) to benzene and H(2)-1,4-benzenedicarboxylate-H(2). To achieve the necessary accuracy, the large Dunning basis sets aug-cc-pVTZ, and aug-cc-pVQZ were used, and the results were extrapolated to the basis set limit. The binding energy results were 4.77 kJ/mol for benzene, 5.27 kJ/mol for H(2)-1,4-benzenedicarboxylate-H(2). We also estimate binding of 5.38 kJ/mol for Li-1,4-benzenedicarboxylate-Li and 6.86 kJ/mol at the zinc oxide corners using second-order M?ller-Plesset perturbation theory. In order to compare our theoretical calculations to the experimental hydrogen storage results, grand canonical Monte Carlo calculations were performed. The Monte Carlo simulations identify a high energy binding site at the corners that quickly saturated with 1.27 H(2) molecules at 78 K. At 300 K, a broad range of binding sites are observed.     

Barrier heights of hydrogen‐transfer reactions with diffusion quantum monte carlo method

Hydrogen‐transfer reactions are an important class of reactions in many chemical and biological processes. Barrier heights of H‐transfer reactions are underestimated significantly by popular exchange–correlation functional with density functional theory (DFT), while coupled‐cluster (CC) method is quite expensive and can be applied only to rather small systems. Quantum Monte‐Carlo method can usually provide reliable results for large systems. Performance of fixed‐node diffusion quantum Monte‐Carlo method (FN‐DMC) on barrier heights of the 19 H‐transfer reactions in the HTBH38/08 database is investigated in this study with the trial wavefunctions of the single‐Slater–Jastrow form and orbitals from DFT using local density approximation. Our results show that barrier heights of these reactions can be calculated rather accurately using FN‐DMC and the mean absolute error is 1.0 kcal/mol in all‐electron calculations. Introduction of pseudopotentials (PP) in FN‐DMC calculations improves efficiency pronouncedly. According to our results, error of the employed PPs is smaller than that of the present CCSD(T) and FN‐DMC calculations. FN‐DMC using PPs can thus be applied to investigate H‐transfer reactions involving larger molecules reliably. In addition, bond dissociation energies of the involved molecules using FN‐DMC are in excellent agreement with reference values and they are even better than results of the employed CCSD(T) calculations using the aug‐cc‐pVQZ basis set. © 2017 Wiley Periodicals, Inc.     

Verification of DFT-predicted hydrogen storage capacity of VC3H3 complex using molecular dynamics simulations

Density functional theory (DFT) and Fourth‐order Möller–Plesset (MP4) perturbation theory calculations are performed to examine the possibility of hydrogen storage in V‐capped VC₃H₃ complex. Stability of bare and H₂ molecules adsorbed V‐capped VC₃H₃ complex is verified using DFT and MP4 method. Thermo‐chemistry calculations are carried out to estimate the Gibbs free corrected averaged H₂ adsorption energy which reveals whether H₂ adsorption on V‐capped VC₃H₃ complex is energetically favorable, at different temperatures. We use different exchange and correlation functionals employed in DFT to see their effect on H₂ adsorption energy. Molecular dynamic (MD) simulations are performed to confirm whether this complex adsorbs H₂ molecules at a finite temperature. We elucidate the correlation between H₂ adsorption energy obtained from density functional calculations and retaining number of H₂ molecules on VC₃H₃ complex during MDs simulations at various temperatures. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

A benchmark quantum chemical study of the stacking interaction between larger polycondensed aromatic hydrocarbons

Large-scale electronic structure calculations were performed for the interaction energy between coronene, C₂₄H₁₂ with circumcoronene, C₅₄H₁₈, and between two circumcoronene molecules, in order to get a picture of the interaction between larger graphene sheets. Most calculations were performed at the SCS-MP2 level but we have corrected them for higher-order correlation effects using a calculation on the coronene-circumcoronene system at the quadratic CI, QCISD(T) level. Our best estimate for the interaction energy between coronene and circumcoronene is 32.1?kcal/mol. We estimate the binding of coronene on a graphite surface to be 37.4 or 1.56?kcal/mol per carbon atom (67.5?meV/C atom). This is also our estimate for the exfoliation energy of graphite. It is higher than most previous theoretical estimates. The SCS-MP2 method which reproduces the CCSD(T) and QCISD(T) values very well for smaller aromatic hydrocarbons, e.g., for the benzene dimer, increasingly overestimates dispersion as the bandgap (the HOMO-LUMO separation) decreases. The barrier to the sliding motion of coronene on circumcoronene is 0.45?kcal/mol, and for two circumcoronene molecules 1.85?kcal/mol (0.018 and 0.034?kcal/mol per C atom, respectively). This means that larger graphenes cannot easily glide over each other.     

The relative energies of polypeptide conformers predicted by linear scaling second-order Møller-Plesset perturbation theory

We describe an implementation of the cluster-in-molecule(CIM) resolution of the identity(RI) approximation second-order M?ller–Plesset perturbation theory(CIM-RI-MP2), with the purpose of extending RI-MP2 calculations to very large systems. For typical conformers of several large polypeptides, we calculated their conformational energy differences with the CIM-RI-MP2 and the generalized energy-based fragmentation MP2(GEBF-MP2) methods, and compared these results with the density functional theory(DFT) results obtained with several popular functionals. Our calculations show that the conformational energy differences obtained with CIM-RI-MP2 and GEBF-MP2 are very close to each other. In comparison with the GEBF-MP2 and CIM-RI-MP2 relative energies, we found that the DFT functionals(CAM-B3LYP-D3, LC-?PBE-D3, M05-2X, M06-2X and ?B97XD) can give quite accurate conformational energy differences for structurally similar conformers, but provide less-accurate results for structurally very different conformers.     

Cooperative enhancement of water binding to crownophane by multiple hydrogen bonds: analysis by high level ab initio calculations.

The intermolecular interaction energy of the model system of the water-crownophane complex was analyzed. The water molecule has four hydrogen bonds, with the two hydrogen-donating phenolic hydroxy groups and two hydrogen-accepting oxygen atoms of the poly-oxyethylene chain of the crownophane in the complex. The MP2/6-311G(2d,2p) level calculations of the model system of the complex (hydrogen donating unit + hydrogen accepting unit + water) indicate that the binding energy of the water is 21.85 kcal/mol and that the hydrogen bond cooperativity increases the binding energy as much as 3.67 kcal/mol. The calculated interaction energies depend on the basis set, while the basis set dependence of the cooperative increment is negligible. Most of the cooperative increment is covered by the HF level calculation, which suggests that the major source of the hydrogen bond cooperativity in this system has its origin in induction. The BLYP/6-311G** and PW91/6-311G** level interaction energies of the model system are close to the MP2/6-311G** interaction energies, which suggests that the DFT calculations with these functionals are useful methods to evaluated the interactions of hydrogen bonded systems.     

含羧酸镁(钙)官能团的多孔芳香骨架材料储氢性能的预测

用MP2方法,TZVPP基组以及基组重叠误差(BSSE)校正计算了氢分子与修饰在多孔芳香骨架(PAF)上的羧酸镁、羧酸钙官能团的相互作用,并建立了描述这一相互作用的分子力学力场.在此基础上用巨正则系综蒙特卡洛(GCMC)模拟预测了氢气在该种新型PAF材料上的吸附等温线.量子化学计算结果表明,每个羧酸镁、羧酸钙官能团分别可以提供13、14个氢分子吸附位点,与每个氢分子的平均结合能在8kJ·mol-1左右.通过比较不同温度和压力下材料的绝对吸附量和超额吸附量发现,在PAF骨架中引入羧酸镁、羧酸钙官能团可以显著提高材料的综合储氢性能,达到并超过了美国能源部提出的2015年储氢标准.同时该工作还揭示了氢吸附量与材料的表面积、空腔体积和分子作用强度间的复杂关系.     

van der Waals corrections to density functional theory calculations: Methane,ethane, ethylene,benzene, formaldehyde,ammonia, water,PBE, and CPMD

Parameters are developed for a practical application of the empirical van der Waals (vdW) correction infrastructure available in the CPMD density functional theory (DFT) code. The binding energy, geometry, and potential energy surface (PES) are examined for methane, ethane, ethylene, formaldehyde, ammonia, three benzene dimer geometries, and three benzene–water geometries. The vdW corrected results compare favorably with MP2 and CCSD(T) calculations near the complete basis set limits, and with experimental results where they are available.