Khóa luận tốt nghiệp: Khảo sát năng lượng tương quan Positron – Electron trong phân tử đồng Oxit

Chia sẻ: Lavie Lavie | Ngày: | Loại File: PDF | Số trang:61

Thêm vào BST

Báo xấu

53
lượt xem 5
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Khóa luận tốt nghiệp: Khảo sát năng lượng tương quan Positron – Electron trong phân tử đồng Oxit gồm có 4 chương. Trong đó, chương 1 - Tổng quan về Positron - Electron, chương 2 - Phương pháp biến phân Monte Carlo lượng tử, chương 3 - Áp dụng phương pháp VMC cho phân tử đồng oxit khi có positron, chương 4 - Các kết quả tính toán.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Khóa luận tốt nghiệp: Khảo sát năng lượng tương quan Positron – Electron trong phân tử đồng Oxit

THƯ TRƯỜNG ĐẠI HỌC SƯ PHẠM THÀNH PHỐ HỒ CHÍ MINH VIỆN KHOA VẬT LÝ CHUYÊN NGÀNH VẬT LÝ HẠT NHÂN  KHOÁ LUẬN TỐT NGHIỆP Đề tài: KHẢO SÁT NĂNG LƯỢNG TƯƠNG QUAN POSITRON – ELECTRON TRONG PHÂN TỬ ĐỒNG OXIT CBHD: ThS. Trịnh Hoa Lăng SVTT: Phạm Thị Phú TP. HỒ CHÍ MINH – 2010
LÔØI CAÛM ÔN Sau khi hoaøn thaønh khoùa luaän toát nghieäp vôùi ñeà taøi "khaûo saùt naêng löôïng töông quan electron - positron trong phaân töû ñoàng oxit". Toâi xin chaân thaønh göûi lôøi caûm ôn ñeán: Thaày Trònh Hoa Laêng ngöôøi ñaõ chæ höôùng, cung caáp taøi lieäu tham khaûo vaø söûa chöõa sai soùt cho toâi trong suoát quaù trình thöïc hieän ñeà taøi. Toaøn theå caùc thaày coâ khoa vaät lyù tröôøng Ñaïi Hoïc Sö Phaïm Tp. Hoà Chí Minh, ñaõ truyeàn ñaït cho toâi nhöõng kieán thöùc quyù baùu trong suoát boán naêm hoïc taïi tröôøng. Ngöôøi thaân vaø taát caû caùc baïn ñaõ ñoäng vieân vaø giuùp ñôõ toâi hoaøn thaønh khoùa luaän. Moät laàn nöõa toâi xin caûm ôn taát caû moïi ngöôøi vaø xin nhaän nôi toâi loøng bieát ôn saâu saéc nhaát. Thaønh phoá Hoà Chí Minh thaùng 5 -2010 Sinh vieân thöïc hieän: Phaïm Thò Phuù - K32.
DANH MUÏC CAÙC KYÙ HIEÄU, CAÙC CHÖÕ VIEÁT TAÉT  Caùc kyù hieäu. n: soá löôïng töû chính. 2: toaùn töû Laplace. l: soá löôïng töû quyõ ñaïo. : delta Dirac. a0: baùn kính Bohr. T : haøm soùng thöû. me: khoái löôïng electron. : haøm soùng cuûa heä.  : soá PI. J: haøm Jastrow. EL: naêng löôïng cuïc boä. N: soá haït electron. Z: ñieän tích hieäu duïng. T: ñoäng naêng.  : haèng soá Planck. E: naêng löôïng. : maät ñoä. E ec  p : naêng löôïng töông quan H: Hamilton. electron – positron. : toaùn töû Gradien. V: theá naêng.  Caùc chöõ vieát taét. Ps: positronium. e+A: hệ liên kết positron - nguyên tử. : gamma. I: theá ion cuûa nguyeân töû CuO: phaân töû ñoàng oxit a.u: Ñôn vò nguyeân töû (atomic unit). e – p: electron – positron. VMC: Variational Monte Carlo.
CAÙC ÑÔN VÒ Caùc ñôn vò tính toaùn trong nguyeân töû Kyù Trong heä Ñaïi löôïng Trong heä SI hieäu ñôn vò nguyeân 5,291 772 108x10-11 Baùn kính Bohr a0 1 m Khoái löôïng electron me 9,109 3826x10-31 kg 1 1,602 176 53x10 -19 Ñieän tích electron. E 1 C 1,054 571 68x10 -34 J Haèng soá Planck  1 s Naêng löôïng tính theo 1 4,359 744 17x10 -18 J 1 Hartree Hartree (27,211 3845 eV) 1,0973731568525x Haèng soá Rydberg Ry 107/m. 0,5 (13,9056923 eV)
LỜI MỞ ĐẦU Positron laø phaûn haït cuûa electron ñöôïc khaùm phaù lyù thuyeát bôûi Paul Dirac naêm 1928, sau ñoù ñöôïc Carl D. Andersen quan saùt thöïc nghieäm naêm 1932. Vaø ngaøy 15 thaùng 3 naêm 1933, sau khi taïp chí khoa hoïc Myõ chuyeân ñeà vaät lyù xuaát baûn baøi "electron mang ñieän döông" bôûi Carl D. Andersen cuûa vieän coâng ngheä California, lòch söû thöïc nghieäm positron baét ñaàu. Trong thöïc nghieäm huyû positron, döïa vaøo ñaëc tính huyû electron - positron maø positron coù phạm vi öùng duïng raát lôùn nhö: Söû duïng positron ñeå phaùt hieän khuyeát taät trong vaät lieäu baèng phöông phaùp ño phoå thôøi gian soáng, CT (Computed Tomography) trong coâng nghieäp ñeå phaùt hieän loã hoång vaät lieäu. Trong y học, positron ñöôïc öùng duïng vaøo coâng ngheä maùy PET (Positron Emission Tomography) duøng phổ biến trong chuẩn đoán và theo dõi bệnh ung thư. Caùc keát quaû thöïc nghieäm ñöôïc ño giaùn tieáp thoâng qua thôøi gian soáng cuûa positron hoaëc xung löôïng huyû caëp electron – positron trong moâi tröôøng khaûo saùt. Vaø nhöõng keát quaû thöïc nghieäm naøy seõ ñöôïc giaûi thích chính xaùc hôn neáu chuùng ta xeùt ñeán töông taùc electron - positron. Töø nhöõng naêm 50, ñaõ coù raát nhieàu nhaø khoa hoïc aùp duïng caùc lyù thuyeát löôïng töû ñeå giaûi thích söï huyû positron trong moâi tröôøng chaát raén nhöng haàu nhö raát phöùc taïp vaø khoâng ñaït keát quaû mong muoán. Ngaøy nay, phöông phaùp bieán phaân Monte Carlo ñang ñöôïc aùp duïng roäng raõi cho caùc heä löôïng töû nhö nguyeân töû, phaân töû, vaät lyù chaát raén. Phöông phaùp naøy tính toaùn ñôn giaûn vaø giaûi quyeát toát caùc vaán ñeà cuûa theá giôùi vi moâ. Vì vaäy trong phaïm vi khoaù luaän naøy, toâi aùp duïng phöông phaùp bieán phaân Monte Carlo trong moâ hình Born - Oppenheimer ñeå khaûo saùt naêng löôïng töông quan electron - positron trong phaân töû ñoàng oxit . Noäi dung khoùa luaän goàm boán chöông: Chöông 1: Toång quan veà heä positron - electron. Chöông 2: Phöông phaùp bieán phaân Monte Carlo löôïng töû. Chöông 3: Aùp duïng phöông phaùp VMC cho phaân töû ñoàng oxit khi coù positron. Chöông 4: Caùc keát quaû tính toaùn.
CHÖÔNG 1 TOÅNG QUAN VEÀ HEÄ ELECTRON - POSITRON Electron laø moät trong nhöõng haït cô baûn caáu thaønh neân vaät chaát, coù phaûn haït laø positron. Positron mang ñieän döông coù cuøng khoái löôïng, cuøng spin, momen töø vaø ñoä lôùn ñieän tích vôùi electron. Neáu positron gaëp electron, chuùng coù xaùc suaát huyû caëp phaàn lôùn taïo ra hai tia gamma naêng löôïng 511 Kev. 1.1.Cô cheá lieân keát positron vôùi vaät chaát Caùc nghieân cöùu gaàn ñaây ñaõ chæ ra raèng moät soá nguyeân töû vaø ion nguyeân töû coù theå lieân keát vôùi moät positron nhö laø liti, heli, natri, canxi, magie, ñoàng, keõm, baïc… Soá nguyeân töû naøy ñöôïc tìm thaáy ngaøy caøng nhieàu. Vaø caáu truùc cuûa caùc heä naøy ñöôïc xaùc ñònh laø lôùn ñeå positron coù theå thaéng löïc ñaåy cuûa haït nhaân vaø lieân keát vôùi caùc electron hoaù trò. Moät soá heä nhö laø e+Be bao goàm moät positron chuyeån ñoäng quanh nguyeân töû beri phaân cöïc. Trong khi trong caùc heä khaùc e+Na, e+Cu, e+O…coù theå ñöôïc moâ taû toát nhaát baèng moät positronium chuyeån ñoäng quanh loõi mang ñieän Na+, Cu+, O+… Ngoaøi ra moät vaøi heä coù khaû naêng coù hai positron ñaõ ñöôïc tìm thaáy. 1.1.1.Positron lieân keát vôùi nguyeân töû Töông taùc giöõa nguyeân töû vaø positron laø ñaåy nhau taïi moïi ñieåm. Ñieàu naøy theå hieän söï khoù khaên trong lieân keát cuûa nguyeân töû vôùi positron. Song ñaùm maây ñieän tích electron cuûa nguyeân töû gaàn positron coù theå töï ñieàu chænh laøm cho töông taùc giöõa positron vaø nguyeân töû laø töông taùc huùt. Khi ñoù ñaùm maây electron khi naøy bò phaân cöïc. Vaø neáu theá huùt phaân cöïc lôùn hôn theá ñaåy cuûa haït nhaân ñoái vôùi positron thì traïng thaùi lieân keát toàn taïi. e+ Hình 1.1. Moâ hình positron lieân keát vôùi nguyeân töû, ñaùm maây ñieän tích cuaû electron ôû gaàn positron thay ñoåi Khi thế ion của nguyên tử I < 0,250 Hartree ( naêng löôïng lieân keát cuûa positronium) thì positron chæ lieân keát vôùi nguyeân töû khi naêng löôïng lieân keát cuûa positron vôùi nguyeân töû lôùn hôn
0,250 - I, ngöôïc laïi heä positron- nguyeân töû seõ taùch thaønh positronium vaø ion döông ( traïng thaùi nguyeân töû khi maát moät electron ) 1.1.2.Lieân keát cuûa positronium vôùi loõi nguyeân töû Caùc electron hoaù trò trong nguyeân töû lieân keát loûng leûo vôùi haït nhaân. Vaø khi thế ion của nguyên tử I < 0.250 Hartree thì moät trong caùc electron hoaù trò naøy seõ bò huùt vaøo positron, taïo thaønh moät positronium ( kí hieäu PS ). Positronium bao goàm moät electron vaø moät positron lieân keát vôùi nhau, noù töông töï nhö nguyeân töû hydro veà maët ñieän tích trong ñoù positron ñöôïc xem nhö proton. Positronium coù khoái löôïng ruùt goïn baèng nöûa khoái löôïng cuûa electron, coù baùn kính baèng hai laàn baùn kính cuûa nguyeân töû hidro. Do positron vaø electron ñeàu coù spin baèng ½ vaø moâ men quyõ ñaïo l=0 neân khi chuùng keát hôïp vôùi nhau coù theå hình thaønh hai traïng thaùi:  Traïng thaùi spin ñoái song coù momen goùc l = 0 vaø spin toaøn phaàn s = 0 ñöôïc goïi laø traïng thaùi singplet, hình thaønh paraPositronium (p-Ps) coù moät traïng thaùi laø 1S0.  Traïng thaùi spin song song coù l = 0 vaø s = 1, ñöôïc goïi laø traïng thaùi triplet, hình thaønh orthoPositronium coù ba traïng thaùi: 3S-1, 3S0, 3S1. Do vaäy, xaùc suaát hình thaønh orthopositronium laø ¾ coøn xaùc suaát hình thaønh para positronium laø ¼. e+ e- Hình 1.2. Moâ hình positron lieân keát vôùi nguyeân töû, Ps bò phaân cöïc do loõi mang ñieän tích döông. Positronium bò phaân cöïc do töông taùc Culong cuûa loõi nguyeân töû vaø coù theå daãn tôùi söï lieân keát cuûa PS vôùi loõi nguyeân töû. Ñeå Positronium lieân keát vôùi nguyeân töû thì naêng löôïng lieân keát cuûa positron vôùi nguyeân töû phaûi nhoû hôn 0,250 - I. 1.2.Söï huyû positron Positron laø phaûn haït cuûa electron, neân theo thuyeát phaûn haït trong cô hoïc löôïng töû thì khi positron gaëp electron seõ huyû caëp vôùi xaùc suaát naøo ñoù. Khi positron ñi vaøo moâi tröôøng vaät chaát seõ
va chaïm vôùi caùc electron, ion nuùt maïng vaø caùc phonon dao ñoäng maïng. Quaù trình va chaïm xaûy ra lieân tieáp laøm cho caùc positron nhanh choùng maát naêng löôïng vaø trôû thaønh positron nhieät. Quaù trình naøy ñöôïc goïi laø quaù trình nhieät hoaù positron. Thôøi gian nhieät hoaù phuï thuoäc vaøo naêng löôïng ban ñaàu cuûa positron. Theo Bergersen, Hautojarvi: positron coù naêng löôïng khoaûng 2 Mev coù thôøi gian nhieät hoaù nhoû hôn 20 ps. Trong khi thôøi gian soáng cuûa positron trong kim loaïi laø khoaûng 200ps. Neân positron bò nhieät hoaù tröôùc khi huyû caëp. Nhieät hoaù e+  1012 s Huyû  1021 e+ Khueách taùn Nguoàn e+ Hình 2 : Quaù trình huyû positron. Khi positron nhieät gaëp electron trong moâi tröôøng vaät chaát thì chuùng coù theå huyû caëp ngay laäp töùc, hoaëc chuùng keát hôïp vôùi electron hoaù trò cuûa nguyeân töû trong moâi tröôøng vaät chaát taïo thaønh PS. Neáu traïng thaùi spin cuûa caëp huyû laø traïng thaùi singlet ( s = 0 ) thì quaù trình huyû chuû yeáu laø huyû 2 . Coøn traïng thaùi spin cuûa caëp huyû laø traïng thaùi triplet thì söï huyû chuû yeáu laø huyû 3 . Neáu coù maët cuûa moät haït thöù ba laø electron hoaëc haït nhaân nguyeân töû cuûa moâi tröôøng thì coù theå xaûy ra quaù trình huyû 1 ( theo ñònh luaät baûo toaøn ñoäng löôïng ). Trong ñoù tröôøng hôïp huûy caëp sinh hai gamma laø chuû yeáu coøn tröôøng hôïp sinh ra ba gamma vaø moät gamma coù xaùc suaát raát nhoû neân coù theå boû qua.
CHÖÔNG 2 PHÖÔNG PHAÙP bieán phaân MONTE CARLO LÖÔÏNG TÖÛ Muïc ñích: trình baøy toång quan veà phöông phaùp bieán phaân Monte Carlo, caùch aùp duïng phöông phaùp naøy cho heä löôïng töû, moâ hình haøm Hamilton vaø haøm soùng thöû cho heä e - p. Caùc noäi dung chính:  Nguyeân lyù bieán phaân.  Phöông phaùp bieán phaân Monte Carlo.  Quy trình thuaät toaùn Monte Carlo löôïng töû.  AÙp duïng phöông phaùp VMC cho heä löôïng töû.  Xaáp xæ Hamilton.  Haøm soùng thöû. 2.1.Nguyeân lyù bieán phaân Vieäc tìm lôøi giaûi gaàn ñuùng cho phöông trình Schodinger cuûa heä nhieàu haït xuaát phaùt töø nguyeân lyù bieán phaân, ñöôïc phaùt bieåu nhö sau: Giaù trò trung bình cuûa Hamilton Hˆ ñöôïc tính vôùi haøm soùng thöû t khoâng bao giôø thaáp hôn giaù trò naêng löôïng traïng thaùi cô baûn 0 cuûa Hamilton Hˆ ñöôïc tính vôùi haøm soùng traïng thaùi cô baûn 0. Töø nguyeân lyù bieán phaân ta thaáy, luoân coù khaû naêng tìm ñöôïc moät giôùi haïn treân ñoái vôùi naêng löôïng traïng thaùi cô baûn. Vaø noù laø moät phöông phaùp xaáp xæ toát nhaát trong vieäc tìm kieám caùc naêng löôïng traïng thaùi cô baûn. Giaû söû haøm soùng thöû t laø haøm soùng traïng thaùi cô baûn cuûa Hamilton thì khi ñoù giaù trò trung bình cuûa Hamilton laø:  * E   (r) H  (r)dr T T (2.1) *   (r) (r)dr T T Vaø khai trieån haøm soùng thöû T laø söï keát hôïp tuyeán tính caùc haøm rieâng cuûa haøm hamilton coù daïng: T (r)  cn n (r) (2.2) n
Thay (2.2) vaøo (2.1) ta ñöôïc  * a m a n   m* ( r ) H  n ( r ) dr m, n E  (2.3)  a m* a n   m* (r ) n (r )dr m,n 2 a E n n  E  n 2 (2.4) a n n Trong ñoù En laø trò rieâng cuûa haøm rieâng  n vaø En   vôùi moïi n suy ra: E  0 (2.5) Phöông phaùp VMC tröïc tieáp aùp duïng nguyeân lyù naøy. Trong ñoù choïn haøm soùng thöû phuï thoäc vaøo caùc tham soá bieán phaân, roài thay ñoåi caùc tham soá bieán phaân naøy ñeå cöïc tieåu hoaù naêng löôïng trung bình. Neáu haøm soùng thöû ñöôïc choïn laø toát, coù ñuû baäc töï do bieán phaân thì keát quaû ñaït ñöôïc coù ñoä chính xaùc raát cao. 2.2.Phöông phaùp bieán phaân Monte Carlo (VMC) Trong vaät lyù chaát raén, phöông phaùp VMC thöôøng moâ taû moâ hình heä höõu haïn nhö moâ taû moät oâ moâ phoûng hôn laø moâ hình heä voâ haïn. Ôû ñaây chuùng ta choïn moâ hình haøm Hamilton vaø haøm soùng thöû phuø hôïp cho heä nhieàu haït, aùp duïng nguyeân taéc bieán phaân leân giaù trò trung bình cuûa haøm hamilton:  * H    (r ) H  (r )dr T T (2.6) *   (r )  ( r )dr T T Ñeå tìm naêng löôïng trung bình, ta phaûi thöïc hieän caùc pheùp tính tích phaân treân (2.6). Tuy nhieân vieäc tính tích phaân naøy laø raát phöùc taïp maø ta khoâng theå tính toaùn baèng giaûi tích toaùn hoïc. Caùc phöông phaùp tích phaân coå ñieån nhö phöông phaùp Gauss-Legendre khoâng coøn phuø hôïp cho heä nhieàu haït. Hieän nay thuaät toaùn lyù töôûng ñeå giaûi quyeát vaán ñeà naøy laø phöông phaùp VMC. Phöông phaùp VMC goàm coù 3 böôùc cô baûn: Böôùc 1: Xaây döïng haøm soùng thöû  T (r , ) cho heä goàm N haït ôû vò trí r = ( r1, r2,...,rN ) vaø phuï thuoäc vaøo boä tham soá  = (1, 2,..., N ) Böôùc 2: Xaùc ñònh giaù trò trung bình cuûa haøm Hamilton * H   T ( r ) H  T ( r )dr (2.7) *  T ( r )( r ) T ( r ) dr
Böôùc 3: Thay ñoåi khoaûng chaïy nhoû nhaát cuûa caùc tham soá bieán phaân trong thuaät toaùn roài quay trôû laïi böôùc moät. Sau khi thöïc hieän moät loaït caùc böôùc chaïy tham soá  , ta ghi nhaän giaù trò naêng löôïng vaø phöông sai töông öùng. Veõ ñoà thò bieåu dieãn söï phuï thuoäc cuûa naêng löôïng, sai soá töông ñoái vaøo caùc tham soá. Sau ñoù tìm giaù trò cöïc tieåu cuûa naêng löôïng vaø sai soá töông ñoái trong daûi caùc tham soá bieán phaân. Keát quaû thu ñöôïc phuï thuoäc vaøo daïng xaáp xæ haøm soùng vaø moâ hình haøm Hamilton. Do ñoù vieäc löïa choïn haøm soùng thöû vaø moâ hình hamilton cho oâ moâ phoûng laø raát quan troïng, quyeát ñònh möùc ñoä thaønh coâng cuûa vieäc moâ taû ñaëc tröng cuûa heä voâ haïn. 2.3. Quy trình thuaät toaùn monte carlo löôïng töû Thieát laäp nhöõng tham soá cuûa heä Thöïc hieän böôùc dòch Thieát laäp chuyeån traïng thaùi caân baèng Chaáp nhaän hoaëc loaïi boû böôùc dòch chuyeån Thöïc hieän böôùc dòch chuyeån Thöïc hieän tính Chaáp nhaän hoaëc loaïi boû toaùn böôùc dòch chuyeån Tích luõy vaø tính trung bình caùc ñaïi löôïng Nhöõng ñaëc tröng cuûa heä Hình 2.1. Sô ñoà quy trình thuaät toaùn Monte Carlo löôïng töû [4].
Quy trình thuaät toaùn Monte Carlo löôïng töû goàm hai böôùc chính:  Thieát laäp traïng thaùi caân baèng Metropolis: Thieát laäp caùc electron ban ñaàu ôû caùc vò trí ngaãu nhieân ñoàng nhaát trong oâ moâ phoûng. Sau ñoù dòch chuyeån, moãi böôùc dòch chuyeån seõ taïo ra caùc böôùc ngaãu nhieân maø chaáp nhaän hay loaïi boû tuaân theo thuaät toaùn Metropolis söû duïng haøm soùng thöû cho heä nhieàu haït. Thöïc hieän moät soá böôùc dòch chuyeån ñuû lôùn ñeå ñaûm baûo heä ñaõ ñaït ñöôïc traïng thaùi caân baèng vaø nhöõng electron ñöôïc phaân 2 boá ngaãu nhieân theo phaân boá T R  . Ñoái vôùi haàu heát nhöõng tính toaùn ñaõ thöïc hieän thì caàn söû duïng 25% cuûa toång soá böôùc dòch chuyeån Monte Carlo cuûa moãi moät electron ñeå taïo ra nhöõng vò trí caân baèng.  Thöïc hieän tính toaùn Monte Carlo: Khi heä ñaït ñöôïc traïng thaùi caân baèng thì caùc ñaïi löôïng quan taâm nhö naêng löôïng cuûa heä ñöôïc tích luyõ. Moãi moät laàn dòch chuyeån ñöôïc taïo ra thì caùc ñaïi löôïng quan taâm ñöôïc tính vaø ñöôïc coäng doàn vaøo quaù trình toång. Khi taïo ra ñuû nhöõng böôùc dòch chuyeån , giaù trò trung bình ñöôïc tính. 2.4.AÙp duïng phöông phaùp VMC cho heä löôïng töû. Ñeå aùp duïng phöông phaùp naøy, ñaàu tieân ta xaây döïng moâ hình cho heä löôïng töû, sau ñoù ta löïa choïn haøm soùng thöû phuø hôïp vôùi moâ hình ñoù. Tieáp ñeán ta tính xaùc suaát phaân boá haøm soùng: 2  T (r ) p(r )  2 (2.8)  t ( r ) dr Naêng löôïng trung bình ñöôïc tính theo coâng thöùc:  * E    ( r ) H  ( r ) dr T T (2.9) *   ( r ) ( r ) dr T T Khi söû duïng haøm soùng thö,û ta ñònh nghóa moät ñaïi löôïng môùi laø naêng löôïng cuïc boä EL:  H  T (r ) EL(r) = (2.10)  T (r ) 2  E    (r ) E dr T L (2.11) 2   (r ) dr T neân: E   p ( r ) EL ( r )dr (2.12)
Phöông trình (2.12) laø böôùc giaûi quyeát ñaàu tieân cuûa phöông phaùp VMC cho heä löôïng töû. Sau khi laøm goïn haøm naêng löôïng E L (r ) ta thöïc hieän treân maùy tính caùc böôùc sau: Caùc böôùc thöïc hieän chöông trình tính toaùn VMC Ñaàu tieân coá ñònh soá böôùc Monte Carlo vaø soá böôùc nhieät hoaù. Choïn giaù trò khôûi taïo cho r, caùc 2 tham soá bieán phaân vaø tính  T (r ) . Xaùc ñònh kích thöôùc böôùc nhaûy ñeå di chuyeån haït töø r ñeán vò trí môùi. Khôûi taïo cuûa EL vaø phöông sai töông öùng. Baét ñaàu tính toaùn Monte Carlo vôùi moät soá chu trình cho tröôùc Böôùc 1: Gieo soá ngaãu nhieân cho vector vò trí r. Böôùc 2: Thöïc hieän böôùc thöû Monte Carlo. Böôùc 3: Di chuyeån vò trí thöû: r’ = r +, vôùi  laø moät soá ngaãu nhieân gieo trong khoaûng [0, 1]. Böôùc 4: Laäp tæ soá p(r)/p(r’), neáu: p (r )   thì thieát laäp böôùc nhaûy môùi ngöôïc laïi ta döøng taïi vò trí naøy. p( r ' ) Böôùc 5: Neáu böôùc nhaûy ñöôïc thieát laäp ta gaùn: r = r’. Böôùc 6: Tính giaù EL vaø phöông sai öùng vôùi giaù trò cuûa r naøy. Khi pheùp thöû Monte-Carlo keát thuùc, ta tính giaù trò trung bình cuûa EL vaø ñoä leäch chuaån töông öùng. 2.5.Xaáp xæ Hamilton Trong heä e -p, caùc haït ñeàu mang ñieän neân chuùng töông taùc vôùi nhau theo nhieàu cô cheá: löïc Coulomb, töông taùc spin-quó ñaïo, töông taùc töø. Ñeå moâ taû caùc töông taùc phöùc taïp ñoù ta söû duïng tröôøng theá töông taùc caëp U (ri  rj ) vôùi i  j, khi ñoù haøm Hamilton coù daïng: N 2 2 N 1 H    i  V (ri )  U (ri  r j ) i 1 2m i 1 2 i j (2.13) Ñeå ñôn giaûn ta söû duïng heä qui chieáu nguyeân töû, trong ñoù :   m  1 vaø haøm Hamilton ñöôïc vieát laïi nhö sau: 1 N 2 N 1 N H   i  2 i 1   i 1 V ( r )   U (ri  r j ) 2 i 1 (2.14) Trong phöông phaùp bieán phaân Monte Carlo, ngöôøi ta thöôøng söû duïng xaáp xæ Born – Oppenheimer [8] thì Hamilton cho heä nhieàu haït goàm caùc electron, positron vaø caùc ion coù daïng:
1 N 1 Z e2 Zne2 Hˆ     i2   2p    n  2 i 1 2 n i ri  d n n rp  d n (2.15) 1 e2 e2 1 Z n Zme2   2 i j ri  r j  i r  r 2 n   m dn  dm i p j i mn ri, rp laàn löôït laø veùc tô toïa ñoä cuûa electron thöù i vaø positron. dn, dm laø veùc tô toïa ñoä cuûa ion thöù n vaø m. Zn, Zm laø ñieän tích cuûa ion thöù n vaø m. Soá haïng ñaàu tieân trong (2.15) laø toång ñoäng naêng cuûa N electron, soá haïng thöù hai laø ñoäng naêng cuûa positron; soá haïng thöù ba laø theá naêng töông taùc giöõa electron vaø ion, soá haïng thöù tö laø theá naêng töông taùc giöõa positron vaø ion; soá haïng thöù naêm laø theá naêng töông taùc giöõa electron – electron, soá haïng thöù saùu laø theá naêng töông taùc iöõa positron – electron vaø soá haïng cuoái cuøng laø theá naêng töông taùc giöõa ion – ion. 2.6.Haøm soùng thöû Heä e - p laø heä caùc haït fermion khoâng ñoàng nhaát, trong ñoù coù söï töông taùc giöõa caùc haït vôùi nhau. Löïc töông taùc giöõa electron vaø haït nhaân nguyeân töû tæ leä nghòch vôùi khoaûng caùch. Neân trong tinh theå ta chæ xeùt caùc nguyeân töû laân caän vôùi nguyeân töû chöùa electron maø ta ñang khaûo saùt. Khi ñoù haøm soùng toaøn phaàn cuûa heä electron-positron coù daïng nhö sau:     i ( ri ) p ( r p )  Je  e ( rij ) Je  p ( rip ) (2.16) i j ji Trong ñoù:  i (ri ) : haøm soùng cuûa electron thöù i trong oâ moâ phoûng. Haøm soùng naøy ñöôïc xaáp xæ nhö laø toå hôïp tuyeán tính caùc haøm soùng quyõ ñaïo:     i r    C nii r  d n n  (2.17)     i r  d n laø haøm soùng ñôn electron trong nguyeân töû thöù n vaø toång ñöôïc laáy treân taát caû caùc nguyeân töû trong oâ moâ phoûng. Cni laø haèng soá chuaån hoùa.  p ( r p ) laø haøm soùng cuûa positron, töông töï haøm soùng cuûa electron nhöng chæ khaùc ôû caùc tham soá trong haøm soùng.  Jee ( rij ), Je p (rip ) Coøn laø heä soá Jastrow theå hieän söï töông quan giöõa electron-electron, electron-positron.
2.6.1. Haøm soùng rieâng phaàn Haøm soùng  (r ) cuûa electron hay positron trong tröôøng nguyeân töû xuyeân taâm coù theå ñöôïc moâ taû baèng moät trong caùc daïng xaáp xæ sau:  Daïng Slater: Haøm soùng kieåu Slater laø haøm soùng thöôøng ñöôïc söû duïng trong nhöõng tính toaùn caáu truùc electron trong nguyeân töû, ñöôïc ñònh nghóa nhö sau:  ( r )  Nr n l e Zr (2.18) Trong ñoù: n>l+1 Töø ñieàu kieän chuaån hoaù haøm soùng ta ñöôïc: 2Z 2n 1 N (2.19) 2n !  Daïng Gauss : Haøm soùng kieåu Gauss ñöôïc ñònh nghóa nhö sau:  ( r )  Nr n  l exp(  r 2 ) (2.20) Töø ñieàu kieän chuaån hoaù ta coù: 2 n   n   3 1 2 2  2 N  (2.21)  2n    !  Daïng Laguerre : Xuaát phaùt töø ña thöùc Laguerre Lk(x) coù daïng: dk Lk ( x )  e x (xk e x ) (2.22) dx k Laáy ñaïo haøm ña thöùc Laguerre s laàn theo bieán x ta thu ñöôïc ña thöùc Laguerre suy roäng: ds Lsk ( x )  Lk ( x) dx s (2.23) Ña thöùc Laguerre suy roäng thoûa phöông trình Laguerre coù daïng: d 2 Lsk ( x ) dLsk ( x ) x  ( s  1  x )  (k  s ) Lsk ( x )  0 (2.24) dx 2 dx Haøm xuyeân taâm ñôn haït ñöôïc bieåu dieãn tröïc tieáp qua ña thöùc Laguerre nhö sau: r l 2 (n  l  1)!  na0  2r  2l 1  2r  nl (r )   2 3 / 2 e     (2.25) n a0 [(n  l )!]2  na  Ln  j  na   0  0
2 trong ñoù a 0  mZe 2 r  2 a0 Traïng thaùi 1s: 10 ( r )  e (2.26) a03 / 2 r 1 r  2a 0 Traïng thaùi 2p:  21 (r )  e (2.27) 2 6a03 a0 Haøm soùng Laguerre coù öu ñieåm hôn haøm soùng Slater vaø Gauss laø noù tröïc giao öùng vôùi moãi moät giaù l vôùi moïi n:  nl  n l   n,n ' ' ' (2.28) 2.6.2. Heä soá Jastrow Heä soá Jastrow thöôøng ñöôïc söû duïng trong phöông phaùp bieán phaân Monte Carlo coù daïng: J r   e J r  (2.29) Trong ñoù J(r) coù nhieàu daïng xaáp xæ khaùc nhau nhö xaáp xæ Page, xaáp xæ Mit, xaáp xæ taâm nguyeân töû, xaáp xæ Willamson… Vôùi moãi tröôøng hôïp ta aùp duïng xaáp xæ phuø hôïp töông öùng. Vaø trong giôùi haïn khoaù luaän naøy ta aùp duïng xaáp xæ Page coù daïng r J r   1  r (2.30) Vôùi , , laø caùc tham soá bieán phaân ñöôïc choïn ñeå thoaû ñieàu kieän chaën vaø r laø khoaûng caùch giöõa hai haït. Xaáp xæ Pade laø moät trong nhöõng xaáp xæ ñôn giaûn vaø phoå bieán nhaát.
CHÖÔNG 3 AÙP DUÏNG PHÖÔNG PHAÙP VMC CHO PHAÂN TÖÛ ÑOÀNG OXIT KHI COÙ POSITRON Muïc ñích: xaùc ñònh naêng löôïng töông quan electron - positron. Caùc noäi dung chính:  Xaây döïng haøm soùng thöû cho heä phaân töû CuO khi coù positron.  Xaây döïng haøm Hamilton.  Giaûi phöông trình Schrodinger.  Xaùc ñònh phöông sai cuûa naêng löôïng. Ñeå xaây döïng haøm soùng thöû cho heä tröôùc heát ta ñi khaûo saùt caáu hình vaø caáu truùc cuûa caùc nguyeân töû oxy vaø ñoàng. 3.1. Khaûo saùt nguyeân töû oxy Oxy laø moät nguyeân toá phi kim, coù tính oxy hoùa maïnh, thuoäc phaân nhoùm chính nhoùm VI vaø naèm ôû vò trí thöù taùm trong baûng heä thoáng tuaàn hoaøn. Neân nguyeân töû oxy coù taùm electron chuyeån ñoäng xung quanh haït nhaân vôùi caáu hình electron:1s22s22p4 Trong ñoù, caùc electron ôû phaân lôùp 1s vaø 2s bò laáp ñaày thöôøng khoâng tham gia vaøo caùc lieân keát hoaù hoïc. Neân trong moâ hình xaáp xæ, caùc electron naøy seõ keát hôïp vôùi haït nhaân taïo thaønh loõi nguyeân töû. Coøn boán electron thuoäc phaân lôùp 2p chöa ñöôïc laáp ñaày neân deã daøng tham gia lieân keát. Loõi nguyeân töû Phaân lôùp 2p Hình 3.1. Moâ hình nguyeân töû oxy, trong ñoù boán electron phaân lôùp 2p ñang chuyeån ñoäng xung quanh loõi nguyeân töû oxy. Ñeå tìm haøm soùng ñôn haït cuûa caùc electron chuyeån ñoäng quanh loõi, ta söû duïng xaáp xæ Slater coù daïng:  nl (r )  Nr n l exp(Zr) (3.1) Boán eletron thuoäc phaân lôùp 2p coù soá löôïng töû chính n=2 vaø soá löôïng töû quó ñaïo l=1 do ñoù haøm soùng ñôn haït coù daïng:  21  N O r exp(  Z O r ) (3.2)
Vôùi r laø khoaûng caùch töø electron ñang xeùt ñeán nguyeân töû oxy. Töø ñieàu kieän chuaån hoùa haøm soùng ta coù: 5 4Z O NO  (3.3) 3 3.2.Khaûo saùt nguyeân töû ñoàng Ñoàng laø moät nguyeân toá kim loaïi chuyeån tieáp, naèm ôû vò trí thöù 29 trong baûng heä thoáng tuaàn hoaøn neân nguyeân töû ñoàng coù 29 electron chuyeån ñoäng xung quanh haït nhaân vôùi caáu hình electron :1s22s22p63s23p64s13d10. Ngoaøi ra ñeå giaûi thích hoaù trò hai cuûa ñoàng ngöôøi ta coøn chaáp nhaän caáu hình thöù hai laø:1s22s22p63s23p64s23d9. Theo moâ hình xaáp xæ, caùc electron thuoäc ba lôùp beân trong keát hôïp vôùi haït nhaân taïo thaønh loõi nguyeân töû. Nhö vaäy trong nguyeân töû Cu chæ xeùt hai electron trong phaân lôùp 4s chuyeån ñoäng xung quanh loõi nguyeân töû. Loõi nguyeân töû Cu Electron phaân lôùp 4s Hình 3.2. Moâ hình nguyeân töû ñoàng, hai electron phaân lôùp 4s ñang chuyeån ñoäng xung quanh loõi nguyeân töû ñoàng. AÙp duïng xaáp xæ Slater cho hai electron phaân lôùp 4s vôùi n = 4, l = 0 ta ñöôïc  1Cu  N Cu r 4 exp(  Z Cu r ) (3.4) Töø ñieàu kieän chuaån hoùa haøm soùng ta coù : 9 4ZCu NCu  (3.5) 315 3.3. Moâ hình phaân töû ñoàng oxit (CuO) Trong phaân töû CuO thì haøm soùng ñôn haït cuûa caùc electron trôû neân phöùc taïp hôn do chuùng chòu taùc duïng cuûa hai tröôøng löïc ñöôïc taïo bôûi hai loõi cuûa ñoàng vaø oxy. Trong ñoù caùc electron laø ñoàng nhaát neân khoâng theå phaân bieät ñöôïc electron naøo laø cuûa ñoàng hay cuûa oxy. Song chuùng ta
coù theå tìm ñöôïc xaùc suaát tìm thaáy haït nhôø caùc haøm soùng ñôn haït vôùi moâ hình phaân töû CuO ñöôïc moâ taû treân hình 3.3 y ei riO riCu ri O Cu d1 x d2 z Hình 3.3. Moâ hình phaân töû CuO, caùc electron ñang chuyeån ñoäng quanh hai loõi nguyeân töû ñoàng vaø oxy. Trong ñoù ta gaén heä truïc Descarses trong khoâng gian 3 chieàu vaøo moâ hình vôùi goác toaï ñoä naèm giöõa hai loõi nguyeân töû. Thì theo moâ hình treân ta coù: Toïa ñoä cuûa loõi nguyeân töû oxy laø (-d1, 0, 0). Toïa ñoä cuûa loõi nguyeân töû ñoàng laø (d2, 0, 0). Toïa ñoä vector vò trí cuûa electron thöù i: ri  ( xi , yi , zi ) (3.6) Toïa ñoä vector khoaûng caùch giöõa electron thöù i vaø haït nhaân oxy: riO  ( xi  d1 , yi , zi ) (3.7) Toïa ñoä vector khoaûng caùch giöõa electron thöù i vaø haït nhaân ñoàng: riCu  ( xi  d 2 , yi , zi ) (3.8) Toaï ñoâ vector khoaûng caùch giöõa electron thöù i vaø electron thöù j: rij  ( xi  x j , yi  y j , zi  z j ) (3.9)
Vaø haøm soùng ñôn haït cuûa caùc electron laø toång hai haøm soùng theo hai tröôøng loõi nguyeân töû Cu vaø O. Suy ra haøm soùng cuûa electron treân quó ñaïo cuûa nguyeân töû oxy coù daïng:  iO  N O riO exp(Z O r )  N O r iCu exp( Z Cu r ) iO iCu (3.10) Vaø haøm soùng cuûa electron treân quó ñaïo cuûa nguyeân töû ñoàng coù daïng:  iCu  N Cu riO4 exp(  Z O r )  N Cu riCu iO 4 exp(  Z Cu r ) iCu (3.11) Khi ñoù haøm soùng cuûa heä caùc haït electron trong phaân töû ñoàng oxit laø tích cuûa caùc haøm soùng ñôn haït vaø heä soá Jastrow: 2 4 6    iCu  iO  e e ij (3.12) i 1 i 1 i 1 j  i 1  ije  e Vôùi laø heä soá töông quan electron - electron, theo xaáp xæ Pade coù daïng: rij  ije  e  exp (3.13) 1  rij vôùi ,  laø caùc haèng soá bieán phaân. Vaäy haøm soùng cuûa heä laø 2     N Cu r14O exp( Z O r1O )  N Cu r14Cu exp(  Z Cu r1Cu )  i 1 4 6 rij (3.14)  N i 1 O riO exp(  Z O riO )  N O r iCu exp(  Z Cu riCu )  exp i 1 1  rij i j Töø phöông trình (3.14) ta thaáy, Haøm soùng cuûa heä trong moâ hình phaân töû CuO ngoaøi chöùa caùc bieán vector vò trí coøn chöùa caùc tham soá bieán phaân ZO, ZCu, , . 3.4. Moâ hình positron trong phaân töû CuO Khi ñöa positron vaøo tröôøng phaân töû CuO, heä taêng theâm moât haït mang ñieän döông. Khi ñoù positron chuyeån ñoäng trong tröôøng taïo bôûi loõi cuûa hai nguyeân töû ñoàng vaø oxy, vôùi haøm soùng ñôn haït coù daïng soùng s laø toång hôïp caùc haøm soùng ñôn haït taïo bôûi hai tröôøng löïc Cu vaø O.