Đề tài Mô phỏng mẫu lớp bằng phương pháp shell model monte carlo (smmc)

ÑAÏI HOÏC QUOÁC GIA TP.HCM TRÖÔØNG ÑAÏI HOÏC KHOA HOÏC TÖÏ NHIEÂN KHOA VAÄT LYÙ CHUYEÂN NGAØNH VAÄT LYÙ HAÏT NHAÂN
 KHOÙA LUAÄN TOÁT NGHIEÄP Ñeà taøi: MOÂ PHOÛNG MAÃU LÔÙP BAÈNG PHÖÔNG PHAÙP SHELL MODEL MONTE CARLO (SMMC) SVTH : LYÙ TUÙ ANH CBHD : Th.S TRÖÔNG THÒ HOÀNG LOAN CN ÑAËNG NGUYEÂN PHÖÔNG CBPB : PGS.TS. CHAÂU VAÊN TAÏO TP. HOÀ CHÍ MINH - 2008 Lôøi caûm ôn “ÔÛ ñôøi khoâng coù con ñöôøng cuøng maø chæ coù nhöõng ranh giôùi ñieàu coát yeáu laø phaûi coù ñuû söùc maïnh ñeå vöôït qua nhöõng ranh giôùi ño ù.” ( Muøa laïc – Nguyeãn Tuaân) Thaät vaäy! Thöïc hieän khoùa luaän naøy coù theå noùi laø moät ranh giôùi trong cuoäc ñôøi cuûa moãi sinh vieân. Ñeå hoaøn thaønh noù, moãi chuùng ta phaûi coù “ñuû söùc maïnh”. Nhöng ñeå coù ñuû söùc maïnh, ngoaøi noå löïc cuûa baûn thaân mình, chính söï taän tuïy cuûa thaày coâ, söï yeâu thöông cuûa gia ñình, söï nhieät tình cuûa beø baïn môùi coù theå giuùp em coù “ ñuû söùc maïnh”. Vaø ñaàu tieân em xin göûi lôøi tri aân ñeán caùc thaày coâ Boä moân Vaät lyù Haït nhaân ñaëc bieät laø Thaïc só Tröông thò Hoàng Loan ñaõ höôùng daãn em hoaøn thaønh toát khoùa luaän naøy. Trong quaù trình hoaøn thaønh khoùa luaän, ñaõ vaáp phaûi raát nhieàu khoù khaên. Do ñoù, em xin caùm ôn ngöôøi ñaõ giuùp cho em raát nhieàu ñeå giaûi quyeát caùc khoù khaên ñoù: Cöû nhaân Ñaëng Nguyeân Phöông. Tieáp theo, em xin göûi lôøi tri aân vaø voâ cuøng bieát ôn ñeán ba maù ñaõ vaát vaû cöïc khoå ñeå em coù theå ñi hoïc vaø tröôûng thaønh nhö ngaøy hoâm nay. Neáu nhö khoâng coù söï vaát vaû, khoå cöïc aáy, em seõ khoâng coù ngaøy hoâm nay! Cuoái cuøng laø lôøi caûm ôn cuûa em ñoái vôùi taát caû ngöôøi baïn cuûa mình. Duø vui hay buoàn caùc baïn ñeàu luoân ôû beân mình, uûng hoä, chia seû cuøng mình. Vaø ñieàu naøy ñaõ goùp phaàn khoâng nhoû ñeå mình coù theå hoaøn thaønh khoùa luaän. Caùm ôn caùc baïn! 1 MUÏC LUÏC Trang MUÏC LUÏC .....................................................................................................................1 CAÙC KYÙ HIEÄU VIEÁT TAÉT ..........................................................................................2 DANH MUÏC CAÙC HÌNH VEÕ - ÑOÀ THÒ - BAÛNG BIEÅU ............................................3 LÔØI MÔÛ ÑAÀU................................................................................................................4 Chöông 1 TỔNG QUAN MAÃU LÔÙP ............................................................................6 1.1 Môû ñaàu ............................................................................................................6 1.2 Maãu lôùp ...........................................................................................................6 Chöông 2 CAÙC PHƯƠNG PHAÙP MOÂ PHOÛNG MONTE CARLO CHO HAÏT NHAÂN ......................................................................................................................................15 2.1 Môû ñaàu ..........................................................................................................15 2.2 Ñònh nghóa phöông phaùp Monte Carlo.........................................................15 2.3 Moâ phoûng Monte Carlo cho haït nhaân ..........................................................15 2.4 Moät soá phöông phaùp Monte Carlo ...............................................................17 2.5 Phöông phaùp SMMC (Shell Model Monte Carlo).......................................21 Chöông 3 MOÂ PHOÛNG MAÃU LÔÙP ............................................................................27 3.1 Môû ñaàu ..........................................................................................................27 3.2 Pheùp nghòch ñaûo ma traän..............................................................................27 3.3 Caùch tính cho maãu lôùp..................................................................................30 3.4 Chöông trình ..........................................................................................................43 3.5 Keát quaû vaø nhaän xeùt.............................................................................................44 KEÁT LUAÄN..................................................................................................................47 KIEÁN NGHÒ .................................................................................................................47 TAØI LIEÄU THAM KHAÛO ...........................................................................................48 Phuï luïc 1 ......................................................................................................................50 Phuï luïc 2 ......................................................................................................................62 Phuï luïc 3 ......................................................................................................................63 2 CAÙC KYÙ HIEÄU VIEÁT TAÉT AFMC Phöông phaùp Auxiliary Field Monte Carlo. DMC Phöông phaùp Monte Carlo khuyeách taùn. GCMC Phöông phaùp Grand Canonical Monte Carlo. GFMC Phöông phaùp Green’s Fucntion Monte Carlo. SMMC Phöông phaùp Shell Model Monte Carlo. PIMC Phöông phaùp Path Integral Monte Carlo. PMC Phöông phaùp Projector Monte Carlo. QMCD Phöông phaùp Quantum Monte Carlo Diagonalization. VMC Phöông phaùp bieán phaân Monte Carlo. 3 DANH MUÏC CAÙC HÌNH VEÕ - ÑOÀ THÒ - BAÛNG BIEÅU Caùc hình veõ Trang Hình 1.1 : Sô ñoà phaân boá caùc lôùp cuûa caùc nucleon trong haït nhaân.13 Hình 2.1 : Söï giao thoa cuûa caùc lôùp trong haït nhaân.18 Hinh2.2 : So saùnh ñoà thò naêng löôïng lieân keát theo soá nguyeân töû baèng caùc phöông phaùp khaùc nhau LDA, DMC, HF vaø thöïc nghieäm....19 Hình 2.3 : Minh hoïa phöông phaùp GFMC baèng moâ phoûng naêng löôïng ôû traïng thaùi cô baûn cuûa moät soá haït nhaân...20 Hình 2.4 : Soá yeáu toá ma traän vôùi caùc khoâng gian maãu khaùc nhau...26 Hình 3.1 : Minh hoïa nhöõng quyõ ñaïo cuûa sô ñoà maãu lôùp..32 Hinh 3.2 : Töông taùc dö trong maãu sô ñoà.33 Hình 3.3 : Nhöõng ñöôøng ñaúng trò cuûa Hamiltonian ñoái vôùi N=8, χ =0.25...40 Hình 3.4 : Nhöõng ñöôøng ñaúng trò Hamiltonian ñoái vôùi N=8, χ =2.541 HìnhPL2.1 : Giao dieän “Chöông trình tính naêng löôïng baèng SMMC”63 HìnhPL2.2: Giao dieän tính thöû naêng löôïng cuûa nhaân He8 vôùi soá voøng laëp N=1000 ...64 Caùc baûng Baûng 3.1 : So saùnh keát quaû thöïc nghieäm, SMMC, GFMC, vaø VMC...45 Caùc ñoà thò Ñoà thò 3.1 : So saùnh keát quaû giöõa thöïc nghieäm, SMMC, GFMC, vaø VMC...46 4 LÔØI MÔÛ ÑAÀU OÂng baø ta coù caâu: Nhaân voâ thập toaøn. Vaø liệu trong thế giới khoa học, treân bước đường chinh phục vaø laøm chủ thế giới ấy điều naøy coù ñuùng khoâng! Chuùng ta đều biết rằng trong thời đại hiện nay, với sự phaùt triển vượt trội vaø khoâng ngừng, dường như khoâng coù gì laø khoâng thể, nhaân loại ñaõ đạt được những thaønh tựu to lớn về khoa học kỹ thuật. Tuy nhieân, beân cạnh đoù, chuùng ta vẫn đang đối mặt những vấn đề nan giải đặc biệt laø trong lĩnh vực vật lyù hạt nhaân. Trong lĩnh vực naøy, điều khoù khăn nhất hiện nay laø xaây dựng một lyù thuyết hoaøn chỉnh về cấu truùc hạt nhaân. Trong thời gian qua ñaõ coù nhiều mẫu cấu truùc ra đời, tuy nhieân, chưa coù một lyù thuyết naøo giải thích toaøn diện vaø chính xaùc. Mỗi mẫu cấu truùc chỉ giải thích được một số tính chất của hạt nhaân hay noùi caùch khaùc laø chuùng bổ sung cho nhau. Trong đoù, Mẫu Lớp laø giải thích được nhiều hiện tượng hạt nhaân nhất. Vaø việc nghieân cứu mẫu lớp trở neân được quan taâm hơn. Chính vì thế, nhiều phương phaùp để nghieân cứu mẫu lớp đaõ ra đời vaø ngaøy caøng chính xaùc hơn. Một trong những phương phaùp đoù laø phương phaùp moâ phỏng Monte Carlo Mẫu Lớp (Shell Model Monte Carlo – SMMC). Phương phaùp Monte Carlo Mẫu Lớp xuất phaùt từ phương phaùp Monte Carlo kết hợp với baøi toaùn lượng tử hoùa lần 2 (caùc toaùn tử sinh – hủy hạt). Phương phaùp naøy cũng ñaõ đạt được những thaønh coâng ñaùng kể như giải cho caùc hệ nhiều hạt (C, O, N,), tính được năng lượng ở caùc trạng thaùi cơ bản vaø kích thích, tính được moment từ của hạt nhaân, Caùc kết quả thu được khaù phuø hợp với thực nghiệm so với caùc phương phaùp khaùc. Việc thực hiện ít rườm raø, ít phức tạp hơn caùc phương 5 phaùp hiện đại hiện nay như DMC, PIMC. Đoù cũng laø lyù do tại sao taùc giaû chọn phương phaùp naøy cho khoùa luận của mình. Khoùa luận đaõ duøng phương phaùp SMMC tính caùc mức năng lượng của hạt nhaân He4, He6, He8. Khoùa luaän goàm 3 chöông: • Chöông 1: Toång quan maãu lôùp. • Chöông 2: Caùc phöông phaùp moâ phoûng Monte Carlo cho haït nhaân. • Chöông 3: Moâ phoûng maãu lôùp. Do thôøi gian hoaøn thaønh khoùa luaän coù haïn cuõng nhö sai soùt laø ñieàu khoâng theå traùnh khoûi, taùc giaû mong nhaän ñöôïc söï caûm thoâng vaø goùp yù cuûa quyù thaày coâ giuùp hoaøn thieän hôn veà ñeà taøi naøy. 6 Chöông 1 TỔNG QUAN MAÃU LÔÙP 1.1 Môû ñaàu Cho ñeán nay, maëc duø khoa hoïc ñaõ raát phaùt trieån nhöng theá giôùi haït nhaân vaãn coøn laø moät theá giôùi ñaày bí aån ñoái vôùi chuùng ta. Moät trong nhöõng noã löïc cuûa caùc nhaø vaät lyù hieän nay laø xaây döïng moät lyù thuyeát hoaøn chænh ñeå giaûi thích toaøn dieän vaø ñuùng ñaén taát caû soá lieäu thöïc nghieäm. Song ñoù vaãn coøn laø moät giaác mô! Do ñoù, caùc maãu haït nhaân laàn löôït ra ñôøi: maãu töông taùc maïnh (ñaïi dieän laø maãu gioït chaát loûng), maãu caùc haït ñoäc laäp (ñieån hình laø maãu lôùp), maãu suy roäng (keát hôïp hai maãu treân). Vaø chöa coù maãu naøo giaûi thích toaøn dieän veà haït nhaân. Nhöng maãu lôùp vaãn laø moâ hình veà caáu truùc haït nhaân ñöôïc aùp duïng ñeå giaûi thích nhieàu hieän töôïng haït nhaân nhaát. Maãu lôùp cuõng chính laø noäi dung ñöôïc trình baøy trong phaàn naøy. 1.2 Maãu lôùp 1.2.1 Nhöõng bieåu hieän toàn taïi maãu lôùp trong haït nhaân. Noäi dung maãu lôùp 1.2.1.1 Nhöõng bieåu hieän toàn taïi cuûa maãu lôùp [1] Thöïc nghieäm cho thaáy naêng löôïng lieân keát cuûa caùc nhaân coù soá neutron vaø soá proton truøng vôùi caùc soá: 2, 8, 20, 28, 50, 82, 126 coù nhöõng tính chaát raát ñaëc bieät. Caùc nhaân naøy goïi laø nhaân magic. Chuùng coù caùc tính chaát sau: Naêng löôïng lieân keát cuûa caùc nhaân magic lôùn hôn so vôùi caùc nhaân beân caïnh. Thöïc nghieäm cho thaáy naêng löôïng lieân keát cuûa caùc nhaân coù soá proton hay neutron truøng vôùi 3, 9, 21, 29, 51, 83 vaø 127 laø ñaëc bieät nhoû. Töùc laø nhaân beân caïnh nhaân magic coù naêng löôïng lieân keát beù hôn nhieàu so vôùi nhaân magic. 7 Tính beàn vöõng cuûa nhaân magic theå hieän ôû ñoä giaøu cao cuûa chuùng trong töï nhieân. Caùc nhaân naøy toàn taïi nhieàu treân traùi ñaát. Tính beàn vöõng cao cuûa nhaân magic coøn theå hieän trong söï giaûm tieát dieän baét neutron. Xaùc suaát chieám neutron cuûa caùc nhaân magic laø raát beù. Nghóa laø caùc nhaân magic raát khoù phaûn öùng baét neutron khi bò chuøm neutron ñaäp vaøo. Moment töù cöïc ñieän cuûa caùc nhaân magic laø nhoû. Nghóa laø caùc nhaân magic coù caáu truùc ñoái xöùng caàu, vì raèng ñoä lôùn moment töù cöïc ñieän ñaùnh giaù ñoä bieán daïng cuûa nhaân khoûi tính ñoái xöùng caàu veà phöông dieän ñieän tích. Theo coâng thöùc baùn thöïc nghieäm Weiszacker, naêng löôïng cuûa α haït phaùt ra bôûi nhöõng haït nhaân phoùng xaï phaûi taêng theo baäc soá Z. Nhöng trong thöïc teá, ñaõ phaùt hieän coù tröôøng hôïp ngoaïi leä cho quy luaät naøy laø haït nhaân polonium coù Z=84 phaùt ra haït α coù naêng löôïng cao hôn naêng löôïng cuûa haït α ñöôïc phaùt ra cuûa haït nhaân ñi sau noù. Noùi chung, nhöõng haït α coù naêng löôïng cao nhaát ñöôïc phaùt ra bôûi nhöõng haït nhaân phoùng xaï coù N=128; Z=84; N=84 ñeå bieán ñoåi thaønh haït nhaân coù N=126; Z=82; N=82. 1.2.1.2 Noäi dung maãu lôùp [2] Noäi dung cuûa maãu lôùp bao goàm: • Moãi nucleon trong nhaân chuyeån ñoäng trong moät tröôøng theá trung bình töï hôïp coù tính ñoái xöùng caàu (taïo ra bôûi caùc nucleon coøn laïi). • Caùc nucleon trong nhaân ñöôïc taäp hôïp thaønh nhoùm theo caùc möùc naêng löôïng. Treân moãi möùc naêng löôïng chæ coù theå chöùa moät soá giôùi haïn caùc nucleon. • Caùc haït nhaân magic laø caùc haït nhaân coù soá nucleon chieám ñaày moät lôùp naøo ñoù. 1.2.2 Nhöõng pheùp tính vôùi maãu lôùp Muïc tieâu cuûa nhöõng tính toaùn trong maãu lôùp laø tìm naêng löôïng vaø phoå naêng löôïng cuûa haït nhaân. Tröôùc tieân, ta seõ tìm hieåu qua veà ñònh thöùc Slater[8]: 8 Trong cô hoïc löôïng töû, ñònh thöùc Slater laø bieåu thöùc mieâu taû haøm soùng cuûa heä fermion maø heä naøy mang tính phaûn xöùng vaø tuaân theo nguyeân lyù ngoaïi tröø Pauli töùc khi moãi caëp femion ñoåi choã heä seõ ñoåi daáu moät laàn. Noù ñöôïc mang teân cuûa ngöôøi khaùm phaù ra noù, John C.Slater. Ông ñaõ tìm ra ñònh thöùc Slater nhö laø phöông tieän dieãn ñaït tính phaûn xöùng cuûa haøm soùng baèng ma traän. Ñònh thöùc Slater phaùt xuaát töø vieäc ngoaïi suy haøm soùng cuûa moät taäp hôïp electron maø moãi haøm soùng ñöôïc coi nhö laø spin quyõ ñaïo χ(x) , x laø vò trí vaø spin cuûa electron ñôn leû vaø khoâng coù haøm soùng cuûa 2 electron cuøng spin quyõ ñaïo. (Nguyeân lyù ngoaïi tröø Pauli) 1.2.2.1 Ñònh thöùc Slater trong tröôøng hôïp 2 haït Caùch ñôn giaûn nhaát xaáp xæ haøm soùng cuûa heä nhieàu haït laø tích cuûa caùc haøm soùng cuûa caùc haït ñoäc laäp. Ñoái vôùi tröôøng hôïp 2 haït ta coù: )(x)χ(xχ)x,ψ(x 221121 = (1.1) Bieåu thöùc naøy ñöôïc goïi laø tích Hatree. Tuy nhieân, noù khoâng thoaû cho haït fermion bôûi haøm soùng treân khoâng phaûn xöùng. Moät haøm soùng phaûn xöùng ñöôïc bieåu dieãn moät caùch toaùn hoïc nhö sau: )x,ψ(x)x,ψ(x 1221 −= (1.2) Vì vaäy, tích Hatree khoâng thoaû nguyeân lyù ngoaïi tröø Pauli. Vaán ñeà naøy ñöôïc khaéc phuïc baèng caùch keát hôïp tuyeán tính caû 2 tích Hatree: )}(x)χ(xχ)(x)χ(x{χ 2 1)x,ψ(x 1221221121 −= )(xχ)(xχ )(xχ)(xχ 2 1 2221 1211= (1.3) Heä soá tröôùc ñònh thöùc coù ñöôïc töø ñieàu kieän chuaån hoaù. Haøm soùng naøy phaûn xöùng vaø hôn nöõa phaân bieät giöõa caùc fermion. Tuy nhieân, noù cuõng tieán veà 0 neáu hai haøm soùng hoaëc hai fermion gioáng nhau. Ñieàu naøy thoaû maõn nguyeân lyù ngoaïi tröø Pauli. 9 1.2.2.2 Ñònh thöùc Slater trong tröôøng hôïp toång quaùt Bieåu thöùc naøy coù theå ñöôïc toång quaùt hoaù ñoái vôùi soá fermion baát kyø baèng caùch vieát noù döôùi daïng ñònh thöùc. =)x,...,x,ψ(x N21 )(xχ)(xχ)(xχ )(xχ)(xχ)(xχ )(xχ)(xχ)(xχ N! 1 nnn2n1 2n2221 1n1211 L MMM L L (1.4) Keát hôïp tuyeán tính caùc tích Hatree trong tröôøng hôïp hai haït, deã daøng thaáy ñöôïc ñoù chính laø ñònh thöùc Slater khi N=2. Ta coù theå chaéc chaén raèng việc söû duïng ñònh thöùc Slater thoaû caùc haøm soùng treân moïi taäp hôïp, coøn caùc haøm soùng ñoái xöùng töï khaéc ñöôïc loaïi boû. Töông töï, caùc ñònh thöùc naøy cuõng thoaû nguyeân lyù Pauli. Thöïc vaäy, ñònh thöùc Slater seõ bieán maát neáu taäp hôïp { }iχ laø ñoäc laäp tuyeán tính. Ñaây chính laø tröôøng hôïp coù hai hay nhieàu hôn spin quyõ ñaïo gioáng nhau. Trong hoaù hoïc, ngöôøi ta dieãn ñaït ñieàu naøy nhö sau: trong cuøng moät traïng thaùi khoâng coù 2 electron cuøng spin quyõ ñaïo. Toång quaùt ñònh thöùc Slater ñöôïc tính baèng khai trieån Laplace. Theo toaùn hoïc, ñònh thöùc Slater laø moät tensor phaûn xöùng. Ñònh thöùc Slater ñôn ñöôïc duøng xaáp xæ haøm soùng electron trong thuyeát Hatree-Fock. Trong nhöõng thuyeát chính xaùc hôn, söï keát hôïp tuyeán tính caùc ñònh thöùc Slater laø ñieàu caàn thieát. Tieáp theo, ta seõ trình baøy caùch tìm naêng löôïng vaø phoå naêng löôïng trong haït nhaân. Ta coù theå laøm ñieàu naøy baèng caùch giaûi phöông trình Schrodinger hoaëc löôïng töû hoaù laàn ΙΙ. 10 1.2.2.3 Giaûi phöông trình Schrodinger (löôïng töû hoaù laàn Ι) Ta giaûi phöông trình Schrodinger nhö sau: Ta coù phöông trình: EφφHˆ = (1.5) Hamiltonian coù daïng: ( )∑ ∑ = < += A 1k A k k k 2 k ξ,ξW 2m pˆ Hˆ l l (1.6) Vôùi { } kk kkkk ipˆ τ,σ,rξ ∇−= ≡ h rrr (1.7) Trong đó: kξ : Vị trí hạt thứ k. kr r : Bán kính hạt thứ k; kσ r : Spin của hạt thứ k; kτ r : Iso spin của hạt thứ k. kpˆ : Xung lượng hạt thứ k ( )lξ,ξW k : Theá töông taùc theo vò trí giöõa haït thöù k vaø thöù l . Ta seõ bieán ñoåi toaùn töû treân thaønh 2 yeáu toá: tröôøng trung bình vaø töông taùc dö ( ) ( ) ( )∑ ∑ ∑ = < = ⎥⎦ ⎤⎢⎣ ⎡ −+⎥⎦ ⎤⎢⎣ ⎡ += A 1k A k A 1k kkk k 2 k ξVξ,ξWξV 2m pˆ Hˆ l l [6] (1.8) Vôùi giaû thieát caùc haït ñoäc laäp, choïn V phuø hôïp vaø loaïi boû töông taùc dö ta coù: ( )∑ = ⎥⎦ ⎤⎢⎣ ⎡ +≈ A 1k k k 2 k ξV 2m pˆHˆ (1.9) Töø ñaây, ta coù phöông trình Schrodinger ñoái vôùi moät haït: (k)φE(k))]φV(ξ 2m pˆ[ iiik k 2 k =+ (1.10) vôùi ( ) ( )kkkii τ,σ,rφkφ rrr≡ (1.11) 11 Do caùc haït ñoäc laäp ta coù: ( ) ( )∏ = = A 1k iiii kφA,1,2,Φ kA21 KK (1.12) Suy ra phöông trình Schrodinger cho heä nhieàu haït: ( ) ( )A,1,2,ΦEA,1,2,ΦHˆ A21kA21 iii A 1k iiii KK KK ⎟⎠ ⎞⎜⎝ ⎛= ∑ = (1.13) Caùc nucleon ñöôïc xem nhö laø caùc haït fermion do ñoù haøm soùng bieåu dieãn noù laø haøm soùng phaûn xöùng töùc ( )A,1,2,Φ A21 iii KK laø phaûn xöùng, vaø coù daïng ñònh thöùc Slater nhö sau: ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )Aφ2φ1φ Aφ2φ1φ Aφ2φ1φ A! 1A)(1,2,...,Φ AAA 222 111 A2...1 iii iii iii iii K MOMM K K = (1. 14) Coâng vieäc coøn laïi laø ta choïn theá V phuø hôïp vaø ñöa vaøo phöông trình giaûi. Choïn ( ) 2kks2k2k ζs.ζrmω2 1 ξV l rrl r lll −−= [6] (1.15) Luùc naøy: ≈Hˆ ∑ = ⎥⎦ ⎤⎢⎣ ⎡ −−+ A 1k 2 kkks 2 k 2 2 k ζs.ζrmω 2 1 2m pˆ lrl r lll (1.16) Vôùi ω : haèng soá dao ñoäng ñieàu hoaø lsζ : ñaëc tröng töông taùc spin_quyõ ñaïo llζ : ñaëc tröng töông taùc quyõ ñaïo_quyõ ñaïo ( 3 heä soá naøy ñöôïc choïn toát nhaát baèng xaáp xæ Hartree-Fock) 12 Vôùi theá treân, ta coù phoå naêng löôïng nhö sau: Hình 1.1[6]: Sô ñoà phaân boá caùc lôùp cuûa caùc nucleon trong haït nhaân Dao ñoäng ñieàu hoaø 13 Caùc thoâng soá [6]: MeV41Aω 3 1−≈h MeV20Aζ 3 2 2 s −≈hl 0.1MeVζ 2 ≈hll Ta ñöôïc caùc soá magic: 2, 8, 20, 28, 50, 82, 126, 1.2.2.4 Löôïng töû hoaù laàn ΙΙ[5] Goïi +αa vaø αa laàn löôït laø toaùn töû sinh vaø huyû haït ôû traïng thaùi α, coù daïng ma traän: ⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛= 00 10 αa ⎟⎟⎠ ⎞ ⎜⎜⎝ ⎛=+ 01 00 αa (1.17) Nhö ta ñaõ bieát moãi traïng thaùi haït nhaân ñöôïc ñaëc tröng bôûi caùc soá löôïng töû ljm (l: moment goùc quyõ ñaïo; j: moment goùc toaøn phaàn; m: hình chieáu cuûa j treân truïc z), ôû ñaây ta duøng α ñeå thay theá cho boä ljm. Moät soá tröôøng hôïp ta coøn xeùt caû isospin τ khi xeùt ñeán traïng thaùi cuûa haït töùc luùc naøy α ≡ ljmτ. Coù theå vieát Hamintonian döôùi daïng: 21 HˆHˆHˆ += (1.18) ∑ += α ααα1 aaεHˆ (1.19) ∑ ++= αβγδ γδβααβγδ2 aaaaV2 1Hˆ (1.20) αβγδV : yeáu toá ma traän khoâng taïo caëp cuûa töông taùc 2 haït. Ñeå baûo ñaûm tính baûo toaøn söï quay vaø caáu truùc lôùp, ngöôøi ta ñaõ vieát laïi Hamiltonian 2 haït nhö sau: ∑∑ ∑ += abcd J M JMJMJ2 (cd)Aˆ(ab)Aˆcd)(ab,V2 1Hˆ ∑∑ ∑ +++= abcd J M JMJM A J 1/2 cdab (cd)Aˆ(ab)Aˆcd)(ab,V)]δ)(1δ[(14 1 (1.21) Toång treân ñöôïc laáy theo taát caû quyõ ñaïo ñôn haït cuûa proton vaø nôtron (ñöôïc kyù hieäu bôûi a, b, c, d) vôùi caùc toaùn töû sinh vaø huyû haït nhö sau: ∑ +++++ ×−== ba baaabb mm JM jjmjmjbbaaJM ]a[aaJM)amjmj((ab)Aˆ (1.22) 14 ∑ ×== ba babbaa mm JM jjmjmjbbaaJM ]a[aaJM)amjmj((ab)Aˆ (1.23) Coøn cd)(ab,VJ laø yeáu toá ma traän moment goùc taïo caëp cuûa theá voâ höôùng )r,rV( 21 rr , ñöôïc xaùc ñònh nhö sau: [ ] [ ]JM2j1j21JM2j1jJ )r(ψ)r(ψ)r,rV()r(ψ)r(ψcd)(ab,V dcba rrrrrr ××= (1.24) Trong tröôøng hôïp 2 haït phaûn ñoái xöùng cd)(ab,VJ ≡ cd)(ab,VAJ ñöôïc xaùc ñònh