Trong việc tính toán cấu trúc, tính chất của phân tử, các chất rắn, lỏng và động lực học các quá trình hóa học bên cạnh hoá lượng tử dựa trên cơ sở hàm sóng đã chứng kiến hai hướng phát triển nổi bật : động lực phân tử và lý thuyết phiếm hàm mật độ. Động lực phân tử (molecular dynamics – MD) là lời giải bằng số của phương trình Newton cổ điển của một hệ N hạt, có khả năng cung cấp các đại lượng vĩ mô tính toán rất gần với thực nghiệm, một bức tranh chi tiết về cấu trúc và các quá trình xảy ra trong thế giới vi mô của một hệ phức tạp các nguyên tử mà bằng các năng lực thực nghiệm hiện nay chưa xác định được.
Bộ giáo dục đào tạo Trường Đại học bách khoa hµ néi - Luận văn thạc sĩ khoa học Nghiên cứu hấp phụ phân ly NO AL(111) phương pháp động lực phân tử bán lượng tử Khiếu Thị hương chi Hà nội 2005 Bộ giáo dục đào tạo Trường Đại học bách khoa hà nội - Luận văn thạc sĩ khoa học Nghiên cứu hấp phụ phân ly NO AL(111) phương pháp động lực phân tử bán lượng tử Nghành : Hoá lý thuyết hoá lý Mà số : 62.44.3101 Khiếu Thị hương chi Người hướng dẫn khoa học : GS TSKH Đặng ứng Vận Hà nội 2005 Lời cảm ơn Sau thời gian làm việc nghiên cứu với bảo, hướng dẫn thầy giáo, giúp đỡ nhiều bạn bè, đồng nghiệp, động viên, khích lệ người thân đà hoàn thành đề tài nghiên cứu Tôi xin bày tỏ lòng biết ơn sâu sắc tới tất người đà bên : Trước tiên, em xin gửi lời cảm ơn tới thầy giáo GS TSKH Đặng ứng Vận, người thầy đà hướng dẫn em từ ngày với phương pháp, hướng nghiên cứu Thầy giáo đà truyền cho em lòng say mê víi khoa häc vµ híng dÉn em hoµn thµnh tèt đề tài nghiên cứu Em xin cảm ơn tất thầy giáo, cô giáo tổ môn Hoá lý trường Đại học Bách Khoa Hà Nội, thầy cô đà giảng dạy, tận tình bảo cho em st thêi gian häc tËp T«i xin gưi lời cảm ơn tới quan tâm , giúp đỡ tạo điều kiện Trung tâm Sau đại học trường Đại học Bách Khoa Hà Nội, khoa Sau đại học trường Đại học Hàng Hải Việt Nam đà dành cho trình học tập thuận lợi, có nhiều điều kiện để phát triển Tôi muốn nói lời cảm ơn tới bạn Trung tâm ứng dụng Tin học Hoá học-Khoa Hoá Học-Đại học Khoa học Tự nhiên-Đại học Quốc Gia Hà Nội, người đà trực tiếp giúp trình thực đề tài Cảm ơn tất bạn đồng nghiệp, người thân đà bên cạnh động viên khích lệ Cuối cùng, xin gửi lời cảm ơn nơi đà dành cho nhiều thời gian, nhiều thuận lợi để tập trung vào việc nghiên cứu đề tài này, trường THPT Đồng Hoà Hải Phòng nơi công tác, cảm ơn bạn đồng nghiệp đà giúp đỡ nhiều cho Hải Phòng, ngày 15 tháng 11 năm 2005 Người viết Khiếu Thị Hương Chi Danh mục bảng, hình Bảng 1: So sánh SQMD với DRC , QMD MD Hình 1.1: Minh hoạ phương trình Newton Hình 1.2: Thuật toán bước nhảy ếch để tính tích phân phương trình Newton Hình 1.3: Mẫu tương tác site-site hai phân tử lưỡng nguyên Hình 1.4: Các loại lượng mẫu MM Hình 1.5: Mặt cắt theo hai tham số khoảng cách Hình 1.6: Đường đẳng lượng theo hai tham số khoảng cách Hình 1.7: Các phương pháp đường tuyến tính đồng thời (LST) đường bậc hai đồng thời(QST) Hình 1.8: Phương pháp điều khiển toạ độ xác định cấu trúc chuyển tiếp Hình 1.9: Sơ đồ biểu diễn trình bề mặt Hình 1.10: Phân tử tiếp cận bề mặt Hình 1.11: Hấp phụ không phân ly Hình 1.12: Hấp phụ phân ly Hình 1.13: Mô tả trình hấp phụ vật lý hấp phụ hoá học Hình 1.14: Sự phụ thuộc hàng rào lượng vào khoảng cách Hình 1.15: Nguyên tử hấp phụ hoá học, hấp phụ vật lý không hoạt hoá Hình 1.16: Nguyên tử hấp phụ hoá học, trạng thái hấp phụ vật lý có hoạt hoá Hình 1.17: Năng lượng xác định khả phản ứng Hình 1.18: Năng lượng dính N, O Hình 1.19: Đường đẳng lượng theo tham số khoảng cách Bảng 2.1: Toạ độ nguyên tử hệ nghiên cứu Hình 2.1 : Cấu hình hộp mô SQMD hệ phân tử NO lớp mỏng Al(111) Hình 3.1 : Cấu hình hệ ban đầu Hình 3.2 : Cấu hình hồi phục trường hợp ghi lại hình 3.2(a,b,c,d) Hình 3.3 : Khảo sát ảnh hưởng thay đổi sè bíc håi phơc (300, 500, 1000 bíc) H×nh 3.4 : Quỹ đạo chuyển dịch phân tử quét b»ng Rasmol H×nh 3.5 : CÊu h×nh ci cđa hƯ sau trình hồi phục Hình 3.6 : Sự biến đổi lượng phân tử NO qua mô 300 bước Hình 3.7 : Sự hấp phụ phân ly phân tử NO(1) Al(111) bước 28 Hình 3.8 : Sự hấp phụ không phân ly phân tử NO(1) Al(111) bước 235 Hình 3.9 : Sù hÊp phơ ph©n ly cđa ph©n tư NO(1) Al(111) bước 242 Hình 3.10 : Sự hấp phụ phân ly phân tử NO(1) Al(111) bước 83 Hình 3.11 : Sự hấp phụ phân ly phân tử NO(1) Al(111) bước 100 Hình 3.12 : Sù hÊp phơ ph©n ly cđa ph©n tư NO(1) Al(111) bước 134 Hình 3.13 : Sự hấp phụ phân ly phân tử NO(1) Al(111) bước 171 Hình 3.14 : Sự hấp phụ phân ly phân tử NO(1) Al(111) bước 205 Hình 3.15 : Sự hấp phụ phân ly phân tử NO(3) Al(111) bước 265 Hình 3.16 : Sự hấp phụ phân ly phân tử NO(3) Al(111) ë bíc 295 H×nh 3.17 : Sù hÊp phơ phân ly phân tử NO(3) Al(111) bước 300 H×nh 3.18 : CÊu h×nh cuèi NO tiÕp cận nguyên tử Al Hình 3.19 : Mặt cắt đường đẳng lượng theo tham số khoảng cách Hình 3.20 : Cấu hình cuối NO tiếp cận nguyên tử Al Hình 3.21 : Mặt cắt đường đẳng lượng theo tham số khoảng cách Hình 3.22 : Cấu hình cuối NO song song Al-Al Danh mục bảng hình Hình 3.23 : Mặt cắt đường đẳng lượng theo tham số khoảng cách Hình 3.24 : Cấu hình cuối NO vuông góc với Al-Al Hình 3.25 : Mặt cắt đường đẳng lượng theo tham số khoảng cách Danh mục bảng hình Các ký hiệu, viết tắt MD-Molecular Dynamics : Động lực phân tử DFT-Density Functional Theory : Lý thuyết phiếm hàm mật độ QMD-Quantum Molecular dynamics : Kỹ thuật mô động lực phân tử lượng tử SQMD-Semi Quantum Molecular Dynamics : Động lực phân tử bán lượng tử MM-Molecular Mechanics : Mẫu học phân tử SQUARED -Semi - QUantum Approximation of REaction Dynamics vi : Vận tốc Fi : Lực tương tác pi : động lượng i K : động hệ : tương tác E : lượng toàn phần hệ s : số mômen KB : số Boltzman AB(R) : tương tác cặp phân tử ij (rij) : tương tác cặp tâm tương tác i-j H : Hamiltonian hệ Qi : toạ độ Pi : moment Vone : hệ OMQD VXC: tương quan trao đổi hệ OMQD EXC : lượng tương quan ~ f i : tổng hợp lực cổ điển ~ f jm : lực tương tác cặp Các ký hiệu, viết tắt f kiint :gradient phân tử k theo toạ độ i : khả hấp phụ rơi vào bẫy : hệ số nhiệt động thuận lợi cho trình hấp phụ C : phân tử khí phân ly hấp phụ lên bề mặt : khả phản ứng SC : hệ số dính : sè ph©n tư che phđ / sè ph©n tư phản ứng , : tham số hồi phục Các ký hiệu, viết tắt -1- Mở đầu I Lý chọn đề tài Trong việc tính toán cấu trúc, tính chất phân tử, chất rắn, lỏng động lực học trình hóa học bên cạnh hoá lượng tử dựa sở hàm sóng ®· chøng kiÕn hai híng ph¸t triĨn nỉi bËt : động lực phân tử lý thuyết phiếm hàm mật độ Động lực phân tử (molecular dynamics - MD) lời giải số phương trình Newton cổ điển hệ N hạt, có khả cung cấp đại lượng vĩ mô tính toán gần với thùc nghiƯm, mét bøc tranh chi tiÕt vỊ cÊu tróc trình xảy giới vi mô hệ phức tạp nguyên tử mà lực thực nghiệm chưa xác định Bởi vì, phương pháp tính nguyên tử xem đối tượng sở (theo nghĩa phân tử nguyên tử tạo thành) kết thu được so sánh với số đại lượng thực nghiệm (mà tất đại lượng tính toán được) nên MD gọi phương pháp mô nguyên tử atomistic simulation Lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT) phương pháp biến phân có nguồn gốc từ vật lý trạng thái rắn để giải toán trạng thái sở hệ n electron tương tác với electron khác với n nhân, đà nhanh chóng trở thành tảng cho tính toán hoá lượng tử Kỹ thuật mô động lực phân tử lượng tử (Quantum Molecular dynamics - QMD) kết hợp hiệu hoá lượng tử, DFT MD với hy vọng có ưu hai phương pháp đồng thời khắc phục Mở đầu -2- nhược điểm phương pháp.Với phát triển gần đây, QMD thực máy tính cá nhân hệ mạng máy tính song song Bằng cách tách việc tính lực tương tác thành phần : phần lượng tử khoảng cách gần học phân tử khoảng cách xa Phương pháp động lực phân tử bán lượng tử (SQMD- Semi Quantum Molecular Dynamics ) đà rút ngắn đáng kể thời gian tính toán so với QMD Vì vậy, QMD SQMD trở thành phương tiện đầy triển vọng nghiên cøu vËt lý c¸c hƯ ngng tơ, khoa häc vËt liệu, hoá lý lý hóa Một phổ rộng ứng dụng hệ phân tử, vật liệu khối, bề mặt đà minh chứng cho khả phương pháp việc phân tích dự báo tính chất cân không cân Những năm gần nghiên cứu trình hấp phụ khí bề mặt kim loại quan tâm thực nghiệm lý thuyết, nghiên cứu trình bề mặt tính toán lý thuyết động lực phân tử lượng tử (QMD) quan tâm Sự hấp phụ phân ly phân tử từ pha khí bề mặt bước phản ứng bản, quan trọng xảy xúc tác, ăn mòn, đặc biệt xảy bình chứa khí kim loại Sự hấp phụ khí NO bề mặt Al kim loại tương đương với trình phân bố phân tử NO lên bề mặt tinh thể Al từ pha khí Những kết công bố luận văn thu nghiên cứu cấu trúc tối ưu kỹ thuật động lực phân tử bán lượng tử tính cho hệ NO Al(111) nhằm khảo sát chi tiết trình hấp phụ phân ly NO bề mặt tinh thể Al kim loại yếu tố ảnh hưởng đến trình Mở đầu 11 CALL ALLFORCE C CHOOSING THE ABSORPT MOLCULE ONLY IATOM=0 SINTERENER=0.D0 IATOM = 1, NSITS(1) SINTERFENER=SINTERENER+(EEX(IATOM)+EEY(IATOM)+EEZ(IATOM)) * /HARTREETOUNITE END DO E(ISTEPZ,ISTEPR) = EPOT+SINTERENER 222 CONTINUE !JOULE C PES CALCULATION 290 CONTINUE WRITE(NQ,'(100F10.3)') (E(ISTEPZ,ISTEPR),ISTEPR=1,NSTEPR) 300 CONTINUE close(NQ) open(UNIT= NQ,FILE= 'PICTURE.XYZ',STATUS = 'unknown') write(NQ,'(I4)')NAE(TASKID) WRITE(NQ,*) 'LAST CONFIGURATION AFTER ',MSTEP,' SIMULATION STEPS' DO I = NAB(TASKID),NAE(TASKID) WRITE(nq,'(A5,3(F12.4),I10)') * NM(NSITE(I)),SX(I),SY(I),SZ(I),NNUM(I) END DO ! theo I close(NQ) RETURN END Phơ lơc 12 Ch¬ng trình thiết lập toạ độ bề mặt Al(111) SUBROUTINE TOADOAl implicit double precision (A-H,O-Z) dimension x(100),y(100),z(100) open(unit=1,file='toadoAl.xyz', status='unknown') d0=4.05 a=d0*sqrt(2.0)/2.0 d=a*sqrt(3.0)/2.0 write(*,*)' nhap so ngtu hang va cot' read(*,*)nh, nc i=0 ih=1,nh ic=1,nc i=i+1 if (mod(ih,2).ne.0) then x(i)=(ic-1)*a else x(i)=(ic-1)*a +a/2.0 endif y(i)= (ih-1)*d z(i)=0.0 end end imax=i i=1,imax write(*,'(A4,3f12.4)')'Al', x(i),y(i),z(i) write(1,'(A4,3f12.4)')'Al', x(i),y(i),z(i) end close(1) end Phụ lục 13 Chương trình hồi qui xử lý kết thu tính mặt PROGRAM HOIQUI implicit double precision (A-H,O-Z) dimension x(100),y(100),a(100) dimension e(100,100) c m2= array(1 10,1 10) of real c m3= array(1 100,1 100) of real c m1= array(0 100) of real c BEGIN (* chuong trinh chinh *) c CLRSCR open(unit=1,file='pes_100d',status='old') read(1,*) x10, x20 read(1,*) nstep1,dstep1,nstep2,dstep2 i=1 , nstep2 c begin read(1,*) (e(i,j),j=1,nstep1) enddo m=5 istep2=1,nstep2 c begin icount=0 istep1=1,nstep1 c begin write(*,'(f10.2)') e(istep1,istep2) if (e(istep1,istep2).ne.0.0) then c begin icount=icount+1 x(icount)= x10+(istep1-1)*dstep1 y(icount)= e(istep1,istep2) endif enddo write(*,*) 'icount=',icount Phô lôc 14 call hoiquy_(ICOUNT,M,X,Y,A) n=m+1 istep1=1,nstep1 c begin if (e(istep1,istep2).eq.0.0) then c begin xx=x10+(istep1-1)*dstep1 ytt=0.0 j=1,n ytt= ytt+ a(j)* xx**(n-j) enddo write(*,'(3F10.4)')XX,E(istep1,istep2),ytt e(istep1,istep2)=ytt endif enddo WRITE(*,*)'istep2 =',istep2 c read(*,*) enddo close(1) open(unit=2,file='out',status='unknown') i=1 , nstep1 c begin write(2,'(100f10.3)') (e(i,k),k=1,nstep2) enddo close(2) end subroutine ghpttt(a,n,xx) implicit double precision (A-H,O-Z) dimension a(10,10),xx(100),aa(10,10) i=1,n k=1 , n+1 Phô lôc 15 aa(i,k)=a(i,k) ! { LUU A BANDAU} end end do i=1,n !* Khu x(i) *) c Begin skk=a(i,i) j=i,n+1 c Begin a(i,j) = a(i,j)/skk end if (i.eq.n) goto 10 k= i+1,n c begin j= i+1,n+1 a(k,j)=a(k,j)-a(k,i)*a(i,j) enddo enddo enddo 10 xx(n) = a(n,n+1) Write(*,*)' Nghiem cua phuong trinh' Write(*,*) ' Xx(',n,')= ', xx(n) i = (n - 1), 1,-1 c Begin ak=0 j = i + , n c Begin ak = ak + xx(j) * a(i,j) enddo XX(i) = a(i,n+1) - ak Write (*,*)' XX(',I,')= ',XX(i) enddo Phô lôc 16 c (* thu lai nghiem cua pt *) write(*,*)' Thu lai nghiem cua phuong trinh ' i=1 , n c begin test=0 j=1 , n c begin test= aa(i,j)*xx(j)+test enddo write(*,'(a10,f15.4,a10,i3,a6,f15.4)') * ' test=',test,' aa(',i,',n+1)=',aa(i,n+1) enddo return end subroutine HOIQUY_(npont,M,x,y,A) implicit double precision (A-H,O-Z) dimension x(100),y(100),a(100) dimension sumx(0:100),sumyx(0:100) dimension hs(10,10) integer twom c BEGIN k= 0,2*m c begin sumx(k)= i=1 , npont sumx(k)=sumx(k)+ x(i)**k enddo enddo k=0,m c begin sumyx(k)= Phô lôc 17 i=1 , npont sumyx(k)=sumyx(k)+ y(i)* x(i)**k enddo enddo n= m+1 twom= 2*m i=1,n c begin j=1,n hs(i,j)= sumx(twom -(i-1)-(j-1)) enddo hs(i,n+1)= sumyx(m-(i-1)) enddo write(*,*) 'OK HS' call GHPTTT(hs,n,a) c (* kiem tra nghiem *) write(*,*) ' ** thu lai phuong trinh hoi quy ** ' write(*,*)' ytt ytn delta ' i=1 , npont c begin ytt= k=1,n ytt= ytt+ a(k)* x(i)**(n-k) enddo write(*,'(3f10.4)') ytt,y(i),ytt-y(i) end return END Phô lôc 18 Chương trình tính trình hồi phục cấu trúc INPUT FILE : *TÝnh lỵng tư $CONTRL SCFTYP=UHF MAXIT=500 RUNTYP=ENERGY COORD=CART NZVAR=0 $END $SYSTEM TIMLIM=6000 $END $BASIS GBASIS=PM3 $END $GUESS GUESS=HUCKEL $END $DATA NO C1 N 7.0 0.0466 -0.8282 0.3483 O 8.0 -0.0466 0.3010 0.3211 Al 13.0 0.0466 -0.8282 1.7483 Al 13.0 0.0466 -0.8282 4.7483 Al 13.0 3.0466 -0.8282 1.7483 $END Inputmd *Tính động lực phân tử # Sample input file for the MD99 program # Simulated system consists of 36 Aluminium molecules and NO molecule # The specified time scale (5ps after 5ps equilibration) is not enogh # for obtaining "true" values Equilibration time for this system # is about 100ps # Lines beginning with "#" are commentaries # Output controll parameter # Base file name for output files: NO_Al # Path to the molecular database /moldata Phô lôc 19 # Read from Dump Check # restart file? restart file? f .t .t only? f # NVT ensemble? Zero counter of cpu time? t NPT ensemble? t Anizotropic NPT? f # Number of molecule types: # which molecules: no al # Number of molecules: 36 # - intramolecular interactions: f .f .f .f .f .f # - intramolecular interactions: f .f .f .f .f f # Scaling factor for 1-4 LJ interactions 1 1 1 # Scaling factor for 1-4 electrostatic interactions 1 1 1 # Intramolecular potential type 0 0 # temperature(K) 800 density(g/cm**3) 0.200 cell type 12.5 # Time step (s) 1.d-15 # Total # MD-steps 1000 pressure (atm) # box size (A) 15.8 12 Small steps in one long Steps for interme- take averages diate averaging each steps 500 # T-termostat param(fs) P-termostat(fs) 30 700 .T dump restart file after steps 5000 NVE+velosity 20 Phô lôc delta T (K) scaling? 20 # Rcutoff(A) Rcut-fast forces 10 check neigbours after steps 10 # Ewald parameters: # alpha fexp 3.14159256 9.81 # which molecules move: t .f .f .f .f .f # recalculate list of intramolec interactions? f # Constrain dynamics? f tolerance parameter 1.d-4 # Initial state (from -1 to 4) zero velocities? f # Keep initial orientation (0/1)? 0 0 0 # Radius of a hole Keep initial # for large molecules 10 Permissible configuration? f # Change temperature at restart? f deviation? Change density? f # Cut large forces? f Parameter 1.d6 Reacdat *TÝnh håi phôc ! kieu tinh toan: Kreaction ! so site dac trung cho moi phan tu a b ! cac site dac trung cua a ! cac site dac trung cua b 3.1 ! khoang cach luong tu ! tiep can be mat 0.2 0.1 ! lamda muy Phơ lơc 21 PES *TÝnh bỊ mỈt thÕ Trường hợp định hướng vuông góc -3.58330 -3.72015 ! XCOOR YCOOR ! IORIENT 0.0 ! ZMAX (khoang cach tu tam khoi den be mat) 0.6 ! RMIN (khoang cach giua nguyen tu N va O) 50 0.05 50 0.05 ! NSTEPZ,DSTEP1,NSTEPR,DSTEP2 ! iapproach Trêng hỵp ®Þnh híng song song -2.1514 -3.72015 ! XCOOR YCOOR ! IORIENT 0.0 ! ZMAX (khoang cach tu tam khoi den be mat) 0.6 ! RMIN (khoang cach giua nguyen tu N va O) 50 0.05 50 0.05 ! NSTEPZ,DSTEP1,NSTEPR,DSTEP2 ! iapproach 10 ! nguyen tu Al duoc xet tinh luong tu Phơ lơc Tµi liƯu tham khảo * Tiếng Việt [1] Trần Thái Hoà (2001), Nghiên cứu trình chuyển dịch nguyên tử phân tử phương pháp Động lực phân tử Hoá lượng tử, Luận án Tiến sĩ Hoá học, Trường ĐHKHTN ĐHQG Hà Nội [2] Trần Thành Huế (2001), Hóa học đại cương, Tập 1, NXB Giáo dục [3] Trung tâm ứng dụng Tin học Hoá học, Khoa Hoá học, Trường ĐHKHTN-ĐHQG Hà Nội (2002), Phần mềm tính động lực phân tử MD99 [4] Đặng ứng Vận (1998), Tin häc øng dơng Hãa häc, NXB Gi¸o dơc [5] Đặng ứng Vận (2001), Kỹ thuật mô vật liệu phương pháp động lực phân tử, NXB ĐHQG HN [6] Đặng ứng Vận (2003), Động lực học phản ứng hoá học, NXB Giáo dục [7] Đặng ứng Vận (2003), Tính toán động lực phân tử bán lượng tử phản ứng đồng phân hoá, Tạp chí hoá học số 2, tr 35-40 [8] Đặng ứng Vận (2003), Phần mềm SQUARED nghiên cứu phản ứng hoá học phương pháp động lực phân tử bán lượng tử, Tuyển tập báo cáo toàn văn hội nghị toàn quốc lần thứ đề tài nghiên cứu khoa học lĩnh vực hoá lý hoá lý thuyết, tr 222-231 [9] Trần Văn Nhân, Nguyễn Thạc Sửu, Nguyễn Văn Tuế (1998), Hoá lý tập II, NXB Giáo dục [10] Lê Chí Kiên(dịch) (1984), Cơ sở hoá học vô phần III, NXB Đại học trung học chuyªn nghiƯp * TiÕng Anh [11].Shailk S S., Schlegel H B., Wolfe S (1976), Theoretical aspects of physical organic chemistry John Wiley & Son, Inc [12] Andrew J Komrowki, at el (2002), Dissociative adsorption of NO upon Al(111): Orientation dependent charge tranfer and chemisorption reaction dynamics J Chem Phys Vol 117, No18 [13] Goldstein (1950), Classical Mechanics Addition Wesley, London [14] Bunday B.D (1984), Basic Optimization method, Edward Amold Press [15] Foresman J.B and Frisch (1993), “ About This book” in Exploring Chemistry with Electronic Structure Method, AGuide to Using Gaussian [16] Scales L.E (1985), Introduction to Non-Linear Optimization, MacMillan [17] Comeau D.C., Zellmer R.J and Shavitt (1986), In Goemetrical Drivatieves of Energy Surfaces and Molecular Properties, Edited by Jorgensen P and Simons J, Reidel., NATO.ASISes., C166,P243 [18] Pople J.A., Head-Gordon M., Fox D.J., Raghavachari K, and Curtiss L.A (1989), “ Gaussian-1 theory : A general procedure for prediction of molecular energies”, J.Chem.Phys, 90, pp.5622 [19] Alya P.Y and Schegel H.B (1997), “ A Combined Method for Detemining Reaction Paths, Minima and Transition State Geometries”, J.Chem.Phys,107,pp 375 [20] Chemical Processing of Electronic Materials Gas Adsorption Kinetics [21] Adsorption.htm [22] J Chem.Phys., Vol 117, No 18, November 2002, Dissociative adsorption of NO on Al(111), p8185, p8189 [23] M Bockstedte, A Kley, J Neugebauer and M Scheffler, Comp Phys.Commun.107,187(1997) [24] M Bellotto, B Rebous, O Clause, J Lynch, D Bazin, E Elkaim, J Phys Chem., 100, 8527 (1996) Môc lôc ... tích NO giống NO2 - Với nghiên cứu tổng quan tính toán động lực phân tử bán lượng tử, tính mặt lượng, tượng hấp phụ, hấp phụ phân ly hướng đà nghiên cứu hấp phụ phân ly với hệ nghiên cứu, luận văn. .. 28 Hình 3.8 : Sự hấp phụ không phân ly phân tử NO( 1) Al(111) bước 235 Hình 3.9 : Sự hấp phụ phân ly phân tử NO( 1) Al(111) bước 242 Hình 3.10 : Sự hấp phụ phân ly phân tử NO( 1) Al(111) bước 83... ly phân tử NO( 1) Al(111) bước 205 Hình 3.15 : Sự hấp phụ phân ly phân tử NO( 3) Al(111) ë bíc 265 H×nh 3.16 : Sù hÊp phơ phân ly phân tử NO( 3) Al(111) bước 295 Hình 3.17 : Sự hấp phụ phân ly phân