Nghiên cứu một số phản ứng hạt nhân cần thiết cho thiên văn học

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ KHOA HỌC VÀ CÔNG NGHỆ VIỆN NĂNG LƯỢNG NGUYÊN TỬ VIỆT NAM ------
------ Nguyễn Ngọc Duy NGHIÊN CỨU MỘT SỐ PHẢN ỨNG HẠT NHÂN CẦN THIẾT CHO THIÊN VĂN HỌC Chuyên ngành: Vật lý nguyên tử và hạt nhân. Mã số: 62 44 05 01 TÓM TẮT LUẬN ÁN TIẾN SĨ VẬT LÝ Hà Nội – 2013 Công trình được hoàn thành tại: - Đại học Tổng hợp Tokyo, Nhật Bản. - Viện Năng lượng Nguyên tử Việt Nam. Người hướng dẫn khoa học: PGS. TS. Lê Hồng Khiêm PGS. TS. Vương Hữu Tấn Phản biện 1:................................................................. Phản biện 2:................................................................. Phản biện 3:................................................................. Luận án sẽ được bảo vệ trước Hội đồng chấm luận án cấp Viện họp tại: . vào hồi ........., ngày ...... tháng ............ năm ... Có thể tìm hiểu luận án tại: - Thư viện Quốc gia Việt Nam. - Viện Năng lượng nguyên tử Việt Nam. 1 Mở đầu Vật lý hạt nhân đóng vai trò quan trọng trong sự cải biến thế giới. Việc nghiên cứu vật lý hạt nhân mang lại nhiều ứng dụng trong nhiều lĩnh vực, chẳng hạn như công nghiệp, nông nghiệp, y học,Hơn nữa, vật lý hạt nhân còn là chiếc chìa khoá để nghiên cứu vũ trụ bao la. Những mô hình của sự hình thành và phát triển sao đã dự đoán về sự tồn tại của các quá trình hạt nhân đang diễn ra trong vũ trụ. Khảo sát những phản ứng hạt nhân cần thiết cho thiên văn học có ý nghĩa rất quan trọng không chỉ trong lĩnh vực thiên văn mà còn cả đối với lĩnh vực cấu trúc hạt nhân. Ngoài ra, những hạt nhân không bền trong môi trường sao được cho là sẽ giúp chúng ta có những kiến thức mới hơn về cấu trúc hạt nhân. Với những thiết bị hiện đại dành cho nghiên cứu hạt nhân hiện nay, chúng ta hoàn toàn có thể tiến hành việc khảo sát các phản ứng xảy ra trong vũ trụ tại các phòng thí nghiệm trên mặt đất. Đây là một điều kiện thuận lợi cho sự khám phá sâu hơn nữa về cấu trúc hạt nhân các hạt không bền và nghiên cứu một cách hiệu quả những phản ứng hạt nhân trong những chuỗi phản ứng của các quá trình tổng hợp nguyên tố trên các sao. Trong đó, một số phản ứng có vai trò hết sức quan trọng trong tiến trình phát triển của sao, cũng như liên quan mật thiết đến những sự bất thường trong quan sát thiên văn. Chúng tôi đặc biệt quan tâm đến phản ứng 22Mg(α,p)25Al vì nó có ý nghĩa quan trọng trong việc nghiên cứu cấu trúc hạt giàu proton 26Si trong vùng năng lượng trên ngưỡng alpha (Ethr = 9.164 MeV). Dữ liệu hạt nhân trong vùng năng lượng này dường như còn bỏ trống, trong khi đối với hạt nhân gương của nó, 26Mg, mật độ mức trong vùng năng lương tương ứng tương đối cao. Ngoài ra, phản ứng này cũng là một trong những mối kết nối quan trọng trong quá trình rp-process tại vị trí hạt nhân 22Mg trong quá trình tổng hợp nguyên tố trong vũ trụ. Suất phản ứng của phản ứng này góp phần giải thích sự bất thường trong việc dò tìm tia gamma năng lượng 1.275 MeV và vấn đề Ne-E hiện nay. Bên cạnh đó, thế “chờ” của 22Mg cũng được làm sáng tỏ dựa trên những kết quả nghiên cứu về phản ứng này và phản ứng 22Mg(p,γ)23Al. 2 Có hai phương pháp cơ bản để khảo sát suất phản ứng: phương pháp sử dụng thông số các mức lượng tử của hạt nhân, có năng lượng tương ứng với vùng nhiệt độ cần khảo sát trong môi trường các sao và phương pháp tính trực tiếp suất phản ứng từ tiết diện phản ứng của phản ứng đó. Chúng tôi đã tiến hành thiết kế đo đạc để thu các mức năng lượng của hạt không bền 26Si và tính toán suất phản ứng của 22Mg(α,p)25Al thông qua các mức cộng hưởng của 26Si trong thực nghiệm tán xạ hạt không bền 22Mg lên bia khí alpha. Trong khuôn khổ luận án này, chúng tôi trình bày những kết quả thu được từ 22Mg(α,α)22Mg cho hạt nhân compound 26Si và tính suất phản ứng của 22Mg(α,p)25Al. Trên thế giới chỉ có hai nghiên cứu về hạt nhân 26Si trên ngưỡng alpha. Nghiên cứu thứ nhất là phân rã beta từ hạt nhân 26P và nghiên cứu thứ hai là phản ứng 28Si(p,t)26Si. Tuy nhiên, các mức năng lượng trong những nghiên cứu này không xa ngưỡng alpha, và do đó, không thoả mãn được vùng nhiệt độ cao (T > 1 GK) trong môi trường Supernova và X-ray Burst. Cho đến trước luận án này, chưa có một công trình nào đo đạc trực tiếp tương tác 22Mg+α. Mục tiêu của luận án - Về vật lý hạt nhân: Xác định các mức năng lượng mới của 26Si trên ngưỡng alpha bao gồm các trạng thái lượng tử spins và chẵn lẻ. Từ những thông tin này, cấu trúc nhóm của hạt giàu proton 26Si sẽ được khảo sát. - Về vật lý thiên văn: Dựa vào các trạng thái kích thích của 26Si từ tán xạ 22Mg+α, xác định suất phản ứng của 22Mg(α,p)25Al, cùng với phản ứng 22Mg(p,γ)23Al và phân rã beta của 22Mg để khảo sát “điểm chờ” 22Mg. Từ đó, đánh giá nguyên nhân của những sự bất thường trong quan sát thiên văn về tia gamma 1.275 MeV và vấn đề Ne-E hiện nay. Thừa số S-factor cần cho việc xác định suất phản ứng không cộng hưởng cũng được tính toán. - Sản xuất chùm hạt không bền 22Mg: Sản xuất chùm hạt không bền 22Mg có độ sạch cao, cường độ lớn và năng lượng thoả mãn vùng trên ngưỡng alpha của 26Si và trong vùng nhiệt độ T = 0.5 - 3 GK. - Kĩ thuật detector: Cải tiến active-gas-target detector để đo đạc đồng thời hạt nặng của chùm hạt tới và các hạt nhẹ bay ra sau phản ứng. Đây là 3 tính chất rất quan trọng và cần thiết đối các loại detector ghi đo phản ứng theo phương pháp động học ngược với kĩ thuật bia dày. Kết quả mới của luận án Về mặt vật lý, luận án đã ghi nhận được sáu mức năng lượng trên ngưỡng alpha của 26Si. Trong đó, có ba mức thấp trùng khớp với hai nghiên cứu trước đó và ba mức mới lần đầu tiên được ghi nhận. Trạng thái spin và chẵn-lẻ của sáu mức này đã được xác định trong luận án này bằng việc làm khớp theo phương pháp R-matrix, dựa trên hàm kích thích tiết diện cộng hưởng của phản ứng 22Mg(α,α)22Mg, những công trình trước đây chưa xác định được trạng thái lượng tử cho các mức này. Đối với thiên văn học, suất phản ứng cộng hưởng của 22Mg(α,p)25Al trong vùng nhiệt độ sao T9 = 0.5 - 2.5 GK được xác định từ trạng thái cộng hưởng trong 22Mg(α,α)22Mg. Bên cạnh đó, hàm thừa số thiên văn S-factor phụ thuộc vào năng lượng cũng được xác định. Thừa số này quan trọng trong việc ngoại suy suất phản ứng không cộng hưởng của 22Mg(α,p)25Al. Về mặt kĩ thuật thực nghiệm, detector khí loại mới dùng để đo đạc các phản ứng thiên văn đã được thiết kế và đưa vào sử dụng. Ngoài ra, chùm hạt không bền 22Mg thoả mãn điều kiện cường độ cao cho việc đo đạc tán xạ lên alpha được tạo thành công. Trong đó, số hạt 22Mg được tạo ra không tỉ lệ tuyến tính với cường độ chùm hạt tới của phản ứng tạo 22Mg. Ứng dụng của kết quả nghiên cứu Kết quả trong nghiên cứu này có thể được sử dụng để khảo sát cấu trúc nhóm trong hạt nhân 26Si trên ngưỡng alpha 9.164 MeV và góp phần giải thích một số hiện tượng bất thường trong quan sát thiên văn, chẳng hạn như tia gamma 1.275 MeV; thế chờ của 22Mg; Những kết quả thu được trong nghiên cứu này đóng góp vào sơ đồ các mức năng lượng chưa từng được ghi nhận của 26Si. Ngoài ra, các mức này góp phần đánh giá trạng thái kích thích của hạt nhân gương của nó, 26Mg vì trạng thái spin và chẵn-lẻ của các mức kích thích của hạt nhân gương này vẫn còn thiếu rất nhiều. Bên cạnh đó, suất phản ứng của phản ứng 22Mg(α,p)25Al có thể được dùng để so sánh sự cạnh tranh của quá trình bắt proton và phân rã beta của 22Mg. Ngoài ra, 4 kết quả suất phản ứng thu được từ thực nghiệm này cũng góp phần đánh giá việc khảo sát suất phản ứng theo lý thuyết thống kê của Hauser-Feshbach. Kĩ thuật tạo chùm hạt không bền 22Mg có thể được áp dụng đối với việc tạo các chùm hạt không bền cho các thực nghiệm khác. Sự bất thường trong sự phụ thuộc vào cường độ chùm vào của cường độ chùm hạt ra cần được khắc phục. Về kĩ thuật phát triển detector, loại detector khí hoạt động (active-gas- target detector) đã được thiết kế để có thể dùng cho các phản ứng trong bia khí theo cơ chế động học ngược. Nội dung của luận án Bố cục luận án, được trình bày bằng tiếng Anh, với nội dung như sau: - Phần mở đầu. - Chương 1: Những kiến thức tổng quan về vật lý hạt nhân, vật lý thiên văn và mục tiêu của luận án. Những chuỗi phản ứng hạt nhân xảy ra trong quá trình hình thành và phát triển sao được khảo sát về mặt lý thuyết theo phương pháp suy luận từ mô hình hình thành sao. Từ đó, xác định một số phản ứng quan trọng để giải thích các hiện tượng thiên văn, trong đó có 22Mg(α,p)25Al và 22Mg(p,γ)23Al. Lý thuyết khảo sát suất phản ứng trong môi trường các sao cũng được trình bày trong chương này. - Chương 2: Đề cập đến thực nghiệm đo trực tiếp 22Mg(α,α)22Mg. Trong chương này, phương pháp và bố trí thực nghiệm được trình bày một cách chi tiết. Kết quả việc tạo chùm hạt không bền 22Mg cũng được chỉ ra trong phần này. - Chương 3: Trình bày việc xử lí số liệu thực nghiệm, những kết quả thu được từ thực nghiệm về hàm tiết diện vi phân theo năng lượng, các mức cộng hưởng của 26Si từ phản ứng 22Mg(α,α)22Mg, suất phản ứng của 22Mg(α,p)25Al so sánh với các phản ứng khác để khảo sát những vấn đề thiên văn. - Kết luận và kiến nghị: Những kết luận quan trọng của luận án và một số vấn đề nghiên cứu tiếp theo. 5 Chương 1. Giới thiệu Trong chương này, chúng tôi đề cập đến những kiến thức cần thiết, quan trọng của vật lý hạt nhân liên quan đến vật lý thiên văn. Từ đó, chúng ta có thể nhận thấy vị trí, vai trò của phản ứng được nghiên cứu của luận án này. 1.1. Nguồn gốc vật chất trong vũ trụ Trình bày những giả thiết và suy luận đưa đến việc xác định vật chất trong vũ trụ được cấu tạo bởi các hạt nucleon, hạt nhân và nguyên tử. 1.2. Quá trình tổng hợp nguyên tố trên các sao Trình bày quá trình tổng hợp nguyên tố trên các sao sau vụ nổ Big Bang diễn ra theo bốn quá trình chính: quá trình cháy hydrogen, quá trình cháy helium, quá trình tổng hợp nguyên tố đến Fe và quá trình tổng hợp nguyên tố nặng hơn Fe. 1.2.1. Quá trình cháy Hydrogen Tất cả quá trình tổng hợp bốn proton thành helium đều được gọi là quá trình cháy hydrogen. Trong đó, có ba quá trình chính là chuỗi proton-proton (pp), chu trình CNO và chu trình NeNa-MgAl. Mỗi quá trình tổng hợp phụ thuộc vào mật độ hạt, nhiệt độ và các hạt nhân xúc tác. Trong đó, vai trò và tầm quan trọng trong các vấn đề thiên văn của các phản ứng ( )22 23,Ne p Naγ , ( )22 22Na Neβ υ+ , ( )22 25,Mg p Alα và 25 26( , )Mg p Alγ được chỉ ra trong phần này. 1.2.2. Quá trình cháy helium Trình bày những nội dung cơ bản trong quá trình cháy Helium. Việc nghiên cứu quá trình cháy helium rất hữu ích cho việc giải thích nguyên nhân của sự không tồn tại hạt nhân bền có số khối A = 5, A = 8 và A = 6 - 11 trong tự nhiên có độ phổ biến thấp. 1.2.3. Quá trình tổng hợp nguyên tố đến Fe Trình bày quá trình cháy carbon, cháy oxygen và cháy silicon. 1.2.4. Quá trình tổng hợp hạt nhân nặng trên Fe Trình bày các quá trình: s- process, r -process và p-process. 1.3. Siêu tân tinh loại II Mục này trình bày những kiến thức cơ bản về siêu tân tinh loại II. Đây là giai đoạn kết thúc của một ngôi sao và bắt đầu hình thành sao mới. 6 1.4. X-ray burst Trong phần này, hiện tượng bùng phát tia X (X-ray Burst) được trình bày. Trong môi trường này, suất phản ứng của các phản ứng được cho là tương đối cao. Ở giai đoạn cuối của X-ray Burst, tất cả các hạt nhân nhẹ đều được biến đổi thành các nguyên tố nặng. 1.5. Tầm quan trọng của số liệu hạt nhân 26Si và phản ứng 22Mg(α,α)22Mg Mục này trình bày lý do của việc nghiên cứu thực nghiệm phản ứng 22Mg+α. Đây cũng chính là mục tiêu chính của đề tài. 1.5.1. Suất phản ứng của 22Mg(α,p)25Al Phần này chỉ ra vị trí và vai trò quan trọng của phản ứng 22Mg(α,p)25Al trong việc giải thích vấn đề Ne-E, quan sát tia gamma 1.275 MeV, làm sáng tỏ thế chờ của 22Mg, như trong hình 1.1 Hình 1.1. Điểm chờ 22Mg và phản ứng 22Mg(α,p)25Al trong quá trình cháy hydrogen thông qua chu trình NeNa-MgAl. 1.5.2. Sự phân bố 26Al trong thiên hà Mục này trình bày tầm quan trọng của phân rã beta của 26Si về 26Al. 26Al tồn tại ở trạng thái cơ bản và trạng thái đồng phân. Nếu 26Si phân rã về 26Al ở trạng thái cơ bản thì sẽ ghi nhận được tia gamma 1.809 MeV từ 26Mg ở trạng thái kích thích. Độ phổ biến của 26Al giúp xác định tuổi thiên hà. 1.5.3. Suất phản ứng của 25Al(p,γ)26Si Mục này trình bày tầm quan trọng của phản ứng 25Al(p,γ)26Si mà suất phản ứng này có thể được xác định bởi thông tin 26Si. 7 1.5.4. Cấu trúc hạt nhân 26Si trên ngưỡng alpha Vai trò và tầm quan trọng của phản ứng 22Mg(α,α)22Mg trong việc nghiên cứu các mức cộng hưởng trong hạt nhân 26Si trên ngưỡng alpha được trình bày trong phần này. 1.6. Mục đích của luận án - Khảo sát về mặt lý thuyết một số phản ứng quan trọng liên quan đến một số hiện tượng thiên văn được chỉ ra, chẳng hạn như phản ứng α+α, 8Be+α, 25Mg(p,γ)26Al, 25Al(p,γ)26Si, 22Mg(α,p)25Al . - Đo đạc thực nghiệm 22Mg+α để ghi nhận số liệu của 26Si trên ngưỡng alpha (9.164 MeV) và từ đó xác định suất phản ứng 22Mg(α,p)25Al, so sánh với suất phản ứng 22Mg(p,γ)23Al và suất phân rã beta của 22Mg để đánh giá hiện tượng thiên văn về vấn đề Ne-E, tia gamma 1.275 MeV, thế chờ của 22Mg. - Đánh giá tính toán suất phản ứng theo mẫu thống kê của Hauser- Feshbach. 1.7. Suất phản ứng trong môi trường sao Việc xác định suất phản ứng được thực hiện dựa trên những lý thuyết sau đây. Suất phản ứng NA vσ của một cặp hạt được đo với đơn vị cm3s-1 đối với phản ứng a A b B+ → + được xác định như sau: ( ) ( ) 0 .A a AN v N N v v v dvσ φ σ ∞ = ∫ (1.1) Trong đó, Na và NA là mật độ các hạt a và A. ( )vσ là tiết diện của phản ứng tương tứng với vận tốc v, năng lượng E. Khi đó: ( ) 1 32 2 0 8 1 exp ,A EN v E E dE kT kT σ σ piµ ∞       = −           ∫ (1.2) với khối lượng rút gọn a A a A m m m m µ = + . 1.7.1. Phản ứng không cộng hưởng Xác suất xuyên rào Pl của hệ phản ứng được xác định bởi: 2 2 ,l l l krP F G = + (1.3) 8 với k là vector sóng, r là khoảng cách phân tách giữa hai hạt, Fl và Gl là nghiệm của hàm Coulomb. Theo cơ học lượng tử chúng ta có: 1/2 2 2 1/2 2(2 ) exp .a Al µ Z Z e P E pi  ∝ −    
(1.4) Biến đổi gần đúng xác suất xuyên rào ta thu được: ( )exp 2 ,lP piη= − (1.5) với η là tham số Sommerfeld xác định bởi: 2 , a AZ Z e v η =
(1.6a) ( )1/21 22 31.29 ,Z Z Eµpiη = (1.6b) E tính theo keV, khối lượng rút gọn µ tính theo amu. Bên cạnh đó, cơ học lượng tử chỉ ra rằng tiết diện phản ứng σ tỉ lệ với bước sóng de-Broglie / p=
như sau: 2 2 1 .p Eσ pi pi − −∝ ∝ ∝
(1.7) Do đó, tiết diện phản ứng sẽ có dạng: ( ) ( ) ( )exp 2 . S E E E σ piη= − (1.8) Trong đó, S(E) được gọi là S-factor, là tham số thiên văn học. Từ (1.2) và (1.8), suất phản ứng được xác định bởi (1.9): ( ) 1/2 3/2 1/2 0 8 1 exp .A E bN v S E dE kT kT E σ piµ ∞       = − −           ∫ (1.9) Trong đó, ( )1/21/21 20.989 .b Z Z MeVµ= (1.10) Kết quả tổng hợp của hai hạng tử trong hàm e-mũ của (1.9) tạo ra một đỉnh phân bố, gọi là đỉnh Gamow. Tiết diện phản ứng trong vùng năng lượng này sẽ trội hơn hẳn so với những vùng khác. Tại đỉnh, suất phản ứng tính bởi: ( ) 1/2 3/2 0 1/2 0 8 1 exp .A E bN v S E dE kT kT E σ piµ ∞       = − −           ∫ (1.11) Với đỉnh: ( ) 2/3 1/32 2 2 0 1 2 61.22 ,2 bkTE Z Z T keVµ = =    (1.12) và độ rộng đỉnh: ( ) ( )1/61/2 2 2 50 1 2 64 0.749 .3 E kT Z Z T keVµ∆ = = (1.13) 9 1.7.2. Suất phản ứng hạt nhân cộng hưởng Cơ chế phản ứng hạt nhân cộng hưởng được cho là xảy ra theo hai bước: hình thành hạt hợp phần compound và phân rã về trạng thái cơ bản. Hình 1.2. Phản ứng cộng hưởng với cơ chế hạt nhân hợp phần. Trong phản ứng cộng hưởng, tiết diện phản ứng đối với các mức năng lượng E được xác định bởi công thức Breit-Wigner: ( ) ( ) 2 22 , 2 a b BW RE E σ pi ω Γ Γ = Γ− +
(1.14) Γa, Γb và Γ là độ rộng riêng phần và toàn phần của kênh vào và kênh ra. ER là năng lượng hạt tới. Thừa số thống kê spin ω có thể tính bởi (1.15): ( ) ( ) 2 1 , 2 1 2 1a A J J J ω + = + + + (1.15) J, Ja và JA là spin của trạng thái cộng hưởng, spin của hạt tới và của bia. Kết hợp những biểu thức trên, suất phản ứng được viết lại như sau: ( ) 1 32 2 0 8 1 exp .A BW EN v E E dE kT kT σ σ piµ ∞       = −           ∫ (1.16) 1.7.2.1. Cộng hưởng hẹp Trong trường hợp Γ/ER < 20% , trong vùng cộng hưởng năng lượng thay đổi rất nhỏ khi đó suất phản ứng sẽ là: ( ) ( )3/29 6 6 116058.08 10 exp iA itot i E N v T T σ µ ωγ−−   = × −    ∑ [cm3.mol-1.s-1] (1.17) Ei có đơn vị MeV, µ có đơn vị amu, T6 có đơn vị MegaKelvin (MK) và hàm lực cộng hưởng được cho bởi: . a bωγ ω Γ Γ = Γ (1.18) 10 1.7.2.2. Cộng hưởng rộng Trong trường hợp Γ /ER > 20% , cộng hưởng xảy ra trong khoảng năng lượng rộng, tiết diện phản ứng cộng hưởng σ(E) tại giá trị năng lượng E xung quanh ER với độ rộng ΓR = Γ(ER) có thể được tính bởi (1.19): ( )2 2 2 / 2( ) ( )( ) ,( ) ( ) ( - ) ( / 2) Ra b R a R b R R E E E E E E E σ σ ΓΓ Γ = Γ Γ + Γ (1.19) Γ và σR là độ rộng toàn phần tương ứng với năng lượng E và tiết diện phản ứng tương ứng với năng lượng cộng hưởng ER. Độ rộng riêng phần của các hạt tương ứng với moment lượng tử quỹ đạo l được tính bởi: ( ) 1/2 22 2 , ,l l R l n E P E E R θ µ   Γ =    
(1.20) 2 lθ là độ rộng rút gọn của trạng thái tương ứng với l và nó có thể được ghi nhận từ thực nghiệm hoặc tính toán thông qua độ rộng Γi và giới hạn Wigner Wi của các kênh của phản ứng, với: 2 ,il iW θ Γ = (1.21) và giới hạn Wigner: ( )2 2 3 , , l n i n P E R W Rµ =
(1.22) với ( )1/3 1/31.45n a AR A A= + . Suất phản ứng cộng hưởng rộng được xác định bởi (1.2) và (1.19). 1.8. Phương pháp R-matrix Nội dung cơ bản của phương pháp R-matrix được trình bày trong phần này. Trong đó, các thông số lượng tử của hạt nhân compound là các yếu tố ma trận R-matrix được xác định từ việc giải phương trình Shrodinger. Luận án này sử dụng code AZURE để làm khớp các giá trị thực nghiệm để thu nhận các trạng thái lượng tử trong phản ứng hạt nhân hợp phần. Chương 2. Thực nghiệm phản ứng 22Mg + α 2.1. Phương pháp thực nghiệm 2.1.1. Xác định vùng năng lượng khảo sát Phản ứng 22Mg + α đã được thiết kế để đo đạc trực tiếp tiết diện và năng lượng của phản ứng 22Mg(α,α)22Mg và 22Mg(α,p)25Al nhằm thu nhận thông tin về trạng thái lượng tử của 26Si trong vùng năng lượng E > Eαthr và 11 suất phản ứng (α,p) trong vùng nhiệt độ sao T9 = 0.5 – 2.5 GK. Cửa sổ Gamow ứng với vùng nhiệt độ môi trường sao được chỉ ra trong bảng 2.1. Bảng 2.1. Cửa sổ Gamow ứng với vùng nhiệt độ và năng lượng quan tâm Nhiệt độ (GK) Đỉnh Gamow (MeV) Độ rộng Gamow (MeV) Năng lượng 26Si (MeV) 0.5 0.960 0.469 9.890 - 10.359 1.5 2.000 1.174 10.575 - 11.749 2.5 2.808 1.796 11.074 - 12.870 2.1.2. Cơ chế phản ứng động học ngược với bia dày Để đảm bảo việc khảo sát năng lượng một cách liên tục trong vùng nhiệt độ T9 = 0.5 – 2.5 GK, chúng ta cần tiến hành thực nghiệm với kỹ thuật bia dày. Hình 2.1 chỉ ra tính chất của kỹ thuật bia dày khi nghiên cứu phản ứng hạt nhân. Hình 2.1. Phản ứng với kỹ thuật bia dày theo cơ chế động học ngược. 2.1.3. Hệ đo CRIB Thực nghiệm 22Mg + α đã được tiến hành tại hệ phổ kế CRIB của đại học tổng hợp Tokyo, đặt tại RIKEN, Nhật Bản. Hình 2.2 là sơ đồ cấu tạo của CRIB (Center for Nuclear Study - CNS - low energy Radioactive Ion Beam) có thể tạo ra các chùm hạt thứ cấp không bền có năng lượng E < 10 MeV/u từ những hạt nhân bền. Các hạt nhân bền (từ 4He đến 40Ca), được tách ra từ nguồn ion ERC và được gia tốc bởi máy gia tốc cyclotron AVF (K = 70) với năng lượng 3 - 10 MeV/u và cường độ từ 100 enA (electron- nanoAmpere ) đến vài eµA (electron-microAmpere). Chùm hạt sơ cấp bắn vào bia sơ cấp F0 để tạo ra chùm hạt không bền thứ cấp, chùm hạt không bền này sẽ được hệ phổ kế lọc hai lần bởi các cặp cực điện, cực từ và Wien filter. Sau khi ra khỏi Wien filter, chùm hạt không bền cần thiết cho phản ứng có độ sạch cao sẽ đi vào buồng bia phản ứng F3 cần khảo sát. 12 Hình 2.2. Sơ đồ hệ phổ kế CRIB của đại học Tokyo Thực nghiệm phản ứng 22Mg + α được thiết kế dựa trên các nguyên lý làm việc của hệ phổ kế nêu trên. Trong đó, chùm hạt sơ cấp được gia tốc từ cyclotron AVF bắn vào bia sơ cấp 3He, được giữ ở nhiệt độ nitrogen lỏng (90K) tại F0, tạo ra chùm 22Mg được lọc bởi các cực điện từ và Wien filter sau đó sẽ bắn vào bia thứ cấp là khí He+CO2(10%) tại F3. Hệ thống detector ghi nhận sản phẩm α và p từ phản ứng 22Mg + α được đặt trong buồng F3. 2.1.4. Detector đo hạt 2.1.4.1. Detector PPAC PPAC là loại detector dùng để xác định vị trí của chùm hạt. Đây là loại detector khí. Nguyên lý hoạt động dựa trên sự ion hoá của chất khí, sinh ra electron và ion di chuyển về các cặp điện cực. Tín hiệu từ các điện cực cho phép chúng ta xác định được vị trí của chùm hạt dùng để bắn vào bia khí trong buồng thực nghiệm F3. Hình 2.3 chỉ ra cấu trúc của PPAC sử dụng trong thực nghiệm này. (a) (b) Hình 2.3. Cấu trúc detector PPAC. Hình 2.4. Hệ detector gồm 3 lớp (a) và 2 lớp (b). 2.1.4.2. Thiết kế hệ thống Silicon-detector Trong thực nghiệm này, chúng tôi sử dụng các detector silicon dạng strip đơn (SSD) được thiết kế bởi công ty Hamamatsu. Trong detector, 13 mặt trước gồm 08 strip và mặt sau gồm 01 pad. Trong phép đo phản ứng 22Mg + α có hai kênh ra cần quan tâm là: 22Mg(α,α)22Mg và 22Mg(α,p)25Al. Do quãng chạy của proton và alpha trong silicon khác nhau, vì vậy cần có sự phân biệt hạt theo phương pháp ∆E-E, nên chúng tôi thiết kế hệ detector gồm 2 hoặc 3 lớp SSD. Hình 2.4 mô tả thiết kế của hệ thống detector silicon dùng trong thực nghiệm này. 2.1.4.3. Thiết kế active-gas-target detector GEM-MSTPC Chúng tôi phát triển loại detector mới, dùng GEM foil và khí sử dụng cho detector cũng chính là bia của phản ứng, gọi là active-gas-target detector GEM-MSTPC. Loại detector này gồm hai vùng chính là vùng dịch chuyển của electron và vùng khuếch đại tín hiệu, như trong hình 2.5. Hình 2.5. Cấu tạo chính của active-target detector GEM-MSTPC. Hình 2.6. Cấu trúc vùng khuếch đại tỉ lệ của detector GEM-MSTPC. Vùng dịch chuyển của electron có kích thước 450 mm x 220 mm x 100 mm (dài x rộng x cao), vùng này chứa bia khí He+CO2(10%). Vùng khuếch đại tỉ lệ được thiết kế thành hai vùng để đo hạt nặng và hạt nhẹ, gọi là vùng low- gain và vùng high-gain. Vùng low-gain sẽ đo chùm hạt nặng tới, có kích thước 200 mm x 110 mm. Vùng high-gain gồm một vùng theo hướng và hai vùng hai bên chùm hạt tới. Trong mỗi vùng, bên dưới các GEM foil, 14 các điện cực anode, gọi là pad readout, được thiết kế dạng tam giác vuông, cạnh 4.2mm, để ghi điện tích của electron từ sự ion hoá ở vùng khí, tín hiệu các hạt sẽ được ghi nhận từ các anode này. Hình 2.6 mô tả cấu trúc các vùng khuếch đại tỉ lệ đo hạt nặng và hạt nhẹ. 2.2. Bố trí thực nghiệm 2.2.1. Bố trí đo đạc phản ứng 22Mg + α Thực nghiệm được tiến hành tại buồng F3 của hệ phổ kế CRIB. Bố trí thực nghiệm được trình bày trong hình 2.7. Hình 2.7. Bố trí hệ thống detector cho thực nghiệm 2.2.2. Hệ thống điện tử và ghi nhận số liệu Hệ điện tử ghi nhận số liệu gồm hai phần cơ bản: tín hiệu thời gian dùng cho trigger và tín hiệu năng lượng dùng để xác định các thông số vật lý của phản ứng, như chỉ ra trong hình 2.8. Hình 2.8. Sơ đồ khối điện tử dùng trong thực nghiệm đo 22Mg+α. Để Trigger được thiết lập bởi tín hiệu từ hai detector PPAC, tín hiệu trên SSD và tín hiệu chia nhỏ của chùm hạt tới được ghi bởi hai PPAC: Trigger = Beam ⊕ SSDs + Beam/n + SSDs OR Beam = PPACa ∧ PPACb 15 Vì thời gian dịch chuyển của electron từ sự ion hoá trong vùng khí về anode chậm hơn so với tín hiệu của alpha và proton trên SSD, do đó chúng ta cần thiết lập cửa sổ trùng phùng cho tín hiệu từ GEM-SMSTPC và SSD. Cửa sổ trùng phùng được xác lập như trình bày trong hình 2.9. Hình 2.9. Sơ đồ thời gian của cửa sổ trùng phùng. 2.3. Ghi đo số liệu Trong thực nghiệm này, hệ thống ghi đo số liệu sử dụng VME và CAMAC. Số liệu từ tín hiệu detector được ghi nhận theo từng sự kiện. Sau đó, thông tin từng sự kiện được lưu dưới dạng tập tin máy tính*.ridf, dạng thập lục phân. Sau khi giải mã file *.ridf , chúng ta có thể lưu số liệu dạng *.dat hoặc *.root. Sau đó, các số liệu cần quan tâm được tách lọc. 2.4. Tạo chùm hạt không bền 22Mg Để phép đo phản ứng 22Mg + α thành công, việc sản xuất chùm hạt 22Mg (T1/2 = 2.875 s) có năng lượng đủ lớn, có cường độ và độ sạch cao rất quan trọng. Phòng thí nghiệm CRIB có thể khắc phục những khó khăn trước đây để tạo chùm 22Mg cho phép đo 22Mg + α với năng lượng trên ngưỡng alpha. Kết quả việc tạo chùm hạt này đã được tác giả công bố trên tạp chí Nuclear Instruments and Methods section A. 2.4.1. Xác định phản ứng tạo chùm hạt 22Mg Sau khi tính toán tiết diện từ các phản ứng khác nhau kết quả cho thấy phản ứng 3He(20Ne,22Mg)n có tiết diện tốt nhất. Do đó, chùm hạt 22Mg được tạo ra theo phản ứng này. 2.4.2. Tạo chùm 22Mg Chùm hạt 22Mg được sản xuất bằng phản ứng 3He(20Ne,22Mg)n. Trong đó, bia khí 3He được làm lạnh ở nhiệt độ nitrogen lỏng (90 K) với áp suất 170 torr. Chùm 22Mg thu được có cường độ 1200 pps với năng lượng tại bia là 0.85 MeV/u. Kết quả tạo chùm hạt không bền 22Mg cho phản ứng 22Mg + α được chỉ ra trong luận án. 16 Chương 3. Phân tích số liệu thực nghiệm và kết quả. Phần này trình bày chi tiết quá trình phân tích số liệu và các kết quả ghi nhận được. Các mức năng lượng của 26Si được ghi nhận và áp dụng tính toán suất phản ứng 22Mg(α,p)25Al. 3.1. Chuẩn năng lượng Việc chuẩn năng lượng được tiến hành cho hệ detector SSD và detector GEM-MSTPC. Đường chuẩn năng lượng có dạng hàm tuyến tính bậc nhất như sau: E = a × Ch + b . (3.1) Các thông số chuẩn năng lượng cho tất cả 96 strip của các SSD được trình bày trong phần phụ lục A của luận án. Kết quả chuẩn tại hình 3.1 và 3.2. Việc chuẩn năng lượng cho GEM-MSTPC được tiến hành theo hai bước: đo độ mất năng lượng của hạt ∆E bằng cách sử dụng nguồn alpha và chuẩn độ mất năng lượng thực tế đo được với tính toán lý thuyết theo công thức bán thực nghiệm của Zigler. Độ mất năng lượng này liên quan đến điện tích Q đo được trực tiếp từ GEM-MSTPC xác định bởi hàm tuyến tính sau: ∆E = αQ + β (3.2) Hình 3.1. Chuẩn năng lượng cho vùng high gain. Hình 3.2. Chuẩn năng lượng cho vùng low gain. 3.2. Nhận biết các hạt trong thực nghiệm 3.2.1. Nhận biết các hạt trong chùm hạt không bền Chúng ta cần nhận biết các hạt thuộc phản ứng 22Mg+α mà chúng ta quan tâm. Do trong chùm nặng hạt tới bia khí chứa hai ion 20Ne và 21Na, nên chúng ta cần phân biệt 22Mg với chúng. Việc phân biệt các hạt này dựa trên độ mất năng lượng dE/dx, đường cong Bragg của chúng đo bởi detector khí, như trong hình 3.3. 17 3.2.2. Nhận biết các hạt nhẹ sinh ra từ phản ứng Đối với các hạt ra, alpha và proton, chúng được phân biệt bằng phương pháp ∆E – E. Trong đó, ∆E được đo bởi high-gain của GEM- MSTPC, trong khi E đo bởi SSD cộng với ∆E. Kết quả nhận biết alpha và proton từ phản ứng được trình bày trong hình 3.4. Hình 3.3. Độ mất năng lượng của 22Mg, 21Na và 20Ne được đo bởi GEM-MSTPC. Hình 3.4. Phân biệt alpha và proton theo phương pháp ∆E-E. 3.2.3. Hiệu chỉnh độ mất năng lượng Năng lượng phản ứng chính là năng lượng còn lại trên pad của GEM-MSTPC. mà ở đó xảy ra phản ứng. Việc đo trực tiếp một cách liên tục độ mất năng lượng trên mỗi pad là điều không thể trong thực tế, do đó chúng ta cần có sự hiệu chỉnh cần thiết cho những kết quả tính toán độ mất năng lượng một cách chính xác. Kết quả việc hiệu chỉnh độ mất năng lượng chi tiết được chỉ trong luận án. Bảng 3.1. Thông số hiệu chỉnh trong tính toán độ mất năng lượng. ∆Eexp = a × ∆ESRIM + b 22Mg 4He R 2 a 0.949 ± 0.003 0.781 ± 0.010 0.999 b 0.466 ± 0.001 0.066 ± 0.002 0.999 3.3. Phân tích số liệu 22Mg(α,α)22Mg 3.3.1. Thuật toán Số liệu được xử lí theo phương pháp từng sự kiện (event - by - event). Mỗi sự kiện chứa một lượng lớn thông số: vị trí chùm hạt trên PPAC, 48 tín hiệu độ mất năng lượng của chùm hạt tới trong vùng low-gain của detector, 8 tín hiệu độ mất năng lượng của hạt alpha trong vùng high- gain, 144 tín hiệu của các SSD (xem hình 3.5). Định dạng file chứa số liệu được chỉ ra trong bảng 3.2. Những sự kiện thoả mãn các điều kiện vật lý của phản ứng mà chúng ta quan tâm sẽ được chọn lọc. 18 Bảng 3.2. Định dạng của fille chứa các thông số trong một sự kiện. Event PPACa PPACb Low-gain Pad1 Pad48 High-gain Pad1 Pad8 Strip 1 X1Y1-X2Y2 X1Y1-X2Y2 P1 P48 P1 P8 Si 2 X1Y1-X2Y2 X1Y1-X2Y2 P1 P48 P1 P8 Si n X1Y1-X2Y2 X1Y1-X2Y2 P1 P48 P1 P8 Si Hình 3.5. Mô tả các thông số ghi nhận trong một sự kiện. Để chọn lọc các hạt 22Mg, chúng ta cần phải thực hiện đặt điều kiện cho nó, gọi là gate beam. Việc này được tiến hành dựa vào sự nhận biết các hạt theo phương pháp dE/dx bằng detector khí và thuật toán xác định đường bay của hạt 22Mg. Việc xác định hướng bay của 22Mg được tiến hành dựa trên độ mất năng lượng khác nhau giữa các anode trong một pad. Nếu hai anode có độ mất năng lượng như nhau thì đó chính là những hạt 22Mg bay thẳng chính giữa bia mà chúng ta quan tâm, như chỉ ra trên hình 3.6. Hình 3.6. Chọn lọc chùm hạt 22Mg dựa trên độ phổ biến của chùm hạt. Bước tiếp theo là chúng ta lọc các sự kiện chứa các hạt alpha và 22Mg và xác định được vị trí của hạt alpha trên các strip của SSD. Bước thứ ba là xác định vị trí phản ứng và góc tán xạ của alpha. Bước thứ tư xác định năng lượng của phản ứng. Khi đó, chúng ta thu được các thông số vật lý của phản ứng: năng lượng phản ứng Ereact, năng lượng alpha Eα, và góc tán xạ θα: Ereact. = Einc. - Ewd - Edl - Edz, (3.3) 19 Eα = Edl1 + EGM + Edl2 + ESSD, (3.4) θα = (PPACa & PPACb) + (R.Point) + (SSD strips). (3.5) Bước thứ năm là xác định hàm kích thích cộng hưởng của tiết diện phản ứng theo năng lượng ( )cmd Edσ Ω . Bước cuối cùng là khớp các giá trị thực nghiệm theo phương pháp R- matrix để xác định các đặc trưng lượng tử của các mức kích thích của 26Si. Những thông số cộng hưởng này tính toán suất phản ứng cho 22Mg(α,p)25Al. 3.3.2. Chương trình máy tính dùng để phân tích số liệu Để xử lí số liệu của thực nghiệm này, một số chương trình máy tính được tác giả luận án lập trình dựa trên các công thức và thuật toán đã trình bày trên. Các chương trình này chạy hệ điều hành Linux và sử dụng kết hợp với chương trình ROOT của trung tâm gia tốc châu Âu, CERN. 3.3.3. Tính toán động học Tính toán động học của phản ứng 22Mg(α,α)22Mg được tiến hành theo mô tả trên hình 3.7. Năng lượng phản ứng Ecm được xác định bởi năng lượng hạt bay vào Elab, 22Mg. Ta có: 4 4 22 4 , 26 He cm lab lab He Mg A E E E A A = = + (3.6) Năng lượng của hạt Elab của 22Mg được xác định bởi: 0 .labE E E= − ∆ (3.7) Hình 3.7. Mô tả hình học việc tính toán động học cho phản ứng. Góc tán xạ được tính toán theo các công thức hình học bên dưới, kết quả tính toán được trình bày trong phụ lục G của luận án. 2 2 tan , cos . OH OP OP OP OH θ θ= = + (3.8) 3.3.4. Độ sai số năng lượng Độ sai lệch năng lượng được đánh giá dựa trên năng lượng phân tán của chùm hạt vào cửa sổ buồng phản ứng, năng lượng phân tán trong bia khí, sai số năng lượng tại vị trí phản ứng, sai số năng lượng của alpha do sai số của góc tán xạ và do phân tán trong bia khí trước khi đến SSD. 20 3.3.5. Góc khối Góc khối của detector được xác định trong phạm vi các góc tán xạ của alpha trên SSD. Góc khối vi phân được xác định như sau: sin .lab lab labd d dθ θ ϕΩ = (3.9) Góc khối xác định trong một khoảng góc tán xạ được xác định bởi: sin .d d dθ θ ϕ∆Ω= Ω=∫ ∫ (3.10) 3.3.6. Số hạt tham gia phản ứng Thời gian tiến hành toàn bộ thực nghiệm là 17 ngày, trong đó mất 5 ngày để sản xuất chùm 22Mg và 1.5 ngày để đo độ mất năng lượng cũng như chuẩn detector bằng chùm hạt không bền và hạt nặng. Phản ứng 22Mg+α được đo trong 9.5 ngày với tổng số hạt vào là 2.36 × 109. 3.3.7. Tiết diện vi phân và trạng thái cộng hưởng Tiết diện vi phân được xác định trong phạm vi góc từ 0 - 5 và 5 -10 độ. Năng lượng phản ứng trong dải từ 0.5 MeV đến 3.0 MeV trong hệ khối tâm. Tiết diện phản ứng thực nghiệm xác định bởi: 0 . . . . . .A I Ad d I N σ τ ε ∆Ω= Ω ∆Ω (3.11) Tiết diện trong hệ khối tâm được xác định thông qua hệ phòng thí nghiệm : ( ) ( )1, , . 4coscm cm lab labcm lab lab d dE E d d σ σ θ θ θ = Ω Ω (3.12) Hình 3.8. Hàm kích thích của tiết diện phản ứng trong hệ khối tâm với góc tán xạ θlab = 0 - 5 độ. Hình 3.9. Hàm kích thích của tiết diện phản ứng trong hệ khối tâm với góc tán xạ θlab = 5 - 10 độ. Kết quả phân tích số liệu các mức cộng hưởng thu được của hạt nhân hợp phần 26Si được chỉ ra trong hình 3.8 và 3.9. Chúng ta thu được các 21 thông số ban đầu cần thiết cho việc làm khớp theo phương pháp R-matrix như trong bảng 3.3. Bảng 3.3. Năng lượng cộng hưởng và độ rộng xác định từ thực nghiệm Năng lượng Ecm(MeV) Độ rộng(keV) 10.360 ± 0.031 101.536 ± 3.116 10.680 ± 0.045 98.103 ± 2.021 10.829 ± 0.050 111.326 ± 0.811 11.242 ± 0.063 103.536 ± 0.927 11.481 ± 0.070 118.821 ± 3.709 11.803 ± 0.079 133.617 ± 1.073 3.4. Phân tích R-matrix cho phản ứng 22Mg(α,α)22Mg Việc phân tích làm khớp R-matrix được tiến hành bằng chương trình máy tính AZURE. Chương trình này đã được kiểm tra với các phản ứng hạt nhân thiên văn 12C+p, 14N+p và 16O+p có độ chính xác cao so với kết quả thực nghiệm. Theo phương pháp R-matrix, ta cần thông số cho hai vùng bên trong và bên ngoài hạt nhân compound. Bán kính kênh phản ứng được xác định bởi: ( )1/3 1/31 21.45R A A= + fm. (3.13) Những thông số cần thiết cho vùng bên trong là các thông số lượng tử. Theo nguyên lý bảo toàn spin và chẵn lẻ, đối với 22Mg + α
26Si, ta có: 1 2J j j l s= + = +     ,l s J l s⇒ − ≤ ≤ + (3.14) ( )1 2 1 . lpi pi pi= − (3.15) Do 22Mg, 4He và26Si là hạt nhân chẵn-chẵn, chúng có spin s = 0 nên: J = l và pi = (-1)l . (3.16) Theo cơ học lượng tử ta có l R p = ×     , do đó, moment quỹ đạo có giá trị cực đại thoả mãn: 2 . R E l µ ≤
(3.17) Với năng lượng cực đại của phản ứng là E = 3.0 MeV thì l ≤ 4, ta thu được các thông số lượng tử khả dĩ của các cộng hưởng như trong bảng 3.4. Đối với vùng ngoài, thông số kênh vào của phản ứng trong bảng 3.5. Sau khi làm khớp bình phương tối thiểu theo phương pháp R-matrix chúng ta thu được kết quả như trình bày trong bảng 3.6 và hình 3.10. Trong thực nghiệm này đã xác định được sáu mức kích thích, trong đó mức đầu và mức 22 cuối có hai giá trị Jpi đều thoả mãn thông số khớp R-matrix. Do đó, chúng ta cần có những nghiên cứu khác để xác định một trong các giá trị này. Bảng 3.4. Thông số cộng hưởng dùng cho AZURE. Bảng 3.5. Thông số ban đầu dùng cho AZURE. Các mức l Jpi 10.36; 10.68; 10.83; 11.24; 11.48; 11.80 0 0+ Projectile 22Mg ( 0J pi += ) 1 1- Target 4He( 0J pi += ) 2 2+ Compound nucleus 26Si( 0J pi += ) 3 3- Separation energy 9.221 (MeV) 4 4+ Channel radius 6.365 (fm) Hình 3.10. Kết quả khớp R-matrix các số liệu thực nghiệm để xác định trạng thái lượng tửcủa các mức kích thích cộng hưởng. Bảng 3.6. Các mức kích thích của 26Si thu được trong nghiên cứu này được so sánh với hai công trình khác với độ tin cậy cao. Mức Nghiên cứu này 28Si(p,t) 26P Er (MeV) Γ (MeV) J pi Er (MeV) Jpi Er (MeV) J pi 1 10.325 ± 0.071 --- 0.218 ± 0.011 --- (2+, 1-) --- --- 10.434 -?- -?- 10.299 10.405 -?- -?- 2 10.678 ± 0.016 0.194 ± 0.006 0+ 10.656 -?- 10.688 -?- 3 10.831 ± 0.113 --- 0.186 ± 0.013 --- 1- --- --- 11.014 -?- --- 10.827 --- -?- --- 4 11.245 ± 0.028 0.208± 0.027 4+ --- --- --- --- 5 11.493 ± 0.216 0.292 ± 0.010 3- --- --- --- --- 6 11.807 ± 0.117 0.156 ± 0.032 (0+, 2+) --- --- --- --- 3.5. Các mức năng lượng cộng hưởng trên ngưỡng alpha của 26Si Theo quy luật của cấu trúc nhóm của hạt nhân, trạng thái kích thích cộng hưởng thường xảy ra gần với năng lượng ngưỡng hình thành nhóm hoặc năng lượng cộng hưởng của nhóm đó, nhất là đối với những hạt nhân 23 có số khối gần gấp bốn lần alpha như 4n, 4n±1, 4n±2, như trường hợp hạt 22Ne, 22Mg. Tương tự với 26Si, chúng ta có thể có cấu trúc nhóm trong nó. Trạng thái 10.325 MeV và 10.831 MeV Vì vùng năng lượng (10.3 - 10.8 MeV) của 26Si cao hơn ngưỡng 3 alpha 3 thrE α (7.2657 MeV) và thấp hơn ngưỡng 4 alpha 4thrE α (14.4382 MeV). Do đó, có thể dự đoán cấu trúc hạt nhân 26Si dạng (p+3α+3α+p). Điều này được dự đoán dựa trên cơ sở là các mức năng lượng của 26Si gần sát với các mức cộng hưởng của 12C và cấu trúc nhóm 3α của 12C. Trạng thái 10.678 MeV, 11.245MeV và 11.493 MeV Các mức này có năng lượng nằm xa các mức cộng hưởng của 12C, do đó, rất khó để có cấu trúc nhóm p+3α+3α+p. Tuy nhiên, các mức này lại rất gần với các mức kích thích cộng hưởng của 24Mg, do đó cấu trúc của 26Si có thể tồn tại dưới dạng lõi là 24Mg liên kết với hai proton “hoá trị”, p+24Mg+p. Các mức cộng hưởng 24Mg, phụ lục E trong luận án, có spin và chẵn-lẻ tương tự như của 26Si nếu spin của hai proton đi kèm bị triệt tiêu. Như đã nói trên, việc xác định trạng thái lượng tử của mức đầu và cuối có hai kết quả cần được kiểm chứng bởi những số liệu khác. Tuy nhiên, cho đến nay chưa có số liệu nào được công bố. Bên cạnh đó, khoảng cách mức và độ rộng mức thu được từ thực nghiệm này rất gần nhau. Do đó, việc tiến hành thực nghiệm xác định các mức năng lượng trên ngưỡng alpha của 26Si vẫn là một việc làm rất cần thiết. Trong hướng nghiên cứu tiếp theo của luận án, kết quả từ việc phân tích số liệu của 22Mg(α,p)25Al sẽ bổ sung thêm những thông tin cần thiết về 26Si trên ngưỡng alpha. 3.6. Suất của phản ứng 22Mg(α,p)25Al Suất phản ứng của 22Mg(α,p)25Al được tính toán bằng cách sử dụng số liệu các mức cộng hưởng của hạt nhân 26Si và theo các công thức trình bày trong chương 1. Trong đó, hàm lực cộng hưởng được xác định như trình bày trong bảng 3.7. Kết quả tính toán suất phản ứng với các thông số hàm lực nêu trên được trình bày trong bảng 3.8. Hình 3.11 trình bày suất phản ứng của 22Mg(α,p)25Al ứng với các mức cộng hưởng của 26Si. So sánh kết quả của nghiên cứu này với nghiên cứu về quá trình bắt proton, hình 3.12, và phân rã beta, kết quả trong bảng 3.9, cho thấy rằng sự cạnh tranh của 24 phản ứng bắt alpha nhỏ hơn so với bắt proton. Tốc độ phân rã beta của 22Mg nhanh hơn so với quá trình bắt alpha nhưng chậm hơn bắt proton. Do đó, tại vùng năng lượng khảo sát trên siêu tân tinh và X-ray Burst, quá trình bắt proton chiếm ưu thế và do đó nó bỏ qua cơ chế phân rã beta của 22Mg. Điều này ảnh hưởng đến việc quan sát tia gamma 1.275 MeV và sự bất thường của tỉ số 22Ne/20Ne trong vấn đề Ne-E hiện nay. Bảng 3.7. Hàm lực công hưởng của các mức kích thích của 26Si. Level Er (MeV) ER = Er - Q (MeV) J pi Γ (MeV) ωγ (MeV) 1 10.325 1.209 2+ 0.218 0.2725 2 10.678 1.512 0+ 0.194 0.0485 3 10.831 1.675 1- 0.186 0.1395 4 11.245 2.057 4+ 0.208 0.4680 5 11.493 2.306 3- 0.292 0.5180 6 11.807 2.640 0+ 0.156 0.0390 Bảng 3.8. Suất phản ứng 22Mg(α,p)25Al tính theo các mức cộng hưởng. T9 1.209 MeV 1.512 MeV 1.675 MeV 2.057 MeV 2.306 MeV 2.640 MeV Tổng 0.5 9.6E-18 1.5E-21 9.8E-23 4.7E-26 1.6E-28 5.1E-33 9.6E-18 1.0 4.2E-12 2.2E-14 9.6E-15 3.8E-16 2.3E-17 3.7E-20 4.2E-12 1.5 2.4E-10 4.2E-12 3.4E-12 6.0E-13 9.7E-14 5.5E-16 2.5E-10 2.0 1.6E-09 5.1E-11 5.7E-11 2.0E-11 5.4E-12 5.9E-14 1.8E-09 2.5 4.8E-09 2.1E-10 2.8E-10 1.6E-10 5.6E-11 9.0E-13 5.5E-09 3.0 9.4E-09 5.1E-10 7.9E-10 6.0E-10 2.5E-10 5.3E-12 1.1E-08 3.5 1.4E-08 9.4E-10 1.5E-09 1.5E-09 7.2E-10 1.8E-11 1.9E-08 4.0 1.9E-08 1.4E-09 2.6E-09 2.8E-09 1.5E-09 4.4E-11 2.8E-08 Hình 3.11. Suất phản ứng của 22Mg(α,p)25Al ứng với các mức cộng hưởng Hình 3.12. Suất phản ứng của 22Mg(p,γ)23Al. 25 Bảng 3.9. Tốc độ phản ứng 22Mg(p,γ)23Al, 22Mg(α,p)25Al và phân rã beta. T9 (GK) Rαp (reaction/s) Rpγ (reaction/s) Rβ+ (decay/s) 1.0 2.99E-06 7.80E+03 0.178881 1.5 1.79E-04 2.47E+04 0.178881 2.0 1.26E-03 4.92E+04 0.178881 2.5 3.90E-03 7.80E+04 0.178881 10 7.50E-02 --- 0.178881 Số liệu suất phản ứng ghi nhận trong nghiên cứu này được so sánh với kết quả tính toán theo lý thuyết thống kê của Hauser-Feshbach với code NON-SMOKERWEB. Kết quả được chỉ ra trong hình 3.13. Hình 3.13. Suất phản ứng trong nghiên cứu này tính toán từ thực nghiệm và từ lý thuyết thống kê thu đượcbởi NON-SMOKERWEB Hình 3.14. S-factor như một hàm theo năng lượng. Thừa số S-factor là một thông số quan trọng trong việc xác định suất phản ứng. Trong nghiên cứu này, thừa số đã được rút ra từ suất phản ứng và các mức cộng hưởng của 26Si như sau: ( ) 1/2 3/2 1/2 8 1 . exp .ii A i E bS E N v kT kT E σ piµ − −      = +           (3.18) T là nhiệt độ trong môi trường sao ứng với các mức năng lượng E. Kết quả S-factor của các mức cộng hưởng đã xác định được chỉ ra trong hình 3.14. và bảng 3.10. Để có thể tính toán ngoại suy cho các mức năng lượng khác, các giá trị trên được làm khớp với hàm đa thức như sau: S(E) = -0.008E2 + 22.915E + 292.143. (3.19) 26 Bảng 3.10. S-factor S(E) tại các mức năng lượng cộng hưởng. E (MeV) Jpi S-factor (keV.barns) 1.209 2+ 1.506×104 ± 626.848 1.512 0+ 1.640×104 ± 542.486 1.675 1- 1.558×104 ± 478.484 2.057 4+ 1.175×104 ± 287.019 2.306 3- 7.008×103 ± 264.092 2.64 0+ 2.746×103 ± 128.035 Từ công thức (3.19) chúng ta có thể tính được suất phản ứng tại các mức năng lượng trong đỉnh Gamow như sau: 1/32 2 1 2 9 1/3 4.2487 9 1 2 02 9 7.83 10 ( ) e . (3.20) Z Z A T A Z ZN v S E AT σ    −       = ⋅     Kết luận và kiến nghị Trong nghiên cứu này chúng tôi đã xác định được 06 mức kích thích trong vùng năng lượng trên ngưỡng alpha của hạt nhân giàu proton 26Si. Trong đó, có 03 mức trùng khớp với những nghiên cứu trước đây và 03 mức hoàn toàn mới lần đầu tiên được ghi nhận. Các trạng thái spin và chẵn- lẻ của 06 mức lần đầu tiên được xác định bởi nghiên cứu này. Cấu trúc nhóm của các mức này đã được đánh giá dựa trên những số liệu ghi nhận được từ hạt nhân 12C và 24Mg. Các mức năng lượng này có thể được bổ sung vào sơ đồ mức của 26Si và góp phần trong việc xác định trạng thái lượng tử của hạt nhân gương 26Mg hiện tại còn rất thiếu. Kết quả này cũng góp phần bổ sung số liệu hạt nhân hiện nay. Tuy nhiên, việc thu thập các số liệu của 26Si bằng các nghiên cứu khác theo cả phương pháp lý thuyết và thực nghiệm là một việc làm rất cần thiết vì đây là những kết quả lần đầu tiên được ghi nhận. Về vấn đề thiên văn học, suất phản ứng của 22Mg(α,p)25Al đã được tính toán và so sánh với các nhánh phân rã beta và bắt proton của 22Mg. Kết quả chứng tỏ rằng, ở vùng nhiệt độ T = 0.5 - 3 GK, quá trình bắt proton diễn ra nhanh hơn bắt alpha. Do đó, khi nhiệt độ tăng cao như trong vùng môi trường của siêu tân tinh và X-ray burst, khả năng cạnh trạnh của phân rã beta không đáng kể. Đây có thể là một trong những lý do mà tia gamma đặc 27 trưng 1.275 MeV bị bỏ qua trong quá trình tổng hợp hạt nhân trong các sao có nhiệt độ cao. Điều này tác động lên độ phổ biến của 22Na và tỉ số 22Ne/20Ne. Thừa số S-factor dùng cho tính toán phản ứng 22Mg(α,p)25Al theo cơ chế trực tiếp cũng được xác định trong nghiên cứu này. Về mặt thực nghiệm, nghiên cứu này đã tạo được chùm 22Mg thoả mãn điều kiện về cường độ và độ sạch để có thể đo trực tiếp tán xạ lên alpha. Một loại detector khí mới được thiết kế và phát triển để dùng trong các thực nghiệm vật lý hạt nhân, nghiên cứu cấu trúc và thiên văn theo cơ chế động học ngược một cách hiệu quả. Hướng phát triển tiếp theo của luận án là xác định tiết diện của 22Mg(α,p)25Al để xác định suất phản ứng từ tiết diện, thay vì tính suất phản ứng bằng các mức cộng hưởng. Đây là một việc làm cần thiết để đánh giá việc khảo sát suất phản ứng thông qua hai phương pháp, dùng các mức cộng hưởng và dùng tiết diện phản ứng. Trong đó, việc phát triển chương trình máy tính ANAROOT để xử lí số liệu thực nghiệm của 22Mg(α,p)25Al cũng sẽ được tiến hành. Ngoài ra, việc áp dụng các mẫu lý thuyết như mẫu quang học, mẫu lớp để đánh giá các mức đã ghi nhận cùng với những kết quả thực nghiệm khác cũng sẽ được thực hiện. Các công bố khoa học Các công trình công bố quốc tế và nước ngoài [1]. N. N. Duy, L. H. Khiem, S. Kubono, D. Kahl et al. (2013), “Low- energy radioactive ion beam production of 22Mg”, Nuclear Instruments and Methods in Physics Research A 723, pp. 99–101. [2]. N. N. Duy, L. H. Khiem, S. Kubono, D. Kahl et al. (2013), “Exotic beam for direct measurement of the 22Mg+α”, Journal of Physical Science and Application (USA), Vol.3, No.2, pp.82 - 85. [3]. N. N. Duy, L. H. Khiem, S. Kubono, D. Kahl et al. (2013), “Role of the 22Mg(α,p)25Al reaction and its first direct measurement”, CNS Annual Report 2011 - the University of Tokyo, pp. 9 - 10. [4]. N. N. Duy, L. H. Khiem, S. Kubono, D. Kahl et al. (2013), “RI beam production in the 22Mg+α experiment at CRIB”, CNS Annual Report 2011 - the University of Tokyo, pp. 7 - 8. 28 [5]. N. N. Duy, L. H. Khiem, S. Kubono, D. Kahl et al. (2012), “Active target measurement of the 22Mg+alpha system in inverse kinematics”, RIKEN Accelerator Process Report 45, p.20. [6]. H. Yamaguchi, S. Kubono, D. Kahl, L. H. Khiem, N. N. Duy et al. (2013), “α-resonance structure in 11C studied via resonant scattering of 7Be+α and with the 7Be(α,p) reaction”, Physics Review C 87, pp. 1 - 12. [7]. H. Yamaguchi, S. Kubono, D. Kahl, L. H. Khiem, N. N. Duy et al. (2012), “Excitation function of 7Be+α elastic and inelastic scattering cross sections”, CNS Annual Report 2010 - the University of Tokyo, pp. 1 - 2. [8]. H. Yamaguchi, S. Kubono, D. Kahl, L. H. Khiem, N. N. Duy et al. (2011), “Measurement of alpha scattering on 7Be”, CNS Annual Report 2009 - the University of Tokyo, pp. 1 - 2. [9]. D. Kahl, A.A.Chen, S. Kubono, H. Yamaguchi, S. Cherubini, N. N. Duy et al. (2012), “Measurement of the 30S+a system for type I X-ray bursts”, XII International Symposium on Nuclei in the Cosmos, pp. 1 - 5. Các công trình công bố trong nước [10]. N.N.Duy and L.H. Khiem (2011), “Feasibility study of direct measurement of stellar reaction 22Mg(α, p)25Al”, Communications in Physics, Vol. 21, No.2, pp. 169 - 178. [11]. N.N.Duy, L.H. Khiem and V. H. Tan (2012), “Investigation of gas gain of GEM-foil used in low energy radioactive beam experiment”, Communications in Physics, Vol. 22, No.3, pp. 283 - 287. [12]. N.N.Duy, L.H. Khiem and V. H. Tan (2011), “Gas gain measurement of GEM-foil in Argon-Carbon dioxide mixture”, Tuyển tập báo cáo hội nghị khoa học và công nghệ hạt nhân toàn quốc lần thứ IX, Ninh Thuận 18- 19/8/2011, pp.209 - 211. [13]. N.N.Duy and L.H. Khiem (2011), “Design of experiment for (α,p) reaction induced by 22Mg radioactive ion beam”, Proceedings of the topical conference on Nuclear Physics, High energy Physics and Astrophysics, Hanoi 9-11/11/2010, pp.90 - 95.