Tính Toán Biến Phân và Ứng Dụng trong Vật Lý

1. Giới thiệu

Tính toán biến phân (Calculus of Variations) là một lĩnh vực toán học nghiên cứu về việc tìm các hàm số tối ưu hóa (cực trị) một hàm chức (functional). Hàm chức là một hàm số nhận đầu vào là một hàm số và trả về một giá trị số thực.

Trong vật lý, tính toán biến phân đóng vai trò quan trọng trong việc thiết lập các nguyên lý vật lý cơ bản, đặc biệt là nguyên lý tác dụng tối thiểu trong cơ học cổ điển và thuyết tương đối.

2. Hàm Chức và Cực Trị Hàm Chức

2.1. Hàm Chức

Một hàm chức, thường ký hiệu là $J$, là một ánh xạ từ một không gian hàm số vào tập số thực $\mathbb{R}$:

$$J[y] = \int_{x_1}^{x_2} L(x, y(x), y'(x)) dx$$

Trong đó:

• $y(x)$ là một hàm số khả vi.
• $y'(x) = \frac{dy}{dx}$ là đạo hàm của $y(x)$.
• $L(x, y, y')$ là một hàm số cho trước, gọi là hàm Lagrangian.
• $x_1$ và $x_2$ là các giới hạn tích phân.

2.2. Cực Trị Hàm Chức

Mục tiêu của tính toán biến phân là tìm hàm số $y(x)$ sao cho hàm chức $J[y]$ đạt cực trị (cực đại hoặc cực tiểu). Để tìm cực trị của hàm chức, ta cần tìm hàm số $y(x)$ thỏa mãn phương trình Euler-Lagrange.

3. Phương trình Euler-Lagrange

3.1. Thiết lập Phương Trình

Giả sử $y(x)$ là hàm số làm cho hàm chức $J[y]$ đạt cực trị. Xét một biến thiên nhỏ $\delta y(x)$ của hàm $y(x)$ sao cho $\delta y(x_1) = \delta y(x_2) = 0$. Khi đó, biến thiên của hàm chức $J$ là:

$$\delta J = J[y + \delta y] - J[y] = \int_{x_1}^{x_2} \left[ L(x, y + \delta y, y' + \delta y') - L(x, y, y') \right] dx$$

Sử dụng khai triển Taylor cho hàm $L$ đến bậc nhất, ta có:

$$L(x, y + \delta y, y' + \delta y') \approx L(x, y, y') + \frac{\partial L}{\partial y} \delta y + \frac{\partial L}{\partial y'} \delta y'$$

Do đó:

$$\delta J \approx \int_{x_1}^{x_2} \left( \frac{\partial L}{\partial y} \delta y + \frac{\partial L}{\partial y'} \delta y' \right) dx$$

Áp dụng tích phân từng phần cho số hạng thứ hai:

$$\int_{x_1}^{x_2} \frac{\partial L}{\partial y'} \delta y' dx = \left[ \frac{\partial L}{\partial y'} \delta y \right]_{x_1}^{x_2} - \int_{x_1}^{x_2} \frac{d}{dx} \left( \frac{\partial L}{\partial y'} \right) \delta y dx$$

Vì $\delta y(x_1) = \delta y(x_2) = 0$, số hạng đầu tiên bằng 0. Vậy:

$$\delta J \approx \int_{x_1}^{x_2} \left[ \frac{\partial L}{\partial y} - \frac{d}{dx} \left( \frac{\partial L}{\partial y'} \right) \right] \delta y dx$$

Để $J$ đạt cực trị, biến thiên $\delta J$ phải bằng 0 với mọi $\delta y$. Điều này dẫn đến phương trình Euler-Lagrange:

$$\frac{\partial L}{\partial y} - \frac{d}{dx} \left( \frac{\partial L}{\partial y'} \right) = 0$$

3.2. Dạng Tổng Quát

Phương trình Euler-Lagrange có thể được viết lại dưới dạng:

$$\frac{\partial L}{\partial y} - \frac{d}{dx}\left(\frac{\partial L}{\partial y'}\right) = \frac{\partial L}{\partial y} - \frac{\partial^2 L}{\partial x \partial y'} - \frac{\partial^2 L}{\partial y \partial y'}y' - \frac{\partial^2 L}{\partial y'^2}y'' = 0$$

4. Ứng Dụng trong Cơ Học Cổ Điển: Nguyên Lý Tác Dụng Tối Thiểu

4.1. Tác Dụng (Action)

Trong cơ học cổ điển, tác dụng $S$ là một hàm chức được định nghĩa là tích phân của hàm Lagrangian $L$ theo thời gian:

$$S = \int_{t_1}^{t_2} L(q(t), \dot{q}(t), t) dt$$

Trong đó:

• $q(t)$ là tọa độ tổng quát của hệ.
• $\dot{q}(t) = \frac{dq}{dt}$ là vận tốc tổng quát.
• $L = T - V$ là hàm Lagrangian, với $T$ là động năng và $V$ là thế năng.

4.2. Nguyên Lý Tác Dụng Tối Thiểu (Principle of Least Action)

Nguyên lý tác dụng tối thiểu phát biểu rằng quỹ đạo thực của hệ cơ học là quỹ đạo mà tác dụng $S$ đạt cực trị (thường là cực tiểu). Điều này có nghĩa là, quỹ đạo thực là nghiệm của phương trình Euler-Lagrange:

$$\frac{\partial L}{\partial q} - \frac{d}{dt} \left( \frac{\partial L}{\partial \dot{q}} \right) = 0$$

Phương trình này tương đương với phương trình Euler-Lagrange đã được thiết lập ở phần trước, chỉ thay biến $x$ bằng $t$ và $y$ bằng $q$.

4.3. Phương Trình Lagrange

Phương trình Euler-Lagrange trong cơ học cổ điển thường được gọi là phương trình Lagrange. Nếu $q$ là một tọa độ tổng quát và $L = T - V$, phương trình Lagrange trở thành:

$$\frac{d}{dt} \left( \frac{\partial L}{\partial \dot{q}} \right) - \frac{\partial L}{\partial q} = 0$$

Ví dụ: Xét một hạt có khối lượng $m$ chuyển động trong trường thế $V(x)$. Động năng là $T = \frac{1}{2} m \dot{x}^2$. Hàm Lagrangian là:

$$L = T - V = \frac{1}{2} m \dot{x}^2 - V(x)$$

Áp dụng phương trình Lagrange:

$$\frac{\partial L}{\partial x} = -\frac{\partial V}{\partial x}$$ $$\frac{\partial L}{\partial \dot{x}} = m \dot{x}$$ $$\frac{d}{dt} \left( \frac{\partial L}{\partial \dot{x}} \right) = m \ddot{x}$$

Phương trình Lagrange trở thành:

$$m \ddot{x} + \frac{\partial V}{\partial x} = 0$$

Đây chính là định luật 2 Newton: $F = -\frac{\partial V}{\partial x} = m \ddot{x}$.

5. Ứng Dụng trong Thuyết Tương Đối

5.1. Tác Dụng trong Thuyết Tương Đối Hẹp

Trong thuyết tương đối hẹp, tác dụng cho một hạt tự do có khối lượng $m$ là:

$$S = -mc \int ds$$

Trong đó:

• $m$ là khối lượng nghỉ của hạt.
• $c$ là vận tốc ánh sáng.
• $ds = \sqrt{c^2 dt^2 - dx^2 - dy^2 - dz^2}$ là phần tử độ dài không thời gian (spacetime interval).

5.2. Hàm Lagrangian Tương Đối Tính

Có thể viết lại tác dụng dưới dạng tích phân theo thời gian:

$$S = \int L dt$$

Với hàm Lagrangian tương đối tính:

$$L = -mc^2 \sqrt{1 - \frac{v^2}{c^2}}$$

Trong đó $v^2 = \dot{x}^2 + \dot{y}^2 + \dot{z}^2$ là bình phương vận tốc của hạt.

5.3. Phương Trình Chuyển Động Tương Đối Tính

Áp dụng phương trình Euler-Lagrange cho từng tọa độ (ví dụ $x$):

$$\frac{\partial L}{\partial x} - \frac{d}{dt} \left( \frac{\partial L}{\partial \dot{x}} \right) = 0$$

Tính các đạo hàm:

$$\frac{\partial L}{\partial x} = 0$$ $$\frac{\partial L}{\partial \dot{x}} = \frac{m \dot{x}}{\sqrt{1 - \frac{v^2}{c^2}}} = \gamma m \dot{x}$$

Trong đó $\gamma = \frac{1}{\sqrt{1 - \frac{v^2}{c^2}}}$ là hệ số Lorentz.

Phương trình Euler-Lagrange trở thành:

$$\frac{d}{dt} (\gamma m \dot{x}) = 0$$

Tương tự cho các tọa độ $y$ và $z$. Kết quả là động lượng tương đối tính được bảo toàn.

6. Tổng Kết

Tính toán biến phân là một công cụ mạnh mẽ trong vật lý, cho phép ta thiết lập các nguyên lý vật lý cơ bản dựa trên nguyên lý cực trị. Phương trình Euler-Lagrange là phương trình quan trọng nhất trong tính toán biến phân, được sử dụng để tìm các hàm số tối ưu hóa hàm chức. Ứng dụng của tính toán biến phân rất đa dạng, từ cơ học cổ điển (nguyên lý tác dụng tối thiểu) đến thuyết tương đối. Hiểu rõ về tính toán biến phân giúp ta có cái nhìn sâu sắc hơn về cấu trúc toán học của các định luật vật lý.