ĐẠI HỌC QUỐC GIA HÀ NỘI
TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN
------------------
BÙI THỊ NHUNG
NGHIÊN CỨU TÍNH CHẤT ĐIỆN TỬ CỦA MỘT SỐ
PEROVSKITE TỪ TÍNH PHA TẠP ĐẤT HIẾM
Chuyên ngành: Vật lý chất rắn
Mã số: 60440104
LUẬN VĂN THẠC SĨ KHOA HỌC
NGƢỜI HƢỚNG DẪN KHOA HỌC:
TS. Bạch Hƣơng Giang
GS.TS. Bạch Thành Công
Hà Nội – Năm 2014
MỤC LỤC
MỞ ĐẦU ......................................................................................................................... 5
CHƢƠNG 1: TỔNG QUAN VỀ PEROVSKITE VÀ VẬT LIỆU CaMnO3 ............ 7
1.1.Tổng quan về vật liệu perovskite ............................................................................ 7
1.1.1.Cấu trúc tinh thể ................................................................................................... 7
1.1.2. Cấu hình điện tử ................................................................................................... 9
1.1.3. Các tƣơng tác vi mô dẫn tới tính chất từ trong hệ vật liệu perovskite ......... 13
1.2. Hệ CaMnO3 pha tạp.............................................................................................. 17
CHƢƠNG 2: LÝ THUYẾT CẤU TRÚC ĐIỆN TỬ CHẤT RẮN VÀ PHƢƠNG
4
MỞ ĐẦU
Vật liệu perovskite được bắt đầu biết đến từ đầu thế kỷ thứ IX. Công thức chung
của loại vật liệu này là AB
trong đó A là kim loại hóa trị 2 (Ca, Sr…), B thuộc nhóm
kim loại chuyển tiếp (Mn, Fe, Ti), loại vật liệu này có độ bền nhiệt rất cao nên có thể
hoạt động trong môi trường nhiệt độ cao. Do có nhiều tính chất điện - từ - hóa đặc biệt
khác nhau nên perovskite có rất nhiều ứng dụng và được coi là một trong những vật
liệu rất lý thú hiện nay. Với tính chất từ điện trở siêu khổng lồ (Colossal
Magnetoresistance- CMR), perovskite từ tính đang được nghiên cứu sử dụng cho các
linh kiện spin tử (spintronics) và các cảm biến từ siêu nhạy. Nhiều Perovskite CMR có
độ dẫn điện đủ lớn và bền vững trong vùng nhiệt độ cao (100oC -1000oC) ... nên là vật
liệu hữu ích để chế tạo các linh kiện điện tử hoạt động trong điều kiện cực đoan. Ngoài
ra, perovskite với các tính chất hấp phụ và xúc tác còn được sử dụng trong các pin
nhiên liệu (fuel cells).
Một trong những loại oxide perovskite được chú ý là CaMnO3 sạch và pha tạp
kim loại đất hiếm hoặc Ytrium (Ca1-xRxMnO3, R= La, Pr, Eu,…Y). Các hợp chất này
thể hiện sự đa dạng về cấu trúc,tính chất từ, đặc biệt là tính chất điện và nhiệt điện
trong vùng nhiệt độ cao.Theo các nghiên cứu thực nghiệm (thí dụ [1]) CaMnO3 là dẫn
điện có điện trở suất khoảng 8.102Ωcm, hệ số Seebeck khá lớn khoảng 200µV/K ở
nhiệt độ phòng. Tuy nhiên các tính toán lý thuyết đều cho kết quả vật liệu là điện môi
phản sắt từ trong trạng thái cơ bản với khe năng lượng xấp xỉ 1.02 eV (xem [17] và tài
liệu trích dẫn).Trong luận văn này, chúng tôi tập trung chủ yếu phân tích tính chất điện
tử của vật liệu khối CaMnO3 pha trực thoi (orthorhombic), cấu trúc không từ tính (ở
nhiệt độ cao vật liệu này là không từ tính) và ảnh hưởng của việc pha tạp một số kim
loại đất hiếm và khuyết oxy dẫn đến tính chất dẫn điện của vật liệu. Tất cả các tính
CHƢƠNG 1: TỔNG QUAN VỀ PEROVSKITE VÀ VẬT LIỆU CaMnO3
1.1.Tổng quan về vật liệu perovskite
1.1.1.Cấu trúc tinh thể
Perovskite có công thức chung là ABO3, trong đó A là kim loại hóa trị 2 (Ca,
Sr…), B thuộc nhóm kim loại chuyển tiếp (Mn, Fe, Ti…). A và B là các ion (cation) có
bán kính khác nhau.
Cấu trúc tinh thể lý tưởng của vật liệu perovskite là dạng lập phương với nhóm
không gian Pm3m. Các ion A2+ chiếm vị trí đỉnh, ion B4+ chiếm vị trí tâm khối và các
ion O2- chiếm vị trí tâm mặt của hình lập phương (xem hình 1.1). Các ion kim loại
chuyển tiếp B4+ tâm khối tạo thành phối trí bát diện (Octahedral coordination) với các
ion Oxy tâm mặt (thường được gọi là bát diện BO6). Giữa các ion này luôn tồn tại song
song hai loại lực ngược chiều là lực hút Coulomb và lực đẩy khoảng cách ngắn (short
range repulsion). Vì vậy ta có thể coi cấu trúc perovskite thực chất là mạng ba chiều
của các bát diện BO6, hay cũng có thể mô tả chúng gồm những hình lập phương xếp
chặt của ion A và Oxy với ion B ở tâm bát diện. Cấu trúc bát diện BO 6 đóng vai trò
quyết định đối với cấu trúc, tính chất điện từ và nhiều tính chất mang ý nghĩa ứng dụng
khác nhau của vật liệu perovskite.
(a)
(b)
7
Hình 1.1.Cấu trúc perovskite (ABO3) lập phương lý tưởng (a), và sự sắp
xếp các bát diện trong cấu trúc perovskite lập phương (b).
Một thông số rất quan trọng của cấu trúc perovskite cần xét để đánh giá sự ổn
định của liên kết giữa các ion A2+, B4+ và O2- đó là thừa số bền vững t:
√
trong một số công trình tính toán có thể coi CaMn
là perovskite lý tưởng không biến
dạng [17]có ô cơ sở giả lập phương (xem hình 1.3)
(b)
(a)
Hình 1.3:Ô đơn vị trực thoi của CaMnO3 (a) và ô đơn vị giả lập phương (b)
1.1.2. Cấu hình điện tử
Đặc trưng tinh thể quan trọng của cấu trúc perovskite ABO3 là sự tồn tại bát diện
BO6. Trên cơ sở cấu trúc bát diện BO6 và sự tương tác tĩnh điện giữa các ion kim loại
chuyển tiếp và ion O2-, chúng ta xét sự ảnh hưởng của trường tinh thể bát diện đến sự
sắp xếp các điện tử trên các mức năng lượng của lớp điện tử d của các ion kim loại
chuyển tiếp.
Như chúng ta đã biết, trong biểu thức Hamiltonian khi thế năng tương tác giữa
hai điện tử là nhỏ hơn nhiều so với động năng chuyển động của điện tử thì hàm thế
năng có thể coi như hàm nhiễu loạn. Trong các vật liệu manganite, thế năng nhiễu loạn
9
được sinh ra do trường tinh thể, sự phá vỡ đối xứng từ lập phương sang bát diện của
trường tinh thể. Kết quả là một quỹ đạo bị suy biến và tách ra thành những quỹ đạo có
mức năng lượng khác nhau.
Hình 1.4 là sơ đồ tách mức năng lượng của ion Mn3+ trong trường tinh thể của
perovskite cấu trúc lập phương. Đối với một nguyên tử tự do, các quỹ đạo có cùng số
Hình 1.4: Mô tả sự tách mức năng lượng của orbital d trong trường tinh thể bát
diện với năng lượng tách mức
(CF- crystal field: trường tinh thể), và tách mức do
méo mạng Jahn-Teller với năng lượng tách mức
(JT – Jahn – Teller)
Theo lý thuyết Jahn – Teller (JT), một phân tử có tính đối xứng cấu trúc cao với
các quỹ đạo điện tử suy biến sẽ phải biến dạng để loại bỏ suy biến, giảm tính đối xứng
và giảm năng lượng tự do. Xét trường hợp của các kim loại chuyển tiếp với cấu trúc
điện tử trên các quỹ đạo 3d không đầy. Hiệu ứng JT xảy ra trong một ion kim loại chứa
số lẻ điện tử trên mức eg (xem hình 1.5). Trong trường hợp của Mn3+ có cấu hình điện
tử 3d4, theo quy tắc Hund thì 4 điện tử này phải nằm trong 4 orbital có mức năng lượng
khác nhau từ thấp lên cao, do đó, sẽ có 3 điện tử nằm trên 3 quỹ đạo của mức t2g. Như
vậy còn 1 điện tử nằm ở mức năng lượng cao eg là mức suy biến bội hai, nên có hai
cách sắp xếp khả dĩ:
và
Nếu theo cách sắp xếp thứ nhất (
.
) thì lực hút tĩnh điện giữa ion O2- với
ion Mn3+ theo trục z sẽ yếu hơn so với trên mặt phẳng xy, điều này dẫn đến độ dài các
liên kết Mn – O không còn đồng nhất như trong trường hợp perovskite lí tưởng: ta sẽ
có 4 liên kết Mn – O ngắn trên mặt xy và 2 liên kết Mn – O dài hơn dọc theo trục z. Ta
12
Lý thuyết JT không cho biết xảy ra loại méo mạng nào, không tiên đoán được
cường độ biến dạng mà chỉ cho thấy méo mạng sẽ làm giảm năng lượng của hệ. Chính
vì thế các điện tử bị định xứ trong các ô mạng cơ sở.
Hiệu ứng JT có ý nghĩa quan trọng trong việc giải thích các bản chất điện và từ
của perovskite, đặc biệt là trạng thái trật tự điện tích trong các perovskite manganite có
pha tạp.
1.1.3. Các tƣơng tác vi mô dẫn tới tính chất từ trong hệ vật liệu perovskite
Trật tự từ trongperovskite được hình thành do sự cạnh tranh nhiều tương tác vi
mô như: siêu trao đổi phản sắt từ, trao đổi kép, trao đổi trực tiếp…Ta xét 2 tương tác
quan trọng nhất trong số đó.
a. Tƣơng tác siêu trao đổi (Super Exchange Interaction- SE)
Như ta đã thấy trong cấu trúc của hợp chất perovskite, các nguyên tử của kim
loại chuyển tiếp ở vị trí B có bán kính nhỏ bị ngăn cách bởi anion O2- ở giữa có bán
kính khá lớn (1.36A0), lớn hơn nguyên tử ở vị trí B nên không có sự xen phủ trực tiếp
giữa các cation kim loại chuyển tiếp. Vì thế cường độ tương tác trao đổi trực tiếp giữa
chúng là rất yếu. Lúc này, các cation kim loại chuyển tiếp chủ yếu tương tác với nhau
thông qua việc trao đổi điện tử với anion O2-. Quá trình trao đổi giữa anion O2- và
cation kim loại chuyển tiếp ở vị trí B là quá trình xen phủ giữa các đám mây điện tử
lai hóa trống eg của cation kim loại chuyển tiếp với đám mây điện tử được chiếm đầy
pб của anion O2- lân cận . Liên kết giữa các ion từ với ion oxy trung gian là liên kết
cộng hoá trị. Tương tác siêu trao đổi là tương tác giữa các spin của ion từ qua nguyên
tử ion oxy trung gian. Hợp chất mà loại tương tác siêu trao đổi chiếm ưu thế đều biểu
hiện tính dẫn điện môi và tương tác phản sắt từ chiếm ưu thế. Bởi vì các điện tử tham
gia liên kết không chuyển động được. Điều này cho biết vật liệu perovskite không pha
tạp và pha tạp hoàn toàn biểu hiện tính điện môi phản sắt từ.
13
14
Hình 1.6: Tương tác siêu trao đổi giữa các ion
Mn+3 qua ion oxy trung gian
b. Tƣơng tác trao đổi kép (Double Exchange Interaction - DE)
Với mô hình tương tác siêu trao đổi đã giải thích được tính chất từ của vật liệu
perovskite không pha tạp và pha tạp hoàn toàn nhưng mô hình này không thể giải thích
được tính chất của vật liệu khi pha tạp một phần. Ví dụ, ta xét hợp chất La1-xSrxMnO3
(0 x 1) : Sr có số oxi hóa là 2, La có số oxi hóa là 3. Khi x = 0 và x = 1 hợp chất này
là phản sắt từ, điện môi. Điều này đã được giải thích thông qua tương tác trao đổi giữa
2 cation Mn3+ gián tiếp qua anion O2-. Còn khi pha tạp hay x khác 0 và 1 thì vật liệu có
tỉ lệ giữa số ion Mn4+ và ion Mn3+ tương ứng là x và 1-x. Khi pha tạp Sr có số oxi hóa
là 2 vào thay thế cho La có số oxi hóa là 3 đã dẫn tới một phần Mn3+ chuyển thành
Mn4+để đảm bảo cân bằng điện tích. Và vật liệu có tính sắt từ ở Tc xung quanh nhiệt độ
phòng, dưới Tc vật liệu có tính kim loại . Để giải thích rõ cơ chế chuyển pha từ phản
15
sắt sang tính sắt từ mạnh, Zener [20] đã đưa ra mô hình cơ chế tương tác trao đổi kép
như sau:
-
Liên kết Hund nội nguyên tử là rất mạnh nên mặc dù cấu hình spin của
các ion thay đổi, do sự trao đổi đồng thời các điện tử của các ion lân cận, nhưng
spin của mỗi điện tử luôn song song với spin định xứ của ion.
-
Ca1-x Rx MnO3 ( R là các nguyên tố đất hiếm như La, Dy, Yb và Y, x~ 0 .1- 0.12) trong
khoảng từ nhiệt độ phòng tới 1000K cho thấy sự pha tạp khác nhau ảnh hưởng phong
phú đến cấu trúc điện tử và các đặc trưng nhiệt điện của vật liệu.
Với việc tăng mức độ pha tạp thì giá trị điện trở suất ρ lúc đầu giảm với
khoảng x~ 0.1-0.12,và sau đó thì tăng.
Điều này được giải thích như sau: sự giảm ρ lúc ban đầu là do sự tăng mật độ
điện tử do pha tạp điện tử bởi các ion đất hiếm hóa trị ba ( R3+ ) thay cho Ca
+2
. Theo
cơ sở cân bằng hóa trị, sự thay thế Ca bởi R3+ sẽ làm xuất hiện các ion Mn3+ , do đó
các điện tử egcó thể nhảy giữa Mn3+ và Mn4+ dễ dàng hơn. Vì vậy điện trở suất ρ giảm.
Tuy nhiên,khi mật độ electron vượt quá một giá trị nhất định(tương ứng với
x~0.1-0.12) trong mạng tinh thể ion có thể xảy ra hiện tượng trật tự điện tích (charge
ordering). Trật tự điện tích này dẫn đến sự giảm độ linh động và tăng sự định xứ của
các electron,do vậy điện trở suất ρ tăng.
Ta thấy rằng sự phụ thuộc của điện trở suất ρ vào mức độ pha tạp không phải là
đơn điệu. Hình 1.8 mô tả kết quả thực nghiệm thu được trong [19 ]
17
Hình 1.8: Sự phụ thuộc của điện trở suất ρ, hệ số Seebeck S , hệ số công suất P
vào nhiệt độ của mẫu Ca1-x Rx MnO3 (R: La, Dy, Yb và Y).
b. Ngoài ra, việc phân tích cấu trúc điện tử và đặc trưng nhiệt điện của hệ
Hình 1.10: Mật độ trạng thái (DOS) của hệ CaMnO3 (trái) và hệ CaMnO3 pha
tạp Sr (phải).
c.
Ngoài ra, một số nghiên cứu về cấu trúc điện tử của CaMnO3 cho thấy
khi pha tạp một số ion kim loại như Eu, Yb hay Y thì thể tích ô cơ sở tăng lên so với
khi không pha tạp, khoảng cách Mn-O tăng lên và góc liên kết Mn-O-Mn giảm đi [9]
[11] .
20
CHƢƠNG 2: LÝ THUYẾT CẤU TRÚC ĐIỆN TỬ CHẤT RẮN VÀ
PHƢƠNG PHÁP LÝ THUYẾT PHIẾM HÀM MẬT ĐỘ
2.1. Phƣơng pháp lý thuyết cấu trúc điện tử chất rắn
Phương pháp chính để nghiên cứu cấu trúc vật liệu đó là phương pháp lý thuyết
cấu trúc điện tử (electron structure theory methods). Phương pháp này thực hiện những
tính toán cơ bản sau:
- Tính toán năng lượng của một cấu trúc phân tử nào đó với một cách bố trí
không gian nhất định các nguyên tử hoặc hạt nhân và các điện tử (thường được gọi là
tính toán năng lượng một điểm – single-point energy calculation). Từ đó, các tính chất
liên quan có thể được dự đoán thông qua một vài phương pháp khác nhau.
- Thực hiện tối ưu hóa cấu trúc (Geometry optimization) nhằm tìm ra trạng thái
cấu trúc có năng lượng nhỏ nhất trong lân cận gần nhất của cấu trúc lý thuyết ban đầu.
Sự tối ưu hóa cấu trúc phụ thuộc chủ yếu vào Gradient của năng lượng (Đạo hàm bậc
nhất của năng lượng theo vị trí của các nguyên tử).
- Tính toán tần số dao động (Frequency calculation) của các phân tử do sự chuyển
động của các nguyên tử trong phân tử. Các tần số này phụ thuộc vào đạo hàm bậc hai
của năng lượng có liên quan tới cấu trúc nguyên tử.
cũng như các thông số thực nghiệm. Thay vào đó, những tính toán ab initio dựa hoàn
toàn vào những định luật cơ bản của cơ học lượng tử và các hằng số vật lý như:
- Vận tốc của ánh sáng c~3.108 m/s;
- Khối lượng me = 9,1.10-31kg và điện tích qe = -1,602.10-19 C của điện tử;
- Hằng số Planck h = 6,626.10-34 J.s;
Trong quá trình thực hiện luận văn này, chúng tôi sử dụng phương pháp abinitio làm phương pháp nghiên cứu khoa học do vậy phần dưới đây sẽ hoàn toàn tập
trung giới thiệu về nhóm phương pháp này.
2.2. Phƣơng pháp phiếm hàm mật độ (DFT).
Lý thuyết phiếm hàm mật độ ( Density Functional Theory) là một lý thuyết
được dùng để mô tả các tính chất của hệ electron trong nguyên tử, phân tử, vật rắn,...
trong khuôn khổ của lý thuyết lượng tử. Trong lý thuyết này, mật độ điện tử được sử
dụng để mô tả trạng thái của hệ thay vì sử dụng hàm sóng.Một hàm sóng mô tả hệ N
điện tử sẽ phải chứa 3N biến tọa độ.Trong khi đó, mật độ điện tử chỉ phụ thuộc vào ba
biến tọa độđộc lập với số điện tử. Vì thế khi gia tăng số điện tử của hệ, hàm sóng sẽ trở
nên phức tạp nhưng mật độ điện tử không thay đổi số biến. Do vậy, lý thuyết phiếm
hàm mật độ có ưu điểm lớn (và hiện nay đang được sử dụng nhiều nhất) trong việc tính
22
toán các tính chất vật lý cho các hệ cụ thể xuất phát từ những phương trình rất cơ bản
của vật lý lượng tử.
2.2.1. Tính chất lƣợng tử của chất rắn
Các nguyên tử hay phân tử liên kết chặt chẽ với nhau và sắp xếp một cách có trật
tự, tuần hoàn trong vật rắn kết tinh. Trong mỗi nguyên tử thì khối lượng của hạt nhân
là chủ yếu, ta bỏ qua khối lượng của các điện tử.
Theo quan điểm lượng tử, Hamiltonian của chất rắn được viết như sau:
̂
trong đó:
Đây là mô hình DFT đầu tiên, trong đó các hàm sóng phức tạp của hệ N electron
bằng một đại lượng mới đơn giản hơn trong việc giải phương trình Schrodinger, đó là
mật độ electron n(r) .Trong phương pháp Thomas-Fermi, động năng của hệ điện tử
được xấp xỉ như là một phiếm hàm của mật độ điện tử và hệ điện tử được lý tưởng hóa
như là các hệ khí đồng nhất gồm các hạt không tương tác. Khi đó mật độ điện tử bằng
với mật độ địa phương tại mọi điểm cho trước.
23
Hạn chế của phương pháp này đó là không tính tới sự trao đổi và tương quan của
các điện tử. Vào năm 1930, Dirac đề ra phương pháp gần đúng địa phương cho sự trao
đổi và vẫn được sử dụng đến ngày nay, khắc phục những hạn chế của Thomas - Fermi.
Phiếm hàm năng lượng cho các điện tử trong một trường thế ngoài Vext(r) có dạng:
1
n(r )n(r ')
(2.3)
EFT [n] C1 d 3rn(r )5/3 d 3rVext (r )n(r ) C2 d 3rn(r )(4/3) C1 d 3rd 3r '
2
r r'
Trong đó thành phần đầu tiên là gần đúng địa phương cho động năng với
C1
3
3 3
(3 2 )2/3 , thành phần thứ 3 là trao đổi địa phương C2 ( )1/3 và thành phần
10
4
-
n(r) là một đại lượng có thể quan sát và có thể đo bằng thực nghiệm, ví dụ như
bằng nhiễu xạ tia X.
-
n(r) sẽ giảm theo tiệm cậm của hàm mũ khi ra xa tất cả các hạt nhân
√
với I là năng lượng ion hóa.
Các định lý Hohenberg – Kohn:
Năm 1964, Hohenberg và Kohn đã làm việc cùng nhau để nghiên cứu khắc phục
các vấn đề cơ bản của mẫu Thomas-Fermi. Họ đã đưa ra và chứng minh hai định lí
quan trọng.
Hình 2.1: Sơ đồ minh hoạ cho định lý Hohenberg-Kohn. Các mũi tên ngắn biểu
thị giải pháp thông thường là giải pháp Schrodinger mà thế Vext(r) xác định tất cả
trạng thái ψi(r), bao gồm trạng thái cơ bản ψ0(r) và mật độ trạng thái cơ bản n0(r).
Các mũi tên dài có kí hiệu “HK” chỉ định lý Hohenberg-Kohn.Chúng liên hệ với nhau
thành một vòng khép kín.
25
Cách tiếp cận của Hohenberg và Kohn được sử dụng để phát biểu lý thuyết phiếm
hàm mật độ như là một lý thuyết chính xác cho hệ nhiều hạt. Phát biểu này được áp
dụng cho mọi hệ hạt tương tác trong một thế ngoài Vext bao gồm tất cả các vấn đề của
điện tử và các nguyên tử đứng yên. Trong đó Hamiltonian được viết dưới dạng:
⟨
⟩
⟨
⟩
⟨
⟩
∫
(2.6)
Tương tự:
⟨
⟩
⟨
|
|
⟩
⟨