p~
~
N~
Ht
N...t.- k.. ~
XKL>. ISSN 1410'76g6
SIMULASI FUNGSI DISTRIBUSI RADIAL PARSIAL GELAS AgPO3 DENGAN METODE METROPOLITAN MONTE CARLO Tri Hardi P., Herrj Mugi Ral1ardjo,Setiawan Puslitbang Iptek Bahan -BAT AN; Kawasan Puspiptek Serpong, Tangerang
ABSTRAK SIMULASI FUNGSI DISTRIBUSIRADIAL PARSIALGELAS AgPOsDENGANMETODEMETROPOLITANMONTECARLO. Telah dilakukansimulasifungsi distribusiradial parsial untuk gelas AgPOJdenganmenggunakanmetodeMetropolitanMonte Carlo Bola Keras. Metode ini merupakan metode Reverse Monte Carlo tanpa parametereksperimen. Dari hasil simulasi dapatdiperolehbahwa puncak amort pactapola hamburanneutrondipengaruhioleh ikatan atom P-P, P-O,0-0 dan a-Ago Gelas AgPOJmemiliki struktur lokal dengan kerapatan atom lokal yang lebih rendahdari kerapatanatomrata-rata.
ABSTRACT SIMULATION OFPARSIALRADIALDISTRIBUTION FUNCTION USINGMETROPOLITAN MONTECARLO.Partialradialdistribution functionof AgP03has beenperformed by the usedof HardSpheresMetropolitan MonteCarlo.Thismethodis a ReverseMonteCarlowithout experimental parameters.Fromthesimulation resultsit canbeseen that~omeamorfpeaksatneutronscattering paternsaredueto P-P,P-O,00 andO-Agatomsinteraction. AgP03glasshaslocalstructures withlocaldensitieslowerthanaverage densities.
TEORI
PENDAHULUAN Pergerakan atom di dalam suatu sistem molekular dapat dianalisis dengan dua sistem yaitu sistem deterministik daD sistem non deterministik (stokastik). Sistem deterministik adalah sistem dimana konfigurasi atom menduduki posisi tertentu. Sedangkan sistem stokastik adalah sistem dimana konfigurasi atom berada pacta konfigurasi acak (random). Contoh sistem deterministik adalah simulasi Dinamika Molekuler. daD sistem stokastik adalah simulasi dengan metode Monte Carlo. Tujuan daTi makalah ini adalah untuk menghitung fungsi distribusi radial parsial pada gelas AgPO3. Sirnulasi ini sangat renting untuk menganalisis puncak amorf daTi pola hamburan neutron atau sinar-x. Dengan mengetahui karakteristik fungsi distribusi radial, maka jarak terdekat antar atom dapat ketahui clan bilangan koordinasi dari ikatan tersebut dapat dihitung. Hal ini secara tidak langsung dapat menentukan faktor struktur parsial daTi bahan yang diteliti, sehingga atomatom yang berinteraksi yang mengakibatkan puncakpuncak amorf dapat dianalisis lebih lanjut. Simulasi dilakukan dengan metode Monte Carlo. Dengan perangkat lunak Reverse Monte Carlo (RMC) clan tanpa memasukkan parameter eksperimen metode RMC mirip denganmetode Monte Carlo.
~,
Metode Metropolitan Monte Carlo Metode Monte Carlo merupakan suatu metode numerik non deterministik dimana konfigurasi suatu sistem tergantung pacta sistem sebelumnya. Dalam analisis struktur bahan, sistem adalah sekumpulan posisi atom dalam tiga dimensi yang menghasilkan konfigurasi atom. Probabilitias konfigurasi atom dihitung menggunakan konsep distribusi energi Boltzman. Konfigurasi atom ini bergantung pacta konfigurasi sebelumnya maka prosedur pembobotan cuplik menggunakan rantai Markov yang memenuhi kriteria berikut: .Variabel x(n) dibangkitkan mengikuti kaidah yang memerlukan elemen ke (n+ J) yang mempunyai distribusi probabilitas x(ntI) yang mana hanya tergantung pactadistribusi x(n)dari elemen ke n .Jika P(x=>y) adalah probabilitas yang dicapai y dari keadaan x, maka P( x=>y) harus diijinkan bergerak untuk setiap keadaan di dalam kumpulan atom. .Mikroreversibilitas harus dipenuhi, yaitu bahwa xP( x=>y) = y P( y=>x) saat sistem dalam keadaan setimbang. Dalam satu kumpulan atom dimana jumlah partikel, volume daD suhu adalah tetap, probabilitas (P) dari konfigurasi atom dapat diperoleh dengan algoritma berikut:
21tJ~ 2000
57
~
F~ ~
~
p~
G~ Afo,~
11~ 11~
11-uc..1l.. T...H~ P..lJ,It.
1
N atom diletakkan di dalam sel dengan syarat batas periodik, secara khayal sel dikelilingi oleh dirinya sendiri. Biasanya digunakan sel dalam bentuk kubus. Untuk kubus dengan sisi L, kerapatan atom p = NIL3 harus sarna dengan kerapatan sistem. Probabilitas untuk konfigurasi ini adalah: PI} -exp
(-U,lkT)
(1)
Pn -exp (-U,/kT)
(7)
-;i=( Probabilitas
persamaan :
dimana v,) adalah energi potensial total, yang dihitung atas dasar bentuk potensial yang spesifik dari potensial interatom. T adalah suhu spesifik. Satu atom dipindahkan secara acak, maka probabilitas konfigurasibarudinyatakansebagai:
2.
dimana m adalahjumlah titik Q; dalam AE dan
PI Po = exp (-(Un -U,j/kT)
=exp (-LJU/ kT)
(3)
3.
Jika .dU < 0 konfigurasi barn akan diterima dan menjadi titik awal[ada konfigurasi selanjutnya. Jika .dU>O konfigurasi akan diterima dengan probabilitas sebesar P/P", Apabila.dU = 0 konfigurasi baru ditolak dan perhitungan selanjutnya akan menggunakankonfigurasi sebelurnnya. 4. Prosedur berulang dari langkah 2. Sebagaimana perhitungan konfigurasi atomatom dilakukan pada suhu tertentu, pada suatu slat U akan menurun sampai mencapai satu nilai setimbang dan kemudian akan berosilasi disekitar titik setimbang. Ukuran maksimum dari perpindahan acak secara normal diatur sehingga rasio pergerakan yang diterima terhadap yang ditolak mendekati satu. Metode ReverseMonte Carlo (RMC) Dalam metode RMC diasumsikan bahwa secara eksperimen faktor struktur yang diukur, AE(QJ, hanya mengandung kesalahan statistik yang mempunyai distribusi normal. Perbedaan antara faktor struktur model, AC(QJ, dan yang diukur secara eksperimental adalah:
e, =Ar(QJ -Ah'(QJ
(4)
p -exp (-iI2)
(5) dimana O'(QJadalah deviasi standardistribusi normal. Probabilitastotal dari A" adalah:
dimana perbedaan antara faktor struktur yang diukur A E(Q) dan yang ditentukan dari konfigurasi A c(Q) adalah
Konfigurasi atom yang baru diperoleh dengan tara memindahkan satu atom secara acak. Kemudian hitung fungsi distribusi radial baru gnc(r) dan faktor struktur total AnC(Q), selanjutnya hitung perbedaan faktor struktur basil perhitungan model dengan faktor struktur eksperimen. Perbedaantersebut dituliskan sebagai
.(~I' 58
exp(-
..
., ~ 11.
.
(6)
~I
)
2
OJ
X,,2 = L (A"C'(Qi)-
AJi(Q,)
/U(Q,Y
(10)
Isl
Jika X: < Xo1maka perpindahan atom diterima dan konfigurasi yang baru berubah menjadi konfigurasi lama. Jika Xn1> Xo1maka perpindahan diterima dengan probabilitas exp (-(Xn1 -Xo1)/2). Jika kedua kondisi tersebut tidak diterima maka konfigurasi baru ditolah dan perhitlmgan menggunakankonfigurasi lama. Perhitunganfungsi distribusi radial parsial (gap(r)) Deskripsi tentang distribusi atom dalam bahan non-kristalin sering menggunakan konsep fungsi distribusi radial. Fungsi distribusi radial (g(r) berkaitan dengan probabilitas untuk menemukan atom lain pads suatu jarak r dari atom asal (misalnya pads posisi r=O). Karena fungsi distribusi radial didapatkan dari data difraksi, informasi eksperimental ini berperan dalam pembahasan struktur dan sifat bahan non-kristalin. Konsep fungsi distribusi radial sellin keadaan bahan kristalin juga dapat diperluas ke sistem non-kristalin Hal yang penting dalam simulasi Monte Carlo adalah menentukan konfigurasi awal posisi atom (x,y,z) secara keseluruhan dengan syarat batas periodik, Andaikan jumlah atom adalah N, posisi N atom dipilih secara acak, Selanjutnya dihitung fungsi distribusi radial untllk konfigurasi awal tersebut, Fungsi distribusi radial dinyatakan sebagai:
...(',~ \ gt~'(r)=
~
(8)
(amort).
dan memiliki probabilitas
P(Ii')
dinyatakan dengan
(2)
Karenaitu:
.H
konfigurasi atom
no2VJ 4Jrr I:lrpo
(11)
dimana n('o(r} adalah jumlah atom antara r din r+b dari atom pusat, yang dirata-ratakan meliputi seluruh atom
2i J~ 2000
~
~
F~ ~
~
terhadap atom pusatnya. Pada prinsipnya ukuran konfigurasi L cukup besar sedemikianrupa sehingg3 tidak ada korelasiyang berpotongandengansel,sehingga g(r>U2) =1. Fungsidistribusi radial g(r) hanyadihitung untuk r
dangan
Dari refs. 3 clan refs. 8 terlihat bahwa X;2/2 dalam metode RMC ekivalen dengan U/kT dalam metode MMC. Hal yang membedakan antarapemodelan dengan metode Metropolis Monte Carlo (MMC) atau simulasi dinamika molekuler dengan metode Reverse Monte Carlo (RMC) adalah pada dua metode yang pertama diperlukan adanya potensial interatomik sedangkanpada metode RMC tidak diperlukan potensial interatomik untuk menentukan konfigurasi atom. Dalam banyak kasus metode MMC lebih menguntungkan, akan tetapi pada kenyataannya jarang sekali potensial secara kuantitatif dihasilkan kembali oleh data difraksi caira,l atau glass, sehingga apa yang dipertimbangkan merupakan suatu hal yang bersifat kualitatif clan belum menghasilkan model yang menjadi perhatian secara rinci dalam analisis struktur amorf [2] karena itu potensial interatomik tidak digunakan dalam analisis bahan amorf
menggunakan RMC. Simulasi MMC dengan model bola keras (hard spheres) merupakan satu metode simulasi Monte Carlo dimana fungsi energi potensial memiliki nilai yang sangat besar (00) untuk jarak terdekat antara dua atom clan bemilai 0 pada jarak yang lain. Hal ini akan menimbulkan gaya interatom yang bersifat diskrit (diskontinyu). Karena itu model ini disebut model bola keras (hard spheres model). Pada model ini gaya akan menggunakan gaya hanya hila mereka melakukan tumbukan. Tumbukan bersifat lenting sempuma. Dengan demikian simulasi RMC tanpa mengikut sertakan data eksperimen dapat dianggap sebagai simulasi MMC model bola keras karena pengenaanjarak terdekat atom adalah ekivalen dengan konsep model bola keras (hard spheres).
p~
AgP03 memiliki struktur acak (amort), dimanapuncak-puncakamorf ditimbulkan oleh interaksi antaratom P,O dan Ag. Padalangkahpertamadilakukan koordinasiantara atom P-P denganbilangankoordinasi 1. Pendekatanjarak atom terdekat (panjang ikatan minimum) ditentukan dari fungsi distribusi radial total yang diperolehdari transformasilangsungfaktor struktur eksperimental [2]. lnteraksi antara ini membentuk rantai P-P-P-P-P Langkah kedua menambahkan satu atom 0 ditambahkan pada pusat ikatan P-P. Maksimalkanproporsi koordinasilipat dua p-o dan O-P antara panjang ikatan minimum Garak pendekatan terdekat)dan 2 A.. Denganmenambahkansatu atom 0
~I
H~ H~
H~ CML., T~H~P.~.
diantara P akan membentuk rantai ikatan ..-O-P-O-P-OP-O.. .Tahap ketiga, tambahkan dua atom 0 untuk setiap atom P. Koordinasi untuk atom O-P daD P-O masingmasing adalah 1 daD 2. Penambahan atom 0 ini alan memberikan struktur jaringan satuan P04 yang membagi dua jenis atom 0, yaitu berjembatan oksigen dan tal berjembatan oksigen. Tahap akhir simulasi adalah mempertahankan koordinasi lipat empat dari P- 0 .Pada tahap ini jaringan dapat mengubah kesesuaiannya sampai pada kondisi yang cocok (conformation) akan tetapi tetap mempertahankanbentuk ikatannya (connectivity).
HASIL DAN PEMBAHASAN Simulasi ini menggunakan 3840 atom AgP03 yang terdiri dari 768 atom P, 768 atom Ag, daD 2304 atom O. Atom-atom tersebut diletakkan pactaposisi acak. Kerapatan atom p= 0,0725 atom! A 3 Atom tersebut berada di dalam kubus dengan sisi L= 18,7774 A Dalarn simulasi ini digunakan jarak terdekat atom rc(P-P}, rc(PO}, rc(O-O} masing masing adalah 2,15, 1,40 daD 2,15A. Kendala koordinasi yang digunakan adalah antara atom P daD 0 dengan bilangan koordinasi lipat 4 pacta jarak maksimum sebesar 4 A. Perbedaaan jarak r untuk menghitung fungsi distribusi radial g(r) adalah sebesar L1r = 0, I A. Perpindahan maksimum atom P, 0 dan Ag masing-masing sebesar loA Dari basil simulasi diperoleh jarak tetangga terdekat pertama daD kedua antar atom. Jarak daD puncak distribusi parsial ditampilkan pada Tabel I. Tabel1. Jarak pendekatan terdekat atom dalam molekul AgPO3 dibandingkan dengan yang dilakukan
oleh M. Tahrez dkk. [3]
j-m
r (A)
Puncakgj-m(r)
r (A) [3]
p-p
2,9
1,6587 1,1306
3,0
5.2
P-Pkedua
P-o
--
152 ,
P-Okedua
250 ,
0-0 O-Okedua
TATA KERJA
Gd A~O,~
I
r\ ~-~g
I
0~~5~e:ua
I
Ag-Ag
Ag-Agkedua
Ag-P Oari pola fungsi distribusi radial parsial atom PP, gp_p(r), memiliki karakteristik osilator teredam dengan puncak osilasi pertama terjadi pada r=2,9 A dengan puncak distribusi 1,6587 clan osilasi kedua pada r=5,2 A denganpuncak distribu'si 1,306 A-I. Untuk r> 8 A fungsi gp_p(r)mencapai kondisi stabil di gp_p(r)=1. Bentuk ini menunjukkan karakteristik bahan dimana bahan memiliki struktur lokal tetapi kerapatan lokal lebih rendah daTi kerapatan rata-rata bahan clan meningkat dengan laju
2CJ J~ 2000
59
~
F~ ~
~
p~
~ At}o.~
H~ H~
H-u ~ T~H~p.m.
yang lebih rendah dari r3. Struktur ini disebut struktur pecahan (fractal structure). Pola harnburan akan dipengaruhi oleh ikatan atom P-P, 0-0 daD O-Ag. lnteraksi P-Ag relatif kurang berpengaruh pada puncak arnorf. Garnbar I menunjukkan fungsi distribusi radial parsial yang dilakukan oleh I. Borjesson untuk interaksi atom P-P, P-O dan 0-0[3]. Garnbar 2 fungsi distribusi radial parsial dari hasil simulasi.
KESIMPULAN Dengantidak memasukkansatu parameterdata eksperimenpada data masukan maka programReverse Monte Carlo (RMC) memiliki kesamaandengansimulasi Metropolitan Monte Carlo (MMC) metode bola keras (hard spheres). Hal ini dapatdigunakankarenafungsi energi potensialpada jarak antara dua atom memiliki nilai tak berhingga dan diluar jarak tersebut energi potensialnya sarna dengan nol dengan demikian parametermasukanbukanlah potensial seperti metode konvensional MMC akan tetapi jarak pendekatan terdekatantar atom yang dinyatakandalam model bola keras(hard spheres). Dari simulasi fungsi distribusi radial parsial dapat disimpulkan bahwa puncak amorf dipengaruhi oleh ikatan atom P-P, P-O, 0-0 dan O-Ag. Bentuk fungsi distribusi radial P-P menunjukkanbahwa bahan mempunyaistruktur lokal dengan kerapatanlokal lebih rendahdari kerapatanrata-rata.
UCAP AN TERIl'.IA KASIH Penulis mengucapkan terima kasih kepada Dr. McGreevy alas kebaikannya memberikan perangkat lunak RMC sehingga simulasi ini dapat dilakukan. Penulsi mengucapkan terima kasih kepada Dr. Agus Purwanto alas diskusinya yang membangun. Penulis juga mengucapkan terima kepada Dr. Abarrul Ikram selaku kepala Balai Spektrometri yang telah memberi dukungan pada penelitian ini
DAFTARPUSTAKA [1]. R.L.MCGREEVY, M.A. HOWE AND J.D.WICKS, A General Purpose Reverse Monte Carlo Code, (1993), p.16 [2]. R.L.MCGREEVY, RMC:progress. problems and prospects, Nuclear Instruments and Methods in Physics ResearchA 354 (1995) 1-16. [3]. M.TACHEZ, R. MERCIER, J.P. MALUGANI AND P.CHIEUX, Structure Determination of AgPOj and (AgPOJJo.5(AgI)o.5 by neutron diffraction and Small Angle Neutron Scattering, Solid State Ionics 25 (1987) 263-270. [4]. L.BORJESSON, R.L.MCGREEVY, W.S HOWELLS, Fractal Aspect of Superionic Glasses from Reverse Monte Carlo simulations, Philosophical Magazine B, (1992). Vol.65, No.2, 261-271.
Gambar 1. Fungsi distribusi parsial gp_p(r), gp-o(r) dan go-o(r) untuk jaringan PO3 yang dihasilkan oleh HSMC pada kerapatan atom yang berhubungan dengan AgPO3 (gambar bawah).
60
~,
2gJ~ 2000
~
F~
~
~
PMui Gd1.4 A$PO, ~
1-1t:tMt 1-1~
1-1-uc...l.. T~H~P.~.
::r::t
v""--
..~.!
;:
--
..-;;--;;-;;;--;:
\,.1""
u
Gambar 2. Fungsi distribusi radial parsial AgPO3. Sebelah kiri, dari atas ke bawah masing-masing adalah gp_p(r), gP-o(r), gp-Ag(r). Sebelah kanan, dari atas ke bawah masing-masing adalah go-o(r), go-Ag(r), gAg-Ag(r)
TANYA-JAWAB Penanya: Sudinnan (P3IB -BATAN) 1. Apa perbedaan Metropolitan MonteCarlo (MMC) denganReverse MonteCarlo (RMC).
Jawaban: 1. Secaraprinsip perbedaanantaraMMC dan RMC adalahsbb. PadametodeMMC diperlukanfungsi energi potensialyang berfungsisebagaipemicu untuk mendapatkankonfigurasiacak dari posisi atom. Sedangkan padaRMC tidak diperlukanfungsienergipotensial.
~, Ke Menu Utama
61
2gJ~ 2000 Ke Daftar Isi