ĐẠI HỌC QUỐC GIA HÀ NỘI
TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN


NGUYỄN TRUNG ĐÔ NGHIÊN CỨU TÍNH CHẤT ĐIỆN TỬ CỦA MỘT SỐ HỢP
CHẤT SỬ DỤNG PHƢƠNG PHÁP PHIẾM HÀM MẬT ĐỘ
LUẬN VĂN THẠC SĨ KHOA HỌC
Hà Nội – 4-2014
ĐẠI HỌC QUỐC GIA HÀ NỘI
TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN

LỜI CẢM ƠN

Trƣớc hết, lời cảm ơn sâu sắc nhất của em xin đƣợc gửi tới thầy giáo hƣớng
dẫn của em, GS.TS.Bạch Thành Công , ngƣời trực tiếp chỉ dẫn và giúp đỡ em nhiều
nhất trong thời gian học tập và hoàn thành luận văn tốt nghiệp của mình.
Em cũng xin gửi lời cảm ơn tới toàn thể các quý thầy cô và tập thể các cán
bộ công nhân viên bộ môn Vật lý Chất rắn cùng gia đình bạn bè , những ngƣời đã
động viên, dạy bảo, chăm sóc và cho em những ý kiến đóng góp quý báu và hết sức
bổ ích giúp em hoàn thành luận này đƣợc dễ dàng và thuận lợi hơn.
Nhân đây, em cũng xin đƣợc gửi lời cảm ơn tới các thầy cô và cán bộ tại
Khoa Vật lý đã hết sức tạo điều kiện thuận lợi cho em trong cả quá trình học tập và
viết luận văn.
Xin cám ơn đề tài QG.12.01 đã hỗ trợ để thực hiện luận văn này.
Hà Nội, ngày tháng năm 2014
Sinh Viên

Nguyễn Trung Đô 4
MỤC LỤC

LỜI CẢM ƠN
MỞ ĐẦU 1
CHƢƠNG I. TỔNG QUAN VỀ VẬT LIỆU PERMALLOY VÀ PEROVSKITE 2
1. Cấu trúc và tính chất của vật liệu Permalloy 2
2.Cấu trúc cơ bản trong vật liệu Perovskite 4

3
47
KẾT LUẬN 50
TÀI LIỆU THAM KHẢO 51

5
DANH MUC HÌNH BẢNG TRONG LUẬN VĂN

Hình 1.1: Sự thay đổi của độ từ thẩm ban đầu của permalloy theo hàm lƣợng Ni
đƣợc chế tạo theo hai phƣơng pháp cán nóng và cán lạnh.
Hình 1.2: Sự phụ thuộc của từ độ bão hòa và hàm lƣợng Ni trong các hợp kim
permalloy.
Hình 1.3: Cấu trúc tinh thể perovskite ABO
3
thuần.
Hình 2.1: Thế Muffin-tin
Hình 3.1: Sự phụ thuộc của năng lƣợng vào thể tích của Fe
Hình 3.2: Sự phụ thuộc của năng lƣợng vào thể tích của permalloy 10 cho hai mô
hình cấu trúc bcc và fcc.
Hình 3.2: Sự phụ thuộc của năng lƣợng vào thể tích của permalloy 10 cho hai mô
hình cấu trúc bcc và fcc.
Hình 3.3: Sự phụ thuộc của năng lƣợng vào thể tích của permalloy 20 cho hai mô
hình cấu trúc bcc và fcc.
Hình 3.4: Sự phụ thuộc của năng lƣợng vào thể tích của permalloy 30 cho hai mô
hình cấu trúc bcc và fcc.
Hình 3.5: Sự phụ thuộc của năng lƣợng vào thể tích của permalloy 40 cho hai mô
hình cấu trúc bcc và fcc.

Hình 3.17: Mật độ trạng thái điện tử cho hai hƣớng spin lên và xuống ( )
của Ni-Fe-20.
Hình 3.18: Mật độ trạng thái điện tử cho hai hƣớng spin lên và xuống ( )
của Ni-Fe-30.
Hình 3.19: Mật độ trạng thái điện tử cho hai hƣớng spin lên và xuống ( )
của Ni-Fe-40.
Hình 3.20: Mật độ trạng thái điện tử cho hai hƣớng spin lên và xuống ( )
của Ni-Fe-50.
Hình 3.21: Mật độ trạng thái điện tử cho hai hƣớng spin lên và xuống ( )
của Ni-Fe-60.
Hình 3.22: Mật độ trạng thái điện tử cho hai hƣớng spin lên và xuống ( )
của Ni-Fe-70.
Hình 3.23: Mật độ trạng thái điện tử cho hai hƣớng spin lên và xuống ( )
của Ni-Fe-80.
Hình 3.24: Mật độ trạng thái điện tử cho hai hƣớng spin lên và xuống ( )
của Ni-Fe-90.
Hình 3.25: Sự phụ thuộc của năng lƣợng và thể tích ô cơ sở của Ni
56.5
Fe
19.0
Ga
24.5

trong hai pha cấu trúc fcc (a) và bcc (b).

7
Hình 3.26: Mật độ trạng thái điện tử cho hai hƣớng spin lên và xuống ( )

1
MỞ ĐẦU Trong thời gian hiện nay permalloy (hợp kim của Niken và sắt) với độ từ
thẩm cao, perovskite có độ dẫn điện tốt nhƣ LaNiO
3
đƣợc sử dụng nhiều trong chế
tạo các cảm biến từ điện, linh kiện đa chức năng. Gần đây permalloy đƣợc dùng nhƣ
một vật liệu thành phần quan trọng để chế tạo các sensor địa từ [1]. LaNiO
3
đƣợc
dùng để làm điện cực cho các linh kiện [2], làm vật liệu xúc tác [3]. Hiện nay để
thiết kế linh kiện với các tham số vật lý cần thiết ngƣời ta có thể dùng phƣơng pháp
ab-initio tính toán thành phần, đặc trƣng của vật liệu đòi hỏi với độ chính xác cao.
Mục tiêu của luận án là sử dụng trên phƣơng pháp phiếm hàm mật độ và gói phần
mềm AKAI-KKR để tính toán các tính chất điện từ của vật liệu permalloy và
LaNiO
3
nhằm góp phần làm sáng tỏ cơ chế vật lý của các kết quả thực nghiệm và
góp phần định hƣớng ứng dụng các vật liệu này cho các mục đích khác nhau.

hợp thức là Ni
1-x
Fe
x
với giá trị x thay đổi từ 20% đến 85%. Trong các tài liệu về từ
học và trong kỹ thuật, ngƣời ta gọi tên của hợp kim này tƣơng ứng với tỉ lệ niken, ví
dụ Permalloy75 là hợp kim permalloy có chứa 75% nguyên tử là niken (hay
Ni
75
Fe
25
). Hợp kim permally là hợp kim có từ tính, thƣờng đƣợc sử dụng trong các
ứng dụng về từ học. Tên gọi permalloy xuất phát từ chữ ghép per (trong chữ
permeability, có nghĩa là từ thẩm), với từ alloy có nghĩa là hợp kim, do permalloy là
hợp kim có độ từ thẩm rất cao
Permalloy với tỉ lệ 75% nguyên tửniken đƣợc gọi là permalloy chuẩn
(Standard permalloy). Permalloy đƣợc pha tạp một số nguyên tố khác (ví dụ nhƣ
Môlipđen - Mo) đƣợc gọi là Supermalloy. Trong kỹ thuật, Permalloy thƣờng đƣợc
viết tắt là Py.
Ở dạng khối và đơn tinh thể, permalloy có cấu trúc lập phƣơng tâm mặt điển
hình [4],[5], cấu trúc tinh thể có thể bị thay đổi tùy theo phƣơng pháp chế tạo (ví dụ
cấu trúc lục giác xếp chặt khi ở dạng màng mỏng chế tạo bằng phƣơng pháp epitaxy
chùm phân tử [6]). Hằng số mạng của permalloy phụ thuộc vào thành phần hợp kim
và quy luật chƣa đƣợc xác định một cách rõ ràng. Thí dụ hợp kim permalloy75 có
cấu trúc lập phƣơng tâm mặt với hằng số mạng a = 0.3555 nm, khối lƣợng riêng
8,57.10
3
kg/m
3
, thuộc nhóm không gian Pm-3m [7], trong khi hợp kim permalloy50

trở suất rất thấp. Hợp kim permalloy có thể cho hiệu ứng từ điện trở khoảng 5% ở
nhiệt độ phòng.
Trong nội dung luận văn này, chúng tôi chỉ tính toán một vài tính chất của
hợp kim Permalloy và ngoài ra còn có một số tính toán một số tính chất của vật liệu
Ni-Fe-Ga
2.Cấu trúc cơ bản trong vật liệu Perovskite
Vật liệu perovskite ABO
3
đƣợc bắt đầu biết đến từ đầu thế kỷ 19.Thời gian
đầu các nhà khoa học cũng chƣa thực sự quan tâm đến những vật liệu này.Trong
thời gian gần đây, đã có rất nhiều nghiên cứu về vật liệu perovskite. Do các vật liệu
perovskite ABO
3
có độ bền nhiệt và thể hiện các tính chất vật lý đặc sắc : sắt điện,
sắt từ, nhiệt, quang … trong vùng nhiệt độ rất rộng nên nếu dùng chúng chế tạo các
linh kiện thì các linh kiện này có miền hoạt động lớn. Ngoài ra, khi pha tạp thay thế
một số nguyên tố nhƣ Ba, Sr, Fe, Ni …. vào vị trí A hoặc B sẽ dẫn đến một số hiệu
ứng vật lý lý thú nhƣ: hiệu ứng nhiệt điện, hiện ứng từ nhiệt, từ trở khổng lồ… Điều
đó đã mở ra những ứng dụng mới về vật liệu perovskite trong một số lĩnh vực công
nghiệp hiện đại nhƣ điện tử, thông tin, làm lạnh mà không gây ô nhiễm môi trƣờng.

Hình 1.3: Cấu trúc tinh thể perovskite ABO
3
thuần
Hợp chất ABO
3
thuần có cấu trúc tinh thể lý tƣởng nhƣ hình 1.3. Ô mạng cở
sở là hình lập phƣơng với các thông số mạng a=b=c và . Ở đây
cation A nằm tại các đỉnh của hình lập phƣơng, còn cation B có bán kính nhỏ hơn
nằm tại tâm của hình lập phƣơng. Cation B đƣợc bao quanh bởi 8 cation A và 6

’
)O
3
với thì vật liệu này đƣợc gọi là vật liệu ABO
3

biến tính. Trong đó, A có thể là các nguyên tố nhƣ La, Nd, Pr thuộc họ đất hiếm còn
A
’
là các kim loại kiềm thổ nhƣ Sr, Ba, Ca … hoặc các nguyên tố nhƣ: Ti, Ag, Bi,
…ion B có thể các nguyên tố nhƣ là Mn, Co trong khi B
’
là các nguyên tố nhƣ Fe,
Ni, Ca…
Trong các perovskite ABO
3
bị biến tính khi pha tạp sẽ xuất hiện trạng thái
hỗn hợp hóa trị và sai lệch cấu trúc.
Sự sai lệch cấu trúc tinh thể đƣợc đánh giá thông qua thừa số dung hạn t do
Goldchmit đƣa ra :

Với R
A
, R
B
, R
O
lần lƣợt là bán kính của các ion A
2+
(A


6
CHƢƠNG II.
PHƢƠNG PHÁP PHIẾM HÀM MẬT ĐỘ VÀ CHƢƠNG TRÌNH AKAI-KKR

Akai-KKR là một gói phần mềm đƣợc sử dụng để tính toán cấu trúc điện tử
của kim loại, bán dẫn và hợp chất dựa trên phép gần đúng mật độ địa phƣơng (LDA)
và phép gẫn đúng tổng quát gradient (GGA) của lý thuyết phiếm hàm mật độ.
Phần mềm này sử dụng phƣơng pháp hàm KKR-Green và có tốc độ tính toán
cao cũng nhƣ độ chính xác. Ngoài ra, phép gần đúng CPA (coherent potential
approximation – gần đúng thế kết hợp) cũng đƣợc đƣa vào phần mềm này cho phép
KKR có thể áp dụng tính toán không chỉ với tinh thể mà còn với các hệ không đồng
nhất nhƣ các hệ pha tạp, hợp kim thế ngẫu nhiên và các tinh thể hỗn hợp.
1. PHƢƠNG PHÁP PHIẾM HÀM MẬT ĐỘ - DFT [10]
Ý tƣởng dùng hàm mật độ hạt để mô tả các tính chất của hệ electron đƣợc
nêu trong các công trình của Llewellyn Hilleth Thomas và Enrico Fermi ngay từ khi
cơ học lƣợng tử mới ra đời. Đến năm 1964, Pierre Hohenberg và Walter Kohn đã
chứng minh chặt chẽ hai định lý cơ bản là nền tảng của lý thuyết phiếm hàm mật
độ. Hai định lý khẳng định năng lƣợng ở trạng thái cơ bản là một phiếm hàm của
mật độ electron, do đó về nguyên tắc có thể mô tả hầu hết các tính chất vật lý của hệ
điện tử qua hàm mật độ. Một năm sau, W. Kohn và Lu Jeu Sham nêu ra qui trình
tính toán để thu đƣợc gần đúng mật độ electron ở trạng thái cơ bản trong khuôn khổ
lý thuyết DFT. Từ những năm 1980 đến nay, cùng với sự phát triển tốc độ tính toán
của máy tính điện tử, lý thuyết DFT đƣợc sử dụng rộng rãi và hiệu quả trong các
ngành khoa học nhƣ: vật lý chất rắn, hóa học lƣợng tử, vật lý sinh học, khoa học vật
liệu, . W. Kohn đã đƣợc ghi nhận những đóng góp của ông cho việc phát triển lý
thuyết phiếm hàm mật độ bằng giải thƣởng Nobel Hóa học năm 1998
1.1. Một số khái niệm cơ bản
a) Phiếm hàm
Hàm số biểu diễn một quy luật để đi từ biến số x đến một giá trị f(x).Phiếm

1
3
(3 )
10
C


,thành phần thứ hai là năng lƣợng điện tử trong trƣờng thế
ngoài thành phần thứ 3 là năng lƣợng trao đổi địa phƣơng với
1/3
2
33
()
4
C


và
thành phần cuối cùng là năng lƣợng Hartree tĩnh điện học.
Năng lƣợng và mật độ trạng thái cơ bản có thể tìm đƣợc băng cực tiểu hóa
phiếm hàm
[]En
trong toàn bộ các mật độ khả dĩ
()nr
liên hệ tới tổng số điện tử
theo biểu thức.

9
Thế ngoài v(r) được xác định hoàn toàn bởi mật độ electron trong trạng
thái cơ bản
Bởi vì xác định số electron nên có thể suy ra rằng cũng xác định
hàm sóng tráng thái cơ bản và tất cả các tính chất của hệ, Lƣu ý rằng v(r) không
chỉ giới hạn trong lực tƣơng tác Coulomb.
Nhƣ vậy, xác định N, v và vì thế xác định tất cả các tính chất ở trạng thái
cơ bản nhƣ động năng T( ), thế năng V( )và tổng năng lƣợng E( ). Ta có phƣơng
trình :
(2.6)
trong đó
+ Số hạng phi cổ điển
là thế năng Hatree. Dạng phi cổ điển là một đại
lƣợng rất khó nắm bắt nhƣng rất quan trọng, là phần chủ yếu của ‘năng lượng trao
đổi – tương quan’ đƣợc xác định ở những phần sau.
Định lý Hohenberg-Kohn thứ hai nói về nguyên lý biến phân năng lƣợng :
với mật độ thử sao cho và ta có :
trong đó, là phiếm hàm năng lƣợng xác định bởi (2.6).
Điều này tƣơng tự với nguyến lý biến phân của hàm sóng .Dựa
trên định lý thứ nhất xác định , Hamilton và hàm sóng , có thể đƣợc dùng
nhƣ hàm thử cho bài toán thế ngoài v. Nhƣ vậy, Giả sử tính khả vi của , nguyên lý biến phân đòi hỏi mật độ trạng thái
cơ bản thỏa mãn điều kiện
(2.8)
Sử dụng phƣơng trình (2.6) ta có phƣơng trình Euler-Lagrange :
(2.9)

với N obitan và với các obitan còn lại.

11
Tƣơng tự với các xác định phiếm hàm tổng hợp , Kohn và Sham đƣa
vào một hệ tham chiếu không tƣơng tác có mật độ electron trạng thái cơ bản chính
xác bằng , tƣơng ứng với Hamilton

Trong đó không có số hạng đẩy electron. Với hệ này, sẽ cố một hàm song định thức
chính xác trạng thái cơ bản:

Trong đó là trạng thái thấp nhất của Hamilton-một electron đƣợc xác định bởi
:

Động năng đƣợc xác định bởi (2.12):

Và mật độ đƣợc phân tích nhƣ trong (2.13).
Ý tƣởng cơ bản của Kohn-Sham là có thể thay bài toán nhiều electron bằng
một tâp tƣơng đƣơng chính xác các phƣơng trình tự hợp cho bài tóan một electron.
Phiếm hàm năng lƣợng tổng cộng của hệ có thể đƣợc viết dƣới dạng tổng của một
số số hạng:

Với một tập cố định các hạt nhân nguyên tử ở . Cả ba số hạng đều là các hàm
của mật độ điện tử . Phƣơng trình này tƣơng đƣơng với phƣơng pháp Hartree
nhƣng số hạng chứa các hiệu ứng tƣơng quan – trao đổi và động năng hạt

trong đó là động năng của hệ các electron không tƣơng tác có mật độ
và là năng lƣợng trao đổi và tƣơng quan của hệ tƣơng tác. 13
Tổng hợp lại, các phƣơng trình obitan Kohn-Sham dƣới dạng chính tắc là: với thế tƣơng quan – trao đổi là:
Trong đó Ψ là các obitan chuẩn trực giao cũng đƣợc gọi là các obitan Kohn-
Sham.Các phƣơng trình này đều là phi tuyến và chỉ có thể giải bằng phƣơng pháp
lặp. Năng lƣợng tổng cộng đƣợc xác định là: trong đó Cũng giống nhƣ lý thuyết Hartree-Fock, năng lƣợng tổng cộng của electron không
bằng tổng năng lƣợng obitan.

trong đó, và lần lƣợt là tích phân Coulomb và tích phân trao đổi
2.CÁC PHƢƠNG PHÁP GẦN ĐÚNG
2.1. Gần đúng mật độ địa phương (LDA - Local Density Approximation)
Phƣơng pháp gần đúng mật độ địa phƣơng đối với phiếm hàm tƣơng quan và
trao đổi dùng để tính toán năng lƣợng tƣơng quan trao đổi trên từng hạt của khí
X a
(hay còn gọi là phƣơng pháp HFS có thông số điều chỉnh

, đặt
trƣớc thế trao đổi). Những phiếm hàm tƣơng quan LDA đầu tiên có chất lƣợng chƣa

15
cao trong những vấn đề thực hiện khai triển nhiễu loạn hệ nhiều hạt ở giới hạn mật
độ thấp.
2.2.Phương pháp gần đúng Gradient suy rộng (Generalized Gradient
Approximation)
Phƣơng pháp gần đúng gradient suy rộng (GGA) là một dạng mở rộng của
phiếm hàm LDA có tính đến Gradient của mật độ điện tử.Trong phƣơng pháp GGA
này phiếm hàm (tổng quát) của năng lƣợng trao đổi tƣơng quan có dạng.
(2.35)
Những bƣớc quan trọng để dẫn tới GGA phần lớn đƣợc thực hiện bởi Perdew
và các cộng sự của ông. Theo ông, GGA có thể đƣợc viết thuận tiện dựa trên một
hàm giải tích đƣợc biết đến nhƣ là một thừa số gia tăng, .
(2.36)
với là năng lƣợng trao đổi của khí điện tử khi không phân cực
Một phiếm hàm đƣợc sử dụng nhiều cho chất rắn là PW91, đƣợc phát triển
bởi Perdew và Wang. Nó đƣợc xây dựng theo lối phi kinh nghiệm, đƣợc xác định từ
các hệ thức cơ học lƣợng tử chính xác. Trong phiếm hàm PW91, thừa số gia tăng
trao đổi có dạng.
(2.37)
s là gradient mật độ không thứ nguyên.
Năng lƣợng tƣơng quan có bổ chính spin trong PW91 có thể viết dƣới dạng.

nồng độ là x
1,
x
2,
x
3
, x
n
. Cho rằng nguyên tử A
i
nằm tại gốc của thế hiệu dụng. Hàm
Green là:

và là hàm Green và ma trận t của môi trƣờng hiệu dụng. CPA là một
phƣơng pháp gần đúng hiệu quả trong việc xác định .
Sử dụng phƣơng trình tự hợp:

Phƣơng trình có nghĩa là xác định hàm Green của môi trƣờng bằng cách lấy
trung bình có trọng số của hàm Green của nguyên tử thành phần đƣợc đặt tại gốc
của môi trƣờng hiệu dụng.
3. PHƢƠNG PHÁP HÀM GREEN [11]
3.1. Bài toán vị trí đơn
Coi sự tán xạ đƣợc gây ra bởi một thế đơn tại gốc. Thế này bằng không tại
một miền cách xa gốc.
3.1.1. Phương trình Shrodinger
Phƣơng trình Shrodinger cho trạng thái dừng của điện tử có dạng
=E . (3.1)
Nếu thế là thế cầu, lời giải có thể đƣợc tách ra thành hai phần khi sử
dụng hệ tọa độ cầu
= (3.2)

dụng hàm Hankel loại một

= + , (3.10) 18

= (3.11)
Kết quả là vấn đề về sự phân tán trở thành vấn đề xác định ma trận t đƣợc định
nghĩa bởi :
= - (3.12)
3.1.3. Wronskian
Ta xây dựng lời giải của phƣơng trình Shrodinger nằm ngoài thế tác dụng từ hệ thức
Wronskian. Từ phƣơng trình Shrodinger cho bán kính:

(3.13)

(3.14)
Các tích phân sau đây đƣợc đánh giá nhƣ sau:

(3.15)

= - - (3.16)
Số hạng đầu tiên trong phƣơng trình (3.15) đƣợc xác định từ dạng tiệm cận của
. Wronskian đƣợc định nghĩa nhƣ sau
= (3.17)
Nhân phƣơng trình (3.15) với và , phƣơng trình (3.16) với