BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH PHAN CƠNG THÀNH NHIỆT ĐỘNG LỰC HỌC CỦA PLASMA Ở TRẠNG THÁI KẾT TINH LUẬN VĂN THẠC SĨ VẬT LÝ Thành phố Hồ Chí Minh - 2011 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH PHAN CÔNG THÀNH NHIỆT ĐỘNG LỰC HỌC CỦA PLASMA Ở TRẠNG THÁI KẾT TINH Chuyên ngành: VẬT LÝ NGUYÊN TỬ HẠT NHÂN VÀ NĂNG LƯỢNG CAO Mã số: 60 44 05 LUẬN VĂN THẠC SĨ VẬT LÝ NGƯỜI HƯỚNG DẪN KHOA HỌC: TS ĐỖ XUÂN HỘI Thành phố Hồ Chí Minh - 2011 LỜI CẢM ƠN Với lịng kính trọng biết ơn xin bày tỏ lời cảm ơn chân thành tới: Phịng Khoa học cơng nghệ sau đại học, khoa Vật lý trường Đại học Sư phạm TP Hồ Chí Minh tạo điều kiện thuận lợi cho tơi suốt q trình học tập thực luận văn Các thầy giúp tiếp cận kiến thức khoa học suốt trình học tập Tiến sĩ Đỗ Xuân Hội, người thầy kính mến nhiệt tình hướng dẫn tơi q trình học tập thực luận văn Tôi thực học tập nhiều điều bổ ích trình làm việc với thầy Vợ gái yêu quý động viên dành cho tình cảm u q Thành phố Hồ Chí Minh, ngày 24 tháng năm 2011 Học viên thực Phan Công Thành report fake MỤC LỤC Phần A MỞ ĐẦU Lí chọn đề tài Mục đích đề tài nghiên cứu Đối tượng phạm vi nghiên cứu Ý nghĩa khoa học thực tiễn đề tài nghiên cứu Phương pháp nghiên cứu Cấu trúc luận văn Phần B NỘI DUNG LUẬN VĂN Chương CÁC KHÁI NIỆM CƠ BẢN 1.1 Mơ hình plasma thành phần 1.2 Khái niệm chắn hàm phân bố xuyên tâm 10 1.2.1.Thế chắn 10 1.2.2 Hàm phân bố xuyên tâm 10 1.2.3 Liên hệ chắn hàm phân bố xuyên tâm Định lí Widom 15 Chương THẾ MÀN CHẮN TRONG MÔI TRƯỜNG 17 2.1 Lịch sử 17 2.1.1 Mô MC cho plasma 17 2.1.2 Phương pháp tham số hóa cực trị địa phương [1] 17 2.1.3 Phương pháp tối thiểu hóa bình phương độ lệch dùng phần mềm tính tốn 19 2.2 Tìm chắn nhờ phương pháp tối thiểu hóa bình phương độ lệch 20 2.2.1 Tìm lại chắn cho OCP ứng với số giá trị tham số tương liên 20 2.2.1.1 Quy trình xác định biểu thức chắn 20 2.2.1.2 Bậc chắn 22 2.2.1.3 Đánh giá sai số biểu thức hàm phân bố xuyên tâm thu 22 2.2.2 Ngoại suy chắn với Γ >160 30 2.2.2.1 Lịch sử 30 2.2.2.2 Ngoại suy chắn cơng trình 32 2.2.2.3 Đánh giá độ tin cậy công thức 2.11 số liệu bảng 2.8 xác định chắn 42 2.2.2.3.1 Đánh giá sai số hàm phân bố xuyên tâm 42 2.2.2.3.2 Đánh giá định tính mối tương quan chắn ngoại suy liệu mô Monte Carlo 45 2.2.2.3.3 Đánh giá định tính mối tương quan cực đại thứ hàm phân bố xuyên tâm 47 2.3 Tính tuyến tính chắn 49 TÓM TẮT CHƯƠNG 59 Chương NHIỆT ĐỘNG LỰC HỌC CỦA PLASMA 63 3.1 Các đại lượng nhiệt động học plasma OCP [4] 63 3.2 Thế chắn mạng tinh thể lập phương tâm khối 65 3.3 Mối tương quan chắn plasma lưu chất chắn mạng tinh thể lập phương tâm khối 67 3.4 Điểm kết tinh plasma OCP 69 3.5 Nội nhiệt dung riêng plasma OCP 73 3.5.1 Nội plasma 73 3.5.2 Nhiệt dung đẳng tích plasma 74 3.5.3 Phần dư nội cho ion u 75 KẾT LUẬN 76 TÀI LIỆU THAM KHẢO 77 Phần A MỞ ĐẦU Lí chọn đề tài Những nghiên cứu vật lý vật chất trạng thái plasma, tạo kiến thức tảng cho việc xây dựng công nghệ tổng hợp hạt nhân có điều khiển, nghiên cứu có tính thời cao giới nhằm hướng đến làm chủ nguồn lượng từ phản ứng tổng hợp hạt nhân Đồng thời, kiến thức plasma quan trọng vật lý thiên văn nghiên cứu lùn trắng nơtron Ở Việt Nam, việc nghiên cứu phổ biến kiến thức vật lý plasma chưa nhiều Do việc tìm hiểu nghiên cứu trạng thái thứ tư vật chất khơng phải cơng việc ý nghĩa Một cản trở lớn phần lớn vật chất vũ trụ tồn trạng thái iơn hóa hồn tồn plasma việc tạo phản ứng tổng hợp hạt nhân plasma phịng thí nghiệm Trái đất điều dễ dàng Do vậy, nhiều thông số plasma thu nhờ tính tốn mơ máy tính, thơng qua việc xử lý số liệu thực nghiệm ảo đó, nhà nghiên cứu rút quy luật chi phối loại vật chất Một điều may mắn cho nghiên cứu từ Việt Nam tiếp cận số liệu mơ kết nghiên cứu cập nhật giới Về plasma OCP liên kết mạnh trạng thái cân lưu chất - tinh thể, nghiên cứu trước trình kết tinh bắt đầu xảy hệ số tương liên Γ =172 , nhiên kết gần cho thấy hệ số tương liên tương ứng phải 175 Do nguyên nhân đó, tồn nhu cầu xác định đặc trưng nhiệt động lực plasma trạng thái chuyển pha ứng với Γ =175 Đồng thời, giá trị chắn hàm phân bố xuyên tâm điểm chuyển pha, số tính chất hàm nhiệt động lực cần xác định.Vì vậy, đề tài ”Nhiệt động lực học plasma trạng thái kết tinh” lựa chọn để nghiên cứu khn khổ luận văn Mục đích đề tài nghiên cứu Đề tài nghiên cứu nhằm đưa biểu thức giải tích cho chắn plasma trạng thái kết tinh ứng với hệ số tương liên Γ =175 đề xuất hệ thức biểu thị đặc trưng nhiệt động lực hệ plasma Đối tượng phạm vi nghiên cứu 3.1 Đối tượng nghiên cứu - Plasma thành phần (OCP - One Component Plasma) cổ điển trạng thái cân lưu chất - tinh thể 3.2 Phạm vi nghiên cứu: - Thế chắn plasma OCP trạng thái cân lưu chất - tinh thể - Một số đặc trưng nhiệt động lực hệ plasma trạng thái kết tinh Ý nghĩa khoa học thực tiễn đề tài nghiên cứu 4.1 Ý nghĩa khoa học - Đề tài dựa liệu mô Monte Carlo – MC – hàm phân bố xuyên tâm số plasma OCP đậm đặc để tìm chắn plasma đậm đặc khác, đặc biệt khai triển phương pháp ngoại suy chắn plasma có tham số tương liên lớn nhằm mục đích so sánh với giá trị đề nghị cơng trình có [10, 5] - Trình bày rõ số khái niệm nhiệt động lực học plasma nhờ biểu thức giải tích cơng trình [9] 4.2 Ý nghĩa thực tiễn Đề tài làm tài liệu tham khảo cho sinh viên chuyên ngành vật lý có học Vật lý Thống kê, để có hội đào sâu kiến thức liên quan đến tương tác hệ nhiều hạt, nhằm cung cấp tư liệu nghiên cứu hữu ích cho nghiên cứu trạng thái plasma đậm đặc vật chất Với thân người thực đề tài, tơi có hội học tập phương pháp nghiên cứu khoa học, học cách xử lí liệu, tìm hiểu thêm chương trình Mathematica dùng để xử lý liệu Monte Carlo Phương pháp nghiên cứu Phương pháp lý thuyết: - Nghiên cứu lý thuyết chắn định lí Widom để xây dựng biểu thức chắn với độ xác cao - Dùng phần mềm Mathematica xử lý số liệu mô Monte Carlo Cấu trúc luận văn Chương Các khái niệm Trình bày mơ hình plasma thành phần, khái niệm chắn hàm phân bố xuyên tâm plasma Thế chắn tinh thể lập phương tâm khối Chương Thế chắn môi trường plasma đậm đặc Dựa liệu mô Monte Carlo – MC – hàm phân bố xuyên tâm số plasma OCP đậm đặc để tìm biểu thức chắn plasma với độ xác cao Trên sở xây dựng cơng thức xác định chắn cho plasma lưu chất có tham số tương liên Γ > 20 Từ đó, ta suy chắn plasma trạng thái kết tinh Chương Các đặc trưng nhiệt động lực học Dựa vào biểu thức fit hàm lượng tự cơng trình [9], tìm đặc trưng nhiệt động lực học plasma OCP nội năng, áp suất, nhiệt dung đẳng áp Kết luận chung luận văn Phần B NỘI DUNG LUẬN VĂN Chương CÁC KHÁI NIỆM CƠ BẢN Những khái niệm vật lý plasma trình bày chương chủ yếu dựa theo từ tài liệu [4] 1.1 Mơ hình plasma thành phần Plasma thành phần cổ điển trạng thái cân gồm ion điện tích dương chuyển động biển electron coi phân bố liên tục có số hạt ion N, thể tích V bình chứa nhiệt độ T xác định Các plasma thường phân loại làm plasma liên kết mạnh plasma liên kết yếu tùy theo tỷ số tương tác Coulomb với lượng chuyển động nhiệt Việc phân loại định lượng cách đưa vào khái niệm khối cầu ion tham số tương liên sau: Xét ion plasma, có điện tích Ze, tạo khối cầu chứa điện tích âm vừa đủ để trung hịa điện tích dương ion Như vậy, gọi ρ = N mật độ ion khối plasma xét, a bán kính khối cầu ion V  4πρ  tính a =     −1/3 Trong vật lý plasma, lượng chuyển động nhiệt quy ước kT Tham số tương liên plasma định nghĩa Γ = ( Ze) tỉ akT số lượng tương tác Coulomb trung bình lượng chuyển động nhiệt trung bình Khi Γ >> , tức lượng Coulomb lớn chuyển động nhiệt, trạng thái plasma gần với trạng thái rắn Ngược lại, Γ