Bài viết Thiết kế nguyên tắc hệ thống chùm positron chậm bằng chương trình mô phỏng Simion trình bày các kết quả nghiên cứu và áp dụng chương trình mô phỏng quỹ đạo hạt mang điện Simion để xây dựng các mô hình thiết kế khả thi cho hệ thống. Một số tính toán mô phỏng thử nghiệm quỹ đạo chuyển động của chùm hạt positron đã được tiến hành nhằm mục đích so sánh giữa các mô hình.
THIẾT KẾ NGUYÊN TẮC HỆ THỐNG CHÙM POSITRON CHẬM BẰNG CHƢƠNG TRÌNH MƠ PHỎNG SIMION CAO THANH LONG1, BÙI XN HUY2 Trung tâm Hạt nhân Thành phố Hồ Chí Minh, 217 Nguyễn Trãi, Phường Nguyễn Cư Trinh, Quận 1, Thành phố Hồ Chí Minh Đại học Sư Phạm Thành phố Hồ Chí Minh, 280 An Dương Vương, Phường 4, Quận 5, Thành phố Hồ Chí Minh ctlong26051993@gmail.com, xuanhuy_vl@yahoo.com Tóm tắt: Hệ thống chùm positron chậm thiết bị quan trọng nghiên cứu vật lý kỹ thuật positron, đặc biệt ứng dụng nghiên cứu vật liệu Báo cáo chúng tơi trình bày kết nghiên cứu áp dụng chương trình mơ quỹ đạo hạt mang điện Simion để xây dựng mơ hình thiết kế khả thi cho hệ thống Một số tính tốn mơ thử nghiệm quỹ đạo chuyển động chùm hạt positron tiến hành nhằm mục đích so sánh mơ hình Từ đó, mơ hình thiết kế ngun tắc đề xuất để làm sở cho việc xây dựng hệ thống chùm positron chậm thực tế thực tương lai Việt Nam Từ khố: Hệ thống chùm positron chậm, chương trình mơ phỏng, Simion, thiết kế nguyên tắc I MỞ ĐẦU Phương pháp phổ kế huỷ cặp positron (Positron Annihilation Spectroscopy – PAS) đóng vai trị quan trọng việc nghiên cứu dạng sai hỏng vật liệu, cấu trúc nano trạng thái kim loại, tạp chất bên vật liệu [1] Do nguồn positron đồng vị truyền thống sử dụng phương pháp PAS có đặc trưng phát positron với phổ phân bố lượng liên tục, phương pháp không cung cấp thông tin chi tiết sai hỏng bề mặt vật liệu nằm độ sâu khác [2] Để khắc phục nhược điểm trên, hệ thống chùm positron chậm chế tạo đưa vào hoạt động số Viện nghiên cứu, phịng thí nghiệm trường Đại học lớn nước phát triển nhằm mục đích tạo chùm positron đơn có lượng thấp để thay cho nguồn phát đồng vị truyền thống Các hệ thống có thiết kế khác tuân theo nguyên lý thiết kế chung Một tỉ lệ positron có lượng cao phát từ nguồn phát làm chậm để trở thành positron chậm có lượng thấp (vài eV) tương đối đơn Các positron chậm sau gia tốc nhẹ nhờ phận tiền gia tốc dẫn qua phận lọc lượng để tách hạt positron có lượng cao khỏi chùm tia Chùm positron chậm tương đối đơn sau lọc khỏi chùm tia ban đầu tiếp tục dẫn qua phận gia tốc để gia tốc đến mức lượng cần thiết dẫn đến buồng chứa mẫu để bắn phá bề mặt vật liệu cần khảo sát Trung tâm Hạt nhân Thành phố Hồ Chí Minh (Trung tâm) bước đầu áp dụng phương pháp PAS bao gồm phương pháp đo thời gian sống positron giãn nở Doppler để nghiên cứu số tính chất vật liệu kim loại, ống nano carbon, zeolite có số kết định [3-6] Tuy nhiên nguồn phát positron chủ yếu sử dụng nguồn 22Na có phổ lượng rộng, gây hạn chế lớn nghiên cứu tính chất bề mặt vật liệu Do đó, nhu cầu xây dựng hệ thống phát chùm positron chậm tương lai Trung tâm để phục vụ cho nghiên cứu ứng dụng chuyên sâu cần thiết Để đảm bảo tính khả thi cho việc xây dựng hệ thống, công việc tiên cần làm giai đoạn đầu tìm hiểu sử dụng chương trình mơ quỹ đạo hạt để phục vụ cho việc thiết kế nguyên tắc cho hệ thống Trong số chương trình mơ tiêu biểu, chương trình Simion với nhiều tính tốt sử dụng rộng rãi với ứng dụng mơ trường điện từ tính tốn quỹ đạo hạt mang điện điện từ trường Chương trình Simion sử dụng cho việc nghiên cứu xây dựng hệ thống chùm positron chậm số phịng thí nghiệm giới Viện nghiên cứu Rossendorf (Đức), Phịng Thí nghiệm quốc gia Lawrence Livermore (Hoa Kỳ), trường Đại học Bath (Anh) nước khác Israel, Trung Quốc [7-10] Bên cạnh đó, số tính tốn mơ thử nghiệm dựa thiết kế tham khảo hệ thống SPONSOR thuộc Viện nghiên cứu Rossendorf chương trình Simion nhóm nghiên cứu Trung tâm thực [11] phù hợp kết mô với kết tính tốn, đo đạc thực nghiệm hệ thống SPONSOR chứng tỏ chương trình Simion cơng cụ tính tốn thích hợp, sử dụng cho việc thiết kế nguyên tắc hệ thống chùm positron chậm Báo cáo trình bày kết nghiên cứu áp dụng chương trình mơ Simion để xây dựng mơ hình thiết kế khả thi hệ thống chùm positron chậm dựa nguyên lý thiết kế số hệ thống tiêu biểu đưa vào hoạt động giới từ đưa mơ hình thiết kế ngun tắc khả thi làm sở cho việc thiết kế chi tiết kỹ thuật, chế tạo, xây dựng cho hệ thống chùm positron chậm thực tương lai Trung tâm II NỘI DUNG II Tổng quan chƣơng trình mơ Simion Simion chương trình mơ chun dụng sử dụng để mơ hình hố, mơ trường tĩnh điện từ trường đồng thời tính tốn quỹ đạo chuyển động hạt mang điện điện từ trường [12] Một cách tổng quan, chương trình Simion giải hai tốn phương pháp số bao gồm tính tốn trường tĩnh điện, từ trường toán xác định quỹ đạo chuyển động hạt mang điện điện từ trường tính tốn Phương pháp sai phân hữu hạn (Finite differential method) tối ưu hoá kỹ thuật lặp over-relaxation chương trình sử dụng để giải tốn tìm nghiệm phương trình đạo hàm riêng phần, cụ thể phương trình Laplace (hoặc phương trình Poisson) nhằm xác định giá trị điểm khơng gian trường Đối với tốn xác định quỹ đạo chuyển động hạt mang điện, chương trình sử dụng phương pháp tích phân số RungeKutta bậc bốn để dự đốn vận tốc vị trí hạt mang điện dựa vào thông số phát hạt ban đầu xác định người dùng dựa trường tính tốn, qua xác định quỹ đạo hạt trường tĩnh điện từ trường xác định trước Một tính hữu ích chương trình Simion việc sử dụng ngơn ngữ lập trình bên chương trình, trình biên dịch ngôn ngữ Lua (Lua 5.1) nhúng trực tiếp bên chương trình Simion từ phiên 8.0 trở Việc sử dụng ngôn ngữ Lua giúp mở rộng khả mơ chương trình cách linh hoạt hiệu II Phƣơng pháp nghiên cứu Dựa thơng tin tiếp cận từ thiết kế hệ thống chùm positron chậm mà nhóm nghiên cứu tham khảo, số mơ hình thiết kế hệ thống xây dựng chương trình mơ Simion Vì hạn chế khả ứng dụng chương trình, mơ hình thiết kế dùng để mô nguyên lý hoạt động hệ thống chùm positron chậm cách đơn giản chưa phải thiết kế chi tiết kỹ thuật Các mơ hình mơ chia làm hai loại sau: - Mơ hình mơ thiết kế theo dạng thẳng, sử dụng phương pháp lọc chùm positron chậm lọc điện từ trường ExB (hình 1) - Mơ hình mơ thiết kế theo dạng cong, sử dụng phương pháp lọc chùm positron chậm đoạn ống cong kết hợp với từ trường chia làm hai mơ hình riêng lẻ bao gồm mơ hình sử dụng đoạn ống cong với góc cong 500 dựa thiết kế hệ thống SPONSOR [13, 14] mơ hình dùng đoạn ống cong với góc cong 900 dựa thiết kế hệ thống trường Đại học Halle [15] (hình 2) Hình Mơ hình hệ thống thiết kế theo dạng thẳng Hình Mơ hình hệ thống thiết kế theo dạng cong với góc cong 500 (trái) góc cong 900 (phải) Sau xây dựng thành cơng mơ hình mơ trên, nhóm nghiên cứu tiến hành số tính tốn thử nghiệm mơ quỹ đạo chuyển động chùm hạt positron nhằm mục đích so sánh mơ hình dựa việc khảo sát chất lượng chùm hạt positron thu bia mẫu cho ba mơ hình thiết kế Chùm hạt positron mô với thông số đầu vào giống cho mơ hình cho trường hợp khảo sát Các trường hợp khảo sát thực bao gồm: - Khảo sát tính tốn quỹ đạo bay chùm hạt positron đơn Khảo sát thực nhằm đánh giá độ hội tụ độ đơn chùm hạt positron thu bia mẫu mơ hình - Khảo sát tính tốn quỹ đạo bay chùm hạt positron đơn trường hợp cuộn solenoid bị lệch khỏi vị trí tối ưu ban đầu Khảo sát thực với mục đích đánh giá độ nhạy chất lượng chùm hạt positron thu bia mẫu phát sinh tình lắp đặt sai vị trí cuộn dây solenoid, dẫn đến sai lệch từ trường tối ưu dọc theo trục mơ hình Nhóm nghiên cứu đánh giá, so sánh kết tính tốn mô đưa lựa chọn thiết kế xem tối ưu số mơ hình Mơ hình tối ưu đề xuất làm mơ hình thiết kế nguyên tắc cho hệ thống chùm positron chậm xây dựng tương lai Bản thơng số thiết kế cuối mơ hình bao gồm thơng số dạng hình học hệ thống, thông số thiết kế cuộn dây điện tạo từ trường thông số điện mơ hình tối ưu lựa chọn II Kết bàn luận Trường hợp 1: Khảo sát mô với chùm hạt positron đơn Chương trình Simion sử dụng để mô quỹ đạo bay chùm hạt bao gồm 1000 hạt positron đơn năng, có động ban đầu eV phát đẳng hướng từ cửa sổ nguồn giả định cho mơ hình thiết kế Giá trị cao cung cấp cho phận tiền gia tốc phận gia tốc cho ba mơ hình cách tương ứng 27 V 30 kV Các kết thống kê chùm hạt positron thu bia mẫu cho mơ hình trình bày bảng Các kết biểu diễn phân bố chùm hạt phân bố lượng chùm hạt bia mẫu cho mơ hình mơ tả hình hình Bảng Kết thống kê bia mẫu trường hợp khảo sát với chùm hạt positron đơn Tổng số positron phát từ nguồn Tổng số positron đến bia mẫu Tỉ lệ positron đến bia mẫu Bán kính tiết diện chùm positron bia mẫu Mơ hình dạng thẳng Mơ hình dạng cong với góc cong 500 Mơ hình dạng cong với góc cong 900 1000 1000 1000 789 807 795 78,9% 80,7% 79,5% 2,75 mm 2,39 mm 2,79 mm Hình Phân bố chùm hạt bề mặt bia mẫu cho mơ hình thiết kế theo dạng thẳng (trái), thiết kế theo dạng cong với góc cong 500 (giữa) góc cong 900 (phải) cho trường hợp khảo sát chùm hạt positron đơn Hình Phân bố lượng chùm hạt bề mặt bia mẫu cho mơ hình thiết kế theo dạng thẳng (trái), thiết kế theo dạng cong với góc cong 500 (giữa) góc cong 900 (phải) cho trường hợp khảo sát chùm hạt positron đơn Từ kết khảo sát cho trường hợp 1, thấy chùm hạt positron thu trường hợp sử dụng mơ hình với góc cong 500 có bán kính tiết diện chùm hạt phân bố bề mặt bia mẫu nhỏ nhất, có độ hội tụ tốt so với hai mơ hình cịn lại Ngồi ra, phổ phân bố lượng chùm hạt positron thu cho ba mơ hình thiết kế có đỉnh phổ nằm xung quanh mức lượng 30030 eV Kết phù hợp với đỉnh lượng dự đoán thu khảo sát chùm hạt positron đơn eV tiền gia tốc với điện áp 27 V gia tốc với cao 30 kV Bên cạnh đó, kết đánh giá phổ phân bố lượng cho thấy độ rộng bán cực đại phổ mơ hình thiết kế với góc cong 500 (FWHM = 1,34 eV) nhỏ so với kết tính tốn cho mơ hình thiết kế theo dạng thẳng (FWHM = 5,50 eV) mô hình thiết kế với góc cong 900 (FWHM = 5,71 eV) Kết cho thấy chùm hạt positron thu bia mẫu mơ hình thiết kế với góc cong 500 tương đối đơn so với hai mơ hình cịn lại Trường hợp 2: Khảo sát mô với chùm hạt positron đơn trường hợp cuộn solenoid bị lệch khỏi vị trí tối ưu ban đầu Quỹ đạo bay chùm hạt đơn mô với thông số đầu giống trường hợp khảo sát trường hợp cuộn solenoid bao quanh phận gia tốc mơ hình thiết kế giả định bị lệch khỏi vị trí tối ưu ban đầu với độ lệch cm cm Các kết so sánh phân bố chùm hạt positron thu bia mẫu cho mơ hình thiết kế mơ tả tương ứng hình 5, Hình Phân bố chùm hạt bề mặt bia mẫu cho mơ hình thiết kế theo dạng thẳng không làm lệch cuộn solenoid (trái) làm lệch cuộn solenoid cm (giữa) cm (phải) Hình Phân bố chùm hạt bề mặt bia mẫu cho mơ hình thiết kế theo dạng cong 500 không làm lệch cuộn solenoid (trái) làm lệch cuộn solenoid cm (giữa) cm (phải) Hình Phân bố chùm hạt bề mặt bia mẫu cho mơ hình thiết kế theo dạng cong 900 không làm lệch cuộn solenoid (trái) làm lệch cuộn solenoid cm (giữa) cm (phải) Từ kết khảo sát cho trường hợp 2, thấy với mơ hình thiết kế theo dạng thẳng sử dụng lọc ExB, độ lệch cuộn solenoid bao quanh phận gia tốc ảnh hưởng nhiều đến tỉ lệ hạt positron đến bia mẫu Với trường hợp cuộn solenoid bị lệch cm, sai lệch từ trường làm cho positron chậm sau khỏi lọc ExB gần bị chặn hoàn toàn khối collimator làm cho số hạt đến bia mẫu Trong trường hợp cuộn solenoid bị lệch cm, chí khơng có hạt positron đến bia mẫu Đối với hai mơ hình thiết kế theo dạng cong, kết so sánh cho thấy chùm hạt positron thu bia mẫu với mơ hình thiết kế theo dạng cong với góc cong 900 bị sai lệch nhiều so với việc sử dụng mơ hình với góc cong 500, chí chùm hạt positron thu cịn bị thay đổi hình dạng trường hợp cuộn solenoid giả định bị lệch cm mô tả hình Các kết so sánh cho thấy chất lượng chùm hạt positron thu với mơ hình thiết kế theo dạng cong với góc cong 500 nhạy với sai số gây lệch vị trí cuộn solenoid so với việc sử dụng hai mơ hình cịn lại III KẾT LUẬN Thơng qua việc thực nghiên cứu này, nhóm nghiên cứu nắm nguyên lý chung thành phần hệ thống chùm positron chậm; đồng thời sử dụng chương trình mơ Simion để xây dựng thành cơng mơ hình thiết kế khả thi cho hệ thống chùm positron chậm Từ việc xem xét tính khả thi mơ hình thiết kế đánh giá kết tính tốn mô thử nghiệm dẫn đến kết luận mô hình thiết kế theo dạng cong với góc cong 500 tối ưu so với hai mơ hình cịn lại Đây mơ hình đề xuất làm thiết kế nguyên tắc làm sở cho việc thiết kế chi tiết kỹ thuật thành phần thiết yếu hệ thống chùm positron chậm xây dựng tương lai Với mơ hình thiết kế ngun tắc lựa chọn, thông số mô hình thay đổi hồn thiện cách linh động qua tối ưu hóa thông số kỹ thuật cách khả thi hợp lý trước tiến tới việc thiết kế, xây dựng hệ thống thực tế TÀI LIỆU THAM KHẢO [1] P K Pujari, K Sudarshan and D Dutta (Ed.), “11th International Workshop on Positron and Positronium Chemistry (PPC-11)”, Journal of Physics: Conference Series, Volume 618, Conference 1, 2015 [2] P.G Coleman (Ed.), “Positron Beams and their applications”, World Scientific, Singapore, 2000 [3] L A Tuyen, Zs Kajcsos, K Lázár, T D Tap, D D Khiem, P T Phuc: “Positron annihilation characteristics in multi-wall carbon nanotubes with different average diameters”, Journal of Physics: Conference Series (IOP), 2013 [4] A T Luu, Zs Kajcsos, N D Thanh, T Q Dung, M V Nhon, K Lázár, K Havancsák, G Huhn, Z E Horváth, T D Tap, L T Son, and P T Phuc “Multi-Wall carbon nanotubes investigated by positron annihilation techniques and microscopies for further production handling” Phys Status Solidi C 6, No 11, 2578–2581, 2009 [5] Nguyen Duc Thanh, Tran Quoc Dung, Luu Anh Tuyen and Khuong Thanh Tuan, “Semiempirical formula for large pore- size estimation from o-Ps annihilation lifetime” Int J Nuclear Energy Science and Technology, Vol 4, No (2008) [6] T Q Dung, K Lázár, K Havancsák, Zs Kajcsos “o-Ps lifetimes in iron containing microand mesoporous media”, Materials Science Forum, Vol 733, pp 197-202, 2013 [7] F A Selim, A.W Hunt, J.A Golovchenko, R H Howell, R Haakenaasen, K.G Lynn, “Improved source and transport of monoenergetic MeV positrons”, Nuclear Instruments and Methods in Physics Research B 171 (2000),182-188, 2000 [8] S May-Tal Beck, D Cohen, E Cohen, A Kelleher, O Hen, J Dumas, E Piasetzky, N Pilip, G Ron, I Sabo-Napadensky, R Weiss-Babai, “Design of the Slow POsitron faciliTy (SPOT) in Israel”, 13th International Workshop on Slow Positron Beam Techniques and Applications, 2014 [9] C K Cheung, P S Naik, C D Beling, S Fung, H M Weng, “Performance of a slow positron beam using a hybrid lens design”, Department of Physics, University of Hong Kong, Pokfulam Road, Hong Kong, PR China, 2006 [10] Xu Hong-Xia, Liu Jian-Dang, Gao Chuan-Bo, Weng Hui-Min, Ye Bang-Jiao, “SIMION simulation of a slow pulsed positron beam”, Department of Modern Physics, University of Science and Technology of China, Hefei 230026, China, 2012 [11] Cao Thanh Long, Nguyen Trung Hieu, Tran Quoc Dung, Huynh Dong Phuong, “Some initial results of simulating a positron beam system by using SIMION”, Nuclear Science and Technology, Vol.7, No 3, 17-24, 2017 [12] David J Manura, David A Dahl, “SIMION Version 8.0/8.1 User Manual”, Document Revision 5, Scientific Instrument Services Inc., 2011 [13] Wolfgang Anwanda, Gerhard Brauer, Maik Butterling, Hans-Rainer Kissener, Andreas Wagner, “Design and Construction of a Slow Positron Beam for Solid and Surface Investigations”, Defect and Diffusion Forum, Vol 331, 2012 [14] Wolfgang Anwanda, Gerhard Brauer, Hans-Rainer Kissener, “Magnetically guided slow positron beam for defect studies”, Positron Group of TU Dresden at Research Centre Rossendorf, Dresden, Germany, 1994 [15] R Krause-Rehberg, “Simple design for a continuous magnetically guided positron beam – and – News from the EPOS project”,Institute of Physics, Martin Luther University HalleWittenberg, Germany, 2010 CONCEPTUAL DESIGN OF A SLOW POSITRON BEAM SYSTEM USING SIMION SIMULATION PROGRAM Abstract: The slow positron beam system is an important device in the study of positron physics and techniques, especially in materials research Our report presents the obtained results in studying and applying a charged particles trajectory simulation program – Simion to build feasible design models for the system Some positron beam trajectory calculation tests have been performed for comparison between the models From that, a conceptual design model has been proposed as a basis for building a real slow positron beam system in the future in Vietnam Keywords: Slow positron beam system, simulation program, Simion, conceptual design ... SPONSOR chứng tỏ chương trình Simion cơng cụ tính tốn thích hợp, sử dụng cho việc thiết kế nguyên tắc hệ thống chùm positron chậm Báo cáo trình bày kết nghiên cứu áp dụng chương trình mơ Simion để xây... mơ hình thiết kế khả thi hệ thống chùm positron chậm dựa nguyên lý thiết kế số hệ thống tiêu biểu đưa vào hoạt động giới từ đưa mơ hình thiết kế ngun tắc khả thi làm sở cho việc thiết kế chi tiết... xây dựng chương trình mơ Simion Vì hạn chế khả ứng dụng chương trình, mơ hình thiết kế dùng để mô nguyên lý hoạt động hệ thống chùm positron chậm cách đơn giản chưa phải thiết kế chi tiết kỹ thuật