Nghiên cứu nhiệt độ nóng chảy của các kim loại dưới áp suất cao

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Nguyễn Ngọc Hà NGHIÊN CỨU NHIỆT ĐỘ NÓNG CHẢY CỦA CÁC KIM LOẠI DƯỚI ÁP SUẤT CAO LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội – 2014 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - Nguyễn Ngọc Hà NGHIÊN CỨU NHIỆT ĐỘ NÓNG CHẢY CỦA CÁC KIM LOẠI DƯỚI ÁP SUẤT CAO Chuyên ngành : Vật lý lý thuyết vật lý toán Mã số : 60440103 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC : TS HỒ KHẮC HIẾU Hà Nội – 2014 Luận văn thạc sỹ LỜI CẢM ƠN Để hoàn thành luận văn em xin bày tỏ lòng biết ơn chân thành, sâu sắc đến TS Hồ Khắc Hiếu tận tâm hướng dẫn, giúp đỡ em suốt trình làm luận văn Em xin cảm ơn thầy cô giáo khoa Vật lí Trường Đại Học Khoa học Tự nhiên nói chung thầy giáo Bộ mơn Vật lý lý thuyết nói riêng hết lịng dạy bảo em suốt hai năm học tập rèn luyện vừa qua Tôi xin cảm ơn thầy ban giám hiệu Trường THPT Chúc ĐộngChương Mỹ- Hà Nội, tồn thể thầy giáo, anh chị bạn đồng nghiệp giúp đỡ, tạo điều kiện thuận lợi để tơi học tập thời gian năm qua hoàn thành luận văn Xin cảm ơn gia đình, anh chị bạn giúp đỡ, ủng hộ, động viên, chia sẻ khó khăn tơi suốt thời gian học tập làm luận văn Tác giả Nguyễn Ngọc Hà Khoa Vật Lý Luận văn thạc sỹ MỤC LỤC MỞ ĐẦU .1 Chương - HIỆN TƯỢNG NÓNG CHẢY VÀ CÁC PHƯƠNG PHÁP NGHIÊN CỨU 1.1 Tổng quan tượng nóng chảy 1.2 Các phương pháp nghiên cứu nóng chảy 1.2.1 Phương pháp ô mạng đế kim cương (Diamond anvil cell - DAC) 1.2.2 Phương pháp thống kê mômen 10 1.2.2.1 Nhiệt độ nóng chảy kim loại P = 11 1.2.2.2 Nhiệt độ nóng chảy kim loại áp suất cao P 12 1.2.3 Phương pháp mô 14 Chương - NGHIÊN CỨU NHIỆT ĐỘ NÓNG CHẢY CỦA VẬT LIỆU DƯỚI ÁP SUẤT CAO 18 2.1 Giới hạn Lindemann nóng chảy hệ số Grüneisen 18 2.1.1 Giới hạn Lindemann nóng chảy 18 2.1.1.1 Mơ hình Debye nghiên cứu tính chất nhiệt động tinh thể 18 2.1.1.2 Hệ số Debye-Waller 20 2.1.1.3 Giả thuyết Lindemann nóng chảy .24 2.1.2 Hệ số Grüneisen áp suất cao 26 2.2 Biểu thức nhiệt độ nóng chảy vật liệu áp suất cao 27 Chương - TÍNH TỐN SỐ VÀ THẢO LUẬN 31 Khoa Vật lý Luận văn thạc sỹ 3.1 Nhiệt độ nóng chảy kim loại Cu, Au Ag áp suất cao 31 3.2 Nhiệt độ nóng chảy kim loại Fe áp suất cao 41 KẾT LUẬN .52 DANH MỤC CÁC CƠNG TRÌNH ĐÃ ĐƯỢC CƠNG BỐ 53 TÀI LIỆU THAM KHẢO 54 Khoa Vật lý Luận văn thạc sỹ DANH MỤC BẢNG BIỂU Tên bảng Nội dung Trang Bảng 3.1 Giá trị áp suất chuyển pha cấu trúc nhiệt độ nóng chảy 31 áp suất P =0 kim loại Bảng 3.2 Các thông số làm khớp q  nhóm Graf theo cơng 32 thức (16) cho hai kim loại Au Cu Bảng 3.3 Giá trị làm khớp K0 K 0 phương pháp bình phương 35 tối thiểu kim loại Ag, Au Cu Bảng 3.4 Sự phụ thuộc thể tích đại lượng vật lý ϵ-Fe nhiệt độ T = 300 K Khoa Vật lý 43 Luận văn thạc sỹ DANH MỤC HÌNH VẼ Tên hình Nội dung Trang Hình 1.1 Phương pháp mạng đế kim cương Hình 3.1 Đồ thị phụ thuộc hệ số nén V/V0 nhiệt độ nóng chảy 34 kim loại Hình 3.2 Đồ thị phụ thuộc áp suất nhiệt độ nóng chảy kim loại Au 37 Hình 3.3 Đồ thị phụ thuộc áp suất nhiệt độ nóng chảy kim loại Au 39 Hình 3.4 Đồ thị phụ thuộc áp suất nhiệt độ nóng chảy Ag 40 Hình 3.5 Giản đồ pha Fe 42 Hình 3.6 Hệ số Grüneisen thực nghiệm lý thuyết kim loại Fe 43 Hình 3.7 Đường cong P-V nhiệt độ T = 300K kim loại Fe 45 Hình 3.8 Đồ thị phụ thuộc hệ số nén V/V0 nhiệt độ nóng chảy kim 46 loại Fe Hình 3.9 Đồ thị phụ thuộc áp suấtcủa nhiệt độ nóng chảy kim loại ϵ-Fe đến áp suất 350 GPa Khoa Vật lý 48 Luận văn thạc sỹ MỞ ĐẦU I Lý chọn đề tài Trong năm gần khoa học công nghệ giới phát triển mạnh, đặc biệt ngành khoa học vật liệu Vì việc nghiên cứu hiểu biết tính chất cơ, nhiệt, điện, quang,… vật liệu yêu cầu cấp thiết, thu hút quan tâm nhà khoa học nói chung nhà vật lý nói riêng Đối với tính chất nhiệt động vật rắn, nhiệt độ nóng chảy áp suất cao hướng nghiên cứu liên ngành ý quan trọng lĩnh vực khoa học vật liệu, địa vật lý, vật lý địa cầu vật lý thiên văn Nói chung khả dẫn điện, dẫn nhiệt cấu trúc điện tử kim loại khác nhau, nhiệt độ nóng chảy chúng khác Ngồi ra, nhiệt độ nóng chảy kim loại chịu ảnh hưởng lớn áp suất bên ngồi Cho đến nay, có nhiều phương pháp nghiên cứu sử dụng để nghiên cứu ảnh hưởng áp suất đến nhiệt độ nóng chảy kim loại phương pháp thực nghiệm (ô mạng đế kim cương), phương pháp phiếm hàm mật độ, phương pháp bán thực nghiệm Tuy nhiên phương pháp nhiều hạn chế như: Các biểu thức toán học cồng kềnh, phức tạp, khó khăn đưa số liệu thực nghiệm, sai số lớn Vì việc nghiên cứu nhiệt độ nóng chảy kim loại áp suất đặc biệt áp suất cao vấn đề thời nhà nghiên cứu thực nghiệm Vì lý đó, chúng tơi chọn tốn “Nghiên cứu nhiệt độ nóng chảy kim loại áp suất cao” làm đề tài luận văn II Đối tượng nghiên cứu Đối tượng nghiên cứu luận văn kim loại chuyển tiếp Cu, Ag, Au, Fe Chúng xác định nhiệt độ nóng chảy áp suất cao Khoa Vật lý Luận văn thạc sỹ kim loại này, từ vẽ đường biểu diễn phụ thuộc nhiệt độ theo áp suất phụ thuộc vào hệ số nén V/V0 nhiệt độ III Mục đích nhiệm vụ nghiên cứu Tiên đốn nhiệt độ nóng chảy kim loại áp suất cao từ biểu diễn đường cong nóng chảy số kim loại điển Cu, Ag, Au, Fe nhiều kim loại khác So sánh giá trị tính tốn lý thuyết số liệu thực nghiệm để kiểm tra lý thuyết IV Phương pháp nghiên cứu Phương pháp nghiên cứu luận văn phương pháp bán thực nghiệm, dựa điều kiện nóng chảy Lindemann phụ thuộc áp suất hệ số Grüneisen để xây dựng biểu thức giải tích nhiệt độ nóng chảy Sử dụng số liệu thực nghiệm biểu thức giải tích lý thuyết thu để xác định đường cong nóng chảy kim loại điển Cu, Ag, Au, Fe V Đóng góp đề tài Đối tượng nghiên cứu luận văn loại vật liệu sử dụng rộng rãi Các kết thu góp phần hồn thiện phát triển ứng dụng điều kiện Lindemann nóng chảy hệ số Grüneisen việc nghiên cứu tính chất vật liệu đặc biệt nhiệt độ nóng chảy áp suất cao Luận văn gợi mở phương pháp để nghiên cứu loại vật liệu khác VI Cấu trúc luận văn Luận văn cấu trúc gồm phần mở đầu, ba chương, phần kết luận tài liệu tham khảo Chương Hiện tượng nóng chảy phương pháp nghiên cứu Khoa Vật lý Luận văn thạc sỹ Nội dung chương trình bày vắn tắt kiến thức tượng nóng chảy kim loại phương pháp để nghiên cứu Mỗi phương pháp có đặc điểm riêng để nghiên cứu tượng nóng chảy, số phương pháp có phương pháp hay dùng phương pháp thống kê momen, nhiên luận văn áp dụng phương pháp khác để nghiên cứu nhiệt độ nóng chảy kim loại áp suất cao giới hạn Lindemann hệ số Grüneisen Chương Nghiên cứu nhiệt độ nóng chảy vật liệu áp suất cao Phần đầu chương tơi trình bày nét Mơ hình Debye nghiên cứu tính chất nhiệt động vật liệu giới hạn Lindemann nghiên cứu tượng nóng chảy kim loại Từ đó, kết hợp với biểu thức phụ thuộc áp suất hệ số Grüneisen chúng tơi thiết lập biểu thức giải tích nhiệt độ nóng chảy Tm hàm thể tích.Dựa phương trình trạng thái Vinet, chúng tơi nghiên cứu ảnh hướng áp suất đến nhiệt độ nóng chảy Tm Chương Tính tốn số thảo luận Trong chương này, sở biểu thức giải tích nhiệt độ nóng chảy Tm thu chương phương trình trạng thái Vinet tinh thể, chúng tơi thực tính tốn số thảo luận cho kim loại chuyển tiếp đồng (Cu), vàng (Au), bạc (Ag) sắt (Fe) Ảnh hưởng thể tích áp suất đến nhiệt độ nóng chảy Tm kim loại nghiên cứu đến giá trị hệ số nén V V0  0,5 đến áp suất tương ứng Từ đồ thị phụ thuộc áp suất nhiệt độ kim loại Cu, Ag, Au, Fe đưa kết thảo luận kết Khoa Vật lý Luận văn thạc sỹ P (GPa) G (Thực (Thực nghiệm) nghiệm) 0.877 P (GPa) G (Lý thuyết) (Lý thuyết) 1.71 1.74 29.96 1.59 28.91 1.58 0.862 35.56 1.57 34.41 1.56 0.847 41.75 1.55 40.51 1.54 0.832 48.62 1.53 47.28 1.52 0.817 56.22 1.51 54.77 1.50 0.802 64.64 1.49 63.09 1.48 0.788 73.99 1.47 71.66 1.47 0.773 84.38 1.45 81.82 1.45 0.758 95.92 1.43 93.10 1.42 0.743 108.76 1.41 105.63 1.40 0.728 123.05 1.39 119.55 1.38 0.713 139.00 1.37 135.04 1.36 0.698 156.80 1.35 152.28 1.34 0.684 176.71 1.32 170.15 1.32 0.669 198.99 1.30 191.42 1.30 0.654 223.98 1.28 215.17 1.28 0.639 252.06 1.25 241.71 1.26 0.624 283.65 1.23 271.41 1.24 0.609 319.26 1.20 304.69 1.22 0.594 359.50 1.18 342.03 1.19 V V0 Về mối liên hệ P-V hay phương trình trạng thái tinh thể, chúng tơi sử dụng phương trình trạng thái Vinet nhiệt độ T = 300 K Phương trình trạng thái thực nghiệm lý thuyết (phương trình trạng thái Vinet) ϵ-Fe Khoa Vật lý 44 Luận văn thạc sỹ biểu diễn hình 3.7 Kết tính tốn hệ số Grüneisen áp suất P theo thực nghiệm lý thuyết chúng tơi đưa bảng 3.4 Từ hình vẽ 3.7 bảng 3.4 kết luận rằng, phương trình trạng thái Vinet mơ tả tốt mối liên hệ áp suất-thể tích kim loại ϵ-Fe Trong phần tiếp theo, sử dụng hệ số làm khớp   1.74232, q  0.72659 T0  1608 K [12, 22] để xác định ảnh hưởng thể tích đến nhiệt độ nóng chảy Tm ϵ-Fe Tiếp đó, phương trình trạng thái Vinet sử dụng để biểu diễn phụ thuộc áp suất đường cong nóng chảy Tm ϵ-Fe Hình 3.7.Đường cong P-V nhiệt độ T = 300 Kcủa kim loại ϵ-Fe ϵ-Fe pha cấu trúc lục giác xếp chặt (Hexagonal close-packed - HCP) kim loại sắt Pha cho tồn ổn định áp suất cao Trong cơng trình [42], Takahashi Bassett rằng, pha sắt với cấu trúc lập phương tâm khối (Body-centered cubic - BCC) (α-Fe) chuyển pha cấu trúc sang ϵ-Fe áp suất 13GPa nhiệt độ phòng Khi áp suất giảm xuống, pha ϵ-Fe nhanh chóng quay trở lại pha α-Fe Đã có nhiều phương pháp thực nghiệm thực để đánh giá nhiệt độ nóng chảy ϵ-Fe áp suất cao Khoa Vật lý 45 Luận văn thạc sỹ như: phương pháp thực nghiệm nén đột ngột [27], nhiễu xạ tia X nén tĩnh [32, 41], thí nghiệm nén tĩnh với nguồn laser đốt nóng [16],thí nghiệm kết hợp nén tĩnh nén đột ngột [47] Tuy nhiên, nay, chưa có lý thuyết hồn chỉnh tiên đốn xác nhiệt độ nóng chảy ϵ-Fe Trên hình 3.8, chúng tơi biểu diễn đồ thị phụ thuộc hệ số nén V/V0 nhiệt độ nóng chảy kim loại ϵ-Fe Dáng điệu đường cong nóng chảy ϵ-Fe có dạng tăng chậm so với kim loại Cu, Au Ag Cụ thể, hệ số nén V V0  0,5 , nhiệt độ nóng chảy ϵ-Fe khoảng 6800 K, giá trị Cu, Au Ag 10000 K, 13000 K 11000 K Hình 3.8.Đồ thị phụ thuộc hệ số nén V/V0của nhiệt độ nóng chảy kim loại ϵ-Fe Sử dụng mối liên hệ P-V từ phương trình trạng thái Vinet, xác định ảnh hưởng áp suất đến nhiệt độ nóng chảy Tm sắt Do ϵ-Fe tồn bền vững áp suất 13 GPa, vậy, hình 3.9, chúng tơi biểu diễn đường cong nóng chảy phụ thuộc áp suất kim loại ϵ-Fe từ 20 GPa đến 350GPa Kết thực nghiệm phương pháp nén đột ngột Nguyen Khoa Vật lý 46 Luận văn thạc sỹ Holmes [36] (ký hiệu x), nhóm Ahrens [12] (ký hiệu +), Brown Mc Queen [18] (ký hiệu *) đưa vào để so sánh Ngồi ra, chúng tơi biểu diễn số liệu thực nghiệm đo nhiễu xạ tia X dựa phương pháp nén tĩnh nhóm Komabayashi Fei [29] (ký hiệu hình tam giác), Ma cộng [32] (ký hiệu hình vng), nhóm Shen [41] (ký hiệu hình kim cương) Có thể thấy từ hình 3.9, đến khoảng áp suất 200 GPa, kết tính tốn lý thuyết nhóm phù hợp tốt với kết thí nghiệm nhiễu xạ tia X nhóm Komabayashi Fei [29], Ma cộng [32], Shen cộng [41] Tuy vậy, sử dụng phương pháp thí nghiệm khác nhau, giá trị nhiệt độ nóng chảy Tm thu khơng thống Lấy ví dụ, áp suất 135 GPa, giá trị lý thuyết nhiệt độ nóng chảy Tm(135GPa)=3649.5 K Trong đó, kết thực nghiệm dựa phương pháp nén đột ngột thực nhóm Ahrens cộng [12] Tm (135GPa)=3400±200 K Dựa chuyển động đối lưu khí nghiệm nén tĩnh với nguồn laser đốt nóng, Boehler xác định giá trị Tm (135GPa)=3200±100 K [16]; ngược lại, nhóm Williams cộng [47] lại đưa số liệu Tm (135GPa)= 4800±200 K kết hợp thí nghiệm nén tĩnh nén đột ngột Trong trường hợp này, giá trị tính tốn lý thuyết có ý nghĩa tiên đốn, giúp thực nghiệm ước lượng độ lớn nhiệt độ nóng chảy áp suất cần nghiên cứu Khoa Vật lý 47 Luận văn thạc sỹ Hình 3.9.Đồ thị phụ thuộc áp suấtcủa nhiệt độ nóng chảy kim loại ϵ-Fe đến áp suất 350 GPa Kết tính tốn lý thuyết nhóm chúng tơi so sánh với giá trị thực nghiệm Nguyen Holmes [36] (dấu x), Ahrens cộng [12] (dấu +), Brown McQueen [18] (dấu *), Komabayashi Fei [29] (dấu◄), Ma cộng [32] (dấu □), Shen cộng [41] (dấu ◊), Boehler [16] (hình sáu cạnh) Ngoài vùng áp suất 200 GPa, độ dốc đường cong nóng chảy lý thuyết có xu hướng giảm dần Nhiệt độ nóng chảy tính tốn nhóm chúng tơi vùng áp suất thấp giá trị thực nghiệm đo nhóm Nguyen Holmes [36], Brown Mc Queen [18]; nhiên lại phù hợp với giá trị tiên đoán Boehler [16] áp suất 330 GPa (được cho áp suất tâm lõi Trái đất) Theo cơng trình cơng bố tạp chí Nature năm 1993, Boehler cho rằng, nhiệt độ nóng chảy ϵ-Fe áp suất 330 GPa vào khoảng 5100 K Tính tốn chúng tơi theo cơng thức lý thuyết (2.34) xác định Tm (330GPa)=5080 K, sai số so với kết Boehler khoảng 0,4% Để kết luận, muốn đưa số nhận xét cách thức tiếp cận tốn nóng chảy dựa điều kiện nóng chảy Lindemann Trong cơng trình [28], Japel cộng vai trò quan trọng lớp vỏ electron đến Khoa Vật lý 48 Luận văn thạc sỹ nhiệt độ nóng chảy kim loại chuyển tiếp Một số kim loại chuyển tiếp có nhiệt độ nóng chảy thấp đáng ngạc nhiên Ni, Mo,… giải thích kể đến cấu hình electron chúng Trong luận văn này, kim loại Ag, Au Cu có cấu hình electron với electron lớp vỏ d lấp đầy (cấu hình electron Ag, Au Cu tương ứng 4d105s1, 5d106s1và 3d104s1) Điều cho phép tiên đốn đường cong nóng chảy ba kim loại có dạng hình học Dựa việc quan sát hình 3.1 (a), hình 3.1 (b) hình 3.1 (c) xác thực kết luận Tuy nhiên, để mơ ta cách xác đường cong nóng chảy áp suất cao kim loại chuyển tiếp, lý thuyết cần xây dựng theo hướng ý đến tính chất điện tử vật liệu.Trong luận văn này, sử dụng giả thiết Lindemann nóng chảy, chúng tơi bỏ qua việc xem xét đến cấu trúc điện tử kim loại Au, Ag, Cu Fe Dầu vậy, điều kiện Lindemann sử dụng để tiên đốn độ lớn nhiệt độ nóng chảy vật liệu áp suất áp suất siêu cao thực nghiệm chưa tiếp cận Lấy ví dụ, quan sát hình 3.3 (b) thấy, giá trị tính tốn lý thuyết nhiệt độ nóng chảy Au 13 GPa khoảng 1800 K Giá trị trùng khớp với giá trị ngoại suy từ kết thực nghiệm tác giả Errandonea [22] Ngồi ra, chúng tơi đưa ý khác là, hàm nhiệt độ nóng chảy Tm khơng thực biến thiên tuyến tính theo áp suất P, đặc biệt vùng áp suất cao Nhận xét dễ dàng kiểm nghiệm thơng qua hình vẽ 3.1: Đường cong nóng chảy kim loại Ag, Au Cu có có dạng hàm phi tuyến áp suất P hệ số nén V V0  0,8 (tương ứng với áp suất P  43 GPa kim loại Ag, P  70 GPa kim loại Au P  50 GPa kim loại Cu [21] Như vậy, luận văn này, đưa cách tiếp cận đơn giản lý thuyết để nghiên cứu tốn nóng chảy áp suất cao kim loại Ag, Au, Cu ϵ-Fe dựa điều kiện Lindemann nóng chảy biểu thức phụ thuộc áp suất hệ số Grüneisen Chúng xây dựng thành công biểu Khoa Vật lý 49 Luận văn thạc sỹ thức giải tích phụ thuộc hệ số nén (tức phụ thuộc thể tích) áp suất nhiệt độ nóng chảy Tm Kết tính tốn số chúng tơi thực đến hệ số nén V V0  0,5 đến áp suất tương ứng với hệ số nén (460 GPa kim loại Ag, 770 GPa kim loại Au 500GPa kim loại Cu) Đối với ϵ-Fe, chúng tơi nghiên cứu tốn đến áp suất 350 GPa (đây giá trị áp suất tiên đoán tồn tâm lõi Trái đất) Trên sở việc so sánh kết tính tốn lý thuyết nhóm chúng tơi với số liệu thực nghiệm, lý thuyết mơ trước chúng tơi kết luận rằng: Cách thức tiếp cận toán nóng chảy dựa giả thuyết Lindemann sử dụng nghiên cứu nhiệt độ nóng chảy kim loại Ag Au tương ứng đến áp suất khoảng 12GPa 6GPa Ở vùng áp suất cao hơn, cách thức tiếp cận Lindemann chủ yếu giúp nghiên cứu định tính tốn nóng chảy, giúp ước lượng độ lớn nhiệt độ nóng chảy áp suất cao Ngoài ra, phương pháp sử dụng để kiểm nghiệm kết thực nghiệm (sử dụng phương pháp đa đế hay ô mạng đế kim cương) kết tính tốn lý thuyết khác tương lai Ví dụ, trường hợp ϵ-Fe, kết lý thuyết phù hợp đáng ngạc nhiên với tiên đoán thực nghiệm áp suất siêu cao P = 330 GPa Cách thức tiếp cận tốn trình bày luận văn sử dụng nghiên cứu ảnh hưởng áp suất đến nhiệt độ nóng chảy kim loại khác Pt, Pd, Riêng kim loại đồng, giá trị lý thuyết TmB (sử dụng thông số làm khớp  , q hệ số Grüneisen xác định gần mơđun nén khối B) sử dụng để tiên tốn nhiệt độ nóng chảy Cu đến áp suất siêu cao (hàng trăm GPa) Chúng cho rằng, để nghiên cứu xác phụ thuộc áp suất nhiệt độ nóng chảy kim loại chuyển tiếp, cần ý đến phát triển lý thuyết có ý đến cấu hình electron Khoa Vật lý 50 Luận văn thạc sỹ Tiểu kết chương 3: Như thông qua đường biểu diễn đường cong nóng chảy số kim loại Fe,Cu,Ag… ta thấy rõ phụ thuộc nhiệt độ nóng chảy vào áp suất dựa điều kiện Lindemann hệ số Grüneisen Từ đường biểu diễn ta thấy với giá trị hệ số định khoảng áp suất, nhệt độ nóng chảy có thay đổi kim loại đường biểu diễn có khác Khoa Vật lý 51 Luận văn thạc sỹ KẾT LUẬN Dựa điều kiện Lindemann nóng chảy phụ thuộc áp suất hệ số Grüneisen nghiên cứu nhiệt độ nóng chảy kim loại chuyển tiếp vùng áp suất cao (hàng trăm GPa) Các kết luận văn bao gồm: Xây dựng biểu thức giải tích tường minh nhiệt độ nóng chảy vật liệu phụ thuộc áp suất thể tích dựa điều kiện Lindemann nóng chảy hệ số Grüneisen Vẽ đường cong nóng chảy phụ thuộc hệ số nén V / V0 phụ thuộc áp suất vật liệu Cu, Ag, Au Fe Kết lý thuyết sai khác so với thực nghiệm không 10% vùng áp suất cao Phương pháp lý thuyết mở rộng nghiên cứu nhiệt độ nóng chảy bán dẫn, siêu mạng áp suất khác Ngoài ra, phương pháp sử dụng để kiểm nghiệm kết thực nghiệm (sử dụng phương pháp đa đế hay ô mạng đế kim cương) kết tính tốn lý thuyết khác tương lai Khoa Vật lý 52 Luận văn thạc sỹ DANH MỤC CÁC CƠNG TRÌNH ĐÃ ĐƯỢC CƠNG BỐ Hồ Khắc Hiếu, Nguyễn Ngọc Hà (2013), “Nghiên cứu nhiệt độ nóng chảy kim loại vàng áp suất cao”, Tạp chí Khoa học Cơng nghệ, Trường Đại học Duy Tân, 3, tr.1-6 Ho Khac Hieu, Nguyen Ngoc Ha (2013), “High pressure melting curves of silver, gold and copper”, AIP Advances, 3, 112125, pp 1-9 Khoa Vật lý 53 Luận văn thạc sỹ TÀI LIỆU THAM KHẢO Phần tiếng Việt [1] Lê Thị Thanh Bình (2010), Nghiên cứu ảnh hưởng vacancy lên nhiệt độ nóng chảy kim loại, Luận văn thạc sĩ khoa học Vật lý, Đại học Sư phạm, Hà Nội [2] Vũ Thị Ngọc Chính (2012), Nghiên cứu nhiệt độ nóng chảy cuả hợp kim đôi phương pháp thống kê momen, Luận văn thạc sĩ khoa học Tự nhiên, Đại học Khoa học Tự nhiên, Hà Nội [3] Nguyễn Thanh Hải, Vũ Văn Hùng (1997), “Nghiên cứu tính chất nhiệt động kim loại áp suất khác không”, Tuyển tập báo cáo hội nghị Vật lý lý thuyết lần thứ XXII , Đồ Sơn, pp 204 [4] Nguyễn Xuân Hãn (1998), Cơ sở lý thuyết trường lượng tử, NXB Đại học Quốc gia, Hà Nội [5] Nguyễn Thị Hằng (2010), Nghiên cứu phương trình trạng thái đường cong nóng chảy kim loại, Luận văn thạc sỹ khoa học vật lý, Trường Đại học sư phạm, Hà Nội [6] Hồ Khắc Hiếu, Nguyễn Ngọc Hà (2013), “Nghiên cứu nhiệt độ nóng chảy kim loại vàng áp suất cao”, Tạp chí Khoa học Cơng nghệ, Trường Đại học Duy Tân, 3, tr.1-6 [7] Nguyễn Văn Hùng (2000), Lý thuyết chất rắn, NXB Đại học Quốc Gia, Hà Nội, tr 189-190 [8] Vũ Văn Hùng (2009), Phương pháp thống kê momen nghiên cứu tính chất nhiệt động đàn hồi tinh thể, NXB Đại học Sư phạm, Hà Nội Khoa Vật lý 54 Luận văn thạc sỹ Phần tiếng Anh [9] Akella J., Kennedy G C (1971), “Melting of gold, silver, and copper-proposal for a new high-pressure calibration scale”, Journal of Geophysical Research, 76, 4969 [10] Anderson O L., Dubrovinsky L., Saxena S K., LeBihan T (2001), “Experimental vibrational Grüneisen ratio values for ϵ-iron up to 330 GPa at 300 K”, Geophysical Research Letters, 28, 399 [11].Anzellini S., Dewaele A., Mezouar M., Loubeyre P., Morard G (2013), “Melting of Iron at Earth’s Inner Core Boundary Based on Fast X-ray Diffraction”, Science, 340, pp 464-466 [12] Ahrens T J., Holland K G., Chen G Q (2002), “Phase diagram of iron, revised-core temperatures”, Geophysical Research Letters, 29(7), pp 54-1– 54-4 [13] Arafin S., Singh R., George A (2013), “Melting of metals under pressure”, Physica B: Condensed Matter, 41, pp 40-44 [14].Bazarov P I., Kotenok V V (1972), “Bogolyubov method in the theory of crystals with three- particle interaction”, Teoreticheskaya Matematicheskaya Fizika, 10(2), pp 275-282 [15].Belonoshko B A., Ahuja R., Eriksson O., Johans-son B (2000), “Quasi ab initio molecular dynamic study of Cumelting”, Physical Review B, 61, pp 3838-384 [16].Boehler R (2000) "High-pressure experiments and the phase diagram of lower mantle and core materials",Review of Geophysics, 38 (2), pp.221–245 [17] Brand H., Dobson P D., Vocadlo L., Wood G I (2006), “Melting curve of copper measured to 16 GPa using a multi-anvil press”, High Pressure Research, 26, pp 185-191 [18] Brown J M., McQueen R G (1986), “Phase transitions, Grüneisen parameter, and elasticity for shocked iron between 77 GPa and 400 GPa”, Journal of Khoa Vật lý 55 Luận văn thạc sỹ Geophysical Research: Solid Earth,91(B7), pp 7485–7494 [19] Burakowsky L., PrestonD L., and SilbarR R (2000), “Analysis of dislocation mechanism for melting of elements: Pressure dependence”, Journal of Applied Physics, 88, 6294 [20] Dewaele A., Torrent M., Loubeyre P., and Mezouar M (2008), “Compression curves of transition metals in the Mbar range: Experiments and projector augmented-wave calculations”, Physical Review B, 78, 104102 [21] Dewaele A., Loubeyre P., and Mezouar M.(2004), “Equations of state of six metals above 94GPa”, Physical Review B, 70, 094112 [22] Errandonea D (2010), “The melting curve of ten metal up to 12 GPa and 1600K”, Journal of Applied Physics, 108(3), pp 033517-033517 [23] Graf J M., Greeff W C., Boettger C J (2004), “High-pressure Debye-Waller and Grüneisen parameters of gold and copper”, AIP Conference Proceedings, 706 , pp 65-68 [24] Girifalco L A (2000), Statistical Mechanics of Solids, Oxford University Press, USA [25] Grüneisen E (1912),“Theorie des festen Zustandes einatomiger Elemente”, An- nalen der Physik, 39, pp 257–306 [26] H K Hieu, N N Ha (2013), “High pressure melting curves of silver, gold and copper”, AIP Advances, 3, 112125, pp 1-9 [27].Huang D., Liu X., Su L., Shao C., Jia R., and Hong S (2007), “Measuring Grüneisen parameter of iron and copper by an improved high pressure-jump method”, Journal of Physics D: Applied Physics, 40, 5327 [28] Japel S., Schwager B., Boehler R., Ros M (2005), “Melting of copper and nickel at high pressure: The role of d electrons”, Physical Review Letters, 95, 167801 [29] Komabayashi T., Fei Y W (2010), “Internally consistent thermodynamic database for iron to the Earth's core conditions”, Journal of Geophysical Research: Solid Earth, 115(B3), 12 Khoa Vật lý 56 Luận văn thạc sỹ [30] Kumari M., Dass N (1990), “An equation of state applied to 50 solids”, Journal of Physics:Condensed Matter, 2(39), pp 7891-7895 [31] Lindemann F (1910), “The calculation of molecular vibration frequencies”, Physica B: Condensed Matter, 419, pp 40-44 [32] Ma Y Z., Somayazulu M., Shen G., Mao H K, Shu J F., Hemley R J (2004), “In situ X-ray diffraction studies of iron to Earth-core conditions”, Physics of the Earth and Planetary Interiors, 143–144, pp 455–467 [33] Metropolis N., Ulam S (1949), "The Monte Carlo Method", Journal of the American Statistical Association, 44(247), pp 335-341 [34] Mirwald P., and Kennedy G C.(1979), “The melting curve of gold, silver, and copper to 60 Kbar pressure: A reinvestigation”, Journal of Geophysical Research: Solid Earth, 84, 6750 [35] N Tang and V V Hung (1990), “Investigation of the Thermodynamic Properties of Anharmonic Crystals by Momentum Method IV.The Limiting of Absolute Stability and the Melting Temperature of Crystals”, Physica Status Solidi B, 162 (2), pp 379–385 [36] Nguyen J H., Neil C., Holmes N C (2004), “Melting of iron at the physical conditions of the Earth's core”, Nature, 427, pp 339-342 [37] Petitgirard S., Daniel I., Dabin Y., Cardon H., Tucoulou R., Susini J (2009), “A diamond anvil cell for x-ray fluorescence measurements of trace elements in fluids at high pressure and high temperature”, The Review of scientific instruments, 80(3), 033906 [38] Ping Y., Coppari F., Hicks D G., Yaakobi B., Fratanduono D E., Hamel S., Eggert J H., Rygg J R., Smith R F., Swift D C., Braun D G., Boehly T R., and Collins G W (2013), “Solid Iron Compressed Up to 560 GPa”, Physical Review Letters, 111, 065501 [39].Pozzo M., and Alfè D (2013), “Melting curve of face-centered-cubic nickel from first-principles calculations”, Physical Review B, 88, 024111 Khoa Vật lý 57 Luận văn thạc sỹ [40].Saxena S K (2004), “Pressure–volume equation of state for solids”, Journal of Physics and Chemistry of Solids, 65, pp.1561–1563 [41] Shen G., Mao H K., Hemley R J., Duffy T S., Rivers M L (1998), “Melting and crystal structure of iron at high pressures and temperatures”, Geophysical Research Letters, 25, 373 [42].Takahashi T., Bassett W A (1964), “High-pressure polymorph of iron”, Science, 145, pp 483-486 [43] Hung V V., and Hai T N (1997), “Investigation of the Melting Temperature of Metal at Various Pressures”, Journal of the Physical Society of Japan, 66, pp 3499-3501 [44].Vinet P., Ferrante J., Rose J., and Smith J.(1987), “Compressibility of solids”, Journal of Geophysical Research: Solid Earth, 92, 9319 [45] Vocadlo L., Alfe D., PriceG D., and Gillan M J (2004), “Ab initio melting curve of copper by the phase coexistence approach”, The Journal of Chemical Physics, 120, 2872 [46].Wang Y., Ahuja R., Johansson B (2001), “ Melting of iron and other metals at earth’s core conditions:A simplified computational approach”, Physical Review B, 65, 014104 [47] Williams Q., Jeanloz R., Bass J., Svendsen B., Ahrens T J (1987), “The melting curve of iron to 250 gigapascals: A constraint on the temperature at Earth's center”, Science, 236, pp 181-182 [48].Wu Y., Wang L., Huang Y., Wang D (2011), “Melting of copper under high pressures by molecular dynamics simulation”, ChemicalPhysics Letters, 515, pp 217-220 Khoa Vật lý 58 ... việc nghiên cứu nhiệt độ nóng chảy kim loại áp suất đặc biệt áp suất cao vấn đề thời nhà nghiên cứu thực nghiệm Vì lý đó, chúng tơi chọn tốn ? ?Nghiên cứu nhiệt độ nóng chảy kim loại áp suất cao? ??... P-V nhiệt độ T = 300K kim loại Fe 45 Hình 3.8 Đồ thị phụ thuộc hệ số nén V/V0 nhiệt độ nóng chảy kim 46 loại Fe Hình 3.9 Đồ thị phụ thuộc áp suấtcủa nhiệt độ nóng chảy kim loại ϵ-Fe đến áp suất. .. thuộc nhiệt độ theo áp suất phụ thuộc vào hệ số nén V/V0 nhiệt độ III Mục đích nhiệm vụ nghiên cứu Tiên đốn nhiệt độ nóng chảy kim loại áp suất cao từ biểu diễn đường cong nóng chảy số kim loại