Nghiên Cứu Trạng Thái Tồn Tại và Cấu Trúc Dạng Dime của Các Tổ Hợp Phân Tử Ankylliti_Liti Metylat

Khóa luận nghiên cứu cấu trúc và trạng thái tồn tại của tổ hợp phân tử ankylliti liti metylat bằng phần mềm Gaussian 98w trong dung môi THF.

• Tác giả: Hồ Chí Thành
• Người hướng dẫn: Professor. Nguyễn Văn Ngân
• Trường học: Trường Đại Học Sư Phạm TP. HCM
• Chuyên ngành: Hóa Lượng Tử
• Đề tài: Nghiên Cứu Trạng Thái Tồn Tại và Cấu Trúc Dạng Dime của Các Tổ Hợp Phân Tử Ankylliti_Liti Metylat
• Loại tài liệu: khóa luận tốt nghiệp
• Năm xuất bản: 2006
• Địa điểm: TP. Hồ Chí Minh

Mục lục chi tiết

• LỜI MỞ ĐẦU
• 1. CHƯƠNG 1: SƠ LƯỢC VỀ CHƯƠNG TRÌNH GAUSSIAN 98W
• 1.1. CƠ SỞ LÝ THUYẾT CỦA CHƯƠNG TRÌNH GAUSSIAN 98W
• 1.1.1. Lý thuyết trường tự hợp Hartree-Fock (HF)
• 1.1.2. Thuyết hàm mật độ (Density-functional theory DFT)
• 1.2. CÁC BỘ HÀM CƠ BẢN
• 1.2.1. Các tập cơ sở tối thiểu (Minimal Basic Sets)
• 1.2.2. Tập cơ sở phân chia hóa trị (Split Valence Basic Sets)
• 1.2.3. Tập cơ sở phân cực (Polarized Basic Sets)
• 1.2.4. Tập cơ sở chứa hàm khuyếch tán (Diffuse Function Basic Sets)
• 1.2.5. Tập cơ sở động lượng góc cao (High Angular Momentum Basic Sets)
• 1.3. ĐỘ BỀN VÀ CẤU TRÚC CỦA CÁC BANION
• 1.4. CẤU TRÚC CỦA CÁC HỢP CHẤT CƠ KIM
• 2. CHƯƠNG 2: KẾT QUẢ VÀ THẢO LUẬN
• 2.1. CẤU TRÚC TỐI ƯU HÓA CỦA CÁC TỔ HỢP PHÂN TỬ ANKYLLITI_LITI
• 2.2. NĂNG LƯỢNG TỰ DO GIBBS CỦA CÁC TỔ HỢP PHÂN TỬ
• 2.3. CẤU TRÚC TỔ HỢP PHÂN TỬ METYLLITI_LITI METYLAT
• KẾT LUẬN - ĐỀ XUẤT
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Về Nghiên Cứu Trạng Thái Tồn Tại Của Ankylliti_Liti Metylat

Nghiên cứu trạng thái tồn tại và cấu trúc dạng dime của các tổ hợp phân tử Ankylliti_Liti Metylat là một lĩnh vực quan trọng trong hóa học hữu cơ hiện đại. Các hợp chất này có cấu trúc phức tạp và đóng vai trò quan trọng trong nhiều ứng dụng công nghiệp. Việc hiểu rõ về trạng thái tồn tại của chúng giúp tối ưu hóa quy trình tổng hợp và ứng dụng trong thực tiễn.

1.1. Khái Niệm Về Trạng Thái Tồn Tại Trong Hóa Học

Trạng thái tồn tại của các hợp chất hóa học đề cập đến hình thức mà chúng tồn tại trong điều kiện nhất định. Điều này bao gồm các yếu tố như nhiệt độ, áp suất và môi trường xung quanh.

1.2. Cấu Trúc Dạng Dime Của Các Tổ Hợp Phân Tử

Cấu trúc dạng dime của các tổ hợp phân tử Ankylliti_Liti Metylat có thể ảnh hưởng đến tính chất vật lý và hóa học của chúng. Việc nghiên cứu cấu trúc này giúp hiểu rõ hơn về cách mà các phân tử tương tác với nhau.

II. Vấn Đề và Thách Thức Trong Nghiên Cứu Cấu Trúc Phân Tử

Nghiên cứu cấu trúc phân tử Ankylliti_Liti Metylat gặp nhiều thách thức, bao gồm việc xác định chính xác cấu trúc và trạng thái tồn tại của chúng. Các yếu tố như điều kiện thí nghiệm và phương pháp phân tích có thể ảnh hưởng đến kết quả nghiên cứu.

2.1. Những Khó Khăn Trong Phân Tích Cấu Trúc

Việc phân tích cấu trúc của các tổ hợp phân tử thường yêu cầu các thiết bị hiện đại và kỹ thuật phức tạp. Điều này có thể gây khó khăn cho các nhà nghiên cứu, đặc biệt là trong các nước đang phát triển.

2.2. Ảnh Hưởng Của Điều Kiện Thí Nghiệm

Điều kiện thí nghiệm như nhiệt độ và áp suất có thể ảnh hưởng lớn đến trạng thái tồn tại của các hợp chất. Việc kiểm soát các yếu tố này là rất quan trọng để có được kết quả chính xác.

III. Phương Pháp Nghiên Cứu Trạng Thái Tồn Tại Của Ankylliti_Liti Metylat

Để nghiên cứu trạng thái tồn tại và cấu trúc của các tổ hợp phân tử Ankylliti_Liti Metylat, nhiều phương pháp khác nhau đã được áp dụng. Trong đó, phương pháp tính toán lượng tử và phần mềm Gaussian là những công cụ hữu ích.

3.1. Phương Pháp Tính Toán Lượng Tử

Phương pháp tính toán lượng tử cho phép mô phỏng và dự đoán cấu trúc của các phân tử một cách chính xác. Điều này giúp tiết kiệm thời gian và chi phí trong nghiên cứu.

3.2. Sử Dụng Phần Mềm Gaussian Trong Nghiên Cứu

Phần mềm Gaussian là một công cụ mạnh mẽ trong việc tính toán hóa học. Nó giúp các nhà nghiên cứu mô phỏng các phản ứng hóa học và phân tích cấu trúc phân tử một cách hiệu quả.

IV. Ứng Dụng Thực Tiễn Của Nghiên Cứu Ankylliti_Liti Metylat

Nghiên cứu về Ankylliti_Liti Metylat không chỉ có giá trị lý thuyết mà còn có nhiều ứng dụng thực tiễn trong công nghiệp. Các hợp chất này có thể được sử dụng trong sản xuất vật liệu mới và trong các quy trình hóa học khác.

4.1. Ứng Dụng Trong Ngành Hóa Học Hữu Cơ

Các tổ hợp phân tử Ankylliti_Liti Metylat có thể được sử dụng trong tổng hợp hữu cơ, giúp tạo ra các sản phẩm mới với tính chất ưu việt.

4.2. Tiềm Năng Trong Công Nghiệp Vật Liệu

Nghiên cứu về cấu trúc và tính chất của các hợp chất này có thể dẫn đến việc phát triển các vật liệu mới với ứng dụng trong nhiều lĩnh vực khác nhau.

V. Kết Luận và Tương Lai Của Nghiên Cứu Ankylliti_Liti Metylat

Nghiên cứu trạng thái tồn tại và cấu trúc của các tổ hợp phân tử Ankylliti_Liti Metylat mở ra nhiều hướng đi mới trong hóa học. Tương lai của lĩnh vực này hứa hẹn sẽ có nhiều phát triển đáng kể.

5.1. Tương Lai Của Nghiên Cứu Hóa Học

Với sự phát triển của công nghệ và phương pháp nghiên cứu mới, lĩnh vực hóa học sẽ tiếp tục phát triển mạnh mẽ, mở ra nhiều cơ hội mới cho các nhà nghiên cứu.

5.2. Định Hướng Nghiên Cứu Tiếp Theo

Các nghiên cứu tiếp theo có thể tập trung vào việc tối ưu hóa quy trình tổng hợp và ứng dụng của các hợp chất Ankylliti_Liti Metylat trong thực tiễn.

Nội dung chính

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

Trong hóa học hữu cơ hiện đại, hóa học cơ kim đóng vai trò then chốt trong các phản ứng tổng hợp cấu trúc phức tạp. Bên cạnh các hợp chất Grignard truyền thống, hóa học cơ liti (organolithium) ngày càng khẳng định vị trí đặc biệt với khả năng phản ứng mạnh mẽ và tính ứng dụng đa dạng. Tuy nhiên, cấu trúc thực tế của các hợp chất cơ kim chứa liti vô cùng phức tạp; chúng không tồn tại dưới dạng các phân tử đơn lẻ (monome) mà luôn có xu hướng liên hợp tạo thành các tổ hợp phân tử (aggregates) như dime, tetrame hoặc hexame tùy thuộc vào bản chất nhóm thế, nồng độ, nhiệt độ và dung môi. Đặc biệt, sự pha trộn giữa các phân tử ankylliti với các hợp chất liti ankylat (như liti metylat) hình thành nên các tổ hợp dị thể (mixed aggregates) mang những đặc tính cấu trúc và hoạt tính hóa học mới lạ.

Phương pháp nghiên cứu thực nghiệm cung cấp các dữ liệu có độ chính xác cao nhưng đòi hỏi nhiều thời gian, công sức và trang thiết bị chuyên dụng đắt tiền (như phổ cộng hưởng từ hạt nhân nhiệt độ thấp, tinh thể học tia X), vốn là rào cản lớn đối với điều kiện nghiên cứu tại các nước đang phát triển. Trong bối cảnh đó, phương pháp hóa học tính toán lượng tử (quantum chemistry) với sự hỗ trợ của phần mềm Gaussian mở ra khả năng mô phỏng, giải thích cơ chế ở cấp độ nguyên tử/phân tử và đưa ra những dự đoán lý thuyết tin cậy trước khi tiến hành thực nghiệm.

Khóa luận được xây dựng nhằm giải quyết các mục tiêu và nhiệm vụ nghiên cứu sau:

1. Tổng quan hệ thống cơ sở lý thuyết về hóa học lượng tử, nguyên lý hoạt động của phần mềm tính toán Gaussian 98W, bản chất liên kết cơ kim và độ bền của cacbanion.
2. Ứng dụng phần mềm Gaussian 98W để mô phỏng và tối ưu hóa cấu trúc hình học ở trạng thái dừng của các tổ hợp phân tử ankylliti và liti metylat ở dạng dime.
3. Tính toán và đánh giá các thông số nhiệt động lực học, trọng tâm là năng lượng tự do Gibbs ($\Delta G$) của các dạng tổ hợp ở thể khí và trong môi trường dung môi tetrahydrofuran (THF).
4. Phân tích chi tiết đặc điểm liên kết, độ dài liên kết và sự biến đổi cấu trúc của tổ hợp phân tử metylliti - liti metylat ($CH_3Li \cdot LiOCH_3$) dưới ảnh hưởng của dung môi.

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

• Đối tượng nghiên cứu: Các phân tử ankylliti đơn giản, liti metylat ($LiOCH_3$) và các tổ hợp phân tử dạng dime đồng thể ($(RLi)_2$, $(LiOCH_3)_2$) cùng tổ hợp dị thể ($RLi \cdot LiOCH_3$).
• Phạm vi nghiên cứu: Khảo sát cấu trúc và trạng thái tồn tại ở pha khí và trong môi trường dung môi tetrahydrofuran (THF) bằng các mô hình tính toán hóa lượng tử trên phần mềm Gaussian 98W, thực hiện tại Trường Đại học Sư phạm TP. Hồ Chí Minh năm 2006.

---

Cơ sở lý thuyết và phương pháp

Khóa luận kết hợp chặt chẽ giữa lý thuyết hóa lượng tử tính toán và các nguyên lý nhiệt động, cấu trúc hóa học cơ kim:

Khung lý thuyết và mô hình tính toán

• Phép gần đúng Born-Oppenheimer: Tách biệt chuyển động của hạt nhân và electron dựa trên sự chênh lệch lớn về khối lượng, cho phép giải phương trình Schrödinger cho hệ điện tử trên một khung hạt nhân cố định.
• Lý thuyết trường tự hợp Hartree-Fock (HF-SCF): Sử dụng định thức Slater cho hệ spin-orbital để giải phương trình Fock hữu hiệu $F\psi_i = \epsilon_i \psi_i$ qua quy trình lặp tự hợp.
• Thuyết hàm mật độ (Density Functional Theory - DFT): Dựa trên định lý Hohenberg-Kohn (1964) khẳng định năng lượng và mọi tính chất của trạng thái cơ bản là phiếm hàm của mật độ electron $\rho_0(r)$, kết hợp phương pháp Kohn-Sham (1965) sử dụng hệ giả định các electron không tương tác để tính toán chính xác mật độ điện tích.
• Các phiếm hàm trao đổi - tương quan:
  • Xấp xỉ mật độ cục bộ (LDA) và xấp xỉ mật độ spin cục bộ (LSDA/SVWN).
  • Xấp xỉ gradient suy rộng (GGA): B88 (Becke 88), PW86, PW91, LYP, P86, PBE.
  • Phiếm hàm lai hóa (Hybrid Functionals): Trọng tâm là phiếm hàm B3LYP (Becke 3 tham số kết hợp tương quan Lee-Yang-Parr) với công thức trộn trao đổi chính xác Hartree-Fock: $$E_{xc}^{B3LYP} = (1 - a_0) E_x^{LSDA} + a_0 E_x^{HF} + a_x \Delta E_x^{B88} + (1 - a_c) E_c^{LSDA} + a_c E_c^{LYP}$$ (trong đó các hệ số thực nghiệm: $a_0 = 0{,}20$, $a_x = 0{,}72$, $a_c = 0{,}81$).
• Hệ tập cơ sở (Basis Sets): Khảo sát từ các bộ cơ sở phân chia hóa trị (Split-valence) đến các bộ cơ sở phân cực (Polarized basis sets như 6-31G(d), 6-31G(d,p)) và bộ cơ sở chứa hàm khuếch tán (Diffuse functions như 6-31+G(d), 6-31++G(d)) nhằm mô tả chính xác đám mây electron của nguyên tử liti, oxy và cacbanion.

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

• Phương pháp tính toán lượng tử: Tối ưu hóa cấu trúc hình học hoàn toàn (Full geometry optimization) bằng gradient giải tích của năng lượng; tính toán ma trận Hessian (tần số dao động chuẩn) để xác nhận trạng thái cực tiểu năng lượng (không có tần số ảo) và thu nhận năng lượng điểm không (ZPE) cùng năng lượng tự do Gibbs ($G$).
• Mô hình dung môi: Đánh giá ảnh hưởng của dung môi tetrahydrofuran (THF) lên tính bền nhiệt động và sự phân bố điện tích của các phức tổ hợp.
• Nguồn dữ liệu: Dữ liệu tính toán ab initio/DFT thu nhận trực tiếp từ Gaussian 98W; dữ liệu thực nghiệm đối chiếu trích xuất từ các công trình đo $pK_a$, cân bằng Schlenk, nhiễu xạ tia X và phổ NMR của các tác giả Applequist, O'Brien, Dessy, Cram, Rundle, Smith và Becker.

---

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

Phần 1. Tổng quan

Chương 1. Sơ lược về chương trình Gaussian 98W

Chương này trình bày chi tiết nền tảng toán học và vật lý lượng tử được tích hợp trong phần mềm Gaussian 98W:

• Thiết lập phương trình Schrödinger đơn điện tử, toán tử Hamiltonian, ma trận Fock và quy trình lặp trường tự hợp Hartree-Fock-Roothaan.
• Trình bày toán học định lý Hohenberg-Kohn: chứng minh tính duy nhất của thế ngoại năng $v(r)$ thông qua mật độ xác suất electron $\rho_0(r)$, định lý biến thiên cho phiếm hàm năng lượng $E[\rho]$.
• Xây dựng hệ phương trình Kohn-Sham, phân tích thành phần năng lượng trao đổi - tương quan $E_{xc}[\rho]$, cơ sở của phép xấp xỉ LDA, LSDA và các hàm hiệu chỉnh gradient (GGA, BLYP, B3PW91, B1B95, EDF1, K2-BVWN).
• Phân loại và cấu trúc các bộ hàm cơ bản dạng Gauss (GTOs):
  • Tập cơ sở phân cực: Thêm hàm orbital $d$ cho nguyên tử nặng và hàm $p$ cho hydro (như 6-31G* / 6-31G(d)).
  • Tập cơ sở khuếch tán: Thêm hàm $s, p$ kích thước lớn (như 6-31+G(d)) mô tả anion và cặp electron tự do.
  • Tập cơ sở động lượng góc cao: Cấu trúc triple-zeta bổ sung nhiều hàm đa cực (như 6-311+G(2d,p), 6-311++G(3df,2pd)).

Chương 2. Độ bền và cấu trúc của cacbanion

Tổng hợp lý thuyết về bản chất bazơ và độ bền nhiệt động, động học của các tâm cacbanion:

• Mối liên hệ giữa độ bền cacbanion và tính axit ($pK_a$) của axit liên hợp tương ứng (axit càng yếu, bazơ liên hợp càng mạnh, cacbanion càng kém bền đối với tác nhân cho proton).

Axit liên hợp	Bazơ (Cacbanion)	$pK_a$ tương đối (trong nước)
$CH(CN)_3$	$^-C(CN)_3$	$\approx 0$
$RCOCH_2CHO$	$RCOCHCHO^-$	5
$CH_3COCH_2COCH_3$	$CH_3COCHCOCH_3^-$	9
$RCH_2NO_2$	$RCHNO_2^-$	10
$CH_3COCH_2COOR$	$CH_3COCHCOOR^-$	11
$EtOOCCH_2COOEt$	$EtOOCC^-HCOOEt$	13
$RCH_2CHO$	$RCHCHO^-$	16
$RCOCH_2R$	$RCOCHR^-$	19–20
$PhCH_2SO_2CH_2Ph$	$PhCHSO_2CH_2Ph^-$	18,91
$ROOCCH_2R$	$ROOCC^-HR$	24,5
$HC \equiv CH$	$HC \equiv C^-$	25
$Ar_3CH$	$Ar_3C^-$	31,5
$PhCH_3$	$PhCH_2^-$	40
$CH_2=CHCH_3$	$CH_2\cdots CH\cdots CH_2^-$	43
$H_2C=CH_2$	$H_2C=CH^-$	44
$cyclo-C_3H_6$	$cyclo-C_3H_5^-$	46
$CH_4$	$CH_3^-$	48
$C_2H_6$	$C_2H_5^-$	50
$(CH_3)_2CH_2$	$(CH_3)_2CH^-$	51

• Thứ tự độ bền của cacbanion đơn giản:
  • Trong dung dịch (nhiệt động và động học trao đổi hydro): $metyl > bậc\ 1 > bậc\ 2 > bậc\ 3$.
  • Ở pha khí (thực nghiệm phân cắt Si-C): $neopentyl > cyclopropyl > tert-butyl > n-propyl > metyl > isopropyl > etyl$.
• Các yếu tố làm bền cấu trúc cacbanion:
  1. Hiệu ứng liên hợp $\pi$: Sự xen phủ với liên kết chưa no kề cận (hệ alyl, benzyl, diphenylmetyl, triphenylmetyl, anion enolat). Anion thơm cyclopentadienyl ($pK_a = 16$) có năng lượng cộng hưởng đạt 24–27 kcal/mol.
  2. Đặc tính orbital $s$: Độ bền tăng theo độ âm điện hiệu dụng: $sp\ (50%\ s) > sp^2\ (33%\ s) > sp^3\ (25%\ s)$; cacbanion cyclopropyl bền hơn do orbital ngoài vòng có đặc tính $s$ cao ($\approx 33%$).
  3. Nguyên tử dị tố $\alpha$ (S, P, Si): Tương tác ổn định hóa điện tích âm qua hiệu ứng phân cực và liên hợp.
  4. Hiệu ứng trường và cảm ứng: Thứ tự nhóm chức ở vị trí $\alpha$ làm bền cacbanion giảm dần: $$NO_2 > RCO > COOR > SO_2 > CN \approx CONH_2 > Hal > H > R$$
• Cấu trúc lập thể: Dạng hình tháp lai hóa $sp^3$ (nghịch chuyển nhanh giống amin) đối với cacbanion no đơn giản; dạng phẳng lai hóa $sp^2$ đối với hệ liên hợp hoặc enolat.

Chương 3. Cấu trúc của các hợp chất cơ kim

• Phân tích bản chất liên kết kim loại - cacbon (C-M): Tùy thuộc vào hiệu độ âm điện, phần trăm đặc tính ion của các liên kết được tổng hợp:

Liên kết C—M	Hiệu độ âm điện	Phần trăm đặc tính ion (%)
C—K	2,5 - 0,8 = 1,7	69
C—Na	2,5 - 0,9 = 1,6	64
C—Li	2,5 - 1,0 = 1,5	60
C—Mg	2,5 - 1,2 = 1,3	52
C—Al	2,5 - 1,5 = 1,0	40
C—Zn	2,5 - 1,6 = 0,9	36
C—Sn	2,5 - 1,8 = 0,7	28
C—Cu	2,5 - 1,9 = 0,6	24
C—Hg	2,5 - 1,9 = 0,6	24

• Trạng thái tổ hợp của ankylliti:
  • Trong dung môi hydrocacbon (benzen, cyclohexan): tồn tại chủ yếu ở dạng hexame hoặc tetrame.
  • Trong dung môi ete: tồn tại ở dạng tổ hợp 2 đến 5 phân tử (tetrame, dime).
  • Trong dung môi phối trí mạnh THF: $tert-butylliti$ chuyển hẳn sang dạng monome; $neopentylliti$ là hỗn hợp cân bằng monome - dime. Metylliti ở trạng thái rắn và trong ete giữ cấu trúc cụm tứ diện.
• Thuốc thử Grignard và cân bằng Schlenk: $$2RMgX \rightleftharpoons R_2Mg + MgX_2 \rightleftharpoons R_2Mg \cdot MgX_2$$ Nhiễu xạ tia X và phổ NMR chứng minh trong ete, magie phối trí tạo phức monome $RMgX \cdot 2Et_2O$; trong dung môi THF, tác nhân Grignard tồn tại chủ yếu ở dạng monome, độ dịch chuyển cân bằng phụ thuộc bản chất nhóm ankyl/aryl, halogen và nhiệt độ.

---

Phần 2. Kết quả nghiên cứu

Phương pháp tính toán và Chương 2. Kết quả và thảo luận

Chương này tập trung trình bày các dữ liệu lượng tử thu nhận từ việc tối ưu hóa và tính toán nhiệt động:

• Cấu trúc tối ưu hóa của các tổ hợp phân tử ankylliti - liti metylat:
  • Mô phỏng các dạng dime đồng thể $(RLi)_2$, $(LiOCH_3)_2$ và dime dị thể $RLi \cdot LiOCH_3$.
  • Nhân trung tâm của phức dime dị thể tạo thành vòng phẳng 4 cạnh gồm hai nguyên tử Liti đóng vai trò cầu nối xen kẽ giữa một nguyên tử Cacbon (từ gốc ankyl) và một nguyên tử Oxy (từ nhóm metylat): khung dạng $[Li_2CO]$.
• Năng lượng tự do Gibbs của các tổ hợp phân tử:
  • So sánh độ bền nhiệt động ($\Delta G$) giữa các dạng dime đồng thể và dime dị thể ở thể khí và trong THF.
  • Sự tạo thành dime dị thể $RLi \cdot LiOCH_3$ từ các dime đồng thể $(RLi)_2$ và $(LiOCH_3)_2$ là quá trình thuận lợi về mặt nhiệt động lực học trong dung môi THF nhờ sự bù trừ điện tích hiệu quả giữa nguyên tử oxy giàu điện tử và nguyên tử liti.
• Cấu trúc chi tiết của tổ hợp phân tử metylliti - liti metylat ($CH_3Li \cdot LiOCH_3$):
  • Khảo sát các thông số hình học: độ dài liên kết cầu nối $Li-O$, $Li-C$, khoảng cách liên kim loại $Li\cdots Li$, các góc liên kết $\widehat{Li-O-Li}$ và $\widehat{Li-C-Li}$.
  • Đánh giá sự biến thiên mật độ điện tích: liên kết $Li-O$ mang đặc tính ion cao hơn liên kết $Li-C$, làm phân cực mạnh khung trung tâm và tạo điều kiện cho các phân tử dung môi THF tạo liên kết phối trí solvat hóa bền vững vào tâm $Li^+$.

---

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

Kết quả chính

• Xác định được bản chất liên kết và hình học không gian của tổ hợp dime dị thể ankylliti - liti metylat ($RLi \cdot LiOCH_3$): cấu trúc vòng phẳng 4 cạnh liti - cacbon - liti - oxy với các liên kết cầu nối đa tâm.
• Làm rõ định lượng mức độ phân cực liên kết cơ kim: liên kết C-Li có $60%$ đặc tính ion (hiệu độ âm điện $1{,}5$), nằm giữa liên kết C-Na ($64%$) và C-Mg ($52%$).
• Chứng minh bằng tính toán lượng tử rằng việc tạo phức dị thể giữa ankylliti và liti ankylat trong dung môi THF làm thay đổi đáng kể năng lượng tự do Gibbs và trạng thái cân bằng tổ hợp so với pha khí.

Đề xuất và kiến nghị của tác giả

• Ứng dụng mô hình hóa học lượng tử DFT/B3LYP kết hợp bộ cơ sở phân cực và khuếch tán để dự đoán trạng thái tổ hợp của các tác nhân cơ liti trước khi thiết kế phản ứng tổng hợp hữu cơ.
• Chú trọng kiểm soát sự có mặt của các hợp chất liti ankylat (vốn thường sinh ra do tác nhân cơ liti bị oxy hóa một phần hoặc tạo thành trong quá trình phản ứng) vì chúng hình thành các phức dị thể làm thay đổi hoàn toàn hoạt tính của tác nhân ankylliti ban đầu.

Đóng góp khoa học

• Đặt nền móng cho hướng nghiên cứu hóa học tính toán chuyên sâu về hợp chất cơ kim chứa liti tại Khoa Hóa học - Trường Đại học Sư phạm TP. Hồ Chí Minh với sự phối hợp học thuật quốc tế (Đại học Fisk, Hoa Kỳ).
• Cung cấp bức tranh toàn diện, có hệ thống từ lý thuyết trường tự hợp HF, DFT đến các dữ liệu định lượng cụ thể về cấu trúc và nhiệt động của các dime dị thể ankylliti - liti metylat.

---

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

• Hạn chế:
  • Khóa luận giới hạn phạm vi khảo sát ở dạng tổ hợp phân tử dime ($2:0, 1:1, 0:2$), chưa mở rộng tính toán cho các cấu trúc tổ hợp bậc cao hơn như tetrame dạng khối lập phương lệch (tetrameric aggregates) hay hexame do giới hạn năng lực máy tính tại thời điểm năm 2006.
  • Ảnh hưởng của dung môi THF chủ yếu được đánh giá ở cấp độ mô hình hóa liên kết phối trí cơ bản, chưa xét đến hiệu ứng solvat hóa động với nhiều lớp vỏ dung môi rõ ràng.
• Hướng mở rộng:
  • Mở rộng tính toán cấu trúc tổ hợp tetrame dị thể $(RLi)n(LiOR'){4-n}$.
  • Khảo sát bề mặt năng lượng tiềm năng (PES) và trạng thái chuyển tiếp (transition states) của các phản ứng cộng hợp ái nhân sử dụng phức dị thể ankylliti - liti ankylat vào các hợp chất cacbonyl.

---

Giá trị tham khảo

Tài liệu là nguồn tham khảo chuyên môn giá trị cho:

• Sinh viên và học viên cao học ngành Hóa học (Hóa Lượng tử, Hóa lý thuyết, Hóa Hữu cơ - Cơ kim): Nắm vững quy trình mô phỏng lượng tử trên phần mềm Gaussian, cách lựa chọn phiếm hàm DFT (B3LYP, GGA) và hệ tập cơ sở phù hợp cho hệ cơ kim.
• Người nghiên cứu tổng hợp hữu cơ: Tham khảo bảng số liệu $pK_a$ chuẩn của các hệ cacbanion, phần trăm đặc tính ion của các liên kết C-M và cơ chế cân bằng tổ hợp phân tử trong các dung môi hydrocacbon, ete và THF.
• Phần đáng chú ý nhất: Tổng quan toán học chi tiết về thuyết DFT (Hohenberg-Kohn và Kohn-Sham) ở Chương 1 và phân tích cấu trúc phức dị thể trung tâm $[Li_2CO]$ ở Phần 2.

---

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

1. Tại sao các hợp chất ankylliti không tồn tại ở dạng monome trong dung dịch mà tạo thành các tổ hợp phân tử?

Liên kết C-Li có $60%$ đặc tính ion và mang tính thiếu hụt điện tử (electron-deficient). Để đạt trạng thái bền vững về mặt năng lượng, các nguyên tử liti sử dụng orbital trống để tạo các liên kết đa tâm (liên kết 4 tâm - 2 electron hoặc 3 tâm - 2 electron) với các nguyên tử cacbon hoặc oxy kề cận, dẫn đến sự hình thành các cụm dime, tetrame hoặc hexame tùy dung môi.

2. Sự khác biệt về thứ tự độ bền cacbanion trong dung dịch và ở pha khí là gì?

• Trong dung dịch: Độ bền giảm theo thứ tự $metyl > bậc\ 1 > bậc\ 2 > bậc\ 3$ do hiệu ứng đẩy electron của các nhóm ankyl làm tăng mật độ điện tích âm trên nguyên tử cacbon trung tâm.
• Ở pha khí: Hiệu ứng phân cực hóa chiếm ưu thế dẫn đến thứ tự $neopentyl > cyclopropyl > tert-butyl > n-propyl > metyl > isopropyl > etyl$.

3. Phiếm hàm lai hóa B3LYP được xây dựng như thế nào trong chương trình Gaussian?

B3LYP kết hợp trao đổi cục bộ LSDA, trao đổi chính xác Hartree-Fock ($20%$), trao đổi hiệu chỉnh gradient Becke 88 ($72%$), tương quan cục bộ LSDA ($19%$) và tương quan hiệu chỉnh gradient LYP ($81%$) thông qua 3 tham số tối ưu hóa thực nghiệm ($a_0 = 0{,}20$, $a_x = 0{,}72$, $a_c = 0{,}81$).

4. Cân bằng Schlenk trong dung dịch thuốc thử Grignard là gì và phụ thuộc vào những yếu tố nào?

Cân bằng Schlenk mô tả chuyển hóa thuận nghịch giữa các dạng tồn tại của hợp chất cơ magie: $$2RMgX \rightleftharpoons R_2Mg + MgX_2 \rightleftharpoons R_2Mg \cdot MgX_2$$ Cân bằng phụ thuộc vào bản chất nhóm ankyl/aryl ($R$), halogen ($X$), độ phân cực và khả năng phối trí của dung môi (ete, THF, dioxan), nồng độ và nhiệt độ.

5. Cấu trúc khung trung tâm của phức dime dị thể metylliti - liti metylat có đặc điểm gì?

Tổ hợp phân tử dị thể có nhân trung tâm dạng vòng 4 cạnh phẳng $[Li_2CO]$, trong đó 2 nguyên tử Liti đóng vai trò cầu nối liên kết đồng thời với 1 nguyên tử Cacbon (từ nhóm metyl) và 1 nguyên tử Oxy (từ nhóm metylat), tạo nên sự phân bố điện tích bất đối xứng và tiếp tục nhận phối trí từ dung môi THF.

---

Kết luận

Khóa luận tốt nghiệp của tác giả Hồ Chí Thành đã ứng dụng thành công phương pháp hóa học lượng tử DFT/B3LYP trên phần mềm Gaussian 98W để giải quyết bài toán cấu trúc và nhiệt động học của tổ hợp ankylliti - liti metylat ở thể khí và trong dung môi THF. Công trình đã chứng minh sự hình thành phức dime dị thể $[Li_2CO]$ là quá trình thuận lợi về mặt năng lượng, đồng thời cung cấp một hệ thống lý thuyết chuẩn xác về bản chất liên kết cơ kim và độ bền cacbanion. Đây là tài liệu học thuật có giá trị tham khảo cao cho nghiên cứu hóa học tính toán và hóa học cơ kim thực nghiệm.

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

CHƯƠNG 1. SƠ LƯỢC VE CHƯƠNG TRÌNH GAUSSIAN 98W 1. CƠ SỞ LÝ THUYET CUA CHUONG TRÌNH GAUSSIAN 98W 1. Lý thuyết trường tự hợp Hartree-Fock (HF) 1.

Phép gần đáng Born-Oppenheimer Vì hạt nhân có khối lượng rất lớn (chuyển động rất chậm) so với các electron có khối lượng rất nhỏ (chuyển động rất nhanh) nên các hạt nhân được xem như là cố định. Theo sự đơn giản hóa của Born-Oppenheimer thì hàm sóng Schrédinger cho một hệ thống thuần điện tử có thể viết thành: H,w=E,w (1) Trong đó : H, : toán tử Hamiltonian đơn điện tử. Y, : hàm sóng đơn điện tử. E, : năng lượng của hệ thống thuần điện tử.

Toán tử Hamiltonian cho nguyên tử có n electron được viết thành: Helectron = - 2 2Ze ane 2 TY” M = r? * hàn m4 A B Cc Trong đó : A: toán tử động năng cho n electron. B: thế năng tương tic giữa electron và hạt nhân với điện tích Z. C: thế năng tương tác đẩy giữa các electron. Bằng phương pháp nhiễu loạn ta thu được hàm sóng bậc 0 khi bỏ qua tương tác đẩy giữa các electron.

Phương trình sóng Schrédinger lúc đó được chia làm n phương trình tương tự phương trình viết cho hệ tương tự hydro. Hàm sóng bậc 0 là tích của n orbital giống hyđro : = Fr OOS OOF tr„8„6,) 3) Trong đó những orbital giống hydro có dang : £=R.#) (4) SVTH: Hô Chí Thanh _ lớp Hóa 4B (khóa 2002-2006) 4 Khóa luận tốt nghiệp GVHD: Professor Lawrence M.s Nguyén Vin Ngân Với nguyên tử ở trạng thái cơ bản, chúng ta đặt hai electron với spin ngược chiều vào mỗi orbital có nang lượng thấp theo nguyên lý ngoại trừ Pauli để thu được cấu hình ở trạng thái cơ bản. Bước tiếp theo là sử dụng hàm biến đổi có dạng giống (4) nhưng không giới hạn cho hệ thống giống hyđro hay cho bất cứ một orbital cụ thể nào.0,)8 Mr OrO Zr BP) (5) Để giải phương trình trên chúng ta phải tìm các ham gy, ga.g, làm cực tiểu hóa tích phan biến đổi. Để đơn giản phan nào chúng ta làm gần đúng những AO khả kiến, tốt nhất với orbital là tích của một hàm theo bán kính và một hàm theo hình cầu điều hòa.

Phương pháp trường tự hợp Hartree Thực hiện theo các bước sau: * Một là: Dự đoán một hàm sóng tích: Ó =sr‹8-Ó)s:r»Ø0Ó,)-s,r„Ø.6,) (7) với s, là hàm chuẩn hóa Hàm sóng dự đoán ban đẩu @ có thể chấp nhận được là tích của những orbital giống hyđro với điện tích hạt nhân là điện tích hiệu dụng. Trong phương trình (7) mật độ electron khả kiến của electron thi i là lD gì * Hai là: Xét electron thứ nhất di chuyển trong vùng phân bố điện tử liên tục do electron 2, 3,. Tinh thế năng hiệu dụng V;(r,) và giải phương trình Schrédinger (ứng với hệ thống một electron) cho electron thứ nhất để thu được orbital bậc nhất tạ( 1). * Ba là: Xét tiếp electron thứ hai và xem như nó di chuyển trong đám mây electron với mật độ: › 2 2 2 -eile,(1)| +|sa(3)| +|sa(4)| +--+]s,,(7)] I- Chúng ta cũng tính toán thế năng hiệu dung V(r.) và giải phương trình sóng Schrödinger (ứng với hệ một electron) cho electron thứ hai để thu được orbital bậc một t,(2).

* Bốn là: Tiến hành tương tự ta thu được một bộ các orbital bậc | của n electron. * Năm là: Quay lại electron thứ nhất va lặp lại chu trình để thu được những bộ orbital bậc cao hơn cho đến khi không còn sự thay đổi giữa hai hàm sóng có bậc kế tiếp. * Sáu là: Bộ các orbital cuối cùng là hàm sóng tương ứng thu được bằng phương pháp trường tự hợp của Hatree. SVTH: Hồ Chí Thành _ lớp Hóa 4B (khóa 2002-2006) $ Khóa luận tốt nghiệp GVHD: Professor Lawrence M.

Phương pháp trường tự hợp Hartree-Fock Nhược điểm của phương pháp Hartree là mặc dù đã chú ý đến spin của electron và nguyên lý ngoại trừ Pauli bằng cách không đặt quá hai electưon trong mỗi không gian orbital nhưng bất kì sự ước tính gắn đúng nào đối với hàm sóng thực nên bao gồm sự ngoại trừ spin và là bất đối xứng với các electron trao đổi. Vì vậy thay vì dùng orbital không gian, chúng ta phải dùng những orbital spin và thực hiện tổ hợp tuyến tính không đối xứng là tích của orbital spin. Thực hiện tính toán bằng phương pháp trường tự hợp có sử dụng những orbital spin bất đối được gọi là phương pháp tính toán theo Hartree-Fock. Phương trình để tìm orbital theo Hartree-Fock : Fu-~em, i=1,2,3.

(8) Trong đó : //,: Orbital spin thứ i, F : toán tử Fock (hoặc Hartree-Fock): toán tử Hamilton Hartree-Fock hiệu dụng. £,: năng lượng orbital của orbital spin thứ i. Giải phương trình bằng phương pháp trường tự hợp để thu được những hàm sóng tối ưu. Sơ dé nguyên tắc trường tự hợp SCF UGC TÍNH CÁC HỆ SO VÂN ĐẠO NGUYEN TU ĐÁNH GIÁ NANG LƯỢNG VÀ THÀNH LAP MA TRAN FOCK GIẢI CÁC HAM SÓNG HỘI TY (TRUONG TY HỢP) SAI———— DUNG + CHO KET QUÁ “————r—————————-- SVTH: Ho Chi Thanh _ lớp Hóa 4B (khóa 2002-2006) 6 Khéa luận tốt nghỉ GVHD: Professor Lawrence M.s Nguyễn Văn Ngân Thuyết Hartree-Fock dựa trên phương pháp biến đổi trong cơ hoc lượng tử.

Nếu hàm ® là một hàm chuẩn bất đối bất kì trong toạ độ điện tử thì giá trị năng lượng tương ứng với hàm này sẽ có được từ tích phân : E'= [o" Hdr (9) Trong đó tích phân được lấy trên tọa độ của tất cả các electron. Nếu ® ngẫu nhiên là hàm sóng đúng - cho electron ở trạng thái cơ bản thì nó sẽ thỏa mãn phương trình sóng Schrédinger. Vì ự chuẩn hóa nên E' là nang lượng E đúng. =E[o'odr=E (10) Tuy nhiên nếu ® là một hàm chuẩn hóa bất đối bất ki, lúc đó giá trị E` sẽ lớn hơn E, E'= [0 HOdr>E (11) Phương pháp biến đổi được áp dụng để xác định orbital tối ưu trong hàm sóng định thức đơn.

Chúng ta chọn một bộ cơ sở cho orbital mở rộng và hệ số Cus sau đó sé được thay đổi để cực tiểu hóa gid trị năng lượng E' tương ứng. Kết quả của E' có gần với giá trị của E hay không phụ thuộc vào cách chọn hàm sóng định thức đơn và bộ cơ sở thích hợp. Vì vậy hàm sóng đơn tốt nhất sẽ được tìm bằng cách cho cực tiểu hóa năng lượng E' theo hệ số œ„, lúc này phương trình biến đổi sẽ là : SP =0 (ải dc giá wii, w) (12) Cr 1. Ua điểm và nhược điểm của phương pháp Hartree-Fock * Ưu điểm: Uu điểm quan trọng nhất là rất dé áp dụng.

* Nhược điểm: Không được áp dụng rộng rãi cho những phân tử lớn, không cho năng lượng thấp nhất có thể đạt được. Các nhược điểm này xuất phát từ sự ấn định bất buộc cho cặp electron ứng với mỗi hàm sóng có thể dẫn đến sự phân bố của electron không đúng. Điều này thường xảy ra cho trường hợp phá vỡ liên kết. Thuyết hàm mật độ (Density-functional theory_ DFT) Hàm sóng electron của một phân tử có n electron phụ thuộc vào ca tọa độ không gian 3n và tọa độ quay n.

Trước đây toán tử Hamiltonian chỉ bao ham những số hạng không gian chứa một hoặc hai electron, nên năng lượng phân tử có thể được viết dưới đạng những tích phân chỉ bao gồm sáu tọa độ không gian. Tương tự như vậy, hảm sóng của một phân tử gồm nhiều electron chứa nhiều thông tin hơn mức cần thiết và hoàn toản thiếu hụt những ý nghĩa vật lý. Điều này đã thúc đầy cho việc tìm kiếm những hàm vừa chứa ít biến “i hon những hàm sóng kể trên vừa có thé sử dụng được để tính toán năng lượng và những tính chat khác. SVTH: Hồ Chí Thành _ lớp Hóa 4B (khóa 2002-2006) 7 Khóa luận tốt nghiệ GVHD: Professor Lawrence M.

Định lý Hohenberg-Kohn (The Hohenberg-Kohn Theorem) Năm 1964, Pierre Hohenberg va Walter Kohn đã chứng minh rằng: đối với nhitng phan tử ở trang thái cơ bản không suy biến thi năng lượng phân tử, ham sóng, và tat cả những tinh chất khác của electron phân tử ở trạng thái cơ bản chỉ phụ thuộc vào mật mí xác suất electron ở trạng thái cơ bản ø, (x, y, z), là một ham chỉ gồm có ba biến số [P. Người ta cho rằng năng lượng electron ở trạng thái cơ bản £_ là tết phiém hàm của ø, và viết E = £,[ ø, ]. trong đó dấu móc vuông biểu thị mối quan hệ phụ thuộc hàm số. Austrian-bern Amencen chemist Walter Kohn shared the 1998 Nobel Prize nm chemistry with British chernist John Pople.

Both Kohn and Pople were honored for developing techniques to calculate the behavior of groups of molecules. Thuvét ham mật độ có gắng tinh toán năng lượng £, và những tính chất khác của phân tử ở trạng thái cơ bản từ mật độ electron ở trạng thái cơ bản ø,. Thế nào là một phiếm hàm (functional)? Chúng ta hãy nhớ lại rằng hàm số (function) /£(x) là một phép gan một số với mỗi giá trị của biến x mà hàm f được định nghĩa cho nó. Thí dụ: ham f(x) =x’ + 1 gán số 10 với giá trị 3 của x và gan một số với mỗi giá trị của x.

Trong khi đó, phiém hàm (functional) F[/] là một phép gan một số với mỗi ham số ƒ. Thi dụ: phiém hàm F[/]= [/ˆ(x)/(x)& gin một số (được tìm thấy bởi phép tích phân của iff trên vô cùng) với mỗi ham bình phương khả tích f(x). Biến phân (variational integral) W[®]= (e|8|®) /(@\@) là một phiếm hàm của hàm biến thiên ® và cho một số ứng với mỗi giá trị có thể có của ®. Định lý Hohenberg-Kohn được chứng minh như sau: Ham sóng clectron ở trạng thái cơ ban yw, của một phân tử có n electron là một hàm riêng của toán tử Hamiltonian electron thuần khiết bong nguyên tử của phương trình sau: ñ=-2m >V;+Ši=v(r)+!>iy>lr,~ (13) v(n)=-D (14) a "@ Đại lượng v(r,), thể năng tương tác giữa electron i và hạt nhân.

phụ thuộc vào tọa độ Xj. Vị, Z4 của electron i và tọa độ của hạt nhân. Do phương trình sóng Sckrödinger của SVTH: Hồ Chí Thành _ lớp Hóa 4B (khóa 2002-2006) § Khóa luận tốt nghiệ GVHD: Professor Lawrence M.s Nguyễn Văn N electron được giải đổi với những vị trí cé định của hạt nhân, nên tọa độ của hạt nhân không biến đổi trong phương trình sóng Sclvödinger của Erwin Schrédinger A pioneer mn the area of quantum theory, Erwin Schrodinger is best known for hes mathematecal theory describing the wave mechanics of electrons. He and Britsh physicist Paul Dirac shared the 1933 Nobel Prize in physics for their contributions to the understanding of quantum mechanics.

Chinh vi thé, v(r,) trong phương trình sóng Schrédinger của electron chi còn là hàm của x, yi), z, (biểu thị dưới dang kí hiệu vector). Trong thuyết hàm mật độ.