IRC in NWChem

Dalam bahasan tentang IRC, kita telah membahas bagaimana menggunakan IRC  di Gaussian dan GAMESS, lantas bagaimana dengan IRC di NWChem? Apakah ada keyword khusus untuk melakukan perhitungan IRC di NWChem? Setelah googling sebentar, akhirnya nemu juga bagaimana menggunakan IRC di NWChem yaitu dengan keyword mepgs.

mepgs merupakan modul yang digunakan untuk mencari dua titik kritis pada potential energy surface yang menghubungkan saddle point dari sebuah molekul. Algoritma yang digunakan dalam mepgs adalah constrained trust region quasi-newton optimization dan update perkiraan energi Hessian.

Beberapa opsi dari mepgs adalah sebagai berikut :

MEPGS
OPTTOL <real opttol default 3e-4>
EPREC <real eprec default 1e-7>
STRIDE <real stride default 0.1>
EVIB <real evib default 1e-4>
MAXMEP <integer maxiter default 250>
MAXITER <integer maxiter default 20>
INHESS <integer inhess default 2>
(BACKWARD || FORWARD) <string default both>
(MSWG || NOMSWG) <string default NOMSWG>
(XYZ || NOXYZ)
END

Penjelasan dari masing-masing bagian dari di atas bisa dilihat di NWChem help hehe, tapi kurang lebih maksudnya sama dengan opsi-opsi yang ada di keyword IRC yang ada di program lain.

Contoh dari perhitungan IRC di NWChem adalah sebagai berikut :

scratch_dir /scr

TITLE “mep test”

echo
start meptest2

GEOMETRY   noautosym
C                    -0.63041882    -0.10565750     0.00000000
N                     0.56689462    -0.05377250     0.00000000
H                    -0.23798723     1.04344672     0.00000000
end

basis
* library DZVP_(DFT_Orbital)
end

task dft freq

freq
reuse meptest2.hess
end

mepgs
maxmep 30
maxiter 20
inhess 2
xyz
evib 0.0005
stride 0.1
print low
backward
end

task dft mepgs

Perhitungan di atas akan melakukan backward IRC. Oke, sekian dulu semoga bermanfa’at!😀

# Diambil dari Forum NWChem
# Perhitungan diambil dari QA/test/mep-test

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s