Current trends in the development of A. M. Butlerov’s theory of chemical structure

The chemical structure concept developed by A. M. Butlerov supplemented by the views on spatial (J. H. van’t-Hoff and J. A. Le Bel) and electronic (G. Lewis) configurations of molecules constitute the basis of the classical theory of chemical structure. The advent of quantum mechanics and development of the computer chemistry extended and enhanced the conceptual basis of theoretical chemistry, which nevertheless retains its independent value and cannot be reduced to direct physical definitions. The review deals with the evolution of the key concepts of the classical theory of chemical structure and introduction of new notions and approaches to analysis of the structure and reactivity problems, which is associated with the advent of the quantum mechanics and quantum chemistry views and methods.     

Quantum mechanical calculations of molecular properties and ideal gas heat capacity of difluoromethane

Ideal gas heat capacities are important for the calculation of caloric properties of real fluids. But as shown in the past, they are not always available with required accuracy or are still lacking. In this paper, we combine quantum mechanical calculations and statistical thermodynamics. In order to find a route to a reliable prediction of ideal gas heat capacity data, we applied various quantum mechanical methods differing by computational effort to difluoromethane (CH₂F₂). Only the structural formula and fundamental physical constants enter into the calculations. It is shown that quantum mechanics leads to accurate molecular data. Reliable experimental heat capacity data reveal that an accuracy of better than 0.5% is obtained for the ideal gas heat capacity of difluoromethane.     

Parity conservation and polarization of differential cross sections in complex-forming chemical reactions

For complex-forming chemical reactions, such as atom-diatom insertion reactions, quantum scattering and quantum statistical calculations usually predict sharp forward/backward peaks in the Differential Cross Sections (DCS). Conversely, the corresponding classical calculations are unable to reproduce these peaks. We show here that the basic reason for such an intriguing failure is that parity conservation is ignored in classical mechanics. A by-product of the analysis is a simple parity-restoring approximation that might significantly increase the ability of classical mechanics to describe DCSs over the whole angular range for the title processes.     

Periodic and high-temperature disordered conformations of polytetrafluoroethylene chains: an ab initio modeling

We report here the main results of a successful attempt to predict some macroscopic properties of representative polymers of technological relevance both in regular and disordered forms by using first principle quantum mechanical approaches at microscopic level. Until now, the prediction of the structural and thermal properties of those polymers has been mostly a domain of molecular mechanics methods. To overcome the limits of those classical computational tools whenever physical properties are significantly influenced by stereoelectronic effects (e.g., electron rich substituents), we employed methods rooted in the Density Functional Theory (DFT). A general computational strategy including the proper choice of periodic boundary conditions (PBC), functional, basis set, and model system size, has been tested and validated for saturated polymers such as polyethylene and isotactic/syndiotactic polypropylenes. On the basis of these results, a comprehensive study of poly(tetrafluoroethylene) (PTFE) chains in both regular periodic and disordered conformations has been performed. A statistical approach has been next applied to obtain the thermal concentration of defects and to reproduce the thermal behavior of the investigated polymer. At the end, a very good agreement with experimental X-ray diffraction and IR results has been achieved, definitely reaching a good understanding of the widely studied disorder phenomenon determining the main technological properties of poly(tetrafluoroethylene) (the trade Teflon) and, at the same time, identifying the proper computational tools to investigate perfluoro-compounds and other complex polymeric systems.     

计算化学实验的课程建设

With the fast development of theoretical chemistry methodologies and computer hardware and software technologies, computational chemistry has become more and more imperative in chemical science. Accordingly, we have initiated "Computational Chemistry Experiments" course based on "Introductory Computational Chemistry" course developed in the Department of Chemistry, Tsinghua University since 2013, to provide basic training and computational chemistry practices for undergraduate students. The popular quantum chemistry software Gaussian is used as the major tool in the course. We focus on the exploratory feature of computational chemistry, creating initial structures and optimizing geometries, to provide fundamental training for students to comprehend the difference of traditional experimental chemistry and modern computational chemistry. Systematic research training is offered to develop the creative thinking capacity of students in using computational chemistry methods in chemistry education and research.     

Theoretical study on aquation reaction of cis-platin complex: RISM-SCF-SEDD, a hybrid approach of accurate quantum chemical method and statistical mechanics

The ligand exchange process of cis-platin in aqueous solution was studied using RISM-SCF-SEDD (reference interaction site model-self-consistent field with spatial electron density distribution) method, a hybrid approach of quantum chemistry and statistical mechanics. The analytical nature of RISM theory enables us to compute accurate reaction free energy in aqueous solution based on CCSD(T), together with the microscopic solvation structure around the complex. We found that the solvation effect is indispensable to promote the dissociation of the chloride anion from the complex.     

Atoms and bonds in molecules and chemical explanations

The concepts of atoms and bonds in molecules which appeared in chemistry during the nineteenth century are unavoidable to explain the structure and the reactivity of the matter at a chemical level of understanding. Although they can be criticized from a strict reductionist point of view, because neither atoms nor bonds are observable in the sense of quantum mechanics, the topological and statistical interpretative approaches of quantum chemistry (quantum theory of atoms in molecules, electron localization function and maximum probability domain) provide consistent definitions which accommodate chemistry and quantum mechanics.     

基于量子化学计算方法的天然气水合物稳定性研究进展

天然气水合物以资源丰富、优质、洁净等特点,被视为21世纪新能源。天然气水合物稳定性的研究对天然气水合物资源勘探开发具有重要意义。本文简述了微观、介观、宏观、矿藏四个尺度天然气水合物稳定性的研究,重点从微观量子尺度介绍了量子化学计算方法对水合物晶体结构及其稳定性以及水合物宏观物理特性微观表征的计算研究。应用量子化学计算方法可以对天然气水合物的晶体结构、电子轨道分布、振动光谱、成键特性及主客体相互作用进行计算研究,其结果能够为天然气水合物在油气储运、水合物成藏、开采及其综合利用等方面的研究提供理论支持。目前,量子化学计算方法的优化与分子动力学模拟、分子力学模拟等方法的结合将有助于水合物形成和分解微观机理研究的发展,提升计算精度和扩大研究体系,为矿场尺度的天然气水合物资源开采利用提供理论支持。     

Early impact of quantum physics on chemistry: George Hevesy’s work on rare earth elements and Michael Polanyi’s absorption theory

After Heitler and London published their pioneering work on the application of quantum mechanics to chemistry in 1927, it became an almost unquestioned dogma that chemistry would soon disappear as a discipline of its own rights. Reductionism felt victorious in the hope of analytically describing the chemical bond and the structure of molecules. The old quantum theory has already produced a widely applied model for the structure of atoms and the explanation of the periodic system. This paper will show two examples of the entry of quantum physics into more classical fields of chemistry: inorganic chemistry and physical chemistry. Due to their professional networking, George Hevesy and Michael Polanyi found their ways to Niels Bohr and Fritz London, respectively, to cooperate in solving together some problems of classical chemistry. Their works on rare earth elements and adsorption theory throws light to the application of quantum physics outside the reductionist areas. They support the heuristic and persuasive value of quantum thinking in the 1920–1930s. Looking at Polanyi’s later oeuvre, his experience with adsorption theory could be a starting point of his non-justificationist philosophy.     

Reaction path potential for complex systems derived from combined ab initio quantum mechanical and molecular mechanical calculations

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.     

The Dirac (Bracket) Notation in the Undergraduate Physical Chemistry Curriculum: A Pictorial Introduction

The Dirac (bracket) notation is ubiquitous in the chemical literature, but it is rarely introduced in the undergraduate chemistry curriculum. In this article we present a pictorial approach to the bracket notation that we have successfully used for the past three years in a junior-senior-level physical chemistry course. We have found that it requires roughly 75 minutes to introduce this topic, and, upon integration into subsequent discussions, it prepares our undergraduate students to routinely use this powerful tool in the study of chemical bonding and spectroscopy. We believe that this approach, when introduced after the traditional integral treatment, enhances student learning of the abstract subject of quantum mechanics.     

Generalization of classical mechanics for nuclear motions on nonadiabatically coupled potential energy surfaces in chemical reactions

Classical trajectory study of nuclear motion on the Born-Oppenheimer potential energy surfaces is now one of the standard methods of chemical dynamics. In particular, this approach is inevitable in the studies of large molecular systems. However, as soon as more than a single potential energy surface is involved due to nonadiabatic coupling, such a naive application of classical mechanics loses its theoretical foundation. This is a classic and fundamental issue in the foundation of chemistry. To cope with this problem, we propose a generalization of classical mechanics that provides a path even in cases where multiple potential energy surfaces are involved in a single event and the Born-Oppenheimer approximation breaks down. This generalization is made by diagonalization of the matrix representation of nuclear forces in nonadiabatic dynamics, which is derived from a mixed quantum-classical representation of the electron-nucleus entangled Hamiltonian [Takatsuka, K. J. Chem. Phys. 2006, 124, 064111]. A manifestation of quantum fluctuation on a classical subsystem that directly contacts with a quantum subsystem is discussed. We also show that the Hamiltonian thus represented gives a theoretical foundation to examine the validity of the so-called semiclassical Ehrenfest theory (or mean-field theory) for electron quantum wavepacket dynamics, and indeed, it is pointed out that the electronic Hamiltonian to be used in this theory should be slightly modified.     

Die Molekulardynamik

This is a brief introduction to some of the basic concepts underlying the modern molecular simulations techniques, in particular the molecular dynamics (MD) method. This field belongs to the statistical mechanics branch of computational chemistry. Some applications in physical chemistry are illustrated through a few examples. The paper ends with epistemological considerations.     

分子间弱相互作用的计算——构建氢键、卤键势能面

介绍一个面向大三下学期本科生的计算化学实验。通过量子化学计算对比研究两种分子间弱相互作用——氢键和卤键的势能面,使学生对势能面的概念及两种弱相互作用的区别有一定的直观认识。通过实际上机操作,初步了解Gaussian、Gaussview及Origin等计算化学相关软件的使用,并深入理解结构化学及计算化学课程中所学的理论知识。     

量子力学和分子力学组合方法及其应用

QM/MM组合方法在研究凝聚态中的化学反应及生物大分子的结构和活性之间的关系等方面已取得重要进展。这一方法的要点在于将大体系配分成几部分,根据需要对不同部分进行不同级别的处理,因此既利用了量子力学的精确性,又利用了分子力学的高效性。对QM/MM组合理论及其一些最新进展作一简单介绍,并以最近进行了几个工作为例说明QM、MM组合方法的应用。     

Evaluation of computational chemistry methods: crystallographic and cheminformatics analysis of aminothiazole methoximes

Cheminformatics is used to validate the capabilities of widely used quantum chemistry and molecular mechanics methods. Among the quantum methods examined are the semiempirical MNDO, AM1, and PM3 methods, Hartree-Fock (ab initio) at a range of basis set levels, density functional theory (DFT) at a range of basis sets, and a post-Hartree-Fock method, local Moller-Plesset second-order perturbation theory (LMP2). Among the force fields compared are AMBER, MMFF94, MMFF94s, OPLS/A, OPLS-AA, Sybyl, and Tripos. Programs used are Spartan, MacroModel, SYBYL, and Jaguar. The test molecule is (2-amino-5-thiazolyl)-alpha-(methoxyimino)-N-methylacetamide, a model of the aminothiazole methoxime (ATMO) side chain of third-generation cephalosporin antibacterial agents. The Ward hierarchical clustering technique yields an insightful comparison of experimental (X-ray) and calculated (energy optimized) bond lengths and bond angles. The computational chemistry methods are also compared in terms of the potential energy curves they predict for internal rotation. Clustering analysis and regression analysis are compared. The MMFF94 force field such as implemented in MacroModel is the best overall computational chemistry method at reproducing crystallographic data and conformational properties of the ATMO moiety. This work demonstrates that going to a higher level of quantum theory does not necessarily give better results and that quantum mechanical results are not necessarily better than molecular mechanics results.     

Quantized Hamilton Dynamics

The paper describes the quantized Hamilton dynamics (QHD) approach that extends classical Hamiltonian dynamics and captures quantum effects, such as zero point energy, tunneling, decoherence, branching, and state-specific dynamics. The approximations are made by closures of the hierarchy of Heisenberg equations for quantum observables with the higher order observables decomposed into products of the lower order ones. The technique is applied to the vibrational energy exchange in a water molecule, the tunneling escape from a metastable state, the double-slit interference, the population transfer, dephasing and vibrational coherence transfer in a two-level system coupled to a phonon, and the scattering of a light particle off a surface phonon, where QHD is coupled to quantum mechanics in the Schrödinger representation. Generation of thermal ensembles in the extended space of QHD variables is discussed. QHD reduces to classical mechanics at the first order, closely resembles classical mechanics at the higher orders, and requires little computational effort, providing an efficient tool for treatment of the quantum effects in large systems.     

Overreact,an in silico lab: Automative quantum chemical microkinetic simulations for complex chemical reactions

Today's demand for precisely predicting chemical reactions from first principles requires research to go beyond Gibbs' free energy diagrams and consider other effects such as concentrations and quantum tunneling. The present work introduces overreact, a novel Python package for propagating chemical reactions over time using data from computational chemistry only. The overreact code infers all differential equations and parameters from a simple input that consists of a set of chemical equations and quantum chemistry package outputs for each chemical species. We evaluate some applications from the literature: gas-phase eclipsed-staggered isomerization of ethane, gas-phase umbrella inversion of ammonia, gas-phase degradation of methane by chlorine radical, and three solvation-phase reactions. Furthermore, we comment on a simple solvation-phase acid–base equilibrium. We show how it is possible to achieve reaction profiles and information matching experiments.