TUGAS APLIKASI KOMPUTER PENENTUAN ENERGI KESTABILAN REAKSI ESTERIFIKASI DENGAN METODE MEKANIKA KLASIK SENYAWA BUTANOL dan ASAM PROPANOAT Disusun Oleh: Dyas Novi Haryanti 083234206 Kimia B 2008 JURUSAN KIMIA FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM UNIVERSITAS NEGERI SURABAYA 2011
94
Embed
083234206_Dyas Novi H._tugas Aplikom Chem Office Revisi1
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
TUGAS APLIKASI KOMPUTER
PENENTUAN ENERGI KESTABILAN REAKSI
ESTERIFIKASI DENGAN METODE MEKANIKA KLASIK
SENYAWA BUTANOL dan ASAM PROPANOAT
Disusun Oleh:
Dyas Novi Haryanti
083234206
Kimia B 2008
JURUSAN KIMIA
FAKULTAS MATEMATIKA DAN ILMU PENGETAHUAN ALAM
UNIVERSITAS NEGERI SURABAYA
2011
A. JUDUL PERCOBAAN
Penentuan Energi Kestabilan Reaksi Esterifikasi Dengan Metode Mekanika Klasik.
B. TUJUAN PERCOBAAN
Mengetahui kestabilan energi dari reaksi antara asam propanoat dan butanol dengan
menggunakan metode mekanika klasik atau MM2.
C. DASAR TEORI
Terjadinya suatu reaksi antar kedua pereaktan dimulai dengan saling bertabrakan.
Biasanya tabrakan antar molekul itu tidak menimbulkan suatu reaksi, molekul-molekul itu
hanyalah terpental kembali. Agar suatu reaksi dapat terjadi, beberapa molekul dan ion yang
bertabrakan dalam suatu wadah itu harus memiliki cukup energi untuk mencapai keadaan
transisi pada waktu bertabrakan. Pada keadaan transisi molekul-molekul mempunyai pilihan
yang sama mudahnya, kembali menjadi perekasi atau terus menjadi produk. Tetapi setelah
melewati puncak, dengan hambatan terkecillah yang dapat menuju ke produk. Selisih antara
energi potensial rata-rata pereaksi dan produk, ialah perubahan entalpi ΔH untuk reaksi
tersebut.
Dibawah ini adalah gambar yang akan menunjukkan diagram energi untuk
berlangsungnya reaksi SN2. Energi potensial yang dibutuhkan untuk mencapai keadaan
transisi membentuk suatu barier energi, dalam gambar barier ini adalah titik energi
maksimum.
Tiap molekul yang bereaksi dan menghasilkan produk harus melewati keadaan transisi,
baik strukturnya maupun eneginya. Karena energi molekul-molekul tersebut tidak sama, maka
diperlukan waktu agar semua molekul itu bereaksi. Persyaratan waktu ini menimbulkan
pengertian dan besaran yang disebut laju reaksi (rate of reaction). Peningkatan konsentasi
pereaksi akan meningkatkan laju terbentuknya produk, karena akan menambah seringnya
tabrakan antar molekul-molekul. Artinya laju reaksi berbanding lurus dengan konsentrasi
kedua pereaksi.
Secara sederhana pengaruh energi pengaktifan terhadap laju relatif reaksi dapat
dirumuskan sebagai; pada kondisi yang sama, reaksi dengan Eakt rendah akan berjalan
dengan lebih cepat. Makin sedikit energi yang diperlukan untuk reaksi, akan makin banyak
molekul yang memiliki cukup energi untuk bereaksi. Bila satu bahan awal dapat mengalami
dua reaksi berbeda yang tak dapat balik (irreversibel) dan yang menghasilkan dua produk
yang berlainan. Bila Eakt dari reaksi balik jauh lebih besar daripada Eakt reaksi, maka reaksi ini
besifat eksotem dan pada hakikatnya tak reversibel. Bila bahan awal dapat mengalami dua
macam reaksi semacam ini, maka produk dari reaksi yang lebih cepat (reaksi dengan Eakt
rendah) akan lebih melimpah.
Eakt ialah energi keadaan tansisi relatif terhadap pereaksi. Oleh karena itu terdapat
hubungan antara laju relatif reaksi dan energi keadaan transisi. Diantara reaksi-reaksi yang
besaing dengan bahan awal sama, reaksi dengan energi keadaan transisi yang rendah adalah
reaksi yang lebih cepat. Suatu spesi dengan energi potensial rendah, akan lebih stabil
dibandingkan dengan spesi berenergi potensial tinggi. Oleh karena itu dapat juga dikatakan,
reaksi dengan struktur keadaan transisi yang lebih stabil adalah reaksi yang lebih cepat.
Gambar 2. Dalam
reaksi – reaksi yang
bersaing, dengan satu
bahan awal, reaksi
dengan Eakt rendah
akan lebih cepat. Jika
reaksi – reaksi itu tak
dapat balik, maka
produk – reaksi yang
lebih cepat akan lebih
banyak.produk
D. PEMBAHASAN
Pada percobaan yang telah dilakukan, tujuan yang diharapkan adalah mengetahui
kestabilan energi dari reaksi antara asam propanoat dan butanol dengan menggunakan metode
mekanika klasik atau MM2. Dengan metode tersebut dapat diketahui perubahan energi
masing-masing senyawa yang dicari baik reaktan, saat keadaan madya (intermediate),
maupun produk.
Reaksi:
HO
butanol +
O
OH
propanoic acid
Reaktan
O
O
O
Madya (intermediet)
O
O
+ H2O
Produk
Langkah yang dilakukan untuk menentukan energi total adalah MM2 dengan cara
mengubah terminate after (molecular dynamics) menjadi 100 steps dan memilih minimize
energy. Langkah tersebut dilakukan sebanyak lima kali untuk membandingkan energi
totalnya. Ternyata ikatan-ikatannya mengalami perubahan setelah kondisi masing-masing
senyawa tersebut distabilkan dengan masing-masing program. Terjadinya perubahan ini
disebabkan ikatan-ikatan yang telah digambarkan belum terstabilkan, selain itu kemampuan
program yang tersedia dapat mengidentifikasi senyawa, sehingga program Chemoffice hanya
mampu menampilkan energi yang dihasilkan. Setelah memperoleh nilai energi total kemudian
menghitung rata-rata masing-masing baik reaktan, madya (intermediet), maupun produk.
Energi rata-rata reaktan =
= 1.60841 kkal/mol
Energi rata-rata madya (intermediet) = 16.7081 kkal/mol
Energi rata-rata produk = 4.34881 kkal/mol
Jadi perubahan entalpi (ΔH) berdasarkan energi total= 4.34881 – 1.60841 = 2.7404 kkal/mol
Energi aktivasi
Ea = energi transisi – energi total reaktan
= 16.7081 kkal/mol – 1.60841 kkal/mol
= 15.09960 kkal/mol
Percobaan ini dilakukan dalam suhu 300 K, Sehingga diperoleh nilai ∆S dan ∆G
sebagai berikut:
∆S = ∆H/∆T = 2.7404/ 300 = 0.00913 kkal/mol K
∆G = ∆H-T∆S = 2.7404- (300 x 0.00913)
= 0.0014 kkal/mol
∆G ≈0 (artinya reaksi dalam keadaan setimbang).
d lnK/ dT = ∆H/ RT2
ln K = ∆H/ RT2
= 2.7404 / ( 8,314 x 300
2) = 3.66 x 10
-6
K = 1.00000366
Berdasarkan reaksi diatas maka dapat diketahui energi dengan program yang dijalankan,
energi stabil dengan menggunakan program MM2 dengan gambar kurva energi adalah
sebagai berikut:
Pada percobaan selanjutnya, yaitu membandingkan nilai yang dihasilkan dari metode
HF dan DFT antara lain energi kinetik, energi potensial, dan energi total. Dari hasil yang
diperoleh, banyak terdapat perbedaan nilai. Setelah memperoleh nilai energinya, dari masing-
masing energi akan dihitung perubahan entalpinya kemudian membuat grafik yang kemudian
juga membandingkan perubahan entalpi dari metode HF dan DFT.
1. Reaktan
a. Butanol
HF DFT
Energi Kinetik 144521.7831 (Kkal/Mol) 144652.3118 (Kkal/Mol)
Energi Potensial -289382.1585 (Kkal/Mol) -290371.8512 (Kkal/Mol)
Energi Total -144860.3754 (Kkal/Mol) -145719.5394 (Kkal/Mol)
b. Asam propanoat
HF DFT
Energi Kinetik 166030.9082 (Kkal/Mol) 166206.6004 (Kkal/Mol)
Energi Potensial -332529.1333 (Kkal/Mol) -333592.9186 (Kkal/Mol)
Energi Total -166498.2251 (Kkal/Mol) -167386.3181 (Kkal/Mol)
Energi reaktan1.60841
Energi transisi intermediet
16.7081
Energi produk4.34881
0
2
4
6
8
10
12
14
16
18
0 0,5 1 1,5 2 2,5 3 3,5
Ene
rgi T
ota
lDiagram Energi
2. Madya atau intermediet
HF DFT
Energi Kinetik 309609.5388 (Kkal/Mol) 309951.3138 (Kkal/Mol)
Energi Potensial -620192.6559 (Kkal/Mol) -622269.6785 (Kkal/Mol)
Energi Total -310583.1171 (Kkal/Mol) -312318.3647 (Kkal/Mol)
3. Produk
a. Butil propanoat
HF DFT
Energi Kinetik 263197.1848 (Kkal/Mol) 263437.6309 (Kkal/Mol)
Energi Potensial -527135.1954 (Kkal/Mol) -528899.0228 (Kkal/Mol)
Energi Total -263938.0106 (Kkal/Mol) -265461.3919 (Kkal/Mol)
b. Air
HF DFT
Energi Kinetik 47396.3303 (Kkal/Mol) 47455.3386 (Kkal/Mol)
Energi Potensial -94825.9378 (Kkal/Mol) -95103.0294 (Kkal/Mol)
Energi Total -47429.6075 (Kkal/Mol) -47647.6908 (Kkal/Mol)
Berdasarkan tabel-tabel diatas, nilai dari metode DFT lebih besar daripada nilai dari
metode HF, baik dari energi kinetik, energi potensial, maupun energi totalnya.
Dari hasil percobaan, juga dapat diperoleh kurva dari masing-masing energi. Berikut
adalah kurva energinya.
1. Energi kinetik
a. Metode HF
Energi kinetik (Kkal/Mol)
Reaktan 155276.34565
Madya 309609.5388
Produk I55296.75755
Sehingga diperoleh nilai perubahan entalpi (∆H) sebesar 20.4119 Kkal/mol.
∆S = ∆H/∆T = 20.4119/300 = 0.0680 kkal/mol K
∆G = ∆H-T∆S = 20.4119- (300 x 0.0680)
= 0.0119 kkal/mol
∆G ≈0 (artinya reaksi dalam keadaan setimbang).
d lnK/ dT = ∆H/ RT2
ln K = ∆H/ RT2
= 20.4119 / ( 8,314 x 300
2) = 2.73 x 10
-6
K = 1.00000273
b. Metode DFT
Energi kinetik (Kkal/Mol)
Reaktan 155429.4561
Madya 309951.3138
Produk 155446.48475
Sehingga diperoleh nilai perubahan entalpi (∆H) sebesar 17.0286 Kkal/mol.
∆S = ∆H/∆T = 17.0286/300 = 0.0568 kkal/mol K
∆G = ∆H-T∆S = 17.0286- (300 x 0.0568)
= 0 kkal/mol
reaktan, 155276.3457
madya, 309609.5388
produk, 155296.76
0
50000
100000
150000
200000
250000
300000
350000
0 1 1 2 2 3 3 4
Kurva Energi Kinetik
∆G ≈0 (artinya reaksi dalam keadaan setimbang).
d lnK/ dT = ∆H/ RT2
ln K = ∆H/ RT2
= 17.0286 / ( 8,314 x 300
2) = 2.27 x 10
-6
K = 1.00000227
2. Energi potensial
a. Metode HF
Energi potensial (Kkal/Mol)
Reaktan 310955.6459
Madya 620192.6559
Produk 310980.5666
Sehingga diperoleh nilai perubahan entalpi (∆H) sebesar 24.9207 Kkal/mol.
∆S = ∆H/∆T = 24.9207/300 = 0.083 kkal/mol K
∆G = ∆H-T∆S = 24.9207- (300 x 0.083)
= 0 kkal/mol
∆G ≈0 (artinya reaksi dalam keadaan setimbang).
reaktan, 155429.4561
madya, 309951.3138
produk, 155446.4848
0
50000
100000
150000
200000
250000
300000
350000
0 1 1 2 2 3 3 4
Kurva Energi Kinetik
d lnK/ dT = ∆H/ RT2
ln K = ∆H/ RT2
= 24.9207 / ( 8,314 x 300
2) = 3.33 x 10
-6
K = 1.00000333
b. Metode DFT
Energi potensial (Kkal/Mol)
Reaktan 311982.8849
Madya 622269.6785
Produk 312002.0261
Sehingga diperoleh nilai perubahan entalpi (∆H) sebesar 19.1412 Kkal/mol.
∆S = ∆H/∆T = 19.1412/300 = 0.0680 kkal/mol K
∆G = ∆H-T∆S = 19.1412- (300 x 0.0680)
= 0 kkal/mol
∆G ≈0 (artinya reaksi dalam keadaan setimbang).
d lnK/ dT = ∆H/ RT2
ln K = ∆H/ RT2
= 19.1412 / ( 8,314 x 300
2) = 2.55 x 10
-6
K = 1.00000255
reaktan, 310955.6459
madya, 620192.6559
produk, 310980.5666
0
100000
200000
300000
400000
500000
600000
700000
0 1 1 2 2 3 3 4
Kurva Energi Potensial
3. Energi total
a. Metode HF
Energi total (Kkal/Mol)
Reaktan 155679.30025
Madya 310583.1171
Produk 155683.80905
Sehingga diperoleh nilai perubahan entalpi (∆H) sebesar 4.5088 Kkal/mol.
∆S = ∆H/∆T = 4.5088/300 = 0.0680 kkal/mol K
∆G = ∆H-T∆S = 4.5088- (300 x 0.0680)
= 0 kkal/mol
∆G ≈0 (artinya reaksi dalam keadaan setimbang).
d lnK/ dT = ∆H/ RT2
ln K = ∆H/ RT2
= 4.5088 / ( 8,314 x 300
2) = 0.60 x 10
-6
K = 1.00000060
reaktan, 311982.8849
madya, 622269.6785
produk, 312002.0261
0
100000
200000
300000
400000
500000
600000
700000
0 1 1 2 2 3 3 4
Kurva Energi Potensial
b. Metode DFT
Energi total (Kkal/Mol)
Reaktan 156552.92875
Madya 312318.3647
Produk 156554.54135
Sehingga diperoleh nilai perubahan entalpi (∆H) sebesar 1.6126 Kkal/mol.
∆S = ∆H/∆T = 1.6126/300 = 0.0680 kkal/mol K
∆G = ∆H-T∆S = 1.6126- (300 x 0.0680)
= 0 kkal/mol
∆G ≈0 (artinya reaksi dalam keadaan setimbang).
d lnK/ dT = ∆H/ RT2
ln K = ∆H/ RT2
= 1.6126 / ( 8,314 x 300
2) = 2.15 x 10
-6
K = 1.00000215
reaktan, 155679.3003
madya, 310583.1171
produk, 155683.8091
0
50000
100000
150000
200000
250000
300000
350000
0 1 1 2 2 3 3 4
Kurva Energi Total
E. Simpulan
Berdasarkan hasil percobaan yang dilakukan, nilai energi total masing-masing akan
berbeda apabila molecular dynamics diubah. Dari percobaan juga diperoleh:
Energi rata-rata reaktan sebesar 1.60841 kkal/mol
Energi rata-rata madya (intermediet) sebesar 16.7081 kkal/mol
Energi rata-rata produk sebesar 4.34881 kkal/mol
Perubahan entalpinya (∆H) berdasakan energi total sebesar 2.7404 kkal/mol
Energi aktivasi sebesar 15.09960 kkal/mol
Perubahan entropi sebesar 0.00913 kkal/mol
Energi Gibbs sebesar 0.0014 kkal/mol, ∆G ≈0 (artinya reaksi dalam keadaan setimbang).
Dengan nilai K sebesar 1.00000366
Energi ikat reaktan sebesar 212.0947 kkal/mol
Energi ikat produk sebesar 224.2237 kkal/mol
Perubahan entalpinya (∆H) berdasarkan energi ikat sebesar -12.129 kkal/mol
reaktan, 156552.9288
madya, 312318.3647
produk, 156554.5414
0
50000
100000
150000
200000
250000
300000
350000
0 1 1 2 2 3 3 4
Kurva Energi Total
Berdasarkan metode HF
Perubahan entalpi menurut energi kinetik sebesar 20.4119 Kkal/mol
Perubahan entalpi menurut energi potensial sebesar 24.9207 Kkal/mol
Perubahan entalpi menurut energi total sebesar 4.5088 Kkal/mol
Berdasarkan metode DFT
Perubahan entalpi menurut energi kinetik sebesar 17.0286 Kkal/mol
Perubahan entalpi menurut energi potensial sebesar 19.1412 Kkal/mol
Perubahan entalpi menurut energi total sebesar 1.6126 Kkal/mol
LAMPIRAN
Gambar senyawa dengan model ball and stick
Butanol
Asam Propanoat
Madya (intermediet)
Produk
dan
Program Mekanika Klasik MM2
1. Energi butanol dengan MM2
a. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 27: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.1782
Bend: 0.7727
Stretch-Bend: 0.0765
Torsion: 0.0095
Non-1,4 VDW: -0.6145
1,4 VDW: 2.5225
Dipole/Dipole: 0.0000
Total: 2.9448
Calculation completed
b. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 26: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.1700
Bend: 0.7612
Stretch-Bend: 0.0731
Torsion: 0.0126
Non-1,4 VDW: -0.6067
1,4 VDW: 2.5391
Dipole/Dipole: 0.0000
Total: 2.9493
Calculation completed
c. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 19: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.1738
Bend: 0.7650
Stretch-Bend: 0.0745
Torsion: 0.0109
Non-1,4 VDW: -0.6080
1,4 VDW: 2.5307
Dipole/Dipole: 0.0000
Total: 2.9470
Calculation completed
d. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 39: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.1687
Bend: 0.7558
Stretch-Bend: 0.0725
Torsion: 0.0099
Non-1,4 VDW: -0.6080
1,4 VDW: 2.5459
Dipole/Dipole: 0.0000
Total: 2.9449
Calculation completed
e. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 26: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.1743
Bend: 0.7670
Stretch-Bend: 0.0748
Torsion: 0.0108
Non-1,4 VDW: -0.6110
1,4 VDW: 2.5300
Dipole/Dipole: 0.0000
Total: 2.9459
Calculation completed
2. Energi asam propanoat dengan MM2
a. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 19: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.0877
Bend: 1.6738
Stretch-Bend: 0.0756
Torsion: -3.2390
Non-1,4 VDW: -0.0907
1,4 VDW: 0.4632
Dipole/Dipole: 1.5859
Total: 0.5565
Calculation completed
b. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 19: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.0835
Bend: 1.6017
Stretch-Bend: 0.0725
Torsion: -3.2141
Non-1,4 VDW: -0.0913
1,4 VDW: 0.4834
Dipole/Dipole: 1.5830
Total: 0.5186
Calculation completed
c. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 129: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.0772
Bend: 1.1378
Stretch-Bend: 0.0553
Torsion: -2.8981
Non-1,4 VDW: -0.3796
1,4 VDW: 0.5128
Dipole/Dipole: 1.5811
Total: 0.0865
Calculation completed
d. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 34: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.0748
Bend: 1.1270
Stretch-Bend: 0.0539
Torsion: -2.8782
Non-1,4 VDW: -0.4056
1,4 VDW: 0.5284
Dipole/Dipole: 1.5792
Total: 0.0796
Calculation completed
e. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 29: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.0731
Bend: 1.1077
Stretch-Bend: 0.0515
Torsion: -2.8390
Non-1,4 VDW: -0.4585
1,4 VDW: 0.5971
Dipole/Dipole: 1.5790
Total: 0.1110
Calculation completed
3. Energi madya (intermediet) dengan MM2
a. Calculation started
Adding lone pairs to O [O(20)]
Adding lone pairs to O [O(5)]
Warning: Some parameters are guessed (Quality = 1).
Iteration 62: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.7850
Bend: 2.4080
Stretch-Bend: 0.3622
Torsion: 0.0507
Non-1,4 VDW: -1.5450
1,4 VDW: 7.4285
Dipole/Dipole: 7.2188
Total: 16.7082
Calculation completed
b. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 1: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.7877
Bend: 2.4049
Stretch-Bend: 0.3627
Torsion: 0.0507
Non-1,4 VDW: -1.5450
1,4 VDW: 7.4293
Dipole/Dipole: 7.2179
Total: 16.7081
Calculation completed
c. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 1: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.7868
Bend: 2.4058
Stretch-Bend: 0.3624
Torsion: 0.0507
Non-1,4 VDW: -1.5449
1,4 VDW: 7.4293
Dipole/Dipole: 7.2180
Total: 16.7081
Calculation completed
d. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 1: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.7876
Bend: 2.4041
Stretch-Bend: 0.3626
Torsion: 0.0507
Non-1,4 VDW: -1.5447
1,4 VDW: 7.4299
Dipole/Dipole: 7.2179
Total: 16.7081
Calculation completed
e. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 1: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.7871
Bend: 2.4044
Stretch-Bend: 0.3624
Torsion: 0.0507
Non-1,4 VDW: -1.5446
1,4 VDW: 7.4300
Dipole/Dipole: 7.2181
Total: 16.7080
Calculation completed
4. Energi butil propanoat dengan MM2
a. Calculation started
Note: Some parameters are not finalized (Quality = 3).
Iteration 68: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.5170
Bend: 1.8333
Stretch-Bend: 0.2448
Torsion: -1.2242
Non-1,4 VDW: -1.3388
1,4 VDW: 6.6806
Dipole/Dipole: 1.9562
Total: 8.6689
Calculation completed
b. Calculation started
Note: Some parameters are not finalized (Quality = 3).
Iteration 109: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.5150
Bend: 1.8341
Stretch-Bend: 0.2446
Torsion: -1.2242
Non-1,4 VDW: -1.3389
1,4 VDW: 6.6803
Dipole/Dipole: 1.9579
Total: 8.6688
Calculation completed
c. Calculation started
Note: Some parameters are not finalized (Quality = 3).
Iteration 128: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.5152
Bend: 1.8333
Stretch-Bend: 0.2447
Torsion: -1.2241
Non-1,4 VDW: -1.3389
1,4 VDW: 6.6805
Dipole/Dipole: 1.9581
Total: 8.6688
Calculation completed
d. Calculation started
Note: Some parameters are not finalized (Quality = 3).
Iteration 129: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.5144
Bend: 1.8346
Stretch-Bend: 0.2446
Torsion: -1.2241
Non-1,4 VDW: -1.3389
1,4 VDW: 6.6802
Dipole/Dipole: 1.9580
Total: 8.6688
Calculation completed
e. Calculation started
Note: Some parameters are not finalized (Quality = 3).
Iteration 73: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.5152
Bend: 1.8337
Stretch-Bend: 0.2447
Torsion: -1.2240
Non-1,4 VDW: -1.3389
1,4 VDW: 6.6804
Dipole/Dipole: 1.9576
Total: 8.6687
Calculation completed
5. Energi air dengan MM2
a. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 6: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.0000
Bend: 0.0288
Stretch-Bend: 0.0000
Torsion: 0.0000
Non-1,4 VDW: 0.0000
1,4 VDW: 0.0000
Dipole/Dipole: 0.0000
Total: 0.0288
Calculation completed
b. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 7: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.0000
Bend: 0.0288
Stretch-Bend: 0.0000
Torsion: 0.0000
Non-1,4 VDW: 0.0000
1,4 VDW: 0.0000
Dipole/Dipole: 0.0000
Total: 0.0288
Calculation completed
c. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 6: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.0000
Bend: 0.0288
Stretch-Bend: 0.0000
Torsion: 0.0000
Non-1,4 VDW: 0.0000
1,4 VDW: 0.0000
Dipole/Dipole: 0.0000
Total: 0.0289
Calculation completed
d. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 6: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.0000
Bend: 0.0288
Stretch-Bend: 0.0000
Torsion: 0.0000
Non-1,4 VDW: 0.0000
1,4 VDW: 0.0000
Dipole/Dipole: 0.0000
Total: 0.0288
Calculation completed
e. Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 7: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 0.0000
Bend: 0.0288
Stretch-Bend: 0.0000
Torsion: 0.0000
Non-1,4 VDW: 0.0000
1,4 VDW: 0.0000
Dipole/Dipole: 0.0000
Total: 0.0288
Calculation completed
ENERGI IKAT
1. Butanol
Calculation started
Adding lone pairs to O [O(5)]
Warning: Some parameters are guessed (Quality = 1).
Iteration 44: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 125.8088
Bend: 1.6154
Stretch-Bend: -0.2088
Torsion: 0.0190
Non-1,4 VDW: 0.2044
1,4 VDW: 168.1119
Total: 295.5507
Note: Due to high VDW interactions, some terms were not computed.
cubic stretch: -2.0000
quartic stretch: 2.3330
p->dielec: 1.5000
p->dieled: 1.5000
Bond Length R(0) K(S) Energy
C-C 1.540 1.5230 4.4000 0.0931 [C(1)-C(2)]
C-H 1.114 1.1130 4.6000 0.0001 [C(1)-H(7)]
C-H 1.543 1.1130 4.6000 34.9016 [C(1)-H(8)]
C-H 1.114 1.1130 4.6000 0.0004 [C(1)-H(6)]
C-C 1.549 1.5230 4.4000 0.2103 [C(2)-C(3)]
C-H 1.116 1.1130 4.6000 0.0033 [C(2)-H(10)]
C-H 1.116 1.1130 4.6000 0.0038 [C(2)-H(9)]
C-C 1.540 1.5140 4.4000 0.2068 [C(3)-C(4)]
C-H 1.116 1.1130 4.6000 0.0036 [C(3)-H(12)]
C-H 1.116 1.1130 4.6000 0.0033 [C(3)-H(11)]
C-O 1.423 1.4080 5.3600 0.0801 [C(4)-O(5)]
C-H 1.112 1.1110 4.6000 0.0004 [C(4)-H(14)]
C-H 1.112 1.1110 4.6000 0.0005 [C(4)-H(13)]
O-H 0.966 0.9610 4.6000 0.0069 [O(5)-H(15)]
O-Lp 1.086 0.6000 4.6000 45.1757 [O(5)-Lp(16)]
O-Lp 1.085 0.6000 4.6000 45.1187 [O(5)-Lp(17)]
Energi ikat butanol = 125.8086 kkal/mol
2. Asam Propanoat
Calculation started
Warning: Some parameters are guessed (Quality = 1).
Iteration 95: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 131.3681
Bend: 1.3171
Stretch-Bend: -0.4075
Torsion: -2.8502
Non-1,4 VDW: -0.4016
1,4 VDW: 184.4035
Total: 313.4293
Note: Due to high VDW interactions, some terms were not computed.
cubic stretch: -2.0000
quartic stretch: 2.3330
p->dielec: 1.5000
p->dieled: 1.5000
Bond Length R(0) K(S) Energy
C-C 1.534 1.5230 4.4000 0.0359 [C(1)-C(2)]
C-H 1.542 1.1130 4.6000 34.8560 [C(1)-H(7)]
C-H 1.543 1.1130 4.6000 34.8866 [C(1)-H(8)]
C-H 1.113 1.1130 4.6000 0.0000 [C(1)-H(6)]
C-C 1.800 1.5090 4.4000 16.4890 [C(2)-C(3)]
C-H 1.115 1.1130 4.6000 0.0009 [C(2)-H(10)]
C-H 1.114 1.1130 4.6000 0.0004 [C(2)-H(9)]
C-O 1.207 1.2080 10.8000 0.0009 [C(3)-O(4)]
C-O 1.344 1.3380 5.0500 0.0139 [C(3)-O(5)]
O-H 0.974 0.9720 7.2000 0.0025 [O(5)-H(11)]
O-Lp 0.600 0.6000 4.6000 0.0000 [O(5)-Lp(12)]
O-Lp 1.085 0.6000 4.6000 45.0820 [O(5)-Lp(13)]
Energi Ikat Asam Propanoat = 86.2861 kkal/mol
3. Butil Propanoat
Calculation started
Note: Some parameters are not finalized (Quality = 3).
Iteration 70: Minimization terminated normally because the gradient norm is less than the
minimum gradient norm
Stretch: 157.1121
Bend: 3.4881
Stretch-Bend: -0.9652
Torsion: -1.2100
Non-1,4 VDW: -1.4493
1,4 VDW: 237.4288
Total: 394.4046
Note: Due to high VDW interactions, some terms were not computed.
cubic stretch: -2.0000
quartic stretch: 2.3330
p->dielec: 1.5000
p->dieled: 1.5000
Bond Length R(0) K(S) Energy
C-C 1.534 1.5230 4.4000 0.0352 [C(1)-C(2)]
C-H 1.113 1.1130 4.6000 0.0001 [C(1)-H(10)]
C-H 1.542 1.1130 4.6000 34.8732 [C(1)-H(11)]
C-H 1.543 1.1130 4.6000 34.8870 [C(1)-H(12)]
C-C 1.801 1.5090 4.4000 16.5979 [C(2)-C(3)]
C-H 1.114 1.1130 4.6000 0.0005 [C(2)-H(13)]
C-H 1.115 1.1130 4.6000 0.0008 [C(2)-H(14)]
C-O 1.209 1.2080 10.8000 0.0005 [C(3)-O(4)]
C-O 1.363 1.3380 5.0500 0.2182 [C(3)-O(5)]
O-C 1.417 1.3890 5.3600 0.2781 [O(5)-C(6)]
O-Lp 0.601 0.6000 4.6000 0.0006 [O(5)-Lp(24)]
O-Lp 0.601 0.6000 4.6000 0.0006 [O(5)-Lp(25)]
C-C 1.534 1.5140 4.4000 0.1248 [C(6)-C(7)]
C-H 1.114 1.1110 4.6000 0.0039 [C(6)-H(15)]
C-H 1.115 1.1110 4.6000 0.0040 [C(6)-H(16)]
C-C 1.550 1.5230 4.4000 0.2154 [C(7)-C(8)]
C-H 1.116 1.1130 4.6000 0.0036 [C(7)-H(17)]
C-H 1.116 1.1130 4.6000 0.0038 [C(7)-H(18)]
C-C 1.540 1.5230 4.4000 0.0919 [C(8)-C(9)]
C-H 1.116 1.1130 4.6000 0.0031 [C(8)-H(19)]
C-H 1.116 1.1130 4.6000 0.0030 [C(8)-H(20)]
C-H 1.114 1.1130 4.6000 0.0002 [C(9)-H(21)]
C-H 1.543 1.1130 4.6000 34.8816 [C(9)-H(22)]
C-H 1.543 1.1130 4.6000 34.8841 [C(9)-H(23)]
Energi Ikat Butil Propanoat = 157.1109 kkal/mol
4. Energi Ikat Air
Bond Length R(0) K(S) Energy
O-H 1.363 0.9420 4.6000 33.5564
O-H 1.363 0.9420 4.6000 33.5564
O-Lp 1.085 0.6000 4.6000 45.1037
O-Lp 0.600 0.6000 4.6000 0.0000
Energi Ikat Air = 67.1128 kkal/mol
Energi Ikat Reaktan = 212.0947 kkal/mol
Energi Ikat Produk = 224.2237 kkal/mol
Perubahan entalpi (∆H) berdasarkan energi ikat = Energi Ikat Reaktan – Energi Ikat Produk