LỜI CẢM ƠN Sau thời gian nghiên cứu với nỗ lực thân giúp đỡ thầy cô gáo bạn học viên hoàn thành đề tài Trước hết xin bày tỏ lòng biết ơn sâu sắc đến thầy giáo, PGS – TS Phạm Đình Tám tận tình giúp đỡ bảo cung cấp cho kiến thức tảng để hoàn thành luận văn truyền cho niềm say mê khoa học suốt thời gian làm việc thầy Tôi xin bày tỏ lòng biết ơn tới Ban giám hiệu, thầy, cô công tác phòng Sau đại học trường Đại học sư phạm Hà Nội thầy cô giáo trực tiếp giảng dạy, truyền đạt cho kiến thức quý báu chuyên môn kinh nghiệm nghiên cứu khoa học thời gian qua Cuối cùng, xin chân thành gửi lời cảm ơn đến người thân gia đình, bạn bè giúp đỡ, động viên tạo điều kiện cho suốt trình học tập hoàn thiện luận văn Học viên thực hiện: Đặng Thị Phương Hải LỜI CAM ĐOAN Tên là: Đặng Thị Phương Hải Học viên: K14 Vật lí lí thuyết Vật lí toán Trường đại học Sư phạm Hà Nội Tôi xin cam kết đề tài “Nghiên cứu lý thuyết trật tự hợp kim β CuZn phương pháp thống kê momen” kết nghiên cứu riêng cá nhân tôi, tìm hiểu thực hướng dẫn bảo tận tình thầy giáo, PGS - TS Phạm Đình Tám Nếu có không trung thực luận văn xin hoàn toàn chịu trách nhiệm trước hội đồng khoa học Học viên thực hiện: Đặng Thị Phương Hải MỤC LỤC Mở đầu…………………………………………………………………… Nội dung…………………………………………………………………… Chƣơng 1: Các phƣơng pháp thống kê nghiên cứu trật tự hợp kim…………………………………………………………………… 1.1 Lý thuyết thống kê trật tự……………………………… 1.2 Phương pháp Kirkwood…………………………………… 1.3 Phương pháp giả hóa……………………………………… 12 Chƣơng 2: Năng lƣợng tự thông số trật tự hợp kim thay AB……………………………………………………………… 17 2.1 Momen biểu thức nhiệt động tinh thể loại nguyên tử…………………………………………………………………… 17 2.1.1 Momen…………………………………………………… 17 2.1.2 Các công thức tổng quát momen……………………… 18 2.1.3 Công thức tổng quát tính lượng tự theo phương pháp momen……………………………………………………………… 23 2.1.4 Biểu thức lượng tự biểu thức nhiệt động tinh thể………………………………………………………………… 23 2.2 Biểu thức lượng tự biểu thức thông số mạng hợp kim thay AB……………………………………………………… 27 2.2.1 Biểu thức lượng tự hợp kim………………… 27 2.2.2 Phương trình trạng thái biểu thức tính thông số mạng hợp kim………………………………………………………………… 29 2.3 Thông số trật tự phương trình xác định nhiệt độ trật tự T0 30 Chƣơng 3: Nghiên cứu trật tự hợp kim β – CuZn phƣơng pháp momen…………………………………………………… 35 3.1 Thế tương tác nguyên tử kim loại hợp kim 35 3.2 Phương trình trạng thái biểu thức thông số mạng hợp kim β – CuZn……………………………………………………………… 37 3.3 Thông số trật tự xa nhiệt độ trật tự T0 hợp kim β – CuZn……………………………………………………………………… 45 3.3.1 Thông số trật tự xa hợp kim β – CuZn……………… 50 3.3.2 Nhiệt độ trật tự hợp kim β – CuZn…………………… 46 3.4 Tính số thảo luận kết quả………………………………… 48 Kết luận…………………………………………………………………… 52 Công trình công bố liên quan đến nội dung luận văn………………… 53 Tài liệu tham khảo……………………………………………………… 54 MỞ ĐẦU Lý chọn đề tài Hợp kim vật thể mang tính kim loại (sáng, dẻo, dẫn điện nhiệt), chứa nhiều nguyên tố, chủ yếu nguyên tố kim loại, nguyên tố lại nguyên tố hợp kim hóa Do hợp kim có nhiều mặt ưu việt kim loại nguyên chất Hợp kim sử dụng rộng rãi nhiều lĩnh vực hàng hải, ứng dụng y tế, quân sự, thương mại, công nghiệp, khu dân cư sản xuất Vì vậy, nhiều ngành khoa học chọn hợp kim làm đối tượng nghiên cứu Có nhiều phương pháp để nghiên cứu trật tự hợp kim phương pháp Bragg – Williams, phương pháp Kirkwood, phương pháp giả hóa… cho phép giải thích nhiều tượng trật tự hợp kim; nhiên, lại không làm rõ thông số trật tự chịu ảnh hưởng tác động nhiệt độ áp suất Do đó, phương pháp thống kê – phương pháp momen đưa vào để giải vấn đề Phương pháp đơn giản cho kết giải tích kết số phù hợp tốt với thực nghiệm nghiên cứu tính chất tinh thể áp dụng thành công để nghiên cứu tính chất nhiệt động tinh thể loại nguyên tử Dựa sở áp dụng phương pháp momen, đưa đề tài “Nghiên cứu lý thuyết trật tự hợp kim β - CuZn phương pháp thống kê momen” Trong đề tài này, tìm hiểu thông tin trật tự tính chất nhiệt động hợp kim β - CuZn, hay gọi đồng thau, hợp kim phổ biến, có vai trò quan trọng ứng dụng rộng rãi thực tế thông qua biểu thức lượng tự Từ tìm phương trình trạng thái, biểu thức tính thông số mạng hợp kim β - CuZn, phương trình xác định thông số trật tự xa cân bằng, phương trình xác định nhiệt độ chuyển pha trật tự ảnh hưởng trật tự lên tính chất nhiệt động hợp kim β - CuZn Mục đích nghiên cứu Tìm hiểu tượng trật tự hợp kim ảnh hưởng nhiệt độ, áp suất tới thông số trật tự hợp kim β - CuZn Nhiệm vụ nghiên cứu - Xác định lượng tự biểu thức tính thông số mạng hợp kim β - CuZn - Xác định thông số trật tự phương trình xác định nhiệt độ trật tự hợp kim β - CuZn - Áp dụng tính số Đối tƣợng nghiên cứu Hợp kim β - CuZn, tính chất nhiệt động tính trật tự hợp kim β CuZn Phƣơng pháp nghiên cứu Phương pháp thống kê momen CHƢƠNG CÁC PHƢƠNG PHÁP THỐNG KÊ NGHIÊN CỨU TRẬT TỰ CỦA HỢP KIM 1.1 Lý thuyết thống kê trật tự Ta tìm biểu thức gần nhiệt động giá trị cân thông số trật tự xa toàn miền biến thiên chúng cách dựa vào sở cho lý thuyết thống kê trật tự xác định mẫu hợp kim đơn giản Trong lý thuyết này, số lượng biểu diễn tường qua lượng tương tác nguyên tử loại khác Do đó, nên lựa chọn mẫu xác định loại nhiệt động đặc biệt điểm chuyển pha loại Mẫu gần tương tác nguyên tử mẫu đơn giản giúp ta giải thích nhiều tượng liên quan đến trật tự nhiều trường hợp nên thường sử dụng lý thuyết thống kê Đối với vật rắn chịu nén nhỏ, coi thể tích V không đổi Khi ta sử dụng lượng tự  thay cho nhiệt động  Bằng phương pháp tổng quát lý thống kê, ta xác định lượng tự Tổng thống kê có dạng sau:  E  Z   exp   n  ,  kT  n (1.1) n số trạng thái hệ Ở đây, ta coi số trạng thái xác định số cấu hình i nguyên tử nút tập hợp số lượng tử (kí hiệu m) đặc trưng cho trạng thái electron dẫn, dao động nhiệt nguyên tử trạng thái liên quan đến bậc tự khác cho cấu hình i Năng lượng En hợp kim mô hình tương tác cặp xác định biểu thức: En = Ei + Em , (1.2) đó, Ei lượng cấu hình Em xem gần không phụ thuộc vào cấu hình i, xác định số lượng tử m Nếu biết trước loại nguyên tử chiếm nút mạng hợp kim xác định cấu hình lượng Ei Thay (1.2) vào (1.1) chuyển tổng theo n sang tổng theo i m, ta được:  E  Em   Em   Ei  Z   exp  i    exp    exp   kT  m   kT  i  kT  i,m (1.3)  E   E  Đặt: Z'   exp  i ; Z''=  exp  m   kT   kT  i m (1.4) Khi đó, (1.3) có dạng: Z  Z'.Z'' (1.5) Mà, lượng tự xác định:  '  kTln Z' Từ ta nhận biểu thức:    '  ''  '  kT ln Z'  ''  kT ln Z'' (1.6)  ' lượng tự cấu hình hợp kim, xác định cấu hình Việc giải toán tìm lượng tự để tính tổng thống kê dựa vào điều kiện cực tiểu từ tính chất cân hợp kim gặp nhiều khó khăn mặt toán học Chỉ giải toán mạng tinh thể chiều hai chiều Như vậy, để xác định thừa số cấu hình Z’ tổng thống kê mạng tinh thể ba chiều cần sử dụng phương pháp gần khác nhau, đơn giản phương pháp không tính tới tương quan, trật tự đặc trưng thông số trật tự xa Lý thuyết thống kê không tính tới tương quan phát triển công trình Gorsky Bragg – Williams Cụ thể sau: Xét hợp kim có nồng độ thành phần c cho, thông số trật tự xa xác định thông số độc lập 1, 2 , , q Khi biểu thức Z’ viết dạng: Z'   1 q Z1 q (1.7) Tổng lấy theo tất giá trị thông số trật tự xa, Z1 q tổng theo tất cấu hình tương ứng với giá trị cho trước thông số trật tự xa 1, 2 , , q xác định: Z1 q   1 q   i  E  exp  i   kT  (1.8) Đối với tinh thể vĩ mô có số lượng nguyên tử lớn thay lnZ’ thành ln Z1 ,2 , ,q với 1, 2 , , q giá trị cân thông số trật tự giá trị Z1 ,2 , ,q đạt cực đại nhiệt độ xác định Tính Z1 q giá trị khác thông số trật tự xa giá trị cân chúng từ điều kiện cực đại Z1 q hay cực tiểu lượng tự cấu hình 1, 2 , , q xác định giá trị cân 1, 2 , , q  '  kTln Z1 q (1.9) Do  '' không phụ thuộc vào 1, 2 , , q nên coi    '  '' với:  '    0; 1 1  '    0; ; 2 2  '   0 q q (1.10) Khi không xét tới tương quan tính gần lượng cấu hình E i tất cấu hình ứng với thông số trật tự xa 1, 2 , , q cho trước Trong trường hợp Ei xác định từ giả thiết hỗn loạn nguyên tử loại nút  E Thay (1.8) vào (1.9), khảo sát gần Ei = E đưa exp   khỏi  kT  dấu tổng theo i ta nhận biểu thức gần  ' sau:  '  E  kTln W , (1.11) W số hoán vị khác nguyên tử theo nút mạng cho biết thông số trật tự xa 1, 2 , , q Tuy nhiên, lý thuyết trật tự hoàn thiện phải tính tới tương quan hợp kim dựa vào việc áp dụng thủ thuật đặc biệt liên quan tới việc xác định gần Phương pháp Kirkwood phương pháp giả hóa thể điều 1.2 Phƣơng pháp Kirkwood Dựa sở khai triển lượng tự thành chuỗi theo lũy thừa  , Kirkwood hoàn thiện phương pháp xác định lượng tự hợp kT kim có tính tương quan Kết thu đạt độ xác cao tỉ số lượng trật tự kT nhỏ, nhiệt độ thấp độ xác giảm lý thuyết cho phép tìm vài số hạng khai triển 39 Mà, y    Cu, Zn  lại tính công thức (2.23): y2   2  A 0 3k 3 (3.12) Trong (3.12) k  ,   có dạng [7]: k      a      a  3a   (3.13) 21  4 8   a      a      a      a  54 9a 9a 9a (3.14) Thế tương tác có dạng: u   a   8  a   6  a  , với a  (3.15) a bán kính cầu phối vị thứ 2   u   a   8  a   6  a   (3.16)   a  tương tác Lennard – Jones có dạng: n m  r0   D   r0  m    n    a0   n  m   a0   a    (3.17) Thay (3.17) (3.16) được: r  8D   r0  m    n   u a0   n  m   a0   a0   n m 6D  3n m +  n  m   2n 3n n  r0  3m n    m m  a0  n m  2D   3n   r0  3m   r0   m   n     n   m       nm    a    a        r0     a0  m   40 Như vậy: n m  2D   3n   r0  3m   r0   m   n     n   m     u a0    nm    a    a       (3.18) Ta có đạo hàm:    2  m n  r0   Dmn  r0        a0   a  n  m   a   a    n m  r0   r0   Dmn  1  n     1  m     a0   a0 n  m   a0   a    m n  r0   r0   Dmn  1  m   m     1  n   n       a   a0  n  m   a0   a     3  4  n r  Dmn [ 1  n   n   n     a0   a0 n  m  a0  m r  - 1  m   m   m    ]  a0  Thay (3.17) vào (3.13): n m  r0   r0   4Dmn   2 1  n  1     m  1    k   a0    a0   3a 3a 02  n  m    a0   a    Đặt A  Dmn ; B1  n  1; B2  m  nm Như vậy: n n m  4AB1  r0   B1  a   k    1    3a  a   B1  r0    (3.19) 41 Thay (3.17) vào (3.14):   21  4  3  2 1   a     a    a    a    0 54 9a 9a 9a 03 n Dmn  21   r0  125  [ n  13n  n  21     9a  n  m    a0  m 125  21  r  -  m3  13m  m  21   ]   a0  Đặt A1  A2  21 125 n  13n  n  21 21 125 m  13m2  m  21 Như vậy: n n m  AA1  r0   A  a       1    9a  a   A1  r0    Ta tính đạo hàm (3.20) u  k  sau: , a a m n  u  2Dmn  3m   r0  3n   r0     m       n       a n  m    a    a      n n m  u  2A  r0   3n   3m   a   Hay:      n     m    a a  a   r 2        m n  r0   r0   k  4A  2  m m2    n n2     a a   a0   a        (3.21) 42 n n m  k  4A  r0   C2  a  , Hay:    1    a a  a   C1  r0    đó: C1  n  n  2; (3.22) C1  m2  m  Chú ý tới công thức gần đúng: 1  x n   nx  x 1 x 1 1 x Thay (3 19), (3.21) (3.12) vào (3 9): n Cu mCu  3n Cu   3mCu   a  [  n     m  ]+  Cu  Cu  r   2      0Cu  n Zn n Zn m Zn 2A Zn  r0Zn   3n Zn   3mZn   a  + ]}   [   n Zn     mZn    a0  a0   r 2      0Zn  2A  r  p a   { Cu  0Cu  12 a  a  n Cu n Cu  r 4A - { C1Cu  0Cu  a 3mCu B1Cu  a0  n Cu mCu   B   a x 1  2Cu    2B1Cu  r0Cu      C  a n Cu mCu  r  4A  x 1  2Cu   C1Zn  0Zn   C1Cu  r0Cu   a 3m Zn B1Zn  a0  n Zn n Zn m Zn   n Zn m Zn   B     a a C  1  2Zn   } x 1  2Zn    2B1Zn  r0Zn    C1Zn  r0Zn   Đặt YCu  a0 a ; YZn  Khi (3.23) có dạng: r0Cu r0Zn x (3.23) 43  3 3n Cu  n Cu  p a  {A Cu YCu   n Cu    n Zn  A Zn YZn  3n Zn   n Zn     3mCu  4  m   Cu     3mZn  4  m   Zn     n Cu mCu   YCu       n Zn m Zn }   YZn     n  B  n Cu mCu  Cu A Cu C  { C1Cu YCu [1   2Cu  2Cu  YCu  4a 3mCu B1Cu 2B C 1Cu   1Cu  + B2Cu C2Cu 2 n Cu mCu  YCu ]+ 2B1Cu C1Cu n  B  Zn A Zn C C1Zn YZn [1   2Zn  2Zn 3m Zn B1Zn  2B1Zn C1Zn   n Zn mZn   YZn  B2Zn C2Zn 2 n Zn mZn  YZn ]} 2B1Zn C1Zn Đối với Cu ta có giá trị: n  9; m=5,5; r0  2,5487A o A  48100; B1  8; B2  4,5 kB A1  6236; A  1664; C1  88; C2  33,75 Đối với Zn ta có giá trị: n  10; m=5,5; r0  2,7622A o A  20552; B1  9; B2  4,5 kB A1  7144; A  1664; C1  108; C2  33,75 Thay giá trị vào biểu thức (3.24): (3.24) 44  3 30 A Cu 39   35,5 23 9 p.10 a 10  { 1,38.10 YCu       5,5 kB      A Zn 310   35,5 23 10   1,38.10 YZn   10     5,5 kB       3,5  YCu         4,5  YZn  }      A Cu 1,38.1023 kB 4,5  { C1Cu YCu x 10 27 4a 10 3M Cu 1,66.10 B1Cu  B x 1   2Cu   2B1Cu A Zn 1,38.1023 C  3,5 B2Cu C2Cu  kB  2Cu  YCu  YCu   x 27 C1Cu  2B1Cu C1Cu 3M 1,66.10 B Zn 1Zn   B C 5  x C1Zn YZn 1   2Zn  2Zn   2B1Zn C1Zn  4,5 B2Zn C2Zn  YCu  }  YCu  2B C  1Zn 1Zn  (3.25) Thực phép tính rút gọn ta thu phương trình cuối có dạng sau: 6 13 4 11,5 3 8,5 5,214.106 pa13  5,12.10 a  0,49.10 a  3,54.10 a   0,008a 80  915,1a 04,5  0,61a 04  223,1a  578,35  (3.26) Thay (3.19), (3.20) vào (3.12) sử dụng công thức gần ta có: a A1  a  y     A 32B13  r0  n n m    3B A2   a  1    ,  =Cu,Zn ,      2B1 2A1   r0   (3.27) với   k BT ta có độ dời nguyên tử khỏi vị trí cân xác định công thức: a A1  a  y    A 32B13  r0  kB n n m    3B A2   a  1   T    2B 2A r      (3.28) 45 Khi đó:   Cu  Zn y  0,5  yCu  y Zn   0,5    3k 3k 3Zn Cu  a0 =0,5T{ A Cu kB a0 A Zn kB A1Cu 32B1Cu A1Zn 32B1Zn  a0     r0Cu   a0     r0Zn  n Zn nCu    nCu  mCu    3B A2Cu   a  2Cu 1       2B 2A r   1Cu  1Cu   0Cu  n Zn m Zn    3B    a A 1   2Zn  2Zn       2B1Zn 2A1Zn   r0Zn   (3.29) A ,A1 ,A 2 ,B1 ,B2 ,r0 kim loại Cu Zn vào kB Thay giá trị (3.29), thực tính toán, ta thu kết quả:  13,5 11 15,5 y  109 T 2,9a10  0,013a  1,14a  0,0033a  (3.30) Như vậy, khoảng cách nguyên tử hợp kim (thông số mạng) xác định sau:   14,5  , a  a 1  109 T 2,9a 90  0,013a12,5  1,14a10 0  0,0033a   (3.31) a0 xác định biểu thức (3.26) 3.3 Thông số trật tự xa nhiệt độ trật tự T0 hợp kim β – CuZn 3.3.1 Thông số trật tự xa hợp kim β – CuZn Thông số trật tự xa cân hợp kim β - CuZn xác định công thức:     0,     T,P,cA (3.32) 46  lượng tự hợp kim β – CuZn Thay  biểu thức (3.5), (3.6) vào (3.32), thực phép tính ta thu phương trình xác định thông số trật tự xa cân hợp kim β – CuZn sau:  3k T  k  k 2  k T 1  Zn Cu  B      B ln k Zn k Cu 1   16  Hay: (3.33)  2 1    k Zn k Cu ln     2  0     k Cu k Zn kT  n Zn n Zn m Zn  4A Zn B1Zn  r0Zn   B2Zn  a     1    a B r 3a   1Zn  0Zn     1   ln  {     4A B  r n Zn  B  a  n Cu mCu  Cu 1Cu 0Cu 2Cu    1    a B r 3a  1Cu  0Cu     (3.34) n Zn  n Cu  mCu  4A Cu B1Cu  r0Cu  B2Cu  a         a B r 3a 02   1Cu  0Cu     2 +  2}  0 n Zn  n Zn m Zn  k T B     4A Zn B1Zn r0Zn B2Zn a    1    a B r 3a   1Zn 0Zn      Thay giá trị A ,B ,r0 Cu Zn ta thu phương trình xác định thông số trât tự xa cân hợp kim β – CuZn:  0,0213a 3,5 1   0,685  0,0052a 04,5 2 ln   0,09a  0,5  0 3,5 4,5  1  a  0,0213a k BT  0,0052a (3.35) 3.3.2 Nhiệt độ trật tự hợp kim β – CuZn Sự chuyển trật tự - không trật tự hợp kim β – CuZn chuyển pha loại nhiệt độ chuyển trật tự - không trật tự tính toán từ điều kiện chuyển pha loại 2: 47   T  TC   T > TC (3.36) Từ (3.36) ta thấy nhiệt độ gần với TC, thông số trật tự  Do đó, điểm gần TC, vế phải phương trình (3.33) có dạng:  k BT   ln  k BT 1  (3.37) Từ (3.33) (3.37), ta tìm phương trình nhiệt độ trật tự TC hợp kim β – CuZn sau: TC  Hay: TC   k B  k Zn  k Cu 2 2 8k Zn k Cu (3.38)  k B  k Zn k Cu    2    k Cu k Zn  n Zn n Zn m Zn  4A Zn B1Zn  r0Zn   B2Zn  a     1    a B r 3a   1Zn  0Zn     2  TC  { (  k B 4A B  r n Cu  B  a n Cu mCu  Cu 1Cu 0Cu 2Cu    1    3a   a   B1Cu  r0Cu   B  a n Cu mCu  1  2Cu     B1Cu  r0Cu   +  2)  2}1 n Zn  n Zn m Zn  4A Zn B1Zn  r0Zn  B2Zn  a         a B r 3a 02  1Zn  0Zn     4A Cu B1Cu  r0Cu    3a 02  a0  n Cu (3.39) Thay giá trị A ,B ,r0 kim loại Cu Zn ta thu nhiệt độ trật tự TC hợp kim β – CuZn: 48   0,0213a 3,5 2  0,685  0,0052a 04,5 TC   0,09a  1,5   k B  a  0,0213a 3,5  0,0052a 04,5  1 (3.40) 3.4 Tính số thảo luận kết Sử dụng kết thu thông số mạng, thông số trật tự hợp kim β – CuZn như: phương trình trạng thái 0oK, áp suất p (3.26); biểu thức tính thông số mạng (3.31), phương trình xác định thông số trật tự nhiệt độ T, áp suất p (3.35), ta tính thông số mạng, thông số trật tự hợp kim β – CuZn nhiệt độ áp suất khác Kết tính số cho bảng biểu diễn hình hình 2: Hình 1: Sự phụ thuộc thông số mạng hợp kim β – CuZn vào nhiệt độ áp suất 49 Bảng 1: Giá trị thông số mạng a, thông số trật tự  nhiệt độ áp suất khác nhau: P(Kbar) T(K) a(A0) TN [1] 1   200 300 400 500 600 700 2,6271 2,6380 2,6488 2,6597 2,6705 2,6814 2,6923 2,7031 a(A0) 100 100  a(A0) 50 0,9991 0,9870 0,9466 0,8616 0,7074 0,3996 0,99 0,95 0,89 0,74 0,41 2,5923 2,6017 2,6110 2,6204 2,6297 2,6391 2,6485 2,6578 1 0,9991 0,9872 0,9472 0,8637 0,7134 0,4216 2,5631 2,5734 2,5818 2,5901 2,5984 2,6068 2,6151 2,6234 1 0,9991 0,9872 0,9476 0,8649 0,7168 0,4331 Hình 2: Sự phụ thuộc thông số trật tự hợp kim β – CuZn vào nhiệt độ áp suất 50 So sánh kết tính lý thuyết với giá trị thực nghiệm áp suất p = (áp suất khí quyển) có phù hợp tốt Trong nhiệt độ, thông số mạng giãn thông số trật tự tăng áp suất tăng Các kết giải thích được: ý tăng áp suất, tinh thể bị co lại, khoảng cách nguyên tử giảm (thông số a giảm), đồng thời khả hoán vị nguyên tử khác loại nút mạng bị cản trở nên khả phá vỡ tính trật tự nhiệt độ giảm (thông số trật tự giảm chậm lại tăng nhiệt độ áp suất cao) Các kết thực nghiệm xác nhận Như vậy, phương trình hệ thức thu để xác định thông số mạng thông số trật tự hợp kim β – CuZn có dạng giải tích đơn giản, dễ tính số kết số phù hợp với thực nghiệm 51 KẾT LUẬN Phương pháp momen phương pháp thống kê mới, áp dụng thành công nghiên cứu tính chất nhiệt động tinh thể loại nguyên tử bước đầu áp dụng tốt tìm hiểu trật tự hợp kim thay AB, cụ thể đề tài đưa ra, phương pháp đưa vào nghiên cứu hợp kim β - CuZn có cấu trúc lập phương tâm khối Đề tài đạt kết sau: - Tìm hiểu số phương pháp thống kê nghiên cứu hợp kim trật tự như: phương pháp Kirkwood, phương pháp giả hóa Đặc biệt, làm quen ứng dụng phương pháp thống kê momen nhà khoa học Việt Nam phát triển ứng dụng có hiệu - Áp dụng phương pháp momen tìm phương trình xác định thông số mạng, phụ thuộc thông số trật tự vào nhiệt độ áp suất hợp kim β CuZn có dạng giải tích đơn giản - Đưa kết tính số nhận xét kết thu hợp kim β CuZn Kết số phù hợp với thực nghiệm khoảng nhiệt độ áp suất lớn Sự phù hợp lý thuyết thực nghiệm nghiên cứu trật tự hợp kim β - CuZn cho phép mở rộng hướng nghiên cứu hợp kim hai thành phần khác hợp kim nhiều thành phần 52 CÔNG TRÌNH ĐÃ CÔNG BỐ LIÊN QUAN ĐẾN NỘI DUNG LUẬN VĂN Pham Dinh Tam, Pham Duy Tan, Nguyen Quang Hoc, Dang Thi Phuong Hai, Nguyen Thi Binh, Order theory of alloy β – CuZn, Proc Natl Conf Theor Phys 36 (2011), pp 195 – 200 53 TÀI LIỆU THAM KHẢO A.A.Cmirnov, Lý thuyết động học phân tử kim loại, Nauka, Mosk., 1966 (Tiếng Nga) K.Butilenco cộng sự, Ảnh hưởng áp suất tới trật tự hợp kim CuZnvà AgZn, Kiev, Nauka, Đumka, 1979 (tiếng Nga) K.Kanniuca, Lý thuyết trật tự hợp kim cấu trúc lập phương tâm khối loại  - đồng thau áp suất tính tới tương quan, .M M , 31, 3(1971)478 (tiếng Nga) P.E.A.Turchi et al., First-Principles Study of Phase Stability in Cu-Zn Substitutional Alloys, Phys.Rev.Let.,67,13(1991)1779 V.F.Degtyareva et al., Stability of Hume-Rothery Phases in Cu-Zn Alloys at Pressures up to 50Gpa, J.Phys : Condens Matter, 17(2005)7955 K.Masuda-Jindo, Vu Van Hung, Pham Dinh Tam, Application of Statistical moment Method to thermodynamic quantities of metals and alloys, Calphad, 26, 1(2002)15 Pham Dinh Tam, The melting temperature for binary alloys AB at various pressures, VNU.Jour.of Sci., 2(1999)35 Pham Dinh Tam, The lattice spacings for binary alloysAB, Comm.in Phys., 2(1998)78 Shuzen, G J Davies, Calculation of Lennard – Jones (nm) potential energy parameter for metals, Phys Stat Sol (a), 78(1983) 595 10 K.Masuda-Jindo, Vu Van Hung, Pham Dinh Tam, Thermodynamic Quantities of Metals investigated by an analytic Statistical Moment Method, Phys.Rev.B,9(2003)094301