Luận văn ứng dụng mathematica trong một số bài toán vật lý phổ thông và vật lý lý thuyết

B 0 GIÂO DyC VÀ DÀO TAO
TRUÔNG BAIHQC SU' PHAM HÀ N 012

M È T IE N M A N H

ÎTNG DUNG MATHEMATICA
TRONG MOT
SÔ BÀI TOÂN VÂT
LŸ
•
•
PHÔ THÔNG YÀ VÂT LŸ LŸ THUYÉT

LUÂN VAN THAC Si KHOA HOC VÂT CHAT
•

•

•

HÀ NQI, 2015

•


B ộ GIÁO DỤC VÀ ĐÀO TẠO
TRƯỜNG ĐẠI HỌC s ư PHẠM HÀ NỘI 2

M È T IẾ N M Ạ N H

ỨNG DỤNG MATHEMATICA

TRONG MỘT SỐ BÀI TOÁN VẬT LÝ
PHỎ THÔNG VÀ VẬT LÝ LÝ THUYẾT
Chuyên ngành : V ật lí lí thuyết và V ật lí toán
M ã số

: 60 44 01 03

LUẬN VĂN THẠC SĨ KHOA HỌC VẬT CHẤT
•

*

*

HÀ NỘI, 2015

•


LỜI CẢM ƠN

Tôi xin chân thành cảm ơn các thầy cô ừong tổ Vật lý lý thuyết, Phòng
Sau đại học Trường Đại học Sư phạm Hà Nội 2, cùng các thầy, cô giáo đã tận
tình giảng dạy quan tâm tạo điều kiện giúp tôi hoàn thành khóa học.
Tôi xin bày tỏ lòng biết ơn sâu sắc tới TS.Trần Thái Hoa đã tận tình
chỉ bảo và giúp đỡ tôi trong suốt quá trình nghiên cứu và hoàn thiện luận văn.
Xin cảm ơn gia đình, bạn bè cùng các học viên lớp KI 7 VLLT & VLT
đã ủng hộ động viên và tạo mọi điều kiện trong thòi gian học tập, nghiên cứu
để hoàn thành luận văn.
Tôi xin chân thành cảm ơn mọi sự giúp đỡ vô cùng quý báu ấy!

Hà Nội, tháng 7 năm 2015
Tác giả

Mè Tiến Mạnh
9


LỜI CAM ĐOAN

Tôi xin cam đoan đề tài “ứng dụng Mathematica trong một số bài
toán vật lý phổ thông và vật lý lý thuyếrYd, đề tài do bản thân tôi nghiên cứu
dưới sự hướng dẫn của thầy giáo, TS. Trần Thái Hoa, Khoa Vật lý trường
Đại học sư phạm Hà Nội 2. Đề tài không hề trùng lặp vói bất kỳ một luận văn
nào, kết quả nghiên cứu không trùng vói tác giả khác.

Hà Nội, tháng 7 năm 2015
Người cam đoan

Mè Tiến Mạnh
«


MỤC LỤC
I. MỞ ĐẦU....................................................................................................... 1
1. Lý do chọn đề tà i....................................................................................... 1
2. Mục đích nghiên cứu.................................................................................2
3. Nhiệm vụ nghiên cứu.................................................................................2
4. Đối tượng và phạm vi nghiên cứu.............................................................2
5. Những đóng góp mới của đề tài.................................................................2

6. Phương pháp nghiên cứu...........................................................................2
II. NỘI DUNG................................................................................................... 3
Chương 1. MỘT VÀI NÉT VỀ PHẦN MỀM MATHEMATICA................... 3
1.1. Giới thiệu sơ bộ về phần mềm Mathematica.......................................... 3
1.2. Khả năng đồ họa của phần mềm Mathematica....................................... 3
1.3. Một số hàm thông dụng của Mathematica............................................. 5
Chương 2. ỨNG DỤNG PHẦN MỀM MATHEMATICA VÀO GIẢI
QUYẾT CÁC BÀI TOÁN VẬT L Ý ................................................................6
2.1. Một số bài toán vật lý phổ thông............................................................6
2.1.1. Bài toán về chuyển động ném xiên .................................................6
2.1.2. Xử lỷ sổ liệu khi làm thực hành ở phổ thông ..................................8
2.2. Một số bài toán về cơ học lượng tử và vật lý thống k ê .........................11
2.2.1. Cơ học lượng tử........................................................................... 11
2.2.2. Vật ỉỷ thống kê ............................................................................ 20
III. KẾT LUẬN...............................................................................................26
IV. DANH MỤC CÁC TÀI LIỆU THAM KHẢO........................................ 27
PHỤ LỤC


1

I. MỞ ĐẦU
1. Lý do chọn đề tàỉ
Để giải quyết nhiều vấn đề ừong vật lý, ừong kỹ thuật, ừong toán học,...
ngưòi ta phải sử dụng các phần mềm toán học. Ngay từ những năm 1960 đã
xuất hiện những bó phần mềm kỹ thuật đầu tiên dựa ừên các hệ đại số tượng
trưng (symbolic algebraic system). Thế hệ thứ nhất của nó là ngôn ngữ
Macsyma và Reduce, chủ yếu dùng cho các bài toán vật lý năng lượng cao,
nhưng nó có nhược điểm là được định hướng chạy chủ yếu ừên các máy tính
lớn (main-frame computer). Thế hệ thứ hai là ngôn ngữ Maple so vói thế hệ

trước có ưu điểm là chạy nhanh hơn và chấp nhận bộ nhớ nhỏ hơn (do vậy có
thể chạy ừên máy tính cá nhân) và được bổ sung nhiều khả năng đại số và đồ
thị hơn. Thế hệ thứ ba của dạng ngôn ngữ này chính là các ngôn ngữ
Mathematica và MatLab (bản có bổ sung phần tính toán đại số tượng trưng).
Trong đó Mathematica có ưu điểm vượt ừội về giao diện thân thiện, về khả
năng vẽ đồ thị siêu việt và khả năng xử lý dữ liệu không thua kém các ngôn
ngữ tính toán khác
Mặc dù lúc đầu ứng dụng của Mathematica chủ yếu trong các lĩnh vực
vật lý, kỹ thuật và toán, tuy nhiên việc ứng dụng của Mathematica ngày càng
được mở rộng ra các lĩnh vực khác như sinh học, các khoa học xã hội nhờ khả
năng mô hình hóa và mô phỏng các hệ lớn, kể cả các hệ động. Hiện nay nó
được sử dụng trong tất cả các công ty có tên trong Fortune 50, trong tất cả 15
bộ của chính phủ Mỹ và được giảng dạy ừong tất cả 50 trường tổng hợp hàng
đầu thế giới. Nó trở thành chương trình ứng dụng lớn nhất được phát triển và
chứa một số lượng lớn các thuật toán và các đổi mới kỹ thuật quan ừọng. Một
trong những sáng kiến kỹ thuật là môi trường phần mềm dựa ừên giao diện
tương tác được biết đến với tên là notebook. Hiện nay đã có vài trăm chương
trình đặc chủng viết ừên Mathematica được thương mại hóa một số đầu tạp
chí chuyên nghành và khoảng 200 đầu sách về ngôn ngữ M athematical]


2

Với những tính năng ưu việt của phàn mềm toán học Mathematica như
khả năng tính toán, khả năng đồ họa, cũng như tính dễ sử dụng của nó ừong
việc xây dựng các mô hình và giải quyết các bài toán vật lý và cũng muốn
mọi người có thêm một công cụ hữu ích để làm việc. Nên bản thân tôi đã
chọn đề tài: “ửng dụng Mathematica trong một sổ bài toán vật lý phổ thông
và vật lý lý thuyết”.
2. Mục đích nghiên cứu

Sử dụng phần mềm Mathematica áp dụng để giải các bài toán vật lý phổ
thông, cơ học lượng tử và vật lý thống kê.
3. Nhiệm vụ nghiên cứu
Lựa chọn các bài toán vật lý và lập trình bằng Mathematica để giải các
bài toán vật lý này.
4. Đổi tượng và phạm vỉ nghiên cứu
Vật lý tính toán và áp dụng phần mềm Mathematica để giải một số bài
toán vật lý phổ thông, vật lý thống kê và cơ học lượng tử.
5. Những đóng góp mói của đề tài
Tìm hiểu rõ hơn về phần mềm Mathematica, cách sử dụng và các tính
năng của phần mềm. Lập trình các bài toán vật lý bằng phần mềm
Mathematica để giúp các bạn học viên, sinh viên có thể dễ dàng sử dụng cộng
cụ Mathematica để thuận lợi hơn ừong quá trình nghiên cứu khoa học cũng
như giải bài tập.
6. Phương pháp nghiên cứu
Sử dụng các phương pháp nghiên cứu của vật lý lý thuyết và vật lý toán
để nghiên cứu các bài toán vật lý phổ thông, vật lý thống kê và cơ học lượng
tử và nghiên cứu các tài liệu về phần mềm Mathematica.


3

II. NỘI DUNG
Chương 1. MỘT VÀI NÉT VÈ PHẦN MÈM MATHEMATICA
1.1. Gỉói thỉệu so bộ về phần mềm Mathematỉca
Mathematica là ngôn ngữ tích hợp đầy đủ các tính toán kỹ thuật
(technical computing), là dạng ngôn ngữ dựa trên nguyên lý xử lý các dữ liệu
tượng trưng (symbolic manipulation). Khỏi thủy của nguyên lý này là ngôn
ngữ LIPS - ngôn ngữ nghiên cứu ừí tuệ nhân (artificial intellect) - nghiên cứu
các vấn đề như xử lý tiếng nói tự nhiên, các hệ chuyên gia (expert system),

các vấn đề logic trong kỹ thuật robot (robotech), điều khiển và tự động hóa.
Thế hệ ngôn ngữ giải tích đầu tiên đó là ngôn ngữ Macsyma và
Reduce, chủ yếu dùng cho các bài toán vật lý năng lượng cao, nhưng chúng
lại có nhược điểm là được định hướng chạy trên các máy tính lớn.
Thế hệ thứ hai là ngôn ngữ Maple so với thế hệ trước là có ưu điểm là
chạy nhanh hơn và chấp nhận bộ nhớ nhỏ hơn và bổ sung nhiều khả năng đại
số và đồ thị hơn và nó thể chạy trên máy tính cá nhân.
Thế hệ thứ ba của dạng ngôn ngữ này chính là các ngôn ngữ
Mathematica và MatLab, trong đó Mathematica có ưu điểm vượt ừội về giao
diện thân thiện, về khả năng vẽ đồ thị siêu việt và khả năng xử lý dữ liệu
không thua kém các môi trường ngôn ngữ tính toán khác.
Nhờ khả năng siêu việt của mình Mathematica không chỉ được ứng dụng
trong các lĩnh vực vật lý, kỹ thuật và toán mà còn mở rộng ừong các lĩnh vực
phức tạp khác như sinh học, khoa học xã hội, ...
Phiên bản đầu tiên của Mathematica được phát hành 23/6/1988. Bản 2.0
được phát hành năm 1991. Bản mới nhất của Mathematica là bản 10.0.2.
1.2. Khả năng đồ họa của phần mềm Mathematỉca
Mathematica cho phép vẽ tất cả các dạng đồ thị của một hàm số vói cấu
trúc lệnh đơn giản nhất như đồ thị hai chiều, đồ thị ba chiều, đồ thị đường
viền, đồ thị mật độ, ... [8]


4

Ví dụ ta sử dụng lệnh sau để vẽ đồ thị của hàm số cos2x+sinx trong đoạn
[0, 20] {Hình 1.1)
Plot [cos [2x] + sin [x], {jt,0,20}]

Hình 1.1
Hoặc ta có thể dùng lệnh sau đây để vẽ đồ thị ba chiều của hàm số

Cos(3xy2) (Hình 1.2)
Plot3D[cos[3xy A2],{jt,0,2},{);,0,2}]

Hình 1.2


5

1.3. Một sổ hàm thông dụng của Mathematỉca
Trong Mathematica

Biểu thức toán

Sqrt[x]
Log[x]

Ln(x)

Sin[x]

Sin(x)

Cos[x]

Cos(x)

Tan[x]

Tan(x)

Log[a,b]

logab

Arcsin[x]

Arcsin(x)

Exp[x]

ex

Factoria[n], n!

n!

Mod[n,m]

Số dư của —
m

Factorlnteger

Phân tích ra thừa sô nguyên sô của n

Abs[x]

Giá trị tuyệt đ ô i của X

xy

xy
Vx

x*y

hoặc X y

Xy

Pi

71

Limit[f(x),x^ x0]

Tính giới hạn

Sum[Function, {U m in,

}

i max , ]

Tính tổng

D[f(x),x]

Tính đạo hàm

Integrate[f(x),x]

Tính nguyên hàm

Integrate [f(x),{x,a,b}]

Tính tích phân xác định

Solve[f(x)==0,x]

Giải phương trình

Solve[fi==0, f == 0,{x,y}]

Giải hệ phương trình

Simplify[f(x),x]

Đơn giản biểu thức

Plot[f(x),{x,a,b}]

Vẽ đô thị


6

Chương 2. ỨNG DỤNG PHẦN MỀM MATHEMATICA
VÀO GIẢI QUYẾT CÁC BÀI TOÁN VẬT LÝ
2.1. Một sổ bàỉ toán vật lý phổ thông

2.1.1. Bài toán về chuyển động ném xiên
A. Lý thuyết
Phương trình chuyển động:
X

= (v0cosa)t, y = (v 0 sin a )

gt

t- —

(2.1)

Phương trình quỹ đạo:
(2.2)

Vận tôc của vật:
v 0x

= v0cosa, v0y = v0sin a, V= y jv l+ v 2y

(2.3)

Thời gian rơi:
(2.4)
Tầm bay xa và độ cao cực đại:
V ổ s in 2 a

(2.5)

g
vổ sin2a

H —-------------

2g

(2.6)

Nếu có lực tác dụng vào vật khi đang bay (ví dụ sức cản của không khí) thì
áp dụng định luật II Newton:
I1—lim

vói 1 ■-

T ±'2T ... T ±n là hợp lực của các lực tác dụng lên vật.

Sau đó chiếu lên các trục tọa độ để tìm phương trình chuyển động. [5]

(2.7)


7

Chiếu lên Ox:

=

F+F

+
m

+F

^

(2.8)

<■
•• E +Ẹ, +... + F
Chiếu lên Oy: y =
^
m

(2.9)

X

B. Bài toán: Một vật có khối lượng m được ném đi với tốc độ v0 hợp với
phương ngang một góc a . Hiệu ứng cản của không khí được mô tả bởi lực
cản 1 c . Lực này tỷ lệ với bình phương tốc độ, được biểu diễn bằng biểu thức:
±"c — - l l l i V V V

Trong đó V và

V

là tốc độ và vận tốc của vật,к là hệ số cảncủa không khí. Vẽ

đồ thị mô tả quá trình chuyển động của vậtkhi hệ sốcản là

k=5,2.10"572 và

k=5,2.10"3 (đồ thị phải mô tả được chuyển động của vật khi không có lực
cản). Từ đó rút ra vai trò của hệ số cản của không khí trong chuyển động ném
xiên với v0=50 m/s và а =45°.
Kết quả thu được khi thực hiện trên Mathematica [Phụ lục 1]: (Hình 2.1)
Chú thích: Đồ thị khi không có lực cản là màu đỏ, khi có hệ số cản kl
là màu xanh lục, khi có hệ số cản k2 là xanh lam.
Đồ thị ném xiên trong các trường hợp

y(t)

Hình 2.1
Từ ừên biểu đồ ta dễ dàng nhận thấy khi hệ số cản к càng lớn thì chuyển


8

động của vật ném xiên sẽ nhanh dừng lại hơn, tầm ném xa và độ cao cực đại
giảm dẫn tới có quỹ đạo nhỏ dần.
2.1.2. Xử lý sổ liệu khi làm thực hành ở phổ thông
2.1.2.1. Thực hành đo gia tốc rơi tự do
A. Lý thuyết
Thả một vật (trong đó trọng lượng của vật rất lớn so với lực cản của
không khí) từ độ cao s, khi đó chuyển động của vật coi như chuyển động rơi
tự do.
Khi một vật chuyển động có vận tốc ban đầu bằng không, chuyển động
nhanh dần đều với gia tốc a, thì quãng đường đi được s sau thời gian t được

xác định bởi công thức:
s = —at2

(2.10)

Đồ thị biểu diễn quan hệ giữa s và t2 có dạng một đường thẳng đi qua
gốc tọa độ và hệ số góc:
tana = —
2

(2.11)

B. Xử lý số liệu
Sau khi làm thí nghiệm về "Khảo sát chuyểnđộng rơi tự do. Xác định
gia tốc roi tự do", người ta đo được một bảnggiá ừị của s vàt như sau (trong
đó s là độ cao vật được thả và t là thời gian rơi):
s(m)

ti(s)

Í2(s)

t3(s>

u(s)

0,2

0,198

0,199

0,199

0,201

0,3

0,247

0,247

0,246

0,247

0,4

0,283

0,284

0,285

0,285

0,5

0,318

0,318

0,319

0,318

0,6

0,349

0,348

0,35

0,349

0,7

0,379

0,379

0,38

0,378


9

Từ bảng số liệu trên:

a, Nghiệm lại công thức s = —gt2.(Chứng minh hệ số« —= 4.905 )
2
t 2
b, Hãy xử lý số liệu để xác định gia tốc rơi tự do g. Và vẽ đồ thị s=s(t2)
và v=v(t), từ đó nhận xét đặc điểm của sự rơi tự do.[1,2]
Kết quả thu được khi thực hiện trên Mathematica [Phụ lục 2]:
Hệ số giữa quãng đường s và bình phương thời gian t2 là: 4.91494
Vậy giá trị của gia tốc rơi tự do là g= 9.88386±0.0696416

ĐỒ thị s=s(f2)
s(m)
0.7 ;
0.6
0.5
0.4
0.3
0.2
0.1

0.02

0.04

0.06

0.08

Hình 2.2

0.10

0.12

0.14

iV >


10

Đồ thị v=v(t)
v(m/s)

1.5

1.0

0.5

— ..................................................................................................
0.05

0.10

0.15

0.20

0.25

0.30

0.35

t(s)

Hình 2.3
Nhận xét: Từ đồ thị s=s(t2), hình 2.2 có dạng 1 đường thẳng. Như vậy,
chuyển động rơi tự do là chuyển động thẳng. Mặt khác đồ thị v=v(t) Hình 2.3
có dạng là đường thẳng hướng lên suy ra vận tốc tăng đều theo thời gian. Vậy
chuyển động rơi tự do là chuyển động thẳng nhanh dần đều.
2.1.2.2. Thực hành đo hệ sổ ma sảt trượt
A. Lý thuyết
Cho một vật nằm trên mặt phẳng nghiêng p với góc nghiêng a so vói
mặt nằm ngang. Khi a nhỏ, vật vẫn nằm yên ừên p không chuyển động. Khi
tăng độ nghiêng vật chuyển động xuống với gia tốc a. Độ lớn của a chỉ phụ
thuộc vào góc ngiêng a và hệ số |at - gọi là hệ số ma sát trượt:
a = g(sin a - |_itcosa)

(2.12)

Bằng cách đo a v à a , ta xác định được hệ số ma sát trượt \xt :
|at = ta n a ----- -—
gcosa

(2.13)


11

Gia tốc a đươc xác đinh theo công thức: a =

t

, trong đó quãng đường

đi được s đo bằng thước, thời gian t đo bằng đồng hồ đo thòi gian hiện số,
điều khiển bằng công tắc và cổng quang điện. Góc nghiêng a có thể đọc
ngay ừên thước đo góc có gắn quả dọi, gắn vào mặt phẳng nghiêng.
B. Xử lý số liệu
Trong thí nghiệm về "Xác định hệ số ma sát", khi tăng độnghiêng của
máng lên 20° thì người ta thấy vật bắt đầu chuyển động xuống với gia tốc a.
Biết độ dài quãng đường đo là s=0.6m, người ta đo được các giá trị thòi gian
vật đi hết quãng đường s là t trong 5 lần đo là t={ 1.018, 1.022, 1.038, 1.036,
1.028}. Hãy xử lý số liệu trên để tính hệ số ma sát trượt giữa mặt phẳng
nghiêng và vật. Lấy g = 9,814m/s2.
Kết quả thu được khi thực hiện trên Mathematica [Phụ lục 3] :
Hệ số ma sát trượt giữa vật và máng lầL\xt = 0,240915 ± 0,00164222.
2.2.

Môt
sổ bài toán về cơ •hoc lương
tử và•vât•/ lý thống
kê
•
• о
Ü

2.2.1. Cơ học lượng tử
2.2.1.1. Cảc bài toán về tính trị trung bình

A. Một sổ kiến thúc để làm bài tập:
Phương trình hàm riêng và trị riêng của toán tửF:
îV „(q^) = f„V„(q,t)

(2.14)

Mô tả phép đô đại lại F để có số đofn. Trong phép đo đó, hệ lượng tò
chuyển từ trạng thái nào đó về nằm ở ừạng thái liên kết \|/ . Trường hợp
phép đo F không thực hiện được thì khi liên kết với hệ máy đo hệ lượng tử
không chuyển về trạng thái \|/ .[4,7]
Hệ hàm riêng {\|/ } của toán tử F là một hệ cơ sở trực chuẩn đủ của
không gian Hilbert các hàm trạng thái, vì vậy có thể khai triển duy nhất theo
hệ cơ sở này:


12

v = Z cnVn

(2.15)

n

Vì hệ {\|/ } là trực chuẩn đủ, dễ dàng rút ra:
cn= J v > d q = (\|/n,\|/)

(2.16)

Như vậy, xác suất để đo được giá trị f xác định của đại lượng F của
hệ lượng tử ừong trạng thái được mô tả bởi hàm sóng \|/(q) nào đó sẽ được

tính bởi
w (fn) = |cn|2 = |JV;(q>|/(q)dq|2

(2.17)

Nếu hệ lượng tử ở trạng thái mô tả bởi hàm sóng \|/ Ф|\|/ } (n=l,2,...),
các số đo f của đại lượng F có xác suất đo là |c |2. Theo lý thuyết xác suất,
các số đo f (n=l,2,..) sẽ có trị trung bình.
(2.18)
Khi \|/ = \|/ thì:
F = J\|/*Fi|/dq = f J\|/*\|/dq = fn; nghĩa là F = fn.

(2.19)

B. Một sổ bài toán về trị trung bình.
Bài 1: Hạt ở trạng thái được mô tả bởi hàm sóng:
\|/(х) = А.е(а“-“2х2). trong đó
trị trung bình:

X, X 2,

-0 0

< X < 00, a và к là các hằng số. Hãy tìm các

p , p2, (Ах)2, (Ap )2 và nghiệm lại hệ thức bất định.

Kết quả thu được khi thực hiện trên Mathematica [Phụ lục 4] :
Giá trị của X= 0
Giá ừi của X2 = —

4a
Giá ừi của Лх2 = —
4a


13

Giá trị của p =kỀ
Giá ừị của p2 = (k2+ a )h
Giá ừị của Лр* = ah
Thỏa mãn hệ thức bất định
Bài 2: Hạt chuyển động ừong giếng thế chữ nhật một chiều có thành cao vô
hạn được mô tả bởi hàm sóng đã chuẩn hóa: [6]

ừong đó d là bề rộng giếng thế và n=l,2,3—
a, Dùng hệ thức bất định ước tính mức năng lượng thấp nhất có thể có
của hạt.
b, Tính các giá trị trung bình của X, X2, Ax2, p , p2 và động năng
trung bình T của hạt.
Kết quả thu được khi thực hiên trên Mathematica [Phụ lục 5] :

Giá tri của E Min=

— ,

2 dm

Giá ừị của X= —

2

n 71 )
n n J
Giá ừị của p = 0


14

G iá ừ ic ủ a ? = í ^ *
' l d!

,

J

Giá ừị của động năng trung bình T = í

n2Ti2h
2d2m )

Bài 3: Hàm sóng chuẩn hóa của nguyên tà đồng dạng Hydro ở trạng thái cơ
bản là:
\|/ = A exp(-pr) ; trong đó A và p là những hằng số
a, Chứng minh rằng
A2 = £
71 ’

B= - m Z
h

e4 m Z 2
h

Trị trung bình của thế năng và động năng lần lượt là 2E và -E
b, Tính r và r2.
Kết quả thu được khi thực hiện trên Mathematica [Phụ lục 6]:
P3

z

= e_ 6m3

71

ry 3

71h

0 = Tíi z
e4^ 2
E = ---- 2h
ũ = 2 E = -^ 5
h
T = -E = ^

2h

3h
r ——- _
2e mZ

3
3/ỉ
r = 4e m 2V2
z

2


15

2.2.1.2. Các bài toán về nhiễu loạn
A. Lý thuyết
Ta sẽ đặt điều kiện hạn chế cho bài toán nhiễu loạn, trước hết ta xét lý
thuyết nhiễu loạn cho các bài toán có phổ gián đoạn.
H\|/j = E\|/1(1 = 1,2,3...)

(2.20)

Giả thiết toán tử H có thể tách ra làm hai phần:
H = H0 + V,

(2.21)

Trong đó, H0 là toán tử Hamilton đã được lý tưởng hóa, còn số thứ hai được
gọi là toán tử nhiễu loạn.
Biểu thị V là nhỏ, để biểu diễn điều đó ta đặt:
V = ẦW

(2.22)

Với X là một thông số nhỏ không thứ nguyên.
Giả sử biết các nghiệm E° và Ф, (1 = 1,2,3...) của phương trình cho hàm
riêng và ừị riêng của toán tử H0.
H qCPj = E “
(2.23)

Và giả thuyết các J ọ > 1dq = SrfG,l, = l,2 ĩ3...)

(2.24)

Với các điều kiện hạn chế đó, việc giải phương trình:

Hcp = Ecp

(2.25)

sẽ quy về việc giải phương trình sau để tìm Ej và \ị/j.
(H 0 + A,W)\|/j = Ejii/j

(2.26)

Nói khácđi chúng ta sẽ hiệu chỉnh cho E° và íp, (1 = 1,2,3...) để sau khi
hiệu chỉnh, các giá trị hiệu chỉnh Ej và \|/j sẽ nghiệm đứng (2.20), (2.25) hay
(2.26).


16

Khi xét bài toán nhiễu loạn dừng sẽ có 2 trường hợp xảy ra:
+Trường hợp bài toán lý tưởng không có suy biến.
+Trường hợp bài toán lý tưởng có suy biến.
B. Một sổ bài toán về lý thuyết nhiễu loạn.
Bài 1 : Xét một elecừon trong hộp một chiều có chiều dài 1 Ả .
a, Hãy tìm 4 hàm sóng đầu tiên (các hàm sóng chuẩn hóa) và vẽ dạng
đồ thị của chứng.
b, Tính 4 mức năng lượng tương ứng và vẽ sơ đồ mức năng lượng.
c, Tại t=0, hạt ở trạng thái n=l. Tại t=0, một giếng thế năng dạng chữ
nhật

v0=-104 eV, có tâm tại a/2 và độ rộng

10"12cm bất ngờ được đưa thêm

vào giếng thế ban đầu và duy trì ừong 5.10"18s, sau đó bị ngắt đi. Sau khi mất
đi sự nhiễu loạn thì xác suất để thấy hệ ở từng trạng thái n=2, n=3, n=4 là bao
nhiêu? (Chiều cao và chiều rộng của giếng thế đặc trưng cho 1 noừon tương
tác với 1 elecừon).
Kết quả khi thực hiện trên Mathematica [Phụ lục 7] :
(*Kết quả của 4 hàm sóng đầu tiên*)
Bốn hàm sóng đầu tiên là:

2

2+

{ { h E ^ ^ -^ - }}(*Vậy ta đã tính được E = ^ -^ - *)
2a m
‘

2a m


17

(*Kết quả 4 mức năng lượng đầu tiên*)
Bốn mức năng lượng đầu tiên tính theo Jun là

Ei= 6.02893 X l0 '18
E2= 2.41157 xlO'17
E3= 5.42604X10'17
E4= 9.64629 X 10'17
(*Vẽ đồ thị hàm sóng và năng lượng. Màu đồ thị n = l: hồng, n=2: đỏ, n=3:
xanh lam, n=4: xanh lục*)

lịí(x)

Đồ thị 4 hàm sóng đầu tiên

Hình 2.4


18

Đ ồ thị 4 m ứ c năng lượng đầu tiên
Enụ)
1.

*

1 0 ~ 16

s. *1er17

6

. X 1 e r 17

4

. *

1(

2

. *

1 С Г 17

r

17

2

.*

1 C T 11

4

.Х

1 С Г 11

6

.Х

1 С Г 11

8

.Х

1 С Г 11

1.

X1 0

Hình 2.5
(*Kết quả thu được*)
Xác suất tìm thấy hệ tại trạng thái 2 sau khi mất đi nhiễu loạn là p2= 0.
Xác suất tìm thấy hệ tại trạng thái 3 sau khi mất đi nhiễu loạn là
p3=0.000583387
Xác suất tìm thấy hệ tại trạng thái 4 sau khi mất đi nhiễu loạn là p4= 0.
Bài 2: Xét một dao động tử điều hòa một chiều có tần sốco0. Kí hiệu các trị
riêng năng lượng bằng n, bắt đầu từ n=0 ứng với giá trị năng lượng thấp nhất.

Một thế nhiễu loạn không phụ thuộc vào thời gian H = V (x) được thêm vào
thế năng dao động tử ban đầu. Thay vì đưa ra dạng của thế nhiễu loạn V(x), ta
sẽ chỉ ra tường minh các phần tử ma trận của nó được tính toán ừong biểu
diễn của các trạng thái riêng không nhiễu loạn. Một phần của ma trận đó được
chỉ ra dưới đây ừong đó 8 là một hằng số nhỏ và không có thứ nguyên.


19

/
1

0

0

0

- i
V2
0

0
0

0

0

0

0

0
0

0

0
V

0

y

a, Hãy tìm các năng lượng mói của năm mức năng lượng thấp nhất tính
đến bậc nhất ừong lý thuyết nhiễu loạn.
b, Tìm các năng lượng mới cho n=0 và n=l tới bậc hai ừong lý thuyết
nhiễu loạn.
Kết quả thu được khi thực hiện trên Mathematica [Phụ lục 8]:
Các mức năng lượng đầu tiên tính tới gần đứng bậc nhất theo lý thuyết nhiễu
loạn là
E' =E„ +H'

Như vậy 5 mức năng lượng đâu tiên là:

E

3ooO h

3ữ)0 h
E 2= ^ 2
E

7coO h

h
V2


20

Hai mức năng lượng đầu tiên tính tới gần đứng bậc 2

3ooO h
E =——

1

2

2.2.2. Vật lý thống kê
Trong khuôn khổ luận văn này chúng ta chỉ xét đến phần các đại lượng
nhiệt động.[3]
A. Tích phân trạng thái và các hàm nhiệt động.
Tích phân trạng thái Z:
(2.27)
Năng lượng tự do F:
F = -kTlnZ

(2.28)

Áp suất của hệ P:
(2.29)
Enừopi của hệ:
S=- ( —

U t Jv

= klnZ + kT

(2.30)

Nội năng u của hệ:
(2.31)
Entanpi của hệ:
H = Ư + PV = kT

fõ \n Z }
— — +kT

(2.32)