Dao động các nguyên tử trong phân tử

Chia sẻ: Trần Bá Trung | Ngày: | Loại File: PDF | Số trang:17

Thêm vào BST

Báo xấu

231
lượt xem 18
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Các nguyên tử trong phân tử dao động như thế nào? Tóm tắt nội dung Sự dao động của phân tử hai nguyên tử rất giống với sự dao động điều hòa của con lắc lò xo cực nhỏ. Trong phần này, chúng ta sẽ giải phương trình Schroedinger cho hệ dao động điều hòa để tìm hàm sóng o và các mức năng lượng được phép, từ đó áp dụng vào phân tử. Đây là một phương trình vi phân khá phức tạp, thường được giải bằng phương pháp chuỗi lũy thừa. Vì vậy, trước hết, ta bàn về phương...

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Dao động các nguyên tử trong phân tử

Các nguyên tử trong phân tử dao động như thế nào? Lý Lê Ngày 6 tháng 8 năm 2009 Tóm tắt nội dung Sự dao động của phân tử hai nguyên tử rất giống với sự dao động điều hòa của con lắc lò xo cực nhỏ. Trong phần này, chúng ta sẽ giải phương trình Schr¨dinger cho hệ dao động điều hòa để tìm hàm sóng o và các mức năng lượng được phép, từ đó áp dụng vào phân tử. Đây là một phương trình vi phân khá phức tạp, thường được giải bằng phương pháp chuỗi lũy thừa. Vì vậy, trước hết, ta bàn về phương pháp chuỗi lũy thừa cho phương trình vi phân. 1 Nghiệm chuỗi lũy thừa của phương trình vi phân Cho đến thời điểm này, chúng ta chỉ mới xét đến những trường hợp mà hàm thế năng V (x) là hằng số; nghĩa là phương trình Schr¨dinger là một phương o trình vi phân thuần nhất tuyến tính bậc hai với hệ số không đổi. Tuy nhiên, thực tế ta sẽ gặp những trường hợp mà thế năng V thay đổi theo tọa độ. Khi đó, phương trình Schr¨dinger sẽ trở nên rất khó tìm nghiệm ở dạng tổ o hợp của các hàm số sơ cấp xác định. Điều này cũng xảy ra ngay cả khi các phương trình vi phân có dạng rất đơn giản. Chẳng hạn phương trình sau y − 3xy + 2y = 0 Đây là phương trình vi phân cấp hai, hệ số hàm nhưng ta không thể tìm được một nghiệm riêng dưới dạng hàm số sơ cấp như đã tiến hành cho hạt trong hộp một chiều. Một trong các phương pháp thông dụng để giải những phương trình vi phân dạng này là ứng dụng lí thuyết chuỗi để tìm nghiệm của phương trình dưới dạng chuỗi lũy thừa y = c0 + c1 x + c2 x2 + · · · + cn xn (1) Sau đây, chúng ta sẽ minh họa bằng cách giải một phương trình vi phân rất đơn giản như sau y (x) = y(x) (2) với điều kiện biên y(0) = 1. 1
Đây là một phương trình vi phân thuần nhất tuyến tính bậc nhất với hệ số không đổi. Ta có phương trình bổ trợ của (2) là s−1=0 ⇒s=1 Vậy nghiệm của (2) là y(x) = ex (3) Bây giờ, ta giải (2) bằng phương pháp chuỗi lũy thừa. Giả sử nghiệm của nó có dạng ∞ y(x) = c0 + c1 x + c2 x2 + c3 x3 + · · · + cn xn = cn xn (4) n=0 Lấy đạo hàm bậc nhất (4) ta được ∞ y (x) = c1 + 2c2 x + 3c3 x2 + · · · + ncn xn−1 = ncn xn−1 (5) n=1 Thế (4) và (5) vào (2), ta được c1 + 2c2 x + 3c3 x2 + · · · + ncn xn−1 = c0 + c1 x + c2 x2 + · · · + cn xn (6) Phương trình (6) đúng khi các hệ số ở hai vế với cùng lũy thừa x bằng nhau. Nghĩa là, ta có c1 = c0 2c2 = c1 3c3 = c2 . . . ncn = cn−1 Từ đó, ta có c1 = c0 c1 c0 c0 c0 c2 = = = = 2 2 1·2 2! c2 c0 c0 c0 c3 = = = = 3 6 1·2·3 3! . . . cn−1 c0 c0 cn = = = n 1 · 2 · 3···n n! Thế các hệ số tìm được ở trên vào phương trình y(x) = c0 + c1 x + c2 x2 + c3 x3 + · · · + cn xn 2
Ta được x2 x3 xn y(x) = c0 1 + x + + + ··· + 2! 3! n! Áp dụng điều kiện y(0) = 1, ta suy ra c0 = 1 và do đó nghiệm của (2) là x2 x3 xn y(x) = 1 + x + + + ··· + (7) 2! 3! n! Từ (3) và (7) ta thấy x2 x3 xn ex = 1 + x + + + ··· + (8) 2! 3! n! Đây chính là công thức khai triển Taylor cho hàm ex theo lũy thừa xn . Kết quả khai triển càng chính xác khi n càng lớn. 2 Dao động điều hòa trong không gian một chiều 2.1 Quan điểm của cơ học cổ điển 1 2 Một hạt khối lượng m chuyển động trong trường thế năng V (x) = kx sẽ 2 chịu một lực tác dụng được xác định như sau dV Fx = − = −kx (9) dx với k là hằng số lực. Theo định luật thứ hai Newton, F = ma, ta có dx2 (t) −kx = m (10) dt2 k với t là thời gian. Đặt ω = , phương trình (10), trở thành m x (t) + ω 2 x(t) = 0 (11) Phương trình bổ trợ của (11) là s2 + ω 2 = 0 ⇒ s2 = −ω 2 = i2 ω 2 ⇒ s = ±iω Như vậy, nghiệm của (11) là x(t) = c1 eiωt + c2 e−iωt (12) 3
Áp dụng phương trình dạng mũ của số phức, ta biến đổi (12) thành x(t) = (c1 + c2 ) cos(ωt) + i(c1 − c2 ) sin(ωt) hay x(t) = C1 cos(ωt) + C2 sin(ωt) (13) Mặt khác, ta có sin(a + b) = sin a cos b + cos a sin b Do đó, ta có thể viết lại (13) như sau x(t) = A sin(ωt + ϕ0 ) (14) Trong đó, A được gọi là biên độ dao động cực đại của x; ϕ0 là hệ số góc. Hai 2π hằng số này được xác định từ điều kiện ban đầu. Đại lượng T = được ω 1 gọi là chu kì dao động (Hz). Đại lượng ν = được gọi là tần số (s−1 ) T 1 ω 1 k ν= = = (15) T 2π 2π m Như vậy tần số dao động điều hòa tỉ lệ thuận với hằng số lực k và tỉ lệ nghịch với khối lượng m. Vận tốc dao động được xác định như sau dx(t) v(t) = = Aω cos(ωt + ϕ0 ) dt Năng lượng của hạt bằng tổng động năng và thế năng 1 1 E = T + V = mv 2 + kx2 (16) 2 2 1 1 E= mA2 ω 2 cos2 (ωt + ϕ0 ) + kA2 sin2 (ωt + ϕ0 ) (17) 2 2 Thế ω = k/m và áp dụng sin2 x + cos2 x = 1, ta được 1 1 E = kA2 [cos2 (ωt + ϕ0 ) + sin2 (ωt + ϕ0 )] = kA2 (18) 2 2 2.2 Quan điểm của cơ học lượng tử 2.2.1 Phương trình Schr¨dinger của dao động điều hòa o Các toán tử động năng và thế năng của dao động điều hòa trong không gian một chiều như sau 2 d2 1 2 T =− ; V = V (x) = kx (19) 2m dx2 2 4
Hamiltonian của dao động điều hòa là 2 d2 1 H =T +V =− + kx2 (20) 2m dx2 2 Như vậy, phương trình Schr¨dinger không phụ thuộc thời gian trong trường o hợp này được viết như sau d22 1 − 2 + kx2 ψ(x) = Eψ(x) (21) 2m dx 2 2m Đơn giản phương trình trên bằng cách nhân với − 2 , ta được d2 ψ(x) mkx2 2mE 2 − 2 ψ(x) = − 2 ψ(x) (22) dx hay d2 ψ(x) mkx2 2mE = − 2 ψ(x) (23) dx2 2 Đây là một phương trình vi phân không tuyến tính. Nó tương tự như một phương trình vi phân đã được giải bởi Charles Hermite, nhà toán học Pháp. 2.2.2 Hàm sóng của dao động điều hòa ở trạng thái cơ bản Giả sử tại một số điểm nào đó nghiệm của (23) có dạng 2 ψ(x) = ce−bx (24) với b, c là những hằng số. Như vậy, ta có dψ(x) 2 = −2bcxe−bx (25) dx d2 ψ(x) 2 2 = −2bce−bx + 4b2 cx2 e−bx (26) dx2 Thế (24) và (26) vào (23), ta được 2 2 mkx2 2 2mE 2 −2bce−bx + 4b2 cx2 e−bx = 2 ce−bx − 2 ce−bx (27) 2 Chia hai vế (27) cho ce−bx ta được mk 2mE 4b2 x2 − 2b = 2 x2 − 2 (28) Từ đó ta suy ra mk 2mE 4b2 = 2 và − 2b = − 2 √ 2 k ⇒b= mk/2 và E=b = m 2 m 5
Mặt khác, ta có tần số dao động điều hòa được tính bởi 1 k k ν= ⇒ = 2πν 2π m m Do đó k h 1 E= = 2πν = πν = hν (29) 2 m 2 2π 2 √ 1 Tóm lại, với b = mk/2 và E = hν, hàm 2 2 ψ(x) = ce−bx thỏa mãn phương trình Schr¨dinger cho hệ dao động điều hòa. Sử dụng giá o trị b đã tìm được, ta viết lại nghiệm ψ(x) như sau √ mk)x2 /2 ψ(x) = ce−( (30) Trong đó c là hằng số chuẩn hóa. Đây chính là hàm sóng của dao động điều hòa ở trạng thái cơ bản. Năng lượng của dao động điều hòa ở trạng thái cơ 1 bản là E = hν. 2 2.2.3 Các mức năng lượng của dao động điều hòa Sau đây, chúng ta sẽ giải phương trình Schr¨dinger của dao động điều hòa o một cách có hệ thống hơn so với cách chúng ta đã tiến hành ở trên. Ta viết lại phương trình sóng (23) như sau d2 ψ 2mE mk 2 2 + 2 − 2 x ψ=0 (31) dx Đặt 2mE mk α= 2 ; β2 = 2 (32) Phương trình (31) trở thành d2 ψ + (α − β 2 x2 )ψ = 0 (33) dx2 Thực hiện đổi biến phương trình (33) bằng cách đặt d2 d2 z= βx ⇒ =β 2 (34) dx2 dz Do đó, (33) tương đương với d2 ψ β + (α − βz 2 )ψ = 0 (35) dz 2 6
hay d2 ψ α + ( − z 2 )ψ = 0 (36) dz 2 β Nghiệm ψ(z) của (36) có thể được viết dưới dạng tích của hai hàm u(z) 2 và e−z /2 như sau1 2 ψ(z) = u(z)e−z /2 (37) Từ đó, ta có các đạo hàm bậc nhất ψ (z) và đạo hàm bậc hai ψ (z) 2 /2 2 /2 ψ = u e−z − uze−z (38) 2 /2 2 /2 2 /2 2 /2 2 /2 ψ = u e−z − u ze−z − u ze−z − ue−z + uz 2 e−z (39) Đơn giản (39), ta được 2 /2 2 /2 2 /2 2 /2 ψ = u e−z − 2u ze−z − ue−z + uz 2 e−z (40) Thế (37) và (40) vào (36), sau khi rút gọn, ta được d2 u du α 2 − 2z + −1 u=0 (41) dz dz β α Nếu ta đặt − 1 = 2n, thì (41) trở thành β d2 u du 2 − 2z + 2nu = 0 (42) dz dz Phương trình vi phân có dạng như (42) được gọi là phương trình Her- mite. Nghiệm của phương trình Hermite được gọi là các đa thức Hermite (Hermite polynomials). Trước khi giải (42), chúng ta xét các mức năng lượng của dao động điều hòa. Ta có α α − 1 = 2n ⇒ = 2n + 1 (43) β β Với 2mE mk α= 2 ; β= 2 α 2E m ⇒ = 2n + 1 = (44) β k Từ đó, ta tính được 1 k E = (n + ) (45) 2 m 2 1 Nghiệm e−z /2 được gọi là nghiệm tiệm cận, khi z → ∞ 7
h 1 k Vì = và ν = nên (45) trở thành 2π 2π m 1 E = (n + )hν (46) 2 Phương trình (46) là biểu thức tính năng lượng của dao động điều hòa. Chúng ta lưu ý vì năng lượng được lượng tử hóa nên số lượng tử n của dao động điều hòa nhận những giá trị nguyên không âm. Khi n = 0, ta có 1 E = hν và được gọi là năng lượng điểm không (zero-point energy). 2 Năng lượng điểm không có thể được xem là năng lượng của những dao động điều hòa ở không độ tuyệt đối. Nếu năng lượng điểm không bằng zero thì cả động năng và thế năng của nó cũng đều bằng zero. Động năng bằng zero nghĩa là Δpx = 0. Thế năng bằng zero nghĩa là hạt sẽ luôn đứng yên tại một điểm, hay Δx = 0. Tuy nhiên, theo nguyên lí bất định thì Δpx và Δx không thể đồng thời bằng zero. Như vậy, sự tồn tại của năng lượng điểm không là phù hợp với nguyên lí bất định Heisenberg. 2.2.4 Hàm sóng của dao động điều hòa Bây giờ chúng ta giải phương trình d2 u du 2 − 2z + 2nu = 0 (47) dz dz bằng phương pháp chuỗi lũy thừa. Nghiệm u(z) được viết dưới dạng chuỗi lũy thừa như sau ∞ 2 u(z) = a0 + a1 z + a2 z + · · · = ak z k (48) k=0 Lần lượt lấy đạo hàm bậc nhất và đạo hàm bậc hai (48), ta được ∞ u (z) = a1 + 2a2 z + 3a3 z 2 + · · · = kak z k−1 (49) k=1 ∞ 2 u (z) = 2a2 + 6a3 z + 12a4 z + · · · = k(k − 1)ak z k−2 (50) k=2 Thế (49) và (50) vào phương trình (47), ta được ∞ ∞ ∞ k(k − 1)ak z k−2 − 2z kak z k−1 + 2n ak z k = 0 (51) k=2 k=1 k=0 ∞ ∞ ∞ k−2 k k(k − 1)ak z − 2kak z + 2nak z k = 0 (52) k=2 k=1 k=0 8
Đặt υ = k − 2, suy ra k = υ + 2. Khi đó, ta có ∞ ∞ k(k − 1)ak z k−2 = (υ + 2)(υ + 1)aυ+2 z υ (53) k=2 υ=0 Theo lí thuyết chuỗi, ta có ∞ ∞ (υ + 2)(υ + 1)aυ+2 z υ = (k + 2)(k + 1)ak+2 z k (54) υ=0 k=0 vì υ và k là những biến số giả (dummy variables). Ví dụ, ta so sánh hai chuỗi 2 2 ci xi = c1 x1 + c2 x2 = cj xj = c1 x1 + c2 x2 i=0 j=0 Mặt khác, ta có ∞ ∞ 2kak z k = 2kak z k (55) k=1 k=0 vì khi k = 0 thì số hạng đầu tiên trong chuỗi là 2 × 0 × a0 × t0 = 0. Do đó, phương trình (52) được viết lại như sau ∞ ∞ ∞ (k + 2)(k + 1)ak+2 z k − 2kak z k + 2nak z k = 0 (56) k=0 k=0 k=0 Từ đó, ta có ∞ (k + 2)(k + 1)ak+2 − 2kak + 2nak z k = 0 (57) k=0 Để (57) đúng với mọi giá trị z thì biểu thức trong dấu móc vuông phải bằng zero (k + 2)(k + 1)ak+2 − 2kak + 2nak = 0 (58) 2n − 2k ak+2 = − ak (k = 0, 1, 2, . . .) (59) (k + 2)(k + 1) Phương trình (59) được gọi là công thức hồi qui (recursion formula). Dựa vào công thức hồi qui, nếu biết a0 ta sẽ tính được các giá trị a2k ; nếu biết a1 ta sẽ tính được các giá trị a2k+1 . Ví dụ 2n − 0 k = 0 : a2 = − a0 = −na0 (0 + 2)(0 + 1) 2n − 2 n−1 k = 1 : a3 = − a1 = − a1 (1 + 2)(1 + 1) 3 2n − 4 n−2 n(n − 2) k = 2 : a4 = − a2 = − a2 = a0 (2 + 2)(2 + 1) 6 6 2n − 6 n−3 (n − 3)(n − 1) k = 3 : a5 = − a3 = − a3 = a1 (3 + 2)(3 + 1) 10 30 9
Sau đây, ta trình bày lại các kết quả trên a2 = −na0 n−1 n−1 2(n − 1) a3 = − a1 = −2 a1 = − a1 3 2·3 3! n(n − 2) n(n − 2) 22 n(n − 2) a4 = a0 = 4 a0 = a0 6 2·3·4 4! (n − 3)(n − 1) (n − 3)(n − 1) 22 (n − 3)(n − 1) a5 = a1 = 4 a1 = a1 30 2·3·4·5 5! Do đó, nghiệm u(z) được viết như sau u(z) = a0 + a1 z + a2 z 2 + a3 z 3 + +a4 z 4 · · · 2(n − 1) 22 n(n − 2) = a0 + a1 z − na0 z 2 − a1 z 3 + a0 z 4 + · · · 3! 4! 22 n(n − 2) 4 2(n − 1) 3 = a0 (1 − nz 2 + z + · · · ) + a1 (z − z + ···) 4! 3! Đặt 22 n(n − 2) 4 u1 (z) = 1 − nz 2 + z + ··· 4! 2(n − 1) 3 u2 (z) = z − z + ··· 3! ⇒ u(z) = a0 u1 (z) + a1 u2 (z) (60) Như vậy, nếu n là số nguyên không âm thì nghiệm của phương trình Hermite là những đa thức, được gọi là đa thức Hermite và có thể được viết như sau n 2 d 2 Hn (z) = (−1)n ez n e−z (61) dz Sau đây là một số đa thức Hermite đầu tiên H0 (z) = 1 H1 (z) = 2z H2 (z) = 4z 2 − 2 H3 (z) = 8z 3 − 12z H4 (z) = 16z 4 − 48z 2 + 12 Hàm sóng của dao động điều hòa có dạng 2 /2 ψn = Nn Hn (z)e−z (n = 0, 1, 2, . . .) (62) với Nn là hằng số chuẩn hóa hàm sóng. Ta thấy, hàm sóng bằng zero khi Hn (z) = 0. Ví dụ, với n = 1 thì H1 (z) = 0 tại z = 0 và hàm sóng có một 10
1 node. Với n = 2 thì H2 (z) = 0 tại z = ± √ và hàm sóng có hai nodes. Một 2 cách tổng quát, số nodes của hàm sóng dao động điều hòa là n. Tóm lại, phương trình Schr¨dinger cho dao động điều hòa cũng được giải o một cách chính xác, dù hơi vất vả. Kết quả, năng lượng và hàm sóng được xác định như sau 1 E = (n + )hν (n = 0, 1, 2, . . .) 2 −z 2 /2 ψn (z) = Nn Hn (z)e (n = 0, 1, 2, . . .) Trong đó, Nn là hằng số chuẩn hóa; Hn (z) là đa thức Hermite bậc n; và mk z= βx = 2 x Như vậy, hàm sóng của dao động điều hòa viết theo biến x có dạng 2 /2 ψn (x) = Nn Hn (x)e−βx (n = 0, 1, 2, . . .) (63) 2 /2 Ta thấy e−βx là một hàm chẵn vì 2 /2 2 /2 e−β(−x) = e−βx Do đó, tính chẵn lẻ của hàm sóng phụ thuộc vào Hn (x). Nếu Hn (x) là hàm chẵn thì ψn (x) cũng là hàm chẵn và ngược lại. Đồ thị của hàm số chẵn thì đối xứng qua trục y. Vì vậy, nếu f (x) là hàm số chẵn thì +a a f (x)dx = 2 f (x)dx (64) −a 0 Hàm g(x) nếu thỏa điều kiện g(−x) = −g(x) thì được gọi là hàm số lẻ. Những hàm lẻ luôn nhận gốc tọa độ làm tâm đối xứng. Vì vậy , nếu g(x) là hàm số lẻ thì +a g(x)dx = 0 (65) −a 3 Phổ dao động của phân tử hai nguyên tử Một phân tử thông thường có các kiểu chuyển động là chuyển động tịnh tiến (translational motion), chuyển động quay (rotational motion) và chuyển động dao động (vibrational motion). Sự tịnh tiến của phân tử dưới dạng khí lý tưởng được mô tả giống như sự chuyển động của hạt trong hộp ba chiều. 11
Trong chuyển động quay, phân tử quay quanh trọng tâm của nó và luôn giữ khoảng cách giữa các nguyên tử không thay đổi. Kiểu chuyển động này chúng ta sẽ xét sau. Trong chuyển động dao động, các nguyên tử dao động gần như điều hòa dọc theo trục liên kết hoặc góc liên kết. Vì vậy, năng lượng dao động của phân tử ở những mức năng lượng thấp có thể được xem như năng lượng của dao động điều hòa. 3.1 Sự dao động của phân tử 1 H 35 Cl Phân tử 1 H 35 Cl gồm một nguyên tử Cl (chứa 17 proton và 18 notron) và một nguyên tử H (chỉ chứa 1 proton). Khoảng cách cân bằng giữa hai nguyên tử là 1, 27 · 10−10 m. Khối lượng của proton và notron theo đơn vị khối lượng nguyên tử (atomic mass units - amu) là mp = 1 amu; mn = 1 amu với 1 amu = 1, 66 · 10−27 kg Như vậy mH = 1 amu; mCl = 35 amu Sự dao động giữa hai nguyên tử H và Cl trong phân tử HCl rất giống với sự dao động của một con lắc lò xo cực nhỏ với hai đầu là những nguyên tử H và Cl. Giá trị thực nghiệm của hằng số đàn hồi (hằng số lực) k trong phân tử 1 H 35 Cl là 513 kg s−2 . Theo cơ học cổ điển, sự chuyển động tương đối của hai hạt có khối lượng m1 và m2 dọc theo một trục nối thì tương đương với sự chuyển động của một hạt có khối lượng rút gọn (reduced mass) là m1 m2 μ= (66) m1 + m2 Trong trường hợp phân tử HCl, ta có khối lượng rút gọn là 1 · 35 μ= 0, 97 amu 1 + 35 Tần số của dao động được tính như sau 1 k ν= = 8, 9 · 1013 s−1 2π μ Năng lượng điểm không của phân tử HCl 1 E0 = hν = 0, 18 eV 2 Và năng lượng ở trạng thái kích thích thứ nhất là 1 E1 = (1 + )hν = E0 + hν = 0, 55 eV 2 12
Khi phân tử HCl ở trạng thái kích thích thứ nhất về trạng thái điểm không, nó sẽ phát ra một photon γ có năng lượng là Eγ = E1 − E0 = 0, 36 eV Photon đó sẽ có tần số Eγ ν= = 8, 9 · 1013 s−1 h tương ứng với bức xạ trong vùng hồng ngoại (Infra-Red − IR). Phổ dao động của phân tử còn được gọi là phổ IR hay phổ hồng ngoại. Ngược lại, để kích thích phân tử HCl từ trạng thái điểm không lên trạng thái kích thích thứ nhất, ta cần cung cấp một năng lượng 0, 36 eV . 12 3.2 Xác định hằng số lực của C 16 O Phổ hồng ngoại của phân tử 12 C 16 O có mũi cực đại tại ν = 2143 cm−1 . ¯ Chúng ta thử xác định hằng số lực k của 12 C 16 O, xem sự dao động của hai nguyên tử C và O dọc theo trục liên kết là một dao động điều hòa. Cũng giống như nhiều phân tử hai nguyên tử khác, ở điều kiện thường thì phân tử CO chủ yếu nằm ở trạng thái điểm không. Khi bị kích thích nó sẽ lên trạng thái thứ nhất, ứng với sự dịch chuyển từ mức năng lượng E0 → E1 . Như vậy, mũi cực đại này tương ứng với sự dịch chuyển từ n = 0 → n = 1. Năng lượng cần cung cấp cho sự dịch chuyển này được tính như sau hc E = E1 − E0 = hν = = hc¯ ν λ Do đó ν = c¯ = (2143 cm−1 )(2, 998 · 10−10 cm/s) = 6, 424 · 1013 s−1 ν Khối lượng rút gọn của CO là m1 m2 12 · 16 μ= = = 6, 857 amu = 11, 383 · 10−27 kg m1 + m2 12 + 16 Để tính hằng số lực k, ta dựa vào biểu thức 1 k ν= 2π m với m là khối lượng rút gọn của CO. Như vậy k = 4π 2 ν 2 μ = 4π 2 (6, 424 · 1013 s−1 )2 (11, 383 · 10−27 kg) = 1855 N/m 13
4 Phương trình Hermite và đa thức Hermite Kí hiệu đa thức Hermite bậc n là Hn (x). Xét hàm g gồm hai biến x, t 2 +2tx g(x, t) = e−t (67) Giả sử g(x, t) được khai triển theo lũy thừa tn với các hệ số phụ thuộc vào Hn (x) như sau ∞ 2 Hn (x) n g(x, t) = e−t +2tx = t (68) n=0 n! Lấy đạo hàm bậc nhất theo t, ta được ∞ ∂ 2 tn−1 g(x, t) = (−2t + 2x)e−t +2tx = Hn (x)n (69) ∂t n! n=1 Ta có 2 +2tx 2 +2tx 2 +2tx (−2t + 2x)e−t = −2te−t + 2xe−t ∞ ∞ tn tn = −2t Hn (x) + 2x Hn (x) n! n! n=0 n=0 ∞ ∞ tn+1 tn = −2 Hn (x) + 2x Hn (x) n! n! n=0 n=0 Do đó (69) trở thành ∞ ∞ ∞ tn+1 tn tn−1 −2 Hn (x) + 2x Hn (x) = Hn (x)n (70) n! n! n! n=0 n=0 n=1 Vì n là biến số giả nên ta có ∞ ∞ ∞ tn−1 tn−1 tn Hn (x)n = Hn (x) = Hn+1 (x) n=1 n! n=1 (n − 1)! n=0 n! và ∞ ∞ ∞ tn+1 tn+1 tn Hn (x) = Hn (x)n = Hn−1 (x) n=0 n! n=0 (n + 1)! n=1 n! Như vậy, (70) tương đương với ∞ ∞ ∞ tn tn tn −2 Hn−1 (x) + 2x Hn (x) = Hn+1 (x) (71) n! n! n! n=1 n=0 n=0 Cho các hệ số trong (71) với cùng số mũ tn bằng nhau, ta được công thức hồi qui Hn+1 (x) = 2xHn (x) − 2nHn−1 (x) (n > 0) (72) 14
Tương tự, lấy đạo hàm (68) theo x, ta được ∞ ∂ 2 tn g(x, t) = 2te−t +2tx = Hn (x) (73) ∂x n! n=0 hay ∞ ∞ tn+1 tn 2 Hn (x) = Hn (x) (74) n! n! n=0 n=0 ∞ ∞ tn tn ⇒2 Hn−1 (x) = Hn (x) (75) (n − 1)! n! n=1 n=1 Từ đó, ta có Hn (x) = 2nHn−1 (x) (n > 0) (76) Kết hợp (72) và (76) ta được Hn+1 (x) = 2xHn (x) − Hn (x) (77) Lấy đạo hàm (77) theo x Hn+1 (x) = 2Hn (x) + 2xHn (x) − Hn (x) 2(n + 1)Hn (x) = 2Hn (x) + 2xHn (x) − Hn (x) Suy ra Hn (x) − 2xHn (x) + 2nHn (x) = 0 (78) Phương trình trên chính là phương trình Hermite. Nghiệm của nó là đa thức Hermite Hn (x). Sau đây, chúng ta xác định một số đa thức Hermite đầu tiên. Ta có ∞ Hn (x) n t2 t = H0 (x) + H1 (x)t + H2 (x) + · · · (79) n! 2! n=0 2 +2tx Mặt khác, áp dụng công thức khai triển Taylor cho hàm e−t theo lũy thừa (−t2 + 2tx)n ta được 2 +2tx (−t2 + 2tx)2 e−t = 1 + (−t2 + 2tx) + + ··· (80) 2! Do đó, ta có t2 (−t2 + 2tx)2 H0 (x) + H1 (x)t + H2 (x) + · · · = 1 + (−t2 + 2tx) + + · · · (81) 2! 2! Để hai vế phương trình trên bằng nhau thì các hệ số của cùng lũy thừa t phải bằng nhau. Vì vậy H0 (x) = 1; H1 (x) = 2x (82) 15
Dựa vào (82) và công thức hồi qui Hn+1 (x) = 2xHn (x) − 2nHn−1(x) ta sẽ xác định các đa thức Hermite có bậc cao hơn. Ví dụ H2 (x) = 2xH1 (x) − 1 × 2H0 (x) = 4x2 − 2 H3 (x) = 2xH2 (x) − 2 × 2H1 (x) = 2x(4x2 − 2) − 8x = 8x3 − 12x H4 (x) = 2xH3 (x) − 3 × 2H2 (x) = 16x4 − 48x2 + 12 Một cách tổng quát, theo Rodrigues các đa thức Hermite được viết như sau 2 dn −x2 Hn (x) = (−1)n ex e dxn Áp dụng công thức trên và kết hợp với một số hiểu biết về đạo hàm ta có thể dễ dàng xác định một đa thức Hermite bậc n bất kì. 16
Bài tập 1. Tìm mối liên hệ hồi qui cho các hệ số cn , từ đó tìm biểu thức tính c4 theo c0 và tính c5 theo c1 trong nghiệm chuỗi lũy thừa của phương trình sau f (x) − 2xf (x) + f (x) = 0 2. Chứng tỏ rằng 2 /2 f (x) = e−x là nghiệm tiệm cận (nghiệm ứng với x → ∞) của phương trình vi phân u (x) + (2 − x2 )u(x) = 0 3. Hàm sóng của dao động điều hòa ở trạng thái cơ bản và ở trạng thái kích thích thứ nhất như sau 2 /2 ψ0 (x) = c0 e−βx 2 /2 ψ1 (x) = c1 xe−βx a. Xác định các hằng số c0 và c1 b. Chứng minh rằng ψ0 (x) và ψ1 (x) trực giao với nhau. 4. Tính các giá trị trung bình T và V của dao động điều hòa ở trạng thái điểm không +∞ ∗ T = ψ0 T ψ0 dx −∞ +∞ ∗ V = ψ0 V ψ0 dx −∞ Chứng tỏ rằng T = V 5. Xác định các nodes của hàm sóng dao động điều hòa ứng với trạng thái n = 3. 6. Phổ IR của phân tử LiH có mũi cực đại tại ν = 1359 cm−1 . Xác định ¯ hằng số lực của liên kết trong phân tử này. 17