Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS.. Ảnh hưởng của hiệu ứng phi điều hoà và áp suất lên các đại lượng nhiệt động của các tinh thể có cấu trúc lập phương trong lý thuyết EXAFS..
Bức xạ tia X
Tia X được R¨ontgen phát hiện năm 1895 đã góp phần quan trọng trong việc nghiên cứu cấu trúc của vật liệu Người ta thường sử dụng bức xạ tia X đóng vai trò là nguồn photon trong các tương tác với vật rắn Kết quả thu được của tương tác này là các phổ có chứa thông tin về cấu trúc của vật rắn [2, 69, 75].
Trong thí nghiệm đối với ống tia X, khi các electron được tạo ra bởi sợi dây Wolfram được nung nóng và chuyển động rất nhanh từ cathode đến anode nhờ một hiệu điện thế cao giữa hai cực trong chân không, do có năng lượng lớn nên các electron xuyên sâu vào lớp vật chất của kim loại chắn và tạo ra bức xạ điện từ có bước sóng rất ngắn gọi là tia
X Theo quy tắc Duane-Hunt, nếu bỏ qua động năng ban đầu và tất cả năng lượng mà electron thu được chuyển thành năng lượng photon của tia X thì bước sóng cực tiểu của tia X sẽ được xác định bởi: c λmin= h ε
, (ε = eV ) (1.1) eV với h là hằng số Planck, c là tốc độ ánh sáng trong chân không, e là điện tích của một electron và λmin là bước sóng cực tiểu ứng với năng lượng cao nhất của tia X có thể nhận được đối với điện thế tăng tốc V Phổ tia X gồm phần liên tục và phần đặc trưng Khi các electron va chạm vào kim loại chắn, chúng truyền năng lượng cho tấm chắn nên trên bề mặt tấm chắn có một điện từ trường biến thiên rất nhanh và tạo ra một sóng điện từ có bước sóng rất ngắn, kết quả là cho một phổ bức xạ hãm liên tục (Bremsstrahlung spectrum) hay phổ tia X liên tục (Hình 1.1) Khi điện thế
Hình 1.1: Phổ liên tục của tia X và phổ bức xạ đặc trưng tăng lên, thì bước sóng λmin giảm và cường độ toàn phần sẽ tăng.
Sự phát sinh tia X đặc trưng có liên quan với sự dịch chuyển electron của các nguyên tử giữa các vùng năng lượng Khi điện thế tăng đến một giá trị nhất định, chùm electron có khả năng xuyên sâu vào vật liệu và va chạm với các nguyên tử làm electron của nguyên tử bị bật ra khỏi lớp vỏ nguyên tử tạo ra ở đó một lỗ trống Theo nguyên lý cực tiểu năng lượng electron từ các mức năng lượng cao hơn trong vùng dẫn nhảy xuống lấp đầy các lỗ trống này và phát ra các bức xạ đặc trưng Trên bảng quang phổ tạo ra một số vạch rõ, gián đoạn, mô tả bằng các đường đặc trưng chồng lên phổ bức xạ hãm liên tục và lớn hơn cường độ của bức xạ hãm cỡ 10 3 lần và được biểu diễn trên Hình (1.1) Cường độ này phụ thuộc vào hai mức năng lượng nguyên tử tham gia vào chuyển dịch Thí dụ các tia K α là do các electon nhảy từ mức năng lượng thứ hai (L) về mức năng lượng thứ nhất (K) phát ra, các tia K β là do sự chuyển dịch của các electron giữa lớp M và K Chúng là các bức xạ đơn sắc và gián đoạn Người ta có thể tạo ra bức xạ này qua sử dụng một thế tăng tốc, các electron trong lượng để làm bật các electron từ trong
Hình 1.2: Tương tác của electron với mô hình đơn giản của nguyên tử. nguyên tử và tạo ra lỗ trống Hình (1.2) mô tả các quá trình vật lý khi dòng electron được phóng qua một nguyên tử [2, 69].
Như vậy, ống tia X có khả năng tạo ra cả phổ tia X liên tục (bức xạ hãm hay bức xạ trắng) và phổ tia X gián đoạn (bức xạ đặc trưng hay bức xạ đơn sắc) Các bức xạ tia X đặc trưng được dùng rộng rãi trong các nghiên cứu nhiễu xạ tia X, các bức xạ tia X liên tục được dùng trong XAFS.
Bức xạ Synchrotron
Bức xạ Synchrotron (còn được gọi là bức xạ Magnetobremsstrahlung) đạt được khi các hạt mang điện như electron hay positron chuyển động với tốc độ gần bằng tốc độ ánh sáng theo một quỹ đạo xoắn ốc trong từ trường và bắn vào vật thử Bức xạ Synchrotron có một số đặc tính cơ bản quan trọng là [75]: (1) Cường độ lớn trong vùng năng lượng rộng, liên tục; (2) Cường độ và vị trí nguồn có độ ổn định cao; (3) Có tính chuẩn trực lớn; (4) Phân cực phẳng, môi trường bức xạ sạch; (5) Cấu trúc thời gian theo xung chuẩn xác, các xung theo micro – giây; (6) Kích thước nguồn nhỏ được xác định qua
Hình 1.3: Hệ số hấp thụ ϵ(E) có phần cấu trúc tinh tế của hấp thụ tia X
Khi chiếu một chùm photon tia X vào vật rắn thì sẽ xảy ra hai quá trình là tán xạ và hấp thụ Quá trình tán xạ là do photon tia X bị phản xạ trở lại sau khi va chạm với electron lõi hoặc nguyên tử bao gồm tán xạ đàn hồi (tán xạ Rayleigh – do photon tia X va chạm hoàn toàn đàn hồi với electron, sau tán xạ bước sóng tia X không thay đổi) và tán xạ không đàn hồi (tán xạ Compton − do photon tia X va chạm với các electron hoá trị và bước sóng tia X bị thay đổi) Quá trình hấp thụ liên quan đến hiệu ứng quang điện là do các electron lõi hấp thụ photon tia
X và chuyển lên mức cao hơn hoặc bắn ra ngoài nguyên tử Nếu electron bắn ra ngoài nguyên tử thì ta có phổ electron quang PES
(Photo Electron Spectrocopy), còn nếu electron quang ở lại trong vật rắn sau khi tán xạ với các nguyên tử lân cận rồi trở lại giao thoa với sóng của quang electron được phát ra từ nguyên tử hấp thụ thì ta thu được phần cấu trúc tinh tế phổ hấp thụ tia X (XAFS). Σ Σ
Lý thuyết về phổ cấu trúc tinh tế hấp thụ tia X
Khi cho một chùm bức xạ điện từ có cường độ I0 đi qua một lớp vật chất, chùm bức xạ sẽ bị hấp thụ một phần, cường độ chùm bức xạ sau khi đi qua là: do đó
I Trong hệ thức (1.2), ϵ là hệ số hấp thụ, x là bề dày lớp vật chất Nếu chùm ánh sáng tới là tia X thì sau cận hấp thụ với năng lượng photon là
ℏω sẽ xuất hiện phổ cấu trúc tinh tế hấp thụ tia X, Hình (1.3).
Như vậy, hệ số hấp thụ ϵ trong trường hợp XAFS bao gồm hai thành phần, đó là hệ số hấp thụ của một nguyên tử biệt lập và phần cấu trúc tinh tế X của phổ tia X, vì thế hệ số hấp thụ toàn phần sẽ là ϵ(E)
= ϵ a (E) 1 + χ(E) Từ đó, phần cấu trúc tinh tế của tia X hay phổ XAFS được viết dưới dạng [2, 62, 63, 68]; χ(E) = ϵ(E) − ϵa(E)
(1.4) ϵ a (E) Để xác định χ(E) trong (1.4) ta phải xác định hệ số hấp thụ ϵ(E) khi vật thể tương tác với sóng điện từ đặc trưng bởi thế vecto A⃗ Hệ số này được tính toán theo quy tắc lọc lựa và có dạng đặc trưng trong gần đúng lưỡng cực [25, 69]. ϵ = 2π Σ
Trong (1.5), các hàm sóng |i⟩ và |f⟩ là hàm riêng của Hamiltonian hiệu dụng ở trạng thái đầu và trạng thái cuối ứng với các mức năng lượng E i và E f Như vậy electron đã hấp thụ năng lượng ℏω của photon tia X với phân cực e và chuyển từ trạng thái đầu |i⟩ với năng lượng E i đến trạng thái cuối |f⟩ với năng lượng E f Khi đó, do tính chất f đối xứng của hàm sóng
− i mà các yếu tố của ma trận dịch chuyển đối với các số lượng tử của trạng thái đầu (l i , m i ) và trạng thái cuối (l f , m f ) sẽ tuân theo quy tắc lọc lựa là l f = l i ± 1, m f = m i , m i = ±1 Từ đây, ta dễ dàng xác định các số lượng tử của trạng thái cuối |f⟩ vào trạng thái đầu |i⟩ mà thu được các cận hấp thụ khác nhau Biến đổi tiếp hệ thức (1.5) và biểu diễn nó qua ma trận mật độ n hay hàm Green G của toàn hệ [42]. ϵ = 2π Σ
E π + ℏω). Đối với cận hấp thụ nhất định thì trạng thái |i⟩ bao giờ cũng được biết trước, cho nờn để đỏnh giỏ hệ số hấp thụ à hay phổ cấu trỳc tinh tế χ, người ta chỉ cần xây dựng trạng thái cuối |f⟩ hay hàm Green G của toàn hệ.
Quang phổ XAFS nói chung xuất hiện trong khoảng 40 − 1000 eV sau cận hấp thụ Trong quá trình nghiên cứu người ta còn phân ra nhiều lý thuyết XAFS khác nhau như EXAFS (Extended XAFS) khi động năng của quang electron E > 50 eV; XANES (X - Ray Absorption Near- Edge Structure) và NEXAFS (Near-Edge XAFS) khi E < 50 eV, tức là cấu trúc gần cận hấp thụ Ngoài ra, đối với vùng mặt tinh thể còn có các lý thuyết như SXANES (Surface XANES) và SEXAFS (Surface EXAFS) [74]. Trong lịch sử nghiên cứu về phương pháp XAFS đã tồn tại hai cách lý luận là mức độ xa (LRO: Long Range Order) và mức độ gần (Short Range Order) [2, 31, 71] Đối với lý luận mức độ xa LRO, các phổ XAFS được đặc trưng bởi mật độ trạng thái của trạng thái cuối, nó được xác định qua cấu trúc vùng năng lượng, quãng đường chuyển động tự do i i xác suất chuyển dịch bị bỏ qua Đối với mức độ gần SRO các phổ XAFS được đặc trưng qua trạng thái cuối, nó bao gồm các hiệu ứng tán xạ bởi các nguyên tử lân cận và tán xạ Σ S
N j Σ S 2 ngược trở lại nguyên tử hấp thụ ban đầu, thời gian sống của quang electron cũng như lỗ trống ở tâm lõi do quang electron để lại được tính qua quãng đường dịch chuyển tự do, các hiệu ứng dao động nhiệt của các nguyên tử được tính qua hệ số Debye – Waller (DWF).
Các lý thuyết LRO và SRO cho các tiên đoán giống nhau về các phổ
XAFS và sự phụ thuộc của chúng vào nhiệt độ vì mật độ trạng thái của trạng thái cuối cũng xuất hiện qua tán xạ của các electron bởi các nguyên tử lân cận Tuy nhiên các phát triển của phương pháp XAFS, lý thuyết SRO có nhiều ưu điểm do việc chuyển hàm Fourier các phổ
XAFS để nhận được các thông tin về cấu trúc nguyên tử của vật rắn.
Ngoài ra, khi tính toán các phổ XAFS người ta sử dụng các tham số của nguyên tử và vật rắn, cho nên khi so sánh các phổ lý thuyết với phổ thực nghiệm thì người ta sẽ nhận thông tin về các tham số này.
Như vậy, trong lý thuyết về quang phổ XAFS hiện đại, phổ cấu trúc tinh tế hấp thụ tia X được coi là hiệu ứng của trạng thái cuối Sóng của electron quang mà nguyên tử phát ra khi hấp thụ photon tia X sẽ bị tán xạ bởi các nguyên tử lân cận rồi quay trở lại nguyên tử hấp thụ Trạng thái cuối là kết quả giao thoa của sóng electron quang bị tán xạ và sóng phát ra ban đầu, vì vậy mà nó chứa thông tin về vị trí của các nguyên tử lân cận.
Phổ XAFS cận K đối với chất đa tinh thể, theo các phương trình (1.4) và (1.6) thường có dạng: χ(k) = o 2 j k
Trong gần đúng điều hòa ở nhiệt độ thấp, khi R i = ⟨r j ⟩ phương trình (1.7) trở thành: χ(k)
(1.8) trong đó, F j (k) là biên độ tán xạ ngược của mỗi nguyên tử lân cận, N j là số nguyên tử lân cận trên lớp nguyên tử thứ j, S 2 là hệ số đặc trưng cho hiệu
H A ứng nhiều hạt, δ(k) là độ dịch pha trong tán xạ, r j là bán kính lớp nguyên tử thứ j, k là số sóng có giá trị được xác định từ hệ thức k .2m (E − E0). Tổng độ dịch pha δ k được xác định theo hệ thức: δ k = 2δ ′ (k) + Ψ, (1.9) trong đó δ ′ (k) đặc trưng cho độ dịch pha của electron quang lúc phát ra ngoài nguyên tử và bằng độ dịch pha của nó so với lúc phản xạ lại nguyên tử ban đầu, còn Ψ là độ dịch pha khi electron quang tán xạ trên nguyên tử lân cận thứ j.
Trong phương trình (1.8), σ 2 là độ dịch chuyển tương đối trung bình bình phương (MSRD: mean-square relative displacement) của khoảng cách giữa hai nguyên tử, nó đóng góp chính vào hệ số Debye-Waller e −2σ 2k2
Ảnh Fourier của phổ cấu trúc tinh tế hấp thụ tia X
Các nhà khoa học nghiên cứu về phổ EXAFS đã phát hiện ra rằng các đỉnh trong ảnh Fourier của phổ EXAFS tương ứng với bán kính của các lớp nguyên tử, vì vậy ảnh Fourier của phổ EXAFS cho thông tin về cấu trúc của vật thể Ở nhiệt độ thấp, bán kính lớp nguyên tử không phụ thuộc nhiệt độ
Hình 1.4: Ảnh Fourier của phổ XAFS đối với tinh thể đồng [47]
(mô hình dao động điều hòa) Tuy nhiên thực nghiệm đã cho thấy khi nhiệt độ tăng, các thông tin về cấu trúc có những sai số đáng kể do ảnh hưởng của dao động phi điều hòa của các nguyên tử [11, 47] Ví dụ trong Hình (1.4) cho thấy ở các nhiệt độ khác nhau ảnh Fourier của các phổ EXAFS cho các thông tin về cấu trúc khác nhau Đỉnh thứ nhất trong ảnh Fourier của phổ XAFS cho thông tin về khoảng cách từ nguyên tử hấp thụ đến nguyên tử tán xạ thuộc lớp thứ nhất, nhưng ở ba nhiệt độ khác nhau (297 0 K, 703 0 K, 973 0 K) ta nhận được ba thông tin về cấu trúc khác nhau.
Theo phương trình (1.8), cấu trúc tinh tế của phổ EXAFS được đặc trưng chủ yếu qua hàm sine, nên ta có thể chuyển hàm EXAFS với biến số là số sóng k trở thành hàm có biến số là toạ độ r thông qua hàm chuyển Fourier như sau [2, 33, 46]
Trong (1.10) ta nhận dược thông tin về toạ độ R = ⟨r⟩, tức là xác định
∫ năng lượng để đánh giá hệ thức (1.10) là rất quan trọng Khi ta mô tả electron được kích thích ở ngoài mặt cầu muffin-tin qua sóng với số sóng bằng k thì năng
∫ ∫ kr 2 j lượng E ∼ k 2 tính từ điểm không của muffin-tin nằm ở cỡ đáy vùng hoá trị, nghĩa là cỡ 10 eV dưới cận hấp thụ Để chuyển Fourier ta cần phải biết sự phụ thuộc của biến số của hàm sine vào số sóng k Sử dụng sự phụ thuộc tuyến tính của pha dao động vào số sóng k dưới dạng:
Thay (1.11) vào các phương trình (1.7), (1.10), thay k bằng k là giá trị trung bình của nó (trừ trong hàm sine) ta nhận được [2, 63];
Như vậy từ các đỉnh của phổ EXAFS được xác định qua phương trình (1.12) ta biết được giá trị của (r j + a), cho nên nếu biết a trong hệ thức (1.11) thì xác định được r j , tức là xác định được cấu trúc nguyên tử của vật rắn.
Khi nhiệt độ tăng, do ảnh hưởng của dao động phi điều hòa nên thông tin về cấu trúc của vật thể sẽ bị thay đổi đáng kể Vì vậy ta cần phải tính đến đóng góp của các nhiễu loạn phi điều hòa tác động lên phổ EXAFS, từ đó sẽ xác định được chính xác cấu trúc của vật thể.
Hàm phân bố hiệu dụng
Trong phổ XAFS, phép khai triển cumulant thường được thực hiện qua hàm phân bố hiệu dụng [12] được định nghĩa dưới dạng:
P(r, γ) ρ1(r)e −2γ , (1.13) r2 với ρ1(k) là phân bố xác xuất của các nguyên tử trên lớp, điều kiện chuẩn hoá cho ta ρ1(r)dr = 1 và ρ1(r) có liên quan đến phân bố ba chiều ρ(r)d 3 r = 1, bởi ρ(r) = 4πr 2 ⟨ρ(r)⟩Ω, ở đây ⟨⟩Ω biểu diễn trung bình góc 4π, ρ1(r) = 0 khi r < 0 Chuyển Fourier, hệ thức (1.13) có dạng: j j
∫ Σ với hàm γ ≡ λ −1 và hàm r là tham số được chọn sau Trong lý thuyết EXAFS, phổ EXAFS cận K thường được viết theo hệ thức của các hàm phân bố dưới dạng: χ(k) = NF(k) ρ1(r)e −2r/λ(k) sin [2kr + δ(k)] dr, (1.15) r2 trong đó N là số nguyên tử trên một lớp nguyên tử, F(k), δ(k) là biên độ và pha của tán xạ, bao gồm tất cả các đóng góp từ nguyên tử hấp thụ, λ(k) là quãng đường tự do trung bình của electron quang và phụ thuộc vào số sóng k Từ các hàm phân bố và hàm chuyển Fourier trong các hệ thức (1.13), (1.14) phương trình của phổ EXAFS cận K có thể viết gọn thành: χ(k) = NF(k)Im Σ e(i(2kr+δ(k))P¯(r¯, γ, k) Σ
Các phương trình (1.15) hoặc (1.16) có thể viết chung theo dạng tổng quát: χ(k) ≡ A(k)sinϕ(k), trong đó biên độ thực của dao động XAFS là:
, (1.17) và pha của dao động là: ϕ(k) = 2kr + δ(k) + argP(r, γ, k) (1.18)
Các biểu thức (1.17), (1.18) tương ứng với các công thức về biên độ và pha của phổ quang electron thu được qua phép lọc Fourier.
Hàm phân bố P(r, γ, k) trong các hệ thức có thể được khai triển theo các mômen dịch chuyển của phân bố hiệu dụng:
Trong các công thức trên, P n là hàm của r và γ Tại các giá trị nhỏ của k thì chỉ có các mômen bậc thấp là quan trọng, nhưng khi k tăng lên, các mômen bậc cao hơn sẽ được lấy theo tất cả các bậc đóng góp Bản chất của khai triển trên là khai triển theo luỹ thừa của 2k∆r, với ∆r là bề rộng đặc trưng trong nhiễu loạn của phân bố. Áp dụng hàm phân bố trên để khai triển các cumulant, người ta thường thực hiện qua giá trị trung bình theo mỗi phân bố của biến x của hệ thức [12],
, (n ≥ 0) (1.21) Ở đây giá trị trung bình ⟨⟩ sẽ triệt tiêu một cách thích hợp rất nhanh ở vô cực Hiển nhiên nếu phân bố là chuẩn hoá Các cumulant sẽ được xác định bởi hệ thức tương quan giữa hàm phân bố hiệu dụng và giá trị trung bình của phân bố:
Khai triển hệ thức trên theo chuỗi Taylor và tách các cumulant bậc chẵn ta sẽ thu được các hệ thức về biên độ dao động:
(2k) 2n σ 2n , (1.23) (2n)! n=0 và các cumulant bậc lẻ sẽ mô tả pha của dao động nguyên tử:
Vì biên độ và pha dao động phụ thuộc vào số sóng k nên chúng ta không thể tùy chọn r, các cumulant σ (n) với (n ≥ 0) là không phụ thuộc vào điểm gốc Điều này cũng được mô tả từ phương trình (1.22) với sự liên quan tới r cùng với việc sử dụng sự phụ thuộc tuyến tính vào luỹ thừa của k Các cumulant bằng hoặc bé hơn các mômen luỹ thừa, nếu r = 0 chúng ta có thể thu được dP n /dq = P n+1 và dσ (n) /dq = σ (n+1) với n ≥ 0 và q ≡ −2γ Theo các phương trình(1.19), (1.20) và (1.22) thì σ (0) (γ) = lnP0(γ) Kết hợp các
∫ hệ thức, ta có thể viết công thức khai triển của các cumulant theo các dạng sau: tương tự ta có: σ(1) σd(0) dq
Nếu ta chọn r từ tâm của phân bố ta sẽ thu được các công thức cumulant dạng rút gọn: σ (0) = lnp0, σ (1) = 0, σ (2) = p2, σ (3) = p3, (1.30) σ (4) = p4 − 3p 2 , σ (5) = p5 − 10p3p2, σ (6) = p6 − 15p4p2 − 10p 2 + 30p 3 ,
Hệ số Debye - Waller
Trong gần đúng dao động điều hòa, hệ số Debye - Waller được mô tả đơn giản qua thừa số exp(−2k 2 σ 2 ) trong phương trình (1.12) Sự giảm dần của quang phổ EXAFS được mô tả qua hàm χ(k) dạng (1.8).
Khi nhiệt độ cao, nhiễu loạn lớn thì χ(k) sẽ được mô tả bởi phương trình tổng quát có dạng [12]: χ(k)
/λ kr 2 Σ sin [2krj + δj(k)] drj, j j j j j trong đó P(r j )dr j là xác xuất tìm thấy nguyên tử thứ j trong vùng từ r j tới (r j + dr j ) Hệ số Debye-Waller được xác định từ việc lấy trung bình công thức EXAFS tán xạ đơn trong hệ nhiều hạt với các cặp nguyên tử lân cận gần nhất với hàm phân bố cặp P(r) Nếu các hệ số khác trong hàm sine của (1.31) có tổng nhận được là hàm dao động nhỏ với r j thì kết quả chính sẽ có dạng:
Thay thế P(r j ) bằng hàm phân bố hiệu dụng P(r j , γ), hàm này kết hợp với biên độ của phổ EXAFS qua hệ thức của hàm phân bố cặp
P(r j , γ) = P(r j )e −2γr j /r 2 , trong đó P(r j ) là phân bố cặp và γ là nghịch đảo của quãng đường tự do trung bình Nếu nhiễu loạn là nhỏ hay có tính đối xứng Gauss, thì chúng ta có thể sử dụng gần đúng
− 2k σ 2 ), (1.33) với σ 2 là trung bình bình phương độ dài liên kết của dao động: σ 2 =
(1.34) Đại lượng σ 2 bao gồm hai yếu tố σ 2 (T ) sinh ra do dao động nhiệt và σ 2 (S) j j j sinh ra do nhiễu loạn cấu trúc và không phụ thuộc nhiệt độ.
Nếu coi gốc toạ độ đặt tại nguyên tử hấp thụ và nguyên tử lân cận tại vị trí j, →−u0 là vectơ dịch chuyển của nguyên tử hấp thụ, thì độ dịch chuyển xuyên tâm bậc 1 của nguyên tử lớp thứ j sẽ là δr j
= (→−u j − →−u 0).→−r j Vì thế trung bình bình phương của biên độ dao động được biểu diễn bằng độ dịch chuyển tương đối trung bình bình phương: σ 2 .Σ (→−u j −
Như vậy, trong EXAFS, σ 2 phụ thuộc vào lớp nguyên tử và chứa độ dịch chuyển trung bình bình phương (MSD - Mean Square Displacement) u 2 và j j Σ Σ Σ Σ(n) (n+m)
. Σ hàm dịch chuyển tương quan (DCF - Displacement Corelation
Biên độ và pha của phổ cấu trúc tinh tế hấp thụ tia X mở rộng
Tính đến sự phụ thuộc vào số sóng của quãng đường tự do trung bình thì trong (1.16) ta thấy các cumulant sẽ biến thiên theo k qua hệ thức P[r, γ(k)] Ở đây coi sự phụ thuộc vào k như một nhiễu loạn, với γ(k) ≡ γ0 + δγ(k) trong đó γ0 là giá trị của γ(k) tại một vài điểm thích hợp và sự phụ thuộc k của hấp thụ trong δγ(k), chú ý khi →−r = 0 thì hệ thức (1.15) sẽ trở thành: exp
Khai triển σ (n) theo γ0 bằng hệ thức
, (1.39) ta sẽ thu được hệ thức biên độ và pha biểu diễn theo tổng các cumulant:Σ ∞ ∞ (2ik) n
Lấy các số hạng bậc một theo δγ(k) ta có:
A(k) lnNF(k) = σ (0) − 2δγ(k)σ (1) − 2k 2 σ (2) + 4k 2 δγ(k)σ 3 , (1.42) ϕ(k) − δ(k) = 2kσ (1) − kδγ(k)σ (2) − 4 k 3 σ (3) + (1.43)Các hệ thức (1.42), (1.43) là các biểu thức tính biên độ và pha của phổ EXAFS qua tổng các cumulant σ (n) tại γ0.
Thế tương tác trong phổ cấu trúc tinh tế hấp thụ tia X mở rộng
Thế năng tương tác có vai trò chính trong năng lượng kết nối giữa các nguyên tử để tạo thành vật rắn Tương tác giữa các nguyên tử theo từng cặp liên kết thường được mô tả qua một số thế tương tác như thế Lennard- Jones, thế Madelung hoặc thế Morse Thế Lennard- Jones thường được sử dụng trong liên kết van-der-Waals và phổ biến đối vật rắn khí trơ Do có cấu trúc của các lớp electron lấp đầy có đối xứng cầu của các nguyên tử khí trơ rất bền vững, ít bị ảnh hưởng khi chúng kết hợp để tạo thành vật rắn, năng lượng tương tác giữa hai nguyên tử chỉ phụ thuộc vào khoảng cách giữa chúng và thường được biểu diễn qua thế Lennard-Jones Thế Madelung thường được dùng khi đánh giá thế tương tác giữa các nguyên tử của các tinh thể ion, thế tương tác này bao gồm thế đẩy giữa các đám mây electron và thế hút Coulomb giữa các ion dương và âm Trong phạm vi nghiên cứu, luận án này sử dụng thế cặp phi điều hòa Morse [17, 48] và xét gần đúng cho các tinh thể có cấu trúc lập phương Thế phi điều hòa Morse có dạng:
, (1.44) trong đó α có thứ nguyên nghịch đảo của khoảng cách ứng với độ rộng của thế (A −1 ) còn D có thứ nguyên của năng lượng (eV) và bằng năng lượng phân ly vì D = U(r0) Vị trí cân bằng của hai nguyên tử được biểu diễn qua r0, như trên (1.5) Khác với các thế khác, thế Morse cho ta độ dịch chuyển của nguyên tử đối với vị trí cân bằng r0 thông qua hiệu r − r0 = x Ta có thể viết biểu thức của thế Morse theo dạng của x:
U(x) = D(e −2αx − 2e −αx ) (1.45) Khai triển (1.45) quanh giá trị nhỏ của x ta có:
Hình 1.5: Hệ số giãn nở nhiệt mạng a mô tả sự bất đối xứng của thế Σ
−2(1 − αx + 1 α 2 x 2 − 1 α 3 x 3 + 1 α 4 x 4 + ) Σ Đối với quang phổ EXAFS ở các bậc cao nhiễu loạn thường là yếu và đóng góp của chúng là không đáng kể Lấy gần đúng tới số hạng bậc
4 khi đó hệ thức của U(x) sẽ có dạng:
Như vậy biểu thức của thế Morse theo độ lệch tức thời x của các nguyên tử sẽ được viết gọn thành
Kết luận Chương 1 23 Chương 2 LÝ THUYẾT VỀ PHỔ CẤU TRÚC TINH TẾ HẤP THỤ
Như vậy, trong nghiên cứu phổ EXAFS chúng ta có thể xác định được cấu trúc của tinh thể Việc quan trọng là chúng ta phải xác định được thế năng tương tác.
MỞ RỘNG PHI ĐIỀU HÒA
Phổ EXAFS phi điều hòa
Tại các nhiệt độ thấp, việc tính toán các phổ EXAFS có thể thực hiện trong gần đúng điều hòa vì các đóng góp phi điều hòa của các dao động nhiệt của nguyên tử là nhỏ nên có thể bỏ qua Khi nhiệt độ tăng cao, thế năng tương tác giữa các nguyên tử trở thành bất đối xứng bởi vì đã xuất hiện các số hạng phi điều hòa. Công thức của phổ EXAFS bao gồm các số hạng phi điều hòa thường được mô tả qua phương pháp gần đúng khai triển cumulant, theo đó hàm dao động EXAFS thường được viết như sau [18]: χ(k) F(k) e −2R/ λ(k) kR2 I m
, (2.1) trong đó, phần thực F(k) biểu diễn biên độ tán xạ nguyên tử ϕ(k) là tổng độ dịch pha của electron quang và λ(k) là quãng đường tự do trung bình của electron quang, σ (n) , (n = 1, 2, 3, ) là các cumulant, chúng xuất hiện do lấy trung bình nhiệt hàm e ikr , trong đó các số hạng bất đối xứng được khai triển theo chuỗi Taylor xung quanh giá trị
R = ⟨r⟩, với r là khoảng cách trung bình giữa nguyên tử hấp thụ và nguyên tử tán xạ tại nhiệt độ T Công thức (2.1) của hàm dao động EXAFS bao gồm các hiệu ứng phi điều hòa có chứa hệ số Debye - Waller do các hiệu ứng dao động nhiệt của các nguyên tử Trong phân tích của các tác giả [44, 72] thì hệ số Debye-Waller của phổ EXAFS sẽ là e ω(k) với: ω(k) = Σ
(T ) + (2.2) trong đó σ (1) là cumulant bậc một hay sự giãn nở nhiệt mạng, σ (2)
= σ 2 là cumulant bậc hai hay độ dịch chuyển tương đối trung bình bình phương (MSRD), σ (3) và σ (4) là các cumulant bậc ba và bậc bốn, các tham số còn lại đã được nêu trong các phần trước Do hiệu ứng phi điều hòa thường là nhỏ nên sự phân tích EXAFS chỉ cần đến các cumulant tới bậc ba hoặc bậc bốn Các cumulant bậc cao hơn ta có thể bỏ qua vì đóng góp của chúng trong dao động nhiệt là rất nhỏ.
Như vậy trong công thức (2.2) các số hạng thứ hai (DWF) và số hạng thứ năm đóng góp vào sự thay đổi biên độ, các số hạng thứ nhất, thứ ba và thứ tư đóng góp vào độ dịch pha của các phổ EXAFS do hiệu ứng phi điều hòa.
Từ các phương trình (2.1) và (2.2), ta có thể viết lại phổ EXAFS phi điều hòa bằng công thức sau: Σ S 2 N j 2 2
Kết hợp với hệ thức (2.2) thì phương trình (2.3) trở thành phổ EXAFS cận K phụ thuộc vào nhiệt độ bao gồm các hiệu ứng phi điều hòa: Σ S 2 N j
Trong các phương trình trên S 2 đặc trưng cho hiệu ứng nhiều hạt, N j là số nguyên tử lân cận thuộc lớp thứ j, các hệ số khác đã được định nghĩa χ(k, T ) χ(k, T )
A ở trên và dấu tổng là tính trên tất cả các lớp nguyên tử.
Chú ý rằng từ phương trình (2.4), thành phần σ 2 (T ) xác định đóng góp phi điều hòa vào biên độ và làm suy giảm biên độ dao động, còn ϕ A (k, T ) là đóng góp phi điều hòa vào độ dịch pha của phổ EXAFS Các đóng góp phi điều hòa này được biểu diễn qua các cumulant Tại nhiệt độ thấp các giá trị này tiến gần tới không và công thức của phổ EXAFS(2.4) sẽ rút về mô hình dao động điều hòa.
Mô hình Debye và Einstein trong EXAFS phi điều hòa
Để xác định biểu thức của các cumulant và tham số nhiệt động trong EXAFS phi điều hòa, có nhiều nghiên cứu đã phát triển, xây dựng và đưa ra các phương pháp và mô hình tính toán khác nhau, mỗi mô hình đều có những điểm mạnh nhưng cũng có một số hạn chế nhất định Trong phạm vi nghiên cứu, luận án sử dụng hai mô hình Debye và Einstein tương quan phi điều hòa.
2.2.1 Mô hình Debye tương quan phi điều hòa
Mô hình Debye tương quan phi điều hòa (ACDM - Anharmonic Correlation Debye Model) được xây dựng để đưa ra một phương pháp tính giải tích cho các tham số nhiệt động và các cumulant phổ EXAFS của các hệ vật liệu, trong đó đã tính tới đóng góp của các thành phần phi điều hòa [45, 49, 52] ACDM được dựa trên các ý tưởng chính là:
(1) Xét sự đóng góp tương quan của các nguyên tử lân cận bao gồm tương tác của nguyên tử hấp thụ với nguyên tử tán xạ cùng các nguyên tử lân cận trong một chùm nhỏ nguyên tử và có tính đến sự tán sắc của các phonon;
(2) Sử dụng hàm thế tương tác hiệu dụng có chứa đóng góp của các thành phần phi điều hòa, phù hợp cho các khai triển đối với độ dịch chuyển mạng nhỏ;
(3) Thành phần phi điều hòa được coi là những nhiễu loạn và là kết quả của tương tác phonon - phonon, trong đó độ dịch chuyển mạng được biểu diễn qua toán tử dịch chuyển phonon và kết quả thu được nhờ phép gần đúng nhiễu loạn hệ nhiều hạt. Σ Σ Σ ℜ Để xem xét ACDM, ta bắt đầu từ phép khai triển gần đúng cumulant của hàm EXAFS phi điều hòa dựa theo hệ thức (1.21) có dạng: ex p⟨2ikR
Hệ thức (2.5) là tổng độ dịch chuyển pha, R = ⟨r⟩ với r là độ dài liên kết tức thời giữa các nguyên tử hấp thụ và tán xạ, σ (n) (n = 1, 2, 3,
4, ) là các cumulant Xét hệ chỉ gồm một loại nguyên tử và khai triển thế năng tương tác hiệu dụng tới bậc 4 ta có:
(x) k 2 eff x 2 + k3 x 3 + k4 x 4 , (2.6) với k eff là hằng số lực hiệu dụng, k3 và k4 là các hệ số phi điều hòa hiệu dụng mô tả sự bất đối xứng của thế phi điều hòa, x là độ lệch của khoảng cách liên kết tức thời giữa hai nguyên tử liền kề khỏi giá trị trung bình của chúng Thế tương tác hiệu dụng được xác định dựa trên việc xem xét dao động tương đối giữa các nguyên tử hấp thụ (0) và tán xạ (1) có tính đến sự tương tác với tất cả các nguyên tử lân cận được xác định nhờ biểu thức:
1.ℜ ˆ ij ), (2.7) trong đó U(x) là thế đơn cặp của nguyên tử hấp thụ và tán xạ Số hạng thứ hai là tổng lấy trên toàn bộ vùng không gian liền kề còn lại trừ nguyên tử hấp thụ và tán xạ, ℜˆ01.ℜˆij là các vectơ đơn vị So sánh phương trình (2.6) với (2.7) ta xác định được các hằng số hiệu dụng k eff , k3, k4 của thế năng hiệu dụng phi điều hòa Biểu thức giải tích của các hằng số này phụ thuộc vào biểu thức thế năng được sử dụng và cấu trúc của mạng tinh thể được nghiên cứu.
Do có sự tán sắc của phonon, chúng ta biểu diễn thông số x qua toán tử dịch chuyển phonon có chứa thành phần tổng thống kê theo các tần số dao động của phonon Xét hệ dao động gồm N dao động tử có tần số thay đổi từ 0 đến tần số Debye cực đại ω D , trong hệ một chiều chỉ gồm một loại nguyên
(e iaq − 1)A , (2.13) tử, hàm phân bố được cho bởi [58]; ω(q) = 2
, (2.8) với q là số sóng, M là khối lượng của cặp nguyên tử hấp thụ và tán xạ, a là hằng số mạng Trong hàm thế liên kết hiệu dụng phi điều hòa, đóng góp của các nguyên tử lân cận được tính theo hình chiếu lên phương của đường nối giữa nguyên tử hấp thụ và tán xạ nên ta hoàn toàn có thể sử dụng hàm phân bố của hệ một chiều ACDM.
Tại biên của vùng Brillouin (BZ) thứ nhất q = ±π/a, tần số nhận giá trị cực đại, ta thu được giá trị tần số Debye tương quan ω D và nhiệt độ Debye θ D có dạng như sau: ω = 2 keff , θ
= ℏD , (2.9) k B với k B là hằng số Boltzmann Biểu diễn thông số x theo độ dịch chuyển của nguyên tử thứ n theo hệ thức: x n = u(n+1) − u n , (2.10) trong đó độ dịch chuyển u n liên hệ với toán tử dịch chuyển phonon A q
Từ (2.11), ta thu được biểu thức x n có dạng:
2 NM q eiqa u n x n q nê Σ n q 2NMω(q) x n = e iaqn f(q)A q , (2.14) q
2NMω(q) Để tính toán các hiệu ứng phi điều hòa, Hamiltonian của hệ được viết lại dưới dạng là tổng của các thành phần điều hòa H0 và phi điều hòa H a , trong đó thành phần phi điều hòa H a bao gồm thành phần bậc ba và bậc bốn:
Thay thế (2.14) vào (2.17) ta có:
. e iq 1anf(q1)f(A q ) Σ. e iq 2anf(q2)f(A q ) Σ. e iq 3anf(q3)f(A q ) ΣΣ
H c = k3 Σ Σ e i(q 1+q2+q3)anf(q1)f(q2)f(q3)A q A q A q Σ (2.18) Đặt ∆(q) = Σ e iqan , ∆(0) = Σ e 0.ian = N, ta có q 1 ,q 2 ,q 3
So sánh (2.19) với (2.16) ta rút ra :
Tương tự, để xác định biểu thức hàm U(q1234), ta thay biểu thức
= k4 Σ Σ Σe iq 1anf(q1)A q Σ Σe iq 2anf(q2)A q Σ Σe iq 3anf(q3)A q Σ Σe iq 4anf(q)A q ΣΣ , n hay q 1 ,q 2 ,q 3 ,q 4 1
So sánh biểu thức (2.17) với biểu thức (2.22), ta rút ra:
Hàm U q (123) và U q (1234) có tính chất đối xứng theo các chỉ số
Mô hình Debye tương quan phi điều hòa được sử dụng nhiều trong phương pháp EXAFS hiện đại và đã thu được nhiều kết quả phù hợp tốt với thực nghiệm [45, 50, 52]. i
2.2.2 Mô hình Einstein tương quan phi điều hòa
Mô hình Einstein tương quan phi điều hòa (ACEM - Anharmomic corre- lated Einstein model) đã đưa ra một phương pháp tính giải tích các tham số nhiệt động và các cumulant phổ EXAFS của các hệ vật liệu, trong đó cũng đã tính tới đóng góp của các thành phần phi điều hòa dựa vào sự đóng góp tương quan của một chùm nguyên tử lân cận gần nhất [18, 44] ACEM được xây dựng tương tự như ACDM, nhưng để đơn giản người ta đã bỏ qua sự tán sắc phonon ACEM cũng sử dụng phép gần đúng khai triển cumulant theo hệ thức (2.5) Sự phát triển quan trọng của phương pháp này là mô hình ACEM đã tính đến sự tương tác giữa nguyên tử hấp thụ và nguyên tử tán xạ với các nguyên tử lân cận trong một chùm nhỏ nguyên tử Để xác định các cumulant, chúng ta khai triển hàm thế tương tác hiệu dụng phi điều hòa đến bậc bốn [44]:
2k eff x 2 + k3 x 3 + k4 x 4 , (2.24) trong đó x = r − r0 là độ dịch chuyển của nguyên tử khỏi vị trí cân bằng với r là khoảng cách giữa hai nguyên tử ở nhiệt độ tuyệt đối T và r0 là khoảng cách cân bằng giữa hai nguyên tử mà tại đó hàm thế đạt giá trị cực tiểu, k eff là hệ số đàn hồi hiệu dụng và nó bao gồm tất cả tương tác của các nguyên tử lân cận, k3 là tham số bậc ba đặc trưng cho tính phi điều hòa được xác định bằng dao động của một liên kết đơn cặp giữa nguyên tử hấp thụ và nguyên tử tán xạ (có khối lượng là M1 và M2 ) Dao động của chúng bị ảnh hưởng bởi các nguyên tử lân cận nên thế tương tác hiệu dụng (2.24) trong mô hình Einstein tương quan phi điều hòa có dạng:
1 2 gọi là khối lượng rút gọn,
ℜˆ là vectơ đơn vị, U(x) đặc trưng cho thế đơn cặp giữa nguyên tử hấp thụ và nguyên tử tánxạ, Σ i j, i=0,1
Phép khai triển cumulant dựa vào ACEM trong phổ EXAFSphi điều hòa
Khi thế năng tương tác giữa các nguyên tử không đối xứng, nghĩa là phải tính đến các hiệu ứng phi điều hòa [2, 18] chúng ta phải sử dụng phép khai triển cumulant trong phổ EXAFS Trong ACEM, phép khai triển này chủ yếu dựa vào công thức của hàm: e⟨2ikr
( n)Σ , n = 1, 2, 3, , (2.32) trong đó σ (n) là các cumulant, r0 là khoảng cách ban đầu giữa các nguyên tử hấp thụ và nguyên tử tán xạ khi chúng ở vị trí cân bằng hay tại vị trí có thế năng tương tác cực tiểu Theo ACEM, r là khoảng cách giữa các nguyên tử tại nhiệt độ T , x = r − r0 là độ lệch liên kết tức thời giữa hai nguyên tử, a(T ) = ⟨r − r0⟩ = σ (1) là độ giãn nở nhiệt mạng, y = x − a và ⟨y⟩ = 0 đã
Σ Σ Σ Σ Σ được định nghĩa. Áp dụng việc khai triển cumulant đã dẫn trong phần trước, phân tích các cumulant dựa vào quan hệ của chúng với các mômen của hàm phân bố ta có: σ (1) = R − r, ⟨y⟩ = 0, σ (2) = σ 2 =
Do các hiệu ứng phi điều hòa thường là nhỏ nên khi sử dụng mô hình Einstein tương quan phi điều hòa trong phổ EXAFS chỉ cần khai triển đến các cumulant tới bậc 3 hoặc bậc 4 Các cumulant bậc cao hơn có thể bỏ qua vì chúng có đóng góp không đáng kể vào phổ EXAFS Trong các hệ thức cumulant (2.33), cumulant bậc hai σ (2) = σ 2 là độ dịch chuyển tương đối trung bình bình phương và đóng góp chủ yếu vào hệ số Debye – Waller nên còn được gọi là hệ số Debye – Waller (DWF).
Sử dụng lý thuyết nhiễu loạn [1, 48] để xây dựng công thức tính giải tích các cumulant và hệ số dãn nở nhiệt và áp dụng cho các tinh thể có cấu trúc lập phương Dao động của các nguyên tử đã được lượng tử hoá là các phonon và tính phi điều hòa là kết quả của tương tác phonon-phonon, đại lượng ⟨y⟩ được mô tả qua các toán tử sinh hạt và huỷ hạt aˆ + , aˆ theo hệ thức [1, 2]:
= e m và sử dụng trạng thái dao động điều hòa |n⟩ như các trạng thái riêng với các năng lượng là các giá trị riêng E n = nℏω E , trong đó, để thuận tiện ta bỏ qua năng lượng điểm không Các toán tử aˆ + , aˆ thoả mãn các tính chất sau: Σaˆ, aˆ + Σ
(2.36) Các cumulant được tính trên cơ sở xác định giá trị trung bình ⟨y m ⟩:
Tr(ρy Z m ), (2.37) trong đó Z là tổng thống kê, ρ = exp(−βH) là ma trận mật độ thống kê với β = 1/k B T Đối với trường hợp không nhiễu loạn thì Z0 Trρ0, với ρ0 = exp(−βH0) Do sử dụng trạng thái dao động tử điều hòa |n⟩ với các giá trị riêng E n = nβℏω E , ta sẽ có:
) = ∞ z = , (2.38) n=0 1 − z trong đó z = exp(−βℏω E ) = exp(−θ E /T ) gọi là biến số nhiệt độ và được xác định bởi nhiệt độ Einstein θ E = ℏω E /k B
Như vậy, nhiễu loạn δU E do hiệu ứng phi điều hòa sẽ dần tới giá trị δρ của ma trận mật độ, ρ = ρ0 + δρ Từ biểu thức của ma trận mật độ [1, 4]: ta nhận được: δρ = Hρ, δρ0 δβ δβ
−βH Áp dụng (2.37) và các công thức trên để tính giải tích các cumulant Ta có các cumulant bậc lẻ: n δρ −
Thay (2.40) vào (2.41) ta nhận được;
Trong phạm vi tính toán, luận án áp dụng tính với các tinh thể có cấu trúc lập phương tâm diện (fcc) Để tính các cumulant bậc lẻ, với m = 1, 3 (tương ứng với ⟨y⟩ và y 3 Σ ), ta sử dụng biểu thức (2.33) và các yếu tố ma trận:
Chú ý rằng từ bước biến đổi thứ hai của dãy hệ thức chúng ta đã thay n ′ = n + 1 Theo (2.37) và (2.41) thì:
Sử dụng các hệ thức của nhóm (2.43) ta có:
3α(σ0) 2 (1 + z) a = (2.44) Đồng thời ta lại có
, 2keff nên (2.44) được viết lại thành a = 3α
Vậy ta thu được cumulant bậc một σ (1) = a 3ℏωE
Tương tự như trên ta sử dụng các yếu tố ma trận trong nhóm (2.43), từ (2.42) ta có σ(3)
Vì các yếu tố ma trận chỉ có tác dụng với y 3 nên ta có Σ e−βnℏω E − e−β(n+3)ℏω E Σ n nℏω E
Thay giá trị của A và B vào (2.47) ta thu được
Ta có cumulant bậc ba viết lại theo (2.49) σ(3) = (ℏω E ) 2
+ Để tính các cumulant bậc chẵn, sử dụng hệ thức (2.37) ta có:
Sử dụng kết quả (2.38), ta có
0 2àω E 2 ke ff eff E hệ thức cumulant bậc 2 trở thành: σ(2) = y2 = ℏω E (1 + z)
Từ hệ thức của hằng số lực đàn hồi hiệu dụng:
Với gần đúng k eff ∼ c1Dα 2 thì ta có hệ thức của cumulant bậc hai được viết dưới dạng: σ(2) = y2Σ
Hệ số Debye-Waller phi điều hòa
H Để xét các đóng góp phi điều hòa vào MSRD hay DWF, luận án đã sử dụng phương pháp Willis và Pryor [60, 75] trong đó sự thay đổi của DWF theo sự thay đổi của nhiệt độ cho bởi hệ thức:
Với σ 2 (T ) là DWF điều hòa và σ 2 (T0) là DWF ở nhiệt độ T0, nhiệt độ này rất thấp để cho σ 2 (T0) là MSRD điều hòa, ta có thể viết:
V trong các hệ thức trên, γ tt là hệ số Gruneisen, ∆V/V là sự thay đổi thể tích tương đối do giãn nở nhiệt, hiện tượng này chỉ xảy ra khi có hiệu ứng dao động phi điều hòa Phát triển tiếp hệ thức (2.61) chúng ta sẽ thu được MSRD tổng cộng: σ 2 (T ) = σ 2 (T ) + β(T ) Σ σ 2 (T ) − σ 2 (T0)Σ
Khi nhiệt độ T0 rất thấp, gần tới không độ tuyệt đối thì MSRD có giá trị rất nhỏ và khi đó σ 2 (T0) → σ 2 Từ đây ta có thể coi như trong phương trình (2.61) MSRD tổng cộng σ 2 (T ) tại một nhiệt độ T là tổng của thành phần điều hòa σ 2 (T ) và phần đóng góp phi điều hòa σ 2 (T ) dưới dạng:
Các biểu thức này xác định phần đóng góp phi điều hòa vào biên độ của phổ EXAFS phi điều hòa.
Hệ số giãn nở nhiệt
Trong gần đúng điều hòa, độ dịch chuyển nguyên tử được coi bằng không tại vị trí cân bằng, năng lượng tự do sẽ có giá trị cực tiểu, vì thế hiện tượng giãn nở nhiệt không xảy ra ở gần đúng điều hòa Khi nhiệt độ tăng thì biên độ dao động của các nguyên tử cũng tăng và độ dịch chuyển nguyên tử trong ô mạng cơ sở cũng tăng, các thành phần phi điều hòa xuất hiện và đóng góp vào năng lượng tự do của các tinh thể nên năng lượng tự do này không bị cực tiểu tại vị trí cân bằng như trong gần đúng điều hòa.
Do ảnh hưởng của hiệu ứng phi điều hòa nên toàn bộ tinh thể bị giãn nở nhiệt và đạt đến thể tích trong đó năng lượng tự do có giá trị cực tiểu.
2.5.1 Hệ số giãn nở khối
Hệ số giãn nở nhiệt toàn phần (giãn nở thể tích hay giãn nở khối) được xác định theo hệ thức [1, 2, 4]; α = 1
, (2.66) trong đó V là thể tích của tinh thể phụ thuộc vào sự biến đổi của nhiệt độ tuyệt đối T dưới áp suất P không đổi Sử dụng phương trình trạng thái của hệ nhiệt động, ta có hay
Từ (2.66) và (2.67) chúng ta thu được α = − 1
. Đại lượng này gọi là modun nén khối đẳng nhiệt và xác định sự thay đổi của thể tích dưới tác dụng của áp suất Ta có: α = − 1
Bỏ qua sự liên kết giữa các dao động và giả thiết rằng năng lượng tự do
T với U là tổng thế năng không phụ thuộc vào nhiệt độ và được tạo ra bởi tương tác giữa các nguyên tử, F q là năng lượng tự do được sinh ra từ các dao động mạng với vectơ sóng q Sự phụ thuộc của áp suất vào thể tích được xác định bởi hệ thức:
Do năng lượng của một dao động tử điều hòa với tần số ω là: ε n = (n + 1/2)ℏω, (n là số nguyên), nên một dao động tử điều hòa đóng góp vào năng lượng tự do của hệ sẽ được xác định bằng hệ thức:
F q = −k B T lnZ, (2.72) trong đó Z là tổng thống kê, với n là các số nguyên nên Tổng thống kê Z có dạng:
Giả sử các mức năng lượng không bị suy biến Thay hệ thức (2.73) vào hệ thức (2.72) ta có hệ thức năng lượng tự do:
Xét trong gần đúng điều hòa, tần số dao động mạng không phụ thuộc thể tích, vì thế các hệ thức (2.74), (2.70), (2.71) cho thấy năng lượng của một dao động tử điều hòa F q , năng lượng tự do Helmholtz F không phụ thuộc thể tích, như vậy các dao động mạng không đóng góp vào áp suất và giãn nở nhiệt Trong hệ thức (2.70), tuy thành
B áp suất nhưng nó lại không phụ thuộc vào nhiệt độ nên nó cũng không có sự đóng góp vào giãn nở nhiệt Như vậy, trong gần đúng điều hòa không có giãn nở nhiệt.
Khi nhiệt độ tăng, hiệu ứng phi điều hòa xuất hiện và hiện tượng quan trọng của hiệu ứng phi điều hòa là sự phụ thuộc của tần số dao động mạng vào thể tích, bởi vì khi có hiệu ứng phi điều hòa hệ sẽ cân bằng ở vị trí mới với thể tích bị giãn nở Bỏ qua sự liên kết giữa các dao động và giả sử năng lượng tự do Helmholtz F vẫn bằng tổng của thế năng tương tác giữa các nguyên tử U và tổng năng lượng tự do của các dao động mạng F q trong hệ thức (2.74) Khi đó áp suất phụ thuộc vào thể tích theo hệ thức (2.71) sẽ là:
Sự phụ thuộc của tần số dao động mạng vào thể tích được thể hiện qua đạo hàm ∂ω q /∂V ở số hạng thứ hai trong (2.75) Giả thiết đơn giản nhất là sự phụ thuộc V của tất cả các tần số dao động mạng là như nhau và có thể biểu diễn qua hệ số Gruneisen dưới dạng; ω ∼ V −γ G ⇒ γ tt d(ln ω )
Hệ số Gruneisen γ tt là đại lượng đặc trưng cho hiệu ứng phi điều hòa và đã được xác định cho một số kim loại có cấu trúc tinh thể [17]
Từ hệ thức (2.76) ta có
Thay hệ thức (2.77) vào hệ thức (2.75) ta sẽ có
Từ các hệ thức (2.72) và (2.73), năng lượng trung bình của một dao động tử điều hòa với tần số ω q sẽ được tính bằng biểu thức: dV exp k T ε¯ = k B T 2 d(lnZ) dT = ℏω q
Sử dụng các hệ thức (2.77), (2.78) và (2.79), ta sẽ có hệ thức của áp suất phụ thuộc vào thể tích dạng:
Thay hệ thức (2.80) vào (2.69) ta sẽ thu được biểu thức của hệ số giãn nở nhiệt khối với C V là nhiệt dung đẳng tích của mạng tinh thể: α V = 1
2.5.2 Hệ số giãn nở tuyến tính
Hệ số giãn nở tuyến tính hay còn gọi là giãn nở dài là sự giãn nở theo một chiều của tinh thể khi nhiệt độ tăng cao, tương tự như hệ thức (2.66), hệ số giãn nở nhiệt tuyến tính được biểu diễn bằng hệ thức: α T = da , (2.82) rdT trong đó r là khoảng cách giữa hai nguyên tử ở nhiệt độ T , a là hệ số giãn nở nhiệt mạng với a(T ) = ⟨r − r0⟩ , r0 là khoảng cách giữa hai nguyên tử ở vị trí cân bằng và thế năng tương tác giữa hai nguyên tử là cực tiểu Trong ACEM, hệ số giãn nở nhiệt mạng a là cumulant bậc một σ (1) Vì vậy theo hệ thức (2.46), chúng ta có thể viết lại hệ thức của hệ số giãn nở dài (2.82) theo dạng:
⇒ α T 40Dαr T 2 (1 − z)2 với nhiệt độ Einstein θ E = ℏω BE , ta có:
Thay lnz = −θE = −ℏωE , ta thu được hệ số giãn nở dài:
Hệ số phi điều hòa và đóng góp của hiệu ứng phi điều hòa vào biên độ của phổ EXAFS
điều hòa vào biên độ của phổ EXAFS
Chúng ta có thể xác định hệ số phi điều hòa từ hệ thức β(T ) 2γ tt ∆V , với độ giãn nở khối tương đối ∆V/V do tính phi điều hòa.
Sử dụng hệ thức R(T ) = R + α(T ) với R là khoảng cách trung bình giữa nguyên tử hấp thụ và tán xạ, ta có thể xác định độ giãn nở khối của tinh thể khi nhiệt độ tăng lên T (K) bởi
Thay a(T ) từ hệ thức (2.46) vào (2.84) ta sẽ có:
B nên hệ thức độ giãn nở khối tương đối của tinh thể có thể viết lại dưới dạng đơn giản:
Theo hệ thức (2.58) ta có hệ số đàn hồi hiệu dụng và tần số Einstein: keff = 5(Dα 2 3 aDα 3 ), 2
Từ hệ thức (2.86) ta có:
(2.87) Đạo hàm theo nhiệt độ hệ thức của hệ số giãn nở nhiệt mạng (2.46), ta lại có: da 3k B θ E Σ(1 − z)dz + (1 + z)dz Σ
Với z = e −θ E /T → dz/dT = z(θ E /T 2 ), biểu thức (2.88) được viết lại thành: da = 3k B θ E Σ−(1 − z)z − (1 + z)z Σ θ E θ E dT
2 Tại nhiệt độ T , thể tích của tinh thể là:
Từ đó ta sẽ có:
∂T 16DαT 2(R + a) (1 − z)2 Thay các hệ thức (2.91) và (2.92) vào (2.76) ta sẽ được hệ số Gruneisen γ tt
Hệ số phi điều hòa (2.93) là hàm của nhiệt độ tuyệt đối T và tỉ lệ nghịch với bán kính R Như vậy nó phản ánh tính chất phi điều hòa như đã nhận được từ thực nghiệm [21, 38] Khi đó, các đóng góp phi điều hòa vào biên độ của phổ EXAFS sẽ tỉ lệ với hàm exp{−β(T )[σ 2 (T )
Pha của phổ EXAFS phi điều hòa
Trong hệ thức (2.2), phần đóng góp vào độ dịch pha của phổ EXAFS phi điều hòa bao gồm số hạng thứ nhất, số hạng thứ ba và số hạng thứ tư Đồng thời từ công thức (2.63) ta thấy đóng góp phi điều hòa vào DWF tại một nhiệt độ đã cho được xác định bởi: σ 2 (T ) = β(T ) Σ σ 2 (T ) − σ 2 (T0) Σ
Theo phân tích đã đưa ra ở trên và từ phương trình (2.2) ta xác định được biểu thức về độ dịch pha ϕ A (T ) của hàm dao động EXAFS phi điều hòa: ϕ (T ) = 2k Σ σ (1) (T ) − 2βσ 2 (T ) 1
Theo hệ thức (2.94) thì (2.95) được viết lại thành: ϕ (T, k) = 2k Σ σ (1) (T ) − 2σ 2 (T ) 1
Biểu thức trên cho thấy, độ dịch pha ϕ A (T ) sẽ giảm nhanh tại các nhiệt độ thấp vì các giá trị phi điều hòa σ (1) , σ (3) , β là nhỏ không đáng kể.
Cumulant bậc 4 trở lên có đóng góp hiệu ứng phi điều hòa làm thay đổi phổ EXAFS nhưng chúng rất nhỏ nên đã bỏ qua Hệ số phi điều hòa phụ thuộc vào nhiệt độ T và hệ số Gruneisen γ tt và thay đổi thể tích tương đối do giãn nở nhiệt ∆V , hai đại lượng này chỉ xuất hiện khi có hiệu ứng phi điều hòa.
Phổ EXAFS cận K phụ thuộc vào nhiệt độ bao gồm hiệu ứng phi điều hòa được mô tả bằng phương trình (2.4), trong đó σ 2 (T ) xác định rõ đóng góp phi điều hòa vào biên độ và làm dao động EXAFS tắt dần, thành phần này
4 nằm trong hàm e mũ chứa hệ số phi điều hòa β(T ) theo hệ thức (2.93) Tức là biên độ phổ EXAFS với đóng góp phi điều hòa tỷ lệ với hàm e −β(T ) σ 2 (T ) Phần đóng góp phi điều hòa vào pha của phổ EXAFS ϕ A (k, T ) trong phương trình (2.4) mô tả bằng hệ thức (2.96). Các đóng góp phi điều hòa này được biểu diễn qua các cumulant, tại các nhiệt độ thấp các giá trị này tiến dần tới 0 và công thức của phổ EXAFS phi điều hòa sẽ rút về mô hình dao động điều hòa theo dạng của hệ thức (1.8)
Hiệu ứng lượng tử ở các nhiệt độ giới hạn
2.8.1 Biểu diễn các tham số nhiệt động qua cumulant bậc 2 Để giảm bớt các phép tính toán số và các phép đo ta có thể biểu diễn biến số nhiệt độ z qua cumulant bậc 2 σ 2 hay DWF dưới dạng: σ 2 − σ 2 z = 0 (2.97) σ2 + σ 2
Thay hệ thức (2.97) vào các biểu thức (2.46), (2.51), (2.59) và (2.83) ta sẽ có các biểu thức đơn giản hơn của các cumulant và hệ số giãn nở nhiệt: σ(1) = σ(1) 1
Trong các hệ thức trên, các đại lượng σ (1) , σ 2 và σ(3) là các đóng góp điểm
0 0 0 không của các cumulant σ (1) , σ 2 và σ (3) Trong phổ EXAFS, người ta thường thiết lập các hệ thức tương quan giữa các cumulant với hệ số giãn nở nhiệt
2 − α T , khoảng cách giữa các nguyên tử r và nhiệt độ tuyệt đối T theo các hệ số cấu trúc và qua cumulant bậc 2 hay hệ số DWF:
Theo cách biểu diễn của các hệ thức (2.98) – (2.103), ngoài các tham số thế Morse đã cho, để tính các tham số nhiệt động ta chỉ cần tính DWF hay cumulant bậc hai σ 2 Như vậy, khi tính toán hoặc đo được σ 2 chúng ta có thể xác định được các tham số nhiệt động và cumulant khác, điều đó làm giảm các phép đo cũng như trong lập trình tính số cho các tham số nhiệt động.
2.8.2 Hiệu ứng lượng tử ở các nhiệt độ giới hạn
Các công thức ở trên thu được từ lý thuyết lượng tử nên có thể áp dụng cho mọi nhiệt độ, trong đó ở nhiệt độ cao, nó bao gồm các kết quả của gần đúng lý thuyết cổ điển và ở giới hạn nhiệt độ thấp nó luôn giữ các hiệu ứng lượng tử thể hiện qua đóng góp của năng lượng điểm không Trong giới hạn nhiệt độ cao (High Temperature – HT), ta có θ E /T