首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Methodology is discussed for mixed ab initio quantum mechanics/molecular mechanics modeling of systems where the quantum mechanics (QM) and molecular mechanics (MM) regions are within the same molecule. The ab initio QM calculations are at the restricted Hartree–Fock level using the pseudospectral method of the Jaguar program while the MM part is treated with the OPLS force fields implemented in the IMPACT program. The interface between the QM and MM regions, in particular, is elaborated upon, as it is dealt with by “breaking” bonds at the boundaries and using Boys-localized orbitals found from model molecules in place of the bonds. These orbitals are kept frozen during QM calculations. Results from tests of the method to find relative conformational energies and geometries of alanine dipeptides and alanine tetrapeptides are presented along with comparisons to pure QM and pure MM calculations. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1468–1494, 1999  相似文献   

2.
3.
Two‐dimensional electronic spectroscopy (2DES) is a cutting‐edge technique for investigating with high temporal resolution energy transfer, structure, and dynamics in a wide range of systems in physical chemistry, energy sciences, biophysics, and biocatalysis. However, the interpretation of 2DES is challenging and requires computational modeling. This perspective provides a roadmap for the development of computational tools that could be routinely applied to simulate 2DES spectra of multichromophoric systems active in the UV region (2DUV) using state‐of‐the‐art ab initio electronic structure methods within a quatum mechanics/molecular mechanics (QM/MM) scheme and the sum‐over‐states (SOS) approach (here called SOS//QM/MM). Multiconfigurational and multireference perturbative methods, such as the complete active space self‐consistent field and second‐order multireference perturbation theory (CASPT2) techniques, can be applied to reliably calculate the electronic properties of multichromophoric systems. Hybrid QM/MM method and molecular dynamics techniques can be used to assess environmental and conformational effects, respectively, that shape the 2D electronic spectra. DNA and proteins are important biological targets containing UV chromophores. We report ab initio simulation of 2DUV spectra of a cyclic tetrapeptide containing two interacting aromatic side chains, a model system for the study of protein structure and dynamics by means of 2DUV spectroscopy. © 2013 Wiley Periodicals, Inc.  相似文献   

4.
The ab initio multiple spawning (AIMS) method has been developed to solve the electronic and nuclear Schrodinger equations simultaneously for application to photochemical reaction dynamics. We discuss some details of the implementation of AIMS in the Molpro program package. A few aspects of the implementation are highlighted, including a new multiple timescale integrator and a scheme for solving the coupled-perturbed multiconfiguration self-consistent field (CP-MCSCF) equations in the context of ab initio molecular dynamics. The implementation is very efficient and we demonstrate calculations on the photoisomerization of ethylene using more than 5000 trajectory basis functions. We have included the capability for hybrid quantum mechanics/molecular mechanics (QM/MM) simulations within AIMS, and we investigate the role of an argon solvent in the photoisomerization of ethylene. Somewhat surprisingly, the surrounding argon has little effect on the timescale of non-adiabatic quenching in ethylene.  相似文献   

5.
A combined ab initio quantum mechanical and molecular mechanical (AI-QM/MM) potential for use in molecular modeling and simulation has been described. In this article, we summarize a procedure for deriving the empirical parameters embedded in a combined QM/MM model and suggest a set of Lennard-Jones parameters for the combined ab initio 3-21G and MM OPLS-TIP3P (AI-3/MM) potential. Interaction energies and geometrical parameters predicted with the AI-3/MM model for over 80 hydrogen-bonded complexes of organic compounds with water were found to be in good accord with ab initio 6-31G(d) results. We anticipate that the AI-3/MM potential should be reasonable for use in condensed phase simulations. © 1996 by John Wiley & Sons, Inc.  相似文献   

6.
The combination of quantum mechanics (QM) and molecular mechanics (MM) methods has become an alternative tool for many applications for which pure QM and MM are not suitable. The QM-MM method has been used for different types of problems, for example, structural biology, surface phenomena, and the liquid phase. In this paper, we have implemented these methods for vitamins, an important kind of biological molecule, and then compared results. The calculations were done by the full ab initio method (HF/3–21 g and HF/6–31 g) and QM-MM (ONIOM) method with HF(3–21 g)/AM1/UFF; then, we found that the geometry obtained by the QM-MM method is very accurate and this rapid method can be used in place of time consuming ab initio methods for large molecules. A comparison of energy values in the QM-MM and QM methods is given. We compare chemical shifts and conclude that the QM-MM method is a perturbed full QM method. The text was submitted by the authors in English.  相似文献   

7.
We present a benchmark study of a combined multipole shielding polarizability/reaction field (MSP/RF) approach to the calculation of both specific and bulk solvation effects on nuclear magnetic shielding constants of solvated molecules. The MSP/RF scheme is defined by an expansion of the shielding constants of the solvated molecule in terms of electric field and field gradient property derivatives derived from single molecule ab initio calculations. The solvent electric field and electric field gradient are calculated based on data derived from molecular dynamics simulations, thereby accounting for solute-solvent dynamical effects. The MSP/RF method is benchmarked against polarizable quantum mechanics/molecular mechanics (QM/MM) calculations. The best agreement between the MSP/RF and QM/MM approaches is found by truncating the electric field expansion in the MSP/RF approach at the linear electric field level which is due to the cancelation of errors. In addition, we investigate the sensitivity of the results due to the choice of one-electron basis set in the ab initio calculations of the property derivatives and find that these derivatives are affected by the basis set in a way similar to the shielding constants themselves.  相似文献   

8.
Born‐Oppenheimer ab initio QM/MM molecular dynamics simulation with umbrella sampling is a state‐of‐the‐art approach to calculate free energy profiles of chemical reactions in complex systems. To further improve its computational efficiency, a mass‐scaling method with the increased time step in MD simulations has been explored and tested. It is found that by increasing the hydrogen mass to 10 amu, a time step of 3 fs can be employed in ab initio QM/MM MD simulations. In all our three test cases, including two solution reactions and one enzyme reaction, the resulted reaction free energy profiles with 3 fs time step and mass scaling are found to be in excellent agreement with the corresponding simulation results using 1 fs time step and the normal mass. These results indicate that for Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, the mass‐scaling method can significantly reduce its computational cost while has little effect on the calculated free energy profiles. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009  相似文献   

9.
10.
This article describes an extension of the quantum supercharger library (QSL) to perform quantum mechanical (QM) gradient and optimization calculations as well as hybrid QM and molecular mechanical (QM/MM) molecular dynamics simulations. The integral derivatives are, after the two‐electron integrals, the most computationally expensive part of the aforementioned calculations/simulations. Algorithms are presented for accelerating the one‐ and two‐electron integral derivatives on a graphical processing unit (GPU). It is shown that a Hartree–Fock ab initio gradient calculation is up to 9.3X faster on a single GPU compared with a single central processing unit running an optimized serial version of GAMESS‐UK, which uses the efficient Schlegel method for ‐ and ‐orbitals. Benchmark QM and QM/MM molecular dynamics simulations are performed on cellobiose in vacuo and in a 39 Å water sphere (45 QM atoms and 24843 point charges, respectively) using the 6‐31G basis set. The QSL can perform 9.7 ps/day of ab initio QM dynamics and 6.4 ps/day of QM/MM dynamics on a single GPU in full double precision. © 2015 Wiley Periodicals, Inc.  相似文献   

11.
12.
13.
14.
15.
We explored water-assisted decompositions of formic acid in supercritical water in terms of local structure near reactant. A hybrid quantum mechanics/molecular mechanics (QM/MM) simulation used in this paper includes QM part as first solvation shell members around the reactant. A present QM/MM approach can simulate supercritical water solution with a reasonable computational load while keeping the simulation preciseness because a density functional theory of B3LYP/6-31+G(d) level was iterated at every 1000 Monte Carlo solute moves. The formic acid converts mainly decarboxylation by water-assisted mechanism, and the coordinated water molecules play an important role for understanding supercritical water density dependence of the reaction. We analyzed a contour map based on the solute–solvent interaction energy along with the reaction pathway. Coordinated water molecule restricted the dehydration pathway by means of hydrogen bond with formic acid, however, the coordinated water promotes the decarboxylation pathway by means of stabilization of the transition state structure with one catalytic water molecule. The contour map of the pair interaction energy along the reaction path elucidates the role of local structure on reactions in supercritical water.  相似文献   

16.
A technique for implementing the integrated molecular orbital and molecular mechanics (IMOMM) methodology developed by Maseras and Morokuma that is used to perform combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations, frequency calculations and simulations of macromolecules including explicit solvent is presented. Although the IMOMM methodology is generalized to any coordinate system, the implementation first described by Maseras and Morokuma requires that the QM and MM gradients be transformed into internal coordinates before they are added together. This coordinate transformation can be cumbersome for macromolecular systems and can become ill-defined during the course of a molecular dynamics simulation. We describe an implementation of the IMOMM method in which the QM and MM gradients are combined in the cartesian coordinate system, thereby avoiding potential problems associated with using the internal coordinate system. The implementation can be used to perform combined QM/MM molecular dynamics simulations and frequency calculations within the IMOMM framework. Finally, we have examined the applicability of thermochemical data derived from IMOMM framework. Finally, we have examined the applicability of thermochemical data derived from IMOMM frequency calculations. Received: 11 May 1998 / Accepted: 14 August 1998 / Published online: 16 November 1998  相似文献   

17.
A system of atom-atom potential functions for computer simulation of aqueous solutions of DNA fragments and counterions was developed. Hydration of Na+, K+, and dimethyl phosphate (DMP) ions was simulated by the Monte Carlo method. The obtained energy and structural characteristics of the solutions reproduce well the experimental data and are in good agreement with the results ofab initio calculations carried out by other authors. Translated fromIzvestiya Akademii Nauk Seriya Khimicheskaya, No. 11, pp. 2166–2173, November, 1998.  相似文献   

18.
Extensive calculations on hydrogen bonded systems were carried out using the improved MM3 directional hydrogen bond potential. The resulting total function was reoptimized. Comparisons of the hydrogen bonding potential function from ab initio calculations (MP2/6-31G**); the original MM3(89); and the reoptimized MM3 force field MM3(96), for a variety of C, N, O, and Cl systems including the formamide dimer and formamide–water complex, are described herein. Hydrogen bonding is shown to be a far more complicated and ubiquitous phenomenon than is generally recognized. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1001–1016, 1998  相似文献   

19.
Combined ab initio quantum mechanical and molecular mechanical calculations have been widely used for modeling chemical reactions in complex systems such as enzymes, with most applications being based on the determination of a minimum energy path connecting the reactant through the transition state to the product in the enzyme environment. However, statistical mechanics sampling and reaction dynamics calculations with a combined ab initio quantum mechanical (QM) and molecular mechanical (MM) potential are still not feasible because of the computational costs associated mainly with the ab initio quantum mechanical calculations for the QM subsystem. To address this issue, a reaction path potential energy surface is developed here for statistical mechanics and dynamics simulation of chemical reactions in enzymes and other complex systems. The reaction path potential follows the ideas from the reaction path Hamiltonian of Miller, Handy and Adams for gas phase chemical reactions but is designed specifically for large systems that are described with combined ab initio quantum mechanical and molecular mechanical methods. The reaction path potential is an analytical energy expression of the combined quantum mechanical and molecular mechanical potential energy along the minimum energy path. An expansion around the minimum energy path is made in both the nuclear and the electronic degrees of freedom for the QM subsystem internal energy, while the energy of the subsystem described with MM remains unchanged from that in the combined quantum mechanical and molecular mechanical expression and the electrostatic interaction between the QM and MM subsystems is described as the interaction of the MM charges with the QM charges. The QM charges are polarizable in response to the changes in both the MM and the QM degrees of freedom through a new response kernel developed in the present work. The input data for constructing the reaction path potential are energies, vibrational frequencies, and electron density response properties of the QM subsystem along the minimum energy path, all of which can be obtained from the combined quantum mechanical and molecular mechanical calculations. Once constructed, it costs much less for its evaluation. Thus, the reaction path potential provides a potential energy surface for rigorous statistical mechanics and reaction dynamics calculations of complex systems. As an example, the method is applied to the statistical mechanical calculations for the potential of mean force of the chemical reaction in triosephosphate isomerase.  相似文献   

20.
Two semiempirical methods (MNDO and AM1), a molecular mechanics technique (MM2) and two ab initio approaches (6–31G* full optimization and 3–21G/6–31G*) were used to calculate the ordering of and energy difference between conformers in 1,3-dichloropropane. The semiempirical methods did not order the conformers properly or predict correct energy differences. Both ab initio methods ordered the conformers and predicted energy differences correctly, with the 6–31G* full optimization performing slightly better. The MM2 results were presented for calculations involving a force field with no hydrogens and a full force field of all atoms. The full force field properly ordered the conformers but did not correctly predict the energy differences. The nonhydrogen field ordered the conformers based on the Cl…Cl nonbonded distance. The data show that conformer stability is not a simple matter of maximizing the Cl…Cl nonbonded distance, but is also related to some other stabilizing interaction(s).  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号