Jump to content

Rep:Mod:seaunicorn12

From ChemWiki

A note on accuracy

Values in summary tables are quoted to the full number of decimal places as given by Gaussian. Any calculated values are given to fewer decimal places to account for errors and inaccuracy in calculations.

Mod 2: Week 1, Day 1

Part 2: BH3 Molecule optimisation 3-21

The BH3 bond lengths were set to 1.5Å as the starting point.

Optimisation file

File:Vm1110 BH3 OPT.LOG
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: 3-21G
Charge: 0
Spin: Singlet
E(RB3LYP): -26.46226338 au
RMS Gradient Norm: 0.00020672 au
Dipole Moment: 0.0000 Debye
Point Group: D3H
Job cpu time: 39.0 seconds.

The symmetry assigned to the molecule by Gaussian corresponds to the actual symmetry of the molecule, no symmetry restrictions required for further steps.

Item table

 
        Item               Value     Threshold  Converged?
 Maximum Force            0.000413     0.000450     YES
 RMS     Force            0.000271     0.000300     YES
 Maximum Displacement     0.001610     0.001800     YES
 RMS     Displacement     0.001054     0.001200     YES
 Predicted change in Energy=-1.071764D-06
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,2)                  1.1935         -DE/DX =    0.0004              !
 ! R2    R(1,3)                  1.1935         -DE/DX =    0.0004              !
 ! R3    R(1,4)                  1.1935         -DE/DX =    0.0004              !
 ! 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



Part 7: BH3 Molecule optimisation 6-31G(d,p)

Optimisation file

File:Vm1110 BH3 OPT2.LOG
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Charge: 0
Spin: Singlet
E(RB3LYP): -26.61532252 au
RMS Gradient Norm: 0.00021672 au
Dipole Moment: 0.0000 Debye
Point Group: D3H
Job cpu time: 31.0 seconds.

The B-H bond length in the BH3 molecule had decreased from the first optimisation from 1.193 to 1.1914 Å. This is a change of 0.02 ± 0.01 Å.

Item table

 
         Item               Value     Threshold  Converged?
 Maximum Force            0.000433     0.000450     YES
 RMS     Force            0.000284     0.000300     YES
 Maximum Displacement     0.001702     0.001800     YES
 RMS     Displacement     0.001114     0.001200     YES
 Predicted change in Energy=-1.189019D-06
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,2)                  1.1914         -DE/DX =    0.0004              !
 ! R2    R(1,3)                  1.1914         -DE/DX =    0.0004              !
 ! R3    R(1,4)                  1.1914         -DE/DX =    0.0004              !
 ! 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


Mod 2: Week 1, Day 2

Part 1: TlBr3 Molecule optimization LanL2DZ

The point group was restricted to D3H on very tight (0.0001) tolerance.
DOI:10042/21611
File:Vm1110 TlBr opt lanL2DZ.log
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: LANL2DZ
Charge: 0
Spin: Singlet
E(RB3LYP): -91.21812851 au
RMS Gradient Norm: 0.00000090 au
Dipole Moment: 0.0000 Debye
Point Group: D3H
Job cpu time: 19.2 seconds.

         Item               Value     Threshold  Converged?
 Maximum Force            0.000002     0.000450     YES
 RMS     Force            0.000001     0.000300     YES
 Maximum Displacement     0.000022     0.001800     YES
 RMS     Displacement     0.000014     0.001200     YES
 Predicted change in Energy=-6.084147D-11
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,2)                  2.651          -DE/DX =    0.0                 !
 ! R2    R(1,3)                  2.651          -DE/DX =    0.0                 !
 ! R3    R(1,4)                  2.651          -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


Optimised bond length: 2.65 Å (literature: 2.5 Å[1])
Optimised Br-Tl-Br bond angle: 120°
The reported bond length from the optimization is 2.65 ± 0.01 Å which is a little different from the literature value of 6.5 Å[1]. The difference in value could be attributed to the bond length being measured in a different molecule. Also, slight disparities between theoretical calculation and experimentally determined values.

Part 2: BBr3 Molecule optimization [Gen: B(6-31G(d,p); Br(LanL2DZ)]

The structure of the BBr3 molecule was based on the optimized BH3 molecule with H atoms replaced with Br atoms. <bt> DOI:10042/21622
File:Vm1110 Bbr3 opt gen2.log
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: Gen
Charge: 0
Spin: Singlet
E(RB3LYP): -64.43645296 au
RMS Gradient Norm: 0.00000383 au
Dipole Moment: 0.0000 Debye
Point Group: D3H
Job cpu time: 10.9 seconds.

         Item               Value     Threshold  Converged?
 Maximum Force            0.000008     0.000450     YES
 RMS     Force            0.000005     0.000300     YES
 Maximum Displacement     0.000046     0.001800     YES
 RMS     Displacement     0.000030     0.001200     YES
 Predicted change in Energy=-5.266345D-10
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,2)                  1.934          -DE/DX =    0.0                 !
 ! R2    R(1,3)                  1.934          -DE/DX =    0.0                 !
 ! R3    R(1,4)                  1.934          -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


Optimised bond length: 1.934 Å
Optimised Br-B-Br bond angle: 120°

Part 3: Comparing BH3, BBr3 and TlBr3 bond lengths

Bond lengths
BH3 BBr3 TlBr3
1.19 1.93 2.65


Bond length discussion

The central atom-ligand bond increases in the order B-H < B-Br < Tl-Br.
From the values for the B-H and B-Br bond lengths, certain information can be deduced regarding the effect of changing the ligand has on bond length. The bond length increases from 1.19 Å to 1.93 Å which is a 0.74 ± 0.01 Å increase. H is a much smaller atom then Br (smaller covalent and van der Waal radii) which means that when it bonds, it bonds much closer to the central atom, hence a shorter bond length. Br and H are similar in that they both require 1 electron to fulfill their octet and doublet respectively. They are also both quite electronegative, Br more so than H, which means that their orbitals are not very diffuse and the bond lengths are shorter than would be for less electronegative atoms.
Comparing the B-Br bond length with the Tl-Br bond length gives the insight on how changing the central atom has on the bond length between central atom and ligand. Boron and Thallium are both in group 13 so they have the same valence electron structure and require the same number of electrons to complete their octet. Boron is slightly more electronegative than Thallium and as such has less diffuse orbitals. However, as Thallium is in group 6 (vs group 2 for Boron) it has many more orbitals and as such this is the largest contributing factor to it's larger size (rather than electronegativity). The increased size of the atom increases the bond length in this case B-Br and Tl-Br bond lengths are 1.93 Å and 2.65 Å respectively and the increase in bond length is 0.72 ± 0.01 Å.

What is a bond?

Gaussview does not always draw bonds between atoms within a structure even if the bonds exist. Gaussview draws in bonds below a pre-defined value for the particular bond in question. For example, in the optimisation of BH3, a bond length of 1.4 Å is not drawn as a visible bond but a bond of 1.2 Å is. This does not mean that there is no bond, Gaussview assigns bonds according to the distance between atoms rather than bonding or non bonding character of the molecular orbitals. <bt> A bond can be defined in one way as the existence of electrons in molecular orbitals of lower energy than the atomic orbitals, with a largely bonding character. In that way if there are equally as many electrons in antibonding orbitals as there are in bonding orbitals, a bond will not form because the molecule is not stabilized. In addition, any electrons residing in non-bonding orbitals do not contribute to the overall bond.

Mod 2: Day 3

Part 2: BH3 Molecule frequency

File:VM1110 BH3 FREQ.LOG
File Type: .log
Calculation Type: FREQ
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Charge: 0
Spin: Singlet
E(RB3LYP): -26.61532252 au
RMS Gradient Norm: 0.00021669 au
Dipole Moment: 0.0000 Debye
Point Group: D3H
Job cpu time: 8.0 seconds.

 Low frequencies ---  -70.3431  -69.5542  -69.5540   -0.0055    0.0695    0.1595
 Low frequencies --- 1161.4019 1212.1084 1212.1111

The low frequencies do not fall into the ±15 range however there are no negative vibrational modes which indicates that the molecule was in fact optimized to a minimum. The large range for the low frequencies may be in part due to the basis set used or the symmetry of the molecule. As it is a very small molecule it is possible to have lower accuracy and therefore larger errors.

Part 3: BH3 Molecule vibrations

BH3 molecule vibrations
no. Form of the vibration Frequency Intensity symmetry D3H point group
1
The H atoms move in and out of the plane of the molecule in a concerted motion. The B atom moves in and out of the plane of the molecule in the opposite direction of the H atoms and isn't displaced as far.
1161.40 92.7103 A2"
2
Two H atoms move sideways towards one another in a concerted motion. The remaining H atom and the B atom move very slightly together, in and out, away from the two H atoms.
1212.11 13.9813 E'
3
Two of the H atoms bend around the central B atom in a concerted motion. The third H atom also stretches around the central B atom but with a greater overall displacement and in a clockwise direction rather than the anticlockwise direction of the other two atoms. The B atom moves very slightly in the opposite direction of the third H atom.
1212.11 13.9849 E'
4
The H atoms move in and out from the B atom in a concerted motion. The B atom remains motionless. All stretching is done in the plane of the molecule.
2587.81 0 A1'
5
Two of the H atoms are moving in and out from the B atom, in the plane of the molecule, but in an asymmetric manner (while one stretches out, the other stretches in). The third H atom remains fixed. The B atoms moves very slightly in the direction of the H atom coming towards it.
2721.52 126.1376 E'
6
All the H atoms move in and out from the B atom, in the plane of the molecule. Two of the H atoms move in a concerted manner and the third moved in the opposite direction. The B atoms moves very slightly towards the H atoms nearest to it.
2721.52 126.1280 E'


IR spectrum


Only 3 peaks are visible in the IR spectrum. This is because vibrations 5 and 6 and, 2 and 3 have the same frequency and therefore show up as a single peak. In addition, vibration 4 has an intensity of 0, hence it has no peak on the spectrum.

Part 4: TlBr3 Molecule frequency

DOI:10042/22056
File:Vm1110 TlBr opt lanL2DZ freq.log
File Type: .log
Calculation Type: FREQ
Calculation Method: RB3LYP
Basis Set: LANL2DZ
Charge: 0
Spin: Singlet
E(RB3LYP): -91.21812851 au
RMS Gradient Norm: 0.00000088 au
Dipole Moment: 0.0000 Debye
Point Group: D3H
Job cpu time: 22.0 seconds.

  Low frequencies ---   -3.4213   -0.0026   -0.0004    0.0015    3.9367    3.9367
  Low frequencies ---   46.4289   46.4292   52.1449

Lowest frequency mode: 46.4289 (E')

Ir Spectrum

Vibrational frequency comparison

Number Boron Thallium
Frequency Symmetry Frequency Symmetry
1 1161.40 A2" 46.4289 E'
2 1212.11 E' 46.4292 E'
3 1212.11 E' 52.1449 A2"
4 2587.81 A1' 165.2685 A1'
5 2721.52 E' 210.6948 E'
6 2721.52 E' 210.6948 E'


By comparison it is possible to see that the frequencies for TlBr3 are much lower than for BH3. This shows that the Tl-Br bonds are much weaker than B-H bonds. In addition to the vast difference in frequencies, there has also been a reordering of the modes. Where as in BH3 the lowest frequency mode was A2" in TlBr3 the lowest frequency mode is E' symmetry and the A2" mode is of the third lowest frequency.
The two spectra are similar in the manner that the both have two peaks at a lower frequency and one large peak at the highest frequency. The peaks on the thallium spectrum are much broader and far less intense. Also there is more of an overlap of the two peaks in the lower frequency region leading to a more spread double pointed peak rather than two distinct peaks. The A1' and E' vibrations lie at a higher frequency for both spectra because they are of a higher energy than the A2" and E' vibrations.

Frequency calculations

Frequency calculations are run using the same basis set as the optimizations because these determine the symmetry of the molecule and the way in which the energy is measured. The lowest energy version of a molecule in one basis set might not be the lowest energy conformation in another basis set. Hence for the frequency analysis to work the energy must be at a minimum and the basis set must be the same.
The frequency analysis gives information regarding the vibrations of the molecule at its lowest energy state. This information can be used to find which vibrations are IR active and hence create a theoretical IR spectrum.
The 6 low frequencies are the motion of the molecule's centre of mass. When these values are low it means the molecule has been fully optimized and the frequency analysis is not affected by the motion of the overall molecule in measuring the molecule's vibrations.

Part 5: BH3 Molecular Orbitals

DOI:10042/22065
File:Vm1110 BH3 MOs.fchk
File Type: .fch
Calculation Type: SP
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Charge: 0
Spin: Singlet
E(RB3LYP): -26.53393782 au
RMS Gradient Norm: 0.00000000 au
Dipole Moment: 0.0000 Debye
Job cpu time: 0.9 seconds.

Molecular orbital Diagram Calculated molecular orbitals

LUMO (a2")

(e')

(e')

(a1')

Discussion

The LCAO orbitals and the calculated orbitals are all quite similar in appearance, the calculated orbitals more combined to give one overall orbital rather rather than distinct overlapped orbitals. However they are still quite similar. This shows that in this case LCAO theory works well to define what is observed through calculation.

Mod 2: Day 4

Part 1: NBO Analysis

NH3 optimisation

DOI:10042/22066
File:VM1110 NH3 OPT 631 GDP 2.log
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Charge: 0
Spin: Singlet
E(RB3LYP): -56.55776856 au
RMS Gradient Norm: 0.00000888 au
Dipole Moment: 1.8464 Debye
Point Group: C1
Job cpu time: 8.7 seconds.

         Item               Value     Threshold  Converged?
 Maximum Force            0.000024     0.000450     YES
 RMS     Force            0.000012     0.000300     YES
 Maximum Displacement     0.000088     0.001800     YES
 RMS     Displacement     0.000056     0.001200     YES
 Predicted change in Energy=-1.759459D-09
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,2)                  1.018          -DE/DX =    0.0                 !
 ! R2    R(1,3)                  1.018          -DE/DX =    0.0                 !
 ! R3    R(1,4)                  1.018          -DE/DX =    0.0                 !
 ! A1    A(2,1,3)              105.7414         -DE/DX =    0.0                 !
 ! A2    A(2,1,4)              105.7486         -DE/DX =    0.0                 !
 ! A3    A(3,1,4)              105.7478         -DE/DX =    0.0                 !
 ! D1    D(2,1,4,3)           -111.8631         -DE/DX =    0.0                 !
 --------------------------------------------------------------------------------
 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

NH3 frequency calculation

DOI:10042/22067
File:VM1110 NH3 freq.log
File Type: .log
Calculation Type: FREQ
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Charge: 0
Spin: Singlet
E(RB3LYP): -56.55776856 au
RMS Gradient Norm: 0.00000888 au
Dipole Moment: 1.8464 Debye
Point Group: C1
Job cpu time: 17.0 seconds.

 Low frequencies ---  -30.7764   -0.0020   -0.0018   -0.0012   20.3142   28.2484
 Low frequencies --- 1089.5557 1694.1237 1694.1868

NH3 population analysis

File:VM1110 NH3 pop.log
File Type: .log
Calculation Type: SP
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Charge: 0
Spin: Singlet
E(RB3LYP): -56.55776856 au
RMS Gradient Norm: au
Dipole Moment: 1.8464 Debye
Point Group: C1
Job cpu time: 15.8 seconds.

NH3 NBO analysis

NH3 charge distribution (-1.0 to +1.0)
N Charge H Charge
-1.125 0.375


Part 2: Association energy

NH3BH3 optimisation

DOI:10042/
File:VM1110 NH3BH3 opt 631 gdp 2.log
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Charge: 0
Spin: Singlet
E(RB3LYP): -83.22468730 au
RMS Gradient Norm: 0.00005971 au
Dipole Moment: 5.5649 Debye
Point Group: C1
Job cpu time: 16.8 seconds.

         Item               Value     Threshold  Converged?
 Maximum Force            0.000122     0.000450     YES
 RMS     Force            0.000057     0.000300     YES
 Maximum Displacement     0.000662     0.001800     YES
 RMS     Displacement     0.000295     0.001200     YES
 Predicted change in Energy=-1.639776D-07
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,8)                  1.21           -DE/DX =   -0.0001              !
 ! R2    R(2,8)                  1.2101         -DE/DX =   -0.0001              !
 ! R3    R(3,8)                  1.2101         -DE/DX =   -0.0001              !
 ! R4    R(4,7)                  1.0186         -DE/DX =   -0.0001              !
 ! R5    R(5,7)                  1.0186         -DE/DX =   -0.0001              !
 ! R6    R(6,7)                  1.0186         -DE/DX =   -0.0001              !
 ! R7    R(7,8)                  1.6681         -DE/DX =   -0.0001              !
 ! A1    A(4,7,5)              107.8722         -DE/DX =    0.0                 !
 ! A2    A(4,7,6)              107.8744         -DE/DX =    0.0                 !
 ! A3    A(4,7,8)              111.0188         -DE/DX =    0.0                 !
 ! A4    A(5,7,6)              107.8673         -DE/DX =    0.0                 !
 ! A5    A(5,7,8)              111.0335         -DE/DX =    0.0                 !
 ! A6    A(6,7,8)              111.0296         -DE/DX =    0.0                 !
 ! A7    A(1,8,2)              113.8814         -DE/DX =    0.0                 !
 ! A8    A(1,8,3)              113.879          -DE/DX =    0.0                 !
 ! A9    A(1,8,7)              104.5977         -DE/DX =    0.0                 !
 ! A10   A(2,8,3)              113.8671         -DE/DX =    0.0                 !
 ! A11   A(2,8,7)              104.5945         -DE/DX =    0.0                 !
 ! A12   A(3,8,7)              104.5926         -DE/DX =    0.0                 !
 ! D1    D(4,7,8,1)            179.9734         -DE/DX =    0.0                 !
 ! D2    D(4,7,8,2)            -60.0193         -DE/DX =    0.0                 !
 ! D3    D(4,7,8,3)             59.9696         -DE/DX =    0.0                 !
 ! D4    D(5,7,8,1)            -60.027          -DE/DX =    0.0                 !
 ! D5    D(5,7,8,2)             59.9803         -DE/DX =    0.0                 !
 ! D6    D(5,7,8,3)            179.9692         -DE/DX =    0.0                 !
 ! D7    D(6,7,8,1)             59.9736         -DE/DX =    0.0                 !
 ! D8    D(6,7,8,2)            179.9809         -DE/DX =    0.0                 !
 ! D9    D(6,7,8,3)            -60.0302         -DE/DX =    0.0                 !
 --------------------------------------------------------------------------------
 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

NH3BH3 frequency

DOI:10042/
File:VM1110 NH3BH3 freq 2.log
File Type: .log
Calculation Type: FREQ
Calculation Method: RB3LYP
Basis Set: 6-31G(d,p)
Charge: 0
Spin: Singlet
E(RB3LYP): -83.22468730 au
RMS Gradient Norm: 0.00005971 au
Dipole Moment: 5.5649 Debye
Point Group: C1
Job cpu time: 1 minute 8.9 seconds.

 Low frequencies ---  -13.5525   -0.0010   -0.0008   -0.0006   16.5718   28.9370
 Low frequencies ---  265.0672  632.1225  638.0187

Energy comparison

Molecule Energy
NH3 -56.55776856 au
BH3 -26.61532252 au
NH3BH3 -83.22468730 au
ΔE -0.052 au (-135.47kJ/mol)

Dissociation energy: -135.47 kJ/mol

Mod 2: Week 2, AlCl2Br dimer isomers

Part 1: AlCl2Br dimer optimisations

Isomer 1


File:Vm1110 AlCl4Br2 isomer 1 opt gen 3.log
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: Gen
Charge: 0
Spin: Singlet
E(RB3LYP): -2352.40630797 au
RMS Gradient Norm: 0.00000327 au
Dipole Moment: 0.0012 Debye
Point Group: C2V
Job cpu time: 28.4 seconds.

         Item               Value     Threshold  Converged?
 Maximum Force            0.000008     0.000450     YES
 RMS     Force            0.000004     0.000300     YES
 Maximum Displacement     0.000324     0.001800     YES
 RMS     Displacement     0.000097     0.001200     YES
 Predicted change in Energy=-1.843622D-09
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,8)                  2.0934         -DE/DX =    0.0                 !
 ! R2    R(2,8)                  2.0934         -DE/DX =    0.0                 !
 ! R3    R(3,7)                  2.0934         -DE/DX =    0.0                 !
 ! R4    R(4,7)                  2.0934         -DE/DX =    0.0                 !
 ! R5    R(5,7)                  2.4893         -DE/DX =    0.0                 !
 ! R6    R(5,8)                  2.4893         -DE/DX =    0.0                 !
 ! R7    R(6,7)                  2.4893         -DE/DX =    0.0                 !
 ! R8    R(6,8)                  2.4893         -DE/DX =    0.0                 !
 ! A1    A(7,5,8)               88.2675         -DE/DX =    0.0                 !
 ! A2    A(7,6,8)               88.2675         -DE/DX =    0.0                 !
 ! A3    A(3,7,4)              121.7613         -DE/DX =    0.0                 !
 ! A4    A(3,7,5)              109.8017         -DE/DX =    0.0                 !
 ! A5    A(3,7,6)              109.8017         -DE/DX =    0.0                 !
 ! A6    A(4,7,5)              109.813          -DE/DX =    0.0                 !
 ! A7    A(4,7,6)              109.813          -DE/DX =    0.0                 !
 ! A8    A(5,7,6)               91.7325         -DE/DX =    0.0                 !
 ! A9    A(1,8,2)              121.7613         -DE/DX =    0.0                 !
 ! A10   A(1,8,5)              109.8017         -DE/DX =    0.0                 !
 ! A11   A(1,8,6)              109.8017         -DE/DX =    0.0                 !
 ! A12   A(2,8,5)              109.813          -DE/DX =    0.0                 !
 ! A13   A(2,8,6)              109.813          -DE/DX =    0.0                 !
 ! A14   A(5,8,6)               91.7325         -DE/DX =    0.0                 !
 ! D1    D(8,5,7,3)           -111.7781         -DE/DX =    0.0                 !
 ! D2    D(8,5,7,4)            111.8043         -DE/DX =    0.0                 !
 ! D3    D(8,5,7,6)              0.006          -DE/DX =    0.0                 !
 ! D4    D(7,5,8,1)            111.7781         -DE/DX =    0.0                 !
 ! D5    D(7,5,8,2)           -111.8043         -DE/DX =    0.0                 !
 ! D6    D(7,5,8,6)             -0.006          -DE/DX =    0.0                 !
 ! D7    D(8,6,7,3)            111.7781         -DE/DX =    0.0                 !
 ! D8    D(8,6,7,4)           -111.8043         -DE/DX =    0.0                 !
 ! D9    D(8,6,7,5)             -0.006          -DE/DX =    0.0                 !
 ! D10   D(7,6,8,1)           -111.7781         -DE/DX =    0.0                 !
 ! D11   D(7,6,8,2)            111.8043         -DE/DX =    0.0                 !
 ! D12   D(7,6,8,5)              0.006          -DE/DX =    0.0                 !
 --------------------------------------------------------------------------------
 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad


Isomer 2


File:Vm1110 AlCl4Br2 isomer 3 opt gen 3.log
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: Gen
Charge: 0
Spin: Singlet
E(RB3LYP): -2352.41632880 au
RMS Gradient Norm: 0.00002607 au
Dipole Moment: 0.1920 Debye
Point Group: CS
Job cpu time: 54.9 seconds.

         Item               Value     Threshold  Converged?
 Maximum Force            0.000083     0.000450     YES
 RMS     Force            0.000031     0.000300     YES
 Maximum Displacement     0.000571     0.001800     YES
 RMS     Displacement     0.000264     0.001200     YES
 Predicted change in Energy=-6.679677D-08
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,3)                  2.0926         -DE/DX =    0.0                 !
 ! R2    R(2,3)                  2.0926         -DE/DX =    0.0                 !
 ! R3    R(3,6)                  2.2958         -DE/DX =    0.0                 !
 ! R4    R(3,7)                  2.2958         -DE/DX =    0.0                 !
 ! R5    R(4,5)                  2.2765         -DE/DX =    0.0                 !
 ! R6    R(4,6)                  2.3006         -DE/DX =    0.0                 !
 ! R7    R(4,7)                  2.3006         -DE/DX =    0.0                 !
 ! R8    R(4,8)                  2.2765         -DE/DX =    0.0                 !
 ! A1    A(1,3,2)              121.7655         -DE/DX =    0.0                 !
 ! A2    A(1,3,6)              110.0697         -DE/DX =    0.0                 !
 ! A3    A(1,3,7)              110.0697         -DE/DX =    0.0                 !
 ! A4    A(2,3,6)              110.0696         -DE/DX =    0.0                 !
 ! A5    A(2,3,7)              110.0696         -DE/DX =    0.0                 !
 ! A6    A(6,3,7)               90.3043         -DE/DX =   -0.0001              !
 ! A7    A(5,4,6)              110.2008         -DE/DX =    0.0                 !
 ! A8    A(5,4,7)              110.2008         -DE/DX =    0.0                 !
 ! A9    A(5,4,8)              121.5012         -DE/DX =    0.0                 !
 ! A10   A(6,4,7)               90.0628         -DE/DX =   -0.0001              !
 ! A11   A(6,4,8)              110.2008         -DE/DX =    0.0                 !
 ! A12   A(7,4,8)              110.2008         -DE/DX =    0.0                 !
 ! A13   A(3,6,4)               89.8164         -DE/DX =    0.0001              !
 ! A14   A(3,7,4)               89.8164         -DE/DX =    0.0001              !
 ! D1    D(1,3,6,4)           -111.5487         -DE/DX =    0.0                 !
 ! D2    D(2,3,6,4)            111.5488         -DE/DX =    0.0                 !
 ! D3    D(7,3,6,4)              0.0001         -DE/DX =    0.0                 !
 ! D4    D(1,3,7,4)            111.5487         -DE/DX =    0.0                 !
 ! D5    D(2,3,7,4)           -111.5488         -DE/DX =    0.0                 !
 ! D6    D(6,3,7,4)             -0.0001         -DE/DX =    0.0                 !
 ! D7    D(5,4,6,3)            111.6137         -DE/DX =    0.0                 !
 ! D8    D(7,4,6,3)             -0.0001         -DE/DX =    0.0                 !
 ! D9    D(8,4,6,3)           -111.6138         -DE/DX =    0.0                 !
 ! D10   D(5,4,7,3)           -111.6137         -DE/DX =    0.0                 !
 ! D11   D(6,4,7,3)              0.0001         -DE/DX =    0.0                 !
 ! D12   D(8,4,7,3)            111.6138         -DE/DX =    0.0                 !
 --------------------------------------------------------------------------------
 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

Isomer 3


File:Vm1110 AlCl4Br2 isomer 4 opt gen 2.log
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: Gen
Charge: 0
Spin: Singlet
E(RB3LYP): -2352.41626675 au
RMS Gradient Norm: 0.00002016 au
Dipole Moment: 0.1669 Debye
Point Group: CS
Job cpu time: 52.6 seconds.

         Item               Value     Threshold  Converged?
 Maximum Force            0.000049     0.000450     YES
 RMS     Force            0.000020     0.000300     YES
 Maximum Displacement     0.000512     0.001800     YES
 RMS     Displacement     0.000174     0.001200     YES
 Predicted change in Energy=-2.899010D-08
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,2)                  2.0939         -DE/DX =    0.0                 !
 ! R2    R(2,5)                  2.2981         -DE/DX =    0.0                 !
 ! R3    R(2,6)                  2.2981         -DE/DX =    0.0                 !
 ! R4    R(2,7)                  2.2747         -DE/DX =    0.0                 !
 ! R5    R(3,4)                  2.2747         -DE/DX =    0.0                 !
 ! R6    R(3,5)                  2.2981         -DE/DX =    0.0                 !
 ! R7    R(3,6)                  2.2981         -DE/DX =    0.0                 !
 ! R8    R(3,8)                  2.0937         -DE/DX =    0.0                 !
 ! A1    A(1,2,5)              109.8477         -DE/DX =    0.0                 !
 ! A2    A(1,2,6)              109.8477         -DE/DX =    0.0                 !
 ! A3    A(1,2,7)              121.4998         -DE/DX =    0.0                 !
 ! A4    A(5,2,6)               90.1702         -DE/DX =    0.0                 !
 ! A5    A(5,2,7)              110.5144         -DE/DX =    0.0                 !
 ! A6    A(6,2,7)              110.5144         -DE/DX =    0.0                 !
 ! A7    A(4,3,5)              110.5164         -DE/DX =    0.0                 !
 ! A8    A(4,3,6)              110.5164         -DE/DX =    0.0                 !
 ! A9    A(4,3,8)              121.4893         -DE/DX =    0.0                 !
 ! A10   A(5,3,6)               90.1678         -DE/DX =    0.0                 !
 ! A11   A(5,3,8)              109.8534         -DE/DX =    0.0                 !
 ! A12   A(6,3,8)              109.8534         -DE/DX =    0.0                 !
 ! A13   A(2,5,3)               89.829          -DE/DX =    0.0                 !
 ! A14   A(2,6,3)               89.829          -DE/DX =    0.0                 !
 ! D1    D(1,2,5,3)            110.7482         -DE/DX =    0.0                 !
 ! D2    D(6,2,5,3)             -0.4771         -DE/DX =    0.0                 !
 ! D3    D(7,2,5,3)           -112.519          -DE/DX =    0.0                 !
 ! D4    D(1,2,6,3)           -110.7482         -DE/DX =    0.0                 !
 ! D5    D(5,2,6,3)              0.4771         -DE/DX =    0.0                 !
 ! D6    D(7,2,6,3)            112.519          -DE/DX =    0.0                 !
 ! D7    D(4,3,5,2)            112.5206         -DE/DX =    0.0                 !
 ! D8    D(6,3,5,2)              0.4771         -DE/DX =    0.0                 !
 ! D9    D(8,3,5,2)           -110.7542         -DE/DX =    0.0                 !
 ! D10   D(4,3,6,2)           -112.5206         -DE/DX =    0.0                 !
 ! D11   D(5,3,6,2)             -0.4771         -DE/DX =    0.0                 !
 ! D12   D(8,3,6,2)            110.7542         -DE/DX =    0.0                 !
 --------------------------------------------------------------------------------
 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

Isomer 4


File:Vm1110 AlCl4Br2 isomer 5 opt gen 2.log
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: Gen
Charge: 0
Spin: Singlet
E(RB3LYP): -2352.41631606 au
RMS Gradient Norm: 0.00001798 au
Dipole Moment: 0.0059 Debye
Point Group: CS
Job cpu time: 54.8 seconds.

         Item               Value     Threshold  Converged?
 Maximum Force            0.000046     0.000450     YES
 RMS     Force            0.000017     0.000300     YES
 Maximum Displacement     0.000619     0.001800     YES
 RMS     Displacement     0.000158     0.001200     YES
 Predicted change in Energy=-2.701528D-08
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,4)                  2.2984         -DE/DX =    0.0                 !
 ! R2    R(1,5)                  2.2984         -DE/DX =    0.0                 !
 ! R3    R(1,7)                  2.0938         -DE/DX =    0.0                 !
 ! R4    R(1,8)                  2.2746         -DE/DX =    0.0                 !
 ! R5    R(2,3)                  2.2746         -DE/DX =    0.0                 !
 ! R6    R(2,4)                  2.2985         -DE/DX =    0.0                 !
 ! R7    R(2,5)                  2.2985         -DE/DX =    0.0                 !
 ! R8    R(2,6)                  2.0937         -DE/DX =    0.0                 !
 ! A1    A(4,1,5)               90.1599         -DE/DX =    0.0                 !
 ! A2    A(4,1,7)              109.8678         -DE/DX =    0.0                 !
 ! A3    A(4,1,8)              110.5187         -DE/DX =    0.0                 !
 ! A4    A(5,1,7)              109.8678         -DE/DX =    0.0                 !
 ! A5    A(5,1,8)              110.5187         -DE/DX =    0.0                 !
 ! A6    A(7,1,8)              121.4683         -DE/DX =    0.0                 !
 ! A7    A(3,2,4)              110.4933         -DE/DX =    0.0                 !
 ! A8    A(3,2,5)              110.4933         -DE/DX =    0.0                 !
 ! A9    A(3,2,6)              121.5473         -DE/DX =    0.0                 !
 ! A10   A(4,2,5)               90.1547         -DE/DX =    0.0                 !
 ! A11   A(4,2,6)              109.8433         -DE/DX =    0.0                 !
 ! A12   A(5,2,6)              109.8433         -DE/DX =    0.0                 !
 ! A13   A(1,4,2)               89.8427         -DE/DX =    0.0                 !
 ! A14   A(1,5,2)               89.8427         -DE/DX =    0.0                 !
 ! D1    D(5,1,4,2)             -0.0017         -DE/DX =    0.0                 !
 ! D2    D(7,1,4,2)           -111.2475         -DE/DX =    0.0                 !
 ! D3    D(8,1,4,2)            112.0414         -DE/DX =    0.0                 !
 ! D4    D(4,1,5,2)              0.0017         -DE/DX =    0.0                 !
 ! D5    D(7,1,5,2)            111.2475         -DE/DX =    0.0                 !
 ! D6    D(8,1,5,2)           -112.0414         -DE/DX =    0.0                 !
 ! D7    D(3,2,4,1)            112.0113         -DE/DX =    0.0                 !
 ! D8    D(5,2,4,1)              0.0017         -DE/DX =    0.0                 !
 ! D9    D(6,2,4,1)           -111.2123         -DE/DX =    0.0                 !
 ! D10   D(3,2,5,1)           -112.0113         -DE/DX =    0.0                 !
 ! D11   D(4,2,5,1)             -0.0017         -DE/DX =    0.0                 !
 ! D12   D(6,2,5,1)            111.2123         -DE/DX =    0.0                 !
 --------------------------------------------------------------------------------
 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

AlCl2Br dimer energy differences

Energies relative to isomer 2
Isomer Energy (au) ΔE (au) ΔE (kJ/mol)
1 -2352.40630797 0.0100 26.26
2 -2352.41632880 0 0
3 -2352.41626675 0.0001 0.26
4 -2352.41631606 0.00001 0.02

Discussion

The isomer of lowest energy is isomer 2 and the isomer of highest energy is isomer 1. Isomer 1 has 2 bridging Br atoms whereas isomers 2, 3 and, 4 have bridging Cl atoms. Hence the dimer is most stable with bridging Cl atoms rather than Br. Isomer 2 has both the terminal Br atoms on one side. Isomers 3 and 4 have the terminal Br atoms on opposite ends in cis and trans conformations respectively. The trans isomer is slightly lower in energy than the cis isomer and both are relatively similar in energy to isomer 2. Hence the molecule is most stable when the Br atoms are terminal and on the same end of the dimer.

Momomer optimization

File:Vm1110 AlCl2Br monomer opt gen 3.log
File Type: .log
Calculation Type: FOPT
Calculation Method: RB3LYP
Basis Set: Gen
Charge: 0
Spin: Singlet
E(RB3LYP): -1176.19013667 au
RMS Gradient Norm: 0.00013439 au
Dipole Moment: 0.1194 Debye
Point Group: C2V
Job cpu time: 9.5 seconds.

         Item               Value     Threshold  Converged?
 Maximum Force            0.000181     0.000450     YES
 RMS     Force            0.000118     0.000300     YES
 Maximum Displacement     0.001031     0.001800     YES
 RMS     Displacement     0.000577     0.001200     YES
 Predicted change in Energy=-2.606668D-07
 Optimization completed.
    -- Stationary point found.
                           ----------------------------
                           !   Optimized Parameters   !
                           ! (Angstroms and Degrees)  !
 --------------------------                            --------------------------
 ! Name  Definition              Value          Derivative Info.                !
 --------------------------------------------------------------------------------
 ! R1    R(1,2)                  2.27           -DE/DX =   -0.0002              !
 ! R2    R(1,3)                  2.0883         -DE/DX =    0.0002              !
 ! R3    R(1,4)                  2.0883         -DE/DX =    0.0002              !
 ! A1    A(2,1,3)              120.0373         -DE/DX =    0.0                 !
 ! A2    A(2,1,4)              120.0373         -DE/DX =    0.0                 !
 ! A3    A(3,1,4)              119.9254         -DE/DX =    0.0                 !
 ! D1    D(2,1,4,3)            180.0            -DE/DX =    0.0                 !
 --------------------------------------------------------------------------------
 GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

Dimer dissociation energy

Energies relative to isomer 2
Molecule Energy (au) Energy (kJ/mol)
monomer -1176.19013667 -
dimer -2352.41632880 -
ΔE -1176.226 -3.088 106

Discussion

Isomer 2 of the dimers is much more stable than the monomer on its own. This is shown by the very large dissociation energy.

Part 2: AlCl2Br dimer frequencies

Isomer 1

DOI:10042/
File:Vm1110 AlCl4Br2 isomer 1 freq.log
File Type: .log
Calculation Type: FREQ
Calculation Method: RB3LYP
Basis Set: Gen
Charge: 0
Spin: Singlet
E(RB3LYP): -2352.40630797 au
RMS Gradient Norm: 0.00000327 au
Dipole Moment: 0.0012 Debye
Point Group: C2V
Job cpu time: 2minutes 1.7 seconds.

 Low frequencies ---   -5.1814   -4.9633   -3.1838    0.0030    0.0038    0.0039
 Low frequencies ---   14.8246   63.2887   86.0876
IR spectrum

Isomer 2

DOI:10042/
File:Vm1110 AlCl4Br2 isomer 3 freq.log
File Type: .log
Calculation Type: FREQ
Calculation Method: RB3LYP
Basis Set: Gen
Charge: 0
Spin: Singlet
E(RB3LYP): -2352.41632880 au
RMS Gradient Norm: 0.00002607 au
Dipole Moment: 0.1920 Debye
Point Group: CS
Job cpu time: 3 minutes 46.9 seconds.

 Low frequencies ---   -3.5964   -2.5820   -1.4168   -0.0020    0.0011    0.0027
 Low frequencies ---   17.7720   51.0000   72.1620
IR spectrum

Isomer 3

DOI:10042/
File:Vm1110 AlCl4Br2 isomer 4 freq.log
File Type: .log
Calculation Type: FREQ
Calculation Method: RB3LYP
Basis Set: Gen
Charge: 0
Spin: Singlet
E(RB3LYP): -2352.41626675 au
RMS Gradient Norm: 0.00002016 au
Dipole Moment: 0.1669 Debye
Point Group: CS
Job cpu time: 3 minutes 11.8 seconds.

 Low frequencies ---   -4.0622   -2.3951   -0.0038   -0.0035   -0.0022    1.0273
 Low frequencies ---   17.1580   50.9238   78.5379
IR spectrum

Isomer 4

DOI:10042/
File:Vm1110 AlCl4Br2 isomer 5 freq.log
File Type: .log
Calculation Type: FREQ
Calculation Method: RB3LYP
Basis Set: Gen
Charge: 0
Spin: Singlet
E(RB3LYP): -2352.41631606 au
RMS Gradient Norm: 0.00001798 au
Dipole Moment: 0.0059 Debye
Point Group: CS
Job cpu time: 3 minutes 40.1 seconds.

 Low frequencies ---   -0.0013    0.0020    0.0021    2.1254    2.5567    4.1893
 Low frequencies ---   18.1395   49.0727   72.9198
IR spectrum

IR spectra discussion

More intense bands appear in the spectrum for isomer 1 because more vibrations change the overall dipole of the molecule. In isomers 1 and 2 the Br atoms are close to one another which means that with every vibration there is more electron density moving, leading to more, intense vibrations which show up as a greater number of intense peaks on the IR spectra. Isomer 1 has 18 vibrations and 10 of those are IR inactive. However, of the ones that are IR active, 4 are of a high intensity. The other three isomers have intense IR peaks at approx. 420 which are more intense than any of the isomer 1 peaks but rather than 4, they only have 3 peaks of the same level of intensity. This means that isomer 1 has more vibrations that make a significant change in the overall dipole of the molecule and isomers 3 and 4 have the fewest.

IR discussion

For a vibration to be IR active it needs to create a change in the overall dipole of the molecule. All vibrations that do not create a change are then IR inactive and show an intensity of 0.

Al-Br bond stretch discussion

Energies relative to isomer 2
Molcule Vibrations
Frequency Images Description
Isomer 1 196.90 0.0 Al sideways movemenmt in opposite directions, slight Br movement.
241.00 99.7720 Al opposite movement towards and from Br.
246.75 0.0 Al concerted inward movement.
341.32 160.6510 Al concerted sideways movement.
467.22 346.5254 Al opposite movement towards and from Br, Cl outward stretch.
493.90 0.0 Al inward movement, Cl outward movement, Br fixed.
Isomer 1 186.76 1.5858 Al concerted inwards movemenmt, Br movement toward Al.
270.68 13.3495 Al movement outwards in opposite directions, Br moment towards Al.
418.55 309.6868 Al between Br atoms outward movement, other atoms relatively fixed.
495.50 133.7192 Al between Br atoms sideways movement, other atoms relatively fixed.
502.68 104.9317 Al between Cl atoms inward movement, other atoms relatively fixed.

Mod 2: Week 2, AlCl2Br lowest energy dimer

Isomer 2 MO calculation

DOI:10042/
File:Vm1110 AlCl4Br2 isomer 3 mo.fchk
File Type: .fchk
Calculation Type: SP
Calculation Method: RB3LYP
Basis Set: Gen
Charge: 0
Spin: Singlet
E(RB3LYP): -2352.41632880 au
RMS Gradient Norm: 0.0000000 au
Dipole Moment: 0.1920 Debye
Point Group: CS
Job cpu time:

MO visualisations

MO AO interaction Weak or strong? AO interaction Bonding or antibonding? Nodes Delocalisation
weak pi antibonding between the two 2 s orbitals on bridging Cl 1 node in the centre of the molecule Highly delocalised
weak pi antibonding between two s orbitals on terminal Cl 1 node between terminal Cl atoms moderately delocalised
weak pi antibonding between two s orbitals on terminal Br 1 node between terminal Br atoms moderately delocalised
Weak interactions Antibonding between two p orbitals on terminal Cl. Bonding interaction between bridging Cl MOs node between and on terminal Cl atoms, nodes on bridging Cl atoms Terminal Cl MOs highly delocalised, bridging Cl MOs not very delocalised
Weak interactions Antibonding between two p orbitals on terminal Br. Bonding interaction between bridging Cl MOs node between and on terminal Br atoms, nodes on bridging Cl atoms Terminal Br MOs highly delocalised, bridging Cl MOs not very delocalised

The above table shows a selection of 5 filled non-core molecular orbitals ranging from the most bonding (top) to very antibonding (bottom).

References