Jump to content

Rep:Mod:syq1234

From ChemWiki

Inorganic Computational Chemistry Lab

Day 1

Optimising a Molecule of BH3

Method

A BH3 molecule was created by GaussView 5.0. The bond length of the three B-H bonds were set to 1.53 angstrom, 1.54 angstrom, 1.55 angstrom separately. The BH3 molecule was optimized using the method: B3LYP and the basis set: 3-21G.

Results

The information of calculation is summarized in the picture. 1.

Picture. 1

















The "Item" table which includes converged forces and distances is displayed below. And It shows that the calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000220     0.000450     YES
RMS     Force            0.000106     0.000300     YES
Maximum Displacement     0.000919     0.001800     YES
RMS     Displacement     0.000447     0.001200     YES
Predicted change in Energy=-1.672480D-07
Optimization completed.
   -- Stationary point found.
                          ----------------------------
                          !   Optimized Parameters   !
                          ! (Angstroms and Degrees)  !
--------------------------                            --------------------------
! Name  Definition              Value          Derivative Info.                !
--------------------------------------------------------------------------------
! R1    R(1,2)                  1.1944         -DE/DX =   -0.0001              !
! R2    R(1,3)                  1.1947         -DE/DX =   -0.0002              !
! R3    R(1,4)                  1.1948         -DE/DX =   -0.0002              !
! A1    A(2,1,3)              119.9983         -DE/DX =    0.0                 !
! A2    A(2,1,4)              119.986          -DE/DX =    0.0                 !
! A3    A(3,1,4)              120.0157         -DE/DX =    0.0                 !
! D1    D(2,1,4,3)            180.0            -DE/DX =    0.0                 !
--------------------------------------------------------------------------------
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

The original output Log. file could be check through the link Optimisation

Picture. 2
Picture. 3

As present in Picture. 2 and Picture. 3, the B-H bond after optimize is 1.19 Å and H-B-H angel is 120 degree.

The energy of each optimization step is displayed in Picture. 4.

Picture. 4



























Day 2

Using a better basis set

Method

A a higher level basis set: 6-31G(d,p)was chose for optimization.

Results

The information of calculation is summarized in the picture. 5.

Picture. 5
















The "Item" table which includes converged forces and distances is displayed below. And It shows that the calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000012     0.000450     YES
RMS     Force            0.000008     0.000300     YES
Maximum Displacement     0.000061     0.001800     YES
RMS     Displacement     0.000038     0.001200     YES
Predicted change in Energy=-1.069047D-09
Optimization completed.
   -- Stationary point found.
                          ----------------------------
                          !   Optimized Parameters   !
                          ! (Angstroms and Degrees)  !
--------------------------                            --------------------------
! Name  Definition              Value          Derivative Info.                !
--------------------------------------------------------------------------------
! R1    R(1,2)                  1.1923         -DE/DX =    0.0                 !
! R2    R(1,3)                  1.1923         -DE/DX =    0.0                 !
! R3    R(1,4)                  1.1923         -DE/DX =    0.0                 !
! A1    A(2,1,3)              120.0055         -DE/DX =    0.0                 !
! A2    A(2,1,4)              119.9938         -DE/DX =    0.0                 !
! A3    A(3,1,4)              120.0007         -DE/DX =    0.0                 !
! D1    D(2,1,4,3)            180.0            -DE/DX =    0.0                 !
--------------------------------------------------------------------------------
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

The original output Log. file could be check through the link Optimisation

Picture. 6
Picture. 7

As present in Picture. 6 and Picture. 7, the B-H bond after optimize is 1.19 Å and H-B-H angel is 120 degree.

The total energy of 3-21G optimised structure is -26.46226429 au. The total energy of 6-31G(d,p) optimised structure is -26.61532360 au. Because two difference basis sets were chose for the optimization, the energies are not related to the stability of structures.

Using pseudo-potentials and larger basis sets

Method

A GaBr3 molecule was created. The symmetry of the molecule was restricted to D3h and the tollerance was adjusted to "Very tight (0.0001)".The GaBr3 molecule was optimized using the method: B3LYP and the basis set: LanL2DZ. The file was uploaded to HPC to run the analysis.

Results

The information of calculation is summarized in the picture. 8.

Picture. 8

















The "Item" table which includes converged forces and distances is displayed below. And It shows that the calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000000     0.000450     YES
RMS     Force            0.000000     0.000300     YES
Maximum Displacement     0.000003     0.001800     YES
RMS     Displacement     0.000002     0.001200     YES
Predicted change in Energy=-1.282679D-12
Optimization completed.
   -- Stationary point found.
                          ----------------------------
                          !   Optimized Parameters   !
                          ! (Angstroms and Degrees)  !
--------------------------                            --------------------------
! Name  Definition              Value          Derivative Info.                !
--------------------------------------------------------------------------------
! R1    R(1,2)                  2.3502         -DE/DX =    0.0                 !
! R2    R(1,3)                  2.3502         -DE/DX =    0.0                 !
! R3    R(1,4)                  2.3502         -DE/DX =    0.0                 !
! A1    A(2,1,3)              120.0            -DE/DX =    0.0                 !
! A2    A(2,1,4)              120.0            -DE/DX =    0.0                 !
! A3    A(3,1,4)              120.0            -DE/DX =    0.0                 !
! D1    D(2,1,4,3)            180.0            -DE/DX =    0.0                 !
--------------------------------------------------------------------------------
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

The HPC output Log. file could be check through the linkDOI:10042/26098

Picture. 9
Picture. 10

As present in Picture. 9 and Picture. 10, the Ga-Br bond length after optimize is 2.35 Å and Br-Ga-Br angel is 120 degree. In literature, The Ga-Br bond length of Gabr3 is 2.249 Å[1]. The optimized Ga-Br bond length is longer than the Ga-Br bond length in the literature. Because the bond length was calculated by Guassian from gaseous Gabr3. However the literature data was measured from crystal solid Gabr3 by experiment. Gabr3 molecules are packed close to each other in solid phase. So the Ga-Br bond length in solid phase is slightly shorter. Meanwhile, the optimized bond length deviates from literature may caused by the basis set chose for the calculation. The basis set was not accurate enough for the calculation. Overall, the bong length difference is 0.101 Å, which is small. It means the optimization is reliable.

Using a mixture of basis-sets and psuedo-potentials

Method

A BBr3 molecule was created. The BBr3 molecule was optimized using 6-31G(d,p) , the basis set: GEN, and the additional keywords section was "pseudo=read gfinput". Finally, the basis sets for all atoms were specified. The file was uploaded to HPC to run the analysis.

Results

The information of calculation is summarized in the picture. 11.

Picture. 11















The "Item" table which includes converged forces and distances is displayed below. And It shows that the calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000020     0.000450     YES
RMS     Force            0.000011     0.000300     YES
Maximum Displacement     0.000077     0.001800     YES
RMS     Displacement     0.000044     0.001200     YES
Predicted change in Energy=-1.117392D-09
Optimization completed.
   -- Stationary point found.
                          ----------------------------
                          !   Optimized Parameters   !
                          ! (Angstroms and Degrees)  !
--------------------------                            --------------------------
! Name  Definition              Value          Derivative Info.                !
--------------------------------------------------------------------------------
! R1    R(1,2)                  1.9339         -DE/DX =    0.0                 !
! R2    R(1,3)                  1.9339         -DE/DX =    0.0                 !
! R3    R(1,4)                  1.934          -DE/DX =    0.0                 !
! A1    A(2,1,3)              119.9967         -DE/DX =    0.0                 !
! A2    A(2,1,4)              119.9968         -DE/DX =    0.0                 !
! A3    A(3,1,4)              120.0064         -DE/DX =    0.0                 !
! D1    D(2,1,4,3)            180.0            -DE/DX =    0.0                 !
--------------------------------------------------------------------------------
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

The HPC output Log. file could be check through the linkDOI:10042/26092

Picture. 12
Picture. 13

As present in Picture. 12 and Picture. 13, the B-Br bond length after optimize is 1.93 Å and Br-B-Br angel is 120 degree. In literature, the the B-Br bond length is 1.892 Å[1].

Comparison of bond distances

The bond distances of BH3, BBr 3, GaBr3 are listed in the table 1.

Table 1

' Bond Length/ Å Bond Length in literature[1]./ Å
B-H of BH3 1.19 1.19
B-Br of BBr3 1.93 1.892
Ga-Br of GaBr3 2.35 2.249

The difference between BH3 and BBr3 is that the -H ligands are replaced by -Br ligands. By the replacement, the bond length increase. The bond length increases form 1.19 Å (B-H) to 1.93 Å (B-Br). It is because that Br is a larger atom and it is more electronegative.

H has the electronic configuration: 1s1. Br has the electronic configuration: [Ar]3d104s24p5. Both of them have a unpaired electron. As a result they are both x-type ligand which donate one unpaired electron to the centre atom to form a sigma bond. However H has a better overlap with B, because the 1s orbital of H has similar size with the sp2 hybridized B orbital. In contrast, the large 4p orbital of Br has a mismatch with the sp2 hybridized B orbital. Though H and Br are both more electronegative than B, Br is more electronegative than H. As a result, B-Br bond is more polar and weaker than B-H bond. In conclusion, B-Br bond is longer, more polar and weaker than B-H bond.

The difference between BBr3 and GsBr3 is that the central B is replaced by Ga. By the replacement, the bond length increase. The bond length increases form 1.93 Å (B-Br) to 2.35 Å (Ga-Br). Because Ga is larger than B.

Ga has the electronic configuration:[Ar]3d104s24p1. B has the electronic configuration: 1s22s22p1. Both of B and Ga are group 3 elements and they can all form sigma bond with Br. However, the large and diffuse (4s,4p)hybridized orbital of Ga has a poor overlap with the large 4p orbital of Br. In contrast (2s,2p)hybridized orbital of B has better overlap with Br. As a results, the Ga-Br bond is weak and long. Meanwhile, Ga is less electronegative than B. As a result, Ga-Br bond is more polar. In conclusion, Ga-Br bond is weaker, longer and more polar than B-Br bond.

In some case, gaussview does not draw bonds as we expect. It does not mean that there is no bond. It is because that the bond length is not in the pre-defined range. "Bond" is defined as the electronic attraction which forms new substance between atoms with opposite charge. In general, there are four types of bonds. They are ionic bond, covalent bond, van der vaals' force and hydrogen bonding[2].

Day 3

Frequency analysis for BH3

Method

The frequency analysis was carried out for the optimized BH3 molecule to analysis the structure of BH3 molecule. The symmetry of the BH3 molecule was restricted to D 3h.

Summary

The information of calculation is summarized in the picture. 14. The summary shows that the " Imaginary frequency" is 0. It means that the molecule was in the stable state.

Picture. 14

















The "Item" table which includes converged forces and distances is displayed below. And It shows that the calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000011     0.000450     YES
RMS     Force            0.000005     0.000300     YES
Maximum Displacement     0.000042     0.001800     YES
RMS     Displacement     0.000021     0.001200     YES
Predicted change in Energy=-6.702220D-10
Optimization completed.
   -- Stationary point found.
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad


The low frequencies lines are displayed below. It shows that the low frequencies are in the range of plus/minus 15cm-1.

Low frequencies ---   -5.0068   -1.2203   -0.0055    0.9880    9.0091    9.0914
Low frequencies --- 1162.9783 1213.1705 1213.1732

The original Log. output could be check through the link Frequency


Animating the vibrations

Method

The vibration models of analyzed BH3 molecule was checked from "Results" - "Vibrations" in Gaussview.

Summary

Six vibrations were found and they are displayed in the Table. 2. Table. 2.

no. form of vibration frequency/ cm-1 Intensity symmetry D3h point group
1
Atom B is stationary. The three H atoms bend in the same direction out of the plane simultaneously
1163 93 A2"
2
The B atom and one H atom are stationary. In the plane of molecule, the other two H atoms bend to each other symmetrically
1213 14 E'
3
The B atom is stationary. In the plane of molecule, the three H atoms bend asymmetrically. Two of the H atoms bend in the same direction. One H atom bends in the opposite direction.
1213 14 E'
4
Atom B is stationary. The three H atoms stretch along the bond symmetrically.
2582 0 A1' totally symmetric
5
The B atom and one H atom are stationary. The other two H atoms stretch along the bond asymmetrically.
2716 126 E'
6
The B atom is stationary. The three H atoms stretch asymmetrically. Two of the H atoms stretch in the same direction. One H atom stretches in the opposite direction.
2716 126 E'
IR

The calculated IR spectrum is showed in Picture. 15

Picture. 15


























In the IR spectrum, there are only three peaks. The three peaks are 1163 cm-1, 1213 cm-1, 2716 cm-1 separately. However, six vibration modes were calculated. The vibration modes 2 and 3 are degenerate and at the frequency 1213 cm-1. The vibration modes 5 and 6 are degenerate and at the frequency 2716 cm-1. The stretch of mode 4 is totally symmetric so this mode cannot be detected by IR.

Frequency analysis for GaBr3

Method

The frequency analysis was run for the optimized GaBr3 to analysis the structure of GaBr3 molecule. The symmetry of the GaBr3 molecule was restricted to D3h. And the file was analysed by HPC system.

Summary

The information of calculation is summarized in the picture. 16.

Picture. 16
















The "Item" table which includes converged forces and distances is displayed below. And It shows that the calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000000     0.000450     YES
RMS     Force            0.000000     0.000300     YES
Maximum Displacement     0.000002     0.001800     YES
RMS     Displacement     0.000001     0.001200     YES
Predicted change in Energy=-6.142862D-13
Optimization completed.
   -- Stationary point found.

The low frequencies lines are displayed below. It shows that the low frequencies are in the range of plus/minus 15cm-1.

Low frequencies ---   -0.5252   -0.5247   -0.0024   -0.0010    0.0235    1.2010
Low frequencies ---   76.3744   76.3753   99.6982

The original Log. output could be check through the linkDOI:10042/26141 The lowest "real" normal mode is 76 cm-1. The optimised structure of GaBr3 is a minimum. Because the imaginary frequency of the structure is zero which is shown in the summary and the low frequencies are in the range of plus/minus 15cm-1.


The calculated IR spectrum is showed in Picture. 17

Picture. 17
Comparison of IR
no. BH3 frequency/ cm-1 Intensity symmetry D3h point group no. GaBr3 frequency/ cm-1 Intensity symmetry D3h point group
1 1163 93 A2" 1 76 3 E'
2 1213 14 E' 2 76 3 E'
3 1213 14 E' 3 100 9 A2"
4 2582 0 A1' 4 197 0 A1'
5 2716 126 E' 5 316 57 E'
6 2716 126 E' 6 316 57 E'

The frequency is proportional to the force constant of the bond and reversely proportional to the reduced mass. BH3 vibrational modes have much larger frequency than GaBr3 vibrational modes. Both B and H are small atoms and have small mass. as a reslut the reduce mass of BH3 is small. In contrast, Ga and Br have larger mass. As a reslut, the reduce mass of GaBr3 is large. It also means that B-H bond is much stronger than Ga-Br bond. Because the overlap between the large 4p orbital of Ga and large 4p orbital of Br is poor. However, the overlap between 1s orbital of H and 2p orbital of B is much better. In the IR spectra of both BH3 and GaBr3, three peaks are recorded. However both BH3 and GaBr3 have six vibrational modes. Because they all have the point group D3h. Two of the E' modes degenerated and the other two E'modes also degenerated. As a result, the four E' modes give two peaks. The A1' is totally symmetric and cannot be detected by IR. Overall, three peaks are observed. In both IR spectra, A2" and E' modes are close in energy. A1' and E' modes are close in energy. In A2" and E' modes, no bond stretch could be observed. So the bond length keeps constant during the vibrations. However, in A1' and E' modes, bond stretch could be observed. Vibrational frequency is higher for the modes with bond stretch. As a result, A1' and E' modes are higher in frequency. For BH3, A2" mode is lower in frequency than E' mode. For GaBr3, A2" mode is higher in frequency than E'mode. It is because than the bond length of Ga-Br is longer and more polar than B-H.

Why must you use the same method and basis set for both the optimisation and frequency analysis calculations? If two difference method and basis are chose for the optimization, two different optimized structures would be obtained. And the two different optimized structure have different total energy. The frequency analysis is used to check if the structure is fully optimized. Same method and basis must be applied, or the frequency analysis would not be applied to the optimized structure.

What is the purpose of carrying out a frequency analysis? The frequency analysis is used to check if the structure is fully optimized. Meanwhile, vibrational modes and IR spectrum could be simulated.

What do the "Low frequencies" represent? The low frequency present the motions on the center of mass of a molecule. For a molecule with 4 atoms, 6 vibrational modes could be find, according to the 3N-6 rule.

Molecular Orbitals of BH3

Method

The molecular orbitals of BH3 was obtained by calculating of the electronic structure. Before the calculation, the method was set to energy. "pop=full" was added to the "additional keywords" section. And "Full NBO" was chose. The file was sent to HPC system to run the analysis.

Summary

The information of calculation is summarized in the picture. 18.

Picture. 18
















The original Log. output could be check through the linkDOI:10042/26146

Picture. 19





























Are there any significant differences between the real and LCAO MOs? There are significant differences between the real and LCAO MOs. LCAO MOS only show the orbitals on each individual atom of a molecule. However, the real MOs delocolized on the whole molecule and it is more diffuse. As a result, real MOs could be used to predict the electron density and distribution of a molecule but LCAO MOs cannot. What does this say about the accuracy and usefulness of qualitative MO theory? Qualitative MO theory could be use to predict the position of orbitals and the their combination. And the nodal planes of the combined orbitals could be predicted. However, qualitative MO theory cannot predict the delocalized electron distribution and the shape of electron density.

NBO Analysis

Method

A NH3 was created. The molecule was carried out 6-31G(d,p) optimization ( key words=nosymm , frequency(key words= int=ultrafine) and MO analysis (key words= pop=full).

Summary
Optimization

The information of calculation is summarized in the picture. 20.

Picture. 20
















The "Item" table which includes converged forces and distances is displayed below. And It shows that the calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000002     0.000015     YES
RMS     Force            0.000001     0.000010     YES
Maximum Displacement     0.000005     0.000060     YES
RMS     Displacement     0.000003     0.000040     YES
Predicted change in Energy=-9.677561D-12
Optimization completed.
   -- Stationary point found.

The original Log. output could be check through the link Optimization


Frequency

The information of calculation is summarized in the picture. 21.

Picture. 21

















The "Item" table which includes converged forces and distances is displayed below. And It shows that the calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000002     0.000450     YES
RMS     Force            0.000001     0.000300     YES
Maximum Displacement     0.000006     0.001800     YES
RMS     Displacement     0.000003     0.001200     YES
Predicted change in Energy=-1.607654D-11
Optimization completed.
   -- Stationary point found.
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

The low frequencies lines are displayed below. It shows that the low frequencies are in the range of plus/minus 15cm-1.

Low frequencies ---   -1.2053   -0.0041   -0.0009    0.0753    3.8651    3.8657
Low frequencies --- 1089.3665 1693.9310 1693.9310

The original Log. output could be check through the link Frequency

MO

The information of calculation is summarized in the picture. 22.

Picture. 22
















The original Log. output could be check through the link MO

The HPC analyzed output file could be checked through the linkDOI:10042/26172

Picture. 23 shows the charge distribution. And the charge range is (-1.0 to +1.0).

Picture. 24 shows the specific NBO chargeS for nitrogen and hydrogen. The charge of nitrogen atom is -1.125. The charge of each hydrogen atom is +0.375.

Picture. 23
Picture. 24

Day 4

Association energies: Ammonia-Borane

Method

A NH3BH3 molecule was created and it was optimized at the b3lyp/6-31G(d,p) level. The structure of NH3BH3 molecule was checked by frequency analysis.

Summary
Optimization

The information of calculation is summarized in the picture. 25.

Picture. 25
















The "Item" table which includes converged forces and distances is displayed below. And It shows that the calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000002     0.000015     YES
RMS     Force            0.000001     0.000010     YES
Maximum Displacement     0.000024     0.000060     YES
RMS     Displacement     0.000010     0.000040     YES
Predicted change in Energy=-8.746364D-11
Optimization completed.
   -- Stationary point found.

The original Log. output could be check through the link OPT


Frequency analysis

The information of calculation is summarized in the picture. 26.

Picture. 26
















The "Item" table which includes converged forces and distances is displayed below. And It shows that the calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000004     0.000450     YES
RMS     Force            0.000001     0.000300     YES
Maximum Displacement     0.000022     0.001800     YES
RMS     Displacement     0.000008     0.001200     YES
Predicted change in Energy=-1.098974D-10
Optimization completed.

The low frequencies lines are displayed below. It shows that the low frequencies are in the range of plus/minus 15cm-1.

Low frequencies ---   -5.5395   -0.3588   -0.0541   -0.0011    1.0000    1.1033
Low frequencies ---  263.2916  632.9628  638.4627
   -- Stationary point found.

The original Log. output could be check through the link Frequency

Dissociation Energy

The energy of optimized BH3, NH3, NH3BH3 are listed in the Table. 3.

Table. 3.

Molecule Energy/au
BH3 -26.6153236
NH3 -56.55776872
NH3BH3 -83.22468911

ΔE=E(NH3BH3)-[E(NH3)+E(BH3)]= -0.05159679 au. The association energy in kJ/mol is -134.15 kJ/mol. The dissociation energy is +134.15 kJ/mol. So the dissociate process is endothermic and it is not a spontaneous process. It is because that NH3BH3 is a stable molecule. The stability is due to a lone pair on N donates to the empty p orbital on B to form a N-B bond.

Project: Lewis acids and bases

Optimization

Four isomers of Al2Cl4Br2 were created by Gaussview. The four isomers were optimized using the method: B3LYP and the basis set: Gen. Full basis set 6-31G(d,p) was chose for Al and Cl. A PP LANL2DZdp was chose for Br.

The structure of each isomer and information of calculation are listed in the Table. 4. In the pictures, the green atoms, red atoms and pink atoms are chlorine, bromine and aluminum separately. The calculated point group and real point group are not necessarily to be the same. So the real point groups of isomer are listed at the bottom of Table. 4.

Table. 4 Isomer 1 Isomer 2 Isomer 3 Isomer 4
Sturucture
File:1.mol
File:1.mol
File:1.mol
File:1.mol
File Name js-isomer1-opt-out js-isomer2-clean-opt js-isomer3-opt-out js-isomer4-clean-opt
File Type .log .log .log .log
Calculation Type FOPT FOPT FOPT FOPT
Calculation Method RB3LYP RB3LYP RB3LYP RB3LYP
Basis Set Gen Gen Gen Gen
Charge 0 0 0 0
Spin Singlet Singlet Singlet Singlet
E(RB3LYP) -2352.40630792 a.u -2352.41109946 a.u -2352.41629860 a.u -2352.41626681 a.u
RMS Gradient Norm 0.00001239 a.u 0.00000701 a.u 0.00000721 a.u 0.00000585 a.u
Imaginary Frequency
Dipole Moment 0.0013 Debye 0.1393 Debye 0.0004 Debye 0.1665 Debye
Point Group C2v C1 C1 C2v
Job cpu Time 0 days 0 hours 3 minutes 30.6 seconds 0 days 0 hours 7 minutes 34.5 seconds 0 days 0 hours 7 minutes 5.4 seconds 0 days 0 hours 6 minutes 32.5 seconds
Link to D-Space DOI:10042/26278 DOI:10042/26300 DOI:10042/26301 DOI:10042/26302
Real point Group D2h C1 C2h C2v
Isomer 1
File:1.mol

The "Item" table which includes converged forces and distances is displayed below. And It shows that the optimization calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000041     0.000450     YES
RMS     Force            0.000017     0.000300     YES
Maximum Displacement     0.001604     0.001800     YES
RMS     Displacement     0.000746     0.001200     YES
Predicted change in Energy=-3.911645D-08
Optimization completed.
   -- Stationary point found.    

The original Log. output could be check through the linkDOI:10042/26278


Isomer 2
File:1.mol

The "Item" table which includes converged forces and distances is displayed below. And It shows that the optimization calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000014     0.000450     YES
RMS     Force            0.000004     0.000300     YES
Maximum Displacement     0.000121     0.001800     YES
RMS     Displacement     0.000033     0.001200     YES
Predicted change in Energy=-1.771796D-09
Optimization completed.
   -- Stationary point found.

The original Log. output could be check through the linkDOI:10042/26300

Isomer 3
File:1.mol


The "Item" table which includes converged forces and distances is displayed below. And It shows that the optimization calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000012     0.000450     YES
RMS     Force            0.000005     0.000300     YES
Maximum Displacement     0.000601     0.001800     YES
RMS     Displacement     0.000239     0.001200     YES
Predicted change in Energy=-4.823693D-09
Optimization completed.
   -- Stationary point found.


The original Log. output could be check through the linkDOI:10042/26301

Isomer 4
File:1.mol

The "Item" table which includes converged forces and distances is displayed below. And It shows that the optimization calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000009     0.000450     YES
RMS     Force            0.000004     0.000300     YES
Maximum Displacement     0.000370     0.001800     YES
RMS     Displacement     0.000170     0.001200     YES
Predicted change in Energy=-2.673342D-09
Optimization completed.
   -- Stationary point found.

The original Log. output could be check through the linkDOI:10042/26302

Summary

Symmetry Analysis

The symmetry elements and point group of each isomer are listed in the Table. 5.

Table. 5 Isomer 1 Isomer 2 Isomer 3 Isomer 4
Structure
Symmetry elements C2(x),C2(y),C2(z), σh(xy), σh(xz), σh(zy) no symmetry element C2(y), σh(zx) C2(z), σv(zx), σv(zy)
Point group D2h C1 C2h C2v
Relative energy of isomers

In the structures, the green atoms, red atoms and pink atoms are chlorine, bromine and aluminum separately. The energy of the four isomers are listed in the Table. 6 in both a.u and KJ/mol unit. By comparing the energies, the energy of isomer 3 is found to be the most negative. As a result, isomer 3 is the most stable isomer in these four isomers. In the bottom of Table. 6, the relative energies to the most stable isomer 3 are listed in the unit KJ/mol for each isomer.

Table. 6 Isomer 1 Isomer 2 Isomer 3 Isomer 4
Sturucture
File:1.mol
File:1.mol
File:1.mol
File:1.mol
Energy/ a.u -2352.40630792 -2352.41109946 -2352.41629860 -2352.41626681
Energy/ KJ/mol -6176242.76 -6176255.34 -6176268.99 -6176268.91
Relative energy to isomer 3/ KJ/mol 26.23 13.65 0.00 0.08

Isomer 3 is the most stable conformer. In the isomer 3, the two bromine atoms are both in the terminal position to different aluminum and trans to each other. In isomer 4, the two bromine atoms are both in the terminal position to different aluminum and cis to each other. The energy of isomer 4 is higher than the energy of isomer 3 and the energy difference is 0.08 KJ/mol. The energy difference is small. The radius of bromine atom is larger than the radius of chlorine atom. When the two bromine atoms are in the cis position the isomer suffers more steric hindrance than the two bromine atoms are in the trans position. The higher steric hindrance is the reason for the slightly high energy of isomer 4 than isomer 3. The most unstable conformer id isomer 1. The two bromine atoms are both in the bridging position and the relative energy to isomer 3 is 26.23 KJ/mol. The second unstable conformer is isomer 2. One of the bromine atom is in the bridging position and the relative energy to isomer 3 is 13.65 KJ/mol which is half of 26.23 KJ/mol. It means that the conformers with bromine in the bridging position are unfavoured. Because aluminum and chlorine are both in the row 3 in the periodic table. The size of 3s and 3p hybridized orbitals of aluminum are similar to the size of 3p orbitals of chlorine. So the overlap between the aluminum (sp3 orbital) and chlorine (3p orbital) is good. The overlap of the diffuse 4p orbital of bromine with the aluminum is less good, especially when the bromine is in the bridging position of the 3c-2e- system.

Monomer analysis

The structure of the monomer AlBrCl2 is shown below. In the pictures, the green atoms, red atoms and pink atoms are chlorine, bromine and aluminum separately.

File:1.mol

Same method and basis set as previous were used for the optimization and frequency analysis for the monomer. And the information of calculation is listed in the Table. 7

Table. 7 Optimization of the monomer Frequency analysis of the monomer
File Name MONOMER-OPT-HPCOUT MONOMER-FRE-HOC-OUT
File Type .log .log
Calculation Type FOPT FREQ
Calculation Method RB3LYP RB3LYP
Basis Set Gen Gen
Charge 0 0
Spin Singlet Singlet
E(RB3LYP) -1176.19013679 a.u -1176.19013679 a.u
RMS Gradient Norm 0.00004196 a.u 0.00004194 a.u
Imaginary Frequency 0
Dipole Moment 0.1075 Debye 0.1075 Debye
Point Group C2v C2v
Job cpu Time 0 days 0 hours 1 minutes 5.9 seconds 0 days 0 hours 0 minutes 37.4 seconds
Link to D-Space DOI:10042/26321 DOI:10042/26322
Real point Group C2v C2v

The "Item" table of optimization which includes converged forces and distances is displayed below. And It shows that the optimization calculation was finished for the monomer.

        Item               Value     Threshold  Converged?
Maximum Force            0.000136     0.000450     YES
RMS     Force            0.000073     0.000300     YES
Maximum Displacement     0.000681     0.001800     YES
RMS     Displacement     0.000497     0.001200     YES
Predicted change in Energy=-7.984394D-08
Optimization completed.
   -- Stationary point found.

The "Item" table of frequency analysis which includes converged forces and distances is displayed below. And It shows that the frequency calculation was finished for the monomer.

        Item               Value     Threshold  Converged?
Maximum Force            0.000081     0.000450     YES
RMS     Force            0.000042     0.000300     YES
Maximum Displacement     0.001588     0.001800     YES
RMS     Displacement     0.000974     0.001200     YES
Predicted change in Energy=-1.810813D-07
Optimization completed.
   -- Stationary point found.

The low frequencies lines are displayed below. It shows that the low frequencies are in the range of plus/minus 15cm-1. The calculated imaginary frequency is 0. It means that the monomer was fully optimized.

Low frequencies ---   -0.0025   -0.0020    0.0031    1.3569    3.6367    4.2604
Low frequencies ---  120.5042  133.9178  185.8950
Dissociation energy

Isomer 3 is the most stable conformer. So isomer 3 was used for the analysis of dissociation energy.

The energy information of isomer 3 and the monomer are listed in the Table. 8

Table. 8 Energy/ a.u Energy/ KJ/mol
Al2Cl4Br2 (Isomer 3) -2352.41629860 -6176268.99
AlBrCl2 (monomer) -1176.19013679 -3088087.20

The association energy ΔE=E(Al2Cl4Br2)-2*E(AlBrCl2)= -94.58 kJ/mol. The dissociation energy is 94.58 kJ/mol. So the dissociate process is endothermic and it is not a spontaneous process. It means that isomer 3 is more stable than the isolated monomer. Because in the monomer, there are only 6 electrons in the valence orbital of Al. So Al want to get two electrons to get the octet structure. When forming the dimer, the electron deficiency are released.

Frequency

The frequency analysis was carried out for the optimized isomers to confirm the isomers were fully optimized. Same method and basis set as the optimize process were chose for the frequency analysis. The structure of each isomer and information of calculation are listed in the Table. 9. In the pictures, the green atoms, red atoms and pink atoms are chlorine, bromine and aluminum separately. The calculated point group and real point group are not necessarily to be the same. So the real point groups of isomer are listed at the bottom of Table. 9.

Table. 9 Isomer 1 Isomer 2 Isomer 3 Isomer 4
Sturucture
File:1.mol
File:1.mol
File:1.mol
File:1.mol
File Name js-isomer1-fre-out js-isomer2-clean-fre js-isomer3-fre-out js-isomer4-clean-fre
File Type .log .log .log .log
Calculation Type FREQ FREQ FREQ FREQ
Calculation Method RB3LYP RB3LYP RB3LYP RB3LYP
Basis Set Gen Gen Gen Gen
Charge 0 0 0 0
Spin Singlet Singlet Singlet Singlet
E(RB3LYP) -2352.40630788 a.u -2352.41109946 a.u -2352.41629860 a.u -2352.41626681 a.u
RMS Gradient Norm 0.00003052 a.u 0.00000700 a.u 0.00000718 a.u 0.00000588 a.u
Imaginary Frequency 0 0 0 0
Dipole Moment 0.0000 Debye 0.1393 Debye 0.0004 Debye 0.1665 Debye
Point Group D2h C1 C1 C2v
Job cpu Time 0 days 0 hours 1 minutes 12 seconds 0 days 0 hours 3 minutes 29.3 seconds 0 days 0 hours 3 minutes 38.2 seconds 0 days 0 hours 1 minutes 40.7 seconds
Link to D-Space DOI:10042/26279 DOI:10042/26304 DOI:10042/26305 DOI:10042/26306
Real point Group D2h C1 C2h C2v


Isomer 1
File:1.mol

The "Item" table which includes converged forces and distances is displayed below. And It shows that the frequency calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000064     0.000450     YES
RMS     Force            0.000031     0.000300     YES
Maximum Displacement     0.000778     0.001800     YES
RMS     Displacement     0.000448     0.001200     YES
Predicted change in Energy=-9.667699D-08
Optimization completed.
   -- Stationary point found.

The low frequencies lines are displayed below. It shows that the low frequencies are in the range of plus/minus 15cm-1. The calculated imaginary frequency is 0. It means that the structure was fully optimized.

Low frequencies ---   -5.7324   -5.5948   -3.3971   -0.0035   -0.0023   -0.0022
Low frequencies ---   14.5079   63.1928   86.0727

The simulated IR spectrum of isomer 1 is presented in Picture. 27.

Picture. 27




























Isomer 2
File:1.mol

The "Item" table which includes converged forces and distances is displayed below. And It shows that the frequency calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000020     0.000450     YES
RMS     Force            0.000007     0.000300     YES
Maximum Displacement     0.000937     0.001800     YES
RMS     Displacement     0.000453     0.001200     YES
Predicted change in Energy=-4.164690D-09
Optimization completed.
   -- Stationary point found.

The low frequencies lines are displayed below. It shows that the low frequencies are in the range of plus/minus 15cm-1. The calculated imaginary frequency is 0. It means that the structure was fully optimized.

Low frequencies ---   -2.5287   -0.0018    0.0023    0.0036    0.6262    3.0981
Low frequencies ---   17.1094   55.9276   80.0590

The simulated IR spectrum of isomer 2 is presented in Picture. 28.

Picture. 28




























Isomer 3
File:1.mol

The "Item" table which includes converged forces and distances is displayed below. And It shows that the frequency calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000019     0.000450     YES
RMS     Force            0.000007     0.000300     YES
Maximum Displacement     0.000488     0.001800     YES
RMS     Displacement     0.000224     0.001200     YES
Predicted change in Energy=-7.038728D-09
Optimization completed.
   -- Stationary point found.

The low frequencies lines are displayed below. It shows that the low frequencies are in the range of plus/minus 15cm-1. The calculated imaginary frequency is 0. It means that the structure was fully optimized.

Low frequencies ---   -5.2222   -0.0026   -0.0013    0.0011    1.3436    1.9047
Low frequencies ---   18.1339   49.1069   73.0122

The simulated IR spectrum of isomer 3 is presented in Picture. 29.

Picture. 29






























isomer 4
File:1.mol

The "Item" table which includes converged forces and distances is displayed below. And It shows that the frequency calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000018     0.000450     YES
RMS     Force            0.000006     0.000300     YES
Maximum Displacement     0.000326     0.001800     YES
RMS     Displacement     0.000160     0.001200     YES
Predicted change in Energy=-3.153276D-09
Optimization completed.
   -- Stationary point found.

The low frequencies lines are displayed below. It shows that the low frequencies are in the range of plus/minus 15cm-1. The calculated imaginary frequency is 0. It means that the structure was fully optimized.

Low frequencies ---   -5.2222   -0.0026   -0.0013    0.0011    1.3436    1.9047
Low frequencies ---   18.1339   49.1069   73.0122

The simulated IR spectrum of isomer 4 is presented in Picture. 30.

Picture. 30


























Summary

There are 8 atoms in each isomer. According to the 3N-6 rule, each isomer has 18 vibration modes. However, not all the modes are active in IR. Only these vibrations which change the dipole moment of the isomer are active in IR spectrum.

Isomer 2 has 18 bands in the IR spectrum. It means that all the 18 vibration modes are active in IR. It is because that isomer 2 has low symmetry (point group is C1). There is no symmetry element for isomer 2. Meanwhile isomer 2 has a total dipole moment 0.1393 Debye. So the vibration modes of isomer 2 cause the change of dipole moment easier than these isomers with total dipole moment 0 Debye.

Only 8 vibration modes of isomer 1 are active in IR. As a result, IR for isomer 1 has less bands than the IR for isomer 2. Because isomer 1 is highly symmetric (point group is D2h) and has the total dipole moment 0.0000 Debye. 10 vibration modes of isomer 1 cannot change the total dipole moment.

Isomer 3 is slightly less symmetric (point group is D2h) than isomer 1 and has a total dipole moment 0.0004 Debye. As a result, isomer 3 has one more active vibration mode than isomer 1 which is 9.

Isomer 4 is more symmetry (point group is C2v) than isomer 2 but less symmetric than isomer 2 and isomer 4. Meanwhile, isomer 4 has a comparably large dipole moment 0.1665 Debye. As a result, there are 15 active vibration modes in the IR.

Overall, the highly unsymmetrical isomer 2 has the most bands in IR. Because all the vibration modes change the dipole moment.

A change of dipole moment is the requirement for a active vibration mode. And only these active vibration mode could be detected by IR.


isomer1-mode11
isomer2-mode11
These two pictures are mode 11 for isomer 1 and isomer 2 separately. In these two modes, The bridging bonds stretch and the bridging atoms move oppositely. The difference of these two modes is that one bridging Br is replaced by Cl in isomer 2. The frequency of mode 11 for the isomer 1 is 197 cm-1. The frequency of mode 11 for the isomer 2 is 211 cm-1. Isomer 1 has lower frequency than isomer 2. Frequency is proportional to the bond strength. Br-Al bond is longer and weaker than Cl-AL bond, as the poor overlap between Br and Al. Br-Al bond is even weaker in the bridging position than in the terminal position.


isomer2-mode12
isomer3-mode12
These two pictures are mode 12 for isomer 2 and isomer 3 separately. All the bonds stretch away from the attached Al atom. One Br is in the bridging position in isomer 2. Two Br are in the trans terminal positions in isomer 3. The frequency of isomer 2 of mode 12 is 257 cm-1. The frequency of isomer 3 of mode 12 is 280 cm-1. Frequency is proportional to the bond strength. It means that the Br-Al bond is stonger and shorter in the terminal position than in the bridging position. As a result, isomer 3 with 2 Br at terminal position is most stable.


isomer3-mode15
isomer4-mode15
These two pictures are mode 15 for isomer 3 and isomer 4 separately. In these two modes, The bridging bonds stretch and the bridging atoms move in the same direction. The 2 Br are in tran terminal position in isomer 3. The 2 Br are in cis terminal position in isomer 3. The frequency of isomer 3 of mode 15 is 421 cm-1. The frequency of isomer 4 of mode 15 is 420 cm-1. Frequency is proportional to the bond strength. It means that the Br-Al bond strength an length of isomer 3 and isomer 4 are similar. As a result, whether the two Br are in the cis or trans terminal position will not effect the nature of Br-Al bond.

MO

The molecular orbitals of isomer 3 was obtained by calculating of the electronic structure. Before the calculation, the method was set to energy. "pop=full" was added to the "additional keywords" section. And "Full NBO" was chose. The file was sent to HPC system to run the analysis. The information of calculation is summarized below:















The original output could be check through the link DOI:10042/26408

Highly bonding orbital 31

This MO is for the energy level 31. In this MO, only strong through bond bonding could be found. And there is no node in this MO. Meanwhile, the electron density of the bonding orbital is delocalized over the bridge. Overall, the MO has strong bonding interactions.

Highly bonding orbital 38

This MO is for the energy level 38. In this MO, only strong through bond bonding and weak through space bonding could be found. And there are six angular nodes in this MO. Meanwhile, the electron density of the bonding orbital is delocalized over the bridge. Overall, the MO has strong bonding interactions.

Bonding orbital 40

This MO is for the energy level 40. In this MO, strong through bond bonding and weak through space bonding could be found.Meanwhile weak through space antibonding also could be found. And there are four angular nodes and a node plane in this MO . Meanwhile, the electron density of the bonding orbital is delocalized over the bridge. Overall, the MO has bonding interactions, but the bonding interaction is less than MO38.


Bonding orbital 57

This MO is for the energy level 57. In this MO, strong through bond antibonding and weak through space antibonding could be found. Meanwhile weak through space bonding also could be found. An angular node plane is in this MO . Meanwhile, the MO is not delocalized at all. Al has bigger P orbital than Cl and Br. Because Al is less electronegatice than Cl and Br. The MO is in high energy. As a result, the less electronegative Al contributes more to the MO. Overall, the MO has antibonding interactions, but the antibonding interaction is less than MO60.

Highly antibonding orbital 60

This MO is for the energy level 60. In this MO, strong through bond antibonding and weak through space antibonding could be found.No bonding interaction could be found. An angular node plane is in this MO . Meanwhile, the MO is not delocalized at all. Al has bigger P orbital than Cl and Br. Because Al is less electronegatice than Cl and Br. The MO is in high energy. As a result, the less electronegative Al contributes more to the MO. Overall, the MO has the strongest antibonding interactions in the five chose MOs.

Further study

The same method and basis set as previous were used to run optimization and frequency analysis for Al2I2Cl4. I is a bigger atom than Br. In theory, Al2I2Cl4 should be less stable than isomer 1 because of the poor overlap. The purple atoms are I. The calculation information is listed below:

Picture. 31
Picture. 32
Picture. 33


The "Item" table of optimization analysis which includes converged forces and distances is displayed below. And It shows that the calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000014     0.000450     YES
RMS     Force            0.000005     0.000300     YES
Maximum Displacement     0.001317     0.001800     YES
RMS     Displacement     0.000641     0.001200     YES
Predicted change in Energy=-6.065606D-09
Optimization completed.
   -- Stationary point found.

The original output could be check through the link DOI:10042/26409

The "Item" table of frequency analysis which includes converged forces and distances is displayed below. And It shows that the calculation was finished.

        Item               Value     Threshold  Converged?
Maximum Force            0.000012     0.000450     YES
RMS     Force            0.000004     0.000300     YES
Maximum Displacement     0.001635     0.001800     YES
RMS     Displacement     0.000710     0.001200     YES
Predicted change in Energy=-6.317924D-09
Optimization completed.
   -- Stationary point found.

The low frequencies lines are displayed below. It shows that the low frequencies are in the range of plus/minus 15cm-1. The calculated imaginary frequency is 0. It means that the monomer was fully optimized.

Low frequencies ---   -6.2631   -5.7571   -4.4599   -0.0010   -0.0006    0.0017
Low frequencies ---    9.2351   65.0734   76.7447

The original output could be check through the link DOI:10042/26410

The total energy of the optimized Al2I2Cl4 is -2348.82462070 a.u ( -6166839.07 KJ/mol). This molecule can created by changing the bridging Br in isomer 1 into I. And the total energy of isomer 1 is -6176242.76 KJ/mol. ΔE=E(Isomer 1) - E(Al2I2Cl4 )= -9403.69 KJ/mol. It means that isomer 1 much stable than Al2I2Cl4. The result satisfy the theory.

The monomer AlICl2 was optimized and analysed below:

Picture. 34
Picture. 35
        Item               Value     Threshold  Converged?
Maximum Force            0.000101     0.000450     YES
RMS     Force            0.000050     0.000300     YES
Maximum Displacement     0.000610     0.001800     YES
RMS     Displacement     0.000365     0.001200     YES
Predicted change in Energy=-3.978942D-08
Optimization completed.
   -- Stationary point found.

The original output could be check through the link DOI:10042/26464

The energy of the monomer is -1174.402499 a.u (-3083393.76 KJ/mol)

The association energy ΔE=E(Al2Cl4I2)-2*E(AlICl2)= -51.55 kJ/mol. The dissociation energy is 51.55 kJ/mol. So the dissociate process is endothermic and it is not a spontaneous process. It means that Al2Cl4I2 is more stable than the isolated monomer. Because in the monomer, there are only 6 electrons in the valence orbital of Al. So Al want to get two electrons to get the octet structure. When forming the dimer, the electron deficiency are released. Tough compare to isomer 3, the chose Al2Cl4I2 is less stable. The monomers are still want to couple to form the dimer.

References

  1. 1.0 1.1 1.2 W. M. Haynes, D. R. Lide and T. J. Bruno, CRC handbook of chemistry and physics : a ready-reference book of chemical and physical data, 2012, 93, 9–23.
  2. Carl R. Nave (2005). HyperPhysics http://hyperphysics.phy-astr.gsu.edu/hbase/chemical/bond.html Retrieved May 18, 2005.