Points stationary

Mathematically, a stationary point is one at which the first derivative of the potential energy with respect to each geometric parameter is zero  

Inspection of Fig. 2.7 shows that the transition state linking the two minima represents amaximum along the direction ofthe IRC, but along all other directions it is a minimum. This is a characteristic of a saddle-shaped surface, and the transition state is called a saddle point (Fig. 2.8). The saddle point lies at the center of the saddle-shaped region and is, like a minimum, a stationary point, since the PES at that point is parallel to the plane defined by the geometry parameter axes we can see that a marble placed (precisely) there will balance. Mathematically, minima and saddle points differ in that although both are stationary points (they have zero first derivatives Eq. (2.1)), a minimum is a minimum in all directions, but a saddle point is a maximum along the 

The geometric parameter corresponding to the reaction coordinate is usually a composite of several parameters (bond lengths, angles and dihedrals), although for some reactions one or two may predominate. In Fig. 2.7, the reaction coordinate is a composite of the 0-0 bond length and the 0-0-0 bond angle. 

The Concept of the Potential Energy Surface Newman projections 

There is a computational chemistry equivalent to the experimentalists purification of a sample before undertaking any measurements the calculated geometry must be at a stationary point in geometric space. For a ground-state system, this means the geometry is such that the calculated AHf is an irreducible minimum for transition states the sum of all forces acting on each atom in the system must be zero, and the system must have exactly one negative force constant. 

In general, it is possible to refine all geometries so that all calculated are very precise, say within 0.001 kcal/mol. This is, however, normally neither necessary nor desirable. Refining geometries to such a degree is costly. Heats of reaction and barriers to activation should not be considered more accurate than 5 kcal/mol, so working to a precision of 0.01 kcal/mol is likely to result in a false sense of accuracy. 

In terms of the residue curve equation, this means algebraically solving for points where x = y. Generally speaking, these stationary points may be classified into three main types a stable node, an unstable node, and a saddle point, depicted in 

FKiURE 2 7 Topological representation of (a) a stable node, (b) an unstable node, and (c) a saddle point. 

FU URE 2.8 RCM plots generated in Aspen Plus S for the acetone/edianol/medianol system using die NRTL thermodynamic model at (a) P = 1 atm and (b) P= 10 atm. 

So far it has been shown that the derivative gives the gradient of the tangent to a specified curve at any point. At the turning or stationary points, also known as maxima and minima, the gradient of the tangent will be zero. Consequently, if the points at which the derivative is zero can be determined, the turning points of the curve will be known. 

Determining the nature of the stationary point is done by looking at the value of 2- if it is negative the stationary point is a maximum whereas if it is positive the stationary point is a minimum. If the second derivative is zero the stationary point may be a maximum, a minumum, or a point of inflexion as shown in Fig. 37.1. Its nature in this case can only be determined by considering the sign of the gradient each side of the stationary point. 

Since this is positive, the stationary point at x = 3/2 is a miniumum. 

Removing the common factor 3 allows this to be rewritten using brackets as 

Note that it is not always necessary to obtain the final numerical value of in order to determine 

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. 

There are two pieces of information from the output which are critical to characterizing a stationary point  

Imaginary frequencies are listed in the output of a frequency calculation as negative numbers. By definition, a structure which has n imaginary frequencies is an n order saddle point. Thus, ordinary transition structures are usually characterized by one imaginary frequency since they are first-order saddle points. 

If applicable, the program notes that there is an imaginary frequency present just prior to the frequency and normal modes output, and the first frequency value is less than zero. Log files may be searched for this line as a quick check for imaginary frequencies. 

Exploring Chemistry with Electronic Structure Methods 

Potential energy surfaces are important because they aid us in visualizing and understanding the relationship between potential energy and molecular geometry, and in understanding how computational chemistry programs locate and characterize structures 

The distinction is sometimes made between a transition state and a transition structure [4]. Strictly speaking, a transition state is a thermodynamic concept, the species an ensemble of which are in a kind of equilibrium with the reactants in Eyring s2 transition-state theory [5]. Since equilibrium constants are determined by free energy differences, the transition structure, within the strict use of the term, is a free energy maximum along the reaction coordinate (in so far as a single species can 

2Henry Eyring, American chemist. Bom Colonia Juararez, Mexico, 1901. Ph.D. University of California, Berkeley, 1927. Professor Princeton, University of Utah. Known for his work on the theory of reaction rates and on potential energy surfaces. Died Salt Lake City, Utah, 1981. 

The chemically most interesting arrangements of nuclei correspond to 

each stationary point of E is a solution of the nonlinear 

When a gradient vector g(x) is decomposed according to Eg. (1) we get the forces that act on the nuclei of the molecular system. In these terms. Eg.(3) means that the forces (the force vector ) vanish(es) at stationary points. The solution set of Eg.(4) is denoted by g (0), 

This notational convention will also be employed for other functions. 

A stationary point can be classified by means of the points of its neighborhoods. 

---

As is a stationary point, the linear temis FlEZ- h/j ) q" i/iin equation B3.5.1 vanish, and higher-order temis do... 

Let us express the displacement coordinates as linear combinations of a set of new coordinates y >q= Uy then AE = y U HUy. U can be an arbitrary non-singular matrix, and thus can be chosen to diagonalize the synmietric matrix H U HU = A, where the diagonal matrix A contains the (real) eigenvalues of H. In this fomi, the energy change from the stationary point is simply AF. = t Uj A 7- h is clear now that a sufBcient... 

The simplest smooth fiuictioii which has a local miiiimum is a quadratic. Such a function has only one, easily detemiinable stationary point. It is thus not surprising that most optimization methods try to model the unknown fiuictioii with a local quadratic approximation, in the fomi of equation (B3.5.1). 

These methods, which probably deserve more attention than they have received to date, simultaneously optimize the positions of a number of points along the reaction path. The method of Elber and Karpins [91] was developed to find transition states. It fiimishes, however, an approximation to the reaction path. In this method, a number (typically 10-20) equidistant points are chosen along an approximate reaction path coimecting two stationary points a and b, and the average of their energies is minimized under the constraint that their spacing remains equal. This is obviously a numerical quadrature of the integral s f ( (.v)where... 

Banerjee A, Adams N, Simons J and Shepard R 1985 Search for stationary points on surfaces J. Phys. Chem. 89 52... 

Conical intersections, introduced over 60 years ago as possible efficient funnels connecting different elecbonically excited states [1], are now generally believed to be involved in many photochemical reactions. Direct laboratory observation of these subsurfaces on the potential surfaces of polyatomic molecules is difficult, since they are not stationary points . The system is expected to pass through them veiy rapidly, as the transition from one electronic state to another at the conical intersection is very rapid. Their presence is sunnised from the following data [2-5] ... 

The origin of a torsional barrier can be studied best in simple cases like ethane. Here, rotation about the central carbon-carbon bond results in three staggered and three eclipsed stationary points on the potential energy surface, at least when symmetry considerations are not taken into account. Quantum mechanically, the barrier of rotation is explained by anti-bonding interactions between the hydrogens attached to different carbon atoms. These interactions are small when the conformation of ethane is staggered, and reach a maximum value when the molecule approaches an eclipsed geometry. 

Example I hc reaction coordinate for rotation about the central carbon-carbon bond in rt-bulane has several stationary points.. A, C, H, and G are m in im a and H, D, an d F arc tn axirn a. Only the structures at the m in im a represen t stable species an d of these, the art/[ con form ation is more stable th an ihc nauchc. 

HyperChem performs ti vibrational analysisat the molecular geometry shown m the IlyperChem workspace, without any automatic pre-optini i/ation. IlyperChem may thus give unreasonable results when yon perform vibrational analysiscalcnlations woth an nnoptimized molecular system, particularly for one far from optimized. Because the molecular system is not at a stationary point, neither at a local minimum nor at a local maximum, the vibra-... 

Dmbining this with the constraint equation enables us to identify the stationary point, hich is at (-59/72, -23/18). 

The problem of a mass suspended by a spring from another mass suspended by another spring, attached to a stationary point (Kreyszig, 1989, p. 159ff) yields the matrix equation... 

Reading the output for H2 is similar to Hj as well. The optimized bond distance is — stationary point found. 

Use a forced convergence method. Give the calculation an extra thousand iterations or more along with this. The wave function obtained by these methods should be tested to make sure it is a minimum and not just a stationary point. This is called a stability test. 

As mentioned earlier, a potential energy surface may contain saddle points , that is, stationary points where there are one or more directions in which the energy is at a maximum. Asaddle point with one negative eigenvalue corresponds to a transition structure for a chemical reaction of changing isomeric form. Transition structures also exist for reactions involving separated species, for example, in a bimolecular reaction... 

---