Selasa, 26 Maret 2013

PERHITUNGAN ENERGI KEADAAN DASAR MOLEKUL FORMALDEHIDE

Berikut adalah video tutorialnya




PEMBACAAN OUTPUT

Distance matrix (angstroms):
                    1          2          3          4
     1  C    0.000000
     2  O    1.220000   0.000000
     3  H    1.084066   1.632054   0.000000
     4  H    1.084066   1.632054   2.168133   0.000000
 Framework group  C2V[C2(CO),SGV(H2)]
 Deg. of freedom     3
Berikut adalah bentuk output dari inputan yang saya ketikkan pada bagian molekul spesifikasi. Pengisian Z matriks pada molekul spesifikasi berdasarkan geometri formaldehyde berikut pada bidang Cartesian, sesuai gambar berikut;


                         Standard orientation:                        
 ---------------------------------------------------------------------
 Center     Atomic     Atomic              Coordinates (Angstroms)
 Number     Number      Type              X           Y           Z
 ---------------------------------------------------------------------
    1          6             0        0.000000    0.000000   -0.610000
    2          8             0        0.000000    0.000000    0.610000
    3          1             0        0.000000    1.084066   -0.610000
    4          1             0        0.000000   -1.084066   -0.610000
 ---------------------------------------------------------------------
 Rotational constants (GHZ):    213.3488346     45.4544629     37.4711481
    34 basis functions,    64 primitive gaussians,    34 cartesian basis functions
     8 alpha electrons        8 beta electrons
       nuclear repulsion energy        32.1096891643 Hartrees.
 NAtoms=    4 NActive=    4 NUniq=    3 SFac= 2.05D+00 NAtFMM=   60 Big=F
 Harris functional with IExCor=  205 diagonalized for initial guess.
 ExpMin= 1.61D-01 ExpMax= 5.48D+03 ExpMxC= 8.25D+02 IAcc=1 IRadAn=         1 AccDes= 1.00D-06
 HarFok:  IExCor= 205 AccDes= 1.00D-06 IRadAn=         1 IDoV=1
 ScaDFX=  1.000000  1.000000  1.000000  1.000000
 Initial guess orbital symmetries:
       Occupied  (A1) (A1) (A1) (A1) (B2) (A1) (B1) (B2)
       Virtual   (B1) (A1) (B2) (A1) (B1) (A1) (A1) (B2) (A1) (A1)
                 (B2) (B1) (B2) (A1) (B2) (A2) (B1) (A1) (A2) (A1)
                 (A1) (B1) (B2) (A1) (A1) (A1)
Pada rote section diisikan sintaks pop=reg disimi mengirim permintaan outpu 5 orbital molekuler terbaik saat iterasi. Orbital Molekuler dibawah ini hanya sebagian kecil dari yang saya dapatkan. Eigen value berisikan besarnya energi orbital dalam satuan Hartres.
Molecular Orbital Coefficients
                           1         2         3         4         5
                        (A1)--O   (A1)--O   (A1)--O   (A1)--O   (B2)--O
     EIGENVALUES --   -20.55383 -11.31877  -1.40209  -0.82217  -0.75354
   1 1   C  1S          0.00001   0.99558  -0.11335   0.16036   0.00000
   2        2S          0.00044   0.02695   0.21039  -0.33865   0.00000
   3        2PX         0.00000   0.00000   0.00000   0.00000   0.00000
   4        2PY         0.00000   0.00000   0.00000   0.00000   0.40998
   5        2PZ        -0.00004   0.00110   0.17550   0.07678   0.00000
   6        3S         -0.00027  -0.00745   0.07293  -0.31292   0.00000
   7        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
   8        3PY         0.00000   0.00000   0.00000   0.00000   0.18023
   9        3PZ        -0.00035  -0.00035  -0.00901   0.05604   0.00000
  10        4XX        -0.00001  -0.00289  -0.01726   0.01149   0.00000
  11        4YY        -0.00005  -0.00175  -0.00318  -0.04107   0.00000
  12        4ZZ        -0.00068  -0.00171   0.03165   0.01732   0.00000
  13        4XY         0.00000   0.00000   0.00000   0.00000   0.00000
  14        4XZ         0.00000   0.00000   0.00000   0.00000   0.00000
  15        4YZ         0.00000   0.00000   0.00000   0.00000   0.01368
  16 2   O  1S          0.99468  -0.00042  -0.19191  -0.11755   0.00000
  17        2S          0.02105   0.00014   0.42369   0.26410   0.00000
  18        2PX         0.00000   0.00000   0.00000   0.00000   0.00000
  19        2PY         0.00000   0.00000   0.00000   0.00000   0.30586
  20        2PZ        -0.00144  -0.00033  -0.13908   0.22928   0.00000
  21        3S          0.00436   0.00059   0.37658   0.39634   0.00000
  22        3PX         0.00000   0.00000   0.00000   0.00000   0.00000
  23        3PY         0.00000   0.00000   0.00000   0.00000   0.16604
  24        3PZ         0.00007   0.00067  -0.05374   0.11319   0.00000
  25        4XX        -0.00425   0.00004  -0.00293   0.00078   0.00000
  26        4YY        -0.00384  -0.00030  -0.00034   0.00482   0.00000
  27        4ZZ        -0.00363  -0.00017   0.01870  -0.01962   0.00000
  28        4XY         0.00000   0.00000   0.00000   0.00000   0.00000
  29        4XZ         0.00000   0.00000   0.00000   0.00000   0.00000
  30        4YZ         0.00000   0.00000   0.00000   0.00000  -0.02188
  31 3   H  1S         -0.00004  -0.00043   0.04580  -0.15071   0.19509
  32        2S         -0.00010   0.00183   0.00133  -0.06376   0.07320
  33 4   H  1S         -0.00004  -0.00043   0.04580  -0.15071  -0.19509
  34        2S         -0.00010   0.00183   0.00133  -0.06376  -0.07320
.
Perhitungan Singgle point energy dengan metode RHF memberikan nilai sebagai berikut dalam satuan Hartress (1 Hartress = 2,625 kj/mol), dengan iterasi sebanyak 6.
 SCF Done:  E(RHF) =  -113.753198883     A.U. after    6 cycles
             Convg  =    0.5610D-04             -V/T =  2.0005
             S**2   =   0.000
Secara default, pekerjaan Gaussian melakukan analisis populasi Mulliken, yang partisi  muatan total di antara atom-atom dalam molekul. Berikut adalah bagian penting dari output untuk formaldehida. Analisis ini menempatkan muatan negatif sedikit pada atom oksigen dan  menyeimbangkan dengan memberikan muatan positif antara 2 atom tersisa.
Condensed to atoms (all electrons):
 Mulliken atomic charges:
              1
     1  C   -0.109690
     2  O   -0.491578
     3  H    0.300634
     4  H    0.300634
 Sum of Mulliken charges=   0.00000
 Electronic spatial extent (au):  <R**2>=    55.5437
 Charge=     0.0000 electrons
Gaussian juga memprediksi momen dipol dan momen muitipole lebih tinggi (melalui  hexadecapole). Momen dipol adalah turunan pertama dari energi sehubungan dengan  medan listrik diterapkan. Ini adalah ukuran dari asimetri dalam distribusi muatan molekul, dan diberikan sebagai vektor dalam tiga dimensi. Untuk Hartree-Fock  perhitungan, ini setara dengan nilai ekspektasi dari X, Y, dan Z, yang merupakan  jumlah yang dilaporkan dalam output.  Berikut adalah momen dipol dan quadrupole diprediksi untuk formaldehida:
Dipole moment (field-independent basis, Debye):
    X=     0.0000    Y=     0.0000    Z=    -2.3107  Tot=     2.3107


 Test job not archived.
 1|1|UNPC-UNK|SP|RHF|6-31G(d)|C1H2O1|PCUSER|23-Mar-2013|0||#T RHF/6-31G
 (D) POP=FULL TEST||Singgle Point Energy Formaldehide||0,1|C,0,0.,0.,0.
 |O,0,0.,0.,1.22|H,0,0.94,-0.54,0.|H,0,-0.94,0.54,0.||Version=IA32W-G03
 RevC.01|State=1-A1|HF=-113.7531989|RMSD=5.610e-005|Dipole=0.,0.,-0.909
 0972|PG=C02V [C2(C1O1),SGV(H2)]||@
Gaussian selalu menampilkan banyak kalimat motivasi setiap selesai melakukan pekerjaan, kalimat2 ini berasal dari tokoh-tokoh fisaka yang ternama, diantarannya adalah sebagai berikut,
THE SOLUTION TO A PROBLEM CHANGES THE PROBLEM.
                                 -- JOHN PEERS
           PAUL DICKSON'S "THE OFFICIAL RULES"
Gaussian melaporkan waktu CPU yang digunakan dan ukuran file awal mereka pada saat selesai. Berikut adalah data untuk pekerjaan formaldehida saya:
 Job cpu time:  0 days  0 hours  0 minutes 20.0 seconds.
 File lengths (MBytes):  RWF=     11 Int=      0 D2E=      0 Chk=      7 Scr=      1
 Normal termination of Gaussian 03 at Sat Mar 23 21:59:49 2013.



SEMOGA MANFAAT ^_^

Tidak ada komentar:

Posting Komentar