DSpace at VNU: Kiểm nghiệm cơ chế phản ứng H2(k)+ Cl2(k) → 2HCl(k) bằng phương pháp tính hóa học lượng tử

DSpace at VNU: Kiểm nghiệm cơ chế phản ứng H2(k)+ Cl2(k) → 2HCl(k) bằng phương pháp tính hóa học lượng tử tài liệu, giáo...

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - VŨ THỊ NGỌC LAN KIỂM NGHIỆM CƠ CHẾ PHẢN ỨNG H2(k)+ Cl2(k) → 2HCl(k) BẰNG PHƢƠNG PHÁP TÍNH HĨA HỌC LƢỢNG TỬ LUẬN VĂN THẠC SĨ Hà Nội - 2015 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - VŨ THỊ NGỌC LAN KIỂM NGHIỆM CƠ CHẾ PHẢN ỨNG H2(k)+ Cl2(k) → 2HCl(k) BẰNG PHƢƠNG PHÁP TÍNH HĨA HỌC LƢỢNG TỬ Chuyên ngành: Hóa lý thuyết hóa lý Mã số: 60440119 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS Phạm Văn Nhiêu TS Vũ Việt Cường Hà Nội - 2015 PHẦN MỞ ĐẦU I Lí chọn đề tài Hóa học lượng tử (HHLT) ngành khoa học nghiên cứu hệ lượng tử dựa vào phương trình tắc CHLT Schrodinger đưa năm 1926 Với xuất HHLT nhanh chóng làm thay đổi quan niệm giới vi mô nhiều ngành khoa học đặc biệt nghiên cứu hóa học HHLT giúp tìm hiểu sâu, nghiên cứu vấn đề cốt lõi hóa học cấu trúc tính chất hóa lý chất Bên cạnh phát triển mạnh mẽ công nghệ máy tính đời phần mềm máy tính phục vụ việc tính HHLT MOPAC, HYPERCHEM, GAUSSIAN, GAUSSVIEW, VASP….Giúp cho việc nghiên cứu HHLT mở rộng Phần mềm tính tốn cho ta biết nhiều tham số xác định cấu trúc phân tử, tham số bề mặt năng, tham số đại lượng nhiệt động học như: ∆H; ∆G; đại lượng động học E, tốc độ phản ứng… HHLT làm sáng tỏ nhiều chế phản ứng hóa học, giải thích đắn quy luật hóa học, kiểm tra kết thực nghiệm Đặc biệt, cịn thể tính ưu việt nghiên cứu chuyển tiếp chất trung gian mà có thời gian xảy ngắn Khảo sát phản ứng hóa học điều kiện khác nhau, phản ứng khó, phản ứng độc hại hay khơng thể thực được, đồng thời tiết kiệm kinh tế Với mong muốn học hỏi, hiểu biết thêm HHLT, tìm hiểu chất phản ứng nên chọn đề tài nghiên cứu: “ Kiểm nghiệm chế phản ứng H2(k)+ Cl2(k) → 2HCl(k) phƣơng pháp tính hóa học lƣợng tử” II Mục đích nghiên cứu Sử dụng lý thuyết hóa học lượng tử phương pháp tính toán gần đúng, tốt áp dụng cho hệ nghiên cứu nhằm thu tham số cấu trúc, tần số dao động loại lượng… hệ chất tham gia, chất sản phẩm, Vũ Thị Ngọc Lan chất trung gian trạng thái chuyển tiếp qua thiết lập bề mặt đầy đủ để giải thích chế phản ứng Khi nghiên cứu chế phản ứng H2 + Cl2  HCl thấy phản ứng xảy theo giai đoạn Trong đó, có giai đoạn phát triển mạch mà nguyên tử công phân tử Vậy ngun tử cơng phân tử theo góc liên kết cho lợi lượng? Về mặt động học phản ứng, phương trình động học phản ứng biểu diễn nào? Với kết nghiên cứu, hy vọng thông số thu được sử dụng làm thơng tin đầu vào cho việc nghiên cứu nhiệt động học động lực học tiếp theo, đồng thời làm tài liệu tham khảo cho việc nghiên cứu hóa học thực nghiệm III Nhiệm vụ nghiên cứu Tìm hiểu sở lý thuyết, phương pháp tính tốn phần mềm liên quan hóa học lượng tử Sưu tầm báo tài liệu liên quan đến hệ chất nghiên cứu Lựa chọn phương pháp tính tốn tốt để kháo sát hệ chất nghiên cứu Dự đoán, kiểm nghiệm hướng phản ứng, xây dựng bề mặt đầy đủ lựa chọn chế phù hợp phản ứng IV Phƣơng pháp nghiên cứu Sử dụng phần mềm Gaussian 09 phần mềm hỗ trợ Gauss View, Chemcraft, Chemoffice… để nghiên cứu Sử dụng phương pháp DFT để tối ưu hóa cấu trúc tính lượng tương quan với hàm 6-31G cho nguyên tố Cl, H Cuối chúng tơi dùng kết tính tốn lượng tương quan để xây dựng bề mặt (PES) hệ chất nghiên cứu, từ rút kết luận quan trọng cho trình nghiên cứu Cấu trúc khóa luận gồm phần mở đầu, nội dung, kết luận, tài liệu tham khảo phụ lục Phần nội dung bao gồm chương: Chƣơng 1: Tổng quan Vũ Thị Ngọc Lan Chƣơng 2: Hệ chất nghiên cứu phƣơng pháp nghiên cứu Chƣơng 3: Kết thảo luận CHƢƠNG 1: TỔNG QUAN 1.1 Cơ sở lý thuyết hóa học lƣợng tử [1][4][9][10][12] 1.1.1 Phƣơng trình Schrodinger Phương trình Schrodinger đưa nhà vật lý người Áo – Schrodinger vào năm 1926 Sự biến đổi trạng thái hạt vi mô theo thời gian hệ lượng tử mơ tả phương trình Schrodinger có dạng tổng quát: 𝑖ħ 𝜕𝜓 𝜕𝑡 =Ĥ (1.1) Ở đây, Ψ(q,t) – Hàm sóng mơ tả trạng thái hệ lượng tử Đó hàm tọa độ q thời gian t Ĥ = 𝑇+ U = Với ∇2 = 𝜕2 𝜕2 𝜕2 𝜕𝑥 𝜕𝑦 𝜕𝑧 + + ħ2 2𝑚 ∇2 + 𝑈 : toán tử Laplace U: (năng lượng tương tác hạt lượng tử hệ) Trong trường hợp hệ không phụ thuộc thời gian: U = U(q) (q: tọa độ) Trong hệ kín hệ chuyển động mơi trường ngồi khơng đổi, tốn tử Hamilton Ĥ không phụ thuộc vào thời gian trùng với tốn tử lượng tồn phần Ĥ(q) trạng thái hệ gọi trạng thái dừng: ψ(q,t)→ψ(q) Phương trình Schrodinger trạng thái dừng: Ĥ(q)ψ(q) = Eψ(q) (1.2) E trị riêng lượng, ψ(q) hàm sóng (q tọa độ) Nghiệm phương trình (1.1) viết dạng: Ψ(q,t) = ψ(q).𝑒 Vũ Thị Ngọc Lan 𝑖𝐸𝑡 ħ (1.3) Những trạng thái mà hệ lượng tử có giá trị xác định, xác định phương trình (1.3) gọi trạng thái dừng (1.2) phương trình Schrodinger cho trạng thái dừng Các hệ lượng tử (nguyên tử, phân tử,…) xem hệ thống ổn định, bền vững theo thời gian, dùng (1.2) để xác định hàm sóng lượng chúng Với hệ lượng tử gồm M hạt nhân N electron, toán tử Hamilton xác định cho hệ hạt là: 𝑁 𝐻=− 𝑖=1 ∇ − 𝑖 𝑀 𝐴=1 ∇ − 2𝑀𝐴 𝐴 𝑁 𝑀 𝑖=1 𝐴=1 𝑍𝐴 + 𝑟𝑖𝐴 𝑁 𝑁 𝑖=1 𝑗 >1 + 𝑟ị 𝑀 𝑀 𝐴=1 𝐵>𝐴 𝑍𝐴 𝑍𝐵 𝑟𝐴𝐵 (1.4) Với A, B kí hiệu cho hạt nhân i, j kí hiệu cho electron hệ MA: khối lượng hạt nhân A ZA, ZB: điện tích hạt nhân A, B rij: khoảng cách electron thứ i j, riA: khoảng cách electron thứ i hạt nhân A, rAB: khoảng cách hạt nhân A B Số hạng thứ thứ hai phương trình (1.4) tốn tử động electron hạt nhân tương ứng, số hạng thứ ba tương tác hút Coulomb (Culong) electron hạt nhân, số hạng thứ tư thứ năm tương tác đẩy electron hạt nhân tương ứng Khi giải phương trình Schrodinger người ta thu hàm sóng ψ mơ tả trạng thái hệ lượng tử trạng thái đó, hệ lượng tử có lượng E Tuy vậy, nguyên tử, phân tử hệ vô phức tạp, nên thực tế phương trình Schrodinger khơng giải cách xác Để giải phương trình Schrodinger cho hệ lượng tử phức tạp, người ta sử dụng phương pháp gần 1.1.2.Sự gần Born – Oppenheirmer Toán tử Hamilton đầy đủ cho hệ phân tử:       H  T elec (r )  T nucl ( R)  V elecnucl ( R, r )  V elec (r )  V nucl ( R) Vũ Thị Ngọc Lan Phép gần Born – Oppenheirmer sử dụng để đơn giản hóa việc giải phương trình Schrodinger Vì khối lượng hạt nhân lớn nhiều khối lượng electron, hạt nhân chuyển động chậm so với electron, nên coi hạt nhân đứng yên Với gần này, động hạt nhân bỏ qua hạt nhân số Phương trình (1.2) viết lại là: Ĥeψe = Eeψe (1.5) Phép gần Born – Oppenheirmer xem electron chuyển động trường hạt nhân tĩnh điện (động hạt nhân 0), lúc số hạng tương tác tĩnh điện hạt nhân nguyên tử số nên tách số hạng khỏi toán tử Hamilton electron Lúc toán tử Hamilton trở thành: 𝑁 𝐻𝑒 = − 𝑖=1 ∇ − 𝑖 𝑁 𝑀 𝑖=1 𝐴=1 𝑍𝐴 − 𝑟𝑖𝐴 𝑁 𝑁 𝑖=1 𝑗 >1 𝑟𝑖𝑗 1.6 Hàm ψe phụ thuộc tọa độ electron tham số tọa độ hạt nhân Đối với hệ nhiều electron xảy tương tác electron Do vấn đề cốt lõi xử lý tương tác electron Vì khơng thể tính xác đại lượng nên thực tế người ta lấy giá trị trung bình 𝑈𝑒𝑒 nhằm mục đích làm cho phương trình Schrodinger giải mà kết đảm bảo độ xác dùng để giải thích kiện thực nghiệm Vì phân tử khơng có tính đối xứng cầu nên dùng phương pháp HF cho phân tử Roothaan thành cơng việc áp dụng phương trình Hartree – Fock cho MO xây dựng dạng tổ hợp tuyến tính obitan nguyên tử: 𝜓𝑖 = 𝑚 𝑗 =1 𝑐𝑖𝑗 𝜑𝑖𝑗 (1.7) cij hệ số khai triển m kích cỡ tập hàm sở, cij xác định phương pháp biến phân 1.1.3.Phƣơng pháp biến phân Mục đích phương pháp biến phân tìm cij gần với hàm sóng thực tế ψ ứng với lượng cực tiểu theo tập hàm sở chọn Năng lượng tính theo phương trình Schrodinger: Vũ Thị Ngọc Lan 𝐸= 𝜓 ∗Ĥ𝜓𝑑𝜏 (1.8) 𝜓∗𝜓𝑑𝜏 Với d𝜏 vi phân thể tích Nếu hàm ψ hàm chuẩn hóa phương trình (1.8) có dạng: 𝐸= 𝜓 ∗ Ĥ𝜓𝑑𝜏 (1.9) Áp dụng phương pháp biến phân, biểu diễn gần cho hàm sóng ψ dạng MO – LCAO: ψ = c1φ1 + c2φ2 + ….+ cnφn (1.10) Đặt (1.10) vào (1.8) nhận thấy trị số E phụ thuộc vào giá trị cj Theo nguyên lý biến phân, hệ số phải chọn cho trị số E cực tiểu Coi cj biến số E phụ thuộc biến số Điều kiện cực tiểu lượng: dE/dcj = (1.11) Xét với hệ hạt, thực vi phân (1.11) ta hệ phương trình tuyến tính với ψ = c1φ1 + c2φ2 𝐻11 − 𝐸𝑆11 𝑐1 + 𝐻12 − 𝐸𝑆12 𝑐2 = 𝐻21 − 𝐸𝑆21 𝑐1 + 𝐻22 − 𝐸𝑆22 𝑐2 = (1.12) Đối với hệ có n hạt, thực vi phân (1.11) hàm sóng (1.10) ta hệ phương trình tuyến tính có dạng : 𝐻11 − 𝐸𝑆11 𝑐1 + 𝐻12 − 𝐸𝑆12 𝑐2 + ⋯ + 𝐻1𝑛 − 𝐸𝑆1𝑛 𝑐𝑛 = 𝐻21 − 𝐸𝑆21 𝑐1 + 𝐻22 − 𝐸𝑆22 𝑐2 + ⋯ + 𝐻2𝑛 − 𝐸𝑆2𝑛 𝑐𝑛 = …………………………………………………………………… 𝐻𝑛1 − 𝐸𝑆𝑛1 𝑐1 + 𝐻𝑛2 − 𝐸𝑆𝑛2 𝑐2 + ⋯ + 𝐻𝑛𝑛 − 𝐸𝑆𝑛𝑛 𝑐𝑛 = (1.13) Hệ phương trình (1.13) viết gọn thành: 𝐻𝑖𝑗 − 𝐸𝑆𝑖𝑗 𝑐𝑗 = (1.14) Trong i số thứ tự phương trình, j số thứ tự số hạng Hệ phương trình (1.13) có nghiệm khác định thức kỉ lập từ phương trình hệ 0: Vũ Thị Ngọc Lan 𝐻11 − 𝐸𝑆11 𝐻12 − 𝐸𝑆12 … … 𝐻1𝑛 − 𝐸𝑆1𝑛 𝐻21 − 𝐸𝑆21 𝐻22 − 𝐸𝑆22 … … 𝐻2𝑛 − 𝐸𝑆2𝑛 ………………………………………………… ………………………………………………… = ………………………………………………… 𝐻𝑛1 − 𝐸𝑆𝑛1 𝐻𝑛2 − 𝐸𝑆𝑛2 … … 𝐻𝑛𝑛 − 𝐸𝑆𝑛𝑛 Hay: │Hij – ESij│= (1.15) Giải định thức kỉ ta tìm biểu thức lượng E Đặt giá trị E vào hệ phương trình (1.13) ta tìm hệ số cj từ suy hàm sóng cần tìm 1.1.4.Thuyết trƣờng tự hợp Hartree – Fock Tất phương pháp tính obitan phân tử đại (ab initio bán kinh nghiệm) sử dụng phương pháp tính gần Hartree – Fock (HF) để giải gần hàm sóng phân tử Từ quan điểm vật lý trường hiệu dụng trung bình hóa electron hợp hút hạt nhân đẩy trung bình hóa tất electron khác sinh Hàm sóng thích hợp mô tả trạng thái electron hàm obitan – spin: ψ(x) xác định tích hàm khơng gian ψ(q) với hàm spin χ(σ) (hàm spin 1 2 kí hiệu α ms = + β ms = − ) Ψ(x) = φ(q).χ(σ) (1.16) Với x tọa độ obitan – spin; σ tọa độ spin; q tọa độ không gian Hàm sóng phản đối xứng đơn giản sử dụng để mô tả trạng thái hệ N electron định thức Slater: 𝜓 𝑥1 , 𝑥2 , … , 𝑥𝑛 = 𝑁! 𝜓1 𝑥1 𝜓2 𝑥1 … 𝜓𝑁 𝑥1 𝜓1 𝑥2 𝜓2 𝑥2 … 𝜓𝑁 𝑥2 …………………………… 𝜓1 𝑥𝑁 𝜓2 𝑥𝑁 … … 𝜓𝑁 𝑥𝑁 (1.17) Với số i ψi số lượng tử; xi tọa độ obitan – spin i Theo nguyên lý biến phân “Hàm sóng tốt xác định theo định thức Slater hàm sóng ứng với lượng cực tiểu: 𝐸 = 𝜓│𝐻│𝜓 ” Áp dụng nguyên lý biến phân ta hệ phương trình Hartree – Fock (HF) Phương trình HF xác định obitan – spin tối ưu có dạng: Vũ Thị Ngọc Lan ℎ 𝑖 𝜓 𝑥𝑖 = 𝜀𝑖 𝜓 𝑥𝑖 (1.18) Ở h 𝑖 toán tử Hamilton hiệu dụng electron trường trung bình hạt cịn lại, gọi tốn tử Fock: ℎ 𝑖 = − ∇2𝑖 − 𝑀 𝑍𝐴 𝐴=1 𝑟 𝑖𝐴 + 𝑉 𝐻𝐹 𝑖 (1.19) VHF(i): trung bình electron thứ i có mặt electron khác Phương pháp gần HF thay nhiều hệ electron phức tạp N hệ electron, đẩy 2e xử lý trung bình hóa Thế HF (VHF (i)) phụ thuộc vào obitan – spin electron khảo sát Phương trình HF (1.18) khơng tuyến tính theo phương pháp “lặp lặp lại” hay gọi phương pháp trường tự hợp SCF (Self Consistent field) Trong phương pháp này, từ hàm obitan – spin dự đoán ban đầu, electron khác (VHF(i)) Tiếp theo giải phương trình (1.18) để xác định obitan – spin Sau đó, dùng obitan để tính VHF(i) Lặp lại q trình trường VHF(i) obitan – spin không cịn thay đổi Kết tìm là: 𝑛 𝑛 𝐸𝜋 = 𝑛 𝜀𝑖 − 𝑖=1 2𝐽𝑖𝑗 − 𝐾𝑖𝑗 (1.20) 𝑖=1 𝑗 =1 Trong đó: Jij Kij tích phân electron Jij = 𝜓𝑖∗ 𝜓𝑗 𝑟12 𝜓𝑖 𝜓𝑗 𝑑𝜏1 𝑑𝜏2 lượng HF, tương đương với tương tác Coulomb cổ điển Kij = 𝜓𝑖∗ 𝜓𝑗 𝑟12 𝜓𝑖 𝜓𝑗 𝑑𝜏1 𝑑𝜏2 gọi lượng trao đổi HF Chất lượng kết HF phụ thuộc vào kích cỡ chất lượng hệ hàm sở Tuy nhiên độ xác HF bị hạn chế việc sử dụng hiệu ứng liên hỗ trợ trung bình 1.1.5.Phƣơng trình Roothaan Phương trình Roothaan áp dụng phân tử khơng có trường Coulomb đối xứng cầu Roothaan (1951) sử dụng tập hàm sở để mở rộng Vũ Thị Ngọc Lan TÀI LIỆU THAM KHẢO Tiếng Việt H Eyring, J Walter, G.E Kimball (1948), Hóa học lượng tử (1976)bản dịch tiếngViệt, Nhà xuất Khoa học Kỹ thuật Việt Nam, Hà Nội Trần Thành Huế (2003), Hóa học đại cương, Nhà xuất giáo dục, Hà Nội Trần Thành Huế (2006), Tư liệu hóa học 10, Nhà xuất giáo dục, Hà Nội Nguyễn Đình Huề, Nguyễn Đức Chuy (2003), Thuyết lượng tử nguyên tử phân tử (Tập 1,2), NXBGD, Hà Nội Đồn Minh Hùng (2015), Khảo sát thơng số nhiệt động, đường phản ứng gốc tự Etinyl (C2H) với phân tử acryonitrin (C3H3N) pha khí phương pháp tính hóa học lượng tử, Luận văn thạc sĩ, Đại học Sư phạm Hà Nội Lê Văn Huỳnh (2014), Hóa học nguyên tố, Nhà xuất Khoa học tự nhiên công nghệ, Hà Nội Nguyễn Hà Mi (2012), Khảo sát số dẫn xuất halogen, ancol, phenol axit cacboxylic phương pháp hóa học lượng tử, Luận văn thạc sĩ, Đại học Khoa học tự nhiên Phạm Thị Thu Ngọc (2014), Nghiên cứu lý thuyết chế phản ứng N2O + H2 pha khí xúc tác cluster Rh5, Luận văn thạc sĩ, Đại học Sư phạm Hà Nội Hồng Nhâm (2000), Hóa vơ (Tập 2), NXBGD, Hà Nội 10 Trần Văn Nhân, Nguyễn Thạc Sửu, Nguyễn Văn Tuế, (2006), Hóa lí (Tập 2,3), NXBGD, Hà Nội 11 Nguyễn Hữu Phú (2006), Hóa lý Hóa keo, NXBKHKT, Hà Nội 12 Lâm Ngọc Thiềm, Phạm Văn Nhiêu, Lê Kim Long (2008), Cơ sở hóa học lượng tử, NXBKHKT, Hà Nội 13 Đào Đình Thức (2005), Cấu tạo nguyên tử liên kết hóa học (Tập 1,2), Nhà xuất giáo dục, Hà Nội Vũ Thị Ngọc Lan 14 Nguyễn Ngọc Trí (2015), Bước đầu nghiên cứu động học phản ứng đơn phân tử phụ thuộc áp suất phương pháp tính hóa học lượng tử, Luận văn thạc sĩ, Đại học Sư phạm Hà Nội 15 Hồ Ngọc Tuấn (2015), Nghiên cứu lí thuyết chế phản ứng gốc etinyl với phân tử etanol pha khí lý thuyết phiếm hàm mật độ, Luận văn thạc sĩ, Đại học Sư phạm Hà Nội 16 Đặng Ứng Vận (1998), Tin học ứng dụng hóa học, Nhà xuất giáo dục, Hà Nội 17 Đào Hữu Vinh, Nguyễn Duy Ái (2014), Tài liệu chuyên Hóa học 10 (Tập 2), Nhà xuất giáo dục, Hà Nội Tiếng Anh 18 Eleen Frisch, Hrantchian, P Hrat Roy D Dennington II, Todd A Keith, John Millam,…(2009), GaussView Reference, Gaussian, Inc 19 Gloria A.A Saracino, Roberto Improta, Vincenzo Barone (2003), “Absolute pKa determination for cacboxylic acids using density functional theory and the polarizable continuum model”, Chemical Physics Letters, 373, pp 411415 20 John A Keith, Emily A Carter (2012), “Quantum Chemical Benchmarking, Validation, and Prediction of Acidity Constants for Substituted Pyridinium Ions and Pyridinyl Radicals”, Journal of Chemical Theory and Computation, 8, pp 3187-3206 21 Kristin S Along, George C Sheilds (2010), “Chapter Theoretical Calculations of Acid Dissociation Constants: A Review Artical”, Annual Reports in Computation Chemistry, Volum 6, pp 113-138 22 Matthew D Liptak, George C Shields (2001), “Accurate pKa Calculation for Carboxylic Acids Using Complete Basis Set and Gaussian-n Models Combined with CPCM Continuum Solvation Methods”, J Am Chem Soc, 123, pp 7314-7319 Vũ Thị Ngọc Lan 23 Vyacheslav S Bryantsev, Mamadou S Diallo, Adri C T van Duin, William A Goddard III (2009), “Evaluation of B3LYP, X3LYP and M06-Class Density Functionals for Predicting the Binding Energies of Neutral, Protonated, and Deprotonated Water Clusters”, Journal of Chemical Theory and Computation, 5, pp 1016-1026 Vũ Thị Ngọc Lan ... tài nghiên cứu: “ Kiểm nghiệm chế phản ứng H2(k)+ Cl2(k) → 2HCl(k) phƣơng pháp tính hóa học lƣợng tử? ?? II Mục đích nghiên cứu Sử dụng lý thuyết hóa học lượng tử phương pháp tính tốn gần đúng,...ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - VŨ THỊ NGỌC LAN KIỂM NGHIỆM CƠ CHẾ PHẢN ỨNG H2(k)+ Cl2(k) → 2HCl(k) BẰNG PHƢƠNG PHÁP TÍNH HĨA HỌC LƢỢNG TỬ Chun ngành: Hóa. .. cứu động học phản ứng đơn phân tử phụ thuộc áp suất phương pháp tính hóa học lượng tử, Luận văn thạc sĩ, Đại học Sư phạm Hà Nội 15 Hồ Ngọc Tuấn (2015), Nghiên cứu lí thuyết chế phản ứng gốc etinyl