Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.
A benchmark set of 28 medium-sized organic molecules is assembled that covers the most important classes of chromophores including polyenes and other unsaturated aliphatic compounds, aromatic hydrocarbons, heterocycles, carbonyl compounds, and nucleobases. Vertical excitation energies and one-electron properties are computed for the valence excited states of these molecules using both multiconfigurational second-order perturbation theory, CASPT2, and a hierarchy of coupled cluster methods, CC2, CCSD, and CC3. The calculations are done at identical geometries (MP26-31G*) and with the same basis set (TZVP). In most cases, the CC3 results are very close to the CASPT2 results, whereas there are larger deviations with CC2 and CCSD, especially in singlet excited states that are not dominated by single excitations. Statistical evaluations of the calculated vertical excitation energies for 223 states are presented and discussed in order to assess the relative merits of the applied methods. CC2 reproduces the CC3 reference data for the singlets better than CCSD. On the basis of the current computational results and an extensive survey of the literature, we propose best estimates for the energies of 104 singlet and 63 triplet excited states.
A hierarchy of three models for combined quantum mechanical (QM) and molecular mechanical (MM) approaches is presented. They simplify the QM description of large molecules by reducing it to the electronically important fragment which interacts with the molecular mechanically treated remainder of the molecule. In the simplest model A, the QM fragments are only mechanically embedded in their MM environment. The more refined models B and C include a quantum mechanical treatment of electrostatic interactions between the fragments and a semiclassical description of polarization. The implementation of models A-C for MNDO type wavefunctions and the MM3 force field is outlined. Selected applications in organic chemistry are discussed, addressing the ability of the proposed models to reproduce substituent effects (MM) on chemical structure and reactivity (QM). These applications include protonations, deprotonations, hydride transfer reactions, nucleophilic additions, and nucleophilic ring cleavage reactions
Time-dependent density functional theory (TD-DFT) and DFT-based multireference configuration interaction (DFT/MRCI) calculations are reported for a recently proposed benchmark set of 28 medium-sized organic molecules. Vertical excitation energies, oscillator strengths, and excited-state dipole moments are computed using the same geometries (MP2/6-31G(*)) and basis set (TZVP) as in our previous ab initio benchmark study on electronically excited states. The results from TD-DFT (with the functionals BP86, B3LYP, and BHLYP) and from DFT/MRCI are compared against the previous high-level ab initio results, and, in particular, against the proposed best estimates for 104 singlet and 63 triplet vertical excitation energies. The statistical evaluation for the latter reference data gives the lowest mean absolute deviations for DFT/MRCI (0.22 eV for singlets and 0.24 eV for triplets) followed by TD-DFT/B3LYP (0.27 and 0.44 eV, respectively), whereas TD-DFT/BP86 and TD-DFT/BHLYP are significantly less accurate. The energies of singlet states with double excitation character are generally overestimated by TD-DFT, whereas triplet state energies are systematically underestimated by the currently investigated DFT-based methods.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.