用第一性原理方法获取周期体系中原子的部分电荷

利用量子化学软件包Crystal计算了立方周期性边界条件下液态水体系的静电势(ESP)和静电场(EF).在此基础上,提出了一种由第一性原理方法获取周期体系中原子的部分电荷的快捷方法.该方法把由周期性边界条件引入的平均静电势准mean作为一个拟合参数,通过对第一性原理静电势与Ewald加和法静电势的最小二乘法拟合而实现.值得说明的是,比较静电势与静电场拟合方法,前者的相对拟合误差仅为2%-3%,比后者小一个数量级.考察了四种电荷限制条件下,静电势、静电场拟合的水分子原子部分电荷及偶极矩的分布情况.     

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.     

Hybrid quantum/classical path integral approach for simulation of hydrogen transfer reactions in enzymes

A hybrid quantum/classical path integral Monte Carlo (QC-PIMC) method for calculating the quantum free energy barrier for hydrogen transfer reactions in condensed phases is presented. In this approach, the classical potential of mean force along a collective reaction coordinate is calculated using umbrella sampling techniques in conjunction with molecular dynamics trajectories propagated according to a mapping potential. The quantum contribution is determined for each configuration along the classical trajectory with path integral Monte Carlo calculations in which the beads move according to an effective mapping potential. This type of path integral calculation does not utilize the centroid constraint and can lead to more efficient sampling of the relevant region of conformational space than free-particle path integral sampling. The QC-PIMC method is computationally practical for large systems because the path integral sampling for the quantum nuclei is performed separately from the classical molecular dynamics sampling of the entire system. The utility of the QC-PIMC method is illustrated by an application to hydride transfer in the enzyme dihydrofolate reductase. A comparison of this method to the quantized classical path and grid-based methods for this system is presented.     

Quasiclassical trajectory calculations to evaluate a kinematic constraint on internal energy in suprathreshold collision energy abstraction reactions

Experimentally observed product quantum state distributions across a wide range of abstraction reactions at suprathreshold collision energies have shown a strong bias against product internal energy. Only a fraction, sometimes quite a small fraction, of the energetically accessible product quantum states are populated. Picconatto et al. [J. Chem. Phys. 114, 1663 (2001)] noted a simple mathematical relationship between the highest-energy rovibrational states observed and the kinematics of the reaction system. They proposed a reaction model based on reaction kinematics that quantitatively explains this behavior. The model is in excellent agreement with measured quantum state distributions. The assumptions of the model invoke detailed characteristics of reactive trajectories at suprathreshold collision energies. Here we test those assumptions using quasiclassical trajectory calculations for the abstraction reactions H+HCl-->H2+Cl, D+HCl-->HD+Cl, and H+DCl-->HD+Cl. Trajectories were run on a potential-energy surface calculated with a London-Eyring-Polyani-Sato function with a localized 3-center term (LEPS-3C) previously shown to accurately reproduce experimentally observed product state distributions for the H+HCl abstraction reaction. The trajectories sample collision energies near threshold and also substantially above it. Although the trajectories demonstrate some aspects of the model, they show that it is not valid. However, the inadequacy of the proposed model does not invalidate the apparent kinematic basis of the observed energy constraint. The present results show that there must be some other molecular behavior rooted in the reaction kinematics that is the explanation and the source of the constraint.     

Convergence of the orbital expansion in a correlated system: A test study on positronium

Positronium, the bound state of an electron and a positron, is an exactly soluble quantum system, similar to a light isotope of hydrogen. It can be studied using the finite basis quantum chemistry codes developed for atoms and molecules. In fact, positronium can be mimicked by two electrons with opposite spins, in the absence of any nucleus and having the sign of the Coulomb interaction reversed. The exact wave function has a cusp in the points of coalescence of the two particles (a “Coulomb peak”), and this fact makes the convergence of the total energy, as a function of the basis set size, extremely slow. For this reason, positronium can be used to test the convergence properties of the quantum chemistry methods used to describe the dynamic correlation. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Prediction of nitroxide hyperfine coupling constants in solution from combined nanosecond scale simulations and quantum computations

We present a combined theoretical approach based on analyzing molecular dynamics trajectories (at the nanosecond scale) generated by use of classical polarizable force fields and on quantum calculations to compute averaged hyperfine coupling constants. That method is used to estimate the constant of a prototypical nitroxide: the dimethylnitroxide. The molecule is embedded during the simulations in a cubic box containing about 500 water molecules and the molecular dynamics is generated using periodic conditions. Once the trajectories are achieved, the nitroxide and its first hydration shell molecules are extracted, and the coupling constants are computed by considering the latter aggregates by means of quantum computations. However, all the water molecules of the bulk are also accounted for during those computations by means of the electrostatic potential fitted method. Our results exhibit that in order to predict accurate and reliable coupling constants, one needs to describe carefully the out-of-plane motion of the nitroxide nitrogen and to sample trajectories with a time interval of 400 fs at least to generate an uncorrelated large set of nitroxide structures. Compared to Car-Parrinello molecular dynamics techniques, our approach can be used readily to compute hyperfine coupling constants of large systems, such as nitroxides of great size interacting with macromolecules such as proteins or polymers.     

Investigating detailed mechanism of hydrogen molecules adsorbing on single-wall carbon nanotubes using fitted force field parameters containing carbon–carbon interactions

The nonbonded and bonded force field parameters for carbon atoms in single-wall carbon nanotubes (SWNT) are fitted by means of quantum chemistry calculations with considering the periodic boundary conditions. The nonbonded parameters between carbon atoms and hydrogen atoms are fitted as well. All the fitted parameters are verified by comparing to quantum chemistry results and by calculating Young's modulus. Adsorption of Hydrogen molecules are then carried out on a bundle of self-assembled SWNTs. The adsorption isotherms are consistent to the Freundlich equation. Both hydrogen molecules adsorbed outside and inside the SWNTs are counted. According to our result, hydrogen molecules adsorbed inside the SWNTs are more stable at a relatively high temperature and are playing an important part in total amount of the adsorbed molecules. While C(10,10) have the highest adsorption capacities in most of the temperatures, hydrogen molecules inside C(5,5) are the most stable of all the four kinds of SWNTs. Thus, balancing adsorption capacities and strength of interaction can be important in choosing SWNT for gas adsorption. Besides, we deduct an equation that can describe the relation between hydrogen pressure and amount of SWNTs based on our simulation results. The hydrogen pressure may decrease by adding SWNTs in the system. The fitting method in our system is valid to SWNTs and can be tested in further studies of similar systems. © 2018 Wiley Periodicals, Inc.     

Computing the density of paths in complex systems

Trajectories of duration tau joining two points q(0) and q(1) in the configuration space of a classical system satisfy Hamilton's principle: they are stationary points of the classical action. The second variation (fluctuations) of the action around the stationary point signals whether the latter is or not a minimum and delivers the density in trajectory space around the points q(0) and q(1). This concept of density of paths is of great importance in semiclassical quantum theory, where it weights the contribution to the propagator from the single classical trajectories. In this paper, two algorithms based on the concepts of molecular dynamics simulation are introduced for computing the density of paths, also called van Vleck [Proc. Natl. Acad. Sci. U.S.A. 14, 178 (1928)] determinant. Examples for realistic systems are presented, together with a suggestion about possible applications in the field of rare events in physics and chemistry.     

Semiclassical initial value series representation in the continuum limit: application to vibrational relaxation

A recently formulated continuum limit semiclassical initial value series representation (SCIVR) of the quantum dynamics of dissipative systems is applied to the study of vibrational relaxation of model harmonic and anharmonic oscillator systems. As is well known, the classical dynamics of dissipative systems may be described in terms of a generalized Langevin equation. The continuum limit SCIVR uses the Langevin trajectories as input, albeit with a quantum noise rather than a classical noise. Combining this development with the forward-backward form of the prefactor-free propagator leads to a tractable scheme for computing quantum thermal correlation functions. Here we present the first implementation of this continuum limit SCIVR series method to study two model problems of vibrational relaxation. Simulations of the dissipative harmonic oscillator system over a wide range of parameters demonstrate that at most only the first two terms in the SCIVR series are needed for convergence of the correlation function. The methodology is then applied to the vibrational relaxation of a dissipative Morse oscillator. Here, too, the SCIVR series converges rapidly as the first two terms are sufficient to provide the quantum mechanical relaxation with an estimated accuracy on the order of a few percent. The results in this case are compared with computations obtained using the classical Wigner approximation for the relaxation dynamics.     

L'Etude Théorique des Systèmes Périodiques. I. La Méthode LCAO?HCO

The LCAO ? HCO method is described. Starting from classical procedures of quantum chemistry, this method enables the theoretical study of periodic systems. The explicit form of some properties is deduced in the Hückel approximation. This method is compared with the tight-binding approximation and applied to the study of polyenes, polyacenes, and graphite. Finally, the general validity of the process is discussed.     

Fast quantum crystallography

Typical contemporary X-ray crystallography delivers the geometries and, at best, the electron densities of molecules or periodic systems in the crystalline phase. Energies, electron momentum densities, and information relating to the pair density such as electron delocalization measures—all crucial to chemistry—are simply missed. Quantum crystallography (QCr) is an emerging line of research aimed at filling this gap by solving the crystallographic problem under the constraints of quantum mechanics. In this way, not only geometries and electron densities become experimentally accessible but also the entire panoply of quantum mechanical properties that are in the output of any quantum chemical software package. However, QCr remains limited to smaller systems (small molecules or small unit cells) due to the exponential bottleneck that plagues quantum mechanical calculations. When combined with a fragmentation technique, termed the “kernel energy method (KEM)”, QCr's reach to larger molecules is extended considerably to almost “any size”, that is, systems of up to many hundreds of thousands of atoms. KEM has made this doable with any chemical model and is capable of providing the entire quantum mechanics of large molecular systems. The smallness of the R-factor adjudicates the accuracy of the quantum mechanics extracted from the crystallography.     

GridMol2.0: Implementation and application of linear-scale quantum mechanics methods and molecular visualization

GridMol is a “one-stop” platform for molecular structure building, scientific computing, and molecular visualization aided by a high-performance computing environment. GridMol version 2.0 introduces two unique features: the first is fragment-based linear-scaling quantum chemistry methods, such as molecular fractionation with conjugate caps and fragment molecular orbital methods; the second is that GridMol enables users to visualize molecular geometries along a geometry optimization and an intrinsic reaction coordinate calculation. Compared with version 1.0, fragment-based linear-scaling quantum chemistry methods implemented in GridMol version 2.0 can be used as a useful tool for performing quantum calculations for large molecular systems to explore the mechanisms involved in protein-ligand or targeted drug interactions.     

Applications of a quantum virial theorem to kronig and penney's model and to a diatomic molecule in static approximation

Kronig and Penney's model is used to test a new formulation of quantum virial theorem for periodic systems. It is shown that the corrective term to the usual formulation of the virial theorem corresponding to boundary conditions is of importance and gives correct energy values for the system. Application of the virial theorem to a diatomic molecule also shows an important boundary condition correction if the distance of the nuclei is slightly different from equilibrium position, in accordance with the above correction.     

Noncovalent interactions in inorganic supramolecular chemistry based in heavy metals. Quantum chemistry point of view

Complexity is a concept that is being considered in chemistry as it has shown potential to reveal interesting phenomena. Thus, it is possible to study chemical phenomena in a new approach called systems chemistry. The systems chemistry has an organization and function, which are regulated by the interactions among its components. At the simplest level, noncovalent interactions between molecules can lead to the emergence of large structures. Consequently, it is possible to go from the molecular to the supramolecular systems chemistry, which aims to develop chemical systems highly complex through intra- and intermolecular forces. Proper use of the interactions previously mentioned allow a glimpse of supramolecular system chemistry in many tasks such as structural properties reflecting certain behaviors in the chemistry of materials, for example, electrical and optical, processes of molecular recognition and among others. In the last time, within this area, inorganic supramolecular systems chemistry has been developed. Those systems have a structural orientation which is defined by certain forces that predominate in the associations among molecules. It is possible to recognize these forces as hydrogen bonding, π-π stacking, halogen bonding, electrostatic, hydrophobic, charge transfer, metal coordination, and metallophilic interactions. The presence of these forces in supramolecular system yields certain properties such as light absorption and luminescence. The quantum theoretical modeling plays an important role in the designing of the supramolecular system. The goal is to apply supramolecular principles in order to understand the associated forces in many inorganic molecules that include heavy metals for instance gold, platinum, and mercury. Relevant systems will be studied in detail, considering functional aspects such as enhanced coordination of functionalized molecular self-assembly, electronic and optoelectronic properties.     

Quantum manifestations of classical trajectories in molecular systems

Quantum chaos, understood as the effect of the underlying classical dynamics on the stationary quantum properties in classically chaotic systems, is examined in two molecular floppy systems. Realistic models of two degrees of freedom for HO₂ and HCN/HNC are considered. The structure of the classical phase space is studied using Poincaré surfaces of section and the dynamical characteristics of the corresponding wave functions analyzed also in phase space with the aid of Husimi functions. Some wave functions show strong localization along periodic orbits. © 2002 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Non-adiabatic molecular dynamics with quantum solvent effects

Three novel approaches extending quantum-classical non-adiabatic (NA) molecular dynamics (MD) to include quantum effects of solvent environments are described. In a standard NA-MD the solute subsystem is treated quantum mechanically, while the larger solvent part of a system is treated classically. The three novel approaches presented here are based on the Bohmian formulation of quantum mechanics, the stochastic Schrödinger equation for the evolution of open quantum systems and the quantized Hamilton dynamics generalization of classical mechanics. The approaches extend the standard NA-MD to incorporate the following quantum effects of the solvent. (1) Branching, i.e. the ability of solvent quantum wave packets to split and follow asymptotically diverging trajectories correlated with different quantum states of the solute. (2) Decoherence, i.e. loss of quantum interference within the solute subsystem induced by the diverging solvent trajectories. (3) Zero point energy that contributes to NA coupling and must be preserved during the energy exchange between solvent and solute degrees of freedom. The Bohmian quantum-classical mechanics, stochastic mean-field and quantized mean-field approximations incorporate the quantum solvent effects into the standard quantum-classical NA-MD in a straightforward and efficient way that can be easily applied to quantum dynamics of condensed phase chemical systems.     

Incorporation of quantum effects for selected degrees of freedom into the trajectory-based dynamics using spatial domains

The approach of defining quantum corrections on nuclear dynamics of molecular systems incorporated approximately into selected degrees of freedom, is described. The approach is based on the Madelung-de-Broglie-Bohm formulation of time-dependent quantum mechanics which represents a wavefunction in terms of an ensemble of trajectories. The trajectories follow classical laws of motion except that the quantum potential, dependent on the wavefunction amplitude and its derivatives, is added to the external, classical potential. In this framework the quantum potential, determined approximately for practical reasons, is included only into the "quantum" degrees of freedom describing light particles such as protons, while neglecting with the quantum force for the heavy, nearly classical nuclei. The entire system comprised of light and heavy particles is described by a single wavefunction of full dimensionality. The coordinate space of heavy particles is divided into spatial domains or subspaces. The quantum force acting on the light particles is determined for each domain of similar configurations of the heavy nuclei. This approach effectively introduces parametric dependence of the reduced dimensionality quantum force, on classical degrees of freedom. This strategy improves accuracy of the quantum force and does not restrict interaction between the domains. The concept is illustrated for two-dimensional scattering systems, where the quantum force is required to reproduce vibrational energy of the quantum degree of freedom.     

Conical intersections and semiclassical trajectories: comparison to accurate quantum dynamics and analyses of the trajectories

Semiclassical trajectory methods are tested for electronically nonadiabatic systems with conical intersections. Five triatomic model systems are presented, and each system features two electronic states that intersect via a seam of conical intersections (CIs). Fully converged, full-dimensional quantum mechanical scattering calculations are carried out for all five systems at energies that allow for electronic de-excitation via the seam of CIs. Several semiclassical trajectory methods are tested against the accurate quantum mechanical results. For four of the five model systems, the diabatic representation is the preferred (most accurate) representation for semiclassical trajectories, as correctly predicted by the Calaveras County criterion. Four surface hopping methods are tested and have overall relative errors of 40%-60%. The semiclassical Ehrenfest method has an overall error of 66%, and the self-consistent decay of mixing (SCDM) and coherent switches with decay of mixing (CSDM) methods are the most accurate methods overall with relative errors of approximately 32%. Furthermore, the CSDM method is less representation dependent than both the SCDM and the surface hopping methods, making it the preferred semiclassical trajectory method. Finally, the behavior of semiclassical trajectories near conical intersections is discussed.