Đề tài Tính toán năng lượng electron bằng phần mềm Mathematica và áp dụng trong giảng dạy hóa học lớp 10

1 MỤC LỤC Trang Trang phụ bìa Lời cam ñoan Lời cảm ơn Mục lục ........................................................................................................ 1 MỞ ðẦU ..................................................................................................... 3 CHƯƠNG 1 CƠ SỞ LÝ THUYẾT 1.1. Hệ một electron, một hạt nhân (nguyên tử hidro)...................................... 5 1.1.1. Mô hình hệ............................................................................................... 5 1.1.2. Phương trình Schrodinger và sơ lược về cách giải.................................... 5 1.1.3. Một số kết quả thu ñược........................................................................... 6 1.1.3.1. Năng lượng của hệ.............................................................................. 6 1.1.3.2. Hàm bán kính ..................................................................................... 6 1.1.3.3. Hàm cầu l,mY ( , )θ φ .............................................................................. 7 1.2. Nguyên tử nhiều electron: ........................................................................... 8 1.2.1. Mô hình hệ............................................................................................... 8 1.2.2. Phương trình Schrodinger và sơ lược về cách giải.................................... 9 1.2.3. Một số kết quả thu ñược........................................................................... 9 1.2.3.1. Sự gần ñúng Slater.............................................................................. 10 1.2.3.2. Hàm sóng ........................................................................................... 11 1.3. Cấu hình electron......................................................................................... 11 1.3.1. Khái niệm ................................................................................................ 11 1.3.2. Những cơ sở ñể viết cấu hình electron..................................................... 11 1.3.2.1. Nguyên lý vững bền hay nguyên lý năng lượng cực tiểu..................... 11 1.3.2.2. Quy tắc Klechkovxki .......................................................................... 12 1.3.2.3. Nguyên lí Pauli................................................................................... 12 1.3.2.4. Qui tắc Hund ...................................................................................... 13 PDF Created with deskPDF PDF Writer - Trial :: 2 1.4. Cấu hình eletron và qui tắc gần ñúng Slater.............................................. 13 1.5. Năng lượng ion hóa và qui tắc gần ñúng Slater ......................................... 14 Chương 2 LẬP TRÌNH TÍNH TOÁN 2.1. Tính tổng năng lượng electron của nguyên tử............................................ 16 2.1.1. Mô hình tính ............................................................................................ 16 2.1.2. Chương trình tính..................................................................................... 17 2.1.2.1. Chương trình tính tổng năng lượng electron cho cấu hình 1 và 2 ........ 17 2.1.2.2. Chương trình tính tổng năng lượng electron cho cấu hình 3 và 4 ........ 21 2.2. Tính năng lượng ion hóa của các nguyên tử............................................... 25 2.2.1. Mô hình tính ............................................................................................ 25 2.2.2. Chương trình tính..................................................................................... 25 2.2.2.1. Tính năng lượng ion hóa thứ nhất ....................................................... 25 2.2.2.2 Tính các trị năng lượng ion hóa có thể của C....................................... 29 Chương 3 KẾT QUẢ VÀ THẢO LUẬN 3.1.Cấu hình electron của các nguyên từ họ d chu kì 3..................................... 34 3.1.1. Các nguyên tố họ d chu kì3 (ngoại trừ Cr và Cu)...................................... 34 3.1.2. Các nguyên tố Cr và Cu ........................................................................... 34 3.2. Năng lượng ion hóa thứ nhất của các nguyên tố chu kì 2 .......................... 35 3.3. Các trị năng lượng ion hóa có thể có của C................................................ 35 3.4. Một số nhận xét............................................................................................ 35 KẾT LUẬN......................................................................................................... 37 TÀI LIỆU THAM KHẢO.................................................................................. 38 PDF Created with deskPDF PDF Writer - Trial :: 3 MỞ ðẦU 1. Lý do chọn ñề tài Cấu hình electron của nguyên tử ñã trở thành một trong nhiều khái niệm trung tâm của hóa học hiện ñại. Từ cấu hình electron của nguyên tử, chúng ta có thể thu ñược rất nhiều thông tin hữu ích như vị trí của nguyên tố trong bảng hệ thống tuần hoàn, số electron hóa trị, số electron ñộc thân, số AO trống,… những thông tin này gần như không thể thiếu khi khảo sát sự hình thành liên kết hóa học giữa các nguyên tử, số hóa trị có thể có của một nguyên tử cũng như sự biến thiên tuần hoàn tính chất hóa học của thế giới vật chất. Trong dạy học hóa học cấu hình electron có vai trò quyết ñịnh ñến việc muốn nghiên cứu hay xác ñịnh các thông tin về một nguyên tố. Vì vậy khái niệm về cấu hình electron ñã ñược ñưa vào giảng dạy ở chương trình hóa học lớp 10 phổ thông trung học. Trên cơ sở này, học sinh có thể tiếp tục lĩnh hội các khái niệm hóa học tiếp theo trong chương trình hóa học phổ thông. Hiện nay khi viết cấu hình electron một số giáo viên, học sinh vẫn còn lúng túng và một số tranh cãi như:
Cách viết cấu hình của các nguyên tố họ d từ chu kì 3 trở ñi là (n+1)s nd hay nd (n+1)s?
Tại sao lại xuất hiện các cấu hình electron theo kiểu 3d54s1, 3d104s1?
Khi tách 1 electron ra khỏi nguyên tử thì năng lương ion hóa sẽ ñược tính như thế nào?
Có thể tách tất cả các electron ra khỏi nguyên tử hay không và năng lượng ion hóa ñược tính như thế nào? ðây chính là nguyên nhân chủ yếu làm cho việc viết cấu hình không chính xác, làm giảm sự hiểu biết, ảnh hưởng ñến chất lượng dạy và học hóa học ở phổ thông. Tất cả những vấn ñề trên có thể ñược giải thích trên cở sở hóa học lượng tử kết hợp với phần mềm toán học và máy tính ñiện tử. PDF Created with deskPDF PDF Writer - Trial :: 4 Với những lí do nêu trên, em quyết ñịnh chọn ñề tài cho khóa luận là: “Tính toán năng lượng electron bằng phần mềm Mathematica và áp dụng trong giảng dạy hóa học lớp 10” 2. Mục tiêu
Xác lập cách viết cấu hình electron hợp lí cho các nguyên tử
Tính các trị năng lượng ion hóa của các nguyên tử 3. Nhiệm vụ
Tìm hiểu các nội dung liên quan ñến việc viết cấu hình electron cho nguyên tử, cách tính mức năng lượng electron và năng lượng ion hóa thứ nhất.
Sử dụng phần mềm Mathematica lập trình tính năng lượng electron
ðưa ra cách viết cấu hình electron hợp lí
So sánh kết quả năng lượng ion hóa tính ñược với thực nghiệm 4. Phương pháp nghiên cứu Nghiên cứu lí thuyết hóa học luợng tử, viết chương trình tính trên cơ sở phần mềm Mathematica và so sánh kết quả tính ñược với thực nghiệm. 5. ðối tượng và phạm vi nghiên cứu 5.1. ðối tượng
Cách viết cấu hình electron cho các nguyên tử lớn
Năng lượng ion hóa cho các nguyên tử nhỏ 5.2. Phạm vi
Lập phương trình tính toán trên phần mềm Mathematica.
Tính năng lượng ion hoá thứ nhất cho các nguyên tố chu kì 1 và 2
Tìm ra cách viết cấu hình electron hợp lí cho các nguyên tố họ d. PDF Created with deskPDF PDF Writer - Trial :: 5 NỘI DUNG CHƯƠNG 1 : CƠ SỞ LÝ THUYẾT 1.1. Hệ một electron, một hạt nhân (nguyên tử hidro) Nguyên tử hidro là mô hình hệ ñơn giản nhất của hóa học lượng tử mà phương trình Schrodinger có thể giải chính xác. Việc giải phương trình Schrodinger cho hệ lượng tử này sẽ ñưa ñến những khái niệm cơ bản nhất của hóa học lượng tử. 1.1.1. Mô hình hệ Hệ lượng tử này bao gồm hạt nhân nguyên tử hidro với ñiện tích là +1 và một electron chuyển ñộng xung quanh hạt nhân này. Xem hạt nhân là gốc của hệ tọa ñộ, vị trí của electron so với gốc ñó ñược xác ñịnh bằng vectơ vị trí r
Hình I.1. Mô hình hệ một electron – một hạt nhân trong hệ tọa ñộ cầu 1.1.2. Phương trình Schrodinger và sơ lược về cách giải Toán tử Hamilton mô tả trạng thái của hệ có dạng: e Zeˆ ˆˆH T U m x y z r 2 2 2 2 2 0 2 2 2  ∂ ∂ ∂ = + =− + + − 2 ∂ ∂ ∂  ℏ Hàm sóng mô tả trạng thái của hệ là (r)Ψ  ðể ñơn giản bài toán ta chuyển hệ tọa ñộ ðecac sang hệ tọa ñộ cầu: z r cos y r sin sin x r sin cos = θ = θ φ = θ φ PDF Created with deskPDF PDF Writer - Trial :: 6 Hàm sóng của hệ ñược viết lại như sau: (r) (r, , ) R(r)Y( , )Ψ = Ψ θ φ = θ φ  R(r) ñược gọi là hàm bán kính Y( , )θ φ ñược gọi là hàm góc hay hàm cầu Phương trình Schrodinger của hệ có dạng: Hˆ EΨ = Ψ Giải phương trình Schodinger cho hệ một lectron, một hạt nhân thu ñược hàm riêng , , ln l m ψ - nói chung là hàm phức – và trị riêng tương ứng là năng lượng En cùng bộ ba số lượng tử . Dùng phép phân li biến số thích hợp (dựa vào ñặc ñiểm của toán tử Hˆ và hàm sóng (r, , )Ψ θ φ ), phương trình Schrodinger sẽ ñược tách thành 2 phương trình riêng rẽ là phương trình góc và phương trình bán kính. 1.1.3. Một số kết quả thu ñược Tiến hành giải riêng lẽ 2 phương trình góc và phương trình bán kính, thu ñược một số kết quả sau: 1.1.3.1. Năng lượng của hệ Năng lượng của hệ khi giải hàm bán kính ñược tính theo biểu thức sau: e n m Z e E n ( ) 2 4 0 2 2 2 0 = − 2 4πεℏ n : số lượng tử chính 0ε : ñộ thẩm từ trong chân không Hàm bán kính n,lR (r): ðể ñơn giản trong việc mô tả hàm sóng, hệ ñơn vị nguyên tử ñược sử dụng: e e m 0 = 1 = 1 0 = 1 4πε = 1 ℏ a0 = 1 1.1.3.2. Hàm bán kính Biểu thức của một số hàm bán kính như sau: PDF Created with deskPDF PDF Writer - Trial :: 7 Zr n 2 ρ = Orbital 1s /,R (r) Z e3 −ρ 21 0 = 2 Orbital 2s /,R (r) Z ( )e3 −ρ 22 0 1 = 2−ρ 8 Orbital 2p /,R (r) Z e3 −ρ 22 1 1 = ρ 24 Orbital 3s /,R (r) Z ( )e3 2 −ρ 23 0 1 = 6−6ρ + ρ 243 Orbital 3p /,R (r) Z ( ) e3 −ρ 23 1 1 = 4− ρ ρ 486 Orbital 3d /,R (r) Z e3 2 −ρ 23 2 1 = ρ 2430 1.1.3.3. Hàm cầu l,mY ( , )θ φ Electron chuyển ñộng trong nguyên tử là một hệ thực, do ñó cần phải chuyển obitan nguyên tử phức thì obitan nguyên tử thực ñể mô tả chuyển ñộng của hệ này. Kết quả giải phương trình Schrodinger cho phần góc sẽ thu ñược các hàm cầu thực hoặc phức: ,Y ( , )0 0 1 θ φ = 4π ,Y ( , ) cos1 0 3 θ φ = 4π θ i ,Y ( , sin) e φ 1 1 3 8π θθ φ = ðể mô tả chuyển ñộng thực của electron ta phải chuyển các hàm cầu phức thành các hàm cầu thực. Kết quả thu ñược các hàm cầu thực như sau: s 1 = 4π z sp co 3 = 4π θ PDF Created with deskPDF PDF Writer - Trial :: 8 xp si c sn o 3 = θ 4π φ yp sin sin 3 = θ 4π φ z d ( cos )2 25= 3 θ−1 16π xzd sin o osc s c 15 = θ 4 θ π φ yzd sin o inc s s 15 = θ 4 θ π φ x y os sin )d sin (c2 2 2 2 − 2φ− 15 = θ 16 φ π xyd sin cos s in 215= θ 4π φ φ 1.2. Nguyên tử nhiều electron ðối với nguyên tử nhiều electron, về nguyên tắc các hàm sóng ψ và năng lượng En tương ứng của cả vỏ nguyên tử ñược tạo ra bởi các electron có thể thu ñược từ việc giải phương trình Schrodinger. 1.2.1. Mô hình hệ Hệ lượng tử này gồm có hạt nhân nguyên tử mang ñiện tích +Z và từ 2 electron trở lên. Trong hệ này, hạt nhân ñược coi là ñứng yên (vì khối lượng hạt nhân rất lớn so với electron). Một ví dụ ñơn giản nhất của loại mô hình này là hệ hạt nhân nguyên tử heli: Hình I.2. Mô hình hệ nguyên tử heli PDF Created with deskPDF PDF Writer - Trial :: 9 1.2.2. Phương trình Schrodinger và sơ lược về cách giải Toán tử Hamilton mô tả trạng thái của hệ có dạng: ij n n n n i i i i i je i Ze e Hˆ m r r 2 2 2 2 =1 =1 < = − ∇ − + 2 ∑ ∑ ∑∑ ℏ Số hạng cuối cùng trong biểu thức trên biểu thị thế năng tương tác tĩnh ñiện giữa các electron. Vì các electron là không thể phân biệt ñược nên số hạng này chỉ có thể xác ñịnh một cách gần ñúng và phương trình Schrodinger do ñó cũng chỉ có thể giải một cách gần ñúng. Theo nguyên lí phản ñối xứng, hàm sóng mô tả trạng thái của các fermion có spin bán nguyên (1/2) phải là hàm phản ñối xứng. Khái niệm phản ñối xứng ở ñây có nghĩa là hàm phải ñổi dấu khi ta hoán vị 2 hạt bất kì trong hệ. ðối với hệ có số chẵn electron (2n, n: nguyên dương) thì ñiều kiện phản ñối xứng có thể thực hiện ñược bằng cách sử dụng hàm sóng dạng ñịnh thức Slater: ( ) ( )... ( ) ( ) ( )... ( ) .................................. ( ) ( )... ( ) n n n n n n x x x x x x x x x 1 1 2 1 2 1 1 2 2 2 2 2 1 2 2 2 2 2 χ χ χ χ χ χ Ψ = χ χ χ Hàm sóng toàn phần hay hàm ASO ( ) i i xχ là tích của hàm sóng không gian kiểu hidro (s, p, d, f, …) thu ñược khi giải phương trình Schrodinger cho nguyên tử hidro và hàm spin ( ,α β ) mô tả chuyển ñộng spin của electron. Phương trình Schrodinger của hệ có dạng: Hˆ EΨ = Ψ 1.2.3. Một số kết quả thu ñược Khi giải phương trình Schrodinger cho hệ nhiều electron thu ñược hai kết quả là hàm riêng là hàm sóng và trị riêng năng lượng. Hàm riêng là tích của hàm cầu và hàm bán kính. Phần hàm cầu giống nguyên tử một electron PDF Created with deskPDF PDF Writer - Trial :: 10 Hàm bán kính ñược xác ñịnh bằng phương pháp lặp trường tự hợp (SCF- self consistent field) và dạng ñịnh thức Slater. Tuy nhiên, công việc này vô cùng phức tạp nên ñể thực hiện cần phải có sự hỗ trợ của phần mềm và máy tính. Luận văn này sẽ không sử dụng phương pháp lặp trường tự hợp và các phần mềm hóa học tính toán cao cấp mà sử dụng những mô hình gần ñúng ñể giải quyết vấn ñề trên. 1.2.3.1. Sự gần ñúng Slater Các electron là những hạt mang ñiện tích âm nên khi chuyển ñộng, chúng sẽ che chắn lẫn nhau khỏi lực hút của hạt nhân nguyên tử. Khi ñó năng lượng của hệ sẽ ñược tính như sau: * n,l * * (Z ) (Z b) E . . (n ) (n ) 2 2 2 2 − = −13 6 = −13 6 b: hằng số chắn n: số lượng tử chính n*: số lượng tử chính hiệu dụng * n ... n . . ... = 1 2 3 4 5 6 = 1 2 3 3 7 4 4 2 Z: số ñiện tích hạt nhân Z*: số ñiện tích hạt nhân hiệu dụng l: số lượng tử phụ ðể tính hằng số chắn, các hàm AO ñược chia thành các nhóm như sau: s/ s p/ s p/ d/ s b/ d/ f /...1 2 2 3 3 3 4 4 4 4 Trị số hằng số chắn ñối với 1 electron ñang xét sẽ bằng tổng các trị số góp của các electron khác. Mỗi electron ở nhóm AO ngoài nhóm AO ñang xét không ñóng góp vào hằng số chắn. Mỗi electron nằm trên cùng một AO (nhóm AO) ñang xét ñóng góp vào hằng số chắn 1 lượng 0.35, riêng 1 electron trên AO-1s chỉ ñóng góp 0.3. Mỗi electron nằm bên trong nhóm AO ñang xét: Ở lớp n có trị số nhỏ hơn lớp ñang xét 1 ñơn vị, ñóng góp 0.85 PDF Created with deskPDF PDF Writer - Trial :: 11 Ở lớp n có trị số nhỏ hơn lớp ñang xét từ 2 ñơn vị trở lên, ñóng góp 1 Nếu nhóm AO ñang xét là AO-d hoặc AO-f thì mỗi electron ở AO trong góp 1 1.2.3.2. Hàm sóng Hàm sóng mô tả chuyển ñộng của electron trong hệ nguyên tử nhiều electron tốt nhất là những hàm thu ñược sau khi thực hiện phép giải lặp trường tự hợp. Thật may mắn, việc sử dụng hàm sóng dạng ñịnh thức Slater mô tả chuyển ñộng của các electron ñã ñưa ñến khả năng sử dụng các hàm sóng kiểu hidro ñể mô tả chuyển ñộng của các electron trong nguyên tử nhiều electron. Tuy nhiên, việc có thêm tương tác giữa các electron với nhau cũng như tương tác giữa hạt nhân và các electron ñã làm thay ñổi thứ tự các mức năng lượng của các hàm sóng kiểu hidro. Vấn ñề này sẽ ñược khảo sát chi tiết ở phần sau. 1.3. Cấu hình electron 1.3.1. Khái niệm Cấu hình electron của nguyên tử nói chung là sơ ñồ biểu thị sự phân bố các electron theo số lượng tử chính và số lượng tử phụ (n và l) hay theo lớp và phân lớp electron. Một số cấu hình electron tiêu biểu như sau: O : s s p hay [He] s p Na : s s p s hay [Ne] s 2 2 4 2 4 2 2 6 1 1 1 2 2 2 2 1 2 2 3 3 Sự sắp xếp các electron vào nguyên tử là một trong những vấn ñề quan trọng khi xét nguyên tử có nhiều electron. Kết quả của sự sắp xếp ñó ñược biểu diễn một cách khái quát bằng cấu hình electron và cơ sở ñể viết cấu hình electron của nguyên tử nhiều electron là: 1.3.2. Những cơ sở ñể viết cấu hình electron 1.3.2.1. Nguyên lý vững bền hay nguyên lý năng lượng cực tiểu Trạng thái hệ lượng tử có năng lượng thấp nhất hay cực tiểu là trạng thái cơ bản, ñó cũng là trạng thái bền vững nhất của hệ. Nguyên lý trên thể hiện một quy luật của thế giới tự nhiên là luôn luôn có xu hướng ñạt tới sự bền vững nhất. Sự sắp xếp các electron vào các AO trong nguyên tử cũng không nằm ngoài qui luật này. Theo ñó, trong nguyên tử, electron chiếm mức năng lượng thấp trước, tiếp ñến các PDF Created with deskPDF PDF Writer - Trial :: 12 mức năng lượng cao hơn. Trạng thái hệ có năng lượng thấp nhất là trạng thái cơ bản. 1.3.2.2. Quy tắc Klechkovxki Việc xuất hiện tương tác giữa các electron trong nguyên tử nhiều electron cũng như tương tác giữa hạt nhân và các electron ñã làm cho năng lượng của các obitan nguyên tử kiểu hidro không còn như trong nguyên tử hidro nữa. Một qui tắc kinh nghiệm ñơn giản nhất mô tả sự thay ñổi này chính là qui tắc Klechkovxki. Nội dung của qui tắc này như sau: Năng lượng của phân mức ,n lε tăng dần theo sự tăng của tổng trị số (n + l), nếu hai phân mức có cùng trị của tổng (n + l) thì ,n lε tăng theo sự tăng của n. Với n là số lượng tử chính, l là số lượng tử phụ, có thể diễn ñạt nội dung của quy tắc ñó bằng sơ ñồ sau: Hình I.3. Biểu diễn qui tắc Klechkovxki 1.3.2.3. Nguyên lí Pauli Nội dung của nguyên li Pauli như sau: “Trong nguyên tử, không thể tồn tại 2 electron giống nhau ở cả 4 số lượng tử n, l, ml,ms.” PDF Created with deskPDF PDF Writer - Trial :: 13 Ta ñã biết rằng AO là hàm không gian (r)Ψ  ñược xác ñịnh bằng một bộ ba số lượng tử n, l, ml (nên ta viết ln,l,m (r)Ψ  ). Khi kể ñến spin electron, ta có hàm sóng toàn phần mô tả ñầy ñủ trạng thái một electron là: l l s l s l nlm s n,l,m ,m n,l,m m nlm s (r). khi m / (r, ) (r). ( ) (r). khi m / ψ β = −1 2 Ψ σ = ψ η σ =  ψ α = +1 2     Như vậy mỗi electron sẽ ñược biểu thị bằng một bộ 4 số lượng tử n, l, ml, ms. Một cách phát biểu khác của nguyên lí Pauli là: “Trong nguyên tử, mỗi AO chỉ có thể bị chiếm bởi tối ña 2 electron” 1.3.2.4. Qui tắc Hund Nội dung cơ bản của qui tắc Hund ñược trình bày như sau: “Các electron phân bố vào các phân lớp sao cho số electron ñộc thân là tối ña và chúng phải cho spin song song” Việc áp dụng chi tiết các nguyên lí, qui tắc trên ñể viết cấu hình electron của nguyên tử nhiều electron ñược trình bày chi tiết trong tài liệu [1] và ở ñây không trình bày lại. 1.4. Cấu hình eletron và qui tắc gần ñúng Slater Việc viết cấu hình electron cho những nguyên tử họ s, p rất ñơn giản. Ở ñây, vấn ñề ñặt ra là ñối với các nguyên tử họ d, f từ chu kì 3 trở ñi, cấu hình electron phải ñược viết như thế nào cho hợp lí. Hiện nay, ñối với các nguyên tố họ d, có 2 cách viết ñang ñược sử dụng là: Cách : ...nd (n )s ( ) trong Ni(Z ) : [Ar] d s Cách : ...(n )s nd ( ) trong Ni(Z ) : [Ar] s d α 2 8 2 2 α 2 8 1 +1 α = 1÷10 = 28 3 4 2 +1 α = 1÷10 = 28 4 3 Hơn nữa, ở một số nguyên tố như Cu, Ag, Cr, Mo cấu hình electron có thể ñược viết như sau: Cách : ...nd (n )s ( ) trong Cu(Z ) : [Ar] d s Cách : ...(n )s nd ( ) trong Cu : [Ar] s d Cách : ...nd (n )s ( , ) trong Cu(Z ) : [Ar] d s Cách : ...(n )s nd ( , ) t α 2 9 2 2 α 2 9 β 1 10 1 1 β 1 +1 α = 1÷10 = 29 3 4 2 +1 α = 1÷10 4 3 3 +1 β = 5 10 = 29 3 4 4 +1 β = 5 10 rong Cu(Z ) : [Ar] s d1 10= 29 4 3 PDF Created with deskPDF PDF Writer - Trial :: 14 Tất nhiên, mỗi cách ñều có cơ sở riêng của nó, ở ñây có thể ñề xuất một tiêu chuẩn là cấu hình nào có năng lượng thấp nhất sẽ ứng với cấu hình electron hợp lí của nguyên tử ở trạng thái cơ bản. Tiêu chuẩn này cũng nghiệm ñúng với nhiều hiện tượng trong tự nhiên mà ở ñó năng lượng sẽ là thước ño của sự vận ñộng vật chất. Một thí dụ ñiển hình là hòn bi ñặt trên ñỉnh ñồi trọc luôn có khả năng chuyển ñộng xuống dốc ñồi dù chỉ bằng một cơn gió thoảng qua. Trên cơ sở của hóa học lượng tử tôi sẽ tiến hành tính tổng năng lượng electron cho từng cấu hình electron. Mức ñộ lí thuyết ñược dùng ở ñây là qui tắc gần ñúng Slater trong mô hình tính năng lượng của nguyên tử nhiều electron. Ở ñây cần lưu ý rằng ñây là một qui tắc gần ñúng rất thô nên những kết quả tính toán mà nó mang lại chỉ mang tính chất ñịnh tính. Việc tính toán ở ñây hoàn toàn có thể thực hiện bằng máy tính cá nhân (calculator), giấy và viết, tất nhiên công việc này sẽ không dễ dàng lắm vì phải thao tác với nhiều con số và nhiều phép toán phức tạp. ðể thuận tiện, việc tính toán sẽ ñược tiến hành trên máy tính thông qua phần mềm Mathematica. Vấn ñề còn lại là viết chương trình tính toán trên máy tính. Mô hình tính toán chi tiết cũng như chương trình tính sẽ ñược trình bày chi tiết trong chương sau. 1.5. Năng lượng ion hóa và qui tắc gần ñúng Slater Nói một cách ñơn giản, năng lượng ion hóa là năng lượng cần cung cấp ñể bứt một electron ra khỏi nguyên tử ở trạng thái khí. Tất nhiên, ñịnh nghĩa này ñang ñề cập ñến năng lượng ion hóa thứ nhất: M I M e+1+ → + Về nguyên tắc, nguyên tử có bao nhiêu electron thì ta có thể tách bấy nhiêu electron ra khỏi nguyên tử ñó. Khi ñó nguyên tử giờ ñây chỉ còn lại một hạt nhân trần trụi (bare nuclear). Như vậy, ta sẽ có các trị năng lượng ion hóa tiếp theo: M I M e+ 2+2+ → + M I M e2+ 3+3+ → + M I M e3+ 4+4+ → + PDF Created with deskPDF PDF Writer - Trial :: 15 Năng lượng ion hóa là một ñại lượng liên quan mật thiết với cấu hình electron của nguyên tử. Việc tách một electron ra khỏi nguyên tử luôn ñòi phải cung cấp năng lượng cho nguyên tử. Trên cơ sở qui tắc gần ñúng Slater và cấu hình electron của nguyên tử, năng lượng ion hóa thứ nhất, thứ hai, thứ 3, ... của các nguyên tử sẽ ñược tính toán. Việc tính năng lượng ion hóa cho các nguyên tử là một trong những vấn ñề còn ñang tranh cãi hiện nay. ðối với các nguyên tử thuộc chu kì nhỏ thì chúng ta có thể tính toán bình thường. Nhưng các nguyên tử từ chu kì ba trở lên khi tính năng lượng ion hóa sẽ gặp trở ngại liên quan ñến cấu hình elcetron, có nghĩa là sự phân bố electron trên obitan. ðể ñơn giản cho việc tính toán, tôi sẽ tiến hành khảo sát các nguyên tử chu kì 1 và 2 từ He (Z=2) ñến Ne (Z=10). Ở ñây không ñề cập ñến việc tính năng lượng ion của nguyên tử H (Z=1) vì công việc này quá ñơn giản. Mô hình tính toán chi tiết cũng như chương trình tính sẽ ñược ñề cập trong chương sau. Kết quả tính ñược cũng sẽ ñược so sánh với các giá trị năng lượng ion hóa thu ñược từ thực nghiệm ñể ñánh giá mức ñộ gần ñúng của qui tắc Slater. PDF Created with deskPDF PDF Writer - Trial :: 16 Chương 2: LẬP TRÌNH TÍNH TOÁN 2.1. Tính tổng năng lượng electron của nguyên tử 2.1.1. Mô hình tính ðể tính ñược tổng năng lượng electron của nguyên tử nhiều electron, phải thực hiện nhiều công ñoạn. Ở ñây các giai ñoạn sẽ ñược thực hiện cho nguyên tử Ni. Viết cấu hình electron cho nguyên tử Giả sử cấu hình electron của Ni là: Ni(Z ) : s s p s p d s2 2 6 2 6 8 2= 28 1 2 2 3 3 3 4 Tính hằng số chắn Các hằng số chắn ñược tính như sau: sb . .1 = 1×0 3 = 0 3 s pb . . .2 2 = 2×0 85+7×0 35 = 4 15 s pb . . .3 3 = 2×1+ 8×0 85+ 7×0 35 = 11 25 db . .3 = 18×1+7×0 35 = 20 45 sb . . .4 = 10×1+16×0 85+1×0 35 = 23 95 Tính năng lượng của 1 electron trên từng nhóm phân lớp Năng lượng của các phân mức ñược tính bằng công thức sau: n,l n,l Z b E . n 2 − = −13 6    Kết quả cụ thể như sau: s . E . . eV 2 1 28−0 3= −13 6 = −10435 10  1 s p . E . . eV 2 2 2  28−4 15= −13 6 = −1934 00  2 s p . E . . eV 2 3 3 28−11 25= −13 6 = −423 96  3 d . E . . eV 2 3  28−20 45= −13 6 = −86 14  3 PDF Created with deskPDF PDF Writer - Trial :: 17 s . E . . eV . 2 4 28−23 95= −13 6 = −16 30  3 7 Tính tổng năng lượng electron của nguyên tử s s p s p d sE E E E E E E ( . ) ( . ) ( . ) ( . ) ( . ) E . eV 1 2 2 3 3 3 4= 2× + 8× + 8× + 8× +2× = 2 −10435 10 + 8 −1934 00 + 8 −423 96 + 8 −86 14 +2 −16 30 = −40455 60 Như vậy qua các giai ñoạn trên, ñã tính ñược tổng năng lượng electron cho nguyên tử Ni với cấu hình Ni(Z ) : s s p s p d s2 2 6 2 6 8 2= 28 1 2 2 3 3 3 4 . Bằng cách thay ñổi các cấu hình khác nhau ta sẽ thu ñược các giá trị năng lượng electron khác nhau của cùng một nguyên tử ứng với các cấu hình electron khác nhau. So sánh các trị năng lượng tìm ñược sẽ tìm ra cấu hình electron hợp lí. 2.1.2. Chương trình tính Trên cơ sở mô hình tính toán ñã trình bày ở trên, tôi sẽ viết chương trình tính tổng năng lượng cho các cấu hình electron khác nhau của các nguyên tử khác nhau. Ở ñây, chương trình chỉ áp dụng cho các nguyên tử họ d chu kì 3 từ Sc (Z=21) ñến Zn (Zn=30). ðể tăng mức ñộ khái quát cho chương trình tính, tôi ñã ñưa các biến ,α β (số eletron trên phân lớp d) vào chương trình. Khi chạy chương trình, chỉ cần nhập các biến α hoặc β là có thể thu ñược kết quả tổng năng lượng electron cho các cấu hình electron khác nhau của cùng một nguyên tử. 2.1.2.1. Chương trình tính tổng năng lượng electron cho cấu hình 1 và 2 Cấu hình electron 1 và 2 ñã ñược xác ñịnh trong chương trước như sau: Cách : ...nd (n )s ( ) trong Ni(Z ) : [Ar] d s Cách : ...(n )s nd ( ) trong Ni(Z ) : [Ar] s d α 2 8 2 2 α 2 8 1 +1 α = 1÷10 = 28 3 4 2 +1 α = 1÷10 = 28 4 3 Chương trình tính tổng năng lượng electron cho hai cấu hình trên ñã ñược viết trên cơ sở phần mềm Mathematica như sau: "Chuong trinh tinh tong nang luong electron cho cau hinh 1 va 2" "So electron tren phan lop 3d" α=? PDF Created with deskPDF PDF Writer - Trial :: 18 "So dien tich hat nhan" Z=20+α "Cau hinh 2: 1s2 2s2 2p6 3s2 3p6 4s2 3dα:" "Gia tri cac hang so chan:" "b1s=" b1s=0.3 "b2s2p=" b2s2p=(7*0.35)+(2*0.85) "b3s3p=" b3s3p=(7*0.35)+(8*0.85)+(2*1) "b4s=" b4s=(1*0.35)+(8*0.85)+(10*1) "b3d=" b3d=((α-1)*0.35)+(20*1) "Nang luong cac nhom phan lop" "E1s=eV" E1s=-13.6*(Z-b1s)2/12 "E2s2p=eV" E2s2p=-13.6*(Z-b2s2p)2/22 "E3s3p=eV" E3s3p=-13.6*(Z-b3s3p)2/32 "E4s=eV" E4s=-13.6*(Z-b4s)2/(3.7)2 "E3d=eV" E3d=-13.6*(Z-b3d)2/(3)2 "Tong nang luong electron cua cau hinh 2:" "E2=eV" EI=(2*E1s)+(8*E2s2p)+(8*E3s3p)+(α*E3d)+(2*E4s) "Cau hinh I: 1s2 2s2 2p6 3s2 3p6 3dα 4s2:" "Gia tri cac hang so chan:" "b3d'=" b3d'=((α-1)*0.35)+(18*1) PDF Created with deskPDF PDF Writer - Trial :: 19 "b4s'" b4s'=(1*0.35)+((8+α)*0.85)+(10*1) "Nang luong cac nhom phan lop" "E4s'=eV" E4s'=-13.6*(Z-b4s')2/(3.7)2 "E3d'=eV" E3d'=-13.6*(Z-b3d')2/(3)2 "Tong nang luong electron cua cau hinh 1:" "E1=eV" EII=(2*E1s)+(8*E2s2p)+(8*E3s3p)+(α*E3d')+(2*E4s') ðể kiểm tra khả năng thực thi của chương trình, tôi ñã tiến hành tính thử cho Ni (Z=28) với α bằng 8. Kết quả thu ñược như sau: Out[720]= Chuong trinh tinh tong nang luong electron cho cau hinh 1 va 2 Out[721]= So electron tren phan lop 3d Out[722]= 8 Out[723]= So dien tich hat nhan Out[724]= 28 Out[725]= Cau hinh 2: 1s2 2s2 2p6 3s2 3p6 4s2 3d\[Alpha]: Out[726]= Gia tri cac hang so chan: Out[727]= b1s= Out[728]= 0.3 Out[729]= b2s2p= Out[730]= 4.15 Out[731]= b3s3p= Out[732]= 11.25 Out[733]= b4s= Out[734]= 17.15 Out[735]= b3d= PDF Created with deskPDF PDF Writer - Trial :: 20 Out[736]= 22.45 Out[737]= Nang luong cac nhom phan lop Out[738]= E1s=eV Out[739]= -10435.1 Out[740]= E2s2p=eV Out[741]= -1934. Out[742]= E3s3p=eV Out[743]= -423.961 Out[744]= E4s=eV Out[745]= -116.949 Out[746]= E3d=eV Out[747]= -46.546 Out[748]= Tong nang luong electron cua cau hinh 2: Out[749]= E2=eV Out[750]= -40340.2 Out[751]= Cau hinh I: 1s2 2s2 2p6 3s2 3p6 3d\[Alpha] 4s2: Out[752]= Gia tri cac hang so chan: Out[753]= b3d'= Out[754]= 20.45 Out[755]= b4s' Out[756]= 23.95 Out[757]= Nang luong cac nhom phan lop Out[758]= E4s'=eV Out[759]= -16.2947 Out[760]= E3d'=eV Out[761]= -86.1371 Out[762]= Tong nang luong electron cua cau hinh 1: Out[763]= E1=eV Out[764]= -40455.6 PDF Created with deskPDF PDF Writer - Trial :: 21 Giá trị -40455.6 là ñúng với giá trị tính bằng tay trong phần mô hình tính ở trên. Chương trình ñã chạy tốt nên sẽ ñược sử dụng ñể tính toán tổng năng lượng electron của các nguyên tử họ d thuộc chu kì 3 từ Sc ñến Zn. 2.1.2.2. Chương trình tính tổng năng lượng electron cho cấu hình 3 và 4 Cấu hình electron 1 và 2 ñã ñược xác ñịnh trong chương trước như sau: Cách : ...nd (n )s ( , ) trong Cu(Z ) : [Ar] d s Cách : ...(n )s nd ( , ) trong Cu(Z ) : [Ar] s d β 1 10 1 1 β 1 10 3 +1 β = 5 10 = 29 3 4 4 +1 β = 5 10 = 29 4 3 Chương trình tính tổng năng lượng electron cho hai cấu hình trên ñã ñược viết trên cơ sở phần mềm Mathematica như sau: "Chuong trinh tinh tong nang luong electron cho cau hinh 3 va 4" "So electron tren phan lop 3d" β=? "So dien tich hat nhan" Z=20+β "Cau hinh 4: la 1s2 2s2 2p6 3s2 3p6 4s1 3dβ:" "Gia tri cac hang so chan:" "b1s=" b1s=0.3 "b2s2p=" b2s2p=(7*0.35)+(2*0.85) "b3s3p=" b3s3p=(7*0.35)+(8*0.85)+(2*1) "b4s=" b4s=(8*0.85)+(10*1) "b3d=" b3d=((β-1)*0.35)+(19*1) "Nang luong cac nhom phan lop" "E1s=eV" E1s=-13.6*(Z-b1s)2/12 "E2s2p=eV" PDF Created with deskPDF PDF Writer - Trial :: 22 E2s2p=-13.6*(Z-b2s2p)2/22 "E3s3p=eV" E3s3p=-13.6*(Z-b3s3p)2/32 "E4s=eV" E4s=-13.6*(Z-b4s)2/(3.7)2 "E3d=eV" E3d=-13.6*(Z-b3d)2/(3)2 "Tong nang luong electron cua cau hinh 4:" "E4=eV" EI=(2*E1s)+(8*E2s2p)+(8*E3s3p)+(β*E3d)+(1*E4s) "Cau hinh 3: 1s2 2s2 2p6 3s2 3p6 3dβ 4s1:" "Gia tri cac hang so chan:" "b3d'=" b3d'=((β-1)*0.35)+(18*1) "b4s'" b4s'=((8+β)*0.85)+(10*1) "Nang luong cac nhom phan lop" "E4s'=eV" E4s'=-13.6*(Z-b4s')2/(3.7)2 "E3d'=eV" E3d'=-13.6*(Z-b3d')2/(3)2 "Tong nang luong electron cua cau hinh 3:" "E3=eV" EII=(2*E1s)+(8*E2s2p)+(8*E3s3p)+(β*E3d')+(1*E4s') ðể kiểm tra khả năng thực thi của chương trình, tôi ñã tiến hành tính thử cho Cr(Z=25) với β bằng 5. Kết quả thu ñược như sau: Out[630]= Chuong trinh tinh tong nang luong electron cho cau hinh 3 va 4 Out[631]= So electron tren phan lop 3d Out[632]= 5 Out[633]= So dien tich hat nhan PDF Created with deskPDF PDF Writer - Trial :: 23 Out[634]= 25 Out[635]= Cau hinh 4: la 1s2 2s2 2p6 3s2 3p6 4s1 3d\[Beta]: Out[636]= Gia tri cac hang so chan: Out[637]= b1s= Out[638]= 0.3 Out[639]= b2s2p= Out[640]= 4.15 Out[641]= b3s3p= Out[642]= 11.25 Out[643]= b4s= Out[644]= 16.8 Out[645]= b3d= Out[646]= 20.4 Out[647]= Nang luong cac nhom phan lop Out[648]= E1s=eV Out[649]= -8297.22 Out[650]= E2s2p=eV Out[651]= -1478.06 Out[652]= E3s3p=eV Out[653]= -285.694 Out[654]= E4s=eV Out[655]= -66.798 Out[656]= E3d=eV Out[657]= -31.9751 Out[658]= Tong nang luong electron cua cau hinh 4: Out[659]= E4=eV Out[660]= -30931.1 Out[661]= Cau hinh 3: 1s2 2s2 2p6 3s2 3p6 3d\[Beta] 4s1: PDF Created with deskPDF PDF Writer - Trial :: 24 Out[662]= Gia tri cac hang so chan: Out[663]= b3d'= Out[664]= 19.4 Out[665]= b4s' Out[666]= 21.05 Out[667]= Nang luong cac nhom phan lop Out[668]= E4s'=eV Out[669]= -15.4999 Out[670]= E3d'=eV Out[671]= -47.3884 Out[672]= Tong nang luong electron cua cau hinh 3: Out[673]= E3=eV Out[674]= -30956.9 Chương trình ñã chạy tốt nên sẽ ñược sử dụng ñể tính toán tổng năng lượng electron của các nguyên tử họ d thuộc chu kì 3 từ Sc ñến Zn. 2.2. Tính năng lượng ion hóa của các nguyên tử 2.2.1. Mô hình tính ðể tính năng lượng ion hóa của các nguyên tử nhiều electron, công việc ñầu tiên là tính tổng lượng electron của các nguyên tử và ion. Vấn ñề này ñã ñược giải quyết trong phần mô hình tính tổng năng lượng electron. Khi ñó năng lượng ion hóa sẽ ñược tính như sau: I E(M ) E(M) I E(M ) E(M ) I E(M ) E(M ) + 1 2+ + 2 3+ 2+ 3 = − = − = − ðể ñơn giản cho việc lập trình, khóa luận sẽ tiến hành tính năng lượng ion hóa thứ nhất cho các nguyên tố chu kì 2 trong bảng hệ thống tuần hoàn. ðối với cacbon, tôi sẽ tiến hành tính tất cả các trị năng lượng có thể có trong nguyên tử này. Các kết quả sau ñó sẽ ñược so sánh với thực nghiệm. 2.2.2. Chương trình tính PDF Created with deskPDF PDF Writer - Trial :: 25 2.2.2.1. Tính năng lượng ion hóa thứ nhất ðể tăng mức ñộ khái quát cho chương trình, tôi chỉ ñưa biến Z (số ñiện tích hạt nhân) và chương trình là có thể thực thi ñược. Nội dung của chương trình năng lượng ion hóa cho He(Z=2) như sau: "TINH NANG LUONG ION HOA THU NHAT CHO He" "So dien tich hat nhan cua nguyen tu" Z=2 "Gia tri cac hang so chan cua nguyen tu" "b1s" b1s=0.3 "Cac muc nang luong electron cua nguyen tu" "E1s" E1s=-13.6 (Z-b1s)2/12 "Tong nang luong electron cua nguyen tu" E0=2*E1s "Tong nang luong electron cua ion" E1=-13.6 (Z)2/12 "Nang luong ion hoa thu nhat" I1=E1-E0 Chương trình cũng ñã ñược chạy thử với kết quả như sau: Out[201]= TINH NANG LUONG ION HOA THU NHAT CHO He Out[202]= So dien tich hat nhan cua nguyen tu Out[203]= 2 Out[204]= Gia tri cac hang so chan cua nguyen tu Out[205]= b1s Out[206]= 0.3 Out[207]= Cac muc nang luong electron cua nguyen tu Out[208]= E1s Out[209]= -39.304 Out[210]= Tong nang luong electron cua nguyen tu Out[211]= -78.608 PDF Created with deskPDF PDF Writer - Trial :: 26 Out[212]= Tong nang luong electron cua ion Out[213]= -54.4 Out[214]= Nang luong ion hoa thu nhat Out[215]= 24.208 Chương trình tính năng lượng ion hóa cho các nguyên tố chu kì 2 như sau: "CHUONG TRINH TINH NANG LUONG ION HOA THU NHAT CHO CAC NGUYEN TO CHU KI 2" "So dien tich hat nhan cua nguyen tu" Z=6 "Gia tri cac hang so chan" "b2s2p" b2s2p=(Z-3)*0.35 +2*0.85 "b1s" b1s=0.3 "Cac muc nang luong electron cua nguyen tu" "E2s2p" E2s2p=-13.6 (Z-b2s2p)2/22 "E1s" E1s=-13.6 (Z-b1s)2/12 "Tong nang luong electron cua nguyen tu" E0=2*E1s+(Z-2)*(E2s2p) "Gia tri cac hang so chan cua ion" "b2s2p'" b2s2p'=(Z-4)*0.35 +2*0.85 "Cac muc nang luong electron cua ion" "E2s2p'" E2s2p'=-13.6 (Z-b2s2p')2/22 "Tong nang luong electron cua ion" E1=2*E1s+(Z-3)*(E2s2p') "Nang luong ion hoa thu nhat" I1=E1-E0 Chương trình cũng ñã ñược chạy thử ñể tính năng lượng ion hóa thứ nhất cho nguyên tử C, kết quả như sau: PDF Created with deskPDF PDF Writer - Trial :: 27 Out[765]= CHUONG TRINH TINH NANG LUONG ION HOA THU NHAT CHO CAC NGUYEN TO CHU KI 2 Out[766]= So dien tich hat nhan cua nguyen tu Out[767]= 6 Out[768]= Gia tri cac hang so chan Out[769]= b2s2p Out[770]= 2.75 Out[771]= b1s Out[772]= 0.3 Out[773]= Cac muc nang luong electron cua nguyen tu Out[774]= E2s2p Out[775]= -35.9125 Out[776]= E1s Out[777]= -441.864 Out[778]= Tong nang luong electron cua nguyen tu Out[779]= -1027.38 Out[780]= Gia tri cac hang so chan cua ion Out[781]= b2s2p' Out[782]= 2.4 Out[783]= Cac muc nang luong electron cua ion Out[784]= E2s2p' Out[785]= -44.064 Out[786]= Tong nang luong electron cua ion Out[787]= -1015.92 Out[788]= Nang luong ion hoa thu nhat Out[789]= 11.458 Như vậy, giá trị năng lượng ion hóa thứ nhất của nguyên tử C tính ñược là 11.458 eV. Chương trình ñã thực thi nên có thể sử dụng cho các nguyên tố còn lại trong chu kì 2 của bảng hệ thống tuần hoàn. 2.2.2.2 Tính các trị năng lượng ion hóa có thể có của C PDF Created with deskPDF PDF Writer - Trial :: 28 Việc tính tất cả các trị năng lượng ion hóa của một nguyên tử là công việc vô cùng khó khăn về mặt lập trình. Do ñó, công việc này chỉ có thể tiến hành riêng rẽ cho từng nguyên tử riêng biệt. Chương trình sau ñây ñã ñược viết ñể tính tất cả các giá trị năng lượng ion hóa có thể có của nguyên tử C (Z=6): "TINH NANG LUONG ION HOA CUA CACBON" "So dien tich hat nhan cua nguyen tu" Z=6 "Tinh tong nang luong cua nguyen tu" "Cau hinh electron cua nguyen tu: 1s2 2s2 2p2" "Gia tri cac hang so chan cua nguyen tu" "b2s2p" b2s2p=(Z-3)*0.35 +2*0.85 "b1s" b1s=0.3 "Cac muc nang luong electron cua nguyen tu" "E2s2p" E2s2p=-13.6 (Z-b2s2p)2/22 "E1s" E1s=-13.6 (Z-b1s)2/12 "Tong nang luong electron cua nguyen tu" E0=2*E1s+(Z-2)*(E2s2p) "Tinh tong nang luong cua ion 1+" "Cau hinh electron cua ion 1+: 1s2 2s2 2p1" "Gia tri cac hang so chan cua ion 1+" "b2s2p'" b2s2p'=(Z-4)*0.35 +2*0.85 "Cac muc nang luong electron cua ion 1+" "E2s2p'" E2s2p'=-13.6 (Z-b2s2p')2/22 "Tong nang luong electron cua ion 1+" PDF Created with deskPDF PDF Writer - Trial :: 29 E1=2*E1s+(Z-3)*(E2s2p') "Nang luong ion hoa thu nhat" I1=E1-E0 "Tinh tong nang luong electron cua ion 2+" "Cau hinh electron cua ion 2+: 1s2 2s2 2p0" "Gia tri cac hang so chan cua ion 2+" "b2s2p''" b2s2p''=(Z-5)*0.35 +2*0.85 "Cac muc nang luong electron cua ion 2+" "E2s2p''" E2s2p''=-13.6 (Z-b2s2p'')2/22 "Tong nang luong electron cua ion 2+" E2=2*E1s+(Z-4)*(E2s2p'') "Nang luong ion hoa thu hai" I2=E2-E1 "Tinh tong nang luong electron cua ion 3+" "Cau hinh electron cua ion 3+: 1s2 2s1 2p0" "Gia tri cac hang so chan cua ion 3+" "b2s2p'''" b2s2p'''=(Z-6)*0.35 +2*0.85 "Cac muc nang luong electron cua ion 3+" "E2s2p'''" E2s2p'''=-13.6 (Z-b2s2p''')2/22 "Tong nang luong electron cua ion 3+" E3=2*E1s+(Z-5)*(E2s2p''') "Nang luong ion hoa thu ba" I3=E3-E2 "Tinh tong nang luong cua ion 4+" "Cau hinh electron cua ion 4+: 1s2 2s0 2p0" "Tong nang luong electron cua nguyen tu" PDF Created with deskPDF PDF Writer - Trial :: 30 E4=2*E1s "Nang luong ion hoa thu tu" I4=E4-E3 "Tinh tong nang luong cua ion 5+" "Cau hinh electron cua ion 4+: 1s1 2s0 2p0" "Tong nang luong electron cua nguyen tu" E5=-13.6 (Z)2/12 "Nang luong ion hoa thu nam" I5=E5-E4 "Nang luong ion hoa thu 6" I6=-E5 Kết quả thu ñược là các trị năng lượng ion hóa của C từ I1 ñến I6: Out[613]= TINH NANG LUONG ION HOA CUA CACBON Out[614]= So dien tich hat nhan cua nguyen tu Out[615]= 6 Out[616]= Tinh tong nang luong cua nguyen tu Out[617]= Cau hinh electron cua nguyen tu: 1s2 2s2 2p2 Out[618]= Gia tri cac hang so chan cua nguyen tu Out[619]= b2s2p Out[620]= 2.75 Out[621]= b1s Out[622]= 0.3 Out[623]= Cac muc nang luong electron cua nguyen tu Out[624]= E2s2p Out[625]= -35.9125 Out[626]= E1s Out[627]= -441.864 Out[628]= Tong nang luong electron cua nguyen tu Out[629]= -1027.38 PDF Created with deskPDF PDF Writer - Trial :: 31 Out[630]= Tinh tong nang luong cua ion 1+ Out[631]= Cau hinh electron cua ion 1+: 1s2 2s2 2p1 Out[632]= Gia tri cac hang so chan cua ion 1+ Out[633]= b2s2p' Out[634]= 2.4 Out[635]= Cac muc nang luong electron cua ion 1+ Out[636]= E2s2p' Out[637]= -44.064 Out[638]= Tong nang luong electron cua ion 1+ Out[639]= -1015.92 Out[640]= Nang luong ion hoa thu nhat Out[641]= 11.458 Out[642]= Tinh tong nang luong electron cua ion 2+ Out[643]= Cau hinh electron cua ion 2+: 1s2 2s2 2p0 Out[644]= Gia tri cac hang so chan cua ion 2+ Out[645]= b2s2p'' Out[646]= 2.05 Out[647]= Cac muc nang luong electron cua ion 2+ Out[648]= E2s2p'' Out[649]= -53.0485 Out[650]= Tong nang luong electron cua ion 2+ Out[651]= -989.825 Out[652]= Nang luong ion hoa thu hai Out[653]= 26.095 Out[654]= Tinh tong nang luong electron cua ion 3+ Out[655]= Cau hinh electron cua ion 3+: 1s2 2s1 2p0 Out[656]= Gia tri cac hang so chan cua ion 3+ Out[657]= b2s2p''' Out[658]= 1.7 Out[659]= Cac muc nang luong electron cua ion 3+ PDF Created with deskPDF PDF Writer - Trial :: 32 Out[660]= E2s2p''' Out[661]= -62.866 Out[662]= Tong nang luong electron cua ion 3+ Out[663]= -946.594 Out[664]= Nang luong ion hoa thu ba Out[665]= 43.231 Out[666]= Tinh tong nang luong cua ion 4+ Out[667]= Cau hinh electron cua ion 4+: 1s2 2s0 2p0 Out[668]= Tong nang luong electron cua nguyen tu Out[669]= -883.728 Out[670]= Nang luong ion hoa thu tu Out[671]= 62.866 Out[672]= Tinh tong nang luong cua ion 5+ Out[673]= Cau hinh electron cua ion 4+: 1s1 2s0 2p0 Out[674]= Tong nang luong electron cua nguyen tu Out[675]= -489.6 Out[676]= Nang luong ion hoa thu nam Out[677]= 394.128 Out[678]= Nang luong ion hoa thu 6 Out[679]= 489.6 Các kết quả thu ñược khi chạy chương trình tính sẽ ñược tổng kết và ñánh giá trong chương sau. Chương 3: KẾT QUẢ VÀ THẢO LUẬN 3.1.Cấu hình electron của các nguyên từ họ d chu kì 3 3.1.1. Các nguyên tố họ d chu kì 3 (ngoại trừ Cr và Cu) PDF Created with deskPDF PDF Writer - Trial :: 33 Các kết quả thu ñược sau khi thực hiện phép tính tổng năng lượng electron của nguyên tử ñược tổng kết như sau: Nguyên tử Sc Ti V Mn ...nd (n )sα 2+1 -20558.3 -22931.7 -25455.1 -30967.1 Năng lượng ...(n )s nd2 α+1 -20557.8 -22926.7 -25441.7 -30924.8 Nguyên tử Fe Co Ni Zn ...nd (n )sα 2+1 -33963.4 -37125.0 -40455.6 -47639.6 Năng lượng ...(n )s nd2 α+1 -33900.7 -37037.9 -40340.2 -47455.6 Phân tích các số liệu từ 2 bảng trên cho thấy ở cả 8 nguyên tố này, cách viết cấu hình electron theo cách 1 ...nd (n )sα 2+1 sẽ có năng lượng thấp hơn cách viết 2 ...(n )s nd2 α+1 nên cấu hình electron viết theo cách 1 sẽ hợp lí hơn. 3.1.2. Các nguyên tố Cr và Cu Các kết quả thu ñược sau khi thực hiện phép tính tổng năng lượng electron của nguyên tử ñược tổng kết như sau: Cấu hình ...nd (n )sα 2+1 ...(n )s nd2 α+1 nd (n )sβ 1+1 ...(n )s nd1 β+1 Cr -28132.3 -28106.4 -30956.9 -30931.1 Năng lượng Cu -43959.2 -43811.5 -47624.0 -47522.8 Từ các số liệu trên có thể thấy rằng chỉ có cách viết số 3 là có năng lượng thấp nhất nên viết theo cách này là hợp lí hơn cả. 3.2. Năng lượng ion hóa thứ nhất của các nguyên tố chu kì 2 Việc thực thi chương trình tính năng lượng ion hóa cho các nguyên tố khác nhau của chu kì 2 thu ñược kết quả sau: PDF Created with deskPDF PDF Writer - Trial :: 34 Nguyên tử Li Be B C N O F Ne Tính ñược 5.75 7.87 9.78 11.46 12.92 14.16 15.18 15.98 I1(eV) Thực nghiệm 5.39 9.32 8.30 11.26 14.53 13.62 17.42 21.56 Các dữ liệu trong bảng trên cho thấy ngoại trừ Ne có sai số giữa trị năng lượng ion hóa giữa lí thuyết tính ñược và thực nghiệm là rất lớn (5.58eV), các trường hợp còn lại sai biệt không quá lớn và có thể chấp nhận ñược. 3.3. Các trị năng lượng ion hóa có thể có của C Các kết quả thu ñược từ việc chạy chương trình tính các trị năng lượng ion hóa có thể có của nguyên tử C ñược tổng kết trong bảng sau: Năng lượng ion hóa (eV) I1 I2 I3 I4 I5 I6 Lí thuyết 11.46 26.10 43.23 62.87 394.13 489.60 Thực nghiệm 11.26 24.39 47.89 64.50 392.14 490.05 Các dữ liệu bảng trên cho thấy giá trị năng lượng ion hóa giữa lí thuyết tính ñược và thực nghiệm sai khác nhau không quá lớn có thể chấp nhận ñược. Các trị năng lượng tính ñược ở trên cho thấy năng lượng ion hóa tăng rất nhanh theo số electron bị tách ra khỏi nguyên tử. ðiều này là hoàn toàn dễ hiểu vì tách một electron ra khỏi một cation bao giờ cũng khó hơn việc tách 1 electron ra khỏi một nguyên tử trung hòa ñiện tích. 3.4. Một số nhận xét Qua các kết quả ñược tổng hợp và phân tích ở trên, có thể rút ra một số nhận xét sau ñây: Về cách viết cấu hình electron cho các nguyên tử họ d chu kì 3, có thể thấy rằng cách viết 1 ...nd (n )sα 2+1 và 3 nd (n )sβ 1+1 có tổng năng lượng electron là nhỏ nhất nên ñây là những cách viết cấu hình electron hợp lí nhất. ðiều này không có gì mới mẻ ñối với những nhà hóa học nhưng ở ñây tôi muốn mang ñến một cách tiếp cận khác về cách viết cấu hình electron cho nguyên tử nhiều electron. PDF Created with deskPDF PDF Writer - Trial :: 35 Về tính toán năng lượng ion hóa cho các nguyên tử nhiều electron ở chu kì 1 và 2 của bảng hệ thống tuần hoàn, các trị năng lượng ion hóa thứ nhất tính ñược cho thấy không có sai khác nhiều giữa các giá trị tính ñược từ lí thuyết và các giá trị thực nghiệm ngoại trừ một số ngoại lệ như Ne. Việc sai khác là tất nhiên nếu chúng ta lưu ý rằng sự gần ñúng Slater là một sự gần ñúng kinh nghiệm nên các tính toán mà nó mang lại chỉ mang tính chất ñịnh tính. Nếu cần các kết quả chính xác hơn, hãy sử dụng các phần mềm hóa học lượng tử thương mại hiện ñại dựa trên cơ sở hóa học lượng tử. Các trị năng lượng ion hóa có thể có của nguyên tử C cũng ñược tiến hành tính toán trên cơ sở phần mềm Mathematica. Kết quả tính ñược cho thấy năng lượng ion hóa sẽ tăng dần khi số electron bị bứt ra khỏi nguyên tử càng lớn. Hóa học hiện ñại ngày nay không chỉ là một khoa học thực nghiệm mà nó còn có một cơ sở lí thuyết vững chắc. Hiểu rõ bản chất về cấu tạo electron của nguyên tử cũng như liên hệ giữa cấu tạo electron và các ñại lượng khác của hóa học như năng lượng ion hóa sẽ mang lại nhiều lợi ích cho việc giảng dạy hóa học sau này. Hơn nữa, qua khóa luận này ñã chứng tỏ ñược rằng Mathematica là một công cụ lập trình toán học không thể thiếu khi giải quyết các bài toán hóa học. Sự kết hợp hoàn hảo giữa phần mềm, máy tính, phòng thí nghiệm và bộ não con người ñã tạo nên ngành hóa học hiện ñại ngày nay. Ở nước ta, một xu hướng hiện ñại ngày nay là ứng dụng công nghệ thông tin vào dạy học ñang bùng phát mạnh mẽ trong nền giáo dục nước nhà. Hy vọng rằng với khóa luận này, tôi sẽ góp phần nhỏ vào phong trào này. KẾT LUẬN PDF Created with deskPDF PDF Writer - Trial :: 36 Qua thời gian tìm hiểu và thực hiện khóa luận này, tôi ñã thực hiện ñược những công việc sau ñây: 1. Tìm hiểu về lí thuyết hóa học lượng tử 2. Tìm hiểu về cách sử dụng phần mềm Mathematica phiên bản 6.0 3. Viết chương trình tính tổng năng lượng electron của nguyên tử bằng qui tắc gần ñúng Slater trên cơ sở phần mềm Mathematica ñể tìm ra cách viết cấu hình electron hợp lí của các nguyên tử họ d chu kì 3 và tính năng lượng ion hóa của các nguyên tố nhiều electron chu kì 1 và 2 Các kết quả tính ñược cho thấy ñối với các nguyên tố họ d chu kì 3, cách viết cấu hình electron hợp lí là ...nd (n )sα 2+1 chứ không phải là ...(n )s nd2 α+1 . Riêng ñối với Cr và Cu thì cách viết cấu hình electron hợp lí là nd (n )sβ 1+1 với ,β = 5 10 . Các trị năng lượng ion hóa thứ nhất của các nguyên tố nhiều electron chu kì 1 và 2 tính ñược không sai khác nhiều so với các trị năng lượng ion hóa thực nghiệm. ðây cũng là một thành công lớn của sự gần ñúng Slater nếu chú ý ñến mức ñộ gần ñúng ban ñầu của sự gần ñúng này. Các trị năng lượng ion hóa có thể có của C cũng ñã ñược tính toán. Tất nhiên, việc tính toán ở ñây sẽ không tuyệt ñối chính xác nhưng cũng ñã mang lại một ñánh giá ñịnh tính về qui luật thay ñổi của các giá trị năng lượng ion hóa có thể trong một nguyên tử. TÀI LIỆU THAM KHẢO PDF Created with deskPDF PDF Writer - Trial :: 37 1. Nguyễn ðức Chuy, (1998), Hóa Học ðại Cương, Nhà Xuất Bản Giáo Dục. 2. Nguyễn ðình Huề, Nguyễn ðức Chuy, (2003), Nhà Xuất Bản Giáo Dục. 3. Trần Thành Huế, (2004), Hóa Học ðại Cương 1 Cấu Tạo Chất, Nhà Xuất bản ðại Học Sư Phạm. 4. Trần Thành Huế, (2006), Tư Liệu Hóa Học 10, Nhà Xuất Bản Giáo Dục. 5. Lâm Ngọc Thềm, Phạm Văn Nhiêu, Lê Kim Long, Cơ Sở Hóa Học Lượng Tử, Nhà Xuất Bản Khoa Học Và Kĩ Thuật 6. Lê Xuân Trọng (Tổng Chủ biên kiêm Chủ Biên), Từ Ngọc Ánh – Lê Mậu Quyền – Phan Quang Thái (2006), Hóa Học 10 nâng cao Nhà Xuất Bản Giáo Dục. 7. Nguyễn Xuân Trường (Tổng Chủ biên kiêm Chủ Biên), Nguyễn ðức Chuy – Lê Mậu Quyền – Lê Xuân Trọng (2006), Hóa Học 10, Nhà Xuất Bản Giáo Dục. 8. PDF Created with deskPDF PDF Writer - Trial ::