Variational optimization

Since molecular HF equations cannot be solved numerically each m.o. is expanded as an LCAO and the expansion coefficients optimized variationally. The equations in the form... 

In the optimized orbitals CCD (00-CCD or OD) model, the orbitals are optimized variationally to minimize the total energy of the 00-CCD wavefunction. This allows one to drop single excitations from the wavefunction. Conceptually, 00-CCD is very similar to the Brueckner CCD (B-CCD) method. Both 00-CCD and B-CCD perform similarly to CCSD in most cases. 

Within the LCAO-MO approximation, the process for finding the MOs that minimize the total HF energy becomes clearer. If the mathematical form of each of the AOs, ( )r, in a molecule is assumed to be known, it is the coefficient, Cjj, of each AO in each filled MO that is optimized variationally. This is accomplished by solving Eq. 3 untU self-consistency is reached. 

This variability in AO sizes is most easily introduced by using more than one independent mathematical function to represent each AO. The contribution of each function to each AO in each MO can then be optimized variationally, by performing a HF calculation. The set of mathematical functions used to represent each of the AOs in a molecule is called the basis set. 

In an interesting application of a gaseous reagent to solid-phase synthesis, Takahashi and co-workers demonstrated hydroformylation of an unactivated alkene using synthetic gas (1 1 H2-CO) and a Rh(I) catalyst (Scheme 8).22 The reaction was typically performed at 40°C in toluene at a pressure of 75 atm. Conversions of 99% were obtained following careful reaction optimization. Variation in the concentration of catalyst could be used to alter the regioselectivity of the reaction. 

Let us consider the first-order properties for the optimized variational energy E (x) in Eq. 8. Using the chain rule, we obtain... 

The HF method treats electron-electron interactions at a mean field level, with the Hartree and exchange interactions exactly written. The method can be implemented either in its spin restricted form (RHF), for closed shell systems, or in the unrestricted form (UHF) for open-shell or strongly correlated systems. In the first case, the one-electron orbitals are identical for electrons of both spin directions, while UHF can account for a non-uniform spin density. The one electron orbitals, which are determined in the course of the self-consistent resolution of the HF equations, are expanded on an over-complete basis set of optimized variational functions. 

Using the same form of an optimized variational wavefunction... 

Numerical values of the interaction energies for these Heitler-London wavefunctions, taken from Magnasco (2008), are given in Tables 1.2 and 1.3. The energies are optimized variationally with respect to the values of the orbital exponents Co of the atomic Is STOs on A and B. 

In ofher words, for afoms and diafomics, in fhe SPSA compufafional mefhodology fhaf was firsf implemenfed in fhe 1970s, a mixfure of numerical and analyfic orbifals is used fo improve convergence and accuracy, for example, [20, 26a-26c]. Specifically, on fhe one hand, fhe accuracy of is secured in ferms of fhe radial defails. On fhe ofher hand, fhe overall calculation, which involves flo P and parts of remains flexible by using analyfic virtual orbitals that are optimized variationally. 

The first is the fact that they facilitate the computations in obtaining converged MCHF solutions for lower lying mixed wavefunctions. In fact, being optimized variationally self-consistently, they absorb some of the contribution from other excited configurations with the same structure. Obviously,... 

Since the width factor Q is in this case the dynamical variable, the centroid frequency can be obtained from Q through the one-dimensional versions of Eqs. (3.74) and (3.75). The optimal variational value for Q is obtained by setting Eq. (3.80) equal to zero and self-consistently solving for Q. For a given centroid initial condition, this value of Q would provide the initial condition for the variable Q in the extended Lagrangian simulation. 

These questions will be discussed separately, starting with the first one, and showing that, within the conventional CO formulation of the magnetic response [58,67], exact cancellation takes place only in the case of optimal variational wavefunctions which satisfy hypervirial constraints [65]. 

To understand the difference between the spurious gauge terms (i), which are removed in the ideal case of optimal variational wavefunctions [65], and terms (ii), which account for the essential origin dependence of the property, let us first discuss die origin dependence of the quadrupole polarizability of magnetic susceptibility (22) within the conventional common-origin representation. 

A more elaborate approach, which involves global functional minimization of the rate constant, is offered by VTST [136, 317]. For outer-sphere ET with a harmonic bath, the reaction coordinate is represented as a linear combination of bath coordinates. The coefficients are optimized variation-ally to minimize the one-way reactive flux at the crossing point of two diabatic surfaces. A clear explanation of this optimization procedure has been presented recently by Benjamin and Poliak [136c]. The basic idea of... 

The exact definition is slightly more complicated, since the wave function has to be properly antisymmetrized and projected onto the actual basis but, for illustration, the above form is sufficient. Such R12 wave functions may then be used in connection with the Cl, MBPT or CC methods described above. Consider for example a Cl calculation with an R12 type wave function. The energy is given by eq. (4.78), where the ayab and bij parameters in (4.77) are optimized variationally. 

---