States Optimal

Gaussian includes a facility for automatically generating a starting structure for a transition state optimization based upon the reactants and products that the transition structure connects, known as the STQN method. This feature is requested with the QST2 option to the Opt keyword. Input files using this option will include two title and molecule specification sections. The facility generates a guess for the transition structure which is midway between the reactants and products, in terms of redundant internal coordinates. 

When this initial guess is poor, you need a more sophisticated—albeit more expensive—means of generating the force constants. This is especially important for transition state optimizations. Gaussian provides a variety of alternate ways of generating them. Here are some of the most useful associated keywords consult the Gaussian User s Reference for a full description of their use ... 

Vinyl alcohol is a good system for discussing transition state optimizations. The (iptimizatior. to a transition state starting from the 9Sf form proceeds easily, making it a suitable introduction to the general topic. 

The transition state optimization (Opt=(TS,CakFC)) of the structure on the right converges in 12 steps. The UHF frequency calculation finds one imaginary frequency. Here is the associated normal mode ... 

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. 

For our initial geometry for the transition structure, we ll detach one hydrogen from the carbon and increase the O-C-H bond angle. We specified the Opt=(TS,CalcFC) keyword in the route section, requesting an optimization to a transition state. The CalcFC option is used to compute the initial force constants, a technique which is generally helpful for transition state optimizations. We ve also included the Freq keyword so that a frequency calculation will automatically be run at the optimized geometry. 

Latour, P., Use of steady-state optimization for computer control in the process industries. In On-line Optimization Techniques in Industrial Control (Kompass, E. J. and Williams, T. J., eds.). Technical Publishing Company, 1979. 

Figure 2-5. Geometries of the 02-bound state optimized using the active-site model (left) and an ONIOM model (right). Note the large differences in geometry of the two calculations, especially the hydrogen bonds donated to O2 in the ONIOM model (marked in grey) (Adapted from Hoffman et al. [25]. Reprinted with permission. Copyright 2004 Wiley Periodicals, Inc.)...

Timm, Gilbert, Ko, and Simmons O) presented a dynamic model for an isothermal, continuous, well-mixed polystyrene reactor. This model was in turn based upon the kinetic model developed by Timm and co-workers (2-4) based on steady state data. The process was simulated using the model and a simple steady state optimization and decoupling algorithm was tested. The results showed that steady state decoupling was adequate for molecular weight control, but not for the control of production rate. In the latter case the transient fluctuations were excessive. 

A structural comparison of the calculated (B3LYP/6-311+G ) ts (transition state in the gas phase), ts-wc (transition state in the cluster of five extra water molecules), ts-CPCM (transition state within the CPCM-solvent model (B3LYP(CPCM)/6-311+G )) and ts-PCM (transition state optimized within the PCM-solvent model (B3LYP(PCM)/6-311+G )), shows no large differences (see Fig. 8), which is also valid for the precursor complexes (see Fig. 9). Modeling solvent effects shrinks in all cases the Be-0 bonds of the entering/leaving water molecules (159). 

The targets for the MPC calculations are generated by solving a steady-state optimization problem (LP or QP) based on a linear process model, which also finds the best path to achieve the new targets (Backx et al., 2000). These calculations may be performed as often as the MPC calculations. The targets and constraints for the LP or QP optimization can be generated from a nonlinear process model using a nonlinear optimization technique. If the optimum occurs at a vertex of constraints and the objective function is convex, successive updates of a linearized model will find the same optimum as the nonlinear model. These calculations tend to be performed less frequently (e.g., every 1-24 h) due to the complexity of the calculations and the process models. 

Optimal rotors - the rotors that are both best suited to perform the run and to achieve the stated optimization criterion,... 

Figure 5.21 Transition state optimized by ab initio DFT/B3LYP calculations for hydride-hydride exchange via the turnstile mechanism. (Reproduced with permission from ref. 32.)...

In scheme (6) the basis set is optimized by invoking the variation principle for each state considered. For the ground state the optimized values of the even-tempered parameters ao and / o given by Schmidt and Ruedenberg [9] are used. We add the subscript 0 to distinguish ground state values. For the excited state optimal ai and / i values for a sequence of Mi values are determined. 

Fig. 2. Relative energy differences between the spin states of [Fe2] 1, for each spin state optimized structures denoted opt as well as for single-point calculations on the optimized high-spin structure denoted as single point, where the energy of the closed-shell spin restricted state has arbitrarily been set to zero.

---