Bối cảnh và vấn đề nghiên cứu

Trong vật lý hạt nhân và vật lý chất rắn, positron (phản hạt của electron, được Paul Dirac tiên đoán lý thuyết năm 1928 và Carl D. Anderson phát hiện thực nghiệm năm 1932) đóng vai trò trung tâm trong nhiều kỹ thuật thực nghiệm hiện đại. Quá trình hủy cặp electron - positron phát ra các photon gamma năng lượng cao mang lại nhiều ứng dụng thực tiễn quan trọng:

  • Trong khoa học vật liệu, phương pháp đo phổ thời gian sống của positron (Positron Annihilation Lifetime Spectroscopy) và kỹ thuật chụp cắt lớp vi tính công nghiệp (Computed Tomography - CT) được ứng dụng rộng rãi nhằm phát hiện các sai hỏng mạng, lỗ hổng vi mô và khuyết tật cấu trúc.
  • Trong y học hạt nhân, kỹ thuật chụp cắt lớp phát xạ positron (Positron Emission Tomography - PET) là công cụ phổ biến phục vụ chẩn đoán hình ảnh và theo dõi tiến triển của bệnh ung thư.

Các thông số thực nghiệm được ghi nhận gián tiếp thông qua thời gian sống của positron hoặc phổ xung lượng của các cặp photon hủy trong môi trường. Để giải thích một cách định lượng và chính xác các kết quả đo đạc này, việc mô tả tương tác và xác định năng lượng tương quan electron - positron là điều kiện tiên quyết. Kể từ thập niên 1950, nhiều nghiên cứu lý thuyết lượng tử đã được triển khai nhằm giải bài toán hủy positron trong môi trường chất rắn, tuy nhiên các phương pháp giải tích truyền thống gặp nhiều trở ngại do tính chất tương tác đa hạt phức tạp. Phương pháp biến phân Monte Carlo lượng tử (Variational Monte Carlo - VMC) nổi lên như một công cụ số hiệu quả, cho phép tính toán trực tiếp các hệ lượng tử vi mô phức tạp mà không bị giới hạn bởi các phép gần đúng giải tích cổ điển.

Mục tiêu và nhiệm vụ nghiên cứu của đề tài gồm 4 nội dung chính:

  1. Tổng quan cơ sở lý thuyết về tương tác giữa electron và positron trong vật chất, các trạng thái liên kết positronium và cơ chế nhiệt hóa, hủy cặp.
  2. Xây dựng mô hình vi mô và hàm sóng thử cho hệ electron trong phân tử đồng oxit (CuO) khi chưa có và khi có mặt positron trong khuôn khổ xấp xỉ Born - Oppenheimer.
  3. Ứng dụng phương pháp biến phân Monte Carlo lượng tử, lập trình thuật toán tính toán số trên máy tính nhằm cực tiểu hóa năng lượng và phương sai để tìm bộ tham số biến phân tối ưu.
  4. Xác định năng lượng trạng thái cơ bản của hệ electron trong phân tử CuO, năng lượng toàn phần khi có thêm positron, từ đó rút ra năng lượng tương quan electron - positron ($E_{corr}^{e-p}$).

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: Năng lượng tương quan positron - electron và cấu hình năng lượng trạng thái cơ bản trong phân tử đồng oxit (CuO).
  • Phạm vi không gian và mô hình: Hệ vi mô phân tử CuO gồm hai lõi nguyên tử (đồng và oxy) đặt trong hệ tọa độ Descartes ba chiều, xét chuyển động của 6 electron hóa trị ngoài lõi (4 electron thuộc phân lớp $2p$ của oxy và 2 electron thuộc phân lớp $4s$ của đồng) cùng một positron chuyển động trong trường thế hiệu dụng.
  • Phạm vi thời gian và công cụ: Đề tài thực hiện năm 2010 tại Khoa Vật lý, Trường Đại học Sư phạm TP. Hồ Chí Minh; công cụ tính toán số được xây dựng bằng ngôn ngữ lập trình Visual C++.

Cơ sở lý thuyết và phương pháp nghiên cứu

Về mặt lý thuyết lượng tử, đề tài xây dựng trên các trụ cột giải tích và mô hình xấp xỉ sau:

  • Nguyên lý biến phân: Năng lượng kỳ vọng $\langle H \rangle$ tính với hàm sóng thử $\Psi$ luôn đóng vai trò là giới hạn trên của năng lượng trạng thái cơ bản $E_0$ của toán tử Hamilton $H$ ($\langle H \rangle \ge E_0$).
  • Xấp xỉ Born - Oppenheimer: Tách biệt chuyển động của các hạt nhân nguyên tử nặng khỏi chuyển động của các hạt nhẹ (electron, positron). Do khối lượng hạt nhân rất lớn so với electron/positron và vận tốc dao động mạng nhỏ, vị trí của hai lõi nguyên tử đồng và oxy được coi là cố định tại các tọa độ xác định.
  • Hàm sóng quỹ đạo đơn hạt dạng Slater: Sử dụng hàm sóng dạng Slater $\psi(r) = N r^{n-1} \exp(-Zr)$ để mô tả các orbital của electron và positron quanh từng tâm lõi nguyên tử.
  • Hệ số tương quan Jastrow dạng Padé: Mô tả tương tác tương quan giữa các cặp hạt tích điện nhằm thỏa mãn điều kiện chặn lượng tử khi khoảng cách giữa hai hạt tiến về 0, có dạng giải tích $f(r) = \exp[a r / (1 + \beta r)]$.

Về phương pháp tính toán số, tác giả sử dụng phương pháp biến phân Monte Carlo lượng tử (VMC) kết hợp thuật toán chọn mẫu Metropolis:

  • Quy trình mô phỏng bao gồm hai giai đoạn: Giai đoạn thiết lập cân bằng nhiệt hóa (sử dụng khoảng 25% tổng số bước dịch chuyển Monte Carlo ban đầu để các hạt đạt phân bố không gian hội tụ theo mật độ $|\Psi|^2$) và giai đoạn lấy mẫu tính toán các đại lượng vật lý (năng lượng cục bộ $E_L(r) = H\Psi/\Psi$ và phương sai năng lượng $\sigma^2$).
  • Kỹ thuật tính đạo hàm vi phân: Toán tử động năng của từng hạt được phân tích chính xác qua các toán tử Gradient ($\nabla$) và Laplace ($\nabla^2$) tác dụng lên logarit tự nhiên của hàm sóng tổng, chuyển đổi việc tính tích phân nhiều chiều phức tạp thành bài toán lấy trung bình thống kê trên tập hợp cấu hình ngẫu nhiên.

Hệ đơn vị và hằng số vật lý: Toàn bộ biểu thức giải tích và thuật toán được chuẩn hóa trong hệ đơn vị nguyên tử (atomic units - a.u.), quy ước $\hbar = m_e = e = 1$. Các hằng số chuyển đổi thực tế được sử dụng gồm:

  • Bán kính Bohr: $a_0 = 5{,}291772108 \times 10^{-11}\text{ m}$.
  • Khối lượng electron: $m_e = 9{,}1093826 \times 10^{-31}\text{ kg}$.
  • Điện tích nguyên tố: $e = 1{,}60217653 \times 10^{-19}\text{ C}$.
  • Hằng số Planck rút gọn: $\hbar = 1{,}05457168 \times 10^{-34}\text{ J}\cdot\text{s}$.
  • Đơn vị năng lượng Hartree: $1\text{ Hartree} = 27{,}2113845\text{ eV} = 4{,}35974417 \times 10^{-18}\text{ J}$.
  • Hằng số Rydberg: $R_\infty = 1{,}0973731568525 \times 10^7\text{ m}^{-1}$ ($13{,}6056923\text{ eV}$).

Nội dung chính theo từng chương

Chương 1: Tổng quan về hệ electron - positron

Chương 1 hệ thống hóa các đặc trưng vật lý cơ bản của positron và các tương tác lượng tử khi hạt này xâm nhập vào môi trường vật chất:

  • Positron là phản hạt của electron, có cùng khối lượng, spin ($s = 1/2$), momen từ nhưng mang điện tích dương ($+e$). Khi positron gặp electron, kênh phân rã chiếm ưu thế là quá trình hủy cặp sinh ra hai tia gamma mang năng lượng đối xứng $511\text{ keV}$.
  • Cơ chế liên kết: Lực Coulomb giữa hạt nhân và positron thuần túy là lực đẩy. Tuy nhiên, đám mây điện tử xung quanh positron có thể tự tái phân bố và phân cực, tạo nên một thế hút phân cực hiệu dụng. Nếu thế hút phân cực thắng được thế đẩy tĩnh điện của hạt nhân, trạng thái liên kết positron - nguyên tử ($e^+A$) sẽ hình thành bền vững (đã được ghi nhận ở các hệ như $e^+\text{Be}, e^+\text{Na}, e^+\text{Cu}, e^+\text{O}$).
  • Hệ liên kết Positronium (Ps): Hình thành khi positron bắt giữ một electron hóa trị của nguyên tử, có năng lượng liên kết $0{,}250\text{ Hartree}$. Ps tồn tại ở hai trạng thái lượng tử: trạng thái para-positronium ($p$-Ps, spin singlet $s = 0, l = 0$, xác suất tạo thành $1/4$, phân rã chủ yếu thành 2 photon gamma) và trạng thái ortho-positronium ($o$-Ps, spin triplet $s = 1, l = 0$, xác suất tạo thành $3/4$, phân rã chủ yếu thành 3 photon gamma). Điều kiện liên kết của Ps với nguyên tử phụ thuộc vào thế ion $I$ của nguyên tử ($I < 0{,}250\text{ Hartree}$).
  • Quá trình nhiệt hóa: Positron mang năng lượng cao ban đầu (khoảng $2\text{ MeV}$) khi đi vào chất rắn va chạm liên tiếp với các electron, ion nút mạng và phonon dao động mạng, nhanh chóng tiêu tán động năng để trở thành positron nhiệt. Dẫn chứng nghiên cứu của Hautojärvi chỉ ra thời gian nhiệt hóa diễn ra dưới $20\text{ ps}$, trong khi thời gian sống trung bình của positron trong kim loại là khoảng $200\text{ ps}$. Do đó, positron luôn bị nhiệt hóa hoàn toàn trước khi xảy ra quá trình hủy cặp.

Chương 2: Phương pháp biến phân Monte Carlo lượng tử

Chương 2 trình bày cơ sở toán học và kiến trúc thuật toán của phương pháp VMC áp dụng cho hệ nhiều hạt lượng tử:

  • Thiết lập nguyên lý biến phân Rayleigh - Ritz và chứng minh toán học biểu thức năng lượng trung bình $\langle H \rangle \ge E_0$ thông qua phép khai triển hàm sóng thử theo hệ trực giao các hàm riêng của toán tử Hamilton.
  • Ba bước cốt lõi của quy trình VMC: (1) Xây dựng hàm sóng thử nhiều hạt $\Psi(r, \alpha)$ phụ thuộc vào bộ tham số biến phân $\alpha$; (2) Tính toán năng lượng trung bình qua mẫu ngẫu nhiên phân bố xác suất $p(r) = |\Psi|^2 / \int |\Psi|^2 dr$; (3) Quét khoảng biến thiên của các tham số biến phân nhằm tìm cực tiểu của năng lượng và phương sai tương ứng.
  • Cấu trúc toán tử Hamilton: Trong hệ nhiều hạt gồm electron, positron và các ion, Hamilton đầy đủ chứa 7 số hạng: tổng động năng của $N$ electron, động năng của positron, thế năng tương tác hút electron - ion, thế năng tương tác đẩy positron - ion, thế năng tương tác đẩy electron - electron, thế năng tương tác hút positron - electron và thế năng tương tác đẩy giữa các ion.
  • Dạng giải tích của hàm sóng thử: Hàm sóng toàn phần được xây dựng từ tích các hàm sóng đơn hạt và các thừa số tương quan Jastrow. Tác giả so sánh ưu thế giữa hàm sóng dạng Slater và đa thức Laguerre suy rộng (với tính chất trực giao theo từng trạng thái lượng tử), đồng thời lựa chọn thừa số Jastrow dạng Padé do tính chất tối ưu về mặt giải tích và khả năng triệt tiêu điểm kỳ dị Coulomb khi hai hạt tích điện tiến sát nhau.

Chương 3: Áp dụng phương pháp VMC cho phân tử đồng oxit khi có positron

Chương 3 chi tiết hóa việc xây dựng mô hình toán - lý cho phân tử CuO đặt trong hệ tọa độ Descartes ba chiều:

  • Phân tích cấu hình điện tử và mô hình lõi: Nguyên tử oxy ($Z = 8$, cấu hình $1s^2 2s^2 2p^4$) được xem như gồm một lõi trơ (kết hợp hạt nhân và các phân lớp $1s, 2s$) cùng 4 electron hóa trị thuộc phân lớp $2p$ ($n = 2, l = 1$). Nguyên tử đồng ($Z = 29$, cấu hình $1s^2 2s^2 2p^6 3s^2 3p^6 4s^2 3d^9$ hoặc $4s^1 3d^{10}$) được xem như gồm lõi trơ và 2 electron hóa trị thuộc phân lớp $4s$ ($n = 4, l = 0$). Hệ phân tử CuO quy về bài toán chuyển động của 6 electron hóa trị trong trường thế của hai lõi.
  • Thiết lập tọa độ hình học: Gốc tọa độ đặt tại trung điểm đoạn nối hai lõi; lõi oxy đặt tại tọa độ $(-d_1, 0, 0)$, lõi đồng đặt tại $(d_2, 0, 0)$. Véctơ vị trí của electron thứ $i$ là $r_i = (x_i, y_i, z_i)$ và của positron là $r_p = (x_p, y_p, z_p)$.
  • Hàm sóng thử toàn phần: Xây dựng dưới dạng tích của các hàm sóng đơn hạt Slater của 6 electron, hàm sóng đơn hạt của positron, hệ số tương quan electron - electron ($f^{ee}$) và hệ số tương quan positron - electron ($f^{ep}$) theo xấp xỉ Padé. Bộ tham số biến phân bao gồm: $Z_{Cu}, Z_O, \alpha, \beta$ (cho hệ electron) và bổ sung thêm $Z_{pCu}, Z_{pO}, \alpha', \beta'$ (cho hệ có positron).
  • Toán tử năng lượng và các đạo hàm giải tích: Biểu thức động năng của electron và positron được rút gọn thông qua việc lấy Gradient và Laplace của logarit hàm sóng. Năng lượng tương quan electron - positron ($E_{corr}^{e-p}$) được xác định thông qua việc tách các số hạng động năng riêng phần và thế năng tương tác Culông $U^{ep}$.
  • Phương sai năng lượng: Sử dụng đại lượng phương sai $\sigma^2 = \langle E_L^2 \rangle - \langle E_L \rangle^2$ làm tiêu chuẩn hội tụ song song với việc cực tiểu hóa năng lượng kỳ vọng; cấu hình tối ưu đạt được khi phương sai tiến gần về 0.

Chương 4: Các kết quả tính toán

Chương 4 công bố kết quả mô phỏng số thu được từ chương trình viết bằng Visual C++. Điều kiện mô phỏng gồm: cấu hình không gian bằng 100, số bước Monte Carlo là 100.000 bước cho mỗi electron và 10.000 bước cho mỗi positron. Quy trình tối ưu hóa tham số được thực hiện qua các bước biến phân độc lập:

  1. Quá trình biến phân hệ electron trong phân tử CuO:
  • Biến phân $Z_{Cu}$: Khảo sát trong khoảng từ 2,0 đến 4,0 với bước nhảy 0,05; năng lượng tổng $E_0$ và sai số tương đối $\sigma/E_0$ đạt cực tiểu tại giá trị tối ưu $Z_{Cu} = 3{,}3$.
  • Biến phân $\alpha$: Cố định $Z_{Cu} = 3{,}3$, khảo sát $\alpha$ từ 0,10 đến 0,30 với bước nhảy 0,01; điểm tối ưu xác định tại $\alpha = 0{,}13$.
  • Biến phân $\beta$: Sử dụng $Z_{Cu} = 3{,}3$ và $\alpha = 0{,}13$, khảo sát $\beta$ từ 3,0 đến 7,0 với bước nhảy 0,10; điểm tối ưu xác định tại $\beta = 4{,}8$.
  • Biến phân $Z_O$: Sử dụng các tham số tối ưu trên, khảo sát $Z_O$ từ 3,0 đến 7,0 với bước nhảy 0,20; điểm tối ưu xác định tại $Z_O = 5{,}4$.
  • Năng lượng tổng kỳ vọng của hệ electron tại điểm tối ưu đạt giá trị: $E_0 = -49{,}8539\text{ a.u.}$ (tương đương $-1356{,}5\text{ eV}$).
  1. Quá trình biến phân hệ electron - positron trong phân tử CuO:
  • Biến phân $Z_{pCu}$: Khảo sát trong khoảng từ 0,0 đến 0,7 với bước nhảy 0,03; điểm tối ưu đạt tại $Z_{pCu} = 0{,}15$.
  • Biến phân $Z_{pO}$: Khảo sát trong khoảng từ 0,2 đến 1,5 với bước nhảy 0,05; điểm tối ưu đạt tại $Z_{pO} = 0{,}85$.
  • Biến phân $\beta'$: Khảo sát trong khoảng từ -3,0 đến -0,1 với bước nhảy 0,10; điểm tối ưu đạt tại $\beta' = -0{,}6$.
  • Biến phân $\alpha'$: Khảo sát trong khoảng từ 0,4 đến 1,8 với bước nhảy 0,04; điểm tối ưu đạt tại $\alpha' = 0{,}92$.
  • Năng lượng tổng kỳ vọng của hệ có positron tại điểm tối ưu đạt: $E = -48{,}6328\text{ a.u.}$ (tương đương $-1323\text{ eV}$).
  1. Xác định năng lượng tương quan: Nạp toàn bộ bộ 8 tham số tối ưu vào mô đun tính năng lượng tương quan, kết quả xuất ra giá trị năng lượng tương quan electron - positron trong phân tử CuO là $E_{corr}^{e-p} = -5{,}06\text{ eV}$.

Kết quả và đóng góp

Bảng tổng hợp các tham số biến phân tối ưu và kết quả năng lượng của phân tử CuO:

Đại lượng / Tham số Ký hiệu Giá trị tối ưu Đơn vị Ghi chú
Điện tích hiệu dụng lõi đồng với electron $Z_{Cu}$ 3,3 Không thứ nguyên Biến phân từ 2,0 đến 4,0
Tham số tương quan electron - electron (1) $\alpha$ 0,13 Không thứ nguyên Biến phân từ 0,10 đến 0,30
Tham số tương quan electron - electron (2) $\beta$ 4,8 Không thứ nguyên Biến phân từ 3,0 đến 7,0
Điện tích hiệu dụng lõi oxy với electron $Z_O$ 5,4 Không thứ nguyên Biến phân từ 3,0 đến 7,0
Năng lượng trạng thái cơ bản hệ electron $E_0$ -49,8539 a.u. (Hartree) Tương đương -1356,5 eV
Điện tích hiệu dụng lõi đồng với positron $Z_{pCu}$ 0,15 Không thứ nguyên Biến phân từ 0,0 đến 0,7
Điện tích hiệu dụng lõi oxy với positron $Z_{pO}$ 0,85 Không thứ nguyên Biến phân từ 0,2 đến 1,5
Tham số tương quan positron - electron (1) $\beta'$ -0,6 Không thứ nguyên Biến phân từ -3,0 đến -0,1
Tham số tương quan positron - electron (2) $\alpha'$ 0,92 Không thứ nguyên Biến phân từ 0,4 đến 1,8
Năng lượng trạng thái cơ bản hệ e - p $E$ -48,6328 a.u. (Hartree) Tương đương -1323 eV
Năng lượng tương quan positron - electron $E_{corr}^{e-p}$ -5,06 eV Xuất từ chương trình mô phỏng

Về mặt bản chất vật lý, kết quả tính toán cho thấy năng lượng toàn phần của hệ khi có mặt positron ($E = -48{,}6328\text{ a.u.}$) có giá trị đại số lớn hơn năng lượng của hệ electron thuần túy ($E_0 = -49{,}8539\text{ a.u.}$). Tác giả khẳng định quy luật này hoàn toàn phù hợp với thực tế vi mô: positron mang điện tích dương khi đưa vào hệ phân tử sẽ đóng góp các thành phần thế năng đẩy dương đối với các lõi nguyên tử mang điện dương cùng động năng dương của chính nó, làm năng lượng liên kết tổng thể của hệ dịch chuyển về phía kém âm hơn.

Về mặt đóng góp so sánh, giá trị năng lượng tương quan electron - positron $E_{corr}^{e-p} = -5{,}06\text{ eV}$ tính toán được nằm trong dải năng lượng tương quan hợp lý khi đối chiếu với các kết quả nghiên cứu trước đó trên các kim loại và oxit kim loại lân cận (như tính toán trên kẽm oxit ZnO của Hoàng Anh Tuấn năm 2007, hoặc các kim loại cấu trúc mạng fcc/bcc của Trịnh Hoa Lăng và Châu Văn Tạo). Công trình đã hoàn thiện một quy trình tính toán số khép kín bằng ngôn ngữ C++ ứng dụng thuật toán VMC cho một đối tượng phân tử oxit cụ thể có cấu trúc vỏ phức tạp.

Hạn chế và hướng nghiên cứu tiếp

Các hạn chế được tác giả nêu rõ trong văn bản bao gồm:

  • Hạn chế về phạm vi đại lượng khảo sát: Do giới hạn về mặt thời gian, luận văn mới dừng lại ở việc xác định năng lượng tương quan electron - positron, chưa mở rộng tính toán các thông số định lượng quan trọng khác của quá trình hủy positron như tốc độ hủy cặp (annihilation rate) và hệ số tăng cường mật độ (enhancement factor).
  • Giới hạn của mô hình xấp xỉ hàm sóng: Khóa luận mới chỉ thử nghiệm dạng hàm sóng quỹ đạo Slater và Gauss, chưa áp dụng các hàm sóng trực giao cao cấp như đa thức Laguerre suy rộng.
  • Giản lược toán tử Hamilton: Mô hình chưa tính đến các hiệu ứng tương tác tinh vi như tương tác spin - quỹ đạo (spin - orbit coupling) và tương tác từ giữa các electron.
  • Quy mô ô mô phỏng: Chỉ khảo sát 6 electron hóa trị của một phân tử CuO đơn lẻ, chưa tính đến sự đóng góp của các lớp electron bên trong lõi cũng như thế tương tác tuần hoàn của mạng tinh thể chất rắn xung quanh.
  • Thông số mô phỏng Monte Carlo: Số cấu hình không gian (100) và số bước lặp (100.000 bước cho electron và 10.000 bước cho positron) bị giới hạn bởi năng lực xử lý tính toán tại thời điểm thực hiện.

Các kiến nghị mở rộng nghiên cứu tiếp theo:

  • Thử nghiệm các dạng hàm sóng thử và hàm Jastrow khác nhau (như xấp xỉ Mitás, xấp xỉ tâm nguyên tử, xấp xỉ Williamson).
  • Đưa số hạng tương tác spin - quỹ đạo vào toán tử Hamilton để tăng độ chính xác của năng lượng tương quan.
  • Bổ sung số electron ở các lớp vỏ bên trong của nguyên tử và tăng số lượng nguyên tử trong ô mô phỏng nhằm mô tả sát hơn trạng thái tinh thể chất rắn.
  • Tăng số bước mô phỏng Monte Carlo và mở rộng số cấu hình không gian để nâng cao độ hội tụ thống kê và giảm thiểu phương sai.
  • Xây dựng thư viện số chuyên biệt lưu trữ dữ liệu năng lượng toàn phần và năng lượng tương quan electron - positron cho các kim loại và oxit kim loại khác nhau.

Giá trị tham khảo

Báo cáo khóa luận cung cấp tài liệu tham khảo có giá trị chuyên môn đối với:

  • Sinh viên, học viên cao học và nghiên cứu sinh chuyên ngành Vật lý hạt nhân, Vật lý lý thuyết và Vật lý tính toán đang tiếp cận phương pháp biến phân Monte Carlo lượng tử (VMC) cho hệ nhiều hạt.
  • Các nhà nghiên cứu trong lĩnh vực vật lý chất rắn, khoa học vật liệu sử dụng phương pháp phổ thời gian sống của positron (PAS) để khảo sát khuyết tật cấu trúc vi mô trong vật liệu gốm oxit và chất bán dẫn.
  • Nghiên cứu viên trong lĩnh vực y học hạt nhân và kỹ thuật xạ trị, chẩn đoán hình ảnh PET cần tài liệu mô tả chi tiết cơ chế tương tác vi mô, quá trình nhiệt hóa và hủy cặp của positron.

Các nội dung có giá trị tham khảo cụ thể trong văn bản:

  • Chương 2 và Chương 3: Trình bày chi tiết các bước khai triển giải tích của toán tử động năng thông qua Gradient và Laplace của logarit hàm sóng, cung cấp công thức toán học tường minh cho việc áp dụng hàm Jastrow Padé trong hệ nhiều fermion.
  • Chương 4 và Phụ lục mã nguồn: Cung cấp thuật toán tối ưu hóa từng bước trên Visual C++, cách thiết lập cấu hình không gian và vòng lặp tính toán năng lượng cục bộ (hàm eLocal, biến tích lũy eSum, eSqdSum), phục vụ trực tiếp cho việc phát triển các chương trình mô phỏng số lượng tử tương tự.

Câu hỏi thường gặp

Tại sao trong phân tử CuO, tác giả chỉ đưa 6 electron vào mô hình tính toán mà không khảo sát toàn bộ 37 electron của hai nguyên tử? Tác giả áp dụng mô hình xấp xỉ lõi nguyên tử để giảm thiểu độ phức tạp tính toán. Đối với oxy ($Z = 8$, cấu hình $1s^2 2s^2 2p^4$), các electron phân lớp $1s$ và $2s$ đã lấp đầy được ghép chung với hạt nhân tạo thành lõi trơ, chỉ 4 electron thuộc phân lớp $2p$ chưa lấp đầy tham gia liên kết. Đối với đồng ($Z = 29$), các electron thuộc 3 lớp bên trong kết hợp với hạt nhân tạo thành lõi trơ, chỉ có 2 electron thuộc phân lớp ngoài $4s$ được tính đến (nhằm giải thích hóa trị II của đồng). Do đó, hệ quy về 6 electron hóa trị chuyển động trong trường hai lõi.

Điều kiện năng lượng nào quyết định việc positron liên kết với nguyên tử hay tạo thành hạt Positronium (Ps) tự do? Theo văn bản, năng lượng liên kết của Positronium là $0{,}250\text{ Hartree}$. Khi thế ion của nguyên tử thỏa mãn điều kiện $I < 0{,}250\text{ Hartree}$, positron chỉ liên kết bền vững với nguyên tử nếu năng lượng liên kết giữa positron và nguyên tử đó lớn hơn hiệu số $(0{,}250 - I)$. Ngược lại, nếu năng lượng liên kết nhỏ hơn hiệu số này, hệ positron - nguyên tử sẽ tách thành hạt Positronium tự do và một ion dương.

Vì sao năng lượng toàn phần của hệ khi có positron ($-48{,}6328\text{ a.u.}$) lại cao hơn năng lượng của hệ khi không có positron ($-49{,}8539\text{ a.u.}$)? Do positron mang điện tích dương ($+e$). Khi positron thâm nhập vào phân tử CuO, hệ xuất hiện thêm thế năng đẩy tĩnh điện Coulomb giữa positron với hai lõi nguyên tử mang điện dương của đồng và oxy, đồng thời hệ phải cộng thêm động năng dương của bản thân positron. Mặc dù có thế năng hút giữa positron và các electron, tổng các đóng góp năng lượng dương làm cho giá trị đại số của năng lượng toàn phần tăng lên (nghĩa là giá trị tuyệt đối của năng lượng âm giảm đi).

Cơ sở khoa học nào cho phép giả định positron đã bị nhiệt hóa hoàn toàn trước khi xảy ra quá trình hủy cặp trong chất rắn? Dựa trên kết quả nghiên cứu của Hautojärvi được trích dẫn trong văn bản: một positron có năng lượng ban đầu khoảng $2\text{ MeV}$ khi đi vào môi trường sẽ liên tục va chạm với electron, ion mạng và phonon, mất dần năng lượng và đạt trạng thái cân bằng nhiệt (positron nhiệt) trong khoảng thời gian nhỏ hơn $20\text{ ps}$. Trong khi đó, thời gian sống trung bình của positron trong kim loại kéo dài khoảng $200\text{ ps}$ (lớn hơn gấp 10 lần). Khoảng chênh lệch thời gian này khẳng định quá trình nhiệt hóa luôn hoàn tất trước khi positron bị hủy.

Các tham số biến phân tối ưu của electron và positron trong phân tử CuO được xác định bằng những giá trị cụ thể nào? Qua quá trình quét cực tiểu hóa năng lượng và phương sai, bộ tham số tối ưu của hệ electron là $Z_{Cu} = 3{,}3$; $\alpha = 0{,}13$; $\beta = 4{,}8$; $Z_O = 5{,}4$ với năng lượng $E_0 = -49{,}8539\text{ a.u.}$. Bộ tham số tối ưu liên quan đến positron là $Z_{pCu} = 0{,}15$; $Z_{pO} = 0{,}85$; $\beta' = -0{,}6$; $\alpha' = 0{,}92$ với năng lượng $E = -48{,}6328\text{ a.u.}$. Năng lượng tương quan electron - positron tương ứng thu được là $E_{corr}^{e-p} = -5{,}06\text{ eV}$.

Kết luận

Khóa luận tốt nghiệp của tác giả Phạm Thị Phú đã áp dụng thành công phương pháp biến phân Monte Carlo lượng tử (VMC) trong gần đúng Born - Oppenheimer để giải quyết bài toán tương tác đa hạt trong phân tử đồng oxit (CuO). Bằng việc thiết lập hàm sóng thử dạng Slater kết hợp hàm tương quan Jastrow Padé và giải thuật toán số trên Visual C++, đề tài đã xác định được bộ tham số biến phân tối ưu cùng năng lượng trạng thái cơ bản của hệ khi có và không có positron. Kết quả tính toán đạt được năng lượng tương quan electron - positron $E_{corr}^{e-p} = -5{,}06\text{ eV}$, phản ánh quy luật vật lý phù hợp với dải số liệu thực nghiệm trong các oxit kim loại chuyển tiếp. Công trình đóng góp một tài liệu học thuật có giá trị tham khảo rõ ràng về mặt phương pháp luận tính toán số lượng tử và mô hình hóa tương tác phản hạt trong vật lý hạt nhân.