Phương pháp Monte Carlo được xem là một phương pháp tính liều chính xác nhất trong xạ trị. Mô phỏng máy gia tốc chính xác là một yêu cầu cần thiết trong tính liều bằng phương pháp này. Trong nghiên cứu này, máy gia tốc Primus M5497 của hãng Siemens tại Bệnh viện Đa khoa Đồng Nai được mô phỏng bằng chương trình EGSnrc cho mức năng lượng photon 6 MV.
103 TẠP CHÍ PHÁT TRIỂN KHOA HỌC & CƠNG NGHỆ: CHUYÊN SAN KHOA HỌC TỰ NHIÊN, TẬP 2, SỐ 4, 2018 Mơ máy gia tốc tuyến tính phương pháp Monte Carlo dùng chương trình EGSnr Dương Thanh Tài, Hoàng Đức Tuân, Lương Thị Oanh, Trương Thị Hồng Loan, Nguyễn Đơng Sơn Tóm tắt—Phương pháp Monte Carlo xem phương pháp tính liều xác xạ trị Mơ máy gia tốc xác yêu cầu cần thiết tính liều phương pháp Trong nghiên cứu này, máy gia tốc Primus M5497 hãng Siemens Bệnh viện Đa khoa Đồng Nai mơ chương trình EGSnrc cho mức lượng photon MV BEAMnrc DOSXYZnrc hai chương trình sử dụng việc mơ tính tốn phân bố liều phantom nước Phần trăm liều theo độ sâu (Percentage depth dose, PDD) phân bố liều theo phương ngang (Beam profiles, OCR) có từ mô so sánh với liệu thực nghiệm để đánh giá độ xác mơ Kết thu có phù hợp tốt mô thực nghiệm với khác biệt PDD 1,26 % OCR nhỏ % Bên cạnh đó, chúng tơi sử dụng phương pháp đánh giá số Gamma chạy Matlab Phần trăm liều theo độ sâu có phần trăm số Gamma đạt 100% phân bố liều theo phương ngang có phần trăm số gamma đạt 98,5 % với yêu cầu sai biệt liều % độ lệch khoảng cách mm Kết cho thấy mô thành công máy gia tốc tuyến tính với độ xác cao Từ khóa—máy gia tốc, mô Monte Carlo, EGSnrc, BEAMnrc, DOSXYZnrc, số gamma Ngày nhận thảo: 29-7-2017; Ngày chấp nhận đăng: 18-12-2017; Ngày đăng:15-10-2018 Tác giả Dương Thanh Tài1,2,*, Hoàng Đức Tuân2,4, Lương Thị Oanh2,4, Trương Thị Hồng Loan2, Nguyễn Đông Sơn3 Bệnh viện Đa khoa Đồng Nai; 2Trường Đại học Khoa học Tự nhiên, ĐHQG-HCM; 3Công ty thiết bị y tế Chí Anh; 4Trường Đại học Nguyễn Tất Thành (Email: thanhtai_phys@yahoo.com) N MỞ ĐẦU gày nay, xạ trị ngành ứng dụng kỹ thuật hạt nhân vào y học mạnh mẽ phương pháp đóng vai trò quan trọng điều trị ung thư Một yêu cầu thiết yếu định đến thành công hay thất bại xạ trị tính tốn liều cho bệnh nhân trước xạ trị Tuy nhiên, việc tính liều có xác hay khơng lại phụ thuộc vào thuật tốn sử dụng để tính tốn Các thuật tốn tính liều phát triển mạnh mẽ từ năm 1950 [1] chia làm nhóm chính: (1) dựa hiệu chỉnh (correction-based), (2) dựa mơ hình hóa (model-based) (3) dựa nguyên lý (principle-based) Thuật tốn tính liều dựa hiệu chỉnh thuật tốn tính liều dựa vào giá trị nội suy ngoại suy từ thực nghiệm như: Phần trăm liều theo độ sâu (PDD) cho kích thước trường khác nguồn bề mặt định (Source to Surface Distance, SSD), phân bố liều theo phương ngang (OCR), tỉ số mô phantom (Tissue Phantom Ratio, TPR)… Sau thuật tốn hiệu chỉnh khác biệt điều kiện điều trị điều kiện đo lường Hiệu chỉnh bao gồm: Hiệu chỉnh suy giảm môi trường không đồng nhất, hiệu chỉnh tán xạ, kích thước trường, … Đối với mơi trường đồng nước, thuật toán cho kết xác Tuy nhiên, mơi trường khơng đồng thể (xương phổi…) thuật tốn xác [1] Thuật tốn tính liều dựa mơ hình hóa ngun lý vật lý sau đơn giản hóa q trình tương tác vật lý nhằm mơ 104 SCIENCE AND TECHNOLOGY DEVELOPMENT JOURNAL: NATURAL SCIENCES, VOL 2, ISSUE 4, 2018 vận chuyển xạ thực tế để đẩy nhanh tốc độ tính tốn Những trình tương tác vật lý đơn giản hóa phương trình tính tốn tích chập (convolution) siêu chồng chập (convolution -superposition) sử dụng môi trường khơng đồng cho độ xác cao Tính liều thuật tốn đưa vào phần mềm lập kế hoạch sử dụng thuật toán tính liều khác gồm: Pencil Beam Convolution (PBC), the Analytical Anisotropic Algorithm (AAA), (Varian Medical System, Inc Palo Alto, CA, USA) Collapse Cone Convolution (CCC) algorithms, (Pinnacle, CMS XiO) [9] Thuật tốn tính liều dựa ngun lý ứng dụng phương pháp Monte Carlo (MC) để mô vận chuyển số lượng lớn hạt photon hạt electron môi trường vật chất Phương pháp MC thường sử dụng công cụ để kiểm tra tính xác cho thuật tốn tính liều khác [2] Phương pháp mơ tả xác chất vật lý tương tác xem xét riêng cho hình học máy gia tốc, phận tạo chùm tia, bề mặt bệnh nhân không đồng mật độ, cho phép xử lý nhiều trường hợp tính liều phức tạp nên kết tính tốn phân bố liều xác Cũng việc tính tốn liều phương pháp MC tốn nhiều thời gian so với thuật tốn tính liều khác Để áp dụng phương việc mơ xác máy gia tốc yêu cầu thiết yếu Mục tiêu nghiên cứu mô hệ máy gia tốc Primus M5497 (của hãng Siemens) sử dụng điều trị Bệnh viện Đa khoa Đồng Nai (gọi tắt Bệnh viện) Hiện nay, chương trình áp dụng Monte Carlo gồm: PENELOPE, MCNP, GEANT4, GATE, EGSnrc… Mỗi chương trình có mạnh riêng ứng dụng nhiều lĩnh vực khác Trong đó, áp dụng nhiều lĩnh vực y khoa phải kể đến chương trình EGSnrc Chương trình EGSnrc phát triển với chương trình linh hoạt BEAMnrc, DOSXYZnrc [3, 4] Vì vậy, nghiên cứu này, chúng tơi sử dụng chương trình EGSnrc để mơ máy gia tốc tuyến tính dùng xạ trị hãng Siemens Bệnh viện Xử lý kết mô thực nghiệm công đoạn quan trọng ảnh hưởng đến kết q trình mơ Các cơng trình nghiên cứu trước mơ máy gia tốc chương trình EGSnrc [5] thường sử dụng sai số tương đối để đánh giá kết Tuy nhiên, việc sử dụng sai số tương đối (sai khác liều điểm) để đánh giá chưa đủ thiếu xác cho trường hợp tính liều vùng có liều cao vùng biến thiên liều [6] Chỉ số gamma đề xuất Low cộng (1998) [6] áp dụng nghiên cứu Ngoài ra, số chất lượng (beam quality) độ phẳng (flatness) chùm tia tính tốn để so sánh mô thực nghiệm VẬT LIỆU VÀ PHƯƠNG PHÁP Mô máy gia tốc chương trình BEAMnrc Sử dụng chương trình BEAMnrc để mô cho chùm photon mức MV hệ Với kích thước trường chiếu 10 × 10 cm2 khoảng cách từ nguồn đến mặt phẳng phantom SSD = 100 cm Tất vật liệu kích thước hệ máy gia tốc cung cấp từ nhà sản xuất Tất thành phần cần thiết đầu máy gia tốc mô chương trình BEAMnrc cách thiết lập số thành phần riêng lẻ gọi mô-đun (component module, CM), vuông góc với hướng chiếu chùm tia Cấu tạo máy gia tốc Bệnh viện gồm thành phần, thành phần mô tả CM BEAMnrc sau (hình 1): Hình Các thành phần máy gia tốc TẠP CHÍ PHÁT TRIỂN KHOA HỌC & CƠNG NGHỆ: CHUYÊN SAN KHOA HỌC TỰ NHIÊN, TẬP 2, SỐ 4, 2018 1) Cửa thoát (exit window) gồm lớp titan (Ti) dày 0,005 cm lớp nước dày 0,066 cm Cửa sổ chân khơng có vị trí bắt đầu -0,424 cm Được khai báo CM: SLABS; 2) Bia (target) gồm lớp theo thứ tự khơng khí dày 0,112 cm; vonfram (W) dày 0,064 cm, hợp kim đồng thau vàng (81,5 %); đồng (16,5 %) nicken (2 %) (Nicoro) dày 0,015 cm, đồng (Cu) dày 0,165 cm, Nicoro dày 0,005 cm, thép khơng rỉ dày 0,102 cm; than chì dày 1,016 cm; thép không rỉ dày 0,004 cm CM: FLATFILT; 3) Bộ lọc phẳng (flattening filter) khai báo CM: FLATFILT; 4) Buồng ion hóa (ionization chamber) gồm lớp: gốm (Al2O3) dày 0,152 cm xen kẽ với lớp nitrogen (N2) dày 0,184 cm CM: CHAMBER; 5) Gương (mirror) cấu tạo SiO2 0,209 cm CM: MIRROR; 6) Ngàm theo trục Y (JAWY) làm vonfram, độ mở theo trục y thay đổi để tạo kích thước trường 10 × 10 cm2 bề mặt Phantom, dày 7,620 cm CM: JAWS; 7) Ngàm theo trục X (JAWX) làm vonfram, độ mở theo trục x thay đổi để tạo kích thước trường 10 × 10 cm2 bề mặt phantom, dày 7,620 cm CM: JAWS; 8) Tấm mica (RECTICLE tray) làm mica dày 0,663 cm CM: SLABS; 9) Lớp khơng khí dày 56,805 cm CM: SLABS Các thông số mô AE = ECUT = 0,700 MeV, AP = PCUT = 0,010 MeV áp dụng nghiên cứu trước [7] Số lịch sử N = 109 hạt electron mô tiến hành chạy vi xử lý 2400 Intel (R) Core i5 Mơ tả xác nguồn electron đập vào bia yêu cầu cần thiết để mô xác máy gia tốc Trong nghiên cứu này, sử dụng nguồn số 19 thư viện nguồn [4] Các thông số nguồn số 19 công bố cơng trình trước [8] Đầu q trình chạy BEAMnrc để mơ máy gia tốc tuyến tính file khơng gian pha, có chứa đầy đủ thơng tin q trình chuyển động hạt lượng, vị trí, góc tới, hướng chuyển động,… File không gian 105 pha ghi nhận CM thứ với khoảng cách SSD = 100 cm tính từ nguồn chiếu Có thể sử dụng file khơng gian pha để phân tích chùm tia nguồn đầu vào cho q trình tính tốn phân bố liều phantom Tính tốn phân bố liều chương trình DOSXYZnrc Chương trình DOSXYZnrc sử dụng để tính tốn phân bố liều cho phantom nước (50 × 50 × 30 cm3) Khai báo nguồn nguồn thư viện DOSXYZnrc Nguồn dùng file không gian pha (*.egsphsp1) tạo từ q trình chạy BEAMnrc Nguồn tới từ phía trước theo phương z nằm mặt phantom Phantom chia thành × 31 × 66 voxels, trình bày hình Phantom nước đặt vị trí cho khoảng cách từ nguồn chiếu đến bề mặt phantom SSD = 100 cm Các electron photon có lượng tối thiểu (ECUT, PCUT) thiết lập 0,700 MeV 0,010 MeV tương ứng Các số liệu mô phân bố liều phantom nước từ trình chạy DOSXYZnrc tiến hành phân tích đánh giá độ sai biệt với thực nghiệm dựa vào sai số tương đối số Gamma tính code Matlab Khảo sát thơng số thực nghiệm Các số liệu thực nghiệm gồm phần đường trăm liều theo độ sâu (PDD), đường phân bố liều theo phương ngang (OCR) độ sâu 1,5 cm, cm, 10 cm, 20 cm thu hệ máy gia tốc tuyến tính Siemens Primus M5497 Hệ đo thiết lập Hình gồm: hai buồng ion hóa CC13 (IBA Dosimetry, Đức), đặt phantom nước lại để khơng khí Đầu dò đặt phantom di chuyển tới vị trí phantom Các đầu dò điều khiển hệ thống phần mềm điều khiển Omni Pro-Accept V7.4c (IBA Dosimetry, Đức) thông qua khối CU500E (IBA Dosimetry, Đức), có nhiệm vụ cung cấp điện áp ± 300 V cho hai đầu dò, điều khiển đầu dò đến vị trí cần đo SCIENCE AND TECHNOLOGY DEVELOPMENT JOURNAL: NATURAL SCIENCES, VOL 2, ISSUE 4, 2018 106 Trong đó: r , r D (r ) Dm (rm ) m c c c + ∆D = liều tính tốn với liều đo sai khác r2 r , r r r m c m c + ∆d = sai khác vị trí giá trị liều + rm vị trí điểm đo, rc vị trí khơng gian phân bố tính tốn ứng với điểm đo Tính thơng số chùm tia photon D10 Gi = Dmax Hình Khai báo phantom cho DOSXYZnrc Tính sai số tương đối Sai số tương đối mô thực nghiệm tính cơng thức sau: D D2 D 100% D1 (1) Trong đó: D1 giá trị liều đo thực nghiệm vị trí theo trục z, D2 giá trị liều mơ vị trí theo trục z (3) Thông số chất lượng chùm tia [9]: Trong đó: D10 liều độ sâu z = 10 cm, Dmax liều độ sâu cực đại - Độ phẳng chùm tia: Độ phẳng chùm tia F định nghĩa thay đổi lớn giá trị phần trăm liều theo phương ngang, vùng chiếm khoảng 80% kích thước trường chiếu tính từ vị trí trung tâm [9] F= Imax Imin 100% Imax + Imin (4) Trong đó: Imax Imin giá trị phần trăm liều hấp thụ lớn nhỏ vùng khoảng 80% kích thước trường chiếu Tính số Gamma Việc sử dụng sai khác liều điểm để đánh giá kết theo cách thơng thường chưa đủ gây sai số vùng liều thấp khu vực có độ biên thiên liều cao [6, 8] Do đó, số gamma phương pháp dùng nghiên cứu để đánh giá kết mô Chỉ số Gamma kết hợp sai số liều lượng ∆D phạm vi khoảng cách cho phép DTA (thông thường ∆D/DTA = 3%/3 mm) Chỉ số gamma tính theo cơng thức sau: rm , rc r rm , rc rm , rc d M2 DM2 (2) KẾT QUẢ VÀ THẢO LUẬN Đường phần trăm liều theo theo độ sâu (depth dose) Đường phân bố liều theo độ sâu chùm photon MV chuẩn hố độ sâu có giá trị liều cực đại 1,5 cm Hình trình bày đường phân bố liều theo độ sâu có từ thực nghiệm mơ chương trình EGSnrc; đường màu đỏ số gamma 107 TẠP CHÍ PHÁT TRIỂN KHOA HỌC & CÔNG NGHỆ: CHUYÊN SAN KHOA HỌC TỰ NHIÊN, TẬP 2, SỐ 4, 2018 Hình Hệ đo thực nghiệm (hình trái) sơ đồ kết nối (hình phải) Hình cho thấy có phù hợp tốt giá trị mô với giá trị thực nghiệm Theo thống kê, sai khác trung bình mô thực nghiệm 1,26 % Trong đó, sai khác lớn phát vùng 0,5 cm từ bề mặt nước 3,15 % vùng từ độ sâu z = 10 cm tới z = 30 cm 1,64 % Sai khác diện (a) vùng từ độ sâu cực đại (dmax) đến vị trí z = 10 cm 0,29 % (hình 2) Bên cạnh tất điểm có số gamma < (đường màu đỏ hình 4, bên trái) Phần trăm số gamma đạt yêu cầu sai biệt liều % độ lệch khoảng cách mm 100 % (b) (a) (b) Hình Kết mơ thực nghiệm đường cong phân bố liều theo độ sâu Chỉ số đánh giá chất lượng chùm tia photon đo tính theo cơng thức (3) cho trường hợp mô là: Qi = 0,668 % Chỉ số phẩm chất chùm tia photon máy LINAC khảo sát theo thực nghiệm Bệnh viện là: Qi = 0,67 % Sai khác tương đối 0,15 % Kết Hình cho thấy đường phân bố liều theo độ sâu mô phù hợp với thực nghiệm SCIENCE AND TECHNOLOGY DEVELOPMENT JOURNAL: NATURAL SCIENCES, VOL 2, ISSUE 4, 2018 108 Hình Sai số tương đối theo độ sâu Phân bố liều theo phương ngang (beam profile) Phân bố liều theo phương ngang tính độ sâu 1,5 cm; cm; 10 cm 20 cm so sánh với phân bố liều thực nghiệm độ sâu tương ứng Phân bố liều theo phương ngang độ sâu 1,5 cm Hình phân bố liều theo phương ngang mức lượng MV độ sâu cm tốt độ sâu 1,5 cm (Hình 6) (có 98% điểm có số gamma < với tiêu chí %/3 mm) Điều lý giải liều độ sâu gần bề mặt (a) nước dao động lớn độ sâu tăng Kết tính tốn độ phẳng chùm tia q trình mơ theo cơng thức (4) 1,12 % Độ phẳng chùm tia F đo thực nghiệm độ sâu tương ứng 1,19 % Vậy độ phẳng chùm tia F mô tương đồng với thực nghiệm Độ phẳng chùm tia mô độ sâu cm tốt so với độ sâu 1,5 cm Bề rộng penumbra có phù hợp tốt Phân bố liều có từ trình mơ vùng umbra có chút sai lệch nhỏ so với thực nghiệm, thấp so với thực nghiệm (chỉ số gamma lớn vị trí từ 22 – 24 cm) (b) (b) (a) Hình Phân bố liều theo phương ngang mức lượng MV độ sâu 1,5 cm 109 TẠP CHÍ PHÁT TRIỂN KHOA HỌC & CƠNG NGHỆ: CHUN SAN KHOA HỌC TỰ NHIÊN, TẬP 2, SỐ 4, 2018 (a) (b) (b) (a) Hình Phân bố liều theo phương ngang mức lượng MV độ sâu cm Phân bố liều theo phương ngang độ sâu 10 cm Kết tính tốn độ phẳng chùm tia q trình mơ theo cơng thức (4) 0,59 % Hình phân bố liều theo phương ngang độ sâu 10 cm theo kết tính độ phẳng thấy (a) có phù hợp tốt mô thực nghiệm, cụ thể số gamma trường hợp đạt 100 % với tiêu chí %/3 mm độ phẳng chùm tia F tốt (b) (b) (a) Hình Phân bố liều theo phương ngang mức lượng MV độ sâu 10 cm Phân bố liều theo phương ngang độ sâu 20 cm Kết tính tốn độ phẳng chùm tia q trình mơ theo cơng thức (4) 1,01 Tương tự trường hợp trên, phân bố liều theo phương ngang độ sâu 20 cm (Hình 9) phù hợp với thực nghiệm (phần trăm số gamma đạt tiêu chí %/3 mm 100 %) SCIENCE AND TECHNOLOGY DEVELOPMENT JOURNAL: NATURAL SCIENCES, VOL 2, ISSUE 4, 2018 110 (b) (a) (b) (a) Hình Phân bố liều theo phương ngang mức lượng MV độ sâu 20 cm Hình 10 thể tổng thể phân bố liều độ sâu khác phantom q trình mơ lý giải lý nhà vật lý y khoa thường chuẩn hóa, khảo sát liều độ sâu Kết cho thấy kết khảo sát liều theo phương ngang độ sâu khác phù hợp tốt với thực nghiệm Kết khảo sát phần trăm liều theo phương ngang độ sau khác cho thấy độ sâu 10 cm phù hơp tốt Kết góp phần (a) (a) (b) (b) (c) (c) (d) (d) Hình 10 Phân bố liều độ sâu khác phantom KẾT LUẬN Chúng mô thành công máy gia tốc tuyến tính Bệnh viện Đa khoa Đồng Nai với độ xác cao: sai số tương đối mô thực nghiệm PDD 0,7 % độ phẳng phân bố liều theo phương ngang nhỏ % Bên cạnh đó, chúng tơi sử dụng phương pháp đánh giá số Gamma chạy code Matlab Phần trăm liều theo độ sâu có phần trăm số Gamma đạt 100 % phân bố liều theo phương ngang có phần trăm số gamma đạt 98,5 % với yêu cầu sai biệt liều % độ lệch khoảng cách mm TẠP CHÍ PHÁT TRIỂN KHOA HỌC & CÔNG NGHỆ: CHUYÊN SAN KHOA HỌC TỰ NHIÊN, TẬP 2, SỐ 4, 2018 Lời cảm ơn: Một phần kết nghiên cứu báo cáo Hội nghị Vật lý y khoa Malaysia (IOS ISMP), ngày 26 tháng 08 năm 2016 Hội nghị Vật lý y khoa giới (ICMP) Thái Lan, ngày 12 tháng 12 năm 2016 TÀI LIỆU THAM KHẢO [1] L Lu, Dose calculation algorithms in external beam photon radiation therapy, Int J Cancer Ther Oncol, 1:01025, 2013 [2] F Verhaegen and J Seuntjens, Monte Carlo modelling of external radiotherapy photon beams, Phys Med Biol 48, R107–R164, 2003 [3] B Walters, I Kawrakow, and D.W.O Rogers, DOSXYZnrc User’s Manual, National Research Council of Canada Report, PIRS–794, 2004 [4] D.W.O Rogers, B Walters, and I Kawrakow, BEAMnrc Users Manual, National Research Council of Canada Report, PIRS –0509a, 2005 111 [5] P.T.T Lý, Tính liều hấp thụ gây chùm tia photon từ máy gia tốc dùng Monte Carlo code EGS, Luận văn Thạc sĩ, Trường Đại học Khoa học Tự nhiên, ĐHQG-HCM, 2009 [6] D.A Low, W.B Harms, S Mutic, J.A Purdy, A technique for the quantitative evaluation of dose distributions, Med Phys 25, 656–661, 1998 [7] Y Shi, L Zhou, X Zhen, S Zhang, The simulation of linear accelerator using BEAMnrc with DOSXYZnrc, IEEE 978, 4244–4713, 2010 [8] D.T Tai, N.D Son, T.T.H Loan, H.D Tuan, A method for determination of parameters of the initial electron beam hitting the target in linac, Journal of Physics: Conference Series, 851, 2017 [9] S Tung, Linac accelerator: Service instruction, Siemens OCS, TH003/02/I (2005), vùng khí hậu nơng nghiệp, 21, 2016 The simulation of a linear accelerator using Monte Carlo method with EGSnrc Program Duong Thanh Tai1,2,*, Hoang Duc Tuan2,4, Luong Thi Oanh2,4, Truong Thi Hong Loan2, Nguyen Dong Son3 Dong Nai General Hospital; 2VNUHCM-University of Science; 3Chi Anh Companny; 4Nguyen Tat Thanh University *Corresponding author: thanhtai_phys@yahoo.com Received: 29-7-2017, Accepted: 18-12-2017, Published:15-10-2018 Abstract—The Monte Carlo method is considered to be the most accurate algorithm for dose calculation in radiotherapy Linear accelerator accurate simulation is required as a prior condition for Monte Carlo dose calculation algorithm In this study, the MV photon beams from a Siemens Primus Linear Accelerator (LINAC) M5497 at the Dong Nai General Hospital was modelled by using EGSnrc The BEAMnrc DOSXYZnrc user code were used for simulation the head of LINAC and the dose distribution in water phantom The percentage depth dose (PDD) and beam profiles (OCR) were calculated and then compared with the measured ones in order to evaluate the simulation accuracy Excellent agreement was found between simulations and measurements with an average difference of 1.26 % for PDD, less than % for OCR In addition, the gamma evaluation method for simulation was also performed using an in-house Matlab code The percentage gamma passing rate was 100% for PDD and 98.5 % for OCR with % dose difference and mm distance to agreement as acceptance criteria The result showed that we had successfully simulated LINAC with excellent agreement Index Terms—linear accelerator, Monte Carlo simulation, EGSnrc, BEAMnrc, DOSXYZnrc, Gamma index ... tơi sử dụng chương trình EGSnrc để mơ máy gia tốc tuyến tính dùng xạ trị hãng Siemens Bệnh viện Xử lý kết mô thực nghiệm công đoạn quan trọng ảnh hưởng đến kết trình mơ Các cơng trình nghiên... GATE, EGSnrc… Mỗi chương trình có mạnh riêng ứng dụng nhiều lĩnh vực khác Trong đó, áp dụng nhiều lĩnh vực y khoa phải kể đến chương trình EGSnrc Chương trình EGSnrc phát triển với chương trình. .. chùm tia tính tốn để so sánh mơ thực nghiệm VẬT LIỆU VÀ PHƯƠNG PHÁP Mô máy gia tốc chương trình BEAMnrc Sử dụng chương trình BEAMnrc để mô cho chùm photon mức MV hệ Với kích thước trường chiếu