Potential energy surfaces determination

Figure 21. Family of reactive trajectories in the ground adiabatic potential energy surface determined by Eq. (13). Crosses indicate the caustics. Taken from Ref. [29].

Fig. 1. Schematic representation of the potential energy surface for the electronic (el) ground state of a molecule existing in two tautomeric forms, A and B. Superscripts exp, HF, CNDO/2, MINDO/3 indicate that energy differences 8 a,b calculated for potential energy surfaces determined either experimentally (exp) or calculated by means of ab initio method in the Hartree-Fock (HF) approximation or by semiempirical methods (CNDO/2, MINDO/3). The symbol eq stands for the geometrical equilibrium of both tautomers, while 2a and Qb indicate nonequilibrium geometries of tautomers A and B, respectively. Note that the theoretical potential surface calculated by sophisticated quantum-mechanical methods ( exact solution of electronic Schrbdinger equation includes electron correlation with geometry optimization) should be the same (or very similar) as that determined experimentally [in this case i>eor) ei

A much more common procedure for reference selection is to accept references whose estimated importance is greater than some given threshold this can involve perturbative estimates of a function s energetic contribution or its coefficient in some preliminary wavefunction. These approaches are more successful at obtaining the best wavefunction at the lowest expense, but they sacrifice the simplicity of the excitation class selection and can become more difficult to implement and to use. One complication is that potential energy surfaces determined using such methods may not be smooth to alleviate this, one may need to determine the important references at each geometry and use the union of these sets at every point. 

The molecule is a van der Waals molecule with no well defined conformation. The value given here corresponds to the minimum of the potential energy surface determined from the lit of IR spectra of Ar...H2 and Ar...D2. At this minimum, the molecule has C2v symmetiy. For further details, refer to the original paper. is the distance between the center of mass of the hydrogen molecule and the argon atom. 

Fig. 3 Potential energy surface determined as a function of an applied force for the electro-cyclic ring opening of cyclobutene. J. Ribas-Arino, M. Shiga and D. Marx. Angew. Chem., Int. Ed., 2009, 48, 4910-4913. Copyright Wiley-VCH Verlag GmbH Co. KGaA. Reproduced with permission.

The full dynamical treatment of electrons and nuclei together in a laboratory system of coordinates is computationally intensive and difficult. However, the availability of multiprocessor computers and detailed attention to the development of efficient software, such as ENDyne, which can be maintained and debugged continually when new features are added, make END a viable alternative among methods for the study of molecular processes. Eurthemiore, when the application of END is compared to the total effort of accurate determination of relevant potential energy surfaces and nonadiabatic coupling terms, faithful analytical fitting and interpolation of the common pointwise representation of surfaces and coupling terms, and the solution of the coupled dynamical equations in a suitable internal coordinates, the computational effort of END is competitive. 

The essence of this analysis involves being able to write each wavefunction as a combination of determinants each of which involves occupancy of particular spin-orbitals. Because different spin-orbitals interact differently with, for example, a colliding molecule, the various determinants will interact differently. These differences thus give rise to different interaction potential energy surfaces. 

If the complete potential energy surface has already been computed, a reaction coordinate can be determined using an adaptation of the IRC algorithm. The IRC computation requires very little computer time, but obtaining the potential energy surface is far more computation-intensive than an ah initio IRC calculation. Thus, this is only done when the potential energy surface is being computed for another reason. 

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. 

Another use of frequency calculations is to determine the nature of a stationary point found by a geometry optimization. As we ve noted, geometry optimizations converge to a structure on the potential energy surface where the forces on the system are essentially zero. The final structure may correspond to a minimum on the potential energy surface, or it may represent a saddle point, which is a minimum with respect to some directions on the surface and a maximum in one or more others. First order saddle points—which are a maximum in exactly one direction and a minimum in all other orthogonal directions—correspond to transition state structures linking two minima. 

You can plot the results of the scan to get a picture of the region of the potential energy surface that you ve explored. By doing so, you may be able to determine the approximate location of the minimum energy structure. However, potential energy surfece scans do not include a geometry optimization. 

Maxima, minima and saddle points are stationary points on a potential energy surface characterized by a zero gradient. A (first-order) saddle point is a maximum along just one direction and in general this direction is not known in advance. It must therefore be determined during the course of the optimization. Numerous algorithms have been proposed, and I will finish this chapter by describing a few of the more popular ones. 

There are many facets to a successful quantum dynamics study. Of course, if comparison with experimental results is a goal, the underlying Bom-Oppenheimer potential energy surface must be known at an appropriately high level of electronic structure theory. For nonadiabatic problems, two or more surfaces and their couplings must be determined. The present chapter, however, focuses on the quantum dynamics of the nuclei once an adequate description of the electronic structure has been achieved. 

The method is composed of the following algorithms (1) transition position is detected along each classical trajectory, (2) direction of transition is determined there and the ID cut of the potential energy surfaces is made along that direction, (3) judgment is made whether the transition is LZ type or nonadiabatic tunneling type, and (4) the transition probability is calculated by the appropriate ZN formula. The transition position can be simply found by... 

---