(Luận văn) cấu trúc và năng lượng của các trạng thái spin của cluster fegen −0 (n=1 2)

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC ĐỒNG THÁP - - NGUYỄN THỊ LAN PHƯƠNG an lu n va gh tn to LUẬN VĂN THẠC SĨ HÓA HỌC p ie d oa nl w oi m ll fu an v an lu CẤU TRÚC VÀ NĂNG LƯỢNG CỦA CÁC TRẠNG THÁI SPIN CỦA CLUSTER FeGen−/0 (n=1-2) nh at z z @ om l.c gm an Lu ĐỒNG THÁP – NĂM 2019 n va a th c si BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC ĐỒNG THÁP - - LUẬN VĂN THẠC SĨ HĨA HỌC an lu n va Chun ngành: Hóa lý thuyết Hóa lý Mã số: 8.44.01.19 p ie gh tn to d oa nl w CẤU TRÚC VÀ NĂNG LƯỢNG CỦA CÁC TRẠNG THÁI SPIN CỦA CLUSTER FeGen−/0 (n=1-2) fu an v an lu Người hướng dẫn khoa học: Nguyễn Thị Lan Phương TS Trần Văn Tân oi m ll Người thực hiện: nh at z z @ om l.c gm an Lu ĐỒNG THÁP – NĂM 2019 n va a th c si i LỜI CAM ĐOAN Tơi xin cam đoan luận văn hồn thành dựa kết nghiên cứu tính tốn tơi Những số liệu, kết chưa dùng cho luận văn khác Đồng Tháp, tháng 07 năm 2019 Người thực luân văn an lu n va p ie gh tn to Nguyễn Thị Lan Phương d oa nl w oi m ll fu an v an lu nh at z z @ om l.c gm an Lu n va a th c si ii LỜI CẢM ƠN Nhà khoa học Isaac Newton nói: “Điều ta biết giọt nước, điều ta chưa biết đại dương” Biển tri thức mênh mông vô tận mà hiểu biết giọt nước bé nhỏ Để tiến xa hành trình tích lũy tri thức, người phải không ngừng sức học tập từ thầy cô, bạn bè,… Với tôi, quãng thời gian học tập, nghiên cứu trường Đại học Đồng Tháp quãng thời gian có ý nghĩa thiêng liêng, cho tơi tiếp cận với tri thức lớn lao nhân loại Và đến ngày hôm nay, luận văn tơi lu an hồn thành niềm vui sướng, mong đợi tin yêu quý thầy cô va Lời đầu tiên, xin bày tỏ lòng biết ơn sâu sắc đến TS Trần Văn Tân - người n gh tn to thầy tận tâm truyền đạt cho tri thức vô giá, kỹ quý báu Bằng tất lịng nhiệt huyết, thầy tận tình hướng dẫn, giúp đỡ hỗ trợ suốt thời gian làm luận văn p ie w Tôi xin gửi lời cảm ơn sâu sắc đến Lãnh đạo Trường Đại học Đồng Tháp, PGS.TS Trần Quốc Trị tất cán Khoa Sư phạm Lý-Hóa-Sinh, người giúp đỡ d oa nl tạo điều kiện thuận lợi cho tơi suốt q trình học tập thực luận văn Và xin chân thành gửi lời tri ân đến tất quý thầy, cô trường Đại học Đồng v an lu Tháp, người tận tình giảng dạy truyền đạt cho tơi kiến thức quý giá fu an suốt thời gian học tập trường Cuối cùng, xin gửi lời biết ơn sâu sắc đến người thân gia đình, bạn bè m ll bạn học viên Cao học Hóa lý thuyết Hóa lý khóa 6-Trường Đại học Đồng oi nh Tháp động viên, giúp đỡ đồng hành với suốt q trình làm luận văn at Tơi xin chúc tất quý thầy cô an khang, ln có nhiều niềm vui, niềm z z hạnh phúc sống gặt hái nhiều thành công nghiệp giáo dục @ gm Đồng Tháp, tháng 07 năm 2019 om l.c Người thực luận văn an Lu n va a th c si iii Nguyễn Thị Lan Phương MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC .iii DANH MỤC CHỮ VIẾT TẮT v an lu DANH MỤC BẢNG .vii n va DANH MỤC HÌNH .viii gh tn to MỞ ĐẦU 1 Lí chọn đề tài p ie Mục tiêu đề tài w Đối tượng nghiên cứu phạm vi nghiên cứu d oa nl Nội dung luận văn Phương pháp nghiên cứu v an lu NỘI DUNG fu an Chương TỔNG QUAN m ll 1.1 Giới thiệu cluster germanium .5 oi nh 1.2 Cluster germanium pha tạp kim loại chuyển tiếp at 1.3 Cluster germanium pha tạp sắt 10 z z @ 1.4 Cơ sở lý thuyết hóa học tính toán 12 gm 1.4.1 Phương trình Schrưdinger .12 om l.c 1.4.2 Sự gần Born – Oppenheimer 13 1.4.3 Phương pháp Hatree – Fock (HF) 14 an Lu n va a th c si C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an iv 1.4.4 Các phương pháp tính đa cấu hình CASSCF/CASPT2 RASSCF/RASPT2 16 1.4.5 Lý thuyết phiếm hàm mật độ 18 1.4.6 Phiếm hàm trao đổi tương quan 21 Chương PHƯƠNG PHÁP NGHIÊN CỨU 22 2.1 Phương pháp tính đa cấu hình CASSCF/CASPT2 RASSCF/RASPT2 22 an lu 2.2 Lý thuyết phiếm hàm mật độ .23 2.3 Các phần mềm tính tốn hóa học lượng tử 24 va n 2.3.1 Phần mềm ORCA 4.0 24 gh tn to 2.3.2 Phần mềm MOLCAS 24 p ie 2.3.3 Phần mềm GABEDIT 24 w Chương KẾT QUẢ VÀ THẢO LUẬN 25 d oa nl 3.1 Kết nghiên cứu cluster FeGe–/0 25 3.1.1 Cấu trúc hình học độ bền cluster FeGe–/0 25 v an lu 3.1.2 Cấu trúc electron trạng thái cluster FeGe–/0 32 fu an 3.2 Kết nghiên cứu đồng phân cluster FeGe2–/0 34 m ll 3.2.1 Cấu trúc độ bền đồng phân cluster FeGe2–/0 34 oi nh 3.2.2 Tần số dao động điều hoà trạng thái spin cluster FeGe2–/0 40 at z 3.2.3 Cấu trúc electron trạng thái spin cluster FeGe2–/0 .44 z @ KẾT LUẬN VÀ KIẾN NGHỊ 49 gm Kết luận 49 l.c Kiến nghị 50 om Lu TÀI LIỆU THAM KHẢO 51 an PHỤ LỤC n va a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an v DANH MỤC CHỮ VIẾT TẮT an lu (Adiabatic Detachment Energy): lượng tách electron ADE B3LYP phiếm hàm trao đổi- tương quan lai ghép B3LYP B3P86 phiếm hàm trao đổi- tương quan lai ghép B3P86 B3PW91 phiếm hàm trao đổi Becke phiếm hàm tương quan PW91 BHLYP phiếm hàm trao đổi- tương quan lai ghép BHLYP BLYP phiếm hàm trao đổi Becke phiếm hàm tương quan LYP n va ADE gh tn to phiếm hàm trao đổi Becke phiếm hàm tương quan tạo BP86 p ie Perdew năm 1986 w CASSCF (Complete Active Space Self-Consistent Field): phương pháp trường d oa nl tự hợp không gian hoạt động đầy đủ (Complete Active Space Second-Order Perturbation Theory): phương v an lu CASPT2 pháp nhiễu loạn bậc hai không gian hoạt động đầy đủ fu an CCSD (Coupled-Cluster with Single and Double excitations): phương pháp CCSD(T) oi m ll coupled-cluster với kích thích đơn, kích thích đơi (Coupled-Cluster with Single and Double and pertubative Triple nh at excitations): phương pháp coupled-cluster với kích thích đơn, kích z thích đơi, nhiễu loạn kích thích ba z @ (Density Functional Theory): lý thuyết phiếm hàm mật độ GGA (Generalized Gradient Approximation) (sự gần gradient suy rộng): phiếm hàm trao đổi-tương quan GGA om l.c gm DFT HSE06 phiếm hàm trao đổi- tương quan lai ghép HSE06 an (Hartree-Fock): Phương pháp Hartree-Fock Lu HF n va a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an vi (Local Spin-Density Approximation) (sự gần mật độ spin địa LSDA phương): phiếm hàm trao đổi-tương quan LSDA NEVPT2 (n-Electron Valence State Second-Order Perturbation Theory): phương pháp nhiễu loạn bậc hai n electron hóa trị RASSCF (Restricted Active Space Self-Consistent Field): phương pháp trường tự hợp không gian hoạt động giới hạn RASPT2 (Restricted Active Space Second-Order Perturbation Theory): phương lu an pháp nhiễu loạn bậc hai không gian hoạt động giới hạn va (Vertical Detachment Energy): lượng tách electron VDE n VDE p ie gh tn to d oa nl w oi m ll fu an v an lu nh at z z @ om l.c gm an Lu n va a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an vii DANH MỤC BẢNG Bảng 3.1 Các thơng số cấu trúc hình học, cấu hình electron lượng tương đối tính cluster anion FeGe– 28 Bảng 3.2 Các thông số cấu trúc hình học, cấu hình electron lượng tương an lu đối tính cluster trung hịa điện FeGe 30 Bảng 3.3 Các thơng số cấu trúc hình học tính theo BP86, cấu hình electron va n lượng tương đối tính theo RASPT2 đồng phân cluster anion FeGe2– 34 gh tn to Bảng 3.4 Các thông số cấu trúc hình học tính theo BP86, cấu hình electron lượng tương đối tính theo RASPT2 cluster trung hòa điện FeGe2 35 p ie w Bảng 3.5 Năng lượng tương đối lượng tách ADE cluster anion FeGe2– tính theo lý thuyết phiếm hàm mật độ phương pháp RASPT2 37 d oa nl Bảng 3.6 Năng lượng tương đối lượng tách ADE cluster trung hịa v an lu FeGe2 tính theo lý thuyết phiếm hàm mật độ phương pháp RASPT2 38 Bảng 3.7 Tần số dao động điều hòa trạng thái electron thuộc cluster fu an FeGe2−/0 tính phiếm hàm BP86 41 oi m ll nh at z z @ om l.c gm an Lu n va a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an viii DANH MỤC HÌNH Hình 1.1 Hình ảnh minh họa cho hàm sóng CASSCF 16 Hình 3.1 Cấu trúc hình học đồng phân cluster FeGe–/0 FeGe2–/0 25 Hình 3.2 Đường cong trạng thái cluster anion FeGe– 26 lu an Hình 3.3 Đường cong trạng thái cluster FeGe 27 n va Hình 3.4 Hình ảnh orbital phân tử số electron chiếm trạng thái gh tn to A2 (2Δ) đồng phân thẳng FeGe– thu từ phép tính CASSCF 32 Hình 3.5 Các kiểu dao động điều hòa trạng thái 4B1 đồng phân vòng p ie cluster anion FeGe2– tính phiếm hàm BP86 42 w Hình 3.6 Các kiểu dao động điều hòa trạng thái 5B2 đồng phân vòng d oa nl cluster trung hịa điện FeGe2 tính phiếm hàm BP86 42 Hình 3.7 Các kiểu dao động điều hòa trạng thái 4A2 đồng phân thẳng v an lu cluster FeGe2– tính phiếm hàm BP86 43 fu an Hình 3.8 Các kiểu dao động điều hòa trạng thái 5A2 (5Δ) đồng phân thẳng m ll cluster FeGe2 tính phiếm hàm BP86 44 oi Hình 3.9 Hình ảnh orbital phân tử số electron chiếm trạng thái nh at B1 đồng phân hình vịng FeGe2– thu từ phép tính RASSCF 45 z Hình 3.10 Hình ảnh orbital phân tử số electron chiếm trạng thái z A2 (4Δ) đồng phân thẳng FeGe2– thu từ phép tính RASSCF 47 @ om l.c gm an Lu n va a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an P3 SYMMETRY ADAPTED BASIS -The coefficients for the symmetry adapted linear combinations (SALCS) of basis functions will now be computed: Number of basis functions 160 Preparing memory done Constructing Gamma(red) done Reducing Gamma(red) done Constructing SALCs done Checking SALC integrity nothing suspicious Normalizing SALCs done Storing the symmetry object: Symmetry file Writing symmetry information c4b1.sym.tmp done an lu - Cấu trúc ban đầu sau thực phép tính FeGe2− (dịng 248 – 256) n va p ie gh tn to ************************************************************* * GEOMETRY OPTIMIZATION CYCLE * ************************************************************* CARTESIAN COORDINATES (ANGSTROEM) Fe 0.000000 0.000000 1.434599 Ge 0.000000 1.176009 -0.551883 Ge 0.000000 -1.176009 -0.551883 + Số vòng lặp SCF: 50 vòng (dòng 462- 527) w -5418.5322709929 -0.000000001529 0.00000695 0.00000012 -5418.5322709930 -0.000000000094 0.00000360 0.00000006 -5418.5322709928 0.000000000160 0.00000191 0.00000003 **** Energy Check signals convergence **** 0.0000025 0.0000 0.0000015 0.0000 0.0000029 0.0000 d oa nl 47 48 49 fu an v an lu ***************************************************** * SUCCESS * * SCF CONVERGED AFTER 50 CYCLES * ***************************************************** + Năng lượng cấu trúc ban đầu tối ưu hóa (dòng 568) : m ll Total Energy -5418.53243224 Eh -147445.76348 eV oi + Năng lượng điểm đơn cấu trúc ban đầu tối ưu hóa (dịng 10117- 10119) nh 5418.532432236496 at z FINAL SINGLE POINT ENERGY - z @ - Cấu trúc hình học FeGe2− sau thực phép tính vịng (dịng 10233 – gm 10241) om l.c ************************************************************* * GEOMETRY OPTIMIZATION CYCLE * ************************************************************* CARTESIAN COORDINATES (ANGSTROEM) Fe 0.000000 -0.000000 1.432297 Ge 0.000000 1.180538 -0.550732 Ge 0.000000 -1.180538 -0.550732 an Lu n va a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an P4 + Số vòng lặp SCF: 10 vòng (dòng 462- 527) -5418.5325482733 -0.000000034517 0.00025300 0.00000217 -5418.5325483636 -0.000000090325 0.00006686 0.00000064 -5418.5325483652 -0.000000001553 0.00003014 0.00000026 **** Energy Check signals convergence **** 0.0001158 0.0000 0.0000179 0.0000 0.0000112 0.0000 ***************************************************** * SUCCESS * * SCF CONVERGED AFTER 10 CYCLES * ***************************************************** + Năng lượng cấu trúc tối ưu hóa vịng (dịng 10377) Total Energy : -5418.53271646 Eh -147445.77121 eV lu an + Năng lượng điểm đơn cấu trúc tối ưu hóa vịng (dòng 10397- 10399) n va 5418.532716458656 gh tn to FINAL SINGLE POINT ENERGY - … p ie - Cấu trúc hình học sau thực phép tính vịng (dịng 10514 – 10522) d oa nl w ************************************************************* * GEOMETRY OPTIMIZATION CYCLE * ************************************************************* CARTESIAN COORDINATES (ANGSTROEM) Fe 0.000000 0.000000 1.392610 Ge 0.000000 1.232408 -0.530888 Ge 0.000000 -1.232408 -0.530888 v an lu + Số vòng lặp SCF: 15 vòng (dòng 462- 527) 0.000000000141 0.00000247 0.00000004 0.000000000093 0.00000093 0.00000002 -0.000000000303 0.00000089 0.00000002 ***DIIS convergence achieved*** 0.0000024 0.0000 0.0000013 0.0000 0.0000014 0.0000 oi m ll -5418.5340015389 -5418.5340015388 -5418.5340015391 fu an 12 13 14 nh ***************************************************** * SUCCESS * * SCF CONVERGED AFTER 15 CYCLES * ***************************************************** at z z @ + Năng lượng cấu trúc tối ưu hóa vịng (dòng 12396) : -5418.53421218 Eh gm Total Energy -147445.81191 eV Lu 5418.534212177907 om FINAL SINGLE POINT ENERGY - l.c + Năng lượng điểm đơn cấu trúc tối ưu hóa vịng (dịng 12416- 12418) an - Phép tính tối ưu hóa hình học chuyển đổi (dịng 12515- 12517) va ***********************HURRAY******************** n a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an P5 *** THE OPTIMIZATION HAS CONVERGED *** ************************************************* - Cấu trúc hình học FeGe2− sau vịng (dịng 12535 – 12544) ******************************************************* *** FINAL ENERGY EVALUATION AT THE STATIONARY POINT *** *** (AFTER CYCLES) *** ******************************************************* CARTESIAN COORDINATES (ANGSTROEM) Fe 0.000000 0.000000 1.393250 Ge 0.000000 1.231500 -0.531208 Ge 0.000000 -1.231500 -0.531208 an lu + Số vòng lặp SCF: 10 vòng (dòng 12739- 12761) n va -5418.5340021062 -5418.5340021150 -5418.5340021154 gh tn to 0.000000004294 0.00007884 0.00000068 -0.000000008778 0.00000982 0.00000012 -0.000000000410 0.00000668 0.00000006 ***DIIS convergence achieved*** 0.0000265 0.0000 0.0000104 0.0000 0.0000025 0.0000 p ie ***************************************************** * SUCCESS * * SCF CONVERGED AFTER 10 CYCLES * ***************************************************** w + Năng lượng cấu trúc tối ưu hóa sau vịng (dịng 12801) Total Energy : -5418.53421267 Eh -147445.81193 eV d oa nl + Năng lượng điểm đơn cấu trúc tối ưu hóa sau vịng kết thúc phép tính v an lu tối ưu hóa hình học (dịng 22350- 22354) 5418.534212672110 -*** OPTIMIZATION RUN DONE *** m ll fu an FINAL SINGLE POINT ENERGY - - Tần số dao động phân tử tối ưu hóa hình học (dịng 22456 - 22468) oi z @ om l.c gm cm**-1 cm**-1 cm**-1 cm**-1 cm**-1 cm**-1 cm**-1 cm**-1 cm**-1 z 0.00 0.00 0.00 0.00 0.00 0.00 192.36 231.68 317.94 at 0: 1: 2: 3: 4: 5: 6: 7: 8: nh VIBRATIONAL FREQUENCIES - Lu - Các kiểu dao động phân tử tối ưu hóa hình học học FeGe2− (dịng an 22456 - 22468); Trong đó, kiểu dao động từ đến (hàng ngang) nguyên tử Fe n va Ge đứng yên a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an P6 -NORMAL MODES an lu n va matrix 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 p ie gh tn to These modes are the cartesian displacements weighted by the diagonal M(i,i)=1/sqrt(m[i]) where m[i] is the mass of the displaced atom Thus, these vectors are normalized but *not* orthogonal 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000002 0.703812 0.000001 -0.551876 0.000002 -0.731841 0.000000 0.000000 0.000000 0.550130 -0.270754 -0.391078 0.212305 0.423103 0.281536 0.000000 0.000000 0.000000 -0.550131 -0.270751 0.391077 0.212302 -0.423105 0.281535 - Cường độ phổ hồng ngoại kiểu dao động (dòng 22501 – 22512) w d oa nl IR SPECTRUM - fu an v an lu Mode freq (cm**-1) T**2 TX TY TZ 6: 192.36 1.685862 ( 0.000000 -0.000210 -1.298407) 7: 231.68 0.013101 ( -0.000000 -0.114461 -0.000040) 8: 317.94 0.064897 ( -0.000000 -0.000227 0.254750) The first frequency considered to be a vibration is The total number of vibrations considered is m ll - Enthalpy (dòng 22569 – 22578) oi nh -ENTHALPY at z z The enthalpy is H = U + kB*T kB is Boltzmann's constant Total free energy -5418.52820598 Eh Thermal Enthalpy correction 0.00094421 Eh 0.59 kcal/mol Total Enthalpy -5418.52726177 Eh @ om l.c gm - Entropy (dòng 22585 – 22602) an Lu ENTROPY - n va The entropy contributions are T*S = T*(S(el)+S(vib)+S(rot)+S(trans) S(el) - electronic entropy a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an P7 S(vib) - vibrational entropy S(rot) - rotational entropy S(trans)- translational entropy The entropies will be listed as mutliplied by the temperature to get units of energy Electronic entropy 0.00130891 Eh 0.82 kcal/mol Vibrational entropy 0.00256601 Eh 1.61 kcal/mol Rotational entropy 0.01217415 Eh 7.64 kcal/mol Translational entropy 0.01986020 Eh 12.46 kcal/mol Final entropy term 0.03590927 Eh 22.53 kcal/mol - Năng lượng tự Gibbs (dòng 22632 – 22641) an lu GIBBS FREE ENTHALPY - n va The Gibbs free enthalpy is G = H - T*S gh tn to Total enthalpy -5418.52726177 Eh Total entropy correction -0.03590927 Eh -22.53 kcal/mol Final Gibbs free enthalpy -5418.56317104 Eh p ie - Phép tính kết thúc bình thường dịng 22655, thời gian tính 29 phút 29 giây w ****ORCA TERMINATED NORMALLY**** TOTAL RUN TIME: days hours 29 minutes 29 seconds 892 msec d oa nl Tính lượng điểm đơn FeGe– với CASSCF/CASPT2 phần v an lu mềm Molcas 3.1 Tập tin đầu vào (INPUT) trạng thái 2A2 fu an Nội dung ý nghĩa tập tin đầu vào: #Bắt đầu vòng lặp input số mục #Bắt đầu vịng lặp input để tối ưu hóa hình học oi *>>> Do while m ll >>> FOREACH R in (2.232) nh at #Chương trình SEWARD #phép tính đối xứng #Đối xứng qua mặt phẳng yz, xz #Từ khóa buộc chương trình phải tn theo u cầu input mà không tự động thay đổi z z @ &SEWARD &END Symmetry x y expert l.c gm # Douglas−Kroll #SEWARD chọn phân tách tích phân hai electron cách sử dụng độ xác cao (ngưỡng 1.0d-8) phân tách cấu hình Mặc định không phân tách om doug high an Lu n va a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an P8 Basis set Fe… / Inline Fe.aug-cc-pwCVTZ-Dk …23s19p11d4f3g.10s9p7d4f3g Fe 0 Angstrom End of basis set Basis Set Ge… / Inline Ge.aug-cc-pVTZ-Dk…21s14p10d2f.7s6p4d2f Ge 0 $R Angstrom End of Basis Set End of Input #Bộ hàm sở #Bộ hàm sở Fe #Hệ tọa độ xyz Fe #kết thúc hàm sở #Bộ hàm sở #Bộ hàm sở Ge #Hệ tọa độ xyz Ge #kết thúc hàm sở #kết thúc input >>copy/home/tan/calc/fege/casvd/2a2/223 #copy file “2232a.rasorb” theo 2a.rasorb INPORB đường dẫn an lu n va #phép tính RASSCF #đối xứng #độ bội spin gh tn to &RASSCF &END Symmetry Spin Nactel 11 0 p ie Frozen 0 0 w Inactive 13 5 d oa nl v an lu Ras2 3 fu an Lumorb oi m ll Linear Orbappear nh at Full End of input #số electron hoạt động (11), số lỗ trống Ras1(0), số e chiếm tối đa Ras3(0) #số lượng orbital không tương quan đối xứng (a1, b2 , b , a ) #số lượng orbital không hoạt động đối xứng a1, b2, b1, a2 13,5,5,1 (tổng 24) #số lượng orbital Ras2 đối xứng a1, b2, b1, a2 6,3,3,2 (tổng 14) #input orbital phân tử định dạng FORTRAN với tên tệp INPORB # phân tử hình đường thẳng #yêu cầu mục với xuất orbital output #chọn hình thức bảng #kết thúc input z z >>copy $Project.RasOrb # copy file từ đường dẫn /home/tan/calc/fege/casvd/2a2/2232a.ras orb >>copy $Project.rasscf.molden /home/tan/calc/fege/casvd/2a2/2232a.mol den @ l.c gm #phép tính CASPT2 #số lượng orbital đóng băng đối xứng a1, b2, b1, a2 8,3,3,0 (tổng số orbital 14) # Số vòng lặp thực tối đa 500 om an Lu &CASPT2 &END Frozen 3 n va Maxiter 500 a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an P9 #thêm đại lượng giả định sau sử dụng giá trị thực kết lượng bậc hai Thay đổi giả định chọn 0,1 #kết thúc input ImaginaryShift 0.1 End of input # Kết thúc vòng lặp >>> ENDDO 2.2 Tập tin đầu (output) Tập tin output bao gồm 1670 dịng với thơng tin quan trọng sau: an lu - Tập tin đầu vào (input) từ dịng 79- 358 - Phép tính SEWARD dòng 364 va - Start Module: seward at Thu Apr 26 15:43:38 2018 - n p ie gh tn to + Các phép tính SEWARD (dòng 384 – 395) d oa nl w SEWARD will generate: Multipole Moment integrals up to order Kinetic Energy integrals Nuclear Attraction integrals (point charge) One-Electron Hamiltonian integrals Velocity integrals Relativistic Douglas-Kroll integrals Cholesky decomposed two-electron repulsion integrals - CD Threshold: 0.10E-07 v an lu Integrals are discarded if absolute value > Do while n va #Bắt đầu vịng lặp input để tối ưu hóa hình học #Chương trình SEWARD #phép tính đối xứng #Đối xứng qua mặt phẳng yz, xz #Từ khóa buộc chương trình phải tuân theo yêu cầu input mà không tự động thay đổi p ie gh tn to &SEWARD &END symmetry x y expert d oa nl w doug high v an lu oi m ll fu an Basis set Fe / Inline Fe.aug-cc-pwCVTZ-Dk 23s19p11d4f3g.10s9p7d4f3g Fe 0.00000000 0.00000000 1.44948825 Angstrom End of basis set Basis Set Ge / Inline Ge.aug-cc-pVTZ-Dk 21s14p10d2f.7s6p4d2f Ge 0.00000000 1.16194804 -0.61391525 Angstrom End of Basis Set End of Input # Douglas−Kroll #SEWARD chọn phân tách tích phân hai electron cách sử dụng độ xác cao (ngưỡng 1.0d-8) phân tách cấu hình Mặc định khơng phân tách #Bộ hàm sở #Bộ hàm sở Fe #Hệ tọa độ xyz Fe nh at #kết thúc hàm sở #Bộ hàm sở z z @ #Bộ hàm sở Ge #Hệ tọa độ xyz Ge gm #kết thúc hàm sở #kết thúc input om l.c >>>copy/home/tan/calc/fege2/ras/c4a2rl.ra sorb INPORB #copy file “c4a2rl.rasorb” theo đường dẫn &RASSCF &END Symmetry #phép tính RASSCF #đối xứng Spin #độ bội spin an Lu n va a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an P14 Nactel 13 #số electron hoạt động (13), số lỗ trống Ras1(0), số e chiếm tối đa Ras3(2) #số lượng orbital đóng băng đối xứng (a1, b2, b1 , a ) #số lượng orbital không hoạt động đối xứng a1, b2, b1, a2 16, 6,13,4 (tổng 39) #số lượng orbital Ras2 đối xứng a1, b2, b1, a2 6,3,3,2 (tổng 14) #input orbital phân tử định dạng FORTRAN với tên tệp INPORB # phân tử hình đường thẳng #yêu cầu mục với xuất orbital output #chọn hình thức bảng #kết thúc input Frozen 0 0 Inactive 16 13 an lu n va gh tn to Ras2 *5 1 Ras3 *2 2 1 Lumorb Orbappear Full Levelshft 1.2 End of input p ie >>> copy $Project.RasOrb /home/tan/calc/fege2/ras/c4a2rl.rasorb >>> copy $Project.rasscf.molden /home/tan/calc/fege2/ras/c4a2rl.molden w # copy file từ đường dẫn d oa nl &CASPT2 &END Frozen 10 Maxiter 500 ImaginaryShift 0.1 End of input nh at *>>> EndDo oi *&SLAPAF &END *End of input m ll *&ALASKA &END *End of input fu an v an lu #phép tính CASPT2 #số lượng orbital khơng tương quan đối xứng a1, b2, b1, a2 8,3,3,0 (tổng số orbital 14) # Số vòng lặp thực tối đa 500 #thêm đại lượng giả định sau sử dụng giá trị thực kết lượng bậc hai #kết thúc input z z # Kết thúc vòng lặp @ gm 3.2 Tập tin đầu (OUTPUT) - Nội dung input từ dịng 79-357 - Phép tính SEWARD dòng 364 Lu - Start Module: seward at Sun Jul om l.c Tập tin output bao gồm 1947 dịng với thơng tin quan trọng sau: 21:56:16 2018 - an + Các phép tính SEWARD (dòng 381 – 393) n va a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an P15 + Các thông tin nhóm đối xứng (từ dịng 394-407) + Các thơng tin hàm sở (từ dòng 414-446) + Các thơng tin cấu trúc phân tử (từ dịng 450-460) + Đối xứng hàm sở (từ dòng 485-487) Basis set specifications : Symmetry species a1 Basis functions 89 b1 50 b2 69 a2 37 + Phép tính SEWARD kết thúc hoàn thành phút 21 giây (từ dòng 527-528) - Stop Module: seward at Sun Jul 21:57:38 2018 /rc=_RC_ALL_IS_WELL_ - Module seward spent minute 21 seconds - lu an - Phép tính RASSCF dịng 532 - Start Module: rasscf at Sun Jul 21:57:38 2018 - va n + Thông tin cấu trúc phân tử (hệ tọa độ Descartes) (dòng 552 – 560) gh tn to + Các thơng tin tối ưu hóa cấu trúc (từ dòng 588-606) ++ Optimization specifications: p ie w RASSCF algorithm: LK Cholesky Maximum number of macro iterations Maximum number of SX iterations Threshold for RASSCF energy Threshold for max MO rotation Threshold for max BLB element Level shift parameter Make Quasi-Newton update d oa nl 200 100 0.100E-07 0.100E+00 0.100E-03 0.500E+00 v an lu m ll fu an The MO-coefficients are taken from the file: INPORB Title: RASSCF average (pseudo-natural) orbitals Total molecular charge -1.00 oi + Phép tính đối xứng (b1) gồm: số vòng lặp RASSCF: thống kê lượng nh at độ hội tụ (từ dòng 631-802); số chiếm orbital hoạt động spin vỏ z mở (từ dòng 812-838) z @ Convergence after 162 iterations 163 -465.53727810 -0.37E-08 -0.26E-04 0.47 1.52 QN YES 0:02:47 40 l.c gm -0.59E-03 om ************************************************************************* Wave function printout: occupation of active orbitals, and spin coupling of open shells (u,d: Spin up or down) an Lu n va ************************************************************************* Note: transformation to natural orbitals has been made, which may change the order of the CSFs a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an P16 printout of CI-coefficients larger than 0.05 for root energy= -5465.537278 conf/sym 11111 22 333 44 11 Coeff Weight 812 u220u 22 200 u0 00 0 0.85686 0.73421 827 ud2uu 22 200 u0 00 0 0.05823 0.00339 Natural orbitals and occupation numbers for root sym 1: 0.998442 1.879021 1.963648 0.087493 1.037515 0.008847 0.007677 sym 2: 1.909276 1.937242 0.020773 sym 3: 1.815077 0.167543 0.030603 0.008584 sym 4: 1.035254 0.086048 0.006958 + Năng lượng cuối thu từ phép tính RASSCF (dòng 906) :: RASSCF root number Total energy: -5465.53727810 lu an + Kết cuối thơng tin orbital, tối ưu hóa cấu trúc (dịng 576- 589) Orbital specifications: - n va ++ to Symmetry species p ie gh tn a1 16 66 89 d oa nl w Frozen orbitals Inactive orbitals Active orbitals RAS1 orbitals RAS2 orbitals RAS3 orbitals Secondary orbitals Deleted orbitals Number of basis functions b1 41 50 b2 13 52 69 a2 30 37 v an lu + Orbital phân tử thu từ phép tính RASSCF loại đối xứng (từ dịng 9161462) điện tích phân tử, độ bội spin (từ dòng 1472-1577) fu an + Phép tính RASSCF kết thúc hồn thành 12 53 phút giây (từ dòng m ll 1601-1602) oi nh - Stop Module: rasscf at Mon Jul 10:50:42 2018 /rc=_RC_ALL_IS_WELL_ - Module rasscf spent 12 hours 53 minutes seconds - at - Phép tính RASPT2 dịng 1608 z 10:50:42 2018 - @ gm + Phép tính RASPT2 (từ dòng 1630 – 1715) z - Start Module: caspt2 at Mon Jul FINAL CASPT2 RESULT: an n va -5465.5372780994 -0.9379521220 -0.0015052479 -0.9394573698 -5466.4767354692 0.0000006564 Lu Reference energy: E2 (Non-variational): Shift correction: E2 (Variational): Total energy: Residual norm: om l.c + Kết cuối phép tính RASPT2 (dòng 1717- 1730) a th c si Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn C.vT.Bg.Jy.Lj.Tai lieu Luan vT.Bg.Jy.Lj van Luan an.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj Do an.Tai lieu Luan van Luan an Do an.Tai lieu Luan van Luan an Do an Stt.010.Mssv.BKD002ac.email.ninhd.vT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.LjvT.Bg.Jy.Lj.dtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn.Stt.010.Mssv.BKD002ac.email.ninhddtt@edu.gmail.com.vn.bkc19134.hmu.edu.vn