Nghiên cứu lý thuyết khả năng phản ứng của axit isocyanic HNCO dựa vào độ mềm cục bộ và hàng rào thế năng

Khả năng phản ứng của axit isocyanic HNCO với các gốc CH3 , NH2 , OH được nghiên cứu bằng phương pháp phiếm hàm mật độ (DFT) B3LYP, với mục đích chúng tôi so sánh khả năng phản ứng theo hai phương pháp là độ mềm cục bộ và hàng rào năng lượng.

LIÊN NGÀNH HĨA HỌC - CƠNG NGHỆ THỰC PHẨM NGHIÊN CỨU LÝ THUYẾT KHẢ NĂNG PHẢN ỨNG CỦA AXIT ISOCYANIC HNCO DỰA VÀO ĐỘ MỀM CỤC BỘ VÀ HÀNG RÀO THẾ NĂNG THEORETICAL STUDY ON THE REACTION ABILITY OF HNCO ISOCYANIC ACID BASED ON LOCAL SOFTNESS AND POTENTIAL ENERGY Vũ Hoàng Phương, Lê Văn Thuỷ Email: hphuong_hp@yahoo.com Trường Đại học Sao Đỏ Ngày nhận bài: 19/01/2018 Ngày nhận sửa sau phản biện: 17/6/2018 Ngày chấp nhận đăng: 28/6/2018 Tóm tắt Khả phản ứng axit isocyanic HNCO với gốc CH3, NH2, OH nghiên cứu phương pháp phiếm hàm mật độ (DFT) B3LYP, với mục đích so sánh khả phản ứng theo hai phương pháp độ mềm cục hàng rào lượng Kết cho thấy sử dụng độ mềm để dự đoán phản ứng trường hợp chưa tối ưu Tính tốn hàng rào lượng phù hợp, với vị trí cơng vào H thuận lợi nhất, O thuận lợi Như vậy, hệ phản ứng khơng có phù hợp tốt hàng rào lượng độ mềm cục bộ, nên dựa vào hàng rào lượng thực nghiệm để dự đoán khả phản ứng tác nhân vào vị trí khác phân tử Từ khóa: Axit isocyanic HNCO; độ mềm cục bộ; thuyết phiếm hàm mật độ (DFT) Abstract The reactivity of isocyanic acid HNCO with CH3, NH2, OH radicals have been studied by the density functional Theory (DFT) using B3LYP functional, for the purpose we compared reaction ability by two methods: local softness and potential energy The results showed that the use of softness to predict the response in this case is not optimal Calculated results indicate that potential energy of reactions well agee, with H atom is most favorable, O atom is less favorable So, reactions are not good fit between the potential energy and the local softness, so using the potential energy and experiment to predict the reactivity of the agents at different positions in the molecule Keywords: Isocyanic acid HNCO; local softness; density functional theory (DFT) GIỚI THIỆU CHUNG Ngày nay, với phát triển khoa học kĩ thuật, đặc biệt công nghệ thơng tin giúp cho mơn hóa học lượng tử phát triển nhanh, xây dựng phần mềm chun dụng dùng tính tốn hóa học lượng tử Thực tế cho thấy kết tính tốn lý thuyết thực nghiệm phù hợp với Từ đó, tính tốn hóa học lượng tử công cụ đắc lực việc nghiên cứu q trình hóa học Song song với thành tựu vô to lớn đạt lĩnh vực khoa học kỹ thuật, người phải đối mặt với nguy to tớn Một nguy nhiễm mơi trường chất thải công nghiệp Tác nhân gây ô nhiễm chất, ngun tố hóa Người phản biện: PGS.TS Ngơ Sỹ Lương TS Hồng Thị Hịa học, mà số oxit nitơ Trong q trình đốt cháy nhiên liệu hóa thạch khơng khí có sinh sản phẩm trung gian axit isocyanic - công thức phân tử HNCO Axit isocyanic (HNCO) tìm thấy tầng trung vũ trụ lần vào năm 1972 hai nhà bác học Lew Snyder David Bulh Trong chổi tìm thấy Hyakutake sau Hale - Bopp [1] HNCO điều chế vào năm 1830 nhà bác học người Đức Justus von Liebig Wohler Là chất không màu, dễ bay độc, sôi 23,5oC Axit isocyanic hợp chất hóa học bền, phân tử nhỏ có chứa bốn nguyên tố thường tìm thấy hợp chất hóa học hữu sinh học HNCO có momen lưỡng cực tương đương với H2O, axit hữu mạnh với giá trị pKa vào khoảng 3,5, nhường proton môi trường nước Tạp chí Nghiên cứu khoa học - Đại học Sao Đỏ, ISSN 1859-4190 Số 2(61).2018 89 NGHIÊN CỨU KHOA HỌC điều kiện định, đặc biệt có mặt bazơ mạnh HNCO lại phản ứng với thành phần khí cháy CH3, NH2, OH Như vậy, HNCO sản phẩm trung gian tồn chế q trình đốt cháy lại giúp giảm bớt nhiễm khí NOx sinh trình đốt cháy nhiên liệu hóa thạch [2], q trình cơng nghiệp hóa Ví dụ phản ứng HNCO với gốc OH, gốc quan trọng phản ứng cháy tính tốn James [3] Sản phẩm phản ứng NH2 CO2 phù hợp với thực nghiệm Hay phản ứng HNCO với NH2 nghiên cứu John [4] f(r) = δρ(r) / δN | v(r) = δµ / δv(r) |v s(r) = f(r) S Hiện nay, có nhiều cơng trình nghiên cứu HNCO [5, 6, 7] Để góp phần cung cấp thêm thơng tin cho trình nghiên cứu sâu sau này, tiến hành nghiên cứu khảo sát khả phản ứng gốc CH3, NH2, OH vào vị trí khác HNCO dựa vào độ mềm cục [8] hàng rào lượng Việc nghiên cứu giúp ta hiểu cách khống chế sản phẩm độc hại gây ô nhiễm môi trường để góp phần bảo vệ mơi trường 2η = 1/S = (I – A) PHƯƠNG PHÁP NGHIÊN CỨU 3.1 Khả phản ứng dựa vào độ mềm cục nguyên tử Dùng phương pháp phiếm hàm mật độ (DFT) B3LYP với hàm sở 6-311++G(d,p) 6-311++G(3df,2p) để tối ưu hóa cấu trúc tính lượng chất tham gia phản ứng, trạng thái chuyển tiếp Trong đó, trạng thái chuyển tiếp kiểm tra việc phân tích tần số dao động tọa độ thực (IRC) phản ứng Tính hàm Fukui theo phương pháp tính điện tích khác độ mềm phân tử độ mềm cục nguyên tử Sau so sánh khả phản ứng HNCO với tác nhân theo hai phương pháp Ứng dụng quan trọng độ mềm dự đốn giải thích khả phản ứng nguyên tử phân tử Thế hóa học (µ), độ âm điện (χ), độ cứng (η), độ mềm(S) tính theo cơng thức: χ = - µ = - δE / δN |v 2η = 1/S = - δ2E /δN2 |v = - δµ / δN|v [13] đó: E, N, v lượng, số lượng điện tử, hệ I, A lượng ion hóa, lực electron 2.2 Phần mềm tính tốn Để nghiên cứu khả phản ứng gốc CH3, NH2, OH vào vị trí khác HNCO phương pháp hóa học lượng tử, sử dụng phần mềm Gaussian 09 [14] phần mềm hỗ trợ Gaussview, Chemoffice KẾT QUẢ VÀ THẢO LUẬN Cấu trúc tối ưu phân tử HNCO tác nhân sau: 3.1.1 Cấu trúc phân tử HNCO tác nhân a HNCO Công thức cấu tạo tương ứng H–N=C=O Tối ưu hóa cấu trúc phần mềm Gauss view sau: 2.1 Hàm Fukui (fo) độ mềm cục (so) Trong hóa học, để xét khả phản ứng phân tử nguyên tử người ta gán cho phân tử nguyên tử trị số biểu khả phản ứng nguyên tử, phân tử Để làm điều này, Yang Mortier [9,10,11,12] đưa hàm Fukui tính nguyên tử : fK + = qK (N+1) – qK(N) f K- = qK (N) – qK (N–1) đó: qK: mật độ điện tích ngun tử K phân tử xét Hàm Fukui f(r) độ mềm cục s(r) có mối liên hệ với theo biểu thức sau: Hình Phân tử HNCO Bảng Độ dài liên kết phân tử Liên kết Tính tốn (Ao) H-N 1,0062 N-C 1,2114 C-O 1,1655 90 Tạp chí Nghiên cứu khoa học - Đại học Sao Đỏ, ISSN 1859-4190 Số 2(61).2018 LIÊN NGÀNH HÓA HỌC - CÔNG NGHỆ THỰC PHẨM b Các tác nhân 3.1.2 Hàm Fukui, độ mềm cục tính theo - Tác nhân OH: Độ dài liên kết O-H 0,976Ao điện tích Mulliken - Tác nhân NH2: Độ dài liên kết N-H 1,0299Ao - Tác nhân CH3: Độ dài liên kết C-H 1,0779Ao Kết tính tốn độ mềm phân tử, độ mềm cục nguyên tử hệ nghiên cứu theo điện tích Mulliken thể bảng Bảng Hàm Fukui (f o), độ mềm cục nguyên tử (so) tính theo điện tích Mulliken Cấu trúc Độ mềm phân tử (S) 2,685 Nguyên tử H N C O C fo 0,963 0,026 -0,131 0,142 0,678 so 1,704 0,046 0,251 1,820 HNCO 1,769 CH3 NH2 OH 2,203 1,846 N O 0,751 0,859 1,654 1,586 Qua bảng ta thấy giá trị độ mềm so giảm theo thứ tự H, O, N, C, từ có kết luận sau: tác nhân CH3, OH, NH2 công vào vị trí phân tử HCNO chúng ưu tiên cơng vào vị trí: Vị trí ưu tiên số 1: vị trí H Vị trí ưu tiên số 2: vị trí O Vị trí ưu tiên số 4: vị trí C Mặt khác, độ mềm S C, N, O ứng với tác nhân CH3, OH, NH2 giảm dần, trị số fo lại tăng dần nên so sánh độ mềm cục C, N, O tương ứng tác nhân, khơng dự đốn khả phản ứng tác nhân với HNCO [11] 3.1.3 Hàm Fukui, độ mềm cục tính theo điện tích ESP Vị trí ưu tiên số 3: vị trí N Bảng Hàm Fukui (fo), độ mềm cục nguyên tử (so) tính theo điện tích ESP Cấu trúc Độ mềm phân tử HNCO 1,769 CH3 NH2 OH 2,685 2,203 1,846 Nguyên tử fo so H N C O C N O 0,637 -0.167 0,144 0,387 1,466 0,816 0,865 1,127 0,255 0,685 4,000 1,796 1,597 Giá trị độ mềm giảm theo thứ tự H, O, C, N, từ Vị trí ưu tiên số 2: vị trí O có kết luận sau: Khi tác nhân CH3, OH, NH2 Vị trí ưu tiên số 3: vị trí C cơng vào vị trí phân tử HCNO Vị trí ưu tiên số 4: vị trí N chúng ưu tiên cơng vào vị trí: Vị trí ưu tiên số 1: vị trí H 3.1.4 Hàm Fukui, độ mềm cục tính theo điện tích NBA Bảng Hàm Fukui (fo), độ mềm cục nguyên tử (so) tính theo điện tích NBA Cấu trúc Độ mềm phân tử Nguyên tử fo so 0,396 0,330 0,035 0,239 0,915 0,700 0,584 0,062 0,423 2,457 HNCO 1,769 CH3 2,685 H N C O C NH2 2,203 N 0,916 2,018 OH 1,846 O 0,937 1,73 Tạp chí Nghiên cứu khoa học - Đại học Sao Đỏ, ISSN 1859-4190 Số 2(61).2018 91 NGHIÊN CỨU KHOA HỌC Giá trị độ mềm giảm theo thứ tự H, N, O, C, từ có kết luận sau: Khi tác nhân CH3, OH, NH2 cơng vào vị trí phân tử HCNO chúng ưu tiên cơng vào vị trí: Vị trí ưu tiên số 1: vị trí H Vị trí ưu tiên số 2: vị trí N Vị trí ưu tiên số 3: vị trí O Vị trí ưu tiên số 4: vị trí C Từ kết quả, xét khả phản ứng dựa vào độ cứng mềm tính phương pháp tính điện tích khác Chúng ta có đồ thị so sánh khả phản ứng nguyên tử phân tử HNCO tác nhân CH3, OH, NH2 cơng sau: Hình Đồ thị so sánh khả phản ứng nguyên tử HNCO theo điện tích Mulliken, ESP, NBA Từ kết đến kết luận sau: Khi tính độ mềm nguyên tử chất hệ phản ứng phương pháp tính điện tích khác (Mulliken, ESP, NBA) nhận khả phản ứng tác nhân (CH3, NH2, OH) vào vị trí phân tử HNCO theo chiều sau: - Khi tính theo Mulliken, vị trí ưu tiên theo thứ tự sau: H, C, N, O - Khi tính theo ESP, vị trí ưu tiên theo thứ tự sau: H, O, C, N - Khi tính theo NBA, vị trí ưu tiên theo thứ tự sau: H, O, N, C Như vậy, với ba cách tính điện tích khác ta thu ba kết khác Điều giải thích HNCO điện tích phân bố cách khơng định vị mà phân bố toàn phân tử Do đó, xét phản ứng dựa vào độ mềm cục ngun tử với phương pháp tính điện tích: TS CH3-H TS CH3-N - Khả phản ứng tác nhân khác vào vị trí khác phân tử HNCO theo chiều định - Mỗi phương pháp tính điện tích khác ta nhận kết khả phản ứng vị trí khác Do xét khả phản ứng dựa vào độ mềm hệ HNCO kết thu có độ xác cao khơng cao Để kiểm nghiệm kết luận này, tiến hành xét khả phản ứng gốc vào vị trí khác phân tử HNCO dựa vào lượng phản ứng 3.2 Khả phản ứng dựa vào hàng rào lượng 3.2.1 Tác nhân CH3 Kết tối ưu hóa hình học theo phương pháp B3LYP/6-311++G(d,p) thu trạng thái chuyển tiếp mà CH3 cơng vào vị trí khác HNCO: TS CH3-C TS CH3-O Hình Trạng thái chuyển tiếp CH3 công vào vị trí khác HNCO (với trạng thái chuyển tiếp kí hiệu TS) 92 Tạp chí Nghiên cứu khoa học - Đại học Sao Đỏ, ISSN 1859-4190 Số 2(61).2018 LIÊN NGÀNH HĨA HỌC - CƠNG NGHỆ THỰC PHẨM Bảng Năng lượng tổng lượng điểm khơng đơn vị (au), lượng tương đối tính theo đơn vị kJ/mol tác nhân CH3 tính theo phương pháp B3LYP/6-311++G(d,p) Cấu trúc Năng lượng tổng ZPE Năng lượng tổng + ZPE Năng lượng tương đối (au) (KJ/mol) (au) (au) HNCO -168,751413 0,021254 CH3 -39,857785 0,029633 HNCO + CH3 -208,609198 0,050887 -208,558311 0,0 TS CH3-H -208,592345 0,050684 -208,541661 43,7 TS CH3-N -208,585805 0,055476 -208,530329 73,5 TS CH3-C -208,581453 0,055713 -208,52574 85,5 TS CH3-O -208,569265 0,055952 -208,513313 118,1 Ta biểu diễn tương quan giá trị ΔE vị trí lên giản đồ sau: trạng thái chuyển tiếp chất phản ứng) thấp chứng tỏ khả hình thành liên kết dễ, thấp ΔE= 43,7 (KJ/mol) TS CH3-H, cao TS CH3-O, thứ tự ưu tiên là: Vị trí số 1: H Vị trí số 2: N Vị trí số 3: C Vị trí số 4: O 3.2.2 Tác nhân NH2 Hình Giản đồ so sánh lương tương quan trạng thái tác nhân CH3 Năng lượng hoạt hóa (năng lượng tương đối TS H2N-H TS H2N-N Kết tối ưu hóa hình học theo phương pháp B3LYP/6-311++G(d,p) thu trạng thái chuyển tiếp mà NH2 công vào vị trí khác HNCO: TS H2N-C TS H2N-O Hình 5.Trạng thái chuyển tiếp NH2 công vào vị trí khác HNCO Từ kết mặt lượng, thiết lập giản đồ so sánh sau: Dựa vào chênh lệch giá trị lương hoạt hóa vị trí NH2 cơng vào phân tử HNCO (thấp ΔE= 9,5 (kJ/mol) TS H2N-H), thứ tự ưu tiên là: Vị trí số 1: H Vị trí số 2: C Vị trí số 3: N Vị trí số 4: O 3.2.3 Tác nhân OH Hình Giản đồ so sánh lương tương quan trạng thái tác nhân NH2 Kết tối ưu hóa hình học theo phương pháp B3LYP/6-311++G(d,p) thu trạng thái chuyển tiếp mà NH2 cơng vào vị trí khác HNCO: Tạp chí Nghiên cứu khoa học - Đại học Sao Đỏ, ISSN 1859-4190 Số 2(61).2018 93 NGHIÊN CỨU KHOA HỌC TS HO-H TS HO-N TS HO-C TS HO-O Hình Trạng thái chuyển tiếp NH2 cơng vào vị trí khác HNCO Từ kết mặt lượng thiết lập giản đồ so sánh sau: Hình Giản đồ so sánh lương tương quan trạng thái tác nhân OH Dựa vào chênh lệch giá trị lượng hoạt hố vị trí tác nhân OH cơng vào phân tử HNCO thứ tự ưu tiên là: Vị trí số 1: H Vị trí số 2: C Vị trí số 3: N Vị trí số 4: O Kết so sánh giá trị ΔE tác nhân công vào vị trí phân tử HNCO trình bày bảng Bảng So sánh giá trị lượng hoạt hóa phản ứng HNCO với gốc Tác nhân Phân tử HNCO H N C O TS CH3 43,7 73,5 85,5 118,1 TS NH2 9,5 105,1 57,6 193,5 TS OH 0,4 127,9 25,2 269,0 Qua bảng phân tích kết tính tốn kiện nhiệt động, nhận thấy tác nhân hầu hết công vào vị trí H Đó H có độ âm điện (2,2) thấp nguyên tố khác phân tử HNCO (N: 3,04; O: 3,44; C: 2,55) nên liên kết dễ dàng bị phân cực hóa tác nhân có độ âm điện lớn hơn, tạo điều kiện thuận lợi cho việc hình thành liên kết Trong đó, tác nhân OH cơng thuận lợi với giá trị ∆E thấp 0,4 (kJ/mol), sau đến NH2 (với giá trị ∆E = 9,5 (kJ/mol)) Điều phân tử HNCO nguyên tử H liên kết với N – nguyên tố có độ âm điện lớn – nên H liên kết chặt chẽ với N, cần phải có tác nhân với trung tâm phản ứng nguyên tố có độ âm điện lớn để làm phân li liên kết H – N Với vị trí N O nguyên tố có độ âm điện lớn, theo chiều tăng dần độ âm điện trung tâm phản ứng (từ CH3, NH2, OH), giá trị ∆E tăng dần nguyên tử C, N, O (của tác nhân) tích điện dấu với nguyên tử N O (của phân tử HNCO) Điện tích nguyên tử (đơn vị C) trạng thái chuyển tiếp sau: TS H3C-N (C: -0,377, N: -0,693); TS H2N-N (N: -0,480, N(HNCO): -0,563); TS HO-N (O: -0,444, N: -0,484) Từ thấy, theo chiều tăng độ âm điện điện tích trung tâm phản ứng N gần nhau, điều làm cho lực đẩy trung tâm phản ứng N tăng dần, tức hình thành trạng thái chuyển tiếp khó dần hay giá trị ∆E tăng dần Đối với vị trí O tương tự Với vị trí C H, giá trị ∆E giảm dần theo thứ tự trung tâm phản ứng tăng dần độ âm điện, khả phân cực hóa tăng dần Khi tác nhân cơng vào vị trí bị che khuất phân tử HNCO ta thấy góc vị trí bị thu hẹp để tạo cấu trúc thuận lợi cho hình thành trạng thái chuyển tiếp Khi trạng thái chuyển tiếp hình thành, liên kết nguyên tử vị trí với ngun tử xung quanh bị kéo giãn Đặc biệt, nguyên tử có độ âm điện bé bị kéo giãn nhiều KẾT LUẬN Đã tính tốn tính khả phản ứng gốc - CH3, -NH2, -OH vào vị trí nguyên tử H, N, C, O phân tử HNCO dựa vào hàng rào lượng phương pháp sử dụng độ mềm cục Việc sử dụng độ mềm để dự đốn phản ứng 94 Tạp chí Nghiên cứu khoa học - Đại học Sao Đỏ, ISSN 1859-4190 Số 2(61).2018 LIÊN NGÀNH HĨA HỌC - CƠNG NGHỆ THỰC PHẨM trường hợp cho kết phù hợp với vị trí H, vị trí N, O, C kết chưa tối ưu Nguyên nhân phân bố điện tích phân tử HNCO khơng định vị mà trải toàn phân tử Đối với hệ khơng có phù hợp tốt hàng rào lượng độ mềm cục bộ, nên dựa vào hàng rào lượng thực nghiệm để dự đoán khả phản ứng tác nhân vào vị trí khác phân tử TÀI LIỆU THAM KHẢO [1] Snyder, L E.; Buhl, D (1972) Interstellar Isocyanic Acid Astrophys J Vol 177, pp 619−623 [2] Yue Ma, Xiaodong Wu, Junyu Zhang, Rui Ran, Duan Weng (2018) Urea-related Reactions and Their Active Sites over Cu-SAPO-34: Formation of NH3 and Conversion of HNCO Applied Catalysis B: Environmental Vol 227, pp 198-208 [3] James A Miller, Craig T Bowman (1991) Kinetic modeling of the reduction of nitric oxide in combustion products by isocyanic acid International Journal of Chemical Kinetics Vol 23, pp 289-313 [4] John D Mertens, Katharina Kohse-Höinghaus, Ronald K Hanson, Craig T Bowman (1991) A shock tube study of H + HNCO → NH2 + CO International Journal of Chemical Kinetics Vol 23, pp 655-668 [5] Luciano Fusina,Ian M Mills (1981) The harmonic force field and rz structure of HNCO Journal of Molecular Spectroscopy Vol 86, pp 488-498 [6] Donghui Quan, Eric Herbst, Yoshihiro Osamura, and Evelyne Roueff (2010) Gas-Grain modeling of isocyanic acid (HNCO), cyanic acid (HOCN), fulminic acid (HCNO), and isofulminic acid (HONC) in assorted interstellar environments The Astrophysical Journal Vol 725, pp 21012109 [7] Valerio Lattanzit, Sven Thorwirth, Carl A Gottlieb and Michael C McCarthy (2012) Two isomers of Protonated Isocyanic Acid J Phys Chem Lett Vol 3, No 23, pp 3420-3424 [8] Yang, W.; Parr, R G (1985). Hardness, softness, and the fukui function in the electronic theory of metals and catalysis Proc Natl Acad Sci. Vol. 82, No 20, pp 6723-6726 [9] Yang, W.; Mortier, W.J (1986) The use of global and local molecular parameters for the analysis of the gas-phase basicity of amines J Am Chem Soc. Vol 108, pp 5708-5711 [10] Fuentealba  and  R Contreras  (2002) Fukui funtion in chemistry Reviews of Modern Quantum Chemistry, pp 1013 - 1052 [11] Asit K Chandra and Minh Tho Nguyen (2002) Use of Local Softness for the Interpretation of Reaction Mechanisms Int J Mol.  Vol 3, No 4, pp 310-323 [12] Miquel Torrent-Sucarrat, Frank De Proft, Paul Geerlings and Paul W Ayers (2008) Do the Local Softness and Hardness Indicate the Softest and Hardest Regions of a Molecule? Chem Eur J Vol 14, pp 8652 - 8660 [13] Parr, R G.; Pearson, R G (1983) Absolute hardness: companion parameter to absolute electronegativity J Am Chem Soc. Vol 105, pp 7512-7516 [14] M J Frisch, G W Trucks, H B Schlegel, J A Pople (2009) Gaussian, Inc Pittsburgh PA [15] Hue Minh Thi Nguyen, Asit K Chandra, Jozef Peeters, and Minh Tho Nguyen (2004) Use of DFT-Based Reactivity Descriptors for Rationalizing Radical Reactions:  A Critical Analysis Journal of Physical Chemistry Vol 108, No 3, pp 484-489 Tạp chí Nghiên cứu khoa học - Đại học Sao Đỏ, ISSN 1859-4190 Số 2(61).2018 95 ... nghiên cứu HNCO [5, 6, 7] Để góp phần cung cấp thêm thơng tin cho q trình nghiên cứu sâu sau này, tiến hành nghiên cứu khảo sát khả phản ứng gốc CH3, NH2, OH vào vị trí khác HNCO dựa vào độ mềm. .. cao Để kiểm nghiệm kết luận này, tiến hành xét khả phản ứng gốc vào vị trí khác phân tử HNCO dựa vào lượng phản ứng 3.2 Khả phản ứng dựa vào hàng rào lượng 3.2.1 Tác nhân CH3 Kết tối ưu hóa hình... tích phân tử HNCO khơng định vị mà trải tồn phân tử Đối với hệ khơng có phù hợp tốt hàng rào lượng độ mềm cục bộ, nên dựa vào hàng rào lượng thực nghiệm để dự đoán khả phản ứng tác nhân vào vị trí