VBTutorial3

<< Return to the program

Valence Bond State correlation diagrams

Exercise 1 : Computation of state correlation Diagrams for a 3 centers / 4 electrons system

In this exercise the <math>\textrm{S}_{\textrm{N}}2</math> reaction Cl<math>{}^{-}</math> + CH3Cl -> ClCH3 + Cl<math>{}^{-}</math> will be studied in both vacuum and solution. Valence Bond State Correlation Diagrams (VBSCD) will be constructed at <math>\pi</math>-D-BOVB level. There are two parts in this exercise: basic part and optional part. The basic part is performed with MCP-DZP basis set in which the inner orbitals in Cl and C are described with MCP pseudo potential. The optional part is performed with 6-31+G* basis set, using the general specification for the xmvb input (expert users). Only reactant and transition state will be computed in this exercise, which is sufficient to build the VBSCD diagrams.

Note:How to perform a VBPCM calculation

A VBPCM calculation is performed in the similar way as the VB calculations in vacuum. One should prepare a GAMESS input file with solvent assigned such as: $PCM SOLVNT=WATER $END The details of PCM calculation in the GAMESS can be found in GAMESS manual. Keyword "VBTYP=XMVB" in CONTRL section is also essential. After the GAMESS input file is prepared, an XMI file with keyword "VBPCM" should be prepared with the same file name as GAMESS input file. In the current XMVB package, VBSCF/PCM and BOVB/PCM calculations are both supported.

Basic part

1. Compute the Energies and Wavefunctions at Reactant and Transition State with Different Sets of VB Structures Write the following VB structure sets for a 3 centers / 4 electrons system : all structures; minimal structures for reactant; minimal structures for product. Perform <math>\pi</math>-D-BOVB calculation for reactant (see "high symmetry cases" here): Perform all-structure <math>\pi</math>-D-BOVB calculation as following: Perform L-VBSCF calculation with "orbtyp=hao frgtyp=sao guess=mo", in which the orbitals are all localized on the Cl and CH<math>{}_3</math> groups; Perform π-D-VBSCF calculation where <math>\pi</math> orbitals are delocalized in the whole system and the <math>\sigma</math> orbitals are kept localized. Use the L-VBSCF orbitals as initial guess; Perform π-D-BOVB calculation with π-D-VBSCF orbitals as initial guess. Perform <math>\pi</math>-D-BOVB calculations with minimal structures for reactant and product with the same procedure as all-structure calculation. Perform <math>\pi</math>-D-BOVB calculation for transition state. The procedure is the same as step 2. Perform <math>\pi</math>-D-BOVB/PCM calculations for reactant: Perform all-structure <math>\pi</math>-D-BOVB/PCM calculation for the reactant in following steps: Perform L-VBSCF/PCM calculation with L-VBSCF orbitals as initial guess; Perform π-D-VBSCF/PCM calculation with L-VBSCF/PCM orbitals as initial guess; Perform π-D-BOVB/PCM calculation with π-D-VBSCF/PCM orbitals as initial guess. Perform <math>\pi</math>-D-BOVB/PCM calculations with minimal structures for reactant and product in the sames steps as all structure calculation. Perform BOVB/PCM calculations for transition state with the same procedure as step 4. 2. Analysis: Wavefunctions and Energies Compare the weights of structures at both reactant and transition state points, find the differences. Compute the Barrier height of the <math>\textrm{S}_{\textrm{N}}2</math> reaction in both vacuum and solution. See the difference of the barrier heights, and find out the reason. Compare the energies of reactant and product structures at reactant and transition state geometries, in both vacuum and solution. What's the difference of the energies at different points? Why? Compute the resonance energies at both reactant and transition state points, see the difference of the resonance energies. Answer VB Structures in the Computations= Total VB Structure Set VB Structure Set of The Reactant VB Structure Set of The Product Weights of Structures Weights of structures at Reactant Geometry S1 S2 S3 S4 S5 S6 VBSCF 0.590 0.002 0.001 0.338 0.070 -0.000 BOVB 0.578 0.022 0.006 0.336 0.059 0.000 VBSCF/PCM 0.604 0.001 0.001 0.316 0.078 0.000 BOVB/PCM 0.596 0.015 0.003 0.318 0.068 0.000 Weights of structures at Transition State Geometry S1 S2 S3 S4 S5 S6 VBSCF 0.239 0.240 0.027 0.496 -0.001 -0.001 BOVB 0.228 0.228 0.041 0.494 0.004 0.004 VBSCF/PCM 0.221 0.221 0.022 0.538 -0.001 -0.001 BOVB/PCM 0.215 0.215 0.035 0.528 0.003 0.003 Barrier of the Reaction Energies (a.u.) and Barriers (kcal/mol) of <math>\textrm{S}_{\textrm{N}}2</math> Reaction VBSCF BOVB VBSCF/PCM BOVB/PCM Reactant -37.03384 -37.05387 -37.13644 -37.15552 Trasition State -36.98034 -37.02691 -37.06980 -37.11506 Barrier 33.6 16.9 41.8 25.4 Resonance Energ of Transition State Energies (a.u.) and Resonance Energies (kcal/mol) of Transition State VBSCF BOVB VBSCF/PCM BOVB/PCM All Structures -36.98034 -37.02691 -37.06980 -37.11506 Reactant -36.95812 -36.97494 -37.05258 -37.06909 Product -36.95812 -36.97494 -37.05258 -37.06910 Resonance Energy 13.9 32.6 10.8 28.8 Resonance Energ of Transition State Energies(a.u.) and Resonance Energies (kcal/mol) of Reactant VBSCF BOVB VBSCF/PCM BOVB/PCM All Structures -37.03384 -37.05387 -37.13644 -37.15552 Reactant -37.03354 -37.05085 -37.13626 -37.15356 Product -36.59694 -36.59694 -36.88408 -36.88750 Resonance Energy 0.2 1.9 0.1 1.2

Optional part (expert users)

In this part, calculations with BFI section are performed with 6-31+G* basis set, which is desired for the experienced users. The inner orbitals are frozen as HF orbitals in all VB calculations and the valence basis functions are reorganized to hybrid basis functions so that the <math>\sigma</math>, <math>\pi_x</math> and <math>\pi_y</math> spaces can be separated well. A D-BOVB calculation is performed in 2 steps: Perform a VBSCF calculation with <math>\pi</math> orbitals delocalized in the whole system and <math>\sigma</math> orbitals localized on the Cl and CH<math>{}_3</math> groups; Perform a BOVB calculation with VBSCF orbitals as initial guess. The VB calculations are the same as the calculations performed above. Try to understand the BFI section, perform the calculations and compare the differences of barrier heights, resonance energies and performances with and without $BFI. Answer Weights of Structures Weights of structures at Reactant Geometry S1 S2 S3 S4 S5 S6 VBSCF 0.632 0.003 0.001 0.299 0.064 0.000 BOVB 0.566 0.061 0.005 0.288 0.081 0.000 VBSCF/PCM 0.645 0.002 0.001 0.279 0.073 0.000 BOVB/PCM 0.599 -0.019 -0.007 0.337 0.089 0.000 Weights of structures at Transition State Geometry S1 S2 S3 S4 S5 S6 VBSCF 0.239 0.240 0.027 0.496 -0.001 -0.001 BOVB 0.263 0.263 0.043 0.416 0.007 0.007 VBSCF/PCM 0.220 0.220 0.020 0.542 -0.001 -0.001 BOVB/PCM 0.251 0.251 0.036 0.452 0.005 0.005 Barrier of the Reaction Energies (a.u.) and Barriers (kcal/mol) of <math>\textrm{S}_{\textrm{N}}2</math> Reaction VBSCF BOVB VBSCF/PCM BOVB/PCM Reactant -958.65818 -958.67599 -958.75774 -958.78203 Trasition State -958.62157 -958.65321 -958.71004 -958.73960 Barrier 23.0 14.3 29.9 26.6 Resonance Energ of Transition State Energies (a.u.) and Resonance Energies (kcal/mol) of Transition State VBSCF BOVB VBSCF/PCM BOVB/PCM All Structures -958.62157 -958.65321 -958.71004 -958.73960 Reactant -958.59188 -958.60816 -958.68547 -958.70142 Product -958.59188 -958.60816 -958.68547 -958.70142 Resonance Energy 18.6 28.3 15.4 23.9 Resonance Energ of Transition State Energies(a.u.) and Resonance Energies (kcal/mol) of Reactant VBSCF BOVB VBSCF/PCM BOVB/PCM All Structures -958.65818 -958.67599 -958.75774 -958.78203 Reactant -958.56752 -958.67272 -958.75734 -958.77250 Product -958.31845 -958.31845 -958.57917 -958.58246 Resonance Energy 0.4 2.1 0.3 6.0

>> general guidelines for BOVB calculations

Optional exercises - Homework

Exercise 2 : computation of H—H + H. -> H. + H—H radical exchange VBSCD diagram 1/ Paper exercise : a/ Considering the following radical exchange process: <math> X^{\bullet} + A-Y \rightarrow X-A + ^{\bullet}Y </math> (X = A = Y = hydrogen atom) Write the HL wave functions for R and R* and derive the value of G using semiempirical VB theory. Hints : 1) write the wave functions of R and R* so that their overlap is positive; 2) neglect the overlap beween the external atoms X and Y, and neglect the overlap between two different determinants. b/ Considering the following reaction: <math> X^{\bullet} + H-X \rightarrow X-H +^{\bullet}X </math> Use semiempirical VB theory to derive the following expression for the avoided crossing term <math> B </math> : <math>B = 0.25\Delta E_{ST} </math>' where <math>\Delta E_{ST} </math>' is the singlet-triplet transition energy of the X-H bond at the geometry of the transitions state. Indication: the resonance energy arising from the mixing of two degenerate VB structures is given by the following expression: <math>RE = [H_{12}-E_{ind}S_{12}]/(1+S_{12}) </math> where <math> E_{ind} </math> is the energy of an individual VB structure, and <math> S_{12} </math> and <math> H_{12} </math> are respectively the overlap and Hamiltonian matrix element between R and R*. c/ It is known that for strong binders, at any given bonding distance the singlet-triplet transition energy is larger than twice the bonding energy of the dimer at equilibrium distance, so that one can write the approximate expression <math>\Delta E_{ST} </math>' <math> = 2 BDE </math>, where <math> BDE </math> is the bonding energy of the dimer at equilibrium distance. Using the latter expression, express the avoided crossing term <math> B </math> as a function of the bonding energy of <math> H_{2}</math>. Answer a - Considering the R and R* states for the radical exchange process, the HL wave functions for R and R* are the following : <math>\psi (R)=(|xa\overline{y}|-|x\overline{a}y|)(1/2)^{-1/2} </math> <math>\psi (R^{*})=(|x\overline{a}y|-|\overline{x}ay|)(1/2)^{-1/2}</math> At the equilibrium distance of R, A and Y are close together while X is far away. Therefore R is stabilized by a covalent bond and its energy is <math> 2 \beta S </math>. In R*, the AY entity has two parallel spins 50% of the time (a repulsive situation), and alternated spins 50% of the time (a non-bonding situation). Therefore the energy of R* is half a triplet repulsion, i.e. <math> - \beta S </math>. As <math>\Delta E_{ST} (A-Y)= 4 \beta S </math>, the final G expression is: <math>G=0.75 \Delta E_{ST} (A-Y)</math> b - Let us express the various energies and matrix elements in terms of the usual <math> \beta </math> and <math> S </math> integrals between the X and H orbitals: <math> E_{ind} = 2 \beta S </math> <math> S_{12} = 0.5 </math> <math> H_{12} = 2 \beta S </math> <math> RE = (2 \beta S - 0.5 \beta S )/1.5 = \beta S </math> Since at the transition state geometry <math>\Delta E_{ST} </math>' = <math> 4 \beta S </math>, we have: <math>B = 0.25\Delta E_{ST} </math>' c - Since <math>\Delta E_{ST} </math>' <math> = 2 BDE </math>, we may re-express the avoided crossing term as: <math>B = 0.5BDE </math> 2/ Computer exercise VBSCD for H—H + H. -> H. + H—H at VBSCF then VBCISD level. In this exercise the VBSCD for H—H + H. -> H. + H—H at VBSCF then VBCISD level will be computed with 6-31G**. Computations for reactant and transition state are requested and other points are optional for advanced users. Computations Compute the Energies and Wavefunctions at Reactant and Transition State with Different Sets of VB Structures Write the following VB structure sets, compare and see the difference: all structures; minimal structures for reactant; minimal structures for product. Perform VBSCF and VBCISD calculations for reactant: Perform a VBSCF calculation with "orbtyp=hao" and "boys"; Perform a VBCISD calculation with VBSCF orbital as initial guess. Perform VBSCF and VBCISD calculations with minimal structures for reactant and product. Perform VBSCF and VBCISD calculations for transition state with the same procedure as in step 2. Analysis: Wavefunctions and Energies Compute the Barrier height of the <math>\textrm{S}_{\textrm{N}}2</math> reaction at VBSCF and VBCISD levels. See the difference of the barrier heights. Compare the energies for reactant and product structures at reactant and transition state, by both VBSCF and VBCISD. What's the difference of the energies at different points? Why? Compute the resonance energies at both reactant and transition state points, see the difference of the resonance energies. Optional : Compute all points and draw the VBSCDs at VBSCF and VBCISD levels. >> Answer Exercise 3 (paper exercise) : Conical intersection in H3• radical (for further reading, see S. Shaik and P.C. Hiberty, "The Chemist's Guide to VB theory", Wiley, Hoboken, New Jersey, 2008, pp. 157-161, exercises 6.11-6.14 pp. 174-176, and answers to the exercises pp. 188-192. Consider three hydrogen atoms Ha, Hb, Hc, with respective atomic orbitals a, b and c, and the two VB structures ] and ] . The Ha-Hb and Hb-Hc distances are equal. By using the thumb rules recalled below, where squared overlap terms are neglected, derive the expression of the energies of R and P, and of the reduced Hamiltonian matrix element between R and P for the 3-orbital/3-electrons reacting system [Ha--Hb--Hc]•. From the sign of this latter integral when θ > 60°, derive the expressions of the ground state Ψ≠ and of the first excited state Ψ* of the H3• system. One may drop the normalization constants for simplicity. What bonding scheme does the excited state represent ? Show that the reduced Hamiltonian matrix element is largest in the collinear transition state geometry, and drops to zero in the equilateral triangular structure. Show that R and P VB structures are degenerate in the equilateral triangular structure, and that Ψ≠ and Ψ* are also degenerate in this geometry. We now extend the above conclusions to the allyl radical. What are the bonding schemes corresponding to the ground state and first excited state ? What geometrical distortion would make these two states degenerate ? What would be the end product of a photochemical excitation of allyl radical to its first excited state ? Appendix : Thumb rules for the calculations of effective Hamiltonian matrix elements between determinants. Energy of a determinant D : <math><D|H|D> = -2 \sum_{i<j}^{ } \beta_{ij} S_{ij}</math> (if orbitals i and j have parallel spins) Matrix element between determinants differing by spin inversion of two spin-orbitals : <math> <D|H|D'>=<|...i\overline{j}...||H||...\overline{i}j...|>= -2 \beta_{ij} S_{ij}</math> (for <math>D</math>, <math>D'</math> differing by spin inversion of two spin-orbitals) >> Answer