Modeling van der Waals interactions between proteins and inorganic surfaces from time-dependent density functional theory calculations

In this work we show how the ab initio determination of van der Waals coefficients within time-dependent density functional theory can be used to build efficient and accurate atomistic models that describe the long-range interactions of proteins with other proteins and of proteins with semi-conducting surfaces. The model parameters are fitted so that they reproduce the ab initio van der Waals coefficients of amino acids and dipeptides. We then assess the quality of our results by comparing ab initio van der Waals coefficients for larger peptides with the coefficients yielded by the models. The different sets of parameters can be easily incorporated in current empirical force field methods, thus providing an essential ingredient for molecular dynamics simulations of proteins close to surfaces.     

Ab initio molecular dynamics of hydrogen dissociation on metal surfaces using neural networks and novelty sampling

We outline a hybrid multiscale approach for the construction of ab initio potential energy surfaces (PESs) useful for performing six-dimensional (6D) classical or quantum mechanical molecular dynamics (MD) simulations of diatomic molecules reacting at single crystal surfaces. The algorithm implements concepts from the corrugation reduction procedure, which reduces energetic variation in the PES, and uses neural networks for interpolation of smoothed ab initio data. A novelty sampling scheme is implemented and used to identify configurations that are most likely to be predicted inaccurately by the neural network. This hybrid multiscale approach, which couples PES construction at the electronic structure level to MD simulations at the atomistic scale, reduces the number of density functional theory (DFT) calculations needed to specify an accurate PES. Due to the iterative nature of the novelty sampling algorithm, it is possible to obtain a quantitative measure of the convergence of the PES with respect to the number of ab initio calculations used to train the neural network. We demonstrate the algorithm by first applying it to two analytic potentials, which model the H2/Pt(111) and H2/Cu(111) systems. These potentials are of the corrugated London-Eyring-Polanyi-Sato form, which are based on DFT calculations, but are not globally accurate. After demonstrating the convergence of the PES using these simple potentials, we use DFT calculations directly and obtain converged semiclassical trajectories for the H2/Pt(111) system at the PW91/generalized gradient approximation level. We obtain a converged PES for a 6D hydrogen-surface dissociation reaction using novelty sampling coupled directly to DFT. These results, in excellent agreement with experiments and previous theoretical work, are compared to previous simulations in order to explore the sensitivity of the PES (and therefore MD) to the choice of exchange and correlation functional. Despite having a lower energetic corrugation in our PES, we obtain a broader reaction probability curve than previous simulations, which is attributed to increased geometric corrugation in the PES and the effect of nonparallel dissociation pathways.     

Theoretical investigation of the He-I2(E3Πg) ion-pair state: ab initio intermolecular potential and vibrational levels

We present a theoretical study on the potential energy surface and vibrational bound states of the E electronic excited state of the HeI(2) van der Waals system. The interaction energies are computed using accurate ab initio methods and large basis sets. Relativistic small-core effective core potentials in conjunction with a quintuple-zeta quality basis set are employed for the heavy iodine atoms in multireference configuration interaction calculations for the (3)A' and (3)A" states. For the representation of the potential energy surface we used a general interpolation technique for constructing potential surfaces from ab initio data based on the reproducing kernel Hilbert space method. The surface presents global and local minima for T-shaped configurations with well-depths of 33.2 and 4.6 cm(-1), respectively. Vibrational energies and states are computed through variational quantum mechanical calculations. We found that the binding energy of the HeI(2)(E) T-shaped isomer is 16.85 cm(-1), in excellent agreement with recent experimental measurements. In lieu of more experimental data we also report our predictions on higher vibrational levels and we analyze the influence of the underlying surface on them. This is the first attempt to represent the potential surface of such a highly excited electronic state of a van der Waals complex, and it demonstrates the capability of the ab initio technology to provide accurate results for carrying out reliable studies to model experimental data.     

Sensitivity of the H + CH3 → CH4 recombination rate constant to the shape of the CH stretching potential

Using a potential-energy surface obtained in part from ab initio calculations, the H + CH₃ → CH₄ bimolecular rate constant at T = 300 K is determined from a Monte Carlo classical trajectory study. Representing the CH stretching potential with a standard Morse function instead ofthe ab initio curve increases the calculated rate constant by an order of magnitude. The experimental recombination rate constant is intermediate of the rate constants calculated with the Morse and ab initio stretching potentials.Two properties of the H + CH₃ α CH₄ potential-energy surface which significantly affect the recombination rate constant are the shape of the CH stretching potential and the attenuation of the H₃CH bending frequencies. Ab initio calculations with a hierarchy of basis sets and treatment of electron correlation indicate the latter is properly described [13]. The exact shape of the CH stretching potential is not delineated by the ab initio calculations, since the ab initio calculations are not converged for bond lengths of 2.0–3.0 Å [12]. However, the form of this stretching potential deduced from the highest-level ab initio calculations, and fit analytically by eq. (2), is significantly different from a Morse function. The experimental recombination rate constant is intermediate of the rate constants calculated with the Morse and ab initio CH stretching potentials. This indicates that the actual CH potential energy curve lies between the Morse and ab initio curves. This is consistent with the finding that potential energy curves for diatomics are not well described by a Morse function [12].     

Atom-atom partitioning of total (super)molecular energy: the hidden terms of classical force fields

Classical force fields describe the interaction between atoms that are bonded or nonbonded via simple potential energy expressions. Their parameters are often determined by fitting to ab initio energies and electrostatic potentials. A direct quantum chemical guide to constructing a force field would be the atom-atom partitioning of the energy of molecules and van der Waals complexes relevant to the force field. The authors used the theory of quantum chemical topology to partition the energy of five systems [H2, CO, H2O, (H2O)2, and (HF)2] in terms of kinetic, Coulomb, and exchange intra-atomic and interatomic contributions. The authors monitored the variation of these contributions with changing bond length or angle. Current force fields focus only on interatomic interaction energies and assume that these purely potential energy terms are the only ones that govern structure and dynamics in atomistic simulations. Here the authors highlight the importance of self-energy terms (kinetic and intra-atomic Coulomb and exchange).     

Investigation of hydrogen and methane adsorption/separation on silicon nanotubes: a hierarchical multiscale method from quantum mechanics to molecular simulation

A combination of ab initio quantum mechanical (QM) calculations and canonical Monte Carlo (CMC) simulations are employed to investigate possible usage of single-walled silicon nanotubes (SWSiNTs) as a novel media for hydrogen and methane adsorption as well as their separation from each other. By fitting the force field, a Morse potential model is selected as an efficient potential to describe the binding energies between both hydrogen-SiNTs and methane-SiNTs obtained from ab initio calculations. Then CMC simulations are performed to evaluate the adsorption and separation behaviors of H₂ and CH₄ on the three different sizes of SiNTs including (5, 5), (7, 7), and (9, 9) SiNTs at ambient temperatures and pressures from 1 up to 10 MPa. As a comparison, the adsorption and separation of H₂ and CH₄ on the (8, 8) CNTs which are isodiameter with (5, 5) SiNTs are also simulated. Results are indicative of remarkable enhancement of H₂ and CH₄ adsorption capacity on the SiNTs compared to the CNTs, which arise from stronger van der Waals (VDW) attractions. In the case of methane adsorption on SiNTs, the stored volumetric energy exceeds the goal of the US Freedom CAR Partnership by 2010, which can not be achieved by methane compression at such low pressures. Moreover, simulation results indicate that SiNTs preferentially adsorb methane relative to hydrogen in their equimolar mixture, which results in efficient separation of these gases from each other at 293 K.     

Van der waals interactions and decrease of the rotational barrier of methyl-sized rotators: a theoretical study

Rotational barriers of methyl-sized molecular rotators are investigated theoretically using ab initio and empirical force field calculations in molecular models simulating various environmental conditions experienced by the molecular rotors. Calculations on neopentane surrounded by methyl groups suggest that the neopentane's methyl rotational potential energy barrier can be reduced by up to an order of magnitude by locating satellite functional groups around the rotator at a geometry that destabilizes the staggered conformation of the rotator through van der Waals repulsive interactions and reduces the staggered/eclipsed relative energy difference. Molecular mechanics and molecular dynamics calculations indicate that this barrier-reducing geometry can also be found in molecular rotators surface mounted on graphite surfaces or carbon nanotube models. In these models, molecular dynamics simulations show that the rotation of methyl-sized functional groups can be catalyzed by van der Waals interactions, thus making very rigid rotators become thermally activated at room temperature. These results are discussed in the context of design of nanostructures and use of methyl groups as markers for microenvironmental conditions.     

Accurate description of van der Waals complexes by density functional theory including empirical corrections

An empirical method to account for van der Waals interactions in practical calculations with the density functional theory (termed DFT-D) is tested for a wide variety of molecular complexes. As in previous schemes, the dispersive energy is described by damped interatomic potentials of the form C6R(-6). The use of pure, gradient-corrected density functionals (BLYP and PBE), together with the resolution-of-the-identity (RI) approximation for the Coulomb operator, allows very efficient computations for large systems. Opposed to previous work, extended AO basis sets of polarized TZV or QZV quality are employed, which reduces the basis set superposition error to a negligible extend. By using a global scaling factor for the atomic C6 coefficients, the functional dependence of the results could be strongly reduced. The "double counting" of correlation effects for strongly bound complexes is found to be insignificant if steep damping functions are employed. The method is applied to a total of 29 complexes of atoms and small molecules (Ne, CH4, NH3, H2O, CH3F, N2, F2, formic acid, ethene, and ethine) with each other and with benzene, to benzene, naphthalene, pyrene, and coronene dimers, the naphthalene trimer, coronene. H2O and four H-bonded and stacked DNA base pairs (AT and GC). In almost all cases, very good agreement with reliable theoretical or experimental results for binding energies and intermolecular distances is obtained. For stacked aromatic systems and the important base pairs, the DFT-D-BLYP model seems to be even superior to standard MP2 treatments that systematically overbind. The good results obtained suggest the approach as a practical tool to describe the properties of many important van der Waals systems in chemistry. Furthermore, the DFT-D data may either be used to calibrate much simpler (e.g., force-field) potentials or the optimized structures can be used as input for more accurate ab initio calculations of the interaction energies.     

A new "spectroscopic" potential energy surface for formaldehyde in its ground electronic state

We report a new "spectroscopic" potential energy surface (PES) of formaldehyde (H(2)(12)C(16)O) in its ground electronic state, obtained by refining an ab initio PES in a least-squares fitting to the experimental spectroscopic data for formaldehyde currently available in the literature. The ab initio PES was computed using the CCSD(T)/aug-cc-pVQZ method at 30 840 geometries that cover the energy range up to 44 000 cm(-1) above equilibrium. Ro-vibrational energies of formaldehyde were determined variationally for this ab initio PES by means of the program TROVE [Theoretical ROtation-Vibration Energies; S. N. Yurchenko, W. Thiel, and P. Jensen, J. Mol. Spectrosc. 245, 126 (2007)]. The parameter values in the analytical representation of the PES were optimized in fittings to 319 ro-vibrational energies with J = 0, 1, 2, and 5. The initial parameter values in the fittings were those of the ab initio PES, the ro-vibrational eigenfunctions obtained from this PES served as a basis set during the fitting process, and constraints were imposed to ensure that the refined PES does not deviate unphysically from the ab initio one in regions of configuration space not sampled by the experimental data. The resulting refined PES, referred to as H(2)CO-2011, reproduces the available experimental J ≤ 5 data with a root-mean-square error of 0.04 cm(-1).     

An ab initio study of van der Waals potential energy parameters for silver clusters

We employ ab initio calculations of van der Waals complexes to study the potential energy parameters (C(6) coefficients) of van der Waals interactions for modeling of the adsorption of silver clusters on the graphite surface. Electronic structure calculations of the (Ag(2))(2), Ag(2)-H(2), and Ag(2)-C(6)H(6) complexes are performed using a coupled-cluster approach that includes single, double, and perturbative triple excitations (CCSD(T)), M?ller-Plesset second-order perturbation theory (MP2), and spin-component-scaled MP2 (SCS-MP2) methods. Using the atom pair approximation, the C(6) coefficients for silver-silver, silver-hydrogen, and silver-carbon atom systems are obtained after subtracting the energies of quadrupole-quadrupole interactions from the total electronic energy.     

Sensitivity analysis of thermodynamic properties of liquid water: a general approach to improve empirical potentials

A sensitivity analysis of bulk water thermodynamics is presented in an effort to understand the relation between qualitative features of molecular potentials and properties that they predict. The analysis is incorporated in molecular dynamics simulations and investigates the sensitivity of the Helmholtz free energy, internal energy, entropy, heat capacity, pressure, thermal pressure coefficient, and static dielectric constant to components of the potential rather than the parameters of a given functional form. The sensitivities of the properties are calculated with respect to the van der Waals repulsive and the attractive parts, plus short- and long-range Coulomb parts of three four site empirical water potentials: TIP4P, Dang-Chang and TTM2R. The polarization sensitivity is calculated for the polarizable Dang-Chang and TTM2R potentials. This new type of analysis allows direct comparisons of the sensitivities for different potentials that use different functional forms. The analysis indicates that all investigated properties are most sensitive to the van der Waals repulsive, the short-range Coulomb and the polarization components of the potentials. When polarization is included in the potentials, the magnitude of the sensitivity of the Helmholtz free energy, internal energy, and entropy with respect to this part of the potential is comparable in magnitude to the other electrostatic components. In addition similarities in trends of observed sensitivities for nonpolarizable and polarizable potentials lead to the conclusion that the complexity of the model is not of critical importance for the calculation of these thermodynamic properties for bulk water. The van der Waals attractive and the long-range Coulomb sensitivities are relatively small for the entropy, heat capacity, thermal pressure coefficient and the static dielectric constant, while small changes in any of the potential contributions will significantly affect the pressure. The analysis suggests a procedure for modification of the potentials to improve predictions of thermodynamic properties and we demonstrate this general approach for modifying potentials for one of the potentials.     

Identification of the active sites and mechanism for partial methane oxidation to methanol over copper-exchanged CHA zeolites

To make methane a suitable energy carrier and transport less costly, it is an urgent and challenging task for us to convert methane to liquid under mild conditions efficiently. In this study, we explored partial methane oxidation to methanol by density functional theory(DFT) calculations using a hybrid functional(HSE06) with van der Waals(vdW) interactions. The stabilities of different active sites over SSZ-13 and SAPO-34, two CHA type zeolites, are thoroughly investigated by ab initio molecular dynamics(AIMD) simulations and ab initio thermodynamics analyses. Four possible active sites, namely [CuOHCu]~(2+), [Cu(OH)_2Cu]~(2+),[CuOCu]~(2+) and [CuOH]~+, are identified stable. Methane-to-methanol reaction mechanisms are further studied upon these most stable active sites, among which [CuOCu]~(2+) and [CuOH]~+ are proved to be reactive. The migration of species among zeolite pores are also discussed, which accounts for the activity on [CuOH]~+ sites. This concept may represent a more complete picture of catalytic reactions over zeolites in general.     

An ab initio investigation of the O(3P)-H2(1sigma(g)+) van der Waals well

We report an ab initio study of the van der Waals region of the O(3P)-H2 potential energy surface based on RCCSD(T) calculations with an aug-cc-pVQZ basis supplemented by bond functions. In addition, an open-shell implementation of symmetry-adapted perturbation theory (SAPT) is used to corroborate the RCCSD(T) calculations and to investigate the relative magnitudes of the various contributions to the van der Waals interaction. We also investigate the effect of the spin-orbit coupling on the position and depth of the van der Waals well. We predict the van der Waals minimum to occur in perpendicular geometry, and located at a closer distance than a secondary well in colinear geometry. The potentials obtained in the present study confirm the previous calculations of Alexander [M. H. Alexander, J. Chem. Phys., 1998, 108, 4467], but disagree with the earlier work of Harding and co-workers [Z. Li, V. A. Apkarian and L. B. Harding, J. Chem. Phys., 1997, 106, 942] as well as with recently refitted surfaces of Brand?o and coworkers [J. Brand?o, C. Mogo and B. C. Silva, J. Chem. Phys., 2004, 121, 8861]. Inclusion of spin-orbit coupling reduces the depth of the van der Waals minimum without causing a change in its position.     

Interpolating moving least-squares methods for fitting potential energy surfaces: an application to the H2CN unimolecular reaction

Classical trajectories have been used to compute rates for the unimolecular reaction H2CN-->H+HCN on a fitted ab initio potential energy surface (PES). The ab initio energies were obtained from CCSD(T)/aug-cc-pvtz electronic structure calculations. The ab initio energies were fitted by the interpolating moving least-squares (IMLS) method. This work continues the development of the IMLS method for producing ab initio PESs for use in molecular dynamics simulations of many-atom systems. A dual-level scheme was used in which the preliminary selection of data points was done using a low-level theory and the points used for fitting the final PES were obtained at the desired higher level of theory. Classical trajectories were used on various low-level IMLS fits to tune the fit to the unimolecular reaction under study. Procedures for efficiently picking data points, selecting basis functions, and defining cutoff limits to exclude distant points were investigated. The accuracy of the fitted PES was assessed by comparing interpolated values of quantities to the corresponding ab initio values. With as little as 330 ab initio points classical trajectory rate constants were converged to 5%-10% and the rms error over the six-dimensional region sampled by the trajectories was a few tenths of a kcal/mol.     

A modification to the COSMIC parameterisation using ab initio constrained potential functions

Summary The H.H non-bonded potential employed in the current COSMIC force field has been contrasted with H.H potentials used in a number of other force fields. Initial conversion of the variety of functions to a Morse format, achieved using a simple graphical fitting procedure, allowed a direct comparison to be made, showing the COSMIC potential to differ considerably from the other potentials. This difference was reflected in the failure of COSMIC to reproduce ab initio and experimental energies for molecules with significant H.H interactions, with particular reference to the energy curves of benzophenone and diphenyl ether. Considerable improvement in these energies is produced by the use of a Morse function originally based on the H.H potential used in MM3.     

An ab initio spin-orbit-corrected potential energy surface and dynamics for the F + CH4 and F + CHD3 reactions

We report an analytical ab initio three degrees of freedom (3D) spin-orbit-correction surface for the entrance channel of the F + methane reaction obtained by fitting the differences between the spin-orbit (SO) and non-relativistic electronic ground state energies computed at the MRCI+Q/aug-cc-pVTZ level of theory. The 3D model surface is given in terms of the distance, R(C-F), and relative orientation, Euler angles ? and θ, of the reactants treating CH(4) as a rigid rotor. The full-dimensional (12D) "hybrid" SO-corrected potential energy surface (PES) is obtained from the 3D SO-correction surface and a 12D non-SO PES. The SO interaction has a significant effect in the entrance-channel van der Waals region, whereas the effect on the energy at the early saddle point is only ～5% of that at the reactant asymptote; thus, the SO correction increases the barrier height by ～122 cm(-1). The 12D quasiclassical trajectory calculations for the F + CH(4) and F + CHD(3) reactions show that the SO effects decrease the cross sections by a factor of 2-4 at low collision energies and the effects are less significant as the collision energy increases. The inclusion of the SO correction in the PES does not change the product state distributions.     

Microwave and ab initio studies of rare gas-methane van der Waals complexes

Rotational spectra of the weakly bound Kr-methane van der Waals complex were recorded using a pulsed molecular beam Fourier transform microwave spectrometer in the range from 3.5 to 18 GHz. Spectra of 25 isotopomers of Kr-methane were assigned and analyzed. For isotopomers containing CH4, 13CH4, and CD4, two sets of transitions with K = 0 and one with K = 1 were recorded, correlating to the j = 0, 1, and 2 rotational levels of free methane, respectively (j is the rotational angular momentum quantum number of the methane monomer). For isotopomers containing CH3D and CHD3, two K = 0 components were recorded, correlating to the j(k) = 0(0) and 1(1) rotational levels of free methane (k corresponds to the projection of j onto the C3 axis of CH3D and CHD3). The obtained spectroscopic results were used to derive van der Waals bond distance R, van der Waals stretching frequency nu(s), and the corresponding stretching force constant k(s). Nuclear spin statistical weights of individual states were obtained from molecular symmetry group analyses and were compared with the observed relative transition intensities. The tentatively assigned j = 2 transitions were more intense than predicted from symmetry considerations. This is attributed to a relatively large effective dipole moment of this state, supported by ab initio dipole moment calculations. Ab initio potential energy calculations of Kr-CH4 and Ar-CH4 were done at the coupled cluster level of theory, with single and double excitations and perturbative inclusion of triple excitations, using the aug-cc-pVTZ basis set supplemented with bond functions. The theoretical results show that the angular dynamics of the dimer does not change significantly when the binding partner of methane changes from Ar to Kr. The dipole moment of Ar-CH4 was calculated at various configurations, providing a qualitative explanation for the unsuccessful spectral searches for rotational transitions of Ar-CH4.     

Ab initio ground and excited state potential energy surfaces for NO-Kr complex and dynamics of Kr solids with NO impurity

The intermolecular potentials for the NO(X 2Pi)-Kr and NO(A 2Sigma+)-Kr systems have been calculated using highly accurate ab initio calculations. The spin-restricted coupled cluster method for the ground 1 2A' state [NO(X 2Pi)-Kr] and the multireference singles and doubles configuration interaction method for the excited 2 2A' state [NO(A 2Sigma+)-Kr], respectively, were used. The potential energy surfaces (PESs) show two linear wells and one that is almost in the perpendicular position. An analytical representation of the PESs has been constructed for the triatomic systems and used to carry out molecular dynamics (MD) simulations of the NO-doped krypton matrix response after excitation of NO. MD results are shown comparatively for three sets of potentials: (1) anisotropic ab initio potentials [NO molecule direction fixed during the dynamics and considered as a point (its center of mass)], (2) isotropic ab initio potentials (isotropic part in a Legendre polynomial expansion of the PESs), and (3) fitted Kr-NO potentials to the spectroscopic data. An important finding of this work is that the anisotropic and isotropic ab initio potentials calculated for the Kr-NO triatomic system are not suitable for describing the dynamics of structural relaxation upon Rydberg excitation of a NO impurity in the crystal. However, the isotropic ab initio potential in the ground state almost overlaps the published experimental potential, being almost independent of the angle asymmetry. This fact is also manifested in the radial distribution function around NO. However, in the case of the excited state the isotropic ab initio potential differs from the fitted potentials, which indicates that the Kr-NO interaction in the matrix is quite different because of the presence of the surrounding Kr atoms acting on the NO molecule. MD simulations for isotropic potentials reasonably reproduce the experimental observables for the femtosecond response and the bubble size but do not match spectroscopic results. A general overall view of the results suggests that, when the Kr-NO interaction takes place inside the matrix, potentials are rather symmetric and less repulsive than those for the triatomic system.     

Rotational spectroscopic and ab initio studies of the Xe-H2O van der Waals dimer

An ab initio potential energy surface of the Xe-H(2)O van der Waals dimer was constructed at the coupled cluster level of theory with single, double, and pertubatively included triple excitations. For the Xe atom, the small-core pseudopotential and augmented correlation-consistent polarized valence quadruple-zeta (aug-cc-pVQZ-PP) basis set was used. Dunning's augmented correlation-consistent polarized valence triple-zeta (aug-cc-pVTZ) basis set was chosen for O and H atoms. Midbond functions were used to supplement the atom-centered basis sets. Rotational spectra of the Xe-H(2)O van der Waals dimer were recorded with a pulsed-nozzle Fourier transform microwave spectrometer. Rotational transitions within two internal rotor states, namely, the 0(00) and 1(01) states, were measured and assigned. Nuclear quadrupole hyperfine structures due to the (131)Xe (I = (3)/(2)), D (I = 1) and (17)O (I = (5)/(2)) nuclei were also observed and analyzed. Information about the molecular structure and the H(2)O angular motions was extracted from the spectroscopic results with the assistance of the ab initio potential.