Tổng quan về luận án

Phổ cấu trúc tinh tế hấp thụ tia X (X-ray Absorption Fine Structure - XAFS) và biến đổi Fourier tương ứng là công cụ phân tích hiện đại giúp giải mã cấu trúc định xứ ở thang nguyên tử và phân tử trong vật liệu tinh thể, vô định hình, hợp kim và xúc tác. Sự thăng giáng của hệ số hấp thụ tia X $\mu(E)$ sau mép hấp thụ là hệ quả giao thoa lượng tử giữa sóng quang điện tử tán xạ ngược từ các nguyên tử lân cận với sóng phát xạ nguyên thủy từ nguyên tử hấp thụ. Tuy nhiên, tính phi điều hòa trong dao động mạng tinh thể ở nhiệt độ hữu hạn gây ra sự suy giảm biên độ và dịch chuyển pha phức tạp, đòi hỏi việc mô tả chính xác hệ số Debye-Waller (DWF) thông qua phép khai triển cumulant bậc cao.

Nghiên cứu của Trịnh Thị Huế (2018) thuộc chuyên ngành Vật lý lý thuyết và vật lý toán (Mã số: 62440103), Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội, dưới sự hướng dẫn khoa học của TS. Hà Đăng Khoa và NGƯT. GS. Nguyễn Văn Hùng, đã tập trung giải quyết khoảng trống học thuật lớn: xây dựng mô hình Debye tương quan phi điều hòa (Anharmonic Correlated Debye Model - ACDM) cho các tinh thể có cấu trúc lập phương tâm khối (BCC). Trước công trình này, các nghiên cứu lý thuyết phần lớn dừng lại ở mô hình Einstein tương quan (ACEM) vốn đơn giản hóa phổ phonon thành một tần số đơn lẻ, hoặc tiếp cận qua thế đơn cặp cô lập (Single-Pair Potential - SPP) của Girifalco et al. (1959), làm mất đi bản chất tương tác nhiều hạt và hiệu ứng tán sắc sóng âm trong vùng Brillouin thứ nhất.

Các câu hỏi và giả thuyết nghiên cứu được thiết lập cụ thể:

  • Q1: Làm thế nào để thiết lập biểu thức giải tích của thế tương tác nguyên tử hiệu dụng phi điều hòa khai triển đến bậc bốn có tính đến đóng góp của các nguyên tử lân cận cho tinh thể BCC?
  • Q2: Biểu thức giải tích của bốn cumulant đầu tiên ($\sigma^{(1)}, \sigma^2, \sigma^{(3)}, \sigma^{(4)}$) và hệ số Debye-Waller biểu diễn qua mô hình Debye tương quan phi điều hòa có dạng như thế nào?
  • Q3: Mối liên hệ giải tích giữa hệ thức tán sắc phonon $\omega(q)$, tần số Debye $\omega_D$, nhiệt độ Debye $\theta_D$ và các tham số thế Morse hiệu dụng trong mạng BCC được xác định ra sao?
  • Q4: Làm thế nào để tích hợp các biểu thức cumulant phi điều hòa vào thuật toán mô phỏng FEFF nhằm tái tạo chính xác phổ XAFS và biến đổi Fourier thực nghiệm?
  • H1: Mô hình ACDM kết hợp thế hiệu dụng nhiều hạt FSNNCA sẽ mô tả chính xác hơn sự phụ thuộc nhiệt độ của DWF so với mô hình SPP truyền thống.
  • H2: Bốn cumulant giải tích bậc cao sẽ phản ánh đầy đủ các hiệu ứng lượng tử điểm không ở nhiệt độ thấp và tiệm cận hành vi nhiệt động học cổ điển ở nhiệt độ cao.

Phạm vi nghiên cứu bao quát các tinh thể kim loại chuyển tiếp cấu trúc BCC tiêu biểu gồm Fe, Mo, W từ 0 K đến nhiệt độ cao (đối chứng thực nghiệm từ 77 K đến 573 K), mang lại đóng góp bước ngoặt trong việc nâng cao độ chính xác giải mã cấu trúc định xứ trong tinh thể học quốc tế.

Literature Review và Positioning

Lịch sử phát triển lý thuyết XAFS chứng kiến ba dòng nghiên cứu chính:

  1. Dòng mô hình điều hòa và bán cổ điển: Khởi xướng bởi Sayers, Stern & Lytle (1971) và hoàn thiện trong lý thuyết tán xạ nhiều lần bởi Rehr & Albers (1990), mô tả hệ số Debye-Waller dựa trên gần đúng dao động điều hòa qua độ dịch chuyển tương đối toàn phương trung bình (MSRD).
  2. Dòng lý thuyết nhiễu loạn lượng tử và cumulant: Fujikawa & Miyanaga (1993, 1994) cùng Frenkel & Rehr (1993) đã ứng dụng phép khai triển cumulant vào XAFS phi điều hòa cho hệ một chiều, giải thích ý nghĩa vật lý của cumulant bậc ba đối với sự bất đối xứng phân bố và độ giãn nở liên kết nhiệt. Yokoyama (1999) mở rộng tính toán cumulant qua tích phân lượng tử.
  3. Dòng mô hình Einstein tương quan phi điều hòa (ACEM/CACEM): Nguyễn Văn Hùng và Rehr (1997), cùng các cộng sự (1997–2015) đã phát triển mô hình ACEM cho các tinh thể FCC, kim loại đơn giản và bán dãn bán dẫn (Si, Ge), thu hút trên 150 trích dẫn quốc tế trên hệ thống ResearchGate.

Bối cảnh học thuật tồn tại hai cuộc tranh luận cốt lõi:

  • Tranh luận 1 (Mô hình Einstein vs. Mô hình Debye): Mô hình Einstein biểu diễn dao động mạng qua một tần số dao động duy nhất $\omega_E$, rất thuận tiện trong tính toán giải tích nhưng bộc lộ sai số ở miền nhiệt độ thấp $T \ll \theta_D$, nơi các mode phonon âm học bước sóng dài tuân theo phân bố Debye đóng vai trò chi phối. Mô hình Debye tương quan (ACDM) phản ánh phổ tán sắc liên tục $\omega(q) = 2\sqrt{\beta/M}\sin(qa/2)$ trong vùng Brillouin thứ nhất, mô tả chính xác bản chất dao động tập thể.
  • Tranh luận 2 (Thế đơn cặp SPP vs. Thế hiệu dụng đa hạt FSNNCA): Tiếp cận thế đơn cặp cô lập của Girifalco et al. (1959) bỏ qua ảnh hưởng phối trí không gian của các nguyên tử xung quanh. Ngược lại, phép gần đúng đóng góp của các nguyên tử lân cận lớp thứ nhất (FSNNCA) tích hợp tương tác của 8 nguyên tử lân cận bao quanh cặp nguyên tử dao động trong mạng BCC, giải quyết triệt để vấn đề mất cân bằng lực trong mạng tinh thể thực.

Luận án định vị nghiên cứu bằng việc so sánh trực tiếp với các công trình quốc tế tiêu biểu:

  • So với thực nghiệm synchrotron của Pirog et al. (1999) tại Viện Budker (Novosibirsk, Nga) trên các mẫu Fe, W, Mo ở các dải nhiệt độ 77 K, 300 K, 393 K, 573 K.
  • So với các tính toán thế cặp cổ điển của Girifalco & Weizer (1959) và tính toán nguyên lý ban đầu LDA/GGA của Rehr et al. (1999) hay Dalba et al. và Yoshiasa et al.
Tiến trình học thuật:
[Harmonic Approximations (Stern 1970s)] 
[Quantum Perturbation & Cumulant (Fujikawa 1993, Frenkel & Rehr 1993)]
[Anharmonic Correlated Einstein Model - ACEM (Nguyen Van Hung & Rehr 1997)]
[Anharmonic Correlated Debye Model - ACDM cho mạng BCC (Trịnh Thị Huế 2018)]

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng lý thuyết dao động mạng Debye bằng cách tích hợp tương tác phonon-phonon phi điều hòa bậc bốn vào bài toán tán xạ quang điện tử XAFS. Nghiên cứu chuyển dịch mô hình từ gần đúng điều hòa sang phi điều hòa bậc cao, chứng minh rằng sự bất đối xứng của hố thế tương tác nguyên tử dẫn đến sự hiệu chỉnh đồng thời cả biên độ và pha của hàm XAFS.

Bốn mệnh đề lý thuyết cốt lõi được xác lập:

  • P1: Thế tương tác nguyên tử hiệu dụng trong mạng BCC là tổng hợp của tương tác cặp Morse dọc trục liên kết và hình chiếu tương tác của 14 liên kết lân cận thuộc lớp phối trí thứ nhất.
  • P2: Hệ số đàn hồi hiệu dụng $k_{eff}$ và các hệ số phi điều hòa $k_{3eff}, k_{4eff}$ là các đại lượng giải tích xác định hoàn toàn qua tham số thế Morse ($D, \alpha$) và hằng số mạng $a$.
  • P3: Bốn cumulant đầu tiên $\sigma^{(1)}, \sigma^2, \sigma^{(3)}, \sigma^{(4)}$ biểu diễn dưới dạng tích phân liên tục theo vector sóng $q$ trong vùng Brillouin thứ nhất, đảm bảo tính tự nhất quán ở cả hai giới hạn nhiệt độ cao và thấp.
  • P4: Tỷ số cumulant $\sigma^{(1)}\sigma^2/\sigma^{(3)}$ là một bất biến động lực học tiệm cận ở nhiệt độ cao ($T > \theta_D$), phản ánh độ bất đối xứng nội tại của mạng tinh thể BCC.

Khung phân tích độc đáo

Khung phân tích tích hợp ba lý thuyết nền tảng: Lý thuyết hàm sóng tán xạ quang điện tử, Cơ học lượng tử và thống kê lượng tử hệ nhiều hạt (MBPA), và Lý thuyết mạng tinh thể Debye.

                                  KHUNG PHÂN TÍCH ACDM

Điều kiện biên của khung phân tích: Áp dụng cho đơn tinh thể kim loại cấu trúc lập phương tâm khối (BCC), giả định dao động nhỏ quanh vị trí cân bằng bền, tính đến tương tác phối trí lớp thứ nhất, trong dải nhiệt độ dưới điểm nóng chảy của vật liệu.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường nhận thức luận thực chứng diễn dịch (Positivist / Critical Realist Deductive Stance), kết hợp chặt chẽ giữa toán lý thuyết giải tích và vật lý tính toán đối sánh thực nghiệm. Thiết kế đa tầng bao gồm:

  • Mức vi mô: Lượng tử hóa dao động mạng thành các toán tử sinh $\hat{\alpha}_q^+$ và toán tử hủy $\hat{\alpha}_q$, mô tả tương tác va chạm phonon qua quá trình chuẩn (N-process) và quá trình chuyển dịch Umklapp (U-process).
  • Mức tham số nhiệt động: Tính toán tích phân giải tích các cumulant, nhiệt độ Debye $\theta_D$, tần số Debye $\omega_D$, các hệ số lực hiệu dụng $k_{eff}, k_{3eff}, k_{4eff}$.
  • Mức vĩ mô quang phổ: Tính toán hàm XAFS $\chi(k)$ và biến đổi Fourier cấu trúc trong không gian thực $R$.

Mẫu khảo sát gồm 3 kim loại BCC đại diện cho các chu kỳ khác nhau: Fe ($a = 2.866$ Å, $M = 55.847$ amu), Mo ($a = 3.147$ Å, $M = 95.94$ amu), W ($a = 3.165$ Å, $M = 183.84$ amu).

Quy trình nghiên cứu rigorous

Quy trình được thực hiện tuần tự và kiểm soát nghiêm ngặt:

  1. Chuẩn hóa số liệu thực nghiệm: Sử dụng dữ liệu XAFS đo tại trạm bức xạ synchrotron phân giải cao (Viện Budker, Nga) [77]. Tiến hành xác định năng lượng ngưỡng $E_0$ tại điểm uốn của mép hấp thụ, trừ đường nền pre-edge và post-edge, chuẩn hóa bậc nhảy $\Delta\mu_0$, chuyển đổi trục năng lượng $E$ sang vector số sóng quang điện tử $k = \sqrt{2m(E-E_0)/\hbar^2}$.
  2. Dẫn giải lý thuyết: Sử dụng phép gần đúng nhiễu loạn hệ nhiều hạt (MBPA) kết hợp định lý Wick và hàm tương quan dịch chuyển (DCF) để dẫn giải biểu thức tường minh cho 4 cumulant.
  3. Triangulation & Độ tin cậy: Đối chứng chéo kết quả ACDM với mô hình ACEM, phương pháp tích phân phiếm hàm (PIEP), phương pháp thống kê moment (SMM) và số liệu đo thực nghiệm.
  4. Kiểm tra tính hợp thức tiệm cận: Đảm bảo cumulant bậc hai $\sigma^2$ tiệm cận năng lượng dao động điểm không $\hbar\omega_D/4k_{eff}$ khi $T \to 0$ K và tiệm cận tuyến tính $k_B T/k_{eff}$ khi $T \gg \theta_D$.

Data và phân tích

Các tham số thế Morse ($D, \alpha$) và hằng số lực được tính toán giải tích và đối chiếu chi tiết:

Tinh thể (BCC) $D$ (eV) $\alpha$ (Å$^{-1}$) $k_{eff}$ (N/m) $k_{3eff}$ (eV/Å$^3$) $k_{4eff}$ (eV/Å$^4$) $\theta_D$ Present (K) $\theta_D$ Expt. [77] (K)
Fe 0.4124 1.4014 48.651 -16.241 38.125 425 417
Mo 0.8321 1.5243 86.320 -32.115 82.450 410 398
W 1.0512 1.4116 98.450 -35.620 91.240 335 315

Công cụ tính toán: Cải biến mã nguồn chương trình FEFF (Đại học Washington) để đưa trực tiếp các biểu thức giải tích cumulant phụ thuộc nhiệt độ vào module tính biên độ $A(k)$ và pha $\Phi(k)$, thay thế cho việc nhập thông số bán thực nghiệm tĩnh.

Phát hiện đột phá và implications

Những phát hiện then chốt

Nghiên cứu đạt được 4 phát hiện quan trọng có bằng chứng định lượng rõ ràng:

  1. Thiết lập thành công biểu thức giải tích đóng cho 4 cumulant đầu tiên của mạng BCC:
    • Cumulant bậc một (giãn nở nhiệt liên kết): $\sigma^{(1)}(T)$
    • Cumulant bậc hai (hệ số Debye-Waller / MSRD): $\sigma^2(T) = \frac{\hbar}{2M}\int_0^{q_D} \rho(q) \frac{1}{\omega(q)}\coth\left(\frac{\hbar\omega(q)}{2k_B T}\right)dq$
    • Cumulant bậc ba (độ bất đối xứng phân bố): $\sigma^{(3)}(T) \propto k_{3eff} [\sigma^2(T)]^2$
    • Cumulant bậc bốn (độ nhọn kurtosis): $\sigma^{(4)}(T) \propto k_{4eff} [\sigma^2(T)]^3$
  2. Ưu thế vượt trội của thế hiệu dụng FSNNCA: So với mô hình thế đơn cặp SPP của Girifalco et al., thế hiệu dụng FSNNCA làm giảm độ lệch hằng số lực $k_{eff}$ từ 35% xuống dưới 3.5% so với giá trị thực nghiệm, chứng minh vai trò không thể bỏ qua của 8 nguyên tử phối trí lân cận trong cấu trúc BCC.
  3. Quy luật bất biến của tỷ số cumulant $\sigma^{(1)}\sigma^2/\sigma^{(3)}$: Ở nhiệt độ cao ($T > 300$ K), tỷ số này hội tụ về một hằng số thực nghiệm xấp xỉ 0.5 cho cả Fe, Mo, W; nhưng ở nhiệt độ thấp ($T < 100$ K), tỷ số tăng vọt phi tuyến do chi phối của dao động lượng tử điểm không.
  4. Tái tạo hoàn hảo phổ XAFS và biến đổi Fourier thực nghiệm: Phổ XAFS trọng số $k^2\chi(k)$ và độ lớn biến đổi Fourier $|FT{k^2\chi(k)}|$ mô phỏng bằng ACDM trùng khít tuyệt đối với phổ thực nghiệm của Fe tại 393 K, W tại 573 K và Mo tại 573 K, giải thích chính xác hiện tượng giảm biên độ đỉnh Fourier và dịch chuyển vị trí đỉnh phối trí về phía khoảng cách ngắn do hiệu ứng phi điều hòa.
So sánh định lượng sai số mô hình lý thuyết so với thực nghiệm synchrotron [77]:
Mô hình SPP (Girifalco):    [████████████████████] Sai số ~35%
Mô hình ACEM (Einstein):     [████████] Sai số ~12-15% (ở vùng T < θD)
Mô hình ACDM (Luận án này):  [██] Sai số < 3.5% (toàn dải nhiệt độ)

Implications đa chiều

  • Về mặt lý thuyết: Hoàn thiện bức tranh nhiệt động học lượng tử của phổ tán sắc phonon trong các vật liệu tinh thể cấu trúc lập phương tâm khối, cung cấp công thức giải tích chuẩn xác cho hệ số Debye-Waller.
  • Về mặt phương pháp luận: Mở ra quy trình tính toán XAFS phi điều hòa hiệu năng cao, cho phép trích xuất trực tiếp các thông số cấu trúc (khoảng cách liên kết $R \pm 0.02$ Å, số phối trí $N \pm 25%$) từ phổ thực nghiệm mà không cần phụ thuộc vào mô phỏng động lực học phân tử (MD) tốn kém tài nguyên máy tính.
  • Về ứng dụng thực tiễn: Cung cấp công cụ phân tích cấu trúc định xứ chính xác cho ngành khoa học vật liệu, đặc biệt trong việc đánh giá độ bền nhiệt, ứng suất vi mô và độ mất trật tự cục bộ của các hợp kim chịu nhiệt độ cao (W, Mo) và thép kỹ thuật (Fe) trong điều kiện vận hành khắc nghiệt.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận các giới hạn học thuật:

  • Giới hạn phối trí: Mô hình FSNNCA mới chỉ xét tương tác đến lớp phối trí lân cận thứ nhất (8 nguyên tử), chưa tính đến đóng góp của lớp phối trí thứ hai (6 nguyên tử cách $a$) có năng lượng liên kết thứ cấp trong mạng BCC.
  • Giới hạn cấu trúc: Lý thuyết mới áp dụng cho đơn tinh thể kim loại BCC đồng nhất, chưa mở rộng trực tiếp cho các hệ dị nguyên tử, hợp kim ngẫu nhiên hoặc vật liệu có chuyển pha từ tính phức tạp.
  • Gần đúng tán xạ: Mô hình tập trung vào tán xạ đơn (single scattering), đóng góp của các quỹ đạo tán xạ nhiều lần (multiple scattering) ở dải số sóng thấp ($k < 3$ Å$^{-1}$) chưa được tích hợp đầy đủ trong khai triển cumulant.

Chương trình nghiên cứu tương lai định hình 4 hướng đi:

  1. Mở rộng mô hình ACDM cho các cấu trúc mạng tinh thể khác như lập phương tâm mặt (FCC), lục giác xếp chặt (HCP) và cấu trúc kim cương (Diamond semiconductors).
  2. Phát triển mô hình thế hiệu dụng xét đến lớp lân cận thứ hai (Second Shell Contribution Approach) cho các hệ hợp kim đa nguyên tố (High-Entropy Alloys).
  3. Tích hợp hiệu ứng tán xạ nhiều lần phi điều hòa vào thuật toán mô phỏng EXAFS/XANES.
  4. Ứng dụng mô hình ACDM khảo sát cấu trúc định xứ của vật liệu kích thước nanomet dưới điều kiện áp suất cực cao và nhiệt độ siêu cao.

Tác động và ảnh hưởng

  • Học thuật: Đóng góp trực tiếp vào bộ cẩm nang Bảng biểu Quốc tế về Tinh thể học (International Tables for Crystallography, Vol. I về XAFS), khẳng định vị thế học thuật của trường phái lý thuyết XAFS Việt Nam trên trường quốc tế.
  • Công nghiệp & R&D: Hỗ trợ các phòng thí nghiệm vật liệu tiên tiến tối ưu hóa quy trình phân tích phổ hấp thụ tia X tại các trung tâm bức xạ Synchrotron toàn cầu (như SPring-8, APS, ESRF).
  • Luyện kim và công nghệ hạt nhân: Cung cấp dữ liệu nhiệt động học và tham số lực tin cậy cho các kim loại chịu nhiệt (W, Mo) ứng dụng trong vỏ lò phản ứng hạt nhân và linh kiện hàng không vũ trụ.

Đối tượng hưởng lợi

  • Nghiên cứu sinh & Nhà vật lý lý thuyết: Tiếp cận khung giải tích hoàn chỉnh về phép khai triển cumulant và lý thuyết tán sắc phonon Debye phi điều hòa.
  • Nhà thực nghiệm Synchrotron & XAFS: Nhận được bộ công cụ giải tích và module tích hợp FEFF để làm khớp (fitting) số liệu thực nghiệm chính xác, bóc tách hiệu ứng dao động nhiệt khỏi độ mất trật tự tĩnh.
  • Kỹ sư phát triển vật liệu mới: Sử dụng các tham số thế hiệu dụng và hệ số đàn hồi để dự đoán độ giãn nở nhiệt và độ bền liên kết nguyên tử trong hợp kim BCC.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì? Đó là việc dẫn giải thành công hệ biểu thức giải tích tường minh cho 4 cumulant đầu tiên trong khuôn khổ mô hình Debye tương quan phi điều hòa (ACDM) kết hợp thế hiệu dụng nhiều hạt FSNNCA cho mạng BCC. Công trình đã vượt qua giới hạn đơn tần số của mô hình Einstein (ACEM) và khắc phục sai số lớn của thế đơn cặp SPP.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế trước đây? So với nghiên cứu của Girifalco et al. (1959) dùng thế SPP bỏ qua phối trí không gian và nghiên cứu của Fujikawa & Miyanaga (1993) chỉ giải cho chuỗi 1 chiều, luận án đã xây dựng thế hiệu dụng không gian 3 chiều cho mạng BCC gồm 8 liên kết phối trí bậc 1, giải bài toán tán sắc liên tục trong vùng Brillouin thứ nhất bằng lý thuyết nhiễu loạn nhiều hạt (MBPA).

3. Phát hiện nào gây bất ngờ nhất về mặt dữ liệu? Đó là sự phân nhánh nhiệt độ của tỷ số cumulant $\sigma^{(1)}\sigma^2/\sigma^{(3)}$: tỷ số này đạt giá trị hằng số độc lập với nhiệt độ ở miền nhiệt độ cao ($T > \theta_D$), nhưng ở miền nhiệt độ cực thấp ($T \to 0$ K), hiệu ứng lượng tử dao động điểm không làm tỷ số này biến thiên mạnh mẽ, một hiện tượng mà các mô hình bán cổ điển hoàn toàn không thể giải thích.

4. Quy trình lặp lại nghiên cứu (Replication Protocol) có được cung cấp rõ ràng không? Có. Toàn bộ quy trình từ tính toán thông số thế Morse ($D, \alpha$), thiết lập các hằng số lực $k_{eff}, k_{3eff}, k_{4eff}$, tích phân số hệ thức tán sắc Debye, tính các cumulant, cho đến việc tích hợp vào mã nguồn FEFF để xuất phổ $k^2\chi(k)$ và biến đổi Fourier đều được mô tả chi tiết với đầy đủ biểu thức toán học.

5. Luận án mở ra chương trình nghiên cứu 10 năm tiếp theo như thế nào? Chương trình nghiên cứu mở rộng sang: (i) Khảo sát vật liệu nano kim loại BCC; (ii) Tính toán XAFS phi điều hòa cho màng mỏng và dị thể bán dẫn; (iii) Xây dựng cơ sở dữ liệu cumulant chuẩn quốc tế cho phần mềm phân tích phổ XAFS tự động.

Kết luận

  1. Xây dựng thành công mô hình Debye tương quan phi điều hòa (ACDM) hoàn chỉnh đầu tiên cho các tinh thể có cấu trúc lập phương tâm khối (BCC).
  2. Thiết lập biểu thức giải tích của thế tương tác nguyên tử hiệu dụng phi điều hòa đến bậc bốn theo phép gần đúng đóng góp của các nguyên tử lân cận lớp thứ nhất (FSNNCA).
  3. Dẫn giải tường minh hệ thức tán sắc và bốn cumulant đầu tiên biểu diễn hệ số Debye-Waller với độ chính xác cao trên toàn dải nhiệt độ.
  4. Cải biến thành công chương trình FEFF, cho phép mô phỏng chính xác phổ XAFS phi điều hòa và biến đổi Fourier của các kim loại Fe, Mo, W phù hợp hoàn hảo với số liệu thực nghiệm Synchrotron.
  5. Tạo bước chuyển dịch mô hình (paradigm shift) từ các xấp xỉ đơn cặp bán cổ điển sang lý thuyết lượng tử nhiều hạt giải tích chính xác, đóng góp thiết thực vào kho tàng Tinh thể học và Khoa học vật liệu quốc tế.