Tạp chí Khoa học ĐHQGHN: Khoa học Tự nhiên Công nghệ, Tập 32, Số (2016) 38-46 Khảo sát cấu trúc dẫn xuất Brom phân tử Bisphenol A lý thuyết phiếm hàm mật độ Vũ Văn Đạt1,*, Lê Kim Long2, Nguyễn Hoàng Trang3 Đại học Quốc gia Hà Nội, 144 Xuân Thủy, Hà Nội, Việt Nam Trung tâm Ứng dụng Tin học Hóa học, Khoa Hoá học, Trường Đại học Khoa học Tự nhiên, ĐHQGHN, 19 Lê Thánh Tông, Hà Nội, Việt Nam Trường Đại học Giáo dục, ĐHQGHN, 144 Xuân Thủy, Hà Nội, Việt Nam Nhận ngày 09 tháng 11 năm 2015 Chỉnh sửa ngày 28 tháng 12 năm 2015; Chấp nhận đăng ngày 10 tháng năm 2016 Tóm tắt: Tiếp sau báo “Khảo sát cấu trúc dẫn xuất clo phân tử Bisphenol A lý thuyết phiếm hàm mật độ” [1] chúng tơi tiếp tục trình bày kết nghiên cứu số dẫn xuất chứa brom Bisphenol A pha khí Trong nghiên cứu tiếp tục sử dụng Lý thuyết phiếm hàm mật độ (DFT) với hàm M06 Hybridmeta - GGA sử dụng hàm sở DGAUSSDZVP TZVP kèm theo hàm DGA1 TZVPfit mật độ tương thích Các tính tốn tập trung tìm cách tối ưu hóa cấu trúc phân tử, tần số dao động giá trị lượng obitan với độ xác phù hợp Kết thu phương pháp M06 mơ tả xác cấu hình phân tử dẫn xuất brom bisphenol A Phân bố mật độ điện tích tính cho thấy rõ ràng thay đổi mật độ điện tích phân tử đưa nhóm brom vào phân tử so với chưa thay Các kết thu sử dụng để bàn luận chế phản ứng thay brom vào phân tử bisphenol A giải thích kết thu từ nghiên cứu thực nghiệm Từ khóa: Bisphenol A, Lý thuyết phiếm hàm mật độ (DFT), hàm M06 Hybridmeta – GGA BPA thường suốt, cứng sử dụng nhiều sản xuất hàng hóa tiêu dùng thơng thường bình sữa cho trẻ em, thiết bị thể thao, đĩa CD DVD, v.v…[1,2] Tuy nhiên thời gian gần có nhiều nghiên cứu khuyến cáo tác hại BPA thể người [3,4] Một số nghiên cứu BPA chất gây ức chế nội tiết hay đóng vai trị chất giả hocmon dẫn đến rối loạn nội tiết thể, góp phần gia tăng rủi ro bệnh liên quan đến tim mạch, béo phì [5], ảnh hưởng khả phát triển trí não trẻ em [6]; chí làm giảm độ nhạy cảm Đặt vấn đề∗ Bisphenol A (BPA) hợp chất hữu chứa hai nhóm phenol có cơng thức (CH3)[4,4'-dihydroxy-2,22C(C6H4OH)2 diphenylpropane] sử dụng rộng rãi sản xuất polycarbonate nhựa epoxy BPA sử dụng chất làm ổn định hay chống oxi hóa cho nhiều loại chất dẻo (plastics) PVC, v.v Các sản phẩm nhựa làm từ _ ∗ Tác giả liên hệ ĐT: 84-934277732 Email: datvv@vnu.edu.vn 38 V.V Đạt nnk / Tạp chí Khoa học ĐHQGHN: Khoa học Tự nhiên Cơng nghệ, Tập 32, Số (2016) 38-46 bệnh nhân sử dụng biện pháp hóa trị liệu, ảnh hưởng hoạt động tuyến tiền liệt, gây ung thư vú, u nang buồng trứng …[7] BPA tự thơi nhiễm thực phẩm từ vật dụng làm nhựa có nguồn gốc từ BPA tác động nhiệt độ, thời gian [4] Khi tan nước, BPA dễ dàng phản ứng với brom tự nước (tồn dạng hipobromua) tạo thành dẫn xuất brom BPA [8,9] Chính có nhiều tác dụng phụ sản phẩm sử dụng BPA mà nhiều nghiên cứu gần tìm chất thay BPA hãng sản xuất đồ đựng phải cơng bố sản phẩm khơng có BPA tự Trong khuôn khổ báo, tiến hành khảo sát cấu trúc phân tử dẫn xuất brom BPA tính tốn lí thuyết hóa 2-bromo-4-[2-(4-hydroxyphenol)propan-2-yl]phenol Phân tử chứa Br (I) 2,6-dibromo-4-[2-(3-bromo-4hydroxyphenyl)propan-2-yl]phenol Phân tử chứa Br (III) 39 lượng tử Việc xem xét cấu trúc BPA mạng tinh thể [10] giúp xây dựng cấu trúc ban đầu cho phân tử BPA dạng tự Q trình tính tốn lượng tử thực theo bước: tối ưu hóa cấu hình xác định tần số dao động ứng với cấu hình tối ưu xác định Tối ưu hóa cấu hình nhằm xác định cấu hình ứng với trạng thái cực tiểu lượng phân tử để thu thông số cấu trúc, lượng obitan phân tử mật độ điện tích, lượng obitan phân tử HOMO- LUMO, momen lưỡng cực, lượng toàn phần Mục đích nghiên cứu tiền đề để tiếp tục thực nghiên cứu sâu để xây dựng mối liên hệ định lượng cấu trúc phân tử mơ hoạt tính sinh học BPA dẫn xuất 2,6-dibromo-4-[2-(4-hydroxyphenyl)propan-2yl]phenol Phân tử chứa Br (II) 2,6-dibromo-4-[2-(3,5-dibromo-4hydroxyphenyl)propan-2-yl]phenol Phân tử chứa Br (IV) Hình Các dẫn xuất chứa BPA Phương pháp nghiên cứu Các tính tốn thực dựa lý thuyết phiếm hàm mật độ (Density Functional Theory- DFT) Việc tối ưu hóa cấu hình xác định cực tiểu toàn cục thực phiếm hàm lai hóa Meta GGA – phiếm hàm M06 [11] Bộ hàm sở DGAUSS 40 V.V Đạt nnk / Tạp chí Khoa học ĐHQGHN: Khoa học Tự nhiên Công nghệ, Tập 32, Số (2016) 38-46 hiệu chỉnh DZVP (Double-Zeta Valance Polarization) [12], TZVP (Triple-Zeta Valance Polarization) [13] sử dụng nhằm làm giảm thời gian tính tốn mà cho kết tối ưu Các hàm sở DZVP TZVP hiệu chỉnh cách đưa thêm vào hàm gốc tệp mật độ sở mô DGA1 [11] (với hàm sở DZVP) TZVPfit [14] (với hàm sở TZVP) tích hợp phiên Gaussian 09 Tất cấu hình phân tử tối ưu ứng với trạng thái cực tiểu tồn cục xác định việc phân tích tần số dao động điều hòa tiêu chuẩn hội tụ tích hợp sẵn phần mềm Gaussian 09 [15] Kết thảo luận Hình Phân tử monobrom bisphenol A (I) thông số cấu trúc ứng với cấu hình tối ưu phân tử khảo sát (độ dài liên kết, góc liên kết) trình bày Bảng Các giá trị so sánh với thông số cấu trúc phân tử BPA phân tử dẫn xuất BPA chứa Clo tính tốn phương pháp [16] Theo kết thu phương pháp M06/TZVP phân tử BPA dẫn xuất brom có cấu trúc khơng gian hình "cánh quạt" với hai vịng benzen phân bố gần so le (Hình 1) Các góc xoắn hai mặt phẳng vịng benzen là: C15 C10 C2 C4 = 131.1° C5 C4C2 C10 = 133.8° phân tử monobrom bisphenol A (Hình 1) Các Bảng Một số thơng số lựa chọn từ cấu hình tối ưu BPA dẫn xuất trình brom hóa với BPA M06/TZVP Độ dài liên kết/Góc liên kết BPA [16] r(Cphenyl – Cphenyl)* r(Cmetyl – Cphenyl )* r(Cmetyl– Cmetyl)* r(Cphenyl– O)* r(Cmetyl – H)* r(Cphenyl – H)* ∠Cphenyl – Cphenyl – O 1.394 1.510 1.528 1.376 1.116 1.098 116.6 I II III IV 1.395 1.528 1.532 1.363 1.096 1.085 117.7 1.393 1.527 1.534 1.351 1.096 1.086 117.9 1.394 1.528 1.533 1.351 1.096 1.089 119.0 1.394 1.528 1.533 1.345 1.095 1.089 119.3 Phân tử chứa clo [16] 1.394 1.527 1.534 1.358 1.095 1.088 118.1 V.V Đạt nnk / Tạp chí Khoa học ĐHQGHN: Khoa học Tự nhiên Công nghệ, Tập 32, Số (2016) 38-46 ∠Cphenyl – O – H 107.8 r(Cphenyl – Br(1) ) r(Cphenyl – Br(2)) r(Cphenyl – Br(3)) r(Cphenyl – Br(4)) ∠Cphenyl – Cphenyl – Br(1) ∠Cphenyl – Cphenyl – Br(2) ∠Cphenyl – Cphenyl – Br(3) ∠Cphenyl – Cphenyl – Br(4) 108.7 1.725 — — — 119.0 — — — 108.9 108.8 108.8 1.908 1.908 — — 119.5 118.8 — — 1.895 1.894 1.908 — 119.5 118.7 118.2 — 1.894 1.906 1.893 1.906 119.5 118.2 118.9 118.2 41 108.9 1.756 — — — 118.5 — — — * Giá trị trung bình liên kết (Å) Từ bảng nhận thấy khuynh hướng thay đổi thông số cấu trúc so với phân tử BPA dẫn xuất brom tương tự với dẫn xuất clo [16] Độ dài liên kết r(Cphenyl – Cphenyl), r(Cphenyl – H) vịng benzen thay đổi khơng đáng kể nguyên tử brom vào vòng benzen phân tử BPA chứng tỏ vịng benzen có cấu trúc ổn định khơng bị ảnh hưởng có mặt nhóm khác Đối với nhóm CH3 q trình tối ưu hóa cấu hình giữ ngun đối xứng không gian C3v Các thông số cấu trúc Cmetyl – Cphenyl, Cmetyl – H thay đổi không đáng kể không phụ thuộc vào số lượng phân tử brom Liên kết Cphenyl – O giảm theo chiều hướng tăng nguyên tử brom phân tử khảo sát (giảm 0.031 Å phân tử chứa brom (IV)) Sự thay đổi mật độ phân bố electron phân tử tăng xuất đẩy lớp vỏ electron nguyên tử brom (4s23d104p5) với lớp vỏ electron nguyên tử oxy (2s22p4) làm cho liên kết C - O phân tử khảo sát bị "ngắn" lại, góc liên kết Cphenyl – Cphenyl – O mà thay đổi tăng dần theo số lượng brom vòng benzen (1.100 – 2.710) Cấu trúc phân tử BPA xác định phương pháp M06/TZVP gần với cấu trúc BPA mạng tinh thể xác định phương pháp X-ray nhiệt độ 115 K [10] Điểm khác biệt góc xoắn hai mặt phẳng vịng benzen Đối với BPA tồn trạng thái rắn [10], phân tử chịu tương tác phân tử khác mạng tinh thể nên cấu trúc không gian bị thu hẹp lại dẫn đến hai góc xoắn C15C10C2 C4 C5 C4C2 C10 có giá trị nhỏ so với phân tử tự xác định phương pháp M06/TZVP Để so sánh với phân tử dẫn xuất clo, Bảng có đưa vào thông số cấu trúc phân tử chứa clo Nhận thấy thông số cấu trúc dẫn xuất brom gần với dẫn xuất clo (Bảng 1) Các khuynh hướng thay đổi cấu trúc dẫn xuất brom dẫn xuất clo giống Như từ phân tích kết thu từ tính tốn lý thuyết nói cấu trúc mạch BPA gần với cấu trúc mạch dẫn xuất brom clo hay nói cách khác việc nguyên tử H ngun tử clo, brom vào vịng phenyl khơng ảnh hưởng đến cấu trúc khơng gian BPA nói chung Cấu trúc dẫn xuất clo brom gần Điều quan trọng tương đồng cấu trúc định phần lớn tương đồng tính chất hóa học tính chất vật lý chất Đây tiền đề để thảo luận tiếp chế phản ứng tạo dẫn xuất clo brom dung dịch báo Trong khuôn khổ báo này, đánh giá thay đổi hai mức lượng obitan phân tử HOMO LUMO (Bảng 2) Có thể thấy mức lượng HOMO LUMO dẫn xuất brom giảm so với phân tử BPA số lượng brom vòng benzen tăng lên (Hình 2) Khuynh hướng giống với trường hợp dẫn xuất clo [16] Khi mức lượng HOMO giảm giá trị MIP (Molecular Ionization Potental) tăng Như xuất 42 V.V Đạt nnk / Tạp chí Khoa học ĐHQGHN: Khoa học Tự nhiên Công nghệ, Tập 32, Số (2016) 38-46 nguyên tử nhóm halogen (clo, brom) làm tăng khả ion hóa dẫn xuất halogen BPA Bảng Các thông số lượng tử thu phương pháp M06/TZVP Phân tử EHOMO (eV) ELUMO (eV) BPA I II III IV -5.92 -6.33 -6.49 -6.58 -6.69 -0.51 -0.62 -0.79 -0.94 -1.12 µ (Debye) 2.42 2.57 1.22 3.66 0.36 Theo giá trị tính tốn thu phương pháp M06/TZVP lượng phân tử dẫn xuất brom nhỏ Nhận thấy lượng phân tử giảm số lượng nguyên tử brom có phân tử khảo sát lớn Như chứng tỏ độ bền phân tử khảo sát tăng số lượng nguyên tử brom phân tử tăng lên dẫn xuất phân tử chứa brom bền so với BPA dẫn xuất lại Giá trị mô men lưỡng cực khác biệt rõ với giá trị cao đạt phân tử chứa brom cho thấy cấu trúc bất đối xứng phân tử (Bảng 2) HOMO LUMO (I) (II) (III) (IV) Hình Hình dạng vị trí định xứ HOMO LUMO dẫn xuất BPA Hình trình bày hình dạng vị trí định xứ obitan phân tử HOMO LUMO dẫn xuất brom Các obitan phân tử HOMO có khuynh hướng định chỗ xung quanh nguyên tử brom, nơi có điện tích riêng phần lớn Trong đó, obitan phân tử LUMO lại có khuynh hướng phân bố hướng xung quanh nhóm nguyên tử –OH, hay – CH3, khu vực có điện tích riêng phần nhỏ Sự phân bố ảnh hưởng đến khả phản ứng đồng phân tác nhân cơng có định hướng khác đến tâm phản ứng phân tử Sự khác biệt có nguyên nhân từ thay đổi mật độ điện tích phân tử nhóm brom đưa vào phân tử BPA Bằng phương pháp M06/TZVP thu phân bố điện tích nguyên tử cacbon vòng benzen BPA dẫn xuất halogen (Bảng 3) Có thể thấy xuất nguyên tử nhóm halogen làm cho electron vòng benzen bị phân bố lại dẫn đến thay đổi mật độ điện tích nguyên tử cacbon vịng Thơng thường phản ứng hợp chất chứa nhóm OH với tác nhân hipobromơ (HOBr), HOBr phản ứng theo chế electrophin vào vòng thơm Kết phân tích mật độ điện tích cho thấy phân tử BPA hai vị trí C3 C21 V.V Đạt nnk / Tạp chí Khoa học ĐHQGHN: Khoa học Tự nhiên Công nghệ, Tập 32, Số (2016) 38-46 mang điện tích riêng phần nhỏ 0.216 -0.211 tham gia phản ứng brom hóa hai vị trí ưu tiên brom so với vị trí khác, khả brom hai vị trí tương đương chênh 43 lệch điện tích riêng phần C3 C21 nhỏ Trong phần tính tốn lượng tử chúng tơi giả định brom vào vị trí C3 để tạo thành monobrom bisphenol A (I) (Hình 4) Bảng Phân bố điện tích ngun tử cacbon vịng benzen BPA dẫn xuất halogen M06/TZVP C2 C3 C4 C5 C6 C17 C18 C19 C20 C21 -0.216 -0.091 C3 HO BPA I II III IV -0.091 -0.216 -0.155 -0.088 -0.083 -0.096 -0.159 -0.211 -0.220 -0.161 -0.502 -0.380 -0.317 -0.323 -0.238 -0.380 -0.507 -0.455 -0.041 -0.156 -0.159 -0.318 -0.105 -0.010 -0.142 -0.144 -0.155 -0.148 -0.091 -0.140 -0.141 -0.163 -0.055 -0.022 -0.0961 -0.080 -0.211 -0.148 -0.045 -0.120 -0.202 -0.110 -0.066 -0.045 -0.066 -0.110 -0.202 -0.120 -0.096 C18 C2 C5 -0.155 C20 C 17 C6 -0.088 -0.083 Br -0.159 C19 C1 C4 Phân tử chứa clo [16] -0.221 -0.095 -0.272 -0.235 -0.183 -0.179 -0.194 -0.254 -0.294 -0.281 -0.220 C3 OH C21 -0.211 HO C2 -0.380 -0.238 C18 C1 C4 C5 -0.380 C6 -0.317 C19 C20 C 17 -0.323 OH C 21 -0.455 4,4 – propane-2,2-diyl)diphenol Phân tử BPA 3-chloro-4-[2-(4-hydroxyphenol)propan-2-yl]phenol Phân tử chứa Clo Hình Mật độ điện tích nguyên tử C vòng phenyl dẫn xuất brom BPA Các phản ứng nguyên tử H nguyên tử brom tạo dibrom bisphenol A (II), tribrom bisphenol A (III) theo chế electronphin tiếp tục xảy vị trí mà nguyên tử cacbon 44 V.V Đạt nnk / Tạp chí Khoa học ĐHQGHN: Khoa học Tự nhiên Công nghệ, Tập 32, Số (2016) 38-46 mang điện tích riêng phần âm Phân tử chứa nguyên tử brom tiếp tục phản ứng với tác nhân hipobromơ, phản ứng electrophin nguyên tử H nguyên tử brom nguyên tử C19 C5 hình thành phân tử tribrom bisphenol A (Phân tử chứa nguyên tử brom) Thực nghiệm [9] cho thấy khả hình thành phân tử tribrom bisphenol A hai vị trí Trong cấu hình phân tử chúng tơi thu được, điện tích riêng phần vị trí nguyên tử C19 nhỏ điện tích riêng phần nguyên tử C5 thế, chúng tơi dự đốn cấu hình sản phẩm phản ứng brom hóa phân tử dibrom bisphenol A sản phẩm phản ứng vị trí C19 dẫn xuất brom cuối q trình brom hóa bisphenol A phân tử tetrabrom bisphenol A (Phân tử chứa nguyên tử brom) So sánh phân bố điện tích vịng benzen monobrom bisphenol A (I) với monoclo bisphenol A thấy nguyên tử cacbon vòng benzen dẫn xuất brom mang điện tích riêng phần âm so với dẫn xuất clo (Bảng 3) Nguyên nhân brom có điện tích âm so với clo nên đẩy electron vòng benzen mạnh làm cho mật độ điện tích phân bố vòng benzen trường hợp dẫn xuất brom lớn so với trường hợp dẫn xuất clo BPA Phân tích phổ IR BPA, tetrabrom bisphenol A, tetraclo bisphenol A thu tính tốn thực nghiệm không quan sát thấy vạch phổ đặc trưng cho nhóm halogen Phổ thu tính toán lý thuyết phù hợp với phổ đo thực nghiệm [17] Sự dịch chuyển vạch phổ thực nghiệm so với tính tốn lý thuyết có tương tác phân tử ảnh hưởng sai số máy đo Phổ BPA dẫn xuất halogen giống Dưới bảng phân tích phổ BPA tetrabrom bisphenol A (IV) thu tính tốn lý thuyết phương pháp M06/DGDZVP Bảng Tần số dao động BPA tetrabrom bisphenol A (IV) (M06/DGDZVP) Tần số dao động (cm-1) OH Dao động vòng brom-benzen Dao động vòng brom-benzen BPA - Kết luận Trong khuôn khổ báo, tập trung nghiên cứu cấu trúc phân tử phân tích thơng số lượng tử BPA bốn dẫn xuất brom bisphenol A thu sau tối ưu cấu hình Phương pháp tính tốn sử dụng M06 với hai hàm sở DGAUSS hiệu chỉnh DZVP TZVP với có mặt tệp mật độ mô DGA1 TZVPfit Các thông số lượng tử cấu hình tối ưu BPA bốn dẫn xuất brom bisphenol A bao gồm: thông số cấu trúc phân tử, lượng obitan phân tử HOMO- LUMO, mô men lưỡng cực, lượng tồn phần xét đến Các thơng số cấu trúc phân tử bốn phân tử phù hợp với số liệu thực nghiệm IV 3733 1660.2-1618.0 1431-1369 BPA [17] 3477 1553 1472 Việc phân tích thơng số lượng tử tạo điều kiện để thực nghiên cứu sâu tính chất hóa học, sinh học dẫn xuất BPA Tài liệu tham khảo [1] M.K Morgan, P.A Jones, A.M Calafat, X Ye, C.W Croghan, J.C Chuang, N.K Wilson, M.S Clifton, Z Figueroa, L.S Sheldon, Environ Sci Technol 2011, 45, 5309 [2] A Schecter, N Malik, D Haffner, S Smith, T.R Harris, O Paepke, L Birnbaum, Environ Sci Technol 2010, 44, 9425 [3] L N Vandenberg, R Hauser, M Marcus, N Olea, W V Welshons Reproductive Toxicology 24 (2007) 139–177 V.V Đạt nnk / Tạp chí Khoa học ĐHQGHN: Khoa học Tự nhiên Công nghệ, Tập 32, Số (2016) 38-46 [4] A V Krishnan, P Stathis, S F Permuth, L Tokes, D Feldman The Endocrine Society 1993, Vol 132, No [5] Rezg R, El-Fazaa S, Gharbi N, Mornagui B (March 2014) "Bisphenol A and human chronic diseases: Current evidences, possible mechanisms, and future perspectives" Environment International 2014, 64, 83–90 [6] D.R Doerge, N.C Twaddle, M Vanlandingham, R.P Brown, J.W Fisher, Toxicol Appl Pharmacol 2011, 255, 261 [7] Ho SM, Tang WY, Belmonte de Frausto J, Prins GS (2006) "Developmental exposure to estradiol and bisphenol A increases susceptibility to prostate carcinogenesis and epigenetically regulates phosphodiesterase type variant 4" Cancer Res 66 (11): 5624–32 [8] Fukazawa H, Hoshino K, Shiozawa T, Matsushita H, Terao Y (2001) "Identification and quantification of chlorinated bisphenol a in wastewater from wastepaper recycling plants".Chemosphere 44 (5) 973–9 [9] Hu JY; Aizawa T.; Ookubo S., Products of Aqueous Chlorination of Bisphenol A and Their Estrogenic Activity Environmental Science & Technology 2002,36 (9), 1980-1987 [10] C F Lim, J M Tanski J Chem Crystallogr (2007) 37:587–595 [11] Zhao, Y.; Truhlar, D., The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic [12] [13] [14] [15] [16] [17] 45 testing of four M06-class functionals and 12 other functionals Theor Chem Account 2008,120 (1-3), 215-241 Godbout, N.; Salahub, D R.; Andzelm, J.; Wimmer, E., Optimization of Gaussian-type basis sets for local spin density functional calculations Part I Boron through neon, optimization technique and validation Canadian Journal of Chemistry 1992,70 (2), 560-571 Schafer, A.; Huber, C.; Ahlrichs, R., Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr The Journal of Chemical Physics 1994,100 (8), 58295835 (a) Eichkorn, K.; Treutler, O.; Öhm, H.; Häser, M.; Ahlrichs, R., Auxiliary basis sets to approximate Coulomb potentials Chemical Physics Letters 1995,240 (4), 283-290;(b) Eichkorn, K.; Weigend, F.; Treutler, O.; Ahlrichs, R., Auxiliary basis sets for main row atoms and transition metals and their use to approximate Coulomb potentials Theor Chem Account 1997,97 (1-4), 119-124 Frisch, M.J.T., G.W et al , Gaussian 09, Revision D.01 Gaussian, Inc., Wallingford CT, 2009 V.V Đạt, L.K Long Tạp chí Khoa học Đại học Quốc gia Hà Nội (Khoa học tự nhiên công nghệ) Tập 31, Số 3, 2015 M P Luda, A I Balabanovich1y, A Hornung, G Camino Polym Adv Technol 14, 741–748 (2003) Study on the Structure Brome Derivatives of Bisphenol a: Density Functional Theory Approach Vũ Văn Đạt1, Lê Kim Long2, Nguyễn Hoàng Trang3 Vietnam National University, Hanoi, 144 Xuân Thủy, Hanoi, Vietnam Center for Applied Informatics in Chemistry, Faculty of Chemistry, VNU University of Science, 19 Lê Thánh Tơng, Hồn Kiếm, Hanoi, Vietnam VNU University of Education, 144 Xuân Thủy, Hanoi, Vietnam Abstract: In this article it is presented a theoretical investigation of some brome derivatives of Bisphenol A in gas phase The Density Functional Theory (DFT) is utilized, using the M06 46 V.V Đạt nnk / Tạp chí Khoa học ĐHQGHN: Khoa học Tự nhiên Công nghệ, Tập 32, Số (2016) 38-46 hybridmeta - GGA functionals with the DGAUSSDZVP and TZVP basis sets accompanying with DGA1 and TZVPfit density fitting sets These calculations focuse on finding the optimized molecular structures, vibrational frequencies, the molecular orbital energies with reasonable accuracy Results obtained have shown that M06 method describes precisely configurations of all chlorine derivatives of bisphenol A When BPA has been chlorinated there are some changes in electric density of molecules Calculation results have been agreed well to those of experiments Keywords: Bisphenol A, Density Functional Theory, M06 hybridmeta - GGA functional ... phân tử dẫn xuất brom nhỏ Nhận thấy lượng phân tử giảm số lượng nguyên tử brom có phân tử khảo sát lớn Như chứng tỏ độ bền phân tử khảo sát tăng số lượng nguyên tử brom phân tử tăng lên dẫn xuất. .. nguyên tử C5 thế, dự đốn cấu hình sản phẩm phản ứng brom h? ?a phân tử dibrom bisphenol A sản phẩm phản ứng vị trí C19 dẫn xuất brom cuối q trình brom h? ?a bisphenol A phân tử tetrabrom bisphenol A (Phân. .. thông số cấu trúc phân tử BPA phân tử dẫn xuất BPA ch? ?a Clo tính tốn phương pháp [16] Theo kết thu phương pháp M06/TZVP phân tử BPA dẫn xuất brom có cấu trúc khơng gian hình "cánh quạt" với hai vịng