scholarly journals Constructing Potential Energy Surfaces for Polyatomic Systems: Recent Progress and New Problems

2012 ◽  
Vol 2012 ◽  
pp. 1-19 ◽  
Author(s):  
J. Espinosa-Garcia ◽  
M. Monge-Palacios ◽  
J. C. Corchado

Different methods of constructing potential energy surfaces in polyatomic systems are reviewed, with the emphasis put on fitting, interpolation, and analytical (defined by functional forms) approaches, based on quantum chemistry electronic structure calculations. The different approaches are reviewed first, followed by a comparison using the benchmark H + CH4 and the H + NH3 gas-phase hydrogen abstraction reactions. Different kinetics and dynamics properties are analyzed for these reactions and compared with the available experimental data, which permits one to estimate the advantages and disadvantages of each method. Finally, we analyze different problems with increasing difficulty in the potential energy construction: spin-orbit coupling, molecular size, and more complicated reactions with several maxima and minima, which test the soundness and general applicability of each method. We conclude that, although the field of small systems, typically atom-diatom, is mature, there still remains much work to be done in the field of polyatomic systems.

Author(s):  
Lionel Raff ◽  
Ranga Komanduri ◽  
Martin Hagan ◽  
Satish Bukkapatnam

When the system of interest becomes too complex to permit the use of ab initio methods to obtain the system potential-energy surfaces (PES), empirical potential surfaces are frequently employed to represent the force fields present in the system under investigation. In most cases, the functional forms present in these potentials are selected on the basis of chemical and physical intuitions. The parameters of the surface are frequently adjusted to fit a very small set of experimental data that comprise bond energies, equilibrium bond distances and angles, fundamental vibrational frequencies, and perhaps measured barrier heights to reactions of interest. Such potentials generally yield only qualitative or semiquantitative descriptions of the system dynamics. Several research groups have significantly improved the accuracy of the values of the experimental properties computed using empirical potential surfaces by fitting the chosen functional form for the potential to the force fields obtained from trajectories using ab initio Car-Parrinello molecular dynamics simulations. The fitting to the force fields is usually done using a least-squares fitting approach. This method has been employed by Izvekov et al. to obtain effective non-polarizable three-site force fields for liquid water. Carré et al. have employed such a procedure to obtain a new pair potential for silica. In their investigation, the vector of potential parameters was fitted using an iterative Levenberg-Marquardt algorithm. Tangney and Scandolo have also developed an interatomic force field for liquid SiO2 in which the parameters were fitted to the forces, stresses, and energies obtained from ab initio calculations. Ercolessi and Adams have used a quasi-Newtonian procedure to fit an empirical potential for aluminum to data obtained from first-principals computations. Empirical potentials can be improved by making the parameters parameterized functions of the coordinates defining the instantaneous positions of the atoms of the system. This approach has been successfully employed by numerous investigators The difficulty with this procedure is that the number of parameters that must be adjusted increases rapidly. Appropriate fitting of these parameters requires a much more extensive database. Finally, the actual fitting process can often be tedious, difficult, and time-consuming.


2019 ◽  
Author(s):  
Marc Riera ◽  
Eric Yeh ◽  
Francesco Paesani

<div> <div> <div> <p>In this study, we extend the scope of the many-body TTM-nrg and MB-nrg potential energy functions (PEFs), originally introduced for halide ion–water and alkali-metal ion–water interactions, to the modeling of carbon dioxide (CO<sub>2</sub>) and water (H<sub>2</sub>O) mixtures as prototypical examples of molecular fluids. Both TTM-nrg and MB-nrg PEFs are derived entirely from electronic structure data obtained at the coupled cluster level of theory and are, by construction, compatible with MB-pol, a many-body PEF that has been shown to accurately reproduce the properties of water. Although both TTM-nrg and MB-nrg PEFs adopt the same functional forms for describing permanent electrostatics, polarization, and dispersion, they differ in the representation of short-range contributions, with the TTM-nrg PEFs relying on conventional Born-Mayer expressions and the MB-nrg PEFs employing multidimensional permutationally invariant polynomials. By providing a physically correct description of many-body effects at both short and long ranges, the MB-nrg PEFs are shown to quantitatively represent the global potential energy surfaces of the CO<sub>2</sub>–CO<sub>2</sub> and CO<sub>2</sub>–H<sub>2</sub>O dimers and the energetics of small clusters as well as to correctly reproduce various properties in both gas and liquid phases. Building upon previous studies of aqueous systems, our analysis provides further evidence for the accuracy and efficiency of the MB-nrg framework in representing molecular interactions in fluid mixtures at different temperature and pressure conditions. </p> </div> </div> </div>


By more complex systems we mean systems containing on the order of hundreds or thousands of atoms, or molecules with less atoms but with “complicated” motions, the latter being the case when considering collisions between polyatomic molecules. In the present chapter we deal with quantum-classical methods for treating energy transfer in collisions involving polyatomic molecules, molecule surface scattering, reactions in polyatomic systems and solution. We will assume that it is possible to construct realistical potential energy surfaces for the systems. Obviously, these surfaces will be of empirical or semi-empirical nature. In some of the methods, as for instance the reaction path method, one tries to minimize the information needed on potential energy surfaces. Chemical reactions and energy transfer processes in the gas phase are often studied using just a single adiabatic Born-Oppenheimer potential energy surface. However non-adiabatic effects, that is, coupling between different electronic states, is an important aspect in chemistry. If the coupling between the various electronic states can be neglected, the “electronic” effect reduces to that of a statistical degeneracy factor ge [180].


2019 ◽  
Author(s):  
Marc Riera ◽  
Eric Yeh ◽  
Francesco Paesani

<div> <div> <div> <p>In this study, we extend the scope of the many-body TTM-nrg and MB-nrg potential energy functions (PEFs), originally introduced for halide ion–water and alkali-metal ion–water interactions, to the modeling of carbon dioxide (CO<sub>2</sub>) and water (H<sub>2</sub>O) mixtures as prototypical examples of molecular fluids. Both TTM-nrg and MB-nrg PEFs are derived entirely from electronic structure data obtained at the coupled cluster level of theory and are, by construction, compatible with MB-pol, a many-body PEF that has been shown to accurately reproduce the properties of water. Although both TTM-nrg and MB-nrg PEFs adopt the same functional forms for describing permanent electrostatics, polarization, and dispersion, they differ in the representation of short-range contributions, with the TTM-nrg PEFs relying on conventional Born-Mayer expressions and the MB-nrg PEFs employing multidimensional permutationally invariant polynomials. By providing a physically correct description of many-body effects at both short and long ranges, the MB-nrg PEFs are shown to quantitatively represent the global potential energy surfaces of the CO<sub>2</sub>–CO<sub>2</sub> and CO<sub>2</sub>–H<sub>2</sub>O dimers and the energetics of small clusters as well as to correctly reproduce various properties in both gas and liquid phases. Building upon previous studies of aqueous systems, our analysis provides further evidence for the accuracy and efficiency of the MB-nrg framework in representing molecular interactions in fluid mixtures at different temperature and pressure conditions. </p> </div> </div> </div>


Sign in / Sign up

Export Citation Format

Share Document