Electronic structure methods dynamics

PAW is a recent addition to the all-electron electronic structure methods whose accuracy appears to be similar to that of the general potential LAPW approach. The implementation of the molecular dynamics fonnalism enables easy stmcture optimization in this method. 

The electronic structure method used to provide the energies and gradients of the states is crucial in photochemistry and photophysics. Ab initio electronic structure methods have been used for many years. Treating closed shell systems in their ground state is a problem that, in many cases, can now be solved routinely by chemists using standardized methods and computer packages. In order to obtain quantitative results, electron correlation (also referred to as dynamical correlation) should be included in the model and there are many methods available for doing this based on either variational or perturbation principles [41],... 

The topic of this review, reactions at metal surfaces, has been in general treated in a similar way to gas-phase reactivity. High level ab initio electronic structure methods are used to construct potential energy surfaces of catalytically important surface reactions in reduced dimensions. Once a chemically accurate potential surface is available, quantum or classical dynamics may be carried out in order to more deeply understand the microscopic nature of the reaction. 

The photolysis experiments have also received theoretical attention, with electronic structure methods used to calculate the nature of the excited states (44), as well as the potential energy curves for loss of CO (45). Theoretical models for the excited-state dynamics leading to dissociation have also been proposed (46). 

The Car-Parrinello approach combines an electronic structure method with a classical molecular dynamics scheme and thus unifies two major fields of computational chemistry, which have hitherto been essentially orthogonal. Through this unification a... 

Another major, future advance in the quantum chemical computation of potential energy surfaces for reaction dynamics will be the ability to routinely compute the energies of molecular systems on the fly . The tedious and time-consuming process of fitting computed quantum chemical values to functional forms could be avoided if it were possible to compute the PES as needed during a classical trajectory or quantum dynamics calculation. For many chemical reactions, it should be practical in the near future to prudently select a sufficiently rapid and accurate electronic structure method to facilitate dynamics computations on the fly. 

Although the theoretical roots of this technique are very well established, it is more often used as a flexible surface which can be adjusted to lit either exprimental data or data established by better electronic-structure methods. The LEPS formalism has also been extensively used to explore the relationships between the potential energy surface and the details of chemical dynamics . Because of the widespread use of this potential for studying gas-phase reactions, the specific form of the equations will not be discussed here. The interested reader is instead referred to references which discuss this approach in more detail . ... 

The major drawback for employing the Car-Parrinello approach in dynamics simulations is that since a variational wavefunction is required, the electronic energy should in principle be minimized before the forces on the atoms are calculated. This greatly increases the amount of computer time required at each step of the simulation. Furthermore, the energies calculated with the electronic structure methods currently used in this approach are not exceptionally accurate. For example, it is well established that potential energy barriers, which are of importance to chemical reactivity, often require sophisticated methods to be accurately determined. Nonetheless, the Tirst-principles calculation of the forces during the dynamics is an appealing idea, and will continue to be developed as computer resources expand. 

One solution to such problems is provided by ab initio molecular dynamics (AIMD) simulations (11,53-57). These are still very expensive, but they include dynamical aspects and they are thus one of the most powerful tools to investigate bioinorganic systems if many conformations and unexpected events play a role. A necessary prerequisite will always be that inexpensive electronic structure methods such as DFT can be applied to the bioinorganic system of interest. 

At several points in diis book, it has been emphasized that the prevalence of classical MC and MD simulations derives from die impracticality of carrying out fully QM dynamics. While diis is largely true, for systems of only modest size where short trajectories may be profitably analyzed, fully QM MD simulations using the so-called Car-Parrinello technique are a viable option (Car and Parrinello 1985). In its most widely used formulation, the Car-Parrinello method employs DFT as the electronic-structure method of choice. In... 

For purposes of discussion, we divide applications of CPL spectroscopy into three categories (1) efforts to develop reliable CPL "sector rules", (2) use of comparative CD and CPL studies to probe excited state geometry changes, and (3) the specific use of the selectivity and sensitivity of CPL to probe details of molecular and electronic structure, and dynamics. Since in this book we are primarily concerned with "analytical" applications of these chiroptical methods, we will emphasize here the last of these categories. 

The main purpose of this chapter is to present the basics of ab initio molecular dynamics, focusing on the practical aspects of the simulations, and in particular, on modeling chemical reactions. Although CP-MD is a general molecular dynamics scheme which potentially can be applied in combination with any electronic structure method, the Car-Parinello MD is usually implemented within the framework of density functional theory with plane-waves as the basis set. Such an approach is conceptually quite distant from the commonly applied static approaches of quantum-chemistry with atom-centered basis sets. Therefore, a main... 

Quantum dynamics effects for hydride transfer in enzyme catalysis have been analyzed by Alhambra et. al., 2000. This process is simulated using canonically variational transition-states for overbarrier dynamics and optimized multidimensional paths for tunneling. A system is divided into a primary zone (substrate-enzyme-coenzyme), which is embedded in a secondary zone (substrate-enzyme-coenzyme-solvent). The potential energy surface of the first zone is treated by quantum mechanical electronic structure methods, and protein, coenzyme, and solvent atoms by molecular mechanical force fields. The theory allows the calculation of Schaad-Swain exponents for primary (aprim) and secondary (asec) KIE... 

Quantum chemical methods are valuable tools for studying atmospheric nucle-ation phenomena. Molecular geometries and binding energies computed using electronic structure methods can be used to determine potential parameters for classical molecular dynamic simulations, which in turn provide information on the dynamics and qualitative energetics of nucleation processes. Quantum chemistry calculations can also be used to obtain accurate and reliable information on the fundamental chemical and physical properties of molecular systems relevant to nucleation. Successful atmospheric applications include investigations on the hydration of sulfuric acid and the role of ammonia, sulfur trioxide and/or ions... 

Unfortunately, quantitatively reliable quantum chemical calculations of nucleation rates for atmospherically relevant systems would require the application of both high-level electronic structure methods and complicated anharmonic thermochemical analysis to large cluster structures. Such computations are therefore computationally too expensive for currently available computer systems, and will likely remain so for the foreseeable future. Instead, a synthesis of different approaches will probably be necessary. In the future, successful nucleation studies are likely to contain combinations of the best features of both classical (Monte Carlo and molecular dynamics) and quantum chemical methods, with the ultimate objective being a chemically accurate, complete configurational sampling. 

The introduction of the Cm-Parrinello method has not only extended the range of classical MD simulations based on empirical potentials but at the same time, it has also significantly increased the capabilities of conventional electronic structure calculations. Through the combination with a MD method a generalization to finite temperature and condensed phase systems was achieved. Furthermore, a whole set of simulation tools based on statistical mechanics can be apphed in this way in the context of an electronic structure method. Consequently, many dynamic as well as thermodynamic properties can be described within the accuracy of a first-principles method. 

---