Potential energy surface molecular geometry

By taking as a reference the calculation in vacuo, the presence of the solvent introduces several complications. In fact, besides the direct effect of the solvent on the solute electronic distribution (which implies changes in the solute properties, i.e. dipole moment, polarizability and higher order responses), it should be taken into account that indirect solvent effects exist, i.e. the solvent reaction field perturbs the molecular potential energy surface (PES). This implies that the molecular geometry of the solute (the PES minima) and vibrational frequencies (the PES curvature around minima in the harmonic approximation) are affected by the presence of a solvating environment. Also, the dynamics of the solvent molecules around the solute (the so-called nonequilibrium effect ) has to be... 

This expansion of the total electronic wavefunction is very compact, and provides a great deal of physical and chemical visuality. The spin-coupled structure (18) by itself reproduces with very reasonable accuracy all the features of a ground-state molecular potential energy surface. For example the spin-coupled function typically yields 85% of the observed binding energy, and equilibrium internuclear separations are accurate to 0.01 A. This function consequently dominates expansion (17) for all nuclear geometries. The various excited structures provide angular and other types of correlation as an extra quantitative refinement but do not alter the qualitative picture. 

For molecules of chemical interest it is not possible to calculate an exact many-electron wave function. As a result, we have to make certain approximations. The most commonly made approximation is the molecular orbital approximation, which is outlined in the next section. Within such a framework, it is useful to define various levels of computational method, each of which can be applied to give a unique wave function and energy for any set of nuclear positions and number of electrons. If such a model is clearly specified and if it is sufficiently simple to apply repeatedly, it can be used to generate molecular potential energy surfaces, equilibrium geometries, and other physical properties. Each such theoretical model can then be explored and the results compared in detail with experiment. If there is sufficient consistent success, some confidence can then be acquired in its predictive power. Each such level of theory therefore should be thoroughly tested and characterized before the significance of its prediction is assessed. 

At present, reaction path methods represent the best approach for utilizing ab initio electronic structure theory directly in chemical reaction dynamics. To study reaction dynamics we need to evaluate accurately the Born-Oppenheimer molecular potential energy surface. Our experience suggests that chemical reaction may take place within in a restricted range of molecular configurations (i.e., there is a defined mechanism for the reaction). Hence we may not need to know the PES everywhere. Reaction path methods provide a means of evaluating the PES for the most relevant molecular geometries and in a form that we can use directly in dynamical calculations. 

Geometry optimization of the proposed mimetic is included as part of the design analysis to ensure the feasibility of the desired molecular conformation. MM and semiempirical quantum mechanical methods have been used most extensively for these purposes. Conformational analysis of the proposed mimetic allows the determination of an energy profile for the molecule under consideration. This has been used by researchers to assess where the desired conformation for the mimetic resides on the molecular potential energy surface. Monte Carlo, MD, and distance geometry-based conformational search techniques have been employed extensively to sample conformational space. Computational methods that attempt to approximate the efifects of aqueous solvation on the conformational profile of the mimetic are being used more frequently as part of these efforts. 

Recently, Jasien k Shepard (11) proposed a new approach for representing molecular potential energy surfaces. They expand the energy (and derivatives, if available) about a line which connects the reactants with the products through the saddle point. This line can be a simple interpolation between the geometries of the above points or it can be an approximate reaction... 

To calculate the properties of a molecule, you need to generate a well-defined structure. A calculation often requires a structure that represents a minimum on a potential energy surface. HyperChem contains several geometry optimizers to do this. You can then calculate single point properties of a molecule or use the optimized structure as a starting point for subsequent calculations, such as molecular dynamics simulations. 

HyperChem provides three types of potential energy surface sampling algorithms. These are found in the HyperChem Compute menu Single Point, Geometry Optimization, and Molecular Dynamics. 

Geometry optimizations usually attempt to locate minima on the potential energy surface, thereby predicting equilibrium structures of molecular systems. Optimizations can also locate transition structures. However, in this chapter we will focus primarily on optimizing to minima. Optimizations to minima are also called minimizations. 

An IRC calculation examines the reaction path leading down from a transition structure on a potential energy surface. Such a calculation starts at the saddle point and follows the path in both directions from the transition state, optimizing the geometry of the molecular system at each point along the path. In this way, an IRC calculation definitively connects two minima on the potential energy surface by a path which passes through the transition state between them. 

There exist a series of beautiful spectroscopy experiments that have been carried out over a number of years in the Lineberger (1), Brauman (2), and Beauchamp (3) laboratories in which electronically stable negative molecular ions prepared in excited vibrational-rotational states are observed to eject their extra electron. For the anions considered in those experiments, it is unlikely that the anion and neutral-molecule potential energy surfaces undergo crossings at geometries accessed by their vibrational motions in these experiments, so it is believed that the mechanism of electron ejection must involve vibration-rotation... 

---