ĐẠI HỌC QUỐC GIA HÀ NỘI
TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN
Vũ Thị Thu Trang
TÁN XẠ HẠT NHÂN CỦA CÁC NƠTRON
PHÂN CỰC TRÊN MẶT TINH THỂ
CÓ CÁC HẠT NHÂN PHÂN CỰC
TRONG ĐIỀU KIỆN CÓ PHẢN XẠ
Chuyên ngành: Vật lý lý thuyết và vật lý toán
Mã số : 60.44.01
LUẬN VĂN THẠC SĨ KHOA HỌC
NGƯỜI HƯỚNG DẪN KHOA HỌC
PGS.TS NGUYỄN ĐÌNH DŨNG
Hà Nội - 2011
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
2
MỤC LỤC
Mở đầu: 3
Chương 1 - Lý thuyết tán xạ của nơtron chậm trong tinh thể . . . . . . . . . . . . . . 5
1.1. Cơ sở lý thuyết tán xạ của nơtron chậm trong tinh thể . . . . . . . . . . . . . . 5
1.2. Thế tương tác của nơtron chậm trong tinh thể . . . . . . . . . . . . . . . . . . . . 8
1.2.1. Yếu tố ma trận của tương tác hạt nhân . . . . . . . . . . . . . . . . . . . . . 8
1.2.2. Yếu tố ma trận của tương tác từ. . . . . . . . . . . . . . . . . . . . . . . . . . . 9
Chương 2 - Tiến động hạt nhân của spin của các nơtron trong
môi trường phân cực . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.1. Tính góc tiến động bằng phương pháp toán tử. . . . . . . . . . . . . . . . . . . 11
2.2. Tính góc tiến động bằng phương pháp hàm sóng. . . . . . . . . . . . . . . . . 13
2.3. Sử dụng bảo toàn năng lượng để tính góc tiến động. . . . . . . . . . . . . . . 16
Chương 3 - Phản xạ gương của nơtron phân cực trên mặt biên gồ ghề
giữa chân không và vật chất có các hạt nhân phân cực . . . . . . . . 18
3.1. Ảnh hưởng của sự gồ ghề mặt biên “chân không – vật chất” có các
hạt nhân phân cực lên phản xạ gương của các nơtron phân cực. . . . . . 18
3.2. Véctơ phân cực của nơtron phản xạ gương trên mặt biên gồ ghề
giữa chân không và vật chất có các hạt nhân phân cực. . . . . . . . . . . . . 23
Chương 4 - Tán xạ hạt nhân của các nơtron phân cực trên mặt tinh thể
có các hạt nhân phân cực trong điều kiện có phản xạ. . . . . . . . . . 26
4.1. Tiết diện hiệu dụng của tán xạ không đàn hồi của các nơtron
trên tinh thể có các hạt nhân phân cực. . . . . . . . . . . . . . . . . . . . . . . . . 26
4.2. Tiết diện tán xạ bề mặt hiệu dụng của các nơtron trong
trường hợp có phản xạ toàn phần. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
Kết luận: . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
Tài liệu tham khảo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
3
MỞ ĐẦU
Trong những năm gần đây, sự tán xạ của nơtron chậm đã được sử dụng rộng
rãi để nghiên cứu vật lý các chất đông đặc.
Các nơtron chậm là một công cụ độc đáo trong việc nghiên cứu động học của
các nguyên tử vật chất và các cấu trúc từ của chúng [14,15,19,20,21]
Hiện nay, để nghiên cứu cấu trúc tinh thể, đặc biệt là cấu trúc từ của tinh thể,
phương pháp quang nơtron đã được sử dụng rộng rãi. Chúng ta dùng chùm nơtron
chậm phân cực bắn vào bia (năng lượng cỡ dưới 1 MeV và không đủ để tạo ra quá
trình sinh hủy hạt ). Nhờ nơtron có tính trung hòa điện, đồng thời môment lưỡng
cực điện vô cùng nhỏ (gần bằng 0) nên nơtron không tham gia tương tác điện dẫn
đến độ xuyên sâu của chùm nơtron vào tinh thể là rất lớn, và bức tranh giao thoa
của sóng tán xạ sẽ cho ta thông tin về cấu trúc tinh thể và cấu trúc từ của bia.
Điều đó giúp ta hiểu rõ hơn về sự tiến động spin của các nơtron trong bia có các
hạt nhân phân cực [2,9,17,18,25]
Các nghiên cứu và tính toán về tán xạ phi đàn hồi của các nơtron phân cực
trong tinh thể phân cực cho phép chúng ta nhận được các thông tin quan trọng về
tiết diện tán xạ của các nơtron chậm trong tinh thể phân cực, hàm tương quan spin
của các hạt nhân [11,12,13,25]. Ngoài ra các vấn đề về nhiễu xạ bề mặt của các
nơtron trong tinh thể phân cực đặt trong trường ngoài biến thiên tuần hoàn và sự
thay đổi phân cực của nơtron trong tinh thể cũng đã được nghiên cứu [9,11,13].
Trong bài luận văn này, chúng tôi nghiên cứu: Tán xạ hạt nhân của các nơtron
phân cực trên mặt tinh thể có các hạt nhân phân cực trong điều kiện có phản
xạ
Một phần kết quả của luận văn đã được báo cáo tại hội nghị vật lý lý thuyết
toàn quốc lần thứ 36 tổ chức tại thành phố Quy Nhơn tháng 8 năm 2011.
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
4
Nội dung của luận văn được trình bày trong 4 chương:
Chương 1 - Lý thuyết tán xạ của nơtron chậm trong tinh thể
Chương 2 - Tiến động hạt nhân của spin của các nơtron trong môi trường
phân cực.
Chương 3 - Phản xạ gương của nơtron phân cực trên mặt biên gồ ghề
giữa chân không và vật chất có các hạt nhân phân cực
Chương 4 - Tán xạ hạt nhân của các nơtron phân cực trên mặt tinh thể có các
hạt nhân phân cực trong điều kiện có phản xạ
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
5
CHƯƠNG 1 – LÝ THUYẾT TÁN XẠ CỦA NƠTRON CHẬM
TRONG TINH THỂ
1.1. Cơ sở lý thuyết tán xạ của nơtron chậm trong tinh thể
Trong trường hợp khi bia tán xạ cấu tạo từ số lớn các hạt (ví dụ như tinh thể),
để tính toán tiết diện tán xạ một cách thuận tiện ta đưa vào lý thuyết hình thức luận
thời gian
Giả sử ban đầu bia được mô tả bởi hàm sóng
n
, là hàm riêng của toán tử
Hamilton của bia
H
n
=E
n
n
(1.1.1)
Sau khi tương tác với nơtron sẽ chuyển sang trạng thái
'n
. Còn nơtron có thể
thay đổi xung lượng và spin của nó. Giả sử ban đầu trạng thái của nơtron được mô
tả bởi hàm sóng
p
. Ta đi xác định xác suất mà trong đó nơtron sau khi tương tác
với hạt nhân bia sẽ chuyển sang trạng thái
'p
và hạt bia chuyển sang trạng thái
'n
Xác suất W
n’p’|np
của quá trình đó được tính theo lý thuyết nhiễu loạn trong gần
đúng bậc nhất sẽ bằng :
( )
2
' '| ' '
2
' '
n p np n p n p
W n p V np E E E E
= + − −
(1.1.2)
Trong đó:
V là toán tử tương tác của nơtron với hạt nhân bia.
' '
, , ,
n p n p
E E E E
là các năng lượng tương ứng của hạt bia và nơtron trước và sau
khi tán xạ.
( )
' 'n p n p
E E E E
+ − −
- hàm delta Dirac.
( )
( )
' '
' '
1
2
n p n p
i
E E E E t
n p n p
E E E E e dt
+∞
− + − −
−∞
+ − − =
∫
(1.1.3)
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
6
Chúng ta quan tâm tới xác suất toàn phần W
p’|p
của quá trình trong đó nơtron
sau khi tương tác với bia sẽ chuyển sang trạng thái
p
; nó nhận được bằng cách
tổng hóa các xác suất W
n’p’|np
theo các trạng thái cuối của bia và lấy trung bình theo
các trạng thái đầu. Bởi vì bia không luôn ở trạng thái cố định do đó ta phải tổng
quát hóa đối với trường hợp khi nó ở trong trạng thái hỗn tạp với xác suất của trạng
thái
n
là
n
. Theo đó ta có:
( )
2
'| ' '
'
2
W ' '
p p n n p n p
nn
n p V np E E E E
= + − −
∑
( )
2
' ' '
'
2
'
n p p n p n p
nn
n V n E E E E
= + − −
∑
(1.1.4)
Ở đây chúng ta đưa vào kí hiệu hỗn hợp để cho các yếu tố ma trận
'
' ' '
p p
n p V np n V n=
(1.1.5)
Như vậy là các yếu tố ma trận của toán tử tương tác của nơtron với hạt bia lấy
theo các trạng thái của nơtron và V
p’p
là toán tử tương đối với các biến số hạt bia
Thay phương trình (1.1.3) vào (1.1.4) ta được:
( )
( )
' '
*
'| ' ' '
2
'
1
W ' '
p p n n
i i
E E t E E t
p p nn p p p p
nn
e dt n V n n V n e
+∞
− −
−∞
=
∑
∫
(1.1.6)
E
n
, E
n’
là các trị riêng của toán tử Hamilton H với các hàm riêng là
n
,
'n
, từ
đó ta viết lại trong biểu diễn Heisenberg:
( )
( )
'
' '
' '
n n
i
E E t
p p p p
n V n e n V t n
−
=
(1.1.7)
Ở đây:
( )
' '
i i
Ht Ht
p p p p
V t e V e
−
=
là biểu diễn Heisenberg của toán tử V
p’p
với toán
tử Hamilton.
Thay (1.1.7) vào (1.1.6), chú ý rằng trong trường hợp này ta không quan tâm
tới sự khác nhau của hạt bia trước và hạt bia sau tương tác, vì vậy công thức lấy
tổng theo n’, n chính là vết của chúng và được viết lại:
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
7
( )
( )
'
'| ' ' '
2
'
1
W '
p p
i
E E t
p p nn p p p p
nn
e dt n V V t n
+∞
−
+
−∞
=
∑
∫
( )
( )
{ }
'
' '
2
1
p p
i
E E t
p p p p
dte Sp V V t
+∞
−
+
−∞
=
∫
(1.1.8)
Ở biểu thức cuối, biểu thức dưới dấu vết có chứa toán tử thống kê của bia
,
các phần tử đường chéo của ma trận của nó chính là xác suất
n
.
Theo qui luật phân bố Gibbs nếu hạt bia nằm ở trạng thái cân bằng nhiệt động
ta có hàm phân bố trạng thái là:
{ }
H
H
e
Sp e
−
−
=
Với:
1
z
k T
=
z
k
- hằng số Boltmann
T - Nhiệt độ
Giá trị trung bình thống kê của đại lượng Vật lý được tính theo các hàm phân
bố là:
{ }
{ }
H
n
H
n
Sp e A
A A
Sp e
−
−
= =
∑
(1.1.9)
Kết hợp (1.1.8) và (1.1.9) ta được:
( )
( )
{ }
'
'| ' '
2
1
W
p p
i
E E t
p p p p p p
dte Sp V V t
+∞
−
+
−∞
=
∫
( )
( )
{ }
{ }
'
' '
2
1
p p
H
i
E E t
p p p p
H
Sp e V V t
dte
Sp e
− +
+∞
−
−
−∞
=
∫
( )
( )
'
' '
2
1
p p
i
E E t
p p p p
dte V V t
+∞
−
+
−∞
=
∫
(1.1.10)
Nếu chuẩn hóa hàm sóng của nơtron trên hàm đơn vị ( trên hàm
) thì tiết
diện tán xạ hiệu dụng được tính trên một đơn vị góc cầu và một khoảng đơn vị năng
lượng
2
d
d dE
Ω
, sẽ liên quan tới xác suất này bởi biểu thức sau:
( ) ( )
( )
( )
'
2 2 2
'|
' '
3 3
5
'
' '
W
2 2
p p
i
E E t
p p
p p p p
p
d m p m p
dte V V t
d dE p p
+∞
−
+
−∞
= =
Ω
∫
(1.1.11)
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
8
Gạch trên đầu là trung bình theo các trạng thái spin của nơtron trong chùm các
nơtron ban đầu và tổng hóa các trạng theo các trạng thái spin trong chùm tán xạ
m - khối lượng nơtron
Trong công thức (1.1.11) đưa vào toán tử mật độ spin của nơtron tới
và sử
dụng công thức:
{ }
L Sp L
=
(1.1.12)
Do đó dạng tường minh của công thức (1.1.11) được viết lại là:
( )
( )
( )
{ }
'
2 2
' '
3
5
'
'
2
p p
i
E E t
p p p p
p
d m p
dte Sp V V t
d dE p
+∞
−
+
−∞
=
Ω
∫
(1.1.13)
Trong đó:
- ma trận mật độ spin nơtron
1.2. Thế tương tác của nơtron chậm trong tinh thể
Tán xạ của nơtron chậm khi đi vào mạng tinh thể sẽ chịu tác động của tương
tác hạt nhân và tương tác từ.
1.2.1. Yếu tố ma trận của tương tác hạt nhân
Ta xây dựng thế hạt nhân của nơtron và hạt nhân bia dưới dạng sau:
( ) ( )n nV r r R
= −
(1.2.1)
Trong đó
( )A B sJ
= +
(1.2.2)
nr
- vị trí của nơtron
R
- Vị trí của hạt nhân
,A B
- là các hằng số
J
- Spin của hạt nhân
s
- Spin của nơtron
Do đó thế tương tác của nơtron với hạt nhân thứ l là:
( ) ( )
l n n l
V r r R
= −
(1.2.3)
Lấy tổng công thức (1.2.3) theo l từ 1 đến số hạt nhân trong bia ta sẽ tìm được
thế tương tác của nơtron với toàn bộ bia:
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
9
'
1
l
N
iqR
p p l
l
V e
=
=
∑
(1.2.4)
Các yếu tố ma trận
'p p
V
thuộc toán tử tương tác hạt nhân V từ xung lượng
p
đến
'p
được ghi nhận trên cơ sở (1.2.3) có dạng:
( )
( )
'
l
iqR
l l l
p p
l
V A B sJ e= +
∑
(1.2.5).
Trong đó
'q p p= −
: Véctơ tán xạ của nơtron
1.2.2. Yếu tố ma trận của tương tác từ.
Tương tác từ của nơtron trong mạng tinh thể xuất hiện do các điện tử tự do
chuyển động. Và bản thân nơtron cũng có mômen từ sinh ra. Thế đặc trưng cho
tương tác này được cho bởi biểu thức [21]
( )
2
' 0
4
( ) , ( )
j
iqR
p p j j
j
V r F q e S s es e
m
= − −
∑
(1.2.6).
Trong đó:
2
0
2
0
e
r
m c
=
: là véctơ bán kính điện từ của electron
m - khối lượng nơtron
1.913
= −
- độ lớn mômen từ hóa trên manhêton Bohr hạt nhân
j
S
- Spin của nguyên tử thứ j
l
R
- là véctơ tọa độ vị trí hạt nhân thứ l
'q p p= −
- véctơ tán xạ của nơtron
q
e
q
=
- véctơ tán xạ đơn vị
s
- spin của nơtron tới
Biểu thức
( )
( )
*
( )
1
j
iqr
Z
j
j j j j
j j
e s S
F q d
S S
=
+
∑
∫
(1.2.7)
Với
j
z
j
S s
=
∑
là toán tử spin của nguyên tử thứ j
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
10
j
S
là đại lượng spin của nguyên tử thứ j
j
là hàm sóng của điện tử thứ j
( )
j
F q
đặc trưng cho sự phân bố của mật độ spin trong nguyên tử
Như vậy khi xét bài toán của một chùm nơtron chậm không phân cực tán xạ
trong tinh thể, ngoài tương tác hạt nhân chúng còn tương tác từ. Do đó trong biểu
thức tiết diện tán xạ vi phân sẽ gồm đóng góp hai phần được đặc trưng bởi hai loại
tương tác ở trên
2 2
2
' ' '
n m
p p p
d d
d
d dE d dE d dE
= +
Ω Ω Ω
(1.2.8)
Do vậy đại lượng
'p p
V
được viết lại dưới dạng sau:
( )
( )
( )
'
2
0
4
( ) , ( )
l
j
iqR
iqR
l l l j j
p p
l j
V A B sJ e r F q e S s es e
m
= + − −
∑ ∑
(1.2.9)
Từ đó ta đi tính được tiết diện tán xạ vi phân
( )
( )
( )
{ }
'
2 2
' '
3
5
'
'
2
p p
i
E E t
p p p p
p
d m p
dte Sp V V t
d dE p
+∞
−
+
−∞
=
Ω
∫
(1.2.10)
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
11
CHƯƠNG 2 – TIẾN ĐỘNG HẠT NHÂN CỦA SPIN CỦA
CÁC NƠTRON TRONG MÔI TRƯỜNG PHÂN CỰC
2.1. Tính góc tiến động bằng phương pháp toán tử.
Giả sử hạt tới và bia đều có spin. Chúng ta xem xét quá trình chuyển động
của nơtron chậm qua vật chất.
Trong trường hợp này, hàm sóng mô tả quá trình va chạm đàn hồi của nơtron
với hạt nhân được gắn ở điểm
i
R
có dạng :
( )
i
i
ik r R
ik r ik R
n nuc n nuc
i
e
r e f e
r R
−
= +
−
(2.1.1)
Trong đó:
n
là hàm sóng spin của nơtron tới,
nuc
là hàm sóng spin của hạt
nhân
Trong trường hợp các nơtron chậm, bước sóng lớn hơn nhiều so với kích
thước của hạt nhân, vì vậy biên độ tán xạ không phụ thuộc vào góc tán xạ và có thể
được viết dưới dạng:
f J
= +
(2.1.2)
Trong đó:
2S
=
,
S
là toán tử spin của nơtron
là toán tử ma trận được tạo bởi các ma trận Pauli
J
là toán tử spin của hạt nhân
1
2 1 2 1
J J
a a
J J
+ −
+
= +
+ +
và
2 1
a a
J
+ −
−
=
+
a
+
là biên độ tán xạ trong trạng thái ứng với mômen tổng cộng của nơtron và
hạt nhân là
1
2
J +
a
−
là biên độ tán xạ trong trạng thái ứng với mômen tổng cộng của nơtron và
hạt nhân là
1
2
J −
Trong quá trình tán xạ trên hệ hạt nhân, hàm sóng viết lại có dạng sau:
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
12
( )
. .
i
i
ik r R
ik r ik R
n nuc m n nuc m
i
i
m m
i
e
r e f e
r R
−
= +
−
∑
∏ ∏
(2.1.3)
Trong đó
.nuc m
m
∏
là hàm sóng spin của các hạt nhân với giả thiết rằng các hạt
nhân không tương tác với nhau.
Để tìm sóng kết hợp trong trường hợp này, chúng ta làm trung bình cộng công
thức (2.1.3) theo phân bố của các hạt nhân bia và theo các trạng thái spin của chúng
Sự trung bình hóa đó dẫn đến biểu thức sau của hàm sóng:
( )
i
i
ik r R
ik r ik R
n m n
i
i
e
r e f e
r R
−
= +
−
∑
(2.1.4)
Trong đó:
f J pJ
= + = +
J
p
J
=
: Véctơ phân cực của hạt nhân
J: spin của hạt nhân
Nếu các hạt nhân được phân bố hỗn loạn trong mặt phẳng
0
z z=
thì chúng ta
sẽ nhận được biểu thức sau cho sóng kết hợp đi qua mặt phẳng trên:
( )
2 '
1 ( )
ik r
n
z
i
r J p e
k
= + +
(2.1.5)
Trong công trình [16], toán tử
1
2
p
i
B
= +
(2.1.6)
được gọi là toán tử spin quay xung quanh một trục đặc trưng bởi vectơ đơn vị
p
một góc
;
<<1
So sánh (2.1.5) và (2.1.6) ta có thể kết luận: Sau khi đi qua mặt phẳng phân
cực, spin của nơtron đã quay đi 1 góc:
( )
4 '
Re
z
Jp
k
=
(2.1.7)
Nếu hàm sóng của các nơtron đi qua m mặt thì góc quay tổng cộng là:
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
13
( )
4 '
Re
m
z
Jpm
k
=
(2.1.8)
Hay, khi đi qua một tấm bia có độ dày l xác định, chúng ta sẽ thu được: Khi
nơtron đi qua bia phân cực này, spin của nó sẽ quay đi góc:
( )
4
Re
z
Jpl
k
=
(2.1.9)
Kết quả này có thể nhận được bằng các phương pháp khác
2.2. Tính góc tiến động bằng phương pháp hàm sóng
Chọn trục lượng tử song song với véctơ phân cực của hạt nhân
p
. Nếu nơtron
tới mặt phẳng có spin song song với véctơ
p
(
1
0
n
=
), thì sóng kết hợp
( )r
+
có dạng:
( )
1
2 '
1
0
ik r
z
i
r f e
k
+
+
= +
(2.2.1)
Trong đó:
f Jp
+
= +
là biên độ tán xạ kết hợp đàn hồi dưới góc bằng 0 của
nơtron với spin song song với véctơ phân cực của hạt nhân p
Đối với nơtron có spin ngược lại thì sóng kết hợp đàn hồi
( )
r
−
có dạng:
( )
0
2
1
1
ik r
z
i
r f e
k
−
−
= +
(2.2.2)
Trong đó:
f Jp
−
= −
là biên độ tán xạ kết hợp đàn hồi dưới góc bằng 0 của nơtron với
spin phản song song với véctơ phân cực của hạt nhân p
Nếu hàm sóng đi qua một lớp vật chất có độ dày xác định thì lặp lại tất cả các
lý luận dẫn đến biểu thức của hệ số khúc xạ đối với bia phân cực mà ta đã biết thì
chúng ta sẽ nhận được hệ số khúc xạ của các nơtron có spin song song với véctơ p
như sau:
( )
2 2
2 2
1 1
z z
n f Jp
k k
+ +
= + = + +
(2.2.3)
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
14
Đối với các nơtron với sự phân cực ngược lại thì:
( )
2 2
2 2
1 1
z z
n f Jp
k k
− −
= + = + −
(2.2.4)
Hiệu số
( )
2
2
z
n n n f f
k
+ − + −
∆ = − = −
(2.2.5)
được xác định bởi hiệu các biên độ tán xạ của sóng kết hợp tương ứng và khác 0 chỉ
trong bia phân cực
Như vậy, trong hạt nhân bia phân cực, nơtron có 2 hệ số khúc xạ
Xét trường hợp nơtron có véctơ phân cực tạo thành một góc tương đối với
hướng của véctơ phân cực hạt nhân. Chọn một hướng của p tạo thành một góc
tương đối với trục z. Véctơ phân cực của hạt nhân bia có phương vuông góc với bề
mặt.
Hàm sóng cơ sở có dạng:
( )
ik r
n
r e
=
,
1
2
n
c
c
=
(2.2.6)
Hay:
1 2
1 0
( )
0 1
ikr ikr
r c e c e
= +
Trạng thái spin
1
0
có liên quan tới chỉ số khúc xạ
n
+
Trạng thái spin
0
1
có liên quan tới chỉ số khúc xạ
n
−
Hàm sóng của nơtron trong trạng thái phân cực thay đổi theo chiều sâu xác
định theo biểu thức sau:
1
1 2
2 _
( )
1 0
( )
0 1
( )
z z
ik n z ik n z
ik r ik r
c r
r c e e c e e
c r
⊥ ⊥
⊥ ⊥
+ −
+
= = +
(2.2.7)
Véctơ phân cực của nơtron là :
n
P
=
(2.2.8)
có các thành phần là :
* *
1 2
2Re
nx
P c c
+ −
=
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
15
* *
1 2
2Im
ny
P c c
+ −
=
(2.2.9)
2 2
1 2nz
P c c
+ −
= −
Giả thiết rằng spin của nơtron có phương vuông góc với vecto phân cực của
hạt nhân và có phương song song với trục x, ta được
1 2
1
2
c c= =
Sử dụng các đẳng thức (2.2.9) ta có :
( )
( )
Im
os Re
z
k n n z
nx z
P c k n n z e
+ −
− +
+ −
= −
( )
( )
Im
sin Re
z
k n n z
ny z
P k n n z e
+ −
− +
+ −
= − −
(2.2.10)
2 Im 2 Im
z z
k n z k n z
nz
P e e
+ −
− −
= −
Suy ra, vectơ phân cực của nơtron hợp với vectơ phân cực của hạt nhân một
góc :
( ) ( )
2
Re Re
z
z
k n n z f f z
k
+ − + −
= − = −
(2.2.11)
Biểu thức của (2.2.11) phù hợp với (2.2.10).
Trong trường hợp tổng quát, vectơ phân cực của hạt nhân không xác định. Để
mô tả hiệu ứng quay của spin nơtron ta dùng toán tử quay spin đi một góc
nào đó.
Sử dụng (2.1.5) ta có : Sau khi đi qua m mặt phẳng phân cực, hàm sóng của
nơtron là :
( ) ( )
2 '
1
m
z
i
r J p
k
= + +
ik r
n
e
(2.2.12)
Sau khi nơ tron đi qua lớp vật chất có bề dày là z = ma (a- bề dày của 1 lớp)
thì
( )
r
được viết như sau :
( )
z
ik r ik nz
n
r e e
⊥ ⊥
=
(2.2.13)
Với :
( )
2
2
1 0
z
n f
k
= +
(2.2.14)
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
16
( )
0f J p
= +
là biên độ tán xạ đàn hồi trên hạt nhân với một góc bằng 0.
So sánh với việc mô tả bằng toán tử quay spin của nơtron đi một góc trong [16]
2
i n
B e
=
(2.2.15)
ta thấy, trong trường hợp này, toán tử quay spin nơtron được mô tả bởi :
( )
2
exp Re
n p
z
B i Jp n z
k
=
(2.2.16)
Trong đó
p
p
n
p
=
Ngoài ra, sự quay spin của nơtron trong bia phân cực có thể nhận được bằng
cách khác.
2.3. Sử dụng bảo toàn năng lượng để tính góc tiến động.
Gọi năng lượng của sóng kết hợp là
'
kh
E
Năng lượng của sóng tự do trong chân không là
ck
E
Theo định luật bảo toàn năng lượng thì thế năng có dạng :
( )
2 2
2
' 2
2
(1 ) 0
2
z
ck kh
k
U E E n f
m m
= − = − = −
(2.3.1)
Như vậy trong hạt nhân bia phân cực, sóng nơtron có khả năng khúc xạ với
các mức năng lượng là :
( )
2
2 2
2
2
(1 ) 0
2
z
k
U n f
m m
± ± ±
= − = −
(2.3.2)
So sánh với (2.2.14) ta viết lại năng lượng dưới dạng toán tử :
( )
2 2
2 2
0 ( )U f J p
m m
= − = − +
(2.3.3)
Khi nơtron chuyển động trong từ trường, năng lượng tương tác của thành phần
spin song song với
H
được tính theo công thức :
W H
+
= −
Tương tự với thành phần spin ngược lại ta có năng lượng bằng
W H
−
=
Hiệu năng lượng là :
W W 2 H
+ −
− = −
Giới hạn của tần số chuyển động tiến động của nơtron trong từ trường H là :
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
17
2
R
H
=
Hoàn toàn tương tự, trong từ trường tồn tại hiệu số thế
U U
+ −
−
, spin của
nơtron chuyển động tiến động quanh trục song song với vectơ phân cực của hạt
nhân với tần số :
4
Re Re
U U
Jp
m
+ −
−
= =
(2.3.4)
Trong khoảng thời gian t, spin của nơtron quay đi một góc
t
=
.
Nếu phần có từ trường có độ dài l, thời gian để nơtron đi qua là :
z
l
t
v
=
Vậy spin của nơtron quay đi một góc :
2
Re
z z z
l m l
Jpl
v k k
= = =
Điều này hoàn toàn phù hợp với công thức (2.1.9).
Trong từ trường thì tương tác giữa spin của nơtron với hạt nhân có từ trường
hiệu dụng :
e ff
2
H
=
Tương tự, nếu như bia phân cực có từ trường phụ thuộc vào thời gian B(t) và
vectơ phân cực của hạt nhân cũng phụ thuộc vào thời gian
( )p p t=
thì từ trường
hiệu dụng tổng hợp là :
( ) ( ) ( )
eff
G t B t H t= +
Trong đó
2
2
( ) ( )
eff
J
H t p t
m
=
Như vậy, năng lượng tương tác spin trong từ trường hiệu dụng là :
( )
eff
( ) ( )V G B t H t
= − = − +
(2.3.5)
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
18
CHƯƠNG 3 – PHẢN XẠ GƯƠNG CỦA CÁC NƠTRON
PHÂN CỰC TRÊN MẶT BIÊN GỒ GHỀ GIỮA “CHÂN
KHÔNG – VẬT CHẤT” CÓ CÁC HẠT NHÂN PHÂN CỰC
3.1. Ảnh hưởng của sự gồ ghề mặt biên “chân không – vật chất” có các hạt
nhân phân cực lên phản xạ gương của các nơtron phân cực
Phản xạ gương của nơtron trên mặt biên giữa vật chất và chân không đã được
nghiên cứu [19]. Sự xuất hiện gồ ghề của mặt biên giới hạn đã dẫn tới sự phụ thuộc
của hệ số phản xạ vào hệ số Debye-Waller [15]. Sự khác biệt giữa công thức mô tả
sự phản xạ gương trên mặt biên phẳng với công thức trong trường hợp có sự gồ ghề
cho phép phán đoán trạng thái bề mặt
Khi xem xét phản xạ gương của các nơtron phân cực trên biên thực tế giữa vật
chất và chân không, chúng ta cần tính đến sự gồ ghề của mặt biên. Sự gồ ghề của
mặt biên thực xuất hiện là do sự gồ ghề của các vị trí của các hạt nhân trong quá
trình dao động nhiệt hoặc là do sự thăng giáng vị trí của biên đến cỡ vài chục
0
A
Giả sử chùm nơtron phân cực tiến đến bề mặt của vật chất có các hạt nhân
phân cực nằm chiếm nửa không gian x >0
Trong bia phân cực như chúng ta biết [18] từ trường tổng cộng hiệu dụng
eff
G
sẽ tác động lên chùm nơtron
eff
G B=
+
eff
nuc
H
(3.1.1)
Ở đó
B
- vectơ cảm ứng từ.
eff
nuc
H
- từ trường hiệu dụng hạt nhân
Chúng ta giả thiết rằng trong nửa không gian x>0, trong vật chất có các hạt
nhân phân cực có từ trường hiệu dụng đồng nhất có dạng:
eff eff
0;
x y
G G= =
eff eff
( )
z
G G x=
Trục z có hướng song song với mặt của bia
Trong trường hợp này quá trình phản xạ, khúc xạ của các nơtron phân cực trên
bia được xác định bởi Hamiltonien
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
19
H=
2
2
p
m
+
eff
( ) ( )
z
V x G x
−
(3.1.2)
Ở đó, p, m- là toán tử xung lượng và khối lượng của nơtron
- moment từ của nơtron
( )V x
: Thành phần thế hạt nhân hiệu dụng không phụ thuộc vào spin
1 0
0 1
z
=
−
- ma trận Pauli.
Ta viết lại (3.1.2) dưới dạng:
H =
0
( , )
z
H x
+
(3.1.3)
0
H
=
2
2
p
m
+
[ ]
0 eff
( )
z
V G x
−
Trong đó:
0
V
và
eff
G
- là các giá trị của
( )V x
và
eff
( )G x
ở sâu trong bia cách xa
biên
( , )
z
x
=
[ ]
0
( ) ( )V x V x
−
[ ]
eff eff
( ) ( )
z
G x G x
− −
ở đó
( )x
=
1
0
, 0
, 0
x
x
>
<
( , )
z
x
- nhiễu loạn xuất hiện khi ta tính đến sự gồ ghề của mặt vật chất
Chúng ta sẽ đi thu nghiệm của phương trình Schrodinger
[ ]
0
( , )
z
H H x E
= + =
(3.1.4)
Dưới dạng sau :
||
||
( )
z
ik r
S
e x
=
Ở đó
z
S
- hàm spin tương ứng với giá trị xác định
z
S
của hình chiếu của spin
của nơtron lên trục z:
1
2
z z
z S z S
S
= ±
,
1
2
z
S =
||
k
và
||
r
- là các thành phần của vectơ sóng và vectơ vị trí của nơtron song
song với bề mặt của vật chất
Đặt (3.1.2) vào (3.1.4) chúng ta sẽ nhận được phương trình để cho
( )x
có
dạng
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
20
[ ]
2
0 eff 1
2
2
( ) ( ) ( ) ( ) 0
x x
m
x k V G x x x
± ± ±
∆ + − + =
(3.1.5)
Ở đó
1/2
2
2
0
x
mE
k
⊥
= >
1
2
2
( ) ( )
m
x x
± ±
= −
[ ] [ ]
0 eff eff
( ) ( ) ( ) ( ) ( )x V x V x G x G x
±
= − −
2 2
0
2
y z
P P
E E
m
⊥
+
= −
- năng lượng chuyển động dọc của nơtron.
Nhờ hàm Green của phương trình Schrodinger mô tả phản xạ gương trên biên
phẳng.
[ ]
2
0 eff
2
2
( , ') ( ) ( , ') ( ')
x x
m
G x x k V G x G x x x x
± ±
∆ + − = −
(3.1.6)
Chúng ta biểu diễn phương trình (3.1.5) trong dạng tích phân:
0 1
( ) ( ) ( , ') ( ') ( ') 'x x G x x x x dx
± ± ± ± ±
= +
∫
(3.1.7)
Ở đó
0
( )x
±
- nghiệm của phương trình thuần nhất xác định phản xạ gương
trên biên phẳng chân không- vật chất:
0
0
, 0
( )
, 0
x x
x
ik x ik x
o
ik x
e A e x
x
B e x
< <
>
±
−
±
±
−
±
+ <
=
>
Ở đó:
2
2
0
x
mE
k
<
⊥
= >
[ ]
0 eff
2
2
x
m
k E V G
>
± ⊥
= − ±
> 0
Từ điều kiện liên tục của hàm sóng và của đạo hàm của hàm sóng trên biên
x= 0 chúng ta xác định được các hệ số của sóng phản xạ và sóng khúc xạ:
0 0
1 A B
± ±
+ =
,
0 0
(1 )
x x
k A k B
< >
± ± ±
− =
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
21
0
x x
x x
k k
A
k k
< >
±
±
< >
±
−
=
+
;
0
2
x
x x
k
B
k k
<
±
< >
±
=
+
(3.1.8)
Để tìm biên độ của sóng phản xạ gương chúng ta cần nghiên cứu tiệm cận của
hàm sóng (3.1.7) khi
x → −∞
. Có thể chỉ ra rằng:
0
lim ( , ') ( ')
x
ik x
x
x
i
G x x e x
k
<
<
→−∞
≈
(3.1.9)
Ở đó
0
( ')x
- nghiệm của phương trình thuần nhất xác định phản xạ gương trên
biên phẳng của chân không - vật chất.
Thay (3.1.9) vào (3.1.7) chúng ta sẽ nhận được biên độ sóng phản xạ có dạng
như sau:
0 0 1
( ') ( ') ( ') '
x
i
A A x x x dx
k
± ± ± ± ±
<
= +
∫
(3.1.10)
Hạn chế ở gần đúng bậc nhất và chú ý đến các công thức (3.1.8) chúng ta sẽ
nhận được:
2 ' 2 '
2
0 0 0 1
2 ( ') '
x x
ik x ik x
x
i
A A e A e A x dx
k
< <
+∞
−
± ± ± ± ±
<
−∞
≈ + + +
∫
( )
2
2
0 0 0 1
1 2 '(1 ) ( ') '
x
x
i
A A ik x A x dx
k
+∞
<
± ± ± ±
<
−∞
≈ + + + −
∫
2
0
0 1 1
( ') ' 2 ' ( ') '
x
x
iB
A x dx ik x x dx
k
+∞ +∞
>
±
± ± ± ±
<
−∞ −∞
≈ + +
∫ ∫
Nếu
1
( ')x
±
là một hàm chẵn thì tích phân thứ hai của biểu thức trên sẽ bằng
không và ta có
2
0
0 1
0
2 ( ') '
x
iB
A A x dx
k
+∞
±
± ± ±
<
≈ +
∫
(3.1.11)
Chúng ta xét một ví dụ khi
1
( ')x
±
có dạng Gauss :
2
2
0
2
1 0
( ')
x
d
x e
−
± ±
=
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
22
Ở đó
0
d
- biên độ đặc trưng của sự gồ ghề. Thay
1
( ')x
±
vào (3.1.11) và tính
tích phân ta sẽ nhận được :
2
2
0
'
2
2
0
0 0
2
0
2 2
'
x
d
x
iB m
A A e dx
k
−
∞
±
± ± ±
<
≈ −
∫
2
0 0 0
0
2
2 2
x
i B d m
A
k
± ±
±
<
= −
0 0
0
2 2
8 2
( )
x
x x
i k md
A
k k
<
±
±
< >
±
= −
+
(3.1.12)
Như vậy cường độ của sóng phản xạ được xác định bởi biểu thức sau :
2
0 0 0
0
2 2
2
16Im
( )
x
x x
A k md
J A
k k
<
± ±
± ±
< >
±
≈ +
+
(3.1.13)
Bây giờ chúng ta đánh giá số hạng bộ xung vào cường độ của sóng phản xạ ở
gần góc tới hạn đặc trưng có sự gồ ghề của bề mặt biên. Để làm được điều đó chúng
ta chọn
k
<
9
10≈
cm
1−
và góc trượt của nơtron
0
0,1
=
.
Trong trường hợp đó
6
10
x
k
<
≈
cm
1−
Theo kết quả của [18] thì
2
0 0
2
(0)V f
m
±
, ở đó
- mật độ hạt nhân, f(0) – biên độ tán xạ về phía trước
của nơtron. Nếu chọn
~
10
22
cm
3−
,
(0)f
10
12−
cm,
7
0
10d cm
−
thì :
2 1
0 0 0
2 2
16 2
10 10
( )
x
x x
A k md
k k
<
− −
± ±
< >
±
−
+
Như vậy chúng ta đã thấy phần đóng góp bổ sung vào cường độ của sóng phản
xạ của nơtron đặc trưng cho sự gồ ghề của bề mặt biên là không nhỏ ngay cả khi d
0
rất nhỏ và bằng
7
10 cm
−
.
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
23
3.2. Vectơ phân cực của nơtron phản xạ gương trên mặt biên gồ ghề giữa chân
không và vật chất có các hạt nhân phân cực.
Ta xét ảnh hưởng của sự gồ ghề của mặt biên tới trạng thái của vectơ phân cực
của nơtron phản xạ.
Véctơ phân cực của nơtron phản xạ được xác định bởi công thức :
px px
px px
P
=
(3.2.1)
Giả sử rằng các nơtron tiến đến bia có các vectơ phân cực hướng theo một
góc nào đó đối với hướng của vectơ phân cực của hạt nhân bia
N
P
.
Trạng thái của nơtron có thể xem như là sự tổ hợp của hai trạng thái phân
cực, phân cực theo vectơ phân cực của hạt nhân bia
N
P
và phân cực theo hướng
ngược lại.
Hàm sóng mô tả trạng thái spin của nơtron tới là :
1 2
1 0
0 1
z
s
c c
= +
Trong đó
2
1
c
và
2
2
c
cho ta xác suất tìm thấy nơtron có trạng thái spin
1
2
z
S =
và
1
2
z
S = −
.
Ta xem xét hàm sóng phản xạ của nơtron có dạng như sau :
( ) ( )
//
//
1
0 0 0 0
0 0
2 2
2 2
2
0
8 2 8 2
0
x
ik x
ik r
x x
px
x x x x
c
i k m d i k m d
e A A e
c
k k k k
<
< <
−
+ −
+ −
< > < >
+ −
= − + −
+ +
(3.2.2)
Đặt :
( )
0
2
2
8 2
x
x x
i k m
k k
<
+
< >
+
=
+
và
( )
0
2
2
8 2
x
x x
i k m
k k
<
−
< >
−
=
+
Thay vào (3.2.2) ta có thể viết sóng phản xạ như sau :
[ ]
[ ]
//
//
1 0 0
1 0 0
x
ik x
ik r
px
c A d
e e
c A d
<
+
−
−
−
=
−
(3.2.3)
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
24
Thay (3.2.3) vào (3.2.1) và lưu ý các ma trận Pauli :
0 1
1 0
x
=
;
0
0
y
i
i
−
=
;
1 0
0 1
z
=
−
Ta có :
px x px
x
px px
P
=
0 1
1 0
px x px px px
= +
( )
[ ]
[ ]
1 0 0
* * * * * *
1 0 0 2 0 0
2 0 0
0 1
,
1 0
x x
ik x ik x
ik r ik r
c A d
e c A d c A d e e e
c A d
< <
+
−
−
+ −
−
−
= − −
−
( )
[ ]
[ ]
2 0 0
* * * * * *
1 0 0 2 0 0
1 0 0
, .
c A d
c A d c A d
c A d
−
+ −
+
−
= − −
−
* * * * * *
1 0 0 2 0 0 2 0 0 1 0 0
. .c A d c A d c A d c A d
+ − − +
= − − + − −
Bỏ qua các số hạng chứa
2
0
d
, ta có :
px x px
( ) ( )
* * * * * * * *
1 2 0 0 0 0 0 0 1 2 0 0 0 0 0 0
c c A A d A A d c c A A d A A d
+ − − + − + − +
= − − + − −
* * * *
1 2 0 0 0 0 0 0
2Re c c A A d A A d
+ − − +
= − −
px px px px
+
=
( )
[ ]
[ ]
1 0 0
* * * * * *
1 0 0 2 0 0
1 0 0
, .
x x
ik x ik x
ik r ik r
c A d
e c A d c A d e e e
c A d
< <
+
−
−
+ −
−
−
= − −
−
( ) ( )
* * * * * * * *
1 1 0 0 0 0 0 0 2 2 0 0 0 0 0 0
c c A A d A A d c c A A d A d A
+ + + + − − − −
= − − + − −
( ) ( )
2 2 2 2
* *
1 0 0 0 2 0 0 0
2Re 2Rec A A d c A A d
+ + − −
= − + −
Vậy :
Luận văn thạc sĩ khoa học
Vũ Thị Thu Trang
25
( )
( ) ( )
* * * *
1 2 0 0 0 0 0
2 2 2 2
* *
1 0 0 0 2 0 0 0
2Re
2Re 2Re
x
c c A A A A d
P
c A A d c A A d
+ − − +
+ + − −
− +
≈
− + −
(3.2.4)
Tính toán tương tự cho
y
P
,
z
P
(bỏ qua các số hạng chứa
2
0
d
), ta có :
px y px
y
px px
P
=
( )
( ) ( )
* * * *
1 2 0 0 0 0 0
2 2 2 2
* *
1 0 0 0 2 0 0 0
2Im
2Re 2Re
c c A A A A d
c A A d c A A d
+ − − +
+ + − −
− +
≈
− + −
(3.2.5)
px z px
z
px px
P
=
( ) ( )
( ) ( )
2 2 2 2
* *
1 0 0 0 2 0 0 0
2 2 2 2
* *
1 0 0 0 2 0 0 0
2Re 2Re
2Re 2Re
c A A d c A A d
c A A d c A A d
+ + − −
+ + − −
− − −
≈
− + −
(3.2.6)
Từ các biểu thức nêu trên, ta thấy trạng thái phân cực của nơtron phản xạ
cũng phụ thuộc vào độ dày
0
d
của lớp chuyển tiếp ở bề mặt gồ ghề. Điều này cho
phép ta dựa vào các kết quả thực nghiệm đo vectơ phân cực của nơtron phản xạ để
nghiên cứu trạng thái gồ ghề của bề mặt vật chất.