Potential energy technique

The above discussion represents a necessarily brief simnnary of the aspects of chemical reaction dynamics. The theoretical focus of tliis field is concerned with the development of accurate potential energy surfaces and the calculation of scattering dynamics on these surfaces. Experimentally, much effort has been devoted to developing complementary asymptotic techniques for product characterization and frequency- and time-resolved teclmiques to study transition-state spectroscopy and dynamics. It is instructive to see what can be accomplished with all of these capabilities. Of all the benclunark reactions mentioned in section A3.7.2. the reaction F + H2 —> HE + H represents the best example of how theory and experiment can converge to yield a fairly complete picture of the dynamics of a chemical reaction. Thus, the remainder of this chapter focuses on this reaction as a case study in reaction dynamics. 

Molecular dynamics consists of the brute-force solution of Newton s equations of motion. It is necessary to encode in the program the potential energy and force law of interaction between molecules the equations of motion are solved numerically, by finite difference techniques. The system evolution corresponds closely to what happens in real life and allows us to calculate dynamical properties, as well as thennodynamic and structural fiinctions. For a range of molecular models, packaged routines are available, either connnercially or tlirough the academic conmuinity. 

In this chapter, we look at the techniques known as direct, or on-the-fly, molecular dynamics and their application to non-adiabatic processes in photochemistry. In contrast to standard techniques that require a predefined potential energy surface (PES) over which the nuclei move, the PES is provided here by explicit evaluation of the electronic wave function for the states of interest. This makes the method very general and powerful, particularly for the study of polyatomic systems where the calculation of a multidimensional potential function is an impossible task. For a recent review of standard non-adiabatic dynamics methods using analytical PES functions see [1]. 

In the following, we shall demonstrate techniques for calculating the electronic potential energy terms up to the second order. For simplicity, we shall study the case of H2 molecule, the simplest multi-electron diatomic molecule. 

The following two techniques were developed to expand such static calculations into a pseudo-dynamic regime by calculating higher derivatives of the potential energy and by introducing an additional degree of freedom. 

For a given potential energy function, one may take a variety of approaches to study the dynamics of macromolecules. The most exact and detailed information is provided by MD simulations in which one solves the equations of motion for the atoms constituting the macromolecule and any surrounding environment. With currently available techniques and methods it is possible... 

Results using this technique are better for force helds made to describe geometries away from equilibrium. For example, it is better to use Morse potentials than harmonic potentials to describe bond stretching. Some researchers have created force helds for a specihc reaction. These are made by htting to the potential energy surface obtained from ah initio calculations. This is useful for examining dynamics on the surface, but it is much more work than simply using ah initio methods to hnd a transition structure. 

The reaction coordinate is one specific path along the complete potential energy surface associated with the nuclear positions. It is possible to do a series of calculations representing a grid of points on the potential energy surface. The saddle point can then be found by inspection or more accurately by using mathematical techniques to interpolate between the grid points. 

This type of calculation does reliably find a transition structure. However, it requires far more computer time than any of the other techniques. As such, this is generally only done when the research requires obtaining a potential energy surface for reasons other than just finding the transition structure. 

If the system can only be modeled feasibly by molecular mechanics, use the potential energy curve-crossing technique or a force held with transition-structure atom types. 

Rather than using transition state theory or trajectory calculations, it is possible to use a statistical description of reactions to compute the rate constant. There are a number of techniques that can be considered variants of the statistical adiabatic channel model (SACM). This is, in essence, the examination of many possible reaction paths, none of which would necessarily be seen in a trajectory calculation. By examining paths that are easier to determine than the trajectory path and giving them statistical weights, the whole potential energy surface is accounted for and the rate constant can be computed. 

If we require similar information regarding the ground state potential energy surface in a polyatomic molecule the electronic emission specttum may again provide valuable information SVLF spectroscopy is a particularly powerful technique for providing it. 

Computer simulation techniques offer the ability to study the potential energy surfaces of chemical reactions to a high degree of quantitative accuracy [4]. Theoretical studies of chemical reactions in the gas phase are a major field and can provide detailed insights into a variety of processes of fundamental interest in atmospheric and combustion chemistry. In the past decade theoretical methods were extended to the study of reaction processes in mesoscopic systems such as enzymatic reactions in solution, albeit to a more approximate level than the most accurate gas-phase studies. 

Chapter 3, Geometry Optimizations, describes how to locate equilibrium structures of molecules, or, more technically, stationary points on the potential energy surface. It includes an overview of the various commonly used optimization techniques and a consideration of optimizing transition strucmres as well as minimizations. 

Though we and others (27-29) have demonstrated the utility and the improved sensitivity of the peroxyoxalate chemiluminescence method for analyte detection in RP-HPLC separations for appropriate substrates, a substantial area for Improvement and refinement of the technique remains. We have shown that the reactions of hydrogen peroxide and oxalate esters yield a very complex array of reactive intermediates, some of which activate the fluorophor to its fluorescent state. The mechanism for the ester reaction as well as the process for conversion of the chemical potential energy into electronic (excited state) energy remain to be detailed. Finally, the refinement of the technique for routine application of this sensitive method, including the optimization of the effi-ciencies for each of the contributing factors, is currently a major effort in the Center for Bioanalytical Research. 

---