Bối cảnh và vấn đề nghiên cứu

Trong lịch sử phát triển của khoa học tự nhiên, toán học và vật lý luôn có mối liên hệ hữu cơ mật thiết. Sự hình thành và phát triển của phương trình toán lý từ thế kỷ XVIII đã cung cấp công cụ toán học để mô tả các quy luật biến thiên của các đại lượng vật lý trong tự nhiên. Khi giải các phương trình vi phân đạo hàm riêng mô tả các trường vật lý (như phương trình Laplace, Poisson, Helmholtz, phương trình sóng và phương trình Schrödinger) trong các hệ tọa độ cong (tọa độ cầu, tọa độ trụ) bằng phương pháp tách biến, nghiệm nhận được thường quy về các phương trình vi phân thường mà nghiệm giải tích của chúng là các hàm đặc biệt như đa thức Legendre, hàm cầu và hàm Bessel.

Trong chương trình đào tạo đại học ngành Sư phạm Vật lý và Vật lý học, sinh viên phải tiếp cận nhiều học phần chuyên ngành đòi hỏi khối lượng tính toán giải tích lớn như Cơ lý thuyết, Điện động lực học, Cơ học lượng tử, Vật lý thống kê, Vật lý nguyên tử và hạt nhân. Mặc dù sự phát triển của công nghệ tính toán hiện nay cho phép giải số trực tiếp các bài toán phức tạp, việc nắm vững bản chất giải tích và tính chất toán học của các hàm đặc biệt vẫn giữ vai trò then chốt trong việc tiên đoán nghiệm gần đúng, phân tích tính đối xứng và hiểu rõ bản chất vật lý của hệ. Xuất phát từ nhu cầu thực tiễn đó, khóa luận tốt nghiệp của tác giả Nguyễn Hữu Phát (ngành Sư phạm Vật lý, Trường Đại học Sư phạm Thành phố Hồ Chí Minh, năm 2022) dưới sự hướng dẫn của thầy Nguyễn Lê Anh đã thực hiện đề tài nhằm hệ thống hóa cơ sở toán học và ứng dụng vật lý của các hàm đặc biệt.

Mục tiêu và nhiệm vụ nghiên cứu của đề tài được xác định cụ thể:

  1. Xây dựng một tài liệu tham khảo có hệ thống về lý thuyết các hàm và đa thức đặc biệt phục vụ học tập và nghiên cứu cho sinh viên chuyên ngành Vật lý.
  2. Khảo sát chi tiết cấu trúc toán học, tính chất giải tích, hệ thức truy hồi, tính trực giao và nghiệm chuỗi của đa thức Legendre, hàm Legendre liên kết và các hàm Bessel (loại một, loại hai, hàm Hankel, hàm cầu Bessel).
  3. Hệ thống hóa các bài toán vật lý điển hình ứng dụng các hàm đặc biệt trong tĩnh điện học, điện động lực học, quang học sóng và cơ học lượng tử.
  4. Xây dựng các thuật toán tính toán số và chương trình hỗ trợ tính toán, vẽ đồ thị phân tích tính chất của các hàm đặc biệt.

Đối tượng nghiên cứu của khóa luận là các hàm đặc biệt thường gặp trong toán lý, trọng tâm gồm đa thức Legendre $P_n(x)$, hàm Legendre liên kết $P_n^m(x)$, hàm Bessel loại một $J_n(x)$, hàm Neumann (Bessel loại hai) $N_n(x)$, hàm Hankel và hàm cầu Bessel $j_n(x)$. Phạm vi nghiên cứu giới hạn trong không gian tọa độ thực và phức liên quan đến các bài toán biên vật lý đối xứng trục, đối xứng trụ và đối xứng cầu trong các điều kiện biên Dirichlet và Neumann tiêu chuẩn.


Cơ sở lý thuyết và phương pháp

Khóa luận dựa trên nền tảng của lý thuyết phương trình toán lý và lý thuyết Sturm-Liouville về bài toán biên tự liên hợp. Khung lý thuyết trọng tâm xoay quanh phương pháp tách biến để giải các phương trình vi phân đạo hàm riêng cấp hai trong hệ tọa độ cực, tọa độ trụ $(r, \phi, z)$ và tọa độ cầu $(r, \theta, \phi)$:

  • Phương trình Laplace và Poisson: $\nabla^2 \phi = 0$ và $\nabla^2 \phi = -\rho/\varepsilon_0$ trong tĩnh điện học.
  • Phương trình Helmholtz: $(\nabla^2 + k^2)\psi = 0$ trong bài toán truyền sóng và tán xạ lượng tử.
  • Phương trình vi phân Legendre: $\frac{d}{dx}\left[(1-x^2)\frac{dP_n(x)}{dx}\right] + n(n+1)P_n(x) = 0$.
  • Phương trình vi phân Legendre liên kết: $(1-x^2)\frac{d^2 P_n^m(x)}{dx^2} - 2x\frac{dP_n^m(x)}{dx} + \left[n(n+1) - \frac{m^2}{1-x^2}\right]P_n^m(x) = 0$.
  • Phương trình vi phân Bessel: $x^2 \frac{d^2 y}{dx^2} + x \frac{dy}{dx} + (x^2 - n^2)y = 0$.

Phương pháp nghiên cứu kết hợp giữa giải tích thuần túy và tính toán số học thuật:

  • Phương pháp giải tích: Sử dụng hàm sinh (generating function) để dẫn xuất các dạng tường minh, hệ thức truy hồi, tính chẵn lẻ và biểu diễn vi phân (công thức Rodrigues); áp dụng tích phân trực giao và khai triển chuỗi trực chuẩn (chuỗi Fourier-Legendre, Fourier-Bessel).
  • Phương pháp bài toán biên: Thiết lập nghiệm tổng quát dưới dạng chuỗi các hàm riêng, áp dụng điều kiện biên vật lý (tại gốc tọa độ, bề mặt vật dẫn, biên vô cực) để xác định các hệ số khai triển đơn trị.
  • Phương pháp tính toán số: Lập trình tính toán các giá trị hàm, xử lý các công thức truy hồi tránh sai số làm tròn, xác định các điểm không (zeros) của hàm Bessel và đạo hàm cấp một, vẽ đồ thị đường biên biểu diễn sự phân bố của các hàm.
Nhóm hàm đặc biệt Phương trình vi phân đặc trưng Khái niệm / Công cụ giải tích cơ bản Miền xác định & Tính trực giao
Đa thức Legendre $P_n(x)$ $\frac{d}{dx}[(1-x^2)y'] + n(n+1)y = 0$ Hàm sinh $g(t,x) = (1-2xt+t^2)^{-1/2}$, công thức Rodrigues $x \in [-1, 1]$, $\int_{-1}^1 P_n P_m dx = \frac{2\delta_{mn}}{2n+1}$
Hàm Legendre liên kết $P_n^m(x)$ $(1-x^2)y'' - 2xy' + [n(n+1) - \frac{m^2}{1-x^2}]y = 0$ Đạo hàm cấp $m$ của $P_n(x)$: $(1-x^2)^{m/2}\frac{d^m P_n}{dx^m}$ $
Hàm Bessel loại một $J_n(x)$ $x^2 y'' + xy' + (x^2 - n^2)y = 0$ Hàm sinh chuỗi lũy thừa, biểu diễn tích phân $x \in [0, \infty)$, trực giao theo hàm trọng số $x$
Hàm Neumann (Bessel loại 2) $N_n(x)$ $x^2 y'' + xy' + (x^2 - n^2)y = 0$ Nghiệm độc lập tuyến tính thứ hai, công thức Wronskian Phân kỳ tại gốc $x \to 0$
Hàm cầu Bessel $j_n(x), y_n(x)$ $x^2 y'' + 2xy' + [x^2 - n(n+1)]y = 0$ Biểu diễn qua hàm Bessel bậc bán nguyên: $\sqrt{\frac{\pi}{2x}}J_{n+1/2}(x)$ Ứng dụng trong giếng thế cầu và tán xạ

Nguồn dữ liệu của đề tài được trích xuất từ các tài liệu kinh điển về phương pháp toán lý, các bài toán chuẩn trong giáo trình cơ học lượng tử, tĩnh điện học và các mô hình mô phỏng dao động cơ học - sóng điện từ.


Nội dung chính theo từng chương

Chương 1: Đa thức Legendre và hàm cầu

Chương 1 trình bày chi tiết cơ sở toán học và ứng dụng vật lý của đa thức Legendre và hàm Legendre liên kết:

  • Cơ sở vật lý và hàm sinh trong tĩnh điện: Đa thức Legendre $P_n(\cos\theta)$ được xây dựng trực tiếp từ bài toán thế tĩnh điện của một điện tích điểm $q$ đặt tại $z = a$ trên trục $z$. Thế tĩnh điện tại điểm $(r, \theta)$ với $r > a$ được khai triển theo định lý cosin thành: $$\phi = \frac{q}{4\pi\varepsilon_0 r}\sum_{n=0}^{\infty} \left(\frac{a}{r}\right)^n P_n(\cos\theta)$$ Từ đó xác định hàm sinh của đa thức Legendre: $g(t, x) = (1 - 2xt + t^2)^{-1/2} = \sum_{n=0}^{\infty} P_n(x)t^n$ với $|t| < 1, x = \cos\theta$.
  • Các tính chất đại số và giải tích:
    • Giá trị đặc biệt: $P_n(1) = 1$, $P_n(-1) = (-1)^n$, $P_{2n+1}(0) = 0$, $P_{2n}(0) = (-1)^n \frac{(2n-1)!!}{(2n)!!}$.
    • Tính chẵn lẻ: $P_n(-x) = (-1)^n P_n(x)$, liên hệ trực tiếp với tính bảo toàn tính chẵn lẻ của toán tử Hamiltonian dưới phép nghịch đảo không gian $r \to -r$ và bảo toàn moment động lượng quỹ đạo.
    • Công thức số hạng tổng quát: $P_n(x) = \sum_{k=0}^{[n/2]} (-1)^k \frac{(2n-2k)!}{2^n k!(n-k)!(n-2k)!} x^{n-2k}$.
    • Công thức Rodrigues: $P_n(x) = \frac{1}{2^n n!} \frac{d^n}{dx^n}(x^2 - 1)^n$, công cụ hữu hiệu để tính toán giải tích và chứng minh tính trực giao.
  • Khai triển đa cực tĩnh điện: Tác giả phân tích khai triển thế của phân bố điện tích $\rho(r_2)$ qua hệ tọa độ vector, dẫn xuất các số hạng đơn cực ($P_0$), lưỡng cực ($P_1$), tứ cực ($P_2$) và bát cực ($P_3$). Số hạng tứ cực điện $\phi_2(r_1)$ được biểu diễn dưới dạng tích phân chứa tensor đối xứng không vết: $\frac{1}{4\pi\varepsilon_0} \frac{1}{2r_1^5}\sum_{i,j} x_{1i}x_{1j}\int [3x_{2i}x_{2j} - \delta_{ij}r_2^2]\rho(r_2)d\tau_2$.
  • Hệ thức truy hồi và phương trình vi phân:
    • Hệ thức truy hồi ba số hạng cơ bản: $(2n+1)xP_n(x) = (n+1)P_{n+1}(x) + nP_{n-1}(x)$.
    • Thuật toán truy hồi ổn định số học cho máy tính nhằm triệt tiêu sai số làm tròn: $P_{n+1}(x) = 2xP_n(x) - P_{n-1}(x) - \frac{xP_n(x) - P_{n-1}(x)}{n+1}$.
    • Dẫn xuất các hệ thức vi phân $P_{n+1}'(x) - P_{n-1}'(x) = (2n+1)P_n(x)$ và chứng minh $P_n(x)$ là nghiệm của phương trình vi phân Legendre cấp hai tự liên hợp: $\frac{d}{dx}[(1-x^2)P_n'(x)] + n(n+1)P_n(x) = 0$.
  • Tính trực giao và bài toán biên:
    • Hệ thức trực chuẩn: $\int_{-1}^1 P_n(x)P_m(x)dx = \frac{2\delta_{mn}}{2n+1}$, cho phép khai triển một hàm $f(x)$ bất kỳ thành chuỗi Legendre $f(x) = \sum a_n P_n(x)$ với $a_n = \frac{2n+1}{2}\int_{-1}^1 f(x)P_n(x)dx$.
    • Khai triển hàm delta Dirac: $\delta(1\mp x) = \sum_{n=0}^\infty \frac{2n+1}{2}(\pm 1)^n P_n(x)$.
    • Bài toán quả cầu dẫn điện trong điện trường đều $E_0$: Bằng cách giải phương trình Laplace $\nabla^2 V = 0$ với điều kiện biên Dirichlet, xác định thế tĩnh điện bên ngoài quả cầu bán kính $r_0$ là $V(r, \theta) = -E_0 r P_1(\cos\theta)\left(1 - \frac{r_0^3}{r^3}\right)$, mật độ điện tích mặt cảm ứng $\sigma = 3\varepsilon_0 E_0 \cos\theta$, và moment lưỡng cực cảm ứng $P = 4\pi\varepsilon_0 r_0^3 E_0$.
    • Bài toán vòng dây dẫn tích điện $q$ bán kính $a$: Xác định thế tĩnh điện tại mọi điểm $r > a$ thông qua so sánh điều kiện biên trên trục đối xứng $z$.
  • Hàm Legendre liên kết $P_n^m(x)$: Xuất hiện khi tách biến phương trình Laplace trong hệ tọa độ cầu khi có sự phụ thuộc vào góc phương vị $\phi$. Định nghĩa qua đạo hàm $P_n^m(x) = (1-x^2)^{m/2}\frac{d^m P_n(x)}{dx^m}$ với điều kiện $|m| \le n$. Tác giả chỉ ra mối liên hệ vật lý lượng tử: điều kiện $m \le n$ tương ứng với hệ thức giá trị kỳ vọng $\langle L_z^2 \rangle \le \langle L^2 \rangle$, trong đó $m$ là trị riêng của $L_z$ và $n(n+1)$ là trị riêng của toán tử bình phương moment động lượng toàn phần $L^2$.

Chương 2: Hàm Bessel

Chương 2 tập trung khảo sát hệ thống hàm Bessel và các ứng dụng trong bài toán đối xứng trụ và đối xứng cầu:

  • Hàm Bessel loại một $J_n(x)$: Định nghĩa qua hàm sinh bậc nguyên $g(x, t) = \exp\left[\frac{x}{2}\left(t - \frac{1}{t}\right)\right] = \sum_{n=-\infty}^\infty J_n(x)t^n$, dạng biểu diễn chuỗi lũy thừa, các hệ thức truy hồi vi phân, biểu diễn tích phân Bessel-Schläfli và tính trực giao trên đoạn hữu hạn theo các điểm không (zeros) của hàm. Tác giả lập bảng số liệu các điểm không của hàm $J_n(x)$ và đạo hàm cấp một $J_n'(x)$.
  • Hàm Neumann (Bessel loại hai) $N_n(x)$ và hàm Hankel: Khảo sát nghiệm độc lập tuyến tính thứ hai của phương trình vi phân Bessel đối với bậc nguyên $n$, công thức Wronskian $W[J_n(x), N_n(x)] = \frac{2}{\pi x}$, và các tổ hợp tuyến tính phức biểu diễn sóng chạy (hàm Hankel loại một $H_n^{(1)}(x)$ và loại hai $H_n^{(2)}(x)$).
  • Khai triển tiệm cận: Phân tích trạng thái tiệm cận của hàm Bessel khi đối số $x \gg 1$ và $x \ll 1$, đánh giá nghiệm gần đúng cho $J_0(x)$ và các đường biên dao động.
  • Hàm cầu Bessel ($j_n(x), y_n(x)$): Khảo sát phương trình vi phân cầu Bessel xuất phát từ phương trình sóng dạng cầu, các hệ thức truy hồi và giá trị giới hạn tại gốc tọa độ.
  • Ứng dụng vật lý điển hình của hàm Bessel:
    • Quang học sóng: Hiện tượng nhiễu xạ Fraunhofer qua khẩu độ tròn dẫn đến hình ảnh vân nhiễu xạ đĩa Airy được mô tả bởi $\left[\frac{2J_1(x)}{x}\right]^2$.
    • Cơ học dao động: Mô hình dao động màng tròn/mặt trống với điều kiện biên cố định tại mép màng $r = a$, xác định các mode dao động riêng tương ứng với các điểm không của $J_m(k_{mn} a) = 0$.
    • Điện từ trường: Phân tích các mode sóng TE và TM trong khoang cộng hưởng trụ kim loại.
    • Cơ học lượng tử và tán xạ: Giải bài toán giếng thế vuông đối xứng cầu, gần đúng WKB, và khai triển sóng phẳng Rayleigh $e^{ikr\cos\gamma} = \sum_{n=0}^\infty i^n (2n+1) j_n(kr) P_n(\cos\gamma)$. Thiết lập công thức tiết diện tán xạ toàn phần trong lý thuyết tán xạ lượng tử: $$\sigma_{tot} = \frac{4\pi}{k^2}\sum_{l=0}^\infty (2l+1)\sin^2\delta_l$$ với $\delta_l$ là độ lệch pha của sóng riêng phần thứ $l$.

Kết quả và đóng góp

Khóa luận đã tổng hợp và đạt được các kết quả cụ thể CÓ TRONG văn bản:

  1. Hệ thống hóa toán học hoàn chỉnh: Thiết lập đầy đủ quy trình giải tích từ hàm sinh, công thức vi phân, hệ thức truy hồi, công thức Rodrigues đến hệ thức trực giao cho cả hai nhóm hàm trọng tâm (Legendre và Bessel).
  2. Giải quyết chính xác các bài toán biên vật lý mẫu mực:
    • Xác định phân bố thế và trường tĩnh điện của hệ đa cực (đơn cực, lưỡng cực, tứ cực tuyến tính, bát cực tuyến tính).
    • Thiết lập nghiệm giải tích cho bài toán quả cầu dẫn điện trung hòa bán kính $r_0$ trong điện trường đều: thế $V(r, \theta) = -E_0 r P_1(\cos\theta)(1 - r_0^3/r^3)$, mật độ điện mặt $\sigma = 3\varepsilon_0 E_0 \cos\theta$, moment lưỡng cực cảm ứng $P = 4\pi\varepsilon_0 r_0^3 E_0$.
    • Chứng minh công thức độ lệch pha và tiết diện tán xạ toàn phần trong cơ học lượng tử qua đa thức Legendre và hàm cầu Bessel.
  3. Đóng góp về phương pháp tính toán số: Cung cấp thuật toán truy hồi cải tiến $P_{n+1}(x) = 2xP_n(x) - P_{n-1}(x) - \frac{xP_n(x) - P_{n-1}(x)}{n+1}$ nhằm triệt tiêu sự tích tụ sai số làm tròn khi tính toán số bậc cao trên máy tính; lập bảng giá trị các điểm không của hàm Bessel và đạo hàm cấp một.
  4. Giá trị học liệu sư phạm: Đóng góp của tác giả so với các tài liệu đơn lẻ trước đây là đã tích hợp lý thuyết toán lý thuần túy với các bài toán vật lý cụ thể, tạo thành tài liệu tham khảo đồng bộ cho sinh viên chuyên ngành Sư phạm Vật lý.
Bài toán vật lý Hàm đặc biệt sử dụng Phương trình vi phân / Công thức cốt lõi Kết quả giải tích chính
Thế quả cầu dẫn trong điện trường đều Đa thức Legendre $P_1(\cos\theta)$ $\nabla^2 V = 0$ (Laplace), ĐK biên Dirichlet $V = -E_0 r P_1(\cos\theta)(1 - r_0^3/r^3)$, $\sigma = 3\varepsilon_0 E_0 \cos\theta$
Khai triển thế đa cực điện Đa thức Legendre $P_n(\cos\theta)$ $g(t, x) = (1-2xt+t^2)^{-1/2} = \sum P_n(x)t^n$ Tách thế đơn cực, lưỡng cực, tứ cực (tensor không vết), bát cực
Khai triển sóng phẳng Rayleigh Đa thức Legendre & Hàm cầu Bessel $e^{ikr\cos\gamma} = \sum a_n j_n(kr)P_n(\cos\gamma)$ Xác định hệ số $a_n = i^n(2n+1)$
Tán xạ lượng tử Đa thức Legendre $P_l(\cos\theta)$ $f(\theta) = \frac{1}{k}\sum (2l+1)e^{i\delta_l}\sin\delta_l P_l(\cos\theta)$ Tiết diện tán xạ toàn phần $\sigma_{tot} = \frac{4\pi}{k^2}\sum (2l+1)\sin^2\delta_l$
Màng rung tròn / Đĩa Airy Hàm Bessel loại một $J_0(x), J_1(x)$ $x^2 y'' + xy' + (x^2 - n^2)y = 0$ Cường độ đĩa Airy $I \propto [2J_1(x)/x]^2$, mode màng $J_m(k_{mn}a) = 0$

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế:
    • Khóa luận tập trung chủ yếu vào hai họ hàm đặc biệt phổ biến nhất là họ Legendre (đa thức Legendre, hàm liên kết, hàm cầu) và họ Bessel (loại một, loại hai, Hankel, cầu Bessel).
    • Chưa mở rộng khảo sát các họ hàm đặc biệt quan trọng khác trong vật lý lượng tử và toán lý như đa thức Hermite (trong bài toán dao động tử điều hòa lượng tử), đa thức Laguerre (trong bài toán nguyên tử hydro), hàm siêu bội (hypergeometric functions) hay hàm Airy.
    • Các chương trình tính toán số dừng lại ở mức độ tính toán bảng số liệu điểm không và hỗ trợ vẽ đồ thị đường biên, chưa xây dựng thành một phần mềm mô phỏng tương tác hoàn chỉnh.
  • Hướng nghiên cứu tiếp:
    • Tiếp tục mở rộng nghiên cứu lý thuyết và ứng dụng của đa thức Hermite và đa thức Laguerre liên kết trong cơ học lượng tử.
    • Xây dựng phần mềm hoặc gói công cụ tính toán số hoàn chỉnh phục vụ giảng dạy trực quan các hàm toán lý và mô phỏng hàm sóng trong chương trình Vật lý đại học và sau đại học.

Giá trị tham khảo

Khóa luận là tài liệu học tập và tham khảo có giá trị chuyên môn cao đối với:

  • Sinh viên ngành Sư phạm Vật lý, Vật lý học và Vật lý kỹ thuật: Là tài liệu hỗ trợ trực tiếp khi học các học phần Phương trình Toán lý, Điện động lực học, Cơ học lượng tử, Cơ học giải tích và Quang học sóng.
  • Sinh viên làm nghiên cứu khoa học, khóa luận tốt nghiệp: Cung cấp phương pháp giải chi tiết, rõ ràng từng bước cho các bài toán biên tĩnh điện và tán xạ sóng.
  • Giảng viên và giáo viên Vật lý: Cung cấp tư liệu bài tập, ví dụ minh họa và hệ thống công thức giải tích chuẩn xác để xây dựng bài giảng chuyên đề toán lý ứng dụng.

Câu hỏi thường gặp

1. Hàm sinh của đa thức Legendre được xây dựng từ bài toán vật lý cụ thể nào?
Hàm sinh đa thức Legendre được xây dựng từ bài toán tính thế tĩnh điện của một điện tích điểm $q$ đặt trên trục $z$ tại tọa độ $z = a$. Khi xét thế tại điểm $(r, \theta)$ trong miền $r > a$, áp dụng định lý cosin cho khoảng cách và khai triển chuỗi nhị thức lũy thừa theo biến $t = a/r$ và $x = \cos\theta$, biểu thức thế quy về dạng hàm sinh $g(t, x) = (1 - 2xt + t^2)^{-1/2} = \sum_{n=0}^{\infty} P_n(x)t^n$.

2. Thế tĩnh điện và mật độ điện tích cảm ứng của quả cầu dẫn điện trung hòa bán kính $r_0$ trong điện trường đều $E_0$ có kết quả ra sao?
Bằng cách giải phương trình Laplace trong hệ tọa độ cầu kết hợp điều kiện biên Dirichlet tại $r = r_0$ ($V = 0$) và tại $r \to \infty$ ($V \to -E_0 r\cos\theta$), nghiệm thế tĩnh điện ngoài quả cầu là $V(r, \theta) = -E_0 r P_1(\cos\theta)\left(1 - \frac{r_0^3}{r^3}\right)$. Mật độ điện tích mặt cảm ứng trên quả cầu là $\sigma = -\varepsilon_0 \left.\frac{\partial V}{\partial r}\right|_{r=r_0} = 3\varepsilon_0 E_0 \cos\theta$, và moment lưỡng cực cảm ứng tương ứng là $P = 4\pi\varepsilon_0 r_0^3 E_0$.

3. Ý nghĩa vật lý lượng tử của điều kiện $|m| \le n$ trong hàm Legendre liên kết $P_n^m(x)$ là gì?
Trong cơ học lượng tử, phương trình Legendre liên kết xuất hiện khi giải bài toán chuyển động trong trường thế xuyên tâm. Chỉ số $n$ liên quan đến trị riêng của bình phương moment động lượng toàn phần $L^2$ (với trị riêng $\hbar^2 n(n+1)$), còn chỉ số $m$ là hình chiếu của moment động lượng lên trục $z$, tức trị riêng của $L_z$ (với trị riêng $m\hbar$). Điều kiện toán học $|m| \le n$ (để đạo hàm cấp $m$ không triệt tiêu hàm) tương ứng chính xác với định luật vật lý: giá trị kỳ vọng bình phương hình chiếu $\langle L_z^2 \rangle$ không bao giờ vượt quá giá trị kỳ vọng của bình phương moment động lượng toàn phần $\langle L^2 \rangle$.

4. Tại sao hệ thức truy hồi của đa thức Legendre cần được viết lại khi tính toán số trên máy tính?
Hệ thức truy hồi chuẩn $(2n+1)xP_n(x) = (n+1)P_{n+1}(x) + nP_{n-1}(x)$ khi thực hiện tính toán số trực tiếp trên máy tính ở bậc $n$ lớn dễ gây tích tụ và khuếch đại sai số làm tròn. Để tăng độ ổn định số học, công thức được biến đổi tương đương thành $P_{n+1}(x) = 2xP_n(x) - P_{n-1}(x) - \frac{xP_n(x) - P_{n-1}(x)}{n+1}$, giúp thuật toán tính toán từ $P_0(x)=1, P_1(x)=x$ đạt độ chính xác cao hơn.

5. Khai triển sóng phẳng Rayleigh liên hệ giữa đa thức Legendre và hàm cầu Bessel như thế nào?
Phương trình Rayleigh biểu diễn một sóng phẳng lan truyền dọc trục thành chồng chất các sóng cầu: $e^{ikr\cos\gamma} = \sum_{n=0}^{\infty} a_n j_n(kr) P_n(\cos\gamma)$. Bằng cách sử dụng tính trực giao của đa thức Legendre $P_n(\cos\gamma)$, lấy đạo hàm cấp $n$ theo $(kr)$ và xét tại gốc $r = 0$, hệ số khai triển được xác định chính xác là $a_n = i^n (2n + 1)$.


Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Hữu Phát đã hoàn thành xuất sắc việc hệ thống hóa toàn diện cơ sở toán học và ứng dụng vật lý của hai họ hàm đặc biệt cơ bản là đa thức Legendre và hàm Bessel. Công trình không chỉ cung cấp các chứng minh giải tích chặt chẽ, các hệ thức truy hồi ổn định và nghiệm giải tích cho các bài toán biên kinh điển trong tĩnh điện học và cơ học lượng tử, mà còn xử lý hiệu quả khía cạnh tính toán số. Đây là một tài liệu học tập và nghiên cứu có giá trị thực tiễn cao, đáp ứng tốt yêu cầu nâng cao chất lượng đào tạo cho sinh viên chuyên ngành Sư phạm Vật lý.