Luận Văn Thạc Sĩ: Tính Toán Phổ Dao Động Của D-Glucose Bằng Phương Pháp DFT

Luận văn thạc sĩ trình bày phương pháp DFT để tính toán phổ dao động của D-glucose, mang lại cái nhìn sâu sắc về cấu trúc phân tử.

• Tác giả: Nguyễn Thị Thủy
• Người hướng dẫn: TS. Hoàng Chí Hiếu
• Trường học: Đại học Quốc gia Hà Nội
• Chuyên ngành: Quang học
• Đề tài: Tính Toán Phổ Dao Động Của D-Glucose Bằng Phương Pháp DFT
• Loại tài liệu: luận văn thạc sĩ
• Năm xuất bản: 2015
• Địa điểm: Hà Nội

Mục lục chi tiết

• LỜI CẢM ƠN
• MỞ ĐẦU
• 1. CHƯƠNG I: TỔNG QUAN VỀ PHƯƠNG PHÁP DFT
• 1.1. Giới thiệu về lý thuyết phiếm hàm mật độ (DFT)
• 1.2. Bài toán của hệ nhiều hạt
• 1.3. Ý tưởng ban đầu về DFT: Thomas-Fermi và các mô hình liên quan
• 1.4. Các định lý Hohenberg-Kohn
• 1.5. Giới thiệu về orbital và hàm năng lượng Kohn-Sham
• 1.6. Phiếm hàm gần đúng mật độ địa phương (LDA – Local Density Approximation). Phương pháp gần đúng gradient suy rộng (GGA)
• 1.7. Mô hình lý thuyết phiếm hàm mật độ trong Dmol3. Nguồn gốc và cấu trúc phổ Raman
• 1.8. Các nguyên tắc chọn lọc cho phổ Hồng ngoại và phổ Raman
• 1.9. Sự dao động của phân tử 2 nguyên tử
• 1.10. So sánh phổ Raman và phổ Hồng ngoại
• 1.11. Ứng dụng của phương pháp phân tích phổ Raman
• 2. CHƯƠNG II: PHƯƠNG PHÁP NGHIÊN CỨU VÀ XÂY DỰNG MÔ HÌNH TÍNH TOÁN
• 2.1. Tổng quan về Glucose. Các phân tử Sacchride
• 2.2. Trạng thái tự nhiên
• 2.3. Tính chất Vật lý
• 2.4. Công thức cấu tạo
• 2.5. Tính chất hóa học
• 2.6. Điều chế và ứng dụng
• 2.7. Tầm quan trọng của Glucose trong đời sống
• 2.8. Một số nghiên cứu quang phổ học dao động của Glucose
• 2.9. Phương pháp tính toán. Mô hình và các thông số tính toán
• 2.9.1. Cấu trúc phân tử của D-Glucose
• 2.9.2. Cấu trúc phân tử của H2O
• 2.9.3. Xây dựng mô hình D-Glucose có thêm phân tử H2O để nghiên cứu ảnh hưởng của H2O lên phổ Raman của D-Glucose
• 3. CHƯƠNG III: KẾT QUẢ VÀ THẢO LUẬN
• 3.1. Đặc trưng phổ Raman của và D-Glucose
• 3.1.1. Vùng phổ từ 0 – 1700 cm-1
• 3.1.2. Vùng phổ từ 2700–3200cm-1
• 3.1.3. Vùng phổ từ 3300–3900cm-1
• 3.2. Đặc trưng phổ Raman của phân tử H2O
• 3.3. Ảnh hưởng của H2O lên phổ Raman của D-Glucose
• 3.3.1. Ảnh hưởng của H2O lên phổ Raman của D-Glucose trong khoảng tần số lân cận 1600cm-1
• 3.3.2. Ảnh hưởng của H2O lên phổ Raman của D-Glucose trong khoảng tần số 3600cm-1 – 3800cm-1
• 4. CHƯƠNG IV: KẾT LUẬN
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Về Tính Toán Phổ Dao Động D Glucose DFT

Nghiên cứu tính toán phổ dao động đóng vai trò quan trọng trong việc xác định cấu trúc và tính chất của các phân tử. Trong đó, D-Glucose, một monosaccharide phổ biến, là đối tượng được quan tâm đặc biệt. Việc sử dụng phương pháp DFT (Density Functional Theory) cho phép mô phỏng và dự đoán phổ dao động của D-Glucose một cách chính xác. Luận văn này tập trung vào việc áp dụng DFT để tính toán phổ dao động của D-Glucose, từ đó làm sáng tỏ các kết quả thực nghiệm và hiểu rõ hơn về cấu trúc phân tử của nó. DFT là một lý thuyết lượng tử được sử 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.

1.1. Giới thiệu về đường D Glucose và ứng dụng thực tiễn

D-Glucose là một monosaccharide quan trọng, đóng vai trò thiết yếu trong nhiều quá trình sinh học và công nghiệp. Nó là nguồn năng lượng chính cho tế bào, tham gia vào cấu trúc của ARN và ADN, và được sử dụng rộng rãi trong thực phẩm và dược phẩm. Nghiên cứu về D-Glucose giúp chúng ta hiểu rõ hơn về các quá trình trao đổi chất và phát triển các ứng dụng mới trong y học và công nghiệp thực phẩm. Với tầm quan trọng như vậy Glucose đã và đang là đối tượng nghiên cứu rộng rãi của rất nhiều nhà khoa học trong nước và trên thế giới.

1.2. Tổng quan về phương pháp tính toán DFT trong hóa học

Phương pháp DFT là một công cụ mạnh mẽ để tính toán cấu trúc điện tử và phổ dao động của các phân tử. DFT dựa trên việc sử dụng mật độ điện tử thay vì hàm sóng, giúp giảm đáng kể chi phí tính toán. Phương pháp này đã được chứng minh là có độ chính xác cao trong nhiều ứng dụng, bao gồm cả việc mô phỏng phổ IR và phổ Raman. Trong những năm gần đây, việc sử dụng các phần mềm hỗ trợ việc mô phỏng và tính toán các tính chất và quá trình biến đổi bên trong các vật chất được quan tâm nhiều.

II. Thách Thức Trong Tính Toán Phổ Dao Động D Glucose DFT

Việc tính toán phổ dao động của D-Glucose bằng phương pháp DFT gặp phải một số thách thức nhất định. Cấu trúc phân tử phức tạp của D-Glucose đòi hỏi các phương pháp tính toán chính xác và tốn kém về mặt thời gian. Ngoài ra, việc lựa chọn phiếm hàm và cơ sở phù hợp cũng ảnh hưởng đáng kể đến độ chính xác của kết quả. Cần có sự cân nhắc kỹ lưỡng để đảm bảo rằng các tính toán DFT có thể tái tạo chính xác các đặc trưng phổ thực nghiệm. Chúng ta đã biết đến việc nghiên cứu các vật liệu bằng phương pháp phổ dao động như: phổ dao động FT-IR, phổ Raman, phổ tần số tổng…Tuy nhiên tất cả các phương pháp đó rất khó để xác định chính xác được các mode dao động cụ thể hoặc có sự nhầm lần giữa những cộng hưởng Fermi với các mode động thực tế.

2.1. Ảnh hưởng của cấu trúc phân tử D Glucose đến tính toán

Cấu trúc vòng của D-Glucose và sự tồn tại của nhiều nhóm hydroxyl (-OH) tạo ra nhiều khả năng tương tác hydro, làm phức tạp việc tối ưu hóa cấu trúc và tính toán năng lượng. Các cấu trúc đồng phân khác nhau (α và β) cũng cần được xem xét để có được bức tranh đầy đủ về phổ dao động của D-Glucose. Cần có tính toán lý thuyết để làm sáng tỏ những kết luận trên.

2.2. Lựa chọn phiếm hàm và cơ sở cho tính toán DFT chính xác

Việc lựa chọn phiếm hàm (ví dụ: DFT/B3LYP, DFT/PBE) và cơ sở ảnh hưởng lớn đến độ chính xác của tính toán DFT. Các phiếm hàm khác nhau có thể cho kết quả khác nhau về tần số dao động và cường độ phổ. Cần thực hiện các kiểm tra và so sánh với dữ liệu thực nghiệm để lựa chọn phiếm hàm và cơ sở phù hợp nhất cho D-Glucose. Trong những năm gần đây, nhiều nghiên cứu về D-Glucose nói riêng và các phân tử saccharide nói chung đã được công bố.

III. Phương Pháp Tính Toán Phổ Dao Động D Glucose Bằng DFT

Để tính toán phổ dao động của D-Glucose bằng phương pháp DFT, cần thực hiện các bước sau: tối ưu hóa cấu trúc, tính toán tần số dao động, và phân tích phương thức dao động chuẩn. Quá trình này thường được thực hiện bằng các phần mềm chuyên dụng như Gaussian, VASP, hoặc Quantum ESPRESSO. Kết quả tính toán cho phép xác định các chế độ dao động đặc trưng của D-Glucose và so sánh với dữ liệu thực nghiệm. Sử dụng phần mềm DMol3 có thể dự đoán được các quá trình xảy ra ở cấp độ phân tử, các cấu trúc tinh thể của vật rắn và các tương tác bề mặt.

3.1. Tối ưu hóa cấu trúc phân tử D Glucose bằng DFT

Tối ưu hóa cấu trúc là bước đầu tiên và quan trọng nhất trong tính toán DFT. Mục tiêu là tìm ra cấu trúc phân tử có năng lượng thấp nhất, tương ứng với cấu trúc ổn định nhất của D-Glucose. Quá trình này đòi hỏi sự hội tụ của các lực và năng lượng, đảm bảo rằng cấu trúc đã được tối ưu hóa hoàn toàn. Đồ thị năng lượng của quá trình tối ưu hóa mô hình D-Glucose.

3.2. Tính toán tần số dao động và phân tích phương thức chuẩn

Sau khi tối ưu hóa cấu trúc, bước tiếp theo là tính toán tần số dao động. Quá trình này cho phép xác định các tần số dao động của các liên kết và nhóm chức trong phân tử D-Glucose. Phân tích phương thức dao động chuẩn giúp xác định loại dao động (ví dụ: kéo giãn, uốn) tương ứng với mỗi tần số. Khi biết được hàm sóng của các điện tử chương trình sẽ cho ta biết được các thông tin khác như mật độ điện tử, mật độ điện tích, tính chất quang, phổ phát quang, phổ dao động, năng lượng của phản ứng, mật độ trạng thái và cấu trúc vùng năng lượng…

IV. Ảnh Hưởng Của Dung Môi Đến Phổ Dao Động D Glucose DFT

Sự hiện diện của dung môi có thể ảnh hưởng đáng kể đến phổ dao động của D-Glucose. Các tương tác giữa D-Glucose và phân tử dung môi (ví dụ: nước) có thể làm thay đổi tần số dao động và cường độ phổ. Việc sử dụng các mô hình dung môi như Polarizable Continuum Model (PCM) cho phép mô phỏng ảnh hưởng của dung môi đến phổ dao động một cách chính xác. Những bằng chứng về sự chuyển hóa một phần từ cấu trúc của D-Glucose sang D-Glucose khi có ảnh hưởng của thành phần H2O lên cấu trúc phân tử D-Glucose đã được phát hiện bằng phương pháp SFG [18].

4.1. Mô phỏng ảnh hưởng của nước bằng mô hình PCM trong DFT

Mô hình PCM là một phương pháp phổ biến để mô phỏng ảnh hưởng của dung môi đến tính toán DFT. PCM coi dung môi như một môi trường liên tục có hằng số điện môi xác định, cho phép tính toán các tương tác tĩnh điện giữa D-Glucose và dung môi. Các vị trí đặt phân tử H2O vào mô hình đã được tối ưu của D- Glucose.

4.2. So sánh phổ dao động D Glucose trong pha khí và dung dịch

So sánh phổ dao động của D-Glucose trong pha khí và dung dịch giúp làm rõ ảnh hưởng của dung môi đến các chế độ dao động của phân tử. Sự thay đổi về tần số dao động và cường độ phổ có thể cung cấp thông tin về các tương tác giữa D-Glucose và dung môi. Phổ Raman thu được của D-Glucose trong khoảng lân cận 1600cm-1 a. Khi có 1 H2O ở vị trí 1 c. Khi có 1 H2O ở vị trí 2 d.

V. Ứng Dụng Kết Quả Tính Toán Phổ Dao Động D Glucose DFT

Kết quả tính toán phổ dao động của D-Glucose bằng phương pháp DFT có nhiều ứng dụng quan trọng. Chúng có thể được sử dụng để giải thích và xác nhận các kết quả thực nghiệm từ phổ IR và phổ Raman, cũng như để nghiên cứu các tương tác giữa D-Glucose và các phân tử khác. Ngoài ra, kết quả tính toán có thể được sử dụng để phát triển các phương pháp phân tích và nhận dạng D-Glucose dựa trên tính chất quang phổ. Với mục đích dùng kết quả thu được từ tính toán lý thuyết để làm sáng tỏ những kết quả thực nghiệm trước đó tôi đã thực hiện đề tài: “Tính toán phổ dao động của D-Glucose bằng phương pháp DFT”.

5.1. So sánh kết quả DFT với phổ IR và Raman thực nghiệm

So sánh kết quả tính toán DFT với phổ IR và phổ Raman thực nghiệm là một cách quan trọng để đánh giá độ chính xác của phương pháp tính toán. Sự phù hợp giữa kết quả tính toán và thực nghiệm cho thấy rằng phương pháp DFT đã mô tả chính xác các chế độ dao động của D-Glucose. Phổ Raman của α-D-Glucose trong khoảng tần số từ 0–4000 cm-1.

5.2. Ứng dụng trong phân tích và nhận dạng D Glucose

Kết quả tính toán phổ dao động có thể được sử dụng để phát triển các phương pháp phân tích và nhận dạng D-Glucose dựa trên tính chất quang phổ. Ví dụ, có thể sử dụng phổ Raman tính toán để xác định sự có mặt của D-Glucose trong một mẫu phức tạp. Trong công nghiệp thực phẩm Glucose được sử dụng làm chất bảo quản. Trong y học, Glucose là nguồn năng lượng chủ yếu và trực tiếp của cơ thể, được dự trữ ở gan dưới dạng glycogen.

VI. Kết Luận Và Hướng Phát Triển Tính Toán Phổ D Glucose DFT

Việc tính toán phổ dao động của D-Glucose bằng phương pháp DFT là một công cụ mạnh mẽ để nghiên cứu cấu trúc và tính chất của phân tử này. Mặc dù có những thách thức nhất định, nhưng với sự phát triển của các phương pháp tính toán và phần mềm, độ chính xác và hiệu quả của tính toán DFT ngày càng được nâng cao. Trong tương lai, có thể kỳ vọng rằng tính toán DFT sẽ đóng vai trò quan trọng hơn nữa trong việc nghiên cứu D-Glucose và các phân tử sinh học khác. Vì vậy, chúng tôi đã sử dụng phương pháp DFT để tính toán phổ dao động của một số vật liệu nhằm giải thích các kết quả thực nghiệm đã thu được trước đó.

6.1. Tóm tắt kết quả và đánh giá độ chính xác của DFT

Luận văn này đã trình bày các kết quả tính toán phổ dao động của D-Glucose bằng phương pháp DFT. So sánh với dữ liệu thực nghiệm cho thấy rằng DFT có thể mô tả chính xác các chế độ dao động chính của D-Glucose, mặc dù vẫn còn một số sai lệch nhỏ. Trên thế giới có nhiều nghiên cứu đã chứng tỏ rằng về việc sử dụng phần mềm DMol3 tính toán dựa trên lý thuyết DFT có độ chính xác cao trong việc nghiên 1 z cứu phổ dao động của các vật liệu [19].

6.2. Hướng phát triển và cải tiến phương pháp tính toán DFT

Trong tương lai, có thể cải tiến phương pháp tính toán DFT bằng cách sử dụng các phiếm hàm và cơ sở tiên tiến hơn, cũng như bằng cách kết hợp với các phương pháp khác như Molecular Dynamics để mô phỏng các hiệu ứng nhiệt độ và dung môi một cách chính xác hơn. Thành phần tham gia vào cấu trúc của tế bào (ARN và ADN) và một số chất đặc biệt khác [38]

Nội dung chính

Tổng quan nghiên cứu

Glucose là một monosaccharide phổ biến và đóng vai trò thiết yếu trong sinh học và công nghiệp. Theo ước tính, glucose là nguồn năng lượng chủ yếu cho tế bào, tham gia cấu trúc của ARN, ADN và nhiều hợp chất sinh học quan trọng khác. Trong công nghiệp thực phẩm, glucose được sử dụng làm chất bảo quản, còn trong y học, nó là nguồn năng lượng dự trữ dưới dạng glycogen tại gan. Nghiên cứu phổ dao động của glucose giúp hiểu rõ hơn về cấu trúc phân tử và các tương tác hóa học, từ đó hỗ trợ phát triển các ứng dụng trong sinh học, y học và công nghiệp.

Vấn đề nghiên cứu tập trung vào việc tính toán phổ dao động của D-Glucose bằng phương pháp lý thuyết phiếm hàm mật độ (DFT), nhằm làm sáng tỏ các kết quả thực nghiệm phổ Raman đã thu được. Mục tiêu cụ thể là xây dựng mô hình tính toán cấu trúc phân tử D-Glucose và ảnh hưởng của phân tử nước (H2O) lên phổ dao động của nó, trong phạm vi nghiên cứu tại Đại học Quốc gia Hà Nội, năm 2015.

Nghiên cứu có ý nghĩa quan trọng trong việc cung cấp dữ liệu lý thuyết chính xác để giải thích các hiện tượng phổ dao động, hỗ trợ phát triển các kỹ thuật phân tích phổ trong nghiên cứu vật liệu sinh học và hóa học phân tử. Các chỉ số phổ Raman được tính toán trong khoảng tần số từ 0 đến 4000 cm⁻¹, bao gồm các vùng phổ đặc trưng như 0–1700 cm⁻¹, 2700–3200 cm⁻¹ và 3300–3900 cm⁻¹, giúp phân tích chi tiết các mode dao động của phân tử.

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

Khung lý thuyết áp dụng

Lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT) là nền tảng chính của nghiên cứu, cho phép mô tả các tính chất của hệ electron trong phân tử dựa trên mật độ điện tử thay vì hàm sóng phức tạp. DFT dựa trên hai định lý Hohenberg-Kohn, khẳng định năng lượng trạng thái cơ bản là một phiếm hàm của mật độ điện tử, và nguyên lý biến phân để tìm mật độ điện tử tối ưu.

Phương pháp Kohn-Sham được sử dụng để tính toán động năng của hệ không tương tác, trong khi năng lượng tương quan trao đổi được xấp xỉ bằng các phiếm hàm như LDA (Local Density Approximation) và GGA (Generalized Gradient Approximation). Trong luận văn, các phiếm hàm GGA phổ biến như PW91, PBE, BLYP được áp dụng để tăng độ chính xác.

Ngoài ra, mô hình phân tử D-Glucose được xây dựng dựa trên cấu trúc dạng mạch vòng, kết hợp với phân tử H2O để nghiên cứu ảnh hưởng của nước lên phổ dao động. Các khái niệm chính bao gồm phổ Raman, phổ hồng ngoại, mode dao động chuẩn tắc, tensor phân cực và các nguyên tắc chọn lọc phổ.

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

Nguồn dữ liệu chính là các mô hình phân tử D-Glucose và H2O được tối ưu hóa bằng phần mềm DMol3 trong bộ công cụ Materials Studio, sử dụng phương pháp DFT với hệ cơ sở DNP (Double Numerical plus Polarization). Cỡ mẫu là các phân tử đơn lẻ và các tổ hợp D-Glucose với một hoặc hai phân tử H2O.

Phương pháp phân tích bao gồm tính toán phổ Raman trong các vùng tần số từ 0 đến 4000 cm⁻¹, phân tích các mode dao động chuẩn tắc và so sánh phổ của D-Glucose khô và khi có sự hiện diện của H2O. Timeline nghiên cứu kéo dài trong năm 2015, với các bước chính: xây dựng mô hình, tối ưu hóa cấu trúc, tính toán phổ dao động, và phân tích kết quả.

Phương pháp chọn mẫu là lựa chọn các cấu hình phân tử tiêu biểu, trong đó vị trí phân tử H2O được đặt tại các vị trí khác nhau để khảo sát ảnh hưởng lên phổ Raman của D-Glucose. Phân tích dữ liệu được thực hiện bằng cách so sánh các đỉnh phổ, sự dịch chuyển số sóng và cường độ phổ.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

1. Đặc trưng phổ Raman của D-Glucose: Phổ Raman của D-Glucose được tính toán trong vùng 0–1700 cm⁻¹ cho thấy các mode dao động hóa trị và biến dạng đặc trưng với các đỉnh phổ rõ ràng. Ví dụ, các mode dao động trong vùng 2700–3200 cm⁻¹ liên quan đến dao động kéo dài liên kết C-H, với cường độ phổ chiếm khoảng 70% tổng cường độ phổ trong vùng này.

2. Ảnh hưởng của phân tử H2O lên phổ Raman: Khi thêm một hoặc hai phân tử H2O vào mô hình, phổ Raman của D-Glucose có sự dịch chuyển số sóng đáng kể, đặc biệt trong vùng tần số lân cận 1600 cm⁻¹ và 3600–3800 cm⁻¹. Sự dịch chuyển này dao động trong khoảng 5–15 cm⁻¹ tùy vị trí phân tử H2O, cho thấy tương tác hydro làm thay đổi cấu trúc phân tử và đặc tính dao động.

3. Phổ Raman của phân tử H2O: Phổ Raman của H2O được tính toán cho thấy các mode dao động chuẩn tắc với tần số đặc trưng ở khoảng 1600 cm⁻¹ (dao động biến dạng) và 3600–3800 cm⁻¹ (dao động kéo dài O-H), phù hợp với dữ liệu thực nghiệm.

4. So sánh phổ Raman của D-Glucose khô và ướt: Kết quả tính toán cho thấy phổ Raman của D-Glucose khi có phân tử H2O có sự thay đổi rõ rệt về cường độ và vị trí đỉnh phổ so với mẫu khô, tương ứng với các kết quả thực nghiệm đã được báo cáo trong ngành.

Thảo luận kết quả

Nguyên nhân của sự dịch chuyển phổ Raman khi có mặt H2O là do tương tác hydro giữa phân tử nước và các nhóm hydroxyl của D-Glucose, làm thay đổi mật độ điện tử và cấu trúc liên kết hóa học. Kết quả này phù hợp với các nghiên cứu thực nghiệm sử dụng phổ SFG và FT-Raman, chứng minh tính chính xác của phương pháp DFT trong mô phỏng phổ dao động.

So sánh với các nghiên cứu khác, phương pháp DFT với phiếm hàm GGA và hệ cơ sở DNP cho độ chính xác cao trong việc dự đoán vị trí và cường độ các mode dao động. Dữ liệu có thể được trình bày qua biểu đồ phổ Raman so sánh giữa các mô hình, hoặc bảng tổng hợp sự dịch chuyển số sóng và cường độ tương ứng.

Ý nghĩa của kết quả là cung cấp cơ sở lý thuyết vững chắc để giải thích các hiện tượng phổ dao động trong các hệ phân tử sinh học, hỗ trợ phát triển các kỹ thuật phân tích phổ trong nghiên cứu vật liệu và sinh học phân tử.

Đề xuất và khuyến nghị

1. Mở rộng nghiên cứu với các phân tử saccharide khác: Thực hiện tính toán phổ dao động cho các phân tử disaccharide và polysaccharide để hiểu rõ hơn về ảnh hưởng cấu trúc lên phổ Raman, nhằm nâng cao độ bao phủ của dữ liệu phổ trong sinh học phân tử. Thời gian thực hiện dự kiến 1-2 năm, do các nhóm nghiên cứu hóa học phân tử và vật lý lý thuyết đảm nhận.

2. Phát triển mô hình tương tác phân tử nước đa phân tử: Xây dựng mô hình tính toán bao gồm nhiều phân tử H2O để mô phỏng môi trường dung dịch thực tế, giúp cải thiện độ chính xác của phổ dao động trong điều kiện sinh học. Khuyến nghị sử dụng các phương pháp mô phỏng động học phân tử kết hợp DFT, thực hiện trong vòng 1 năm.

3. Ứng dụng phổ Raman tính toán trong phân tích thực nghiệm: Kết hợp dữ liệu tính toán với phổ Raman thực nghiệm để xác định chính xác các mode dao động, giảm thiểu nhầm lẫn do cộng hưởng Fermi hoặc các hiệu ứng phức tạp khác. Đề xuất triển khai trong các phòng thí nghiệm quang phổ và nghiên cứu vật liệu sinh học.

4. Phát triển phần mềm hỗ trợ phân tích phổ dựa trên DFT: Tạo công cụ tính toán và phân tích phổ dao động tự động, giúp các nhà nghiên cứu dễ dàng áp dụng phương pháp DFT trong phân tích phổ Raman và IR. Thời gian phát triển dự kiến 1-2 năm, phối hợp giữa các chuyên gia tin học và hóa lý.

Đối tượng nên tham khảo luận văn

1. Nhà nghiên cứu vật lý lý thuyết và hóa lượng tử: Luận văn cung cấp phương pháp và kết quả tính toán phổ dao động bằng DFT, giúp họ phát triển các mô hình tính toán và áp dụng vào nghiên cứu vật liệu phân tử.

2. Chuyên gia quang phổ học và phân tích vật liệu: Thông tin chi tiết về phổ Raman và ảnh hưởng của môi trường lên phổ dao động hỗ trợ họ trong việc giải thích dữ liệu thực nghiệm và thiết kế thí nghiệm.

3. Nhà khoa học sinh học phân tử và hóa sinh: Hiểu rõ cấu trúc và dao động của glucose giúp nghiên cứu các quá trình sinh học liên quan đến carbohydrate, như trao đổi năng lượng và tương tác phân tử.

4. Ngành công nghiệp thực phẩm và dược phẩm: Dữ liệu về phổ dao động hỗ trợ kiểm soát chất lượng sản phẩm chứa glucose và các dẫn xuất, cũng như phát triển các sản phẩm mới dựa trên tính chất phân tử.

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

1. Phương pháp DFT có ưu điểm gì trong tính toán phổ dao động?
   DFT sử dụng mật độ điện tử thay vì hàm sóng phức tạp, giúp giảm đáng kể độ phức tạp tính toán mà vẫn đảm bảo độ chính xác cao. Ví dụ, phần mềm DMol3 dựa trên DFT cho phép mô phỏng phổ Raman với độ tin cậy cao, phù hợp với các hệ phân tử lớn như glucose.

2. Tại sao cần tính toán ảnh hưởng của phân tử nước lên phổ dao động của glucose?
   Phân tử nước tương tác hydro với glucose, làm thay đổi cấu trúc và mật độ điện tử, dẫn đến sự dịch chuyển và biến đổi cường độ các mode dao động trong phổ Raman. Điều này giúp giải thích các kết quả thực nghiệm trong môi trường dung dịch.

3. Phổ Raman và phổ hồng ngoại khác nhau như thế nào?
   Phổ Raman dựa trên sự thay đổi độ phân cực của phân tử trong quá trình dao động, còn phổ hồng ngoại dựa trên sự thay đổi moment lưỡng cực. Do đó, một số dao động chỉ xuất hiện trong phổ Raman hoặc phổ hồng ngoại, cung cấp thông tin bổ sung cho nhau.

4. Hệ cơ sở DNP trong DMol3 có vai trò gì?
   Hệ cơ sở DNP (Double Numerical plus Polarization) cung cấp các hàm cơ sở số nhân đôi kèm theo các hàm phân cực, giúp mô phỏng chính xác hơn các orbital phân tử và tương tác electron, từ đó nâng cao độ chính xác của tính toán phổ dao động.

5. Kết quả tính toán có thể áp dụng trong thực tế như thế nào?
   Kết quả giúp giải thích các phổ Raman thực nghiệm, hỗ trợ phân tích cấu trúc phân tử và tương tác trong các hệ sinh học và công nghiệp. Ví dụ, trong công nghiệp thực phẩm, dữ liệu này giúp kiểm soát chất lượng sản phẩm chứa glucose và các hợp chất liên quan.

Kết luận

• Luận văn đã thành công trong việc sử dụng phương pháp DFT để tính toán phổ dao động của D-Glucose, bao gồm ảnh hưởng của phân tử H2O lên phổ Raman.
• Kết quả tính toán phù hợp với các dữ liệu thực nghiệm, chứng minh tính hiệu quả của phương pháp DFT và phần mềm DMol3 trong nghiên cứu phổ dao động phân tử.
• Phổ dao động được phân tích chi tiết trong các vùng tần số 0–1700 cm⁻¹, 2700–3200 cm⁻¹ và 3300–3900 cm⁻¹, cung cấp thông tin về các mode dao động hóa trị và biến dạng.
• Nghiên cứu mở ra hướng phát triển các mô hình tính toán phổ dao động cho các phân tử sinh học phức tạp hơn và trong môi trường dung dịch thực tế.
• Đề xuất các bước tiếp theo bao gồm mở rộng nghiên cứu sang các phân tử saccharide khác, phát triển mô hình tương tác đa phân tử nước, và ứng dụng kết quả trong phân tích phổ thực nghiệm.

Để tiếp tục phát triển nghiên cứu, các nhà khoa học và chuyên gia trong lĩnh vực quang phổ, hóa lượng tử và sinh học phân tử được khuyến khích áp dụng phương pháp và kết quả luận văn nhằm nâng cao hiểu biết về cấu trúc và tính chất phân tử trong các hệ phức tạp.

Trích đoạn nội dung tài liệu

CHƯƠNG I TỔNG QUAN VỀ PHƯƠNG PHÁP DFT 1. Giới thiệu về lý thuyết phiếm hàm mật độ (DFT) Trong cơ học lượng tử, để nghiên cứu hệ có N điện tử chúng ta phải đi giải phương trình Schrödinger để tìm ra hàm sóng  của hệ là hàm của 3N biến số. Cho đến hiện nay, chúng ta chỉ có lời giải chính xác đối với trường hợp nguyên tử hyđro (bài toán 1 điện tử, N = 1). Đối với phân tử hyđro chúng ta chỉ có thể giải gần đúng phương trình Schrödinger.Về mặt giải tích, hiện tại chưa có phương pháp nào giải được chính xác phương trình Schrödinger của hệ nhiều điện tử.

Lý thuyết phiếm hàm mật độ (Density-functional Theory, DFT) là một cách tiếp cận khác mà có thể hiện thực hóa việc nghiên cứu các hệ nhiều hạt. DFT là một lý thuyết hiện đại dựa trên nền tảng của cơ học lượng tử. DFT có thể được dùng để mô tả các tính chất của hệ điện tử trong nguyên tử, phân tử, vật rắn… Điểm cốt yếu trong lý thuyết này là các tính chất của hệ N điện tử được biểu diễn thông qua hàm mật độ điện tử của hệ (là hàm của ba biến tọa độ không gian) thay vì hàm sóng của 3N biến tọa độ không gian trong cơ học lượng tử. Vì vậy, DFT có ưu điểm lớn (và hiện nay đang được sử dụng nhiều nhất) trong việc nghiên cứu các tính chất của các hệ vật liệu từ nguyên tử, phân tử cho tới chất rắn… Ý tưởng dùng hàm mật độ điện tử để mô tả các tính chất của hệ điện tử được nêu trong các công trình của Llewellyn Hilleth Thomas và Enrico Fermi ngay từ khi cơ học lượng tử mới ra đời.Đến năm 1964, Pierre Hohenberg và Walter Kohn đã chứng minh chặt chẽ hai định lý cơ bản là nền tảng của lý thuyết phiếm hàm mật độ.Hai định lý khẳng định năng lượng ở trạng thái cơ bản là một phiếm hàm của mật độ điện tử, do đó về nguyên tắc có thể mô tả hầu hết các tính chất vật lý của hệ điện tử qua hàm mật độ điện tử.

Một năm sau, Walter Kohn và Lu Jeu Sham nêu ra qui trình tính toán để thu được gần đúng mật độ điện tử ở trạng thái cơ bản trong khuôn khổ lý thuyết DFT. Từ những năm 1980 đến nay, cùng với sự phát triển tốc 3 z độ tính toán của máy tính điện tử, lý thuyết DFT được sử dụng rộng rãi và hiệu quả trong các ngành khoa học như: vật lý chất rắn, hóa học lượng tử, vật lý sinh học, khoa học vật liệu… Walter Kohn đã được ghi nhận những đóng góp của ông cho việc phát triển lý thuyết phiếm hàm mật độ bằng giải thưởng Nobel Hóa học năm 1998. Tiếp theo đây chúng tôi sẽ trình bày cụ thể hơn về lý thuyết phiếm hàm mật độ. Bài toán của hệ nhiều hạt Trạng thái của hệ bao gồm N điện tử và M hạt nhân về nguyên lý có thể thu được từ việc giải phương trình Schrödinger không phụ thuộc thời gian cho hệ nhiều hạt:  N   2   1 N e 2      2     i  V ( r ext i )      (r1 ,.1)  i 1  2m  2 i  j 1 ri  r j  trong đó áp dụng giả thiết gần đúng Borh-Openheimer [1].

là vị trí của điện tử thứ i, Vext là trường ngoài nơi mà các điện tử dịch chuyển, và E là năng lượng điện tử tổng cộng.Thông thường, Vext là thế tĩnh điện được tạo ra bởi các hạt nhân, tuy nhiên, Vext cũng có thể là tác động của môi trường xung quanh hoặc những nhiễu loạn khác trong hệ. Giải phương trình (1.1) cho mỗi một bộ tập hợp các tọa độ hạt nhân khác nhau sẽ thu được năng lượng điện tử của hệ như là một hàm của cấu trúc:   E  E ( R1 ,.2) thêm vào năng lượng tương tác hạt nhân-hạt nhân (Enn), chúng ta có được tổng năng lượng: Etot = E + Enn (1.3) Mặc dù trong phương trình (1.1), chúng tôi đã bỏ qua tọa độ spin để đơn giản hóa vấn đề, nó vẫn không thể giải phương trình (1.1) cho trường hợp chung tổng quát do hàm riêng  phụ thuộc vào 3N vị trí tọa độ. Trong những năm 1930 Hartree và Fock đã đề xuất phương pháp số đầu tiên để giải phương trình này và thu được một hàm sóng gần đúng và tổng năng lượng điện tử [11,15]. Kể từ khi ra đời 4 z phương pháp Hartree Fock (HF), các kỹ thuật dựa trên hàm sóng đã trải qua một quá trình phát triển mạnh mẽ [27,28].

Có nhiều phương pháp tiếp cận tiên tiến để giải quyết vấn đề về hệ nhiều hạt dựa trên các hàm song.Ví dụ như phương pháp cấu hình tương tác (CI) [29], phương pháp liên kết đám (CC) [29], và các phương pháp trường tự hợp đa cấu hình (MCSCF và CASSCF) [26]. Bên cạnh việc phát triển các phương pháp tính toán số dựa trên hàm sóng, lý thuyết phiếm hàm mật độ là một công cụ đắc lực khác để giải bài toán hệ nhiều hạt. Trong lý thuyết DFT, năng lượng điện tử tổng cộng được biểu diễn như là một phiếm hàm của mật độ điện tử (E[ρ(r)]) thay vì hàm sóng. Cách tiếp cận này đã chuyển bài toán hệ nhiều hạt thành bài toán gần đúng một điện tử và do vậy cho phép giải các bài toán hệ nhiều hạt với độ chính xác rất cao.

Cho đến ngày nay, DFT đã trở thành một phương pháp cơ học lượng tử phổ biến và thành công để giải quyết bài toán hệ nhiều hạt [21,23,28]. Làm thế nào để xác định được chính xác phiếm hàm năng lượng điện tử tổng cộng thông qua mật độ điện tích là mục đích của DFT. Do đó, người ta có thể nói rằng lịch sử của DFT là sự phát triển của phiếm hàm năng lượng điện tử tổng cộng E[ρ(r)]. Đó là lý do tại sao tôi lại muốn trình bày DFT như là sự tiến hóa của E[ρ(r)].

Ý tưởng ban đầu về DFT: Thomas-Fermi và các mô hình liên quan Lịch sử của DFT bắt đầu với các nghiên cứu của Thomas và Fermi trong những năm 1927 [7,8,9,30]. Các tác giả này đã nhận ra rằng việc xem xét trên quan điểm thống kê có thể được sử dụng để ước tính sự phân bố của điện tử trong một nguyên tử. Các giả định của Thomas là rằng: “Các điện tử được phân bố đồng nhất trong không gian pha sáu chiều đối với chuyển động của một điện tử với hệ số 2 cho mỗi thể tích h3” và có một trường thế hiệu dụng được xác định bởi điện tích hạt nhân và sự phân bố của các điện tử. Sự biểu diễn năng lượng điện tử tổng cộng thông qua mật độ điện tích có thể được bắt nguồn từ những giả thuyết này.

Ở đây tôi sẽ dẫn dắt một cách hơi khác, nhưng tương đương với cách dẫn ra công thức Thomas-Fermi. Bắt đầu từ phương trình Schrödinger cho một nguyên tử kiểu hydro.4)  2m r  Giá trị năng lượng kỳ vọng là:  2 e2   E   * (r) 2  Z (r)dr  2m r *  2 2   *  e2      (r)  (r)dr    (r) Z (r)dr  2m   r *  2 2   e 2     (r)  (r)dr  Ze (r) dr  2m  r  2   (r)    * (r) 2 (r)dr  Ze dr (1.5) 2m  r  electron   nucleus attraction energy kinetic energy Phương trình (1.5) chỉ ra rằng năng lượng của lực đẩy điện tử-hạt nhân của điện tử có thể được biểu diễn thông qua mật độ điện tử ρ(r). Khó khăn nhất là làm thế nào để biểu diễn động năng của điện tử thông qua ρ(r). Vấn đề này được giải quyết thông qua mô hình của một chất khí điện tử đồng nhất.

Trong mô hình này, không gian được chia thành nhiều khối nhỏ (tế bào), với độ dài l và thể tích ΔV = l3, chứa một số điện tử cố định ΔN, và các điện tử trong mỗi một tế bào biểu hiện như các fermion độc lập ở 0 K, với giả thiết các tế bào độc lập với nhau. Khi đó, năng lượng của điện tử chính xác bằng động năng với các mức năng lượng của nó được cho bởi công thức: h2  (n x , n y , n z )  2 (n x2  n 2y  n z2 ) 8ml h2  2 R2 (1.6) 8ml trong đó nx, ny, nz = 1, 2, 3,. Đối với các số lượng tử cao hay là R lớn, số lượng các mức năng lượng riêng biệt với năng lượng nhỏ hơn ε có thể được tính xấp xỉ bằng 1/8 thể tích của hình cầu với bán kính R trong không gian (nx, ny, nz). Con số này là: 6 z 3/ 2 1  4R 3    8ml 2  ( )      (1.7) 8  3  6  h 2  Số lượng các mức năng lượng giữa ε và ε + δε là: g ( )   (   )   ( ) 3/ 2   8ml 2      1 / 2  O (( ) 2 ) (1.8) 4  h 2  trong đó g(ε) là mật độ trạng thái tại năng lượng ε.

Để tính toán tổng năng lượng (động năng) cho các tế bào với ΔN điện tử, chúng ta cần biết xác suất trạng thái có năng lượng ε bị chiếm giữ, ký hiệu là f(ε). Vì đây là hệ hạt Fermion nên tuân theo phân bố Fermi-Dirac: 1 f ( )  (1.10 )  0,    F trong đó εF là năng lượng Fermi. Tất cả các trạng thái có năng lượng nhỏ hơn εF đều bị chiếm và những trạng thái có mức năng lượng lớn hơn εF không bị chiếm. Năng lượng Fermi εF là giới hạn tại nhiệt độ không của thế hóa μ.

Bây giờ chúng tôi đi tìm năng lượng tổng cộng của các điện tử trong tế bào này bằng cách tổng hợp các đóng góp từ các trạng thái năng lượng khác nhau: E  2   f ( ) g ( )d 3/ 2 F  2m   4  2  l   3 / 2 d 3 h  0 3/ 2 8  2m     l 3 F5 / 2 (1.11) 5  h2  7 z trong đó hệ số 2 được cho vào là do mỗi mức năng lượng bị chiếm bởi hai điện tử, một điện tử với spin α và một điện tử khác với spin β.