( Tạp chí Khoa học Trường ĐHSP TPHCM ) ( Lý Duy Nhất và tgk ) ( TẠP CHÍ KHOA HỌC HO CHI MINH CITY UNIVERSITY OF EDUCATION TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINH JOURNAL OF SCIENCE Tập 17, Số 6 (2020) 1[.]
TẠP CHÍ KHOA HỌCHO CHI MINH CITY UNIVERSITY OF EDUCATION TRƯỜNG ĐẠI HỌC SƯ PHẠM TP HỒ CHÍ MINHJOURNAL OF SCIENCE Tập 17, Số (2020): 1019-1029 ISSN: 1859-3100 Vol 17, No (2020): 1019-1029 Website: Bài báo nghiên cứu* NĂNG LƯỢNG CHÍNH XÁC CAO CHO TRẠNG THÁI CƠ BẢN CỦA NGUYÊN TỬ HYDRO Ở MÔI TRƯỜNG PLASMA TRONG TỪ TRƯỜNG ĐỀU Lý Duy Nhất1,2*, Hoàng Đỗ Ngọc Trầm1 Trường Đại học Sư phạm Thành phố Hồ Chí Minh, Việt Nam Viện Năng lượng Nguyên tử Việt Nam * Tác giả liên hệ: Lý Duy Nhất – Email: nhatld@hcmue.edu.vn Ngày nhận bài: 31-3-2020; ngày nhận sửa: 21-5-2020, ngày chấp nhận đăng: 10-6-2020 TÓM TẮT Gần đây, phương pháp tốn tử FK áp dụng thành cơng cho toán nguyên tử hydro từ trường với cường độ Trong cơng trình này, kết phát triển cho trường hợp tổng qt có tính đến chắn Trong đo đạc vật lí thiên văn, nguyên tử hydro tồn thời gian ngắn mơi trường plasma mà chắn có vai trị quan trọng Phương pháp toán tử FK sử dụng với phép biến đổi tọa độ Kustaanheimo-Stiefel hàm sóng khai triển theo hàm sở dao động tử điều hòa bốn chiều Tuy nhiên, để tính yếu tố ma trận cho thành phần tương tác chắn, cần phát triển phương pháp đại số Điều làm báo với kết biểu thức giải tích cho yếu tố ma trận thu tường minh dạng chuỗi cho phép xây dựng chương trình FORTRAN tính tốn lượng với độ xác cao, lên đến 32 chữ số thập phân Giá trị số đưa minh họa báo với tham số chắn từ đến a.u cho trạng thái 1s Từ khóa: phương pháp đại số; nguyên tử hydro; chắn; hệ nguyên tử ba chiều; toán tử sinh; hủy; hàm sở Giới thiệu Nguyên tử hydro tốn đơn giản, có lời giải xác giải tích, nên thường dùng để kiểm chứng phát triển phương pháp khác trước áp dụng cho toán nguyên tử phức tạp Cho đến nay, có nhiều cơng trình nghiên cứu liên quan đến nguyên tử hydro, nguyên tử hydro từ trường (Cao, Ly, Hoang, & Le, 2019; Sasmal, 2014; Thirumalai, & Heyl, 2009), nguyên tử hydro điện trường (Bahar, & Soylu, 2014; Main, Schwacke, & Wunner, 1998) Trong năm qua, nguyên tử hydro plasma nghiên cứu nhiều tác giả (Bahar, & Soylu, Cite this article as: Ly Duy Nhat, & Hoang Do Ngoc Tram (2020) Highly accurate numerical energy for the ground state of a hydrogen atom in a plasma with the presence of a uniform magnetic field Ho Chi Minh City University of Education Journal of Science, 17(6), 1019-1029 Tạp chí Khoa học Trường ĐHSP TPHCM Tập 17, Số (2020): 10191029 2014, 2015, 2016; Kar, & Ho, 2005; Lin, & Ho, 2011; Paul, & Ho, 2009, 2010, 2011; Saha, Mukherjee, & Diercksen, 2002; Saha, Mukherjee, Mukherjee, Fricke, 2011; Soylu, 2012) Cơng trình Saha (2002) nghiên cứu đặc tính nguyên tử hydro plasma Debye phương pháp nhiễu loạn, mơ hình chắn Coulomb tĩnh (Static Screened Coulomb – SSC) sử dụng với dạng Yukawa cho phần chắn Sau đó, Paul (2009) dùng chắn SSC để tìm trị riêng lượng hàm sóng cho trạng thái phương pháp Ritz (Ritz variation method) Paul (2011) dựa mơ hình chắn Coulomb tổng qt dạng mũ cosine (Generalized Exponential Cosine Screened Coulomb – GECSC) để tính trị riêng lượng phương pháp Ritz Thế GECSC lấy từ mơ hình ECSC Shukla đề xuất năm 2008 (Shukla, & Eliasson, 2008) Sau đó, Soylu (2012) đề xuất mơ hình tổng quát (More Generalized Exponential Cosine Screened Coulomb – MGECSC) để nghiên cứu đặc tính nguyên tử hydro plasma Các kết tiếp tục tác giả phát triển có mặt điện trường cho mơ hình hai chiều (Bahar, & Soylu, 2014, 2015, 2016) Trong cơng trình nêu trên, tác giả giải số phương trình Schrưdinger chắn với độ xác đến chữ số thập phân Trong đó, lượng mức kích thích bậc cao gần nhau, đặc biệt xét đến chắn, cần phải giải số phương trình Schrưdinger với độ xác cao Phương pháp tốn tử FK giới thiệu chuyên khảo (Feranchuk, Ivanov, Le, & Ulyanenkov, 2015) có ưu điểm cho nghiệm hội tụ nhanh có độ xác cao Các ưu điểm thể cơng trình (Cao, Ly, & Hoang, 2016; Cao et al., 2019; Hoang, Nguyen, Hoang, & Le, 2016; Le, Hoang, & Le, 2017; Le, & Hoang, 2017) Trong cơng trình này, chúng tơi áp dụng phương pháp toán tử FK phép biến đổi Kustaanheimo-Stiefel (xem (Bartsch, 2003)) để tính phổ lượng nguyên tử hydro đặt từ trường môi trường plasma theo mơ hình chắn SSC Kết thu nghiệm xác số có độ hội tụ đến 32 chữ số thập phân cho trạng thái 1s Cấu trúc báo gồm ba phần chính: phương trình Schrưdinger cho tốn ngun tử hydro chắn, yếu tố ma trận chắn, kết thảo luận, kết luận Phương trình Schrưdinger cho tốn ngun tử hydro chắn Phương trình Schrưdinger cho ngun tử hydro từ trường có xét đến chắn Yukawa hệ nguyên tử viết dạng Hˆ ψ = ε ψ , Hˆ = 1 − ∂2 ∂x + ∂2 ∂y + −i ∂2 ∂z γ x ∂y ∂ ∂ −y +1 ) ∂x γ − (x + y2 x2 y2 z2 Ze−λ , x2 y2 z2 (1) Tạp chí Khoa học Trường ĐHSP TPHCM Tập 17, Số (2020): 10191029 theo hệ thức D = 1/ với λ tham số chắn Do mối liên hệ với bán kính Debye λ (B Saha et al., 2002) nên tham số λ phụ thuộc vào nhiệt độ mật độ plasma Đơn vị độ Tạp chí Khoa học Trường ĐHSP TPHCM Lý Duy Nhất tgk dài bán kính Bohr a = 4π ε 2 / m e2 = 0.5292 Å; đơn vị lượng Hartree 0 e E = 2 / m a2 = 27.21 eV; đơn vị từ trường B = / ea2 = 2.3505 ×105 h e 0 T Trong cơng trình này, chúng tơi tính toán cho nguyên tử hydro nên số hiệu nguyên tử Z = Phương trình (1) giải thành cơng phương pháp tốn tử FK trường hợp khơng chắn, λ = , cơng bố cơng trình (Cao et al., 2019) Bằng cách tương tự, cơng trình chúng tơi giải phương trình (1) cho trường hợp λ theo bước sau: (i) Qua phép biến đổi Kustaanheimo-Stiefel x2= 2u u + 2vv , y = 2u v − 2u v , z = u2 − u2 + v2 − v , 2 2 1 φ = arctan + (v u γ + 2 1 + arctan u1 phương trình (1) viết thành ∂ m ∂ − + − ε − γ ( + 2+ 2+ ∂ 22 2 u v u v ∂v ∂u ∂u ∂ ∂v 1 2 v1 , u2 ) (2) )( u +) v( u v2 + uv + 2 2 ( ) + v ) − Ze u2 v2 u2 v2 ψ (u, v) = 1 2 Ở đây, góc φ :0 →2π đưa vào biến số thêm, chuyển từ không gian (u, v) 2 1 2 không gian (x, y, z) ta cần điều kiện ∂ ψ / ∂ φ = để xác định (ii) Biểu diễn Hamiltonian phương trình (2) qua toán tử sinh, hủy + bsˆ , bˆ (s = 1, 2) sau aˆ =1 s ˆ b =1 s ( αˆ s ( − iβ ˆ s αˆ + i β ˆ s s ), aˆ+ = s ) α 1ˆ = ω u + ( αˆ + + iβˆ ) , + s s + + ˆ , b1 = αˆ + − i β ˆ , s s s ( đó, tốn tử sinh, hủy ∂ αˆ , αˆ + , , (3) ) + β ˆ , β ˆ định nghĩa s s αˆ + = aˆs , aˆ + , ω u − ∂ , Tạp chí Khoa học Trường ĐHSP TPHCM s s ˆ β = ω v + s s ω ∂u s ∂ = ω ∂v s s + , βˆ s s ∂ s ω v − ω ∂u Lý Duy Nhất tgk s ω ∂v s Việc tổ hợp biểu thức (3) giúp cho moment động lượng theo trục Oz: lˆ trở z thành toán tử chéo để dễ dàng áp dụng tính chất bảo tồn đại lượng vật lí Ta thu Hamiltonian dạng đại số Tạp chí Khoa học Trường ĐHSP TPHCM Tập 17, Số (2020): 10191029 ( + + ˆ H = − ω aˆ bˆ + aˆ bˆ + aˆ + bˆ + aˆ + bˆ −aˆ + aˆ − − + + bˆ bˆ − aˆ + aˆ − bˆ bˆ ) 1 2 2 m 1 ˆ − ε − γ aˆ b + + aˆ + bˆ + aˆ + aˆ 2ω 11 1 + γ aˆ bˆ + + aˆ + bˆ + aˆ + aˆ ( 1 1 ( 1 + bˆˆ bˆ + aˆ b 1 1 + bˆ bˆ + aˆ bˆ + )( 1 2 ) + + bˆ + +ˆ+ ˆ + + aˆ + aˆ aˆ b b 2 2 2 + ˆ + aˆ + b + + + aˆ + aˆ bˆ bˆ +1 + ( 16ω3 ) 2 2 1 − ( 1 1 1 2 2 1 1 1 1 2 2 2 2 2ω 2 +2 2 + + + + aˆ bˆ + aˆ + bˆ + aˆ + aˆ +bˆ bˆ + aˆ bˆ + aˆ + bˆ + aˆ + aˆ +bˆ bˆ + − Ze ) + + × aˆ bˆ + aˆ + bˆ + aˆ + aˆ + bˆ bˆ + aˆ bˆ + + + aˆ + bˆ + aˆ + aˆ + bˆ bˆ λ (4) ) Ở đây, biểu thức (4), nhờ vào tính chất bảo toàn moment động lượng theo trục Oz, ta thay tốn tử lˆ trị riêng theo phương trình lˆ ψ (u, v) = mψ (u, v) z z (iii) Sau đó, chúng tơi sử dụng hàm sở hàm dao động tử điều hòa dạng đại số thực tương tác đại số để tính yếu tố ma trận cách dễ dàng (xem chi tiết hàm sở (Cao et al., 2019)) Tuy nhiên, phương trình (4) có thành phần mà tốn tử sinh, hủy dạng mũ việc tính yếu tố ma trận tầm thường Trongλ bước này, tập trung tính yếu tố ma trận − ( aˆ bˆ + aˆ + bˆ+ + aˆ + aˆ +bˆ+ bˆ + aˆ bˆ + aˆ + bˆ+ + aˆ + aˆ +bˆ+ bˆ + ) Vk1k2 , j1 j2 ( ω, λ ) = −Z k1k2m 1 1 1 1 2 2 2 2 j1 j2m (5) e 2ω (iv) Giải phương trình hàm riêng, vector riêng để thu phổ lượng hàm sóng Yếu tố ma trận chắn 3.1 Dạng chuẩn tốn tử dạng e mũ Chúng tơi sử dụng toán tử sinh, hủy + aˆ (ω), aˆ + (ω), bˆ (ω), bˆ (ω) (s = 1, 2) , s s s định nghĩa biểu thức (3), với giao hoán tử scủa chúng cho việc tính tốn yếu tố ma trận phương pháp đại số Trước hết, xét toán tử λ − ( aˆ bˆ + aˆ+ bˆ+ + aˆ+ aˆ +bˆ+ bˆ + aˆ bˆ + aˆ+ bˆ+ + aˆ+ aˆ +bˆ+ bˆ + 2) 1 1 1 1 2 2 2 2 Sˆ = e , (6) 2ω ta thấy tách thành − λ Sˆ = Sˆ với Sˆ ( aˆ bˆ + aˆ + bˆ+ + aˆ + aˆ +bˆ+ b1ˆ +12) − λ ( aˆ bˆ + aˆ + bˆ + + + aˆ + aˆ +bˆ bˆ +1 ) Tạp chí Khoa học Trường ĐHSP TPHCM Sˆ1 = e 1 1 1 Sˆ2 = e 1 2ω 2 2 2 Tập 17, Số (2020): 10191029 2 2ω Điều dễ dàng thấy toán tử sinh, hủy với số toán tử với số giao hốn với Để tính tốn, tốn tử Sˆ , Sˆ đưa dạng chuẩn theo nghĩa tích, tốn tử sinh bên phải, tốn tử hủy bên trái tốn tử trung hịa Do Sˆ S có dạng nên ta trình bày chi tiết bước ˆ chuẩn hóa tốn tử Sˆ Để tinh gọn toán tử hàm e mũ, chúng tơi đặt tốn tử + + + Mˆ = aˆ bˆ , Mˆ = aˆ + bˆ , Nˆ = aˆ + aˆ + bˆ bˆ + , 1 1 1 1 1 (7) Tạp chí Khoa học Trường ĐHSP TPHCM −η M + M qua viết lại toán tử ˆ = e ( + Nˆ S 1) ˆ+ ˆ Lý Duy Nhất tgk với η = λ / 2ω Dễ dàng kiểm chứng toán tử (7) tạo nên đại số kín thể qua hệ thức giao hốn: Mˆ , Mˆ + = , = Nˆ +, Mˆ + Nˆ , Nˆ 2Mˆ , = 2Mˆ Mˆ 1 1 1 1 (8) Sử dụng giao hoán tử (8) thực bước chuẩn hóa tốn tử trình bày chi tiết biểu thức (6.41 ̶ 6.44) trang 232 tài liệu chuyên khảo (Feranchuk et al., 2015), ta thu kết + −η Mˆ + Mˆ + ˆ( η η η Nˆ () Mˆ Mˆ − − + 1 ˆ S1 = η (9) 1+ − ln 1+ ) N1 e = e e 1+ η e Tương tự vậy, ta thu ˆ − η ( Mˆ + ˆ η + ˆM + N ) − + Mˆ 2 − η Mˆ 2 − ln 1+ ) N2 e = Se ˆ2 = e e 1+η 1+ η 3.2 Tính yếu tố ma trận Một cách để giải phương trình Schrưdinger tìm hàm sóng dạng khai triển tuyến tính theo hàm sở cho trước Khi thay giải phương trình vi phân đạo hàm riêng bậc hai ta giải hệ phương trình tuyến tính với hệ số khai triển hàm sóng lượng ẩn số cần tìm Trong này, chúng tơi sử dụng hàm sở hàm sóng nguyên tử hydro đối xứng trụ nguyên tử hydro từ trường đối xứng cầu giữ nguyên đối xứng trụ theo trục dọc theo véc-tơ từ trường Bộ hàm sở qua phép biến đổi Kustaanheimo-Stiefel hàm sóng dao động tử điều hịa bốn chiều viết dạng đại số qua toán tử sinh, hủy xây dựng (Feranchuk et al., 2015) dạng: k,k ,m = Mˆ k1 0, 0, m (10) (ω) (ω) Mˆ k2 Ở đây, ta có ( ) k1 !k2 !k1 m !k2 m ! 0, 0, m (ω) với m < Trạng thái chân không aˆ1 (ω) = 0, ( a (ω) , ˆ2 ) + = (aˆ + bˆ (ω) với m ≥ )m 0, 0, m (ω) + = (aˆ + bˆ (ω) )m (ω) định nghĩa từ phương trình: ˆ (ω) = 0, ˆ (ω) = (11) Tạp chí Khoa học Trường ĐHSP bTPHCM Lý Duy Nhất tgk b2 với điều kiện chuẩn hóa 0 (ω) = Dùng biểu thức quen thuộc toán tử sinh, hủy tác dụng lên véc-tơ trạng thái, có sách giáo khoa lượng tử, ta thu biểu thức tác dụng sau: Mˆ k1 , k2 , m k1 −1, k2 , , k k1 m m = +1, k , m Mˆ + k , k , m k1 1k1 m 1 k , = 1 (12) Nˆ k1 , k2 , = ( 2k1 + m + 1) k1 , k2 , m m Ta có biểu thức tương tự cho tác dụng toán tử Mˆ , Mˆ , Nˆ + 2 Tạp chí Khoa học Trường ĐHSP TPHCM Tập 17, Số (2020): 10191029 Ta tính yếu tố ma trận Vk k , j j cách khai triển thành phần dạng mũ 2 toán tử Sˆ S theo chuỗi Taylor Tiếp theo, tương tác đại số toán tử ˆ + Mˆ , Mˆ s s , Nˆ ( s = 1, với 2) lên hàm sở j1 j2m theo biểu thức (12) Các bước tính tốn chi tiết tương tự trình bày phụ lục (Cao et al., 2019) Kết thu ∞ k− k−j jˆ η (i) S j j m = ( η ) k, j m (13) A , −1 ( ) ∑ k =0 Aklm ( η) = (ii) j1 + m ( 1+ η ) +k +1 l ∑ p=max( 0,l −k ) ˆ ˆ SS jjm = ∑2 (1−11)2 ∞ k− η k =0 k − j1 A ( 1+ η ) +k +1 η k − j +l −j ) k ,l =0 l l m k k m p p k l p k l p η p j = ∞ ( −∑ kj1m j1 + m kj1m (15) ( η) k − j1 +l − k, j m j2 1 j k j l 2 m 2 (η ˆ ) S Akj m (iii) (14) ( η) Alj k, l, m m Vk k ,j j (η) = −Z k1k2m j1 j2m Sˆ Sˆ 2 = −Z k k m ( −1) 2 ∑ ∞ k − j1 +l − j2 η k − j +l − j2 1 j k j l 2 m 2 k ,l =0 = −Z )k jm 1 1 2 (η) kj1m 2 (η Trong trường hợp riêng, không chắn, Vk k , j j ( 0) = −Zδ k , j δ k ,j 10 ( η) m k, l, (16) lj2m k1 − j1 + k2 − j2 ( −η ) A A 1 j k j k 2 m 2 k j m ( η ) A A λ = , ta có kết Tạp chí Khoa học Trường ĐHSP TPHCM Khi 2 khơng có từ trường, 1 Tập 17, Số (2020): 10191029 toán túy nguyên tử 2 hydro với nghiệm giải tích Các hàm sóng Rn ( r ) Ylm ( ϕ,θ ) E = −1 / 2n2 , ( n = 1, 2, ) ứng với mức lượng n Các yếu tố ma trận cịn lại hồn tồn tốn ngun tử hydro từ trường khơng chắn đưa Phụ lục (Cao et al., 2019) Ví dụ với yếu tố ma trận ω + + K = j j m Nˆ + Nˆ + − Mˆ − Mˆ − Mˆ − Mˆ kk m , ( j1 j2 ,k1k2 ) ta dễ dàng tính sau: 11 2 Tạp chí Khoa học Trường ĐHSP TPHCM K + j k1k2 , j1 j2 ( ω) = ω Lý Duy Nhất tgk 2( j + m +1) δ δ − j1 j1 m 4 − k1 , j1 2 δ k1 , j1 −1 k2 , j2 k2 , j2 δ j1 1 j1 m 1 j2 j2 m δ k , j −1δ k , j δk , j +1δk , j − j2 1− j2 m 1 2 1 1 δ k , j +1δ k , (17) j Kết thảo luận Với yếu tố ma trận tính theo biểu thức (16), (17), chương trình giải số phương pháp tốn tử FK viết ngơn ngữ FORTRAN có sử dụng thư mục “subroutine DSYGVX” gói LAPACK (Netlib.org LAPACK: Linear Algebra PACKage, n.d.) với biến số thực viết lại REAL*16 tương đương độ xác 32 chữ số thập phân Chương trình cung cấp miễn phí theo yêu cầu người dùng công bố tạp chí tương ứng Một điểm bật phương pháp tốn tử FK có đưa vào tham số tự ω qua biểu thức (3) định nghĩa toán tử sinh, hủy Chọn giá trị thích hợp cho tham số ta thu độ xác cao tính lượng Điều tương đương với tốc độ hội tụ cao tính lượng với độ xác cho trước Tuy tính chất tham số tự quan sát nhiều hệ khác nhau, xem ví dụ (Feranchuk et al., 2015), với trường hợp ta cần kiểm chứng lại số Với tốn xét ta tìm miền giá trị tham số cho tốc độ hội tụ cao Trên Hình 1, thể cho ví dụ với tham số chắn λ = (Hình 1a) λ = (Hình 1b) cho trường 0.01 0.1 hợp khơng có từ trường lượng trạng thái 1s tính với độ xác đến 15 chữ số thập phân (đường màu đỏ) 10 chữ số thập phân (đường màu đen) Trục tung biểu diễn số bậc gần cần phải tính để đạt độ xác cho ứng với giá trị tham số tự khác nhau, biểu diễn trục hồnh Hình cho thấy vùng tham số tự tối ưu tồn gần với độ xác 10 hay 15 chữ số thập phân Thử cho nhiều độ xác khác ta kết luận miền tham số tối ưu khơng phụ thuộc vào độ xác cần tính Điều cho phép ta tính lượng với độ xác thấp, ví dụ đến 10 chữ số thập phân, với tham số tự khác qua xác định miền tham số tối ưu Sau sử dụng tham số để tính lượng với độ xác lên đến 32 chữ số thập phân Trong chương trình tính tốn, giá trị tham số tối ưu tính trước cho số trường hợp ứng với cường độ từ trường khác từ γ = →100 với bước nhảy ∆γ = 0.1×(γ +1) tham số chắn λ = →1 với bước nhảy ∆λ = 0.05 Các trường hợp khác thu nội suy Như vậy, chương trình tính tốn tích hợp giá trị tối ưu tham số tự 12 Tạp chí Khoa học Trường ĐHSP TPHCM Tập 17, Số (2020): 10191029 Hình Độ hội tụ thể qua bậc gần s theo tham số tự ω cho nghiệm tính với độ xác đến 15 chữ số thập phân (đường màu đỏ) 10 chữ số thập phân (đường màu đen) γ = hai trường hợp tham số chắn: λ= (b) cho (a) 0.01 λ = 0.1 Chương trình kiểm chứng cho giá trị lượng 1s xác đến 32 chữ số thập phân với tham số chắn λ có giá trị từ đến 1.0 a.u ứng với từ trường γ có giá trị từ đến 100 Để minh họa, Bảng đưa số kết cho trường hợp khơng có từ trường, kết so sánh với cơng trình khác (Paul, & Ho, 2009; Saha et al., 2002; Soylu, 2012) Bảng cho trường hợp từ trường γ = Trong trường hợp λ > 0.02 bậc gần hàm sóng phải chọn lên đến 180, tương đương với (s) số lượng phần tử ma trận lên tới 8281× 8281 hàm sở đối xứng Chương trình phát triển cho trạng thái kích thích cho dạng chắn khác GECSC công bố cơng trình Bảng Các giá trị lượng mức 1s theo tham số chắn λ hệ đơn vị nguyên tử trường hợp từ trường γ = 0.000 0.010 -0.5000000000000000000000000000000 -0.4900745067466941718092082523330 -0.490075 -0.490075 Tài liệu (Saha et al., 2002) -0.50000 -0.49008 0.050 -0.4518164285245055438562170757554 -0.451816 -0.451816 -0.45182 0.100 -0.4070580306134031567545070703611 -0.407058 -0.405705 -0.40705 1.000 -0.0102857899900176968047742153111 λ Tài liệu (Soylu, 2012) ε 13 Tài liệu (Paul, & Ho, 2009) 0.16598 Tạp chí Khoa học Trường ĐHSP TPHCM Lý Duy Nhất tgk Bảng Các giá trị lượng mức 1s theo tham số chắn λ hệ đơn vị nguyên tử trường hợp từ trường γ = λ ε 0.000 0.010 0.050 0.100 1.000 -0.3311688967331580356102982 -0.3212281942984668084487882 -0.2826224109660852632933711 -0.2368473261576955381981450 0.2593299390848640488750535 Kết luận Qua phép biến đổi Kustaanheimo-Stiefel, chúng tơi áp dụng thành cơng phương pháp tốn tử FK cho tốn ngun tử hydro mơi trường plasma với chắn Yukawa Chúng giải phương trình Schrưdinger cho ngun tử hydro từ trường có xét chắn thu lượng với độ xác đến 32 chữ số thập phân cho trạng thái 1s so sánh với kết tính số trước cho thấy vượt trội độ xác Một chương trình FORTRAN viết để tự động hóa hồn tồn tính tốn cho người dùng Chương trình phát triển kiểm chứng tiếp cho trạng thái kích thích cao, độ thay đổi từ trường rộng miền tham số chắn mở rộng Ngồi ra, phát triển cho toán heli nguyên tử cao Tuyên bố quyền lợi: Các tác giả xác nhận hồn tồn khơng có xung đột quyền lợi Lời cảm ơn: Bài báo nghiên cứu tài trợ Trường Đại học Sư phạm Thành phố Hồ Chí Minh qua Đề tài trọng điểm cấp sở mã số CS.2019.19.43TĐ Các tác giả chân thành cám ơn GS TSKH Lê Văn Hoàng hướng dẫn, đặt vấn đề góp ý, chỉnh sửa cho báo TÀI LIỆU THAM KHẢO Bahar, M K., & Soylu, A (2014) The hydrogen atom in plasmas with an external electric field Physics of Plasmas, 21(9), 092703 https://doi.org/10.1063/1.4894684 Bahar, M K., & Soylu, A (2015) Confinement effects of magnetic field on two-dimensional hydr-ogen atom in plasmas Physics of Plasmas, 22(5), 052701 https://doi.org/10.1063/1.4919613 Bahar, M K., & Soylu, A (2016) Probe of hydrogen atom in plasmas with magnetic, electric, and Aharonov-Bohm flux fields Physics of Plasmas, 23(9), 092712 https://doi.org/10.1063/1.4963772 Bartsch, T (2003) The Kustaanheimo–Stiefel transformation in geometric algebra Journal of Physics A: Mathematical and General, 36(25), 6963-6978 https://doi.org/10.1088/03054470/36/25/305 Cao, H T X., Ly, D N., & Hoang, D N T (2016) Nang luong trang thai co ban cua nguyen tu hydro tu truong deu co cuong bat ki [Ground state energy of a hydrogen atom in a 14 Tạp chí Khoa học Trường ĐHSP TPHCM Lý Duy Nhất tgk uniform magnetic field with arbitrary strength] Ho Chi Minh City University of Education Journal of Science, 12, 39-51 15 Tạp chí Khoa học Trường ĐHSP TPHCM Tập 17, Số (2020): 10191029 Cao, T X H., Ly, D N., Hoang, N.T D., & Le, V H (2019) High-accuracy numerical calculations of the bound states of a hydrogen atom in a constant magnetic field with arbitrary strength Computer Physics Communications, 240, 138-151 https://doi.org/10.1016/j.cpc.20 19.02.013 Feranchuk, I., Ivanov, A., Le, V.-H., & Ulyanenkov, A (2015) Non-perturbative Description of Quantum Systems (Vol 894) Cham: Springer International Publishing https://doi.org/10.1007/978-3-319-13006-4 Hoang, N T D., Nguyen, D A P., Hoang, V H., & Le, V H (2016) Highly accurate analytical energy of a two-dimensional exciton in a constant magnetic field Physica B: Condensed Matter, 495, 16-20 https://doi.org/10.1016/j.physb.2016.04.038 Kar, S., & Ho, Y K (2005) Electron affinity of the hydrogen atom and a resonance state of the hydrogen negative ion embedded in Debye plasmas New Journal of Physics, 7, 141-141 https://doi.org/10.1088/1367-2630/7/1/141 Le, D N., Hoang, N.T D., & Le, V.H (2017) Exact analytical solutions of a two-dimensional hydrogen atom in a constant magnetic field Journal of Mathematical Physics, 58(4), 042102 https://doi.org/10.1063/1.4979618 Le, D N., & Hoang, D N T (2017) Binding Energy of Exciton in Monolayer Semiconductor Ws With Yukawa-Like Screening Potential Ho Chi Minh City University of Education Journal of Science, 14(9), 43-50 Retrieved from http://journal.hcmue.edu.vn/index.php/hcmuejos/issue/archive Lin, C Y., & Ho, Y K (2011) The photoionization of excited hydrogen atom in plasmas Computer Physics Communications, 182(1), 125-129 https://doi.org/10.1016/j.cpc.20 10.06.013 Main, J., Schwacke, M., & Wunner, G (1998) Hydrogen atom in combined electric and magnetic fields with arbitrary mutual orientations Physical Review A, 57(2), 1149-1157 https://doi org/10.1103/PhysRevA.57.1149 Netlib.org LAPACK: Linear Algebra PACKage (n.d.) Subroutine DSYGVX.f Retrieved from http://www.netlib.org/lapack/explore-ht-ml/d2/d97/dsyevx-8f.html Paul, S., & Ho, Y K (2009) Hydrogen atoms in Debye plasma environments Physics of Plasmas, 16(6), 063302 https://doi.org/10.1063/1.3152602 Paul, S., & Ho, Y K (2010) Two-colour three-photon transitions in a hydrogen atom embedded in Debye plasmas Journal of Physics B: Atomic, Molecular and Optical Physics, 43(6), 065701 https://doi.org/10.1088/0953-4075/43/6/065701 Paul, S., & Ho, Y K (2011) Solution of the generalized exponential cosine screened Coulomb potential Computer Physics Communications, 182(1), 130-133 https://doi.org/10.1016/j cpc.2010.06.014 Saha, B., Mukherjee, P K., & Diercksen, G H F (2002) Energy levels and structural properties of compressed hydrogen atom under Debye screening Astronomy & Astrophysics, 396(1), 337–344 https://doi.org/10.1051/0004-6361:20021350 Saha, J K., Mukherjee, T K., Mukherjee, P K., & Fricke, B (2011) Hyperpolarizability of hydrogen atom under spherically confined Debye plasma The European Physical Journal D, 62(2), 205-211 https://doi.org/10.1140/epjd/e2011-10668-4 16 Tạp chí Khoa học Trường ĐHSP TPHCM Lý Duy Nhất tgk Sasmal, G P (2014) On computation for a hydrogen atom in arbitrary magnetic fields using finite volume method Journal of Atomic and Molecular Sciences, 5(3), 187-205 https://doi.org/10 4208/jams.110813.021414a Shukla, P K., & Eliasson, B (2008) Screening and wake potentials of a test charge in quantum plasmas Physics Letters A, 372(16), 2897-2899 https://doi.org/10.1016/j.physleta.2007.12 067 Soylu, A (2012) Plasma screening effects on the energies of hydrogen atom Physics of Plasmas, 19(7), 072701 https://doi.org/10.1063/1.4736947 Thirumalai, A., & Heyl, J S (2009) Hydrogen and helium atoms in strong magnetic fields Physical Review A, 79(1), 012514 https://doi.org/10.1103/PhysRevA.79.012514 HIGHLY ACCURATE NUMERICAL ENERGY FOR THE GROUND STATE OF A HYDROGEN ATOM IN A PLASMA WITH THE PRESENCE OF A UNIFORM MAGNETIC FIELD Ly Duy Nhat1,2*, Hoang Do Ngoc Tram2 Ho Chi Minh City University of Education, Vietnam Vietnam Atomic Energy Institute, Vietnam * Corresponding author: Ly Duy Nhat – Email: nhatld@hcmue.edu.vn Received: March 31, 2020; Revised: May 21, 2020; Accepted: June 10, 2020 ABSTRACT The FK operator method has been recently applied successfully to the problem of a hydrogen atom in a magnetic field with any intensity In this work, this method is applied more generally, considering the screening effect Screening effect plays an important role in practical astronomy research on hydrogen atom in a plasma environment In this paper, the FK operator method is combined with the Kustaanheimo-Stiefel transformation, so the wavefunction is expanded in terms of the eigenfunctions of a four-dimensional harmonic oscillator There is an improvement in algebraic calculations to obtain the matrix elements of the screening term in the Hamiltonian, which results in the explicit expressions of the matrix elements in the series expansion These results help to build a FORTRAN program finding highly accurate numerical energies with the precision up to 32 decimal places For an illustration purpose, energies for the ground state with the screening parameter ranging from to a.u are shown in the paper Keywords: algebraic method; hydrogen atom; screened Coulomb potential; annihilation and creation operators; basic set 17 ... dụng thành cơng phương pháp tốn tử FK cho tốn ngun tử hydro mơi trường plasma với chắn Yukawa Chúng giải phương trình Schrưdinger cho ngun tử hydro từ trường có xét chắn thu lượng với độ xác đến... Trong cơng trình này, chúng tơi tính toán cho nguyên tử hydro nên số hiệu nguyên tử Z = Phương trình (1) giải thành cơng phương pháp tốn tử FK trường hợp khơng chắn, λ = , cơng bố cơng trình (Cao. .. tơi sử dụng hàm sở hàm sóng nguyên tử hydro đối xứng trụ nguyên tử hydro từ trường đối xứng cầu giữ nguyên đối xứng trụ theo trục dọc theo véc-tơ từ trường Bộ hàm sở qua phép biến đổi Kustaanheimo-Stiefel