VBTutorial2

<< Return to the program

How to modify this page :

• first : log in (top right of this page) ;
• click on [edit] (far right) to edit a section of the page ;
• write your text directly in the wiki page, and click on the "Save page" button (bottom left) to save your modifications

Pictures : how to insert a picture in your text

See also this page for an introduction to the basics of the wiki syntax

VB applications on PI systems

Remarks

>> Initial remarks for tutors

Exercices

Remark : in all the following exercises, <math>\pi</math> the system will be taken as active, and the <math>\sigma</math> system as inactive. In all VB calculations, the <math>\sigma</math> orbitals shall be described by MOs delocalized onto the whole molecule.

Exercise 1 : The allyl radical and cation

***** INPUT FILES TO BE WRITTEN *****

1. Structures of the allyl radical (paper exercise) :
   1. What are the three possible covalent structures for the allyl radical molecule ? Use qualitative VB theory to compute their energy, and show that two of them are degenerate.
   2. Show that the third structure can be expressed as a linear combination of the first two structures, and thus that only two of the three covalent structures form a complete basis of non-redundant structures (Rumer basis).
   3. What are the possible ionic structures for the allyl radical ? Based on your chemical knowledge, propose a selection subset of the most chemically meaningful covalent + ionic structures.
2. Structures of the allyl cation, and weights from Hückel theory (paper exercise) :
   1. Find a complete basis of covalent structures for the allyl cation.
   2. What are the possible ionic structures for the allyl cation ? Based on your chemical knowledge, propose a selection subset of the most chemically meaningful covalent + ionic structures.
   3. Use the HuLis software to retrieve the Hückel MOs for allyl cation. Write a single-determinant MO wave function based on Hückel orbitals. Develop it into the basis of atomic orbitals, to get an expression in terms of VB structures. Compute by hand the weights of the different structures (neglecting all overlaps for simplicity).
3. Computation of allyl radical and cation weights :
   1. Compute a VBSCF wave function for allyl radical (6-31G* basis set) using your selected set of structures (questions 1.3 and 2.2). Compute a VBSCF wave function including the complete set of VB structures, using the "str=all" keyword. Compare the weights and energies for both wave functions to validate your selection of structures.
   2. Compute for allyl radical a BOVB wave function which include only your selected set of most chemically meaningful structures (use the guess orbitals obtained at the VBSCF level). Compare the weights obtained at the VBSCF and BOVB level.
   3. Repeat questions 3.1 and 3.2 for the allyl ion. Compare the BOVB weights to the weights predicted by simple MO theory.
4. Computation of resonance energies :
   1. We want to build a wave-function corresponding to only one Lewis structures for the allyl radical. To do so, we will include in the wave-function only one covalent structure, and the ionic structures associated to this covalent bond. Propose a selection of VB structures which would describe one Lewis structure for the allyl radical.
   2. Compute this wave function at VBSCF then BOVB levels. Deduce what is the resonance energy of the allyl radical at the BOVB level.
   3. Repeat questions 4.1 and 4.2 for the allyl cation.
5. Stability of the diabetic quantities with the basis set :
   1. Recompute the resonance energy of allyl
6. I would just add here an explicit expression of of the covalent structure: 1/6^1/2[ |abc| + |abc| - 2|abc|] and showing the origins of spin polarization (see Book p 198, and Ex. 7.1. In fact, it is good to stretch to Ex 7.2 by jusy drawing the spin alternant determinant and predicting the patterns of spin polrization)

>> Answer

Exercise 2 : Radical character of ozone

1. Paper exercice :
   1. Propose a complete basis of non-redondant VB structures for the ozone molecule. Based on your chemical knowledge, propose a selection subset of the most chemically meaningful structures.
   2. Use the HuLis software to retrieve the Hückel MOs for ozone. Write a single-determinant MO wave function based on Hückel orbitals. Develop it into the basis of atomic orbitals, to get an expression in terms of VB structures.
   3. Compute by hand the weights of the different structures (neglecting all overlaps for simplicity). What is the radical character of ozone according to simple MO theory ?
2. Computer exercise :
   1. Compute a VBSCF wave function for ozone (6-31G* basis set) using your selected set of structures (question 1.1). Compute a VBSCF wave function including the complete set of VB structures, using the "str=all" keyword. Compare the weights and energies for both wave functions to validate your selection of structures.
   2. Compute a BOVB wave function for ozone which include only your selected set of most chemically meaningful structures (use the guess orbitals obtained at the VBSCF level).
   3. Compare the weights obtained at the VBSCF and BOVB level. Compare the BOVB diradical weight to the value predicted by simple MO theory.

>> Answer

Exercise 3 : Resonance energy of Benzene

***** EXERCISE AND INPUT FILES TO BE WRITTEN *****

ideas :

• Covalent only description :
  • Ask to compute the benzene with str=cov (5 covalent structures),
  • then 1 kékulé strucure and the VRE => too low (reference : 90kcal/mol)
• Full description :
  • Ask to compute the benzene with str=full (175 structures)
  • Ask to compute one corresponding Kékulé (27 structures which have to be hand-written !) => good VRE, but cumbersome... We'll see in tutorial 4 that it is more simple with BLW
• Delocalization problem :
  • Ask people to compute 5 covalent struc. (keyword : "str=cov") first using strictly localized orbitals (VBSCF), then using OEO ('orbtyp=oeo).
  • Compute then 1 kékulé with OEO, and decuce the VRE. It will be much too low, as using OEO implicitly include other struc. than the kékulé, i.e. the individual structures loose their identity => OEO shall not be used to compute diabatic quantities, and more generally if active orbitals are too extended (too large basis set with diffuse functions for instance) diabatic quantities will be wrong (underestimated).

>> general guidelines for BOVB calculations