scispace - formally typeset
Search or ask a question
Journal ArticleDOI

Calculation of ground and excited state potential surfaces of conjugated molecules. I. Formulation and parametrization

01 Aug 1972-Journal of the American Chemical Society (American Chemical Society)-Vol. 94, Iss: 16, pp 5612-5625
TL;DR: In this article, Levitt and Levitt developed a method for the consistent calculation of ground and excited state potential surfaces of conjugated molecules, which is based on the formal separation of u and 7r electrons, the former being represented by an empirical potential function and the latter by a semi-empirical model of the Pariser-Parr-Pople type corrected for nearest-neighbor orbital overlap.
Abstract: A formulation is developed for the consistent calculation of ground and excited state potential surfaces of conjugated molecules. The method is based on the formal separation of u and 7r electrons, the former being represented by an empirical potential function and the latter by a semiempirical model of the Pariser-Parr-Pople type corrected for nearest-neighbor orbital overlap. A single parameter set is used to represent all of the molecular properties considered; these include atomization energies, electronic excitation energies, ionization potentials, and the equilibrium geometries and vibrational frequencies of the ground and excited electronic states, and take account of all bond length and bond angle variations. To permit rapid determination of the potential surfaces, the u potential function and SCF-MO-CI energy of the r electrons are expressed as analytic functions of the molecular coordinates from which the first and second derivatives can be obtained. Illustrative applications to 1,3butadiene, 1,3,5-hexatriene, a,w-diphenyloctatetraene, and 1,3-cyclohexadiene are given. detailed interpretation of electronic transitions and A concomitant photochemical processes in conjugated molecules requires a knowledge of the ground and excited state potential surfaces. The determination of such surfaces has long been a goal of theoretical chemistry. Difficulties in a reliable a priori approach to the problem for a system as simple as ethylene2 are such that calculations for more complicated molecules are prohibitive at present. Consequently, a variety of methods that utilize experimental data have been introduced. Completely empirical treatments, in which the energy surface is expressed as a function of potential parameters fitted to the available information (1) Supported in part by Grant EY00062 from the National Institute of Health. (2) U. Kaldor and I. Shavitt, J . Chem. Phys., 48, 191 (1968); R. J. Buenker, S. D. Peyerimhoff, and W. E. Kammer, ibid., 55, 814 (1971). (equilibrium geometry, vibrational frequencies, etc.), have had considerable success in applications to molecules for which a localized electron description is app l i~ab le .~ The great advantage of this type of approach, which leaves open questions of reliability when extended from one class of molecules to another, is the ease and speed of the calculations; this had made possible applications to systems as large as certain nucleic acids and globular proteins. For conjugated molecules, however, the importance of delocalization introduces difficulties into such an empirical treatmenL5 (3) (a) See, for example, J. E. Williams, P. J . Stand, and P. v. R. Schleyer, Annu. Reu. Phys. Chem., 19, 531 (1969); (b) S. Lifson and A. Warshel, J . Chem. Phys., 49, 5116 (1968); A. Warshel and S . Lifson, ibid., 53, 8582 (1970). (4) M. Levitt and S. Lifson, J. Mol. B i d , 46, 269 (1969); M. Levitt, Nature (London), 224, 759 (1969). ( 5 ) C. Tric, J . Chem. Phys., 5 1 , 4778 (1969). Journal of the American Chemical Society 1 94:16 1 August 9, 1972
Citations
More filters
Journal ArticleDOI
TL;DR: An overview of the CHARMM program as it exists today is provided with an emphasis on developments since the publication of the original CHARMM article in 1983.
Abstract: CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecu- lar simulation program. It has been developed over the last three decades with a primary focus on molecules of bio- logical interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estima- tors, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numer- ous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

7,035 citations

Journal ArticleDOI
TL;DR: Electrostatic stabilization is an important factor in increasing the rate of the reaction step that leads to the formation of the carbonium ion intermediate, found in the cleavage of a glycosidic bond by lysozyme.

3,951 citations

Journal ArticleDOI
TL;DR: In this paper, a combined quantum mechanical and molecular mechanical potential has been developed for the study of reactions in condensed phases, where semi-empirical methods of the MNDO and AM1 type are used, while the molecular mechanics part is treated with the CHARMM force field.
Abstract: A combined quantum mechanical (QM) and molecular mechanical (MM) potential has been developed for the study of reactions in condensed phases. For the quantum mechanical calculations semiempirical methods of the MNDO and AM1 type are used, while the molecular mechanics part is treated with the CHARMM force field. Specific prescriptions are given for the interactions between the QM and MM portions of the system; cases in which the QM and MM methodology is applied to parts of the same molecule or to different molecules are considered. The details of the method and a range of test calculations, including comparisons with ab initio and experimental results, are given. It is found that in many cases satisfactory results are obtained. However, there are limitations to this type of approach, some of which arise from the AM1 or MNDO methods themselves and others from the present QM/MM implementation. This suggests that it is important to test the applicability of the method to each particular case prior to its use. Possible areas of improvement in the methodology are discussed.

2,197 citations

Journal ArticleDOI
TL;DR: A new computational scheme integrating ab initio and molecular mechanics descriptions in different parts of the same molecule is presented, which evaluates explicitly the energy derivatives with respect to geometrical parameters and therefore has a straightforward application to geometry optimization.
Abstract: A new computational scheme integrating ab initio and molecular mechanics descriptions in different parts of the same molecule is presented. In contrast with previous approaches, this method is especially designed to allow the introduction of molecular mechanics corrections in full geometry optimizations concerning problems usually studied through ab initio calculations on model systems. The scheme proposed in this article intends to solve some of the systematic error associated with modeling through the use of molecular mechanics corrections. This method, which does not require any new parameter, evaluates explicitly the energy derivatives with respect to geometrical parameters and therefore has a straightforward application to geometry optimization. Examples of its performance on two simple cases are provided: the equilibrium geometry of cyclopropene and the energy barriers on SN2 reactions of alkyl chloride systems. Results are in satisfactory agreement with those of full ab initio calculations in both cases. © 1995 by John Wiley & Sons, Inc.

1,565 citations

Journal ArticleDOI
TL;DR: In this article, an approach to couple ab initio quantum mechanical geometry optimiuzations with molecular mechanical optimizations is presented, with the added capability to carry out molecular dynamics simulations of the systems to earch for new local minima.
Abstract: We present an approach to couple ab initio quantum mechanical geometry optimiuzations with molecular mechanical optimizations, with the added capability to carry out molecular dynamics simulations of the systems to earch for new local minima. The approach is applied to the aqueous solution CH3Cl + Cl− exchange reaction and the gas phase protonation of polyethers.

1,046 citations