Khuếch đại của tốc độ phản ứng tổng hợp hạt nhân trong môi trường plasma OCP đậm đặc

pdf 11 trang ngocly 1880
Bạn đang xem tài liệu "Khuếch đại của tốc độ phản ứng tổng hợp hạt nhân trong môi trường plasma OCP đậm đặc", để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên

Tài liệu đính kèm:

  • pdfkhuech_dai_cua_toc_do_phan_ung_tong_hop_hat_nhan_trong_moi_t.pdf

Nội dung text: Khuếch đại của tốc độ phản ứng tổng hợp hạt nhân trong môi trường plasma OCP đậm đặc

  1. Tạp chí KHOA HỌC ĐHSP TP HCM Đỗ Xuân Hội, Lý Thị Kim Thoa ___ KHUẾCH ĐẠI CỦA TỐC ĐỘ PHẢN ỨNG TỔNG HỢP HẠT NHÂN TRONG MÔI TRƯỜNG PLASMA OCP ĐẬM ĐẶC ĐỖ XUÂN HỘI *, LÝ THỊ KIM THOA TÓM TẮT Thế màn chắn ở khoảng cách liên hạt nhân rất nhỏ đóng vai trò quan trọng là làm tăng tốc độ phản ứng tổng hợp hạt nhân trong những thiên thể đậm đặc như sao lùn trắng hay sao neutron. Trong bài báo này, sau khi đã thảo luận bao quát và chi tiết về những công trình mới nhất liên quan đến việc khuếch đại của tốc độ này theo giá trị của tham số tương liên, chúng tôi đề nghị ba công thức nhằm mục đích cải tiến độ chính xác của các biểu thức xấp xỉ của thế màn chắn ở khoảng cách gần bằng không. Các kết quả thu được có tính đến sự liên tục giữa đặc tính của plasma OCP và plasma hỗn hợp nhiều thành phần cũng như quy tắc hỗn hợp tuyến tính nổi tiếng. ABSTRACT Magnifiying synthetic nuclear reaction rate in the condensed plasma OCP The screening potential at the tiny values of the antinuclear distance plays an important role in increasing the rate of the synthetic nuclear reaction in condensed celestial bodies such as the white dwarfs or the neutron stars. In this paper, after referring the most recent works concerning the magnification of this rate due to the values of the coupling parameter, we propose three formulae with the aim to improve the precision of the approximate expressions for the screening potential at near zero distance. Those results obtained were taken into account the continuity between the characteristics of OCP plasma and the multi component mixture plasma as well as the well known linear mixing rule. 1. Mở đầu Nguồn năng lượng chính được bức xạ từ các sao trong vũ trụ có nguồn gốc là phản ứng tổng hợp hạt nhân. Các phản ứng này ảnh hưởng đến quá trình tiến hóa của những thiên thể tạo bởi plasma có mật độ khối lượng cao như sao lùn trắng hoặc sao neutron. Trong plasma, các hạt nhân có thể vượt qua hàng rào thế Coulomb do hiệu ứng đường ngầm lượng tử để gây ra phản ứng tổng hợp. Tốc độ phản ứng hạt nhân được tính bởi hệ thức tổng quát: R  () v vNi N j , trong đó,  ()v là tiết diện tán xạ, hàm theo tốc độ tương đối v giữa các hạt nhân tương tác, và Ni , N j là các mật độ hạt. Bản thân  ()v được biểu thị qua thừa số vật * TS, Trường Đại học Quốc tế, ĐHQG TP HCM CN, Trường THPT Tân Bình, TP HCM 69
  2. Tạp chí KHOA HỌC ĐHSP TP HCM Số 21 năm 2010 ___ lí thiên văn S() , là hàm theo năng lượng  , đặc trưng cho tương tác hạt nhân giữa hai ion, theo: S()  (  ) exp G .    exp G Ở hệ thức trên, số hạng cho biết xác suất truyền qua hàng rào thế  Coulomb của hai hạt nhân, với  G là năng lượng Gamow [4]. Trong một công trình tiên phong trong lĩnh vực vật lí thiên văn hạt nhân (nuclear astrophysics) [14], Salpeter đã khẳng định rằng khi mật độ khối lượng của hệ plasma đủ lớn để ảnh hưởng màn chắn của các hạt nhân xung quanh hai hạt nhân tương tác trở nên đáng kể, tốc độ phản ứng hạt nhân ở trên được nhân thêm hệ số khuếch đại eh0 / kT , đặc trưng cho tác dụng của môi trường bên ngoài lên tương tác giữa hai hạt nhân này. Hệ số h0 là thế màn chắn ở khoảng cách liên ion gần bằng không (là khoảng cách tại đó phản ứng tổng hợp hạt nhân xảy ra), k và T lần lượt là hằng số Boltzmann và nhiệt độ của môi trường. Việc xác định giá trị của hệ số h0 có tầm quan trọng đặc biệt, liên quan đến lĩnh vực vật lí hạt nhân để tính toán hiệu suất phản ứng tổng hợp hạt nhân trong thiên thể và đồng thời nhằm đưa ra một biểu thức của thế màn chắn cho mọi giá trị của khoảng cách liên ion và cho mọi tương quan giữa thế năng tương tác Coulomb và năng lượng chuyển động nhiệt. Chẳng hạn, trong môi trường plasma liên kết mạnh, hệ số khuếch đại được đánh giá lớn đến cỡ 1074 . Trên thế giới, cho đến thời gian gần đây, đã có nhiều công trình liên quan đến vấn đề tính hệ số h0 (xem, ví dụ như, [2], [3]). Trong bài báo này, dựa trên phương pháp tham số hóa hiệu ứng trật tự địa phương đã được trình bày trong một công trình trước [8], chúng tôi sẽ đề xuất một biểu thức giải tích cho hệ số h0 tương thích với các kết quả mô phỏng Monte Carlo (MC) cho plasma liên kết mạnh, mô hình vật lí phù hợp với quá trình tiến hóa xảy ra trong sao lùn trắng và sao neutron và đồng thời đáp ứng được các yêu cầu về vật lí của bài toán như tính chẵn và luân phiên dấu của hệ thức Widom, tính chính xác của hệ số Jancovici, quy tắc hỗn hợp của plasma đa ion liên kết mạnh cũng như tính liên tục giữa plasma một thành phần và plasma nhiều thành phần. 2. Mô hình khảo sát Trong khuôn khổ của bài báo này, ta sẽ chỉ đề cập đến việc xác định sự khuếch đại của phản ứng hạt nhân trong môi trường plasma một thành phần (OCP – One Component Plasma) cổ điển, gồm các hạt mang điện tích dương Ze như nhau, đặt trong một biển đồng nhất electron mang điện tích e có tác dụng trung hòa điện. Mô hình OCP được xem là thích hợp để diễn tả chẳng hạn như sự đốt cháy carbon: 12CC 12 , phản ứng tổng hợp xảy ra trong sao lùn trắng có mật độ khối 70
  3. Tạp chí KHOA HỌC ĐHSP TP HCM Đỗ Xuân Hội, Lý Thị Kim Thoa ___ lượng lớn đến khoảng 109g / cm 3 và nhiệt độ đạt đến 108 K . Hiệu ứng lượng tử liên quan đến tốc độ phản ứng tổng hợp hạt nhân trong trường hợp đang xét có tác dụng nhất định khi “làm nhòe” (quantum smearing) vùng phân kì của thế Coulomb, sẽ được khảo sát trong một công trình khác. Để đặc trưng cho sự tương quan giữa năng lượng chuyển động nhiệt có độ lớn vào cỡ kT của các ion và tương tác đẩy Coulomb giữa các ion mang điện tích cùng dấu Ze, ta sử dụng tham số tương liên quen thuộc, được định nghĩa bởi: Ze 2  , akT trong đó, a là bán kính khối cầu ion, được tính theo mật độ hạt n như: 1/3 3 a . 4 n Với định nghĩa trên, ta có thể hình dung được tính chất của plasma phụ thuộc vào độ lớn của tham số tương liên : khi chuyển động nhiệt chiếm ưu thế, môi trường plasma sẽ trạng thái lưu chất và ngược lại, nếu tương tác Coulomb quan trọng hơn, ta sẽ có plasma kết tinh. Giá trị ngưỡng của , tại đó có sự chuyển pha từ lưu chất sang tinh thể lập phương tâm khối (bcc) được đánh giá vào cỡ: m 172 [15]. Nếu kí hiệu H(R) là hàm đặc trưng cho tác dụng của màn chắn của môi trường xung quanh lên hai ion tương tác nhau, thì với R là khoảng cách liên ion, hai ion này sẽ chuyển động trong trường thế hiệu dụng Ze 2 VRHR()() , (1) R 2 R Ze 1 hoặc nếu tính theo đơn vị r và , ta viết: V()() r H r . a a r Xác suất tương tác (contact probability) giữa hai ion cho bởi hàm phân bố xuyên tâm [2]: a 2 A ( Ze )2  g( r ) exp  H ( r ) exp P ( ) d  , 3/2 R ()kTh kT trong đó, A là số khối hạt nhân tính theo đơn vị khối lượng nguyên tử (amu) và P() r là đại lượng đặc trưng cho xác suất truyền qua hàng rào thế Coulomb. Nếu gọi t là điểm giật lùi (turning point) trong lí thuyết tán xạ cổ điển, ta có: 2 A rt P()() dr V r  . h 0 71
  4. Tạp chí KHOA HỌC ĐHSP TP HCM Số 21 năm 2010 ___ 2 r ()Ze Giá trị của t được xác định bởi điều kiện:  . rt 3. Khuếch đại tốc độ phản ứng hạt nhân trong plasma OCP liên kết mạnh Trong những sao lùn trắng hoặc sao neutron, độ lớn của tham số tương liên được đánh giá vào khoảng 100 nên vấn đề tính xác định giá trị của hệ số h0 là rất quan trọng. Để tính các giá trị của hệ số này tương ứng với các tham số tương liên khác nhau, ta có thể sử dụng một trong hai phương pháp: Thu thập dữ liệu cho hàm phân bố xuyên tâm g() r từ các mô phỏng MC và aln g ( r ) suy ra giá trị cho thế màn chắn H() r từ hệ thức: H() r . R  Mặt khác, tính đối xứng của bài toán cho phép ta viết thế màn chắn H() r dưới dạng đa thức Widom có bậc chẵn, luân phiên dấu [16]: 2 4 6 H( r ) h0 h 1 r h 2 r h 3 r , (2) trong đó, hệ số Jancovici h1 đã được chứng minh có giá trị chính xác: h1 0.25 [10]. Vì các mô phỏng MC không thể cho được các giá trị của g() r ở các khoảng cách liên ion quá nhỏ nên ta bắt buộc phải sử dụng các phép tính ngoại suy để có được các biểu thức của h0 tương ứng với các giá trị khác nhau của tham số . Khai thác phương pháp tham số hóa hiệu ứng trật tự địa phương [8], các hệ số của đa thức H(r) được biểu thị theo các hàm phụ thuộc tham số  như sau 5 k hi  a k ln  , (3) k 0 rất thuận tiện khi muốn thực hiện các chương trình trên máy tính. Các hệ số ak trong khai triển đa thức hi đều đã được cho trong một bảng. Cách tiếp cận khác để tính được các giá trị của hệ số h0 là sử dụng các hàm nhiệt động lực: Hệ plasma được xem như gồm N 2 ion và một phân tử lưỡng nguyên tử tạo bởi hai ion gần nhau đến một khoảng cách nào đó. Ta có thể chứng minh rằng h0 là hiệu số giữa năng lượng tự do của hệ plasma trước và sau phản ứng. Sử dụng quy tắc hỗn hợp tuyến tính (linear mixing rule), ta có được: 5/3 h0 2 f () f (2 ), với kí hiệu f () là năng lượng tự do của từng ion theo đơn vị của kT, và số hạng 25/3  tương ứng với tham số tương liên của hệ lưỡng nguyên tử [17]. Do giá trị của năng lượng tự do f () được cung cấp trực tiếp từ các mô phỏng MC, ta có thể suy ra các hệ số h0. Sự tương đồng giữa các kết quả thu được cho h0 từ hai phương pháp trên đã được phân tích trong một công trình gần đây [1]. 72
  5. Tạp chí KHOA HỌC ĐHSP TP HCM Đỗ Xuân Hội, Lý Thị Kim Thoa ___ Để có được hệ số h0 cho plasma lưu chất loãng, các tác giả H. DeWitt và W. 8 Slattery [6] đã sử dụng các mô phỏng MC với 5 10 cấu hình cho 1 000 ion với cách tiếp cận nhiệt động lực ở trên để thu được biểu thức:  h h() lm DWS (4a) 0DWS 0 DWS 100 với: 1,020822 1 h( lm ) 1,056349 0,274823ln  1,084312 , (4b) 0DWS 0,676936  DWS 2,7 ln  4,8 . (4c) Trong công thức (4a), các hàm h0DWS () lm xuất phát từ kết quả của quy tắc xấp xỉ hỗn hợp tuyến tính đã trình bày ở trên và hàm  DWS là kết quả của hiệu chính từ phép tính cho plasma nhiều thành phần. Các hệ thức xấp xỉ (4a, b, và c) ở trên tương đối phù hợp với các dữ liệu thu được từ các mô phỏng MC với phương pháp hỗn hợp tuyến tính thực hiện bởi chính các tác giả này, nhưng chỉ trong khoảng  [1, 20]. Một hệ thức tương tự với (4a, b, và c) cũng đã được các tác giả trên nêu ra trong một công trình khác [5]: 1,039957 1 1 h 1,056299 0,274823ln  1,084319 0,0271ln  0,048 ,(4d) 0DWS 0,676936   và cũng cho kết quả có sai số nhỏ khi so sánh với các dữ liệu MC cung cấp bởi chính công trình này. 1/2 Nhằm mục đích có được đồng thời kết quả h0 3  ở giới hạn của chế độ nhiệt hạt nhân cổ điển tương ứng với  1 theo [14] và quy luật h0 const đối với lưu chất Coulomb ( : 1), L. R. Gasque et al đã đề nghị [9]: 1,07541/ 2 h0G 1/4 . (5) 4 1,0754 2  3 Tuy nhiên, hệ thức trên chỉ cho giá trị tương đối chính xác của h0 với  80 , còn đối với các giá trị khác của , sai số là khá đáng kể. Theo các tác giả, các sai số trên là chấp nhận được nếu so sánh với sai số do phép tính thừa số vật lí thiên văn S(). 16 Gần đây hơn, khi quan tâm đến sự tổng hợp của các hạt nhân 12C và O trong các sao, trên cơ sở kết hợp phương pháp gần đúng WKB cho hiệu ứng đường ngầm lượng tử xuyên qua hàng rào thế Coulomb tạo bởi hai hạt nhân tương tác và phương pháp thế của trường trung bình tĩnh (static mean field potential) do các ion lân cận, 73
  6. Tạp chí KHOA HỌC ĐHSP TP HCM Số 21 năm 2010 ___ các tác giả A. I. Chugunov và H. E. DeWitt đã thực hiện 129 mô phỏng MC để thu được các dữ liệu cho plasma BIM được xem như chính xác nhất cho đến nay [2]. Trường hợp plasma OCP được suy ra từ hệ thức tổng quát, cho thấy ở vùng giá trị của tham số tương liên mà ta quan tâm trong khuôn khổ của bài báo này, kết quả thu được là hoàn toàn tương thích với thông báo trước đó thực hiện riêng biệt cho plasma OCP [3], hoặc với kết quả có được từ các mô phỏng Monte Carlo sử dụng tích phân lộ trình (PIMC - Path Integral Monte Carlo) do Militzer và Pollock thực hiện [11]. Theo [3], sự phụ thuộc của hệ số h0 vào độ lớn của tham số tương liên được cho bởi: ABAB   h 1/ 2 1 3 1 3 (6) 0CHU 1  BB  2 A2  2 4 với: A1 2,7822 , A2 98,34, AAA3 3 1 / 2 1,4515, B1 1,7476, B2 66,07 , B3 1,12, và B4 65. Đặc điểm của hệ thức trên, tương tự như trong [9], là ta có thể suy ra giá trị của h0 tiệm cận với lí thuyết Debye-Hückel vận dụng cho plasma liên kết loãng: 1/2 h0CHU 3  (  1). Khảo sát sự tương hợp giữa hệ thức (6) ở trên với các dữ liệu MC của H. DeWitt và W. Slattery [5], ta thấv sai số của h0 là khoảng 0,13% cho các giá trị của  10 , hoàn toàn chấp nhận được so sánh với sai số của mô phỏng MC. Trên cơ sở đó, ta có thể xem như hệ thức (6) cho ta các giá trị của h0 tương ứng với mọi giá trị của tham số . Bảng B.I. So sánh các giá trị h0 trong các bài báo [7] và [5]  h0MC96 [7] h0MC99 [5] 10 1,0962 1,0994 20 1,0962 1,0953 40 1,0901 1,0879 80 1,0828 1,0803 160 1,0770 1,0737 Để có thêm thông tin về mức độ chính xác của các số liệu cung cấp bởi các mô phỏng MC, chúng ta có thể tham khảo bảng B.I: Ở các cột thứ hai và thứ ba lần lượt là các giá trị của hệ số h0 cung cấp bởi các công trình [7] và [5]. Các giá trị của h0MC96 ở cột thứ hai được chúng tôi suy ra bằng phương pháp xấp xỉ trực tiếp từ các dữ liệu cho hàm phân bố xuyên tâm g(r). Ta có thể nhận thấy các giá trị này có độ lệch lớn nhất là 3,3 10 3 tương ứng với  160 . Điều đáng lưu ý là khi thực hiện các phép tính xấp xỉ trên, chúng tôi đã kiểm nghiệm được các hàm thế màn chắn 74
  7. Tạp chí KHOA HỌC ĐHSP TP HCM Đỗ Xuân Hội, Lý Thị Kim Thoa ___ h(r) theo khoảng cách liên ion r là các đa thức bậc chẵn, luân phiên dấu theo đúng đa thức Widom (2), đồng thời, các hệ số h1 trong khai triển (2) có giá trị gần đúng: h0 , 25 1 , như đã chứng minh bởi Jancovici [10]. Trong bài báo này, chúng ta sẽ sử dụng các giá trị của hệ số h0 cho bởi [5], tức là các giá trị ở cột thứ ba của bảng B.I, do sự tương thích giữa các số liệu này với hệ thức (6) như đã trình bày ở trên. Để thấy rõ mối quan hệ trong việc tính toán hệ số khuếch đại phản ứng hạt nhân trong trường hợp tổng quát cho plasma nhiều thành phần (MCP - Multicomponent Plasmas) và plasma OCP, đồng thời để có sự tương thích với phương pháp sử dụng quy tắc hỗn hợp tuyến tính, chúng tôi đề nghị hệ thức sau cho h0:  h h() lm (7a) 0 0 100 trong đó: 1,039957 1 h( lm ) 1,056299 0, 274823ln  1,084319 (7b) 0 0,676936  5 a (ln )k và:   k  , (7c) k 0 Các hệ số ak được cho bởi bảng B.II dưới đây : Bảng B.II. Các hệ số của công thức (7c) a0 a1 a2 a3 a4 a5 6,69370 0,69922 2,80549 1,95369 0,43372 0,03298 ) 0MC99 h - 0 0 h h ( 3 10 H.2. Đồ thị sai số giữa hệ thức (4d) và dữ liệu MC cùng từ công trình [5]. H.1. Biến thiên của h0 theo ln. Chấm tròn là dữ liệu MC [5]. Đường liền nét là hệ thức đề nghị (7a, b, và c). Đường đứt nét là hệ thức (4d) từ [5], và đường gạch chấm là hệ thức (5) từ [9]. 75
  8. Tạp chí KHOA HỌC ĐHSP TP HCM Số 21 năm 2010 ___ Bảng trên cho các giá trị của ak khác với ở [8] và cho giá trị h0 gần với dữ liệu mô phỏng MC chính xác nhất cho tới nay [2, 3, và 5]. Điều này có thể giải thích được là vì trong [8], các hệ số ak được tính theo kết quả mô phỏng MC cung cấp bởi H. DeWitt et al [7]. Như có thể thấy được trên hình H.1, các hệ thức xấp xỉ (7a, b, và c) ưu tiên cho vùng plasma đậm đặc, có tham số tương liên  1, là nội dung chủ yếu của bài báo này. Giá trị cực đại của hệ số h0 là vào khoảng 1,0994 tương ứng với  9,6 . Trên hình H.2, ta thấy sai số phạm phải giữa hệ thức (4d) và các dữ liệu MC [5] lớn nhất tương ứng với  10 và giảm dần khi  tăng, trong khi theo tính toán của chúng tôi, sai số giữa các hệ thức đề nghị (7a, b, và c) và cũng với các dữ liệu MC này là không đáng kể. Đồng thời cũng cần chú ý rằng, theo (7c), dáng điệu biến thiên của hàm  lệch khá xa khỏi dạng tuyến tính (4c) đề nghị bởi [6]. Các giá trị của h0 tương ứng với một số giá trị của tham số  được trình bày trên bảng B.IV: Ở cột thứ hai, ta có các giá trị cho bởi mô phỏng MC từ công trình [5], trong khi ở các cột thứ ba, thứ tư, và thứ năm lần lượt là các giá trị có được từ các hệ thức (4d), (5), và (7a, b, và c). Ở cột thứ sáu, ta có các giá trị của h0 từ công thức xấp xỉ đơn giản sau đây, được thiết lập nhằm mục đích tạo nên sự tương thích với công thức tính thế màn chắn đã thu được trong [8], và cũng để thuận tiện trong việc thực hiện chương trình trên máy tính: 5 h b (ln )k 0  k  , (8) k 0 với các hệ số bk cho bởi bảng B.III: Bảng B.III. Các hệ số của công thức (8) b1 b2 b3 b4 b5 b0 0,9450000 0,1993204 0, 0959109 0,0218715 0,0025140 0,0001177 Bảng B.IV. Các giá trị của hệ số h0 trích từ [5, 9] và từ công trình này  h0MC99 [5] h0DWS (4d) h0G (5) h0 (7a, b, c) h0 (8) 1 0,9450 0,9639 1,0388 0,9450 0,9450 10 1,0994 1,0924 1,0750 1,0994 1,0994 20 1,0953 1,0913 1,0753 1,0953 1,0953 40 1,0879 1,0858 1,0754 1,0879 1,0879 80 1,0803 1,0792 1,0754 1,0803 1,0803 160 1,0737 1,0731 1,0754 1,0737 1,0737 76
  9. Tạp chí KHOA HỌC ĐHSP TP HCM Đỗ Xuân Hội, Lý Thị Kim Thoa ___ Tương tự như với các hệ thức (7a, b, và c), sai số giữa công thức (8) và các dữ liệu MC cho bởi công trình [5] là rất nhỏ như ta có thể nhận xét ở bảng B.IV. ) 2 h 2MC96 h 100 - 2 h ( 3 10 H.4. Sai số giữa hệ thức (10) và dữ liệu MC [7] H.3. Biến thiên của h2 theo tham số tương liên Ở mức độ chính xác cao hơn, thay vì chỉ quan tâm đến giá trị của h0 H( r 0) , ta có thể viết biểu thức của thế màn chắn dưới dạng 2 4 H( r ) h0 0,25 r h 2 r , (9) tức là để khảo sát ảnh hưởng của khoảng cách liên ion r lên hiệu suất phản ứng. Dựa trên số liệu MC, một số tác giả đã đề nghị các giá trị cho hệ số h2. Ví dụ như các tính toán của Ogata et al [12] cho thấy h2 0 , hoặc từ phép tính xấp xỉ tuyến 2 tính H() r C0 C 1 r , Rosenfeld đã tìm thấy giá trị h2 3 10 [13]. Chú ý rằng các giá trị được đề nghị ở trên đều độc lập đối với tham số . Tuy nhiên, trên hình H.3, ta có thể thấy sự biến thiên khá nhạy của h2 theo giá trị của tham số tương liên. Vì vậy, để đạt độ chính xác cao cho thế màn chắn và đồng thời để nhấn mạnh sự phụ thuộc của h2 vào , là đại lượng đặc trưng cho tính chất của plasma, từ các số liệu MC cung cấp bởi các tác giả H. DeWitt và W. Slattery [5], ta có thể sử dụng hệ thức sau: 2 h2 5,61765 0,41802ln  10 . (10) Ta thấy có sự chênh lệch giữa các giá trị của h2 tìm được và giá trị cho bởi [13]. Ngược lại, giá trị bằng không của hệ số h2 đề nghị bởi [12] là quá sơ lược và không phản ánh được sự biến thiên của thế màn chắn theo tính chất của plasma ở khoảng cách liên ion không quá nhỏ. Hình H.3 cho thấy giá trị của hệ số h2 giảm dần theo tham số . Độ lệch giữa hệ thức (10) và số liệu MC cho bởi [7] được trình bày trên hình H.4. Ta thấy sai số phạm phải trong phép tính xấp xỉ có cực đại là 2 khoảng 1, 8 10 tương ứng với  160. Tuy nhiên, sai số này là chấp nhận được r4 nếu ta nhớ lại rằng trong đa thức Widom (2), hệ số h2 này được nhân với và ta đang chú ý đến khoảng cách r nhỏ đối với phản ứng tổng hợp hạt nhân. 77
  10. Tạp chí KHOA HỌC ĐHSP TP HCM Số 21 năm 2010 ___ 4. Kết luận Hệ số khuếch đại của phản ứng tổng hợp hạt nhân trong môi trường plasma đậm đặc như sao lùn trắng hoặc sao neutron chịu ảnh hưởng bởi hai yếu tố chính: thừa số vật lí thiên văn S() và thế màn chắn H(r) do tác dụng của hệ nhiều hạt dưới điều kiện mật độ khối lượng lớn. Trong bài báo này, sau khi tham khảo kết quả của những công trình liên quan mới nhất ở mức độ quốc tế, các tác giả đề nghị các biểu thức (7a, b, và c), (8), và (10) để đánh giá ảnh hưởng của điều kiện thứ hai kể trên. Các biểu thức được trình bày có quan tâm đến tính liên tục của hệ plasma một thành phần OCP với hệ plasma nhiều thành phần MCP. Những điểm tương đồng hay dị biệt với các kết quả thu được gần đây nhất cũng được phân tích, cho thấy các hệ thức đề nghị ở trên hoàn toàn tương thích với những kết quả của mô phỏng MC tin cậy nhất cho đến nay. Ngoài ra, tính đúng đắn của đa thức Widom cũng như của hệ số Jancovici cũng được kiểm nghiệm khi thực hiện các phép tính xấp xỉ trực tiếp trên các số liệu mô phỏng MC của hàm phân bố xuyên tâm. Tính tiệm cận với lý thuyết Debye-Hückel áp dụng cho plasma loãng cũng như ảnh hưởng của hiệu ứng màn chắn điện tử (electron screening) cũng là những hướng được ngắm đến nhằm có bức tranh toàn diện về tác dụng của môi trường lên hiệu suất của phản ứng tổng hợp áp suất hạt nhân (pycnonuclear reactions). TÀI LIỆU THAM KHẢO 1. Chugunov A. I. and DeWitt H. E. (2009), “Corrections to Linear Mixing in Binary Ionic Mixtures and Plasma Screening at Zero Separation”, Contrib. Plasma Phys., 49(10), pp. 696 - 699. 2. Chugunov A.I., DeWitt H.E. (2009), “Nuclear fusion reaction rates for strongly coupled ionic mixtures”, Phys. Rev. C, 80(1), pp. 014611-1- 014611- 12 3. Chugunov A. I., DeWitt H. E., Yakovlev D. G. (2007), “Coulomb tunneling for fusion reactions in dense matter: Path integral Monte Carlo versus mean field”, Phys. Rev. D, 76(2), pp. 025028-1-025028-13. 4. Clayton D. D. (1983), Principles of Stellar Evolution and Nucleosynthesis, The University of Chicago Press. 5. DeWitt H. and Slattery W. (1999), “Screening Enhancement of Thermonuclear Reactions in High Density Stars”, Contrib. Plasma Phys., 39(1–2), pp. 97 -100. 6. DeWitt H. and Slattery W. (1998), PNP-9 Workshop, Rostok Germany. 7. DeWitt H., Slattery W., and Chabrier G., (1996), “Numerical simulation of strongly coupled binary ionic plasmas”, Physica B, 228(1-2), pp. 21-26. 78
  11. Tạp chí KHOA HỌC ĐHSP TP HCM Đỗ Xuân Hội, Lý Thị Kim Thoa ___ 8. Đỗ Xuân Hội (2002), “Thế màn chắn trong plasma với tham số tương liên  [5, 160]”, Tạp chí Khoa học Tự nhiên ĐHSP TP HCM, (28), tr. 55-66. 9. Gasques L. R., Afanasjev A. V., Aguilera E. F., Beard M., Chamon L. C., Ring P., Wiescher M., and Yakovlev D. G. (2005), “Nuclear fusion in dense matter: Reaction rate and carbon burning”, Phys. Rev. C, 72(2), pp. 025806-1-025806- 14. 10. Jancovici B. (1977), “Pair correlation function in a dense plasma and pycnonuclear reactions in stars”, J. Stat. Phys., 17(5), pp. 357-370. 11. Militzer B., Pollock E. L. (2005), “Equilibrium contact probabilities in dense plasmas”, Phys. Rev. B 71(13), pp. 134303-1-134303-10. 12. Ogata S., Iyetomi H., and Ichimaru S. (1991), “Nuclear reaction rates in dense carbon-oxygen mixtures”, Astrophys. J., 372, pp. 259-266. 13. Rosenfeld Y. (1992), “Screening potentials in strongly coupled plasmas: Reanalysis of recent highly accurate simulations”, Phys. Rev. A, 46(2), pp. 1059-1065. 14. Salpeter, E. E. (1954), “Electron Screening and Thermonuclear Reactions”, Australian J. Phys., 7, pp. 373-388. 15. Slattery W. L., Doolen G. D., and De Witt H. E. (1982), “N dependence in the classical one-component plasma Monte Carlo calculations”, Phys. Rev. A, 26(4), pp. 2255–2258. 16. Widom B. (1963), “Some Topics in the Theory of Fluids”, J. Chem. Phys., 39(11), pp. 2808-2812. 17. Xuan Hoi Do (1999), Thèse de Doctorat de l’Université Paris 6 –Pierre et Marie Curie, Paris (Pháp). 79