Mô hình Sherrington–Kirkpatrick khác gì so với mô hình Edwards–Anderson?

Mô hình Sherrington–Kirkpatrick giả thiết tương tác trao đổi vô hạn tầm giữa mọi cặp spin trong hệ thống, trong khi mô hình Edwards–Anderson chỉ xem xét tương tác giữa các spin láng giềng gần trong mạng tinh thể hữu hạn.

Tại sao nghiệm đối xứng bản sao ban đầu của Sherrington và Kirkpatrick bị coi là nghịch lý?

Nghiệm đối xứng bản sao dẫn đến kết quả entropy trên mỗi spin nhận giá trị âm xấp xỉ âm 0,16 hằng số Boltzmann ở nhiệt độ 0 tuyệt đối, vi phạm định luật thứ ba của nhiệt động lực học đối với hệ spin rời rạc.

Đường de Almeida–Thouless có ý nghĩa vật lý gì?

Đường de Almeida–Thouless phân định ranh giới ổn định của nghiệm đối xứng bản sao trong mặt phẳng nhiệt độ và từ trường ngoài. Dưới đường này, mode dao động replicon làm ma trận Hessian có trị riêng âm, khiến nghiệm đối xứng bản sao bị mất ổn định hoàn toàn.

Lý thuyết phá vỡ đối xứng bản sao của Giorgio Parisi đã giải quyết bài toán như thế nào?

Giorgio Parisi phân chia ma trận độ phủ giữa các bản sao thành các khối thứ bậc vô hạn, dẫn đến hàm tham số trật tự liên tục q(x), khôi phục entropy không âm ở nhiệt độ 0 và khám phá cấu trúc không gian trạng thái siêu mêtric.

Mô hình Sherrington–Kirkpatrick (SK model) là gì?

Mô hình Sherrington–Kirkpatrick là mô hình trường trung bình kinh điển trong cơ học thống kê mô tả pha thủy tinh spin với tương tác trao đổi vô hạn tầm và hỗn loạn ngẫu nhiên giữa các spin Ising. Mô hình này được xây dựng nhằm tìm kiếm một hệ thống khả giải chính xác về mặt giải tích để nghiên cứu các đặc tính nhiệt động lực học của trạng thái vô trật tự có tương tác cạnh tranh. Bài viết này trình bày chi tiết về nguồn gốc lịch sử, cấu trúc giải tích của Hamiltonian, nghiệm đối xứng bản sao cùng nghịch lý entropy âm, tiêu chuẩn mất ổn định de Almeida–Thouless, lý thuyết phá vỡ đối xứng bản sao của Giorgio Parisi, chứng minh toán học hiện đại, các ứng dụng liên ngành trong trí tuệ nhân tạo và các nghiên cứu liên quan tại Việt Nam.

Bối cảnh lịch sử và sự ra đời của mô hình

Trong vật lý chất rắn và vật lý vật liệu từ, trạng thái thủy tinh spin đại diện cho một trạng thái từ tính phức tạp xuất hiện trong các hợp kim kim loại loãng, nơi các mômen từ tạp chất tương tác gián tiếp với nhau thông qua cơ chế trao đổi Ruderman–Kittel–Kasuya–Yosida. Do vị trí không gian ngẫu nhiên của các nguyên tử từ tính trong mạng tinh thể kim loại nền, các liên kết tương tác giữa chúng nhận cả giá trị sắt từ (dương) và phản sắt từ (âm), gây ra hiện tượng cạnh tranh tương tác và làm tắc nghẽn trật tự từ tính thông thường.

Vào năm 1975, hai nhà vật lý lý thuyết David Sherrington và Scott Kirkpatrick đã công bố công trình nền tảng trên tập san Physical Review Letters (tập 35, số 26, trang 1792 đến 1796), giới thiệu một mô hình trường trung bình chính xác cho hệ thủy tinh spin. Thay vì giả định tương tác tầm ngắn chỉ giữa các láng giềng gần như mô hình Edwards–Anderson đề xuất trước đó, mô hình Sherrington–Kirkpatrick giả thiết rằng mọi cặp spin trong hệ thống đều tương tác trực tiếp với nhau bất kể khoảng cách hình học, tương tự như cách tiếp cận của mô hình Curie–Weiss đối với mô hình Ising truyền thống.

Mục tiêu ban đầu của Sherrington và Kirkpatrick là cung cấp một lời giải giải tích tường minh cho hiện tượng chuyển pha thủy tinh spin mà không phụ thuộc vào cấu trúc hình học của không gian mạng tinh thể hữu hạn. Tuy nhiên, quá trình phân tích mô hình này đã mở ra một cuộc khủng hoảng lý thuyết sâu sắc kéo dài nhiều năm trong cộng đồng vật lý thống kê quốc tế.

Công thức toán học và cơ sở vi mô

Xét một hệ gồm $N$ spin Ising cổ điển, trong đó mỗi biến số spin $\sigma_i$ tại vị trí nút mạng thứ $i$ chỉ nhận một trong hai giá trị rời rạc là $+1$ hoặc $-1$ . Hamiltonian toàn phần của mô hình Sherrington–Kirkpatrick được mô tả bằng biểu thức:

H = -\sum_{1 \le i < j \le N} J_{ij} \sigma_i \sigma_j - h \sum_{i=1}^N \sigma_i

Trong biểu thức Hamiltonian trên, đại lượng $J_{ij}$ là hằng số tương tác trao đổi giữa cặp spin thứ $i$ và spin thứ $j$ , còn $h$ là cường độ của từ trường ngoài đều tác dụng lên toàn bộ các spin trong hệ thống. Khác với mô hình Ising thuần nhất có hằng số trao đổi cố định, các giá trị tương tác $J_{ij}$ trong mô hình Sherrington–Kirkpatrick là các biến ngẫu nhiên độc lập, cùng tuân theo phân bố xác suất Gauss:

P(J_{ij}) = \sqrt{\frac{N}{2 \pi J^2}} \exp\left( -\frac{N (J_{ij} - J_0/N)^2}{2 J^2} \right)

Trong hàm phân bố xác suất trên, giá trị kỳ vọng trung bình của tương tác là $J_0/N$ và phương sai là $J^2/N$ , với $J$ là thang năng lượng tương tác đặc trưng của hệ. Để nghiên cứu thuần túy pha thủy tinh spin mà không làm xuất hiện trật tự sắt từ cảm ứng, các phân tích tiêu chuẩn thường chọn giá trị kỳ vọng trung bình bằng 0, nghĩa là $J_0 = 0$ .

Quy tắc tỉ lệ hóa phương sai tỷ lệ nghịch với số lượng hạt $N$ là điều kiện tiên quyết mang tính giải tích. Bởi vì mỗi spin liên kết đồng thời với $N-1$ spin khác, tổng số cặp tương tác trong hệ thống là $N(N-1)/2$ . Nếu phương sai tương tác không được giảm theo hệ số $1/N$ , tổng năng lượng của hệ sẽ tăng nhanh hơn bậc nhất theo $N$ , phá vỡ tính chất mở rộng của năng lượng tự do và ngăn cản việc xác lập giới hạn nhiệt động lực học khi $N$ tiến tới vô cùng.

Nghiệm đối xứng bản sao và nghịch lý nhiệt động lực học

Trong các hệ thống vô trật tự bị đóng băng (quenched disorder), việc tính toán năng lượng tự do đòi hỏi phải lấy trung bình của logarit của hàm phân bố thống kê $Z$ trên toàn bộ các cấu hình phân bố tương tác ngẫu nhiên $J_{ij}$ , được ký hiệu là $\mathbb{E}[\ln Z]$ . Để vượt qua độ phức tạp toán học khi lấy tích phân logarit trực tiếp, Sherrington và Kirkpatrick đã áp dụng phương pháp bản sao (replica trick) dựa trên đẳng thức giải tích:

\ln Z = \lim_{n \to 0} \frac{Z^n - 1}{n}

Phương pháp này quy việc tính toán trung bình của logarit về việc tính trung bình cấu hình của tích phân $n$ bản sao độc lập của hệ thống $\mathbb{E}[Z^n]$ , sau đó kéo dài giải tích tham số nguyên $n$ về giới hạn thực bằng 0.

Khi thực hiện phép tính tích phân phiếm hàm qua phương pháp điểm yên ngựa, Sherrington và Kirkpatrick đưa ra giả thiết đối xứng bản sao (replica-symmetric ansatz). Giả thiết này phát biểu rằng tất cả $n$ bản sao của hệ thống có vai trò hoàn toàn tương đương và đối xứng hình học với nhau dưới nhóm hoán vị, nghĩa là ma trận độ phủ $q_{\alpha \beta}$ giữa hai bản sao khác nhau $\alpha$ và $\beta$ nhận một giá trị duy nhất:

q_{\alpha \beta} = q \quad (\text{với mọi } \alpha \neq \beta)

Đại lượng $q$ trong biểu thức trên chính là tham số trật tự Edwards–Anderson, đo lường mức độ đóng băng định hướng của các mômen từ tại cùng một vị trí nút mạng theo thời gian. Khi nhiệt độ $T$ cao hơn nhiệt độ tới hạn $T_c = J$ (với quy ước hằng số Boltzmann bằng 1), hệ thống nằm trong pha thuận từ vô trật tự với tham số trật tự $q = 0$ . Khi nhiệt độ giảm xuống dưới nhiệt độ tới hạn $T_c$ , tham số trật tự $q$ nhận giá trị dương khác 0, đánh dấu sự chuyển pha sang trạng thái thủy tinh spin.

Mặc dù nghiệm đối xứng bản sao mô tả trơn tru đường chuyển pha ở nhiệt độ cao, nó dẫn đến một mâu thuẫn nhiệt động lực học nghiêm trọng ở vùng nhiệt độ thấp. Khi khảo sát giới hạn nhiệt độ tiến về 0 tuyệt đối, nghiệm đối xứng bản sao của Sherrington và Kirkpatrick cho ra kết quả entropy trên mỗi spin nhận giá trị âm, xấp xỉ bằng âm 0,16 hằng số Boltzmann (chính xác là $-1/(2\pi)$ ). Trong cơ học thống kê lượng tử và cổ điển của các hệ spin rời rạc, không gian trạng thái vi mô là hữu hạn nên entropy tuyệt đối bắt buộc phải lớn hơn hoặc bằng 0. Giá trị entropy âm này vi phạm trực tiếp định luật thứ ba của nhiệt động lực học, chứng tỏ rằng giả thiết đối xứng bản sao đã chứa đựng sai lầm căn bản.

Đường de Almeida–Thouless và giới hạn ổn định

Năm 1978, hai nhà vật lý J. R. L. de Almeida và D. J. Thouless đã giải quyết bí ẩn về nguồn gốc của giá trị entropy âm thông qua công trình xuất bản trên Journal of Physics A (tập 11, số 5, trang 983 đến 990). Bằng cách phân tích độ ổn định cục bộ của điểm yên ngựa trong không gian bản sao, de Almeida và Thouless đã tính toán ma trận đạo hàm bậc hai (ma trận Hessian) của phiếm hàm năng lượng tự do.

Phân tích đại số cho thấy ma trận Hessian trong không gian bản sao chứa ba trị riêng độc lập. Điểm yên ngựa đối xứng bản sao chỉ thực sự tương ứng với cực tiểu năng lượng tự do khi toàn bộ các trị riêng này đều dương. Tuy nhiên, de Almeida và Thouless phát hiện ra rằng trị riêng liên quan đến mode dao động phá vỡ tính đối xứng giữa các bản sao (được gọi là mode replicon) sẽ đổi dấu từ dương sang âm tại một ranh giới xác định trong mặt phẳng nhiệt độ và từ trường ngoài $(T, h)$ .

Ranh giới này được gọi là đường de Almeida–Thouless. Dưới đường de Almeida–Thouless, nghiệm đối xứng bản sao không còn là cực tiểu ổn định của năng lượng tự do mà trở thành một điểm yên ngựa không bền vững với hướng dốc âm. Khi từ trường ngoài $h$ có độ lớn nhỏ dần về 0, độ lệch nhiệt độ tới hạn tuân theo quy luật tiệm cận lũy thừa:

1 - \frac{T}{T_c} \propto h^{2/3}

Số mũ phân số 2/3 xuất hiện trong quy luật tiệm cận trên phản ánh đặc tính không giải tích của bề mặt năng lượng tự do dưới tác động phối hợp giữa từ trường ngoài và hiện tượng vô trật tự tương tác. Phát hiện về đường de Almeida–Thouless khẳng định dứt khoát rằng trong toàn bộ pha thủy tinh spin ở nhiệt độ thấp, tính đối xứng hoán vị giữa các bản sao bắt buộc phải bị phá vỡ.

Phương pháp phá vỡ đối xứng bản sao của Parisi

Năm 1979, nhà vật lý lý thuyết người Ý Giorgio Parisi đã công bố một giải pháp cách mạng trên Physical Review Letters (tập 43, số 23, trang 1754 đến 1756), thiết lập nên lý thuyết phá vỡ đối xứng bản sao (replica symmetry breaking). Thay vì gán một giá trị duy nhất cho mọi phần tử ngoài đường chéo của ma trận độ phủ $q_{\alpha \beta}$ , Parisi đề xuất một quy trình phân chia ma trận kích thước $n \times n$ thành các khối ma trận con lồng nhau một cách có thứ bậc.

Trong bước phá vỡ thứ nhất, ma trận $n \times n$ được chia thành các khối con kích thước $m_1 \times m_1$ . Các phần tử bên trong cùng một khối nhận giá trị $q_1$ , trong khi các phần tử nằm ngoài khối nhận giá trị $q_0$ . Trong các bước phá vỡ tiếp theo, mỗi khối con lại tiếp tục được phân tách thành các khối nhỏ hơn kích thước $m_2 \times m_2$ với giá trị độ phủ $q_2$ . Khi lặp lại quá trình phân chia này vô hạn lần và lấy đồng thời giới hạn $n \to 0$ , ma trận rời rạc $q_{\alpha \beta}$ biến đổi thành một hàm số liên tục $q(x)$ được định nghĩa trên đoạn thực $x \in [0, 1]$ .

Nghiệm phá vỡ đối xứng bản sao liên tục của Parisi đã giải quyết triệt để toàn bộ các nghịch lý nhiệt động lực học của mô hình Sherrington–Kirkpatrick:

Entropy trạng thái đáy: Khôi phục hoàn toàn tính không âm của entropy ở nhiệt độ 0 tuyệt đối, triệt tiêu nghịch lý vi phạm nhiệt động lực học của nghiệm cũ.
Năng lượng trạng thái đáy: Xác định giá trị năng lượng trạng thái đáy trên mỗi spin xấp xỉ bằng âm 0,7633 lần thang năng lượng $J$ , chính xác hơn đáng kể so với ước lượng âm 0,798 của nghiệm đối xứng bản sao.
Phân bố độ phủ: Hàm số $q(x)$ mô tả hàm phân bố xác suất của độ phủ giữa hai trạng thái cân bằng độc lập trong hệ thống, chỉ ra sự tồn tại đồng thời của vô số thung lũng năng lượng tự do không tương đương.
Cấu trúc cây siêu mêtric: Không gian các trạng thái cân bằng thuần nhất sắp xếp theo cấu trúc hình học cây siêu mêtric. Với ba trạng thái cân bằng bất kỳ $A$ , $B$ và $C$ , độ đo khoảng cách giữa chúng luôn thỏa mãn tính chất tam giác cân có cạnh đáy ngắn hơn hoặc bằng hai cạnh bên, nghĩa là độ phủ giữa chúng thỏa mãn hệ thức:

q_{AC} \ge \min(q_{AB}, q_{BC})

Chứng minh toán học nghiêm ngặt và ý nghĩa vật lý

Suốt hơn hai thập kỷ sau khi Parisi công bố lý thuyết, nghiệm phá vỡ đối xứng bản sao vẫn bị xem là một phương pháp trực giác đầy bí ẩn nhưng thiếu cơ sở toán học chặt chẽ. Việc kéo dài giải tích số chiều ma trận $n$ về giá trị 0 và đảo ngược thứ tự kích thước khối ma trận $0 \le m_k \le \dots \le m_1 \le n \le 1$ từng bị nhiều nhà toán học nghi ngờ về tính hợp lệ.

Bước ngoặt toán học đã diễn ra khi Francesco Guerra chứng lập được cận trên nghiêm ngặt cho năng lượng tự do thông qua kỹ thuật nội suy ngẫu nhiên, sau đó nhà toán học người Pháp Michel Talagrand đã hoàn thành chứng minh chặt chẽ cho công thức Parisi trên toàn bộ thang nhiệt độ. Công trình đồ sộ của Talagrand được công bố năm 2006 trên tập san Annals of Mathematics (tập 163, số 1, trang 221 đến 263), biến nghiệm Parisi từ một phỏng đoán vật lý trừu tượng thành một định lý toán học xác thực.

Thành tựu giải mã mô hình Sherrington–Kirkpatrick đã được tôn vinh ở cấp độ khoa học cao nhất. Năm 2021, Viện Hàn lâm Khoa học Hoàng gia Thụy Điển đã trao Giải Nobel Vật lý cho Giorgio Parisi vì những khám phá mang tính tiên phong về sự tác động qua lại giữa hỗn loạn và thăng giáng trong các hệ vật lý từ quy mô nguyên tử đến quy mô hành tinh. Trong bài giảng Nobel xuất bản năm 2023 trên tập san Reviews of Modern Physics (tập 95, số 3, mã bài 030501), Parisi nhấn mạnh rằng cấu trúc đa cân bằng và phương pháp luận bản sao được phát triển từ mô hình Sherrington–Kirkpatrick đã trở thành chiếc chìa khóa vạn năng mở ra sự hiểu biết cho lý thuyết các hệ phức tạp hiện đại.

So sánh với các mô hình hệ từ tính liên quan

Để làm rõ vị trí của mô hình Sherrington–Kirkpatrick trong bức tranh tổng thể của vật lý thống kê các hệ spin, bảng dưới đây đối chiếu các đặc tính cấu trúc và nhiệt động lực học then chốt giữa ba mô hình tiêu biểu:

Đặc tính so sánh	Mô hình Ising chuẩn	Mô hình Edwards–Anderson	Mô hình Sherrington–Kirkpatrick
Bản chất tương tác	Đồng nhất, sắt từ hoặc phản sắt từ	Hỗn loạn ngẫu nhiên có dấu thay đổi	Hỗn loạn ngẫu nhiên phân bố Gauss
Phạm vi tương tác	Láng giềng gần trong mạng tinh thể hữu hạn	Láng giềng gần trong mạng tinh thể hữu hạn	Vô hạn tầm, liên kết toàn bộ các cặp spin
Bản chất trường trung bình	Chỉ đúng khi số chiều không gian tiến tới vô cùng	Không phải trường trung bình	Mô hình trường trung bình chính xác
Nghẽn tương tác	Không xuất hiện hiện tượng nghẽn tương tác	Xuất hiện nghẽn tương tác cục bộ	Nghẽn tương tác toàn cục cực đại
Số lượng trạng thái đáy	Hai trạng thái suy biến đối xứng	Chưa xác định chính xác trong không gian 3 chiều	Vô số trạng thái cân bằng tổ chức siêu mêtric
Khả năng giải giải tích	Giải được chính xác ở số chiều 1 và 2	Không thể giải chính xác giải tích	Giải được chính xác qua nghiệm Parisi

Ứng dụng liên ngành và các hướng nghiên cứu tại Việt Nam

Mặc dù xuất phát từ bài toán từ tính chất rắn thuần túy, mô hình Sherrington–Kirkpatrick đã trở thành nền tảng toán học cốt lõi cho hàng loạt lĩnh vực khoa học liên ngành hiện đại:

Mạng thần kinh nhân tạo: Mô hình mạng Hopfield do John Hopfield đề xuất về bộ nhớ liên kết thực chất là một biến thể trực tiếp của mô hình Sherrington–Kirkpatrick với ma trận liên kết Hebbian. Các kỹ thuật bản sao của Parisi cho phép tính toán chính xác dung lượng lưu trữ cực đại của mạng thần kinh trước khi xảy ra hiện tượng quá tải thông tin.
Tối ưu hóa tổ hợp và khoa học máy tính: Các bài toán tối ưu hóa tổ hợp khó thuộc lớp NP-đầy đủ như bài toán người bán hàng hoặc bài toán phân vùng đồ thị có bề mặt hàm mục tiêu tương đương với cảnh quan năng lượng tự do gồ ghề của thủy tinh spin. Thuật toán tôi luyện mô phỏng (simulated annealing) khai thác trực tiếp cơ chế chuyển động nhiệt trên các thung lũng năng lượng của mô hình này để tìm nghiệm gần tối ưu.
Lý thuyết học máy và mạng nơ-ron sâu: Trong học sâu đương đại, cảnh quan hàm mất mát của các mạng nơ-ron sâu tham số lớn biểu hiện các tính chất tương đồng kỳ lạ với cảnh quan năng lượng của mô hình Sherrington–Kirkpatrick, bao gồm sự tồn tại của vô số điểm yên ngựa và các thung lũng phẳng giúp việc huấn luyện hội tụ hiệu quả.

Tại Việt Nam, các hướng nghiên cứu lý thuyết và mô phỏng số về hệ vô trật tự, mô hình Ising và thủy tinh spin luôn thu hút sự quan tâm của các nhóm nghiên cứu vật lý lý thuyết và vật lý tính toán. Tiêu biểu vào năm 2013, nhóm nghiên cứu gồm Ngô Văn Thanh, Hoàng Đình Tiến và Điệp Thế Hùng tại Viện Vật lý (Viện Hàn lâm Khoa học và Công nghệ Việt Nam) phối hợp cùng Đại học Cergy-Pontoise đã công bố công trình chuyên sâu trên tập san Communications in Physics (tập 23, số 3, trang 227 đến 234). Công trình đã khảo sát hành vi chuyển pha và hiện tượng trật tự thủy tinh spin trong các hệ Ising bị pha tạp và thất vọng từ tính trên mạng tinh thể lập phương tâm diện (FCC) thông qua phương pháp mô phỏng Monte Carlo, đóng góp vào việc làm sáng tỏ cơ chế vi mô của các trạng thái từ tính cạnh tranh trong vật liệu từ thực tế.

Hạn chế và các câu hỏi mở

Dù mô hình Sherrington–Kirkpatrick là một chuẩn mực thành công của lý thuyết trường trung bình, nó vẫn tồn tại những ranh giới và giới hạn vật lý rõ rệt:

Giả định tương tác vô hạn tầm: Trong vật liệu từ thực nghiệm, tương tác trao đổi suy giảm theo lũy thừa khoảng cách không gian. Do đó, mô hình Sherrington–Kirkpatrick không phản ánh các thăng giáng không gian cục bộ xảy ra trong các hệ ba chiều thực tế.
Tranh luận giữa lý thuyết trường trung bình và mô hình giọt: Trong suốt nhiều thập kỷ qua, câu hỏi liệu tính chất phá vỡ đối xứng bản sao của nghiệm Parisi có thực sự tồn tại trong các hệ thủy tinh spin không gian ba chiều hữu hạn có từ trường hay không vẫn là một chủ đề tranh luận học thuật chưa có hồi kết giữa trường phái trường trung bình Parisi và trường phái mô hình giọt (droplet model) của Daniel Fisher và David Huse.
Động lực học thủy tinh spin: Các phân tích giải tích của mô hình chủ yếu tập trung vào trạng thái cân bằng nhiệt động lực học. Quá trình hồi phục nhiệt động học phi cân bằng, hiện tượng lão hóa (aging) và trí nhớ từ tính trong thang thời gian vĩ mô vẫn là những thách thức tính toán phức tạp đang tiếp tục được nghiên cứu sâu rộng.