Luận án Tiến sĩ Vật lý: Chuyển pha smectic-isotropic, tương tác vi mô tinh thể lỏng

Tổng quan về luận án

Nghiên cứu về tinh thể lỏng (Liquid Crystals - LCs) và các quá trình chuyển pha nhiệt động học đóng vai trò trung tâm trong vật lý vật chất ngưng tụ và khoa học vật liệu hiện đại. Tinh thể lỏng thermotropic sở hữu các pha trung gian độc đáo nằm giữa trạng thái rắn tinh thể và chất lỏng đẳng hướng, tiêu biểu là pha nematic, pha cholesteric và pha smectic. Trong đó, pha smectic đặc trưng bởi sự tự tổ chức của các phân tử hình khuẩn que (calamitic) thành các cấu trúc phân lớp song song, đồng thời sở hữu cả trật tự định hướng tầm xa (long-range orientational order) và trật tự vị trí một chiều (one-dimensional positional order). Tuy nhiên, cơ chế vi mô của hiện tượng chuyển pha từ pha smectic sang pha đẳng hướng (smectic - isotropic) vẫn là một bài toán mở đầy thách thức. Luận án tiến sĩ với tiêu đề "Nghiên cứu tính chất chuyển pha smectic - isotropic sử dụng tương tác vi mô trong cấu trúc tinh thể lỏng" của nghiên cứu sinh Nguyễn Thị Phương Thúy, chuyên ngành Vật lý lý thuyết và Vật lý toán tại Học viện Khoa học và Công nghệ - Viện Hàn lâm Khoa học và Công nghệ Việt Nam (2024), đã giải quyết trực diện khoảng trống học thuật này thông qua việc xây dựng mô hình vi mô và phát triển phương pháp tính toán số tiên tiến.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ sự không thống nhất giữa các quan sát thực nghiệm và các mô hình lý thuyết hiện hành. Thực nghiệm trên các màng mỏng smectic tự do (free-standing films) ghi nhận hiện tượng tan chảy từng lớp (layer-by-layer melting) rất đặc trưng theo quy luật hàm số mũ, nhưng các lý thuyết trường trung bình liên tục kinh điển như lý thuyết Landau-de Gennes không phản ánh được cấu trúc gián đoạn vi mô và động học phân tử. Ngược lại, các mô hình spin tĩnh truyền thống (như mô hình Ising, mô hình Lebwohl-Lasher) lại bỏ qua tính linh động tịnh tiến của các phân tử chất lỏng. Nghiên cứu tiên phong của GS. H. T. Diep (2020) đã bước đầu ứng dụng mô hình Potts linh động 6 trạng thái để mô phỏng pha smectic, nhưng thuật toán Metropolis tiêu chuẩn khi đó gặp phải hiện tượng bẫy tại các trạng thái siêu bền cục bộ (metastable states), khiến bản chất nhiệt động học tại lân cận nhiệt độ chuyển pha chưa được làm sáng tỏ hoàn toàn.

Để giải quyết triệt để vấn đề trên, luận án tập trung giải quyết hai câu hỏi nghiên cứu trọng tâm:

1. Câu hỏi nghiên cứu 1 (RQ1): Cơ chế tương tác vi mô dị hướng nào (giữa tương tác trao đổi nội lớp, tương tác liên lớp và thế tương tác hút - đẩy tầm ngắn) quyết định trật tự phân lớp và động học tan chảy từng lớp của màng smectic khi tăng nhiệt độ?
2. Câu hỏi nghiên cứu 2 (RQ2): Bằng cách nào có thể vượt qua rào cản năng lượng tự do trong chuyển pha loại 1 để tính toán chính xác hàm mật độ trạng thái (Density of States - DOS) và các hàm nhiệt động học của hệ tinh thể lỏng linh động?

Luận án thiết lập hai giả thuyết khoa học tương ứng:

• Giả thuyết 1 (H1): Việc tích hợp tương tác vi mô dị hướng ($J_k, J_\perp$) cùng thế năng Lennard-Jones ($V_0$) vào mô hình Potts linh động 6 trạng thái sẽ mô tả chính xác sự phân tách giữa quá trình mất trật tự vị trí phân lớp và mất trật tự định hướng phân tử.
• Giả thuyết 2 (H2): Sự kết hợp giữa thuật toán Metropolis cải tiến cơ chế dịch chuyển vị trí và kỹ thuật lấy mẫu vi chính tắc Wang-Landau sẽ triệt tiêu hiện tượng trễ nhiệt động, xác định chính xác các điểm chuyển pha bậc 1 thông qua cấu trúc phân bố năng lượng hai đỉnh ($P_L(E)$).

Phạm vi nghiên cứu bao quát các hệ mạng lập phương kích thước từ $N_s = 12^3$ ($1.728$ nút mạng), $18^3$ ($5.832$ nút mạng) đến $24^3$ ($13.824$ nút mạng), với nồng độ lấp đầy phân tử $c = N_s/N_L$ thay đổi linh hoạt từ $30%, 50%, 60%, 80%$ đến $100%$, khảo sát dưới các tỷ số tương tác trao đổi $J_k = 1.0, 1.5$, $J_\perp = \pm 0.5$ và độ sâu thế Lennard-Jones $V_0 = 0.2, 0.5, 1.0$. Ý nghĩa của nghiên cứu không chỉ dừng lại ở việc hoàn thiện bức tranh lý thuyết vi mô của vật lý chất mềm, mà còn cung cấp cơ sở định lượng để tối ưu hóa thiết kế vật liệu hiển thị quang điện tử, cảm biến nhiệt độ tinh thể lỏng và các hệ vi lưu tự lắp ghép.

Literature Review và Positioning

Lịch sử nghiên cứu tinh thể lỏng khởi đầu từ phát hiện của Friedrich Reinitzer (1888) trên hợp chất cholesteryl benzoate với hai điểm nóng chảy: điểm vẩn đục ($T_{mp} = 145.5^\circ\text{C}$) và điểm trong suốt ($T_{cp} = 178.5^\circ\text{C}$), được Otto Lehmann (1889) và Georges Friedel (1922) phân loại cấu trúc quang học thành các pha nematic, smectic và cholesteric. Nền tảng lý thuyết nhiệt động học tiếp tục được củng cố bởi lý thuyết biến dạng đàn hồi Frank (1958), lý thuyết trường phân tử Maier-Saupe (1958–1960) cho trật tự định hướng, và công trình đột phá của Pierre-Gilles de Gennes (Nobel Vật lý 1991) khi mở rộng lý thuyết chuyển pha Ginzburg-Landau để xây dựng phiếm hàm năng lượng tự do Landau-de Gennes cho hệ tinh thể lỏng.

Dòng nghiên cứu thực nghiệm về chuyển pha smectic - isotropic ghi nhận dấu mốc then chốt từ công trình của T. Stoebe và các cộng sự (1994) trên tạp chí Physical Review Letters. Bằng phương pháp đo hệ số phản xạ quang học $R = cN^2$ trên màng mỏng smectic tự do của hợp chất $H_{10}F_5MOPP$ và $H_8F_7MOPP$, nhóm nghiên cứu đã phát hiện quy luật tan chảy từng lớp độc đáo:

"Tại nhiệt độ $85^\circ\text{C}$, 10 lớp phân tử bị tan chảy, còn lại 15 lớp phân tử. Nhiệt độ hệ tiếp tục tăng lên, số lượng lớp phân tử còn lại lần lượt là 11, 9, 8, 7, 6, 5, 4, 3 và 2. Tại nhiệt độ $112^\circ\text{C}$, cấu trúc màng của pha smectic hoàn toàn bị phá vỡ."

Các nghiên cứu tiếp theo của Johnson et al. (1997) và Pankratz et al. (1999) đã xác lập quy luật phụ thuộc hàm số mũ phổ quát giữa số lớp $N$ và nhiệt độ chuyển pha:

"Số lớp phân tử $N$ tỉ lệ với nhiệt độ $t$ theo quy luật hàm số mũ: $N \propto t^{-\nu}$ với $\nu \approx 0.4$, trong đó $t = (T_0 - T_c(N))/T_0$."

Tuy nhiên, trong y văn quốc tế tồn tại những cuộc tranh luận lý thuyết sâu sắc chưa có lời giải cuối cùng:

1. Bản chất của cơ chế tan chảy lớp: Pankratz et al. (1999) cùng Géminard et al. (1997) giải thích sự tan chảy từng lớp thông qua cơ chế tạo mầm nhiệt động học (nucleation) và sự lan truyền của các vòng lệch mạng sai hỏng cấu trúc (dislocation loops) trong lòng khối. Ngược lại, Dolganov et al. (2022) đưa ra quan điểm đối lập rằng sự mất ổn định hình học tại mặt khum bề mặt (meniscus) liên kết với khung đỡ mới là nguồn gốc kích hoạt quá trình phá hủy lớp.
2. Tính phổ quát của hiện tượng: Một số tác giả từng giả định rằng hiện tượng tan chảy từng lớp chỉ tồn tại riêng biệt ở các hợp chất fluor hóa (fluorinated liquid crystals) do lực liên kết nội phân tử đặc thù. Giả thuyết này đã bị bác bỏ bởi Mirantsev et al. (1997) khi quan sát thấy quá trình tan chảy lớp tương tự trên hợp chất phi-fluor $54COOBC$, đòi hỏi phải có một mô hình vi mô mang tính quy luật phổ quát.
3. Định nghĩa tham số trật tự: Các nghiên cứu lý thuyết mở rộng Landau-de Gennes của Mukherjee et al. (2001, 2015) gặp khó khăn trong việc tách biệt độc lập giữa tham số trật tự định hướng $S$ và tham số trật tự vị trí $\psi_0$, dẫn đến việc không tái hiện được trạng thái trung gian linh động giữa các lớp.

Luận án này định vị chính xác tại điểm nghẽn của y văn: chuyển dịch từ cách tiếp cận hiện tượng học vĩ mô sang phương pháp mô hình hóa vi mô dựa trên mạng spin linh động. Bằng cách so sánh đối chiếu trực tiếp với dữ liệu thực nghiệm của Stoebe et al. (1994) và mô hình hạt cứng trụ cầu (spherocylinders) của Bolhuis & Frenkel (1997), luận án chứng minh rằng việc kết hợp Hamiltonian Potts dị hướng với kỹ thuật tính toán mật độ trạng thái Wang-Landau cho phép giải thích trọn vẹn cả tính chất chuyển pha loại 1 lẫn cơ chế động học tan chảy lớp mà không cần phụ thuộc vào các giả định tham số vĩ mô tùy ý.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng nền tảng lý thuyết mô hình Potts kinh điển (Ashkin-Teller 1943, Potts 1952, Wu 1982) từ trạng thái spin tĩnh sang hệ chất lỏng linh động có cấu trúc không gian ba chiều. Trong mô hình Potts tiêu chuẩn $q$ trạng thái, các spin cố định tại các nút mạng; khi $q > 4$ trong không gian hai chiều hoặc $q \ge 3$ trong không gian ba chiều, hệ chuyển pha trật tự - mất trật tự theo loại 1. Luận án đã kế thừa và nâng tầm mô hình Potts linh động 6 trạng thái do GS. H. T. Diep (2020) đề xuất bằng cách thiết lập một khung Hamiltonian vi mô hoàn chỉnh:

$$\mathcal{H} = - \sum_{\langle i,j \rangle} J_{ij} \delta(\sigma_i, \sigma_j) + \mathcal{H}_{ext}$$

Trong đó, mỗi phân tử được gán một biến trạng thái định hướng $\sigma_i \in {1, 2, 3, 4, 5, 6}$ tương ứng với 6 hướng trực giao trong không gian 3D. Tương tác trao đổi $J_{ij}$ được phân rã thành hai thành phần dị hướng rõ rệt: $J_\parallel = J_k$ biểu diễn liên kết giữa hai phân tử nằm trong cùng một mặt phẳng lớp $(xy)$, và $J_\perp$ biểu diễn liên kết phân tử giữa hai mặt phẳng kế cận dọc theo trục $z$.

Đóng góp đột phá về mặt lý thuyết của luận án thể hiện ở hai mệnh đề khoa học:

• Mệnh đề 1 (P1): Trật tự pha smectic được duy trì bền vững khi và chỉ khi năng lượng trao đổi nội lớp chiếm ưu thế ($J_k > J_\perp$). Tương tác liên lớp sắt từ ($J_\perp = +0.5 > 0$) bảo toàn trật tự đồng nhất giữa các lớp, trong khi tương tác phản sắt từ ($J_\perp = -0.5 < 0$) tạo ra sự phân cực ngược chiều, tái hiện trạng thái smectic phản sắt điện (antiferroelectric smectic phase).
• Mệnh đề 2 (P2): Khi tích hợp thế năng Lennard-Jones vào Hamiltonian vi mô: $$\mathcal{H}{LJ} = \sum{\langle i,j \rangle} 4\varepsilon \left[ \left(\frac{\sigma_{LJ}}{r_{ij}}\right)^{12} - \left(\frac{\sigma_{LJ}}{r_{ij}}\right)^6 \right]$$ sự cạnh tranh giữa lực đẩy lõi cứng tầm ngắn và lực hút van der Waals tầm xa điều khiển trực tiếp khoảng cách cân bằng giữa các lớp, giải thích cơ chế tan chảy bề mặt trước khi cấu trúc tinh thể lỏng sụp đổ hoàn toàn.

+-------------------------------------------------------------------------+
|                  KHUNG PHÂN TÍCH VI MÔ CỦA LUẬN ÁN                      |
+-------------------------------------------------------------------------+
|  [Cấu trúc pha vi mô]       [Hamiltonian tích hợp]    [Động học chuyển pha]
|  - Trật tự định hướng (q=6) -> H_Potts = -J*delta    -> Chuyển pha loại 1
|  - Trật tự phân lớp (z)     -> H_aniso = J_k + J_perp-> Tan chảy từng lớp
|  - Độ linh động vị trí (c)  -> H_LJ = 4e[(s/r)^12-(s/r)^6]-> Hóa lỏng Isotropic
+-------------------------------------------------------------------------+

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp đồng thời ba lý thuyết vật lý cốt lõi: Lý thuyết chuyển pha Landau (phân loại bậc chuyển pha qua tính gián đoạn của tham số trật tự), Cơ học thống kê vi chính tắc (sử dụng mật độ trạng thái để vượt qua rào cản năng lượng tự do), và Lý thuyết thế liên phân tử Lennard-Jones. Cách tiếp cận này tạo nên một cấu trúc phân tích đa tầng:

1. Tầng trật tự định hướng: Đo lường qua tham số trật tự chuẩn hóa $M = \frac{q \cdot \max(N_1, \dots, N_q)/N_s - 1}{q - 1}$ với $q = 6$.
2. Tầng trật tự phân lớp vị trí: Khảo sát sự phân bố mật độ phân tử theo từng mặt phẳng $z$ dọc theo trục vuông góc.
3. Tầng vi nhiệt động: Khảo sát năng lượng nội vi mô $U = \langle \mathcal{H} \rangle / N_s$, nhiệt dung riêng $C_v = \frac{\langle E^2 \rangle - \langle E \rangle^2}{N_s k_B T^2}$, và độ cảm thăng giáng tham số trật tự $\chi = \frac{N_s (\langle M^2 \rangle - \langle M \rangle^2)}{k_B T}$.

Điều kiện biên (boundary conditions) được xác lập chặt chẽ: áp dụng điều kiện biên tuần hoàn (Periodic Boundary Conditions - PBC) theo các phương $x, y$ để mô phỏng hệ vô hạn trong mặt phẳng lớp, và điều kiện biên tự do (Free Boundary Conditions) hoặc tuần hoàn theo phương $z$ để phân tích hiệu ứng màng mỏng và hiệu ứng kích thước thể tích.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án tuân thủ lập trường bản thể luận hiện thực phản biện (Critical Realism) kết hợp phương pháp luận thực chứng thống kê (Positivist Statistical Mechanics). Bản chất của các hiện tượng nhiệt động vĩ mô được xác định là hệ quả tất yếu xuất phát từ các quy luật tương tác và phân bố xác suất của các cấu hình vi mô. Để khắc phục triệt để các hạn chế của phương pháp mô phỏng Monte Carlo cổ điển đối với các chuyển pha có bước nhảy entropy lớn, luận án đã xây dựng một quy trình tính toán số hai giai đoạn:

+-------------------------------------------------------------------------+
|                  SƠ ĐỒ QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU                 |
+-------------------------------------------------------------------------+
|  [Khởi tạo mô hình]                                                    |
|  -> Thiết lập mạng L x L x L (L = 12, 18, 24)                          |
|  -> Phân bố N_s phân tử theo nồng độ c = 30%, 50%, 60%, 80%, 100%      |
|                                                                         |
|  [Giai đoạn 1: Thuật toán Metropolis cải tiến]                          |
|  -> Cải tiến hoán vị vị trí trống (hopping move) song song đổi spin q   |
|  -> Khảo sát cân bằng nhiệt động & hệ số khuếch tán D                   |
|                                                                         |
|  [Giai đoạn 2: Kỹ thuật Wang-Landau vi chính tắc]                      |
|  -> Lấy mẫu không định kiến tính hàm mật độ trạng thái g(E)            |
|  -> Kiểm tra tiêu chí phẳng Histogram: H(E) >= 80% <H(E)>              |
|  -> Tinh chỉnh hệ số biến đổi: f_{i+1} = sqrt(f_i) đến ln(f) ~ 10^-8   |
|                                                                         |
|  [Xử lý dữ liệu & Phân tích nhiệt động]                                |
|  -> Tích phân thống kê suy ra U(T), C_v(T), M(T), chi(T)               |
|  -> Phân tích hàm phân bố năng lượng P_L(E) tại nhiệt độ tới hạn       |
+-------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình mô phỏng được chuẩn hóa với độ tin cậy và tính lặp lại tuyệt đối:

1. Cải tiến thuật toán Metropolis: Trong thuật toán Metropolis truyền thống cho mô hình Potts linh động, việc thử lật spin và di chuyển vị trí diễn ra ngẫu nhiên độc lập thường dẫn đến xác suất chấp nhận cực thấp ở mật độ cao. Luận án đã tối ưu hóa quy trình: tại mỗi bước Monte Carlo (MCS), hệ thực hiện đồng thời phép thử thay đổi trạng thái định hướng $\sigma_i \rightarrow \sigma_i'$ với xác suất $W = \min(1, \exp(-\beta \Delta E))$ và phép thử nhảy vị trí (spatial displacement) sang một trong các nút mạng trống lân cận gần nhất (Nearest Neighbours - NN).
2. Kỹ thuật Wang-Landau tính mật độ trạng thái: Thay vì lấy mẫu theo phân bố Boltzmann truyền thống vốn dễ bị "bẫy" tại các cực tiểu cục bộ của năng lượng tự do, thuật toán Wang-Landau thực hiện bước đi ngẫu nhiên trong không gian năng lượng $E$ để xây dựng trực tiếp hàm mật độ trạng thái $g(E)$.
   • Ban đầu khởi tạo $g(E) = 1$ và hệ số biến đổi $f = f_0 = e^1 \approx 2.71828$.
   • Tại mỗi bước nhảy từ cấu hình có năng lượng $E_1$ sang $E_2$, xác suất chuyển dịch là: $$p(E_1 \rightarrow E_2) = \min\left(1, \frac{g(E_1)}{g(E_2)}\right)$$
   • Cập nhật hàm mật độ trạng thái: $\ln g(E_2) \leftarrow \ln g(E_2) + \ln f$ và tích lũy biểu đồ năng lượng $H(E_2) \leftarrow H(E_2) + 1$.
   • Tiêu chuẩn kiểm soát độ phẳng: Khi biểu đồ năng lượng đạt độ phẳng tối thiểu $H(E) \ge 80% \cdot \langle H(E) \rangle$ trên toàn bộ phổ năng lượng khảo sát, hệ số $f$ được tinh chỉnh theo quy tắc $f_{new} = \sqrt{f_{old}}$, đồng thời biểu đồ $H(E)$ được làm mới về 0. Quá trình lặp dừng lại khi $\ln f \le 10^{-8}$.

Data và phân tích

Dữ liệu mô phỏng được phân tích thông qua các kỹ thuật thống kê chuẩn xác:

• Sau khi thu được hàm mật độ trạng thái $g(E)$ chính xác, tổng thống kê chính tắc $Z(T)$ được tính toán liên tục tại mọi nhiệt độ mà không cần chạy lại mô phỏng: $$Z(T) = \sum_E g(E) \exp\left(-\frac{E}{k_B T}\right)$$
• Hàm phân bố xác suất năng lượng $P_L(E)$ tại nhiệt độ $T$ được trích xuất bằng biểu thức: $$P_L(E, T) = \frac{g(E) \exp(-E / k_B T)}{Z(T)}$$
• Phân tích hiệu ứng kích thước hữu hạn (Finite-Size Scaling - FSS) được thực hiện trên ba kích thước $L = 12, 18, 24$ để xác định nhiệt độ chuyển pha trong giới hạn nhiệt động học ($L \rightarrow \infty$).

Phát hiện đột phá và implications

Những phát hiện then chốt

Kết quả mô phỏng số chi tiết trong luận án đã mang lại 5 phát hiện mang tính đột phá về cấu trúc và nhiệt động học của chuyển pha smectic - isotropic:

1. Minh chứng thuyết phục về bản chất chuyển pha loại 1: Tại nhiệt độ chuyển pha, tham số trật tự $M$ suy giảm đột ngột từ giá trị bão hòa về tiệm cận 0, đi kèm với bước nhảy gián đoạn của năng lượng nội $U$. Đồ thị nhiệt dung riêng $C_v$ và độ cảm $\chi$ thể hiện các đỉnh nhọn phân kỳ đối xứng. Đặc biệt, hàm phân bố xác suất năng lượng $P_L(E)$ tại nhiệt độ tới hạn hiển thị cấu trúc hai đỉnh (bimodal energy distribution) rõ rệt với một rào cản năng lượng tự do ngăn cách ở giữa—bằng chứng nhiệt động học tuyệt đối của sự cùng tồn tại hai pha (phase coexistence) đặc trưng cho chuyển pha bậc 1.

2. Phát hiện hiện tượng phân tách đa chuyển pha theo nồng độ: Khi nồng độ phân tử giảm dần từ $c = 100%$ xuống $c = 30%$, hệ thể hiện sự phân rã của đỉnh nhiệt dung đơn lẻ thành chuỗi nhiều đỉnh nhiệt dung kế tiếp nhau:

   • Tại $c = 100%$ ($N_s = 12^3$): Hệ ghi nhận chuyển pha tập trung tại $T = 0.8790$.
   • Tại $c = 80%$: Chuyển pha xuất hiện đỉnh sắc nét tại $T_1 = 0.8320$.
   • Tại $c = 60%$: Đỉnh chuyển pha chính dịch chuyển về $T_1 = 0.6010$.
   • Tại $c = 50%$: Xuất hiện rõ rệt hai bước chuyển pha tại $T_1 = 0.7641$ (bắt đầu phá vỡ trật tự phân lớp) và $T_2 = 0.8796$ (phá hủy hoàn toàn trật tự định hướng).
   • Tại $c = 30%$: Hệ thể hiện cấu trúc 3 đỉnh chuyển pha phân tách tại các nhiệt độ $T_1 = 0.7253$, $T_2 = 0.7916$ và $T_3 = 0.8791$. Đây chính là sự tái hiện thành công về mặt vi mô của quá trình tan chảy tuần tự từng lớp từng được quan sát trong thực nghiệm của Stoebe et al. (1994).

3. Vai trò định hình cấu trúc của loại tương tác liên lớp ($J_\perp$): Khảo sát so sánh giữa tương tác hút (ferromagnetic $J_\perp = +0.5$) và tương tác đẩy (antiferromagnetic $J_\perp = -0.5$) tại $c = 50%$ cho thấy: với $J_\perp > 0$, các phân tử duy trì sự định hướng đồng pha giữa các lớp, hệ tan chảy tuần tự từ nhiệt độ bắt đầu $T = 0.747$ đến khi mất trật tự hoàn toàn ở $T = 0.949 > T_c$. Ngược lại, với $J_\perp < 0$, sự đối kháng định hướng giữa các lớp liền kề thúc đẩy quá trình trượt lớp và làm giảm nhiệt độ chuyển pha tổng thể.

4. Tác động điều chế của thế năng Lennard-Jones trong mô hình Potts mở rộng: Khi bổ sung đóng góp của thế năng Lennard-Jones vi mô ($V_0 = 1.0, 0.5, 0.2$), nhiệt độ chuyển pha của hệ bị dịch chuyển mạnh về vùng nhiệt độ thấp hơn. Cụ thể tại nồng độ $c = 30%$, mô hình Potts mở rộng ghi nhận chuỗi 6 điểm chuyển tiếp nhiệt động tinh vi tại các nhiệt độ: $$T_1 = 0.531146,; T_2 = 0.604178,; T_3 = 0.635270,; T_4 = 0.656481,; T_5 = 0.681227,; T_6 = 0.700268$$ Phát hiện này chứng minh rằng lực hút tầm xa van der Waals có xu hướng làm mềm mạng tinh thể và tạo điều kiện cho các lớp ngoài cùng tan chảy trước ở mức kích thích nhiệt thấp hơn.

5. Xác lập ngưỡng nồng độ giới hạn: Luận án chỉ ra giới hạn nồng độ vật lý tối thiểu để duy trì pha smectic:

   "Với nồng độ phân tử thấp $c = 20%$, các phân tử bay hơi hoàn toàn, nên không xác định được nhiệt độ chuyển pha." Dưới ngưỡng $c_{crit} \approx 20%$, mật độ liên kết vi mô không đủ để bù đắp dao động nhiệt entropy, khiến hệ chuyển thẳng từ trạng thái cụm phân tử phân tán sang pha khí-lỏng loãng mà không hình thành pha trung gian smectic.

Implications đa chiều

• Về mặt lý thuyết: Kết quả nghiên cứu cung cấp một cầu nối hoàn chỉnh giữa cơ học thống kê vi mô và lý thuyết hiện tượng học Landau-de Gennes. Mô hình đã chứng minh rằng hiện tượng tan chảy từng lớp không nhất thiết phải bắt nguồn từ các khuyết tật mạng phức tạp hay cấu trúc meniscus bề mặt, mà hoàn toàn có thể được kích hoạt thuần túy từ sự bất đối xứng giữa năng lượng tương tác nội lớp ($J_k$) và liên lớp ($J_\perp$).
• Về phương pháp luận: Quy trình kết hợp giữa thuật toán Metropolis cải tiến động học dịch chuyển và kỹ thuật Wang-Landau đã thiết lập một tiêu chuẩn mới cho việc mô phỏng các hệ chất lỏng phức tạp, polyme lỏng và các màng sinh học lipid có rào cản năng lượng tự do cao.
• Về ứng dụng thực tiễn: Các phát hiện về sự phụ thuộc của nhiệt độ chuyển pha vào độ dày lớp và nồng độ phân tử cung cấp thông số nền tảng để tinh chỉnh dải nhiệt độ hoạt động của các thiết bị hiển thị tinh thể lỏng (LCD), công tắc quang học (optical switches), và nâng cao độ nhạy của nhiệt kế tinh thể lỏng hoạt động trong dải $+18^\circ\text{C}$ đến $+34^\circ\text{C}$.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra 3 giới hạn nghiên cứu chính:

1. Giới hạn của mạng không gian rời rạc: Việc mô phỏng các phân tử trên mạng lập phương cố định ($L^3$) chưa phản ánh đầy đủ bậc tự do quay và tịnh tiến liên tục vô hạn của các phân tử calamitic thực tế trong không gian $O(3)$.
2. Xấp xỉ trạng thái định hướng rời rạc: Số lượng trạng thái định hướng trong mô hình Potts linh động được giới hạn ở $q = 6$ (trực giao), là một sự đơn giản hóa so với sự phân bố góc phương vị liên tục ($\theta, \phi$) trong các pha smectic C thực tế.
3. Quy mô kích thước hữu hạn: Do rào cản tài nguyên tính toán khi thực hiện thuật toán Wang-Landau đa thông số, kích thước mạng tối đa khảo sát đạt $N_s = 24^3$.

Chương trình nghiên cứu tương lai được định hình qua 4 hướng triển khai cụ thể:

• Mở rộng mô hình Potts linh động sang không gian tọa độ liên tục (off-lattice Monte Carlo simulation) kết hợp thế tương tác dị hướng Gay-Berne.
• Khảo sát đáp ứng động học của pha smectic dưới tác dụng của các trường ngoài định hướng (điện trường $\mathbf{E}$ và từ trường $\mathbf{B}$) nhằm tối ưu hóa thời gian đáp ứng milli-giây trong màn hình LCD thế hệ mới.
• Nghiên cứu hiện tượng chuyển pha trong các hệ tinh thể lỏng chiral smectic ($SmC^*$) sở hữu tính chất sắt điện (ferroelectricity) ứng dụng trong bộ nhớ quang-điện tử.
• Ứng dụng mô hình vi mô mở rộng để mô phỏng các màng kép phospholipid và các cấu trúc micelle trong tinh thể lỏng lyotropic phục vụ dẫn truyền thuốc hướng đích trong y sinh.

Tác động và ảnh hưởng

Nghiên cứu của luận án tạo ra những tác động khoa học và công nghệ rõ nét:

• Ảnh hưởng học thuật: Cung cấp bộ dữ liệu chuẩn về mật độ trạng thái và các điểm chuyển pha vi mô cho cộng đồng vật lý lý thuyết và vật lý tính toán. Các bài báo công bố từ luận án trên các tạp chí chuyên ngành quốc tế uy tín (ISI/Scopus) đóng góp trực tiếp vào việc giải quyết tranh luận kéo dài ba thập kỷ về cơ chế tan chảy lớp.
• Đổi mới công nghiệp quang điện tử: Dữ liệu về nhiệt độ chuyển pha và sự ổn định cấu trúc màng mỏng cung cấp cơ sở kỹ thuật cho các tập đoàn công nghệ sản xuất tấm nền hiển thị LCD cao cấp, giúp nâng cao độ tương phản và mở rộng dải nhiệt độ vận hành của màn hình từ môi trường cực lạnh đến môi trường nhiệt độ cao.
• Lợi ích xã hội và kinh tế: Việc tối ưu hóa quy trình tổng hợp vật liệu tinh thể lỏng dựa trên mô phỏng tính toán giúp giảm thiểu đáng kể chi phí thử nghiệm hóa học thực nghiệm, hạn chế phát thải hóa chất độc hại ra môi trường trong quá trình nghiên cứu phát triển vật liệu mới.

Đối tượng hưởng lợi

• Nghiên cứu sinh và Giảng viên Vật lý: Tiếp cận một khung lý thuyết và phương pháp tính toán số hoàn chỉnh về mô hình Potts linh động và kỹ thuật Wang-Landau cho các hệ chuyển pha phức tạp.
• Các nhà nghiên cứu Vật lý lý thuyết & Vật lý tính toán: Sử dụng thuật toán Metropolis cải tiến và các dữ liệu hàm phân bố năng lượng $P_L(E)$ làm hệ quy chiếu để kiểm chứng các mô hình thống kê mới.
• Kỹ sư R&D trong ngành công nghiệp vật liệu quang học: Khai thác các quy luật tương tác vi mô $J_k, J_\perp$ và thế Lennard-Jones để thiết kế các hợp chất tinh thể lỏng có dải pha smectic ổn định theo yêu cầu kỹ thuật.
• Các nhà khoa học công nghệ sinh học: Vận dụng cơ chế chuyển pha phân lớp để phân tích tính linh động và độ bền nhiệt của màng tế bào sinh học lyotropic.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?
Trả lời: Đóng góp độc đáo nhất là việc mở rộng thành công mô hình Potts $q$ trạng thái kinh điển của Ashkin-Teller (1943) và Potts (1952) thành Mô hình Potts Linh động Mở rộng 6 trạng thái tích hợp thế năng Lennard-Jones và tương tác trao đổi dị hướng ($J_k \ne J_\perp$). Khung lý thuyết này đã giải thích được bản chất vi mô của hiện tượng tan chảy từng lớp mà lý thuyết vĩ mô Landau-de Gennes không thể mô tả chi tiết.

2. Đổi mới phương pháp luận của luận án vượt trội hơn các nghiên cứu trước như thế nào?
Trả lời: So với nghiên cứu của GS. H. T. Diep (2020) chỉ sử dụng thuật toán Metropolis tiêu chuẩn dễ bị bẫy trạng thái cục bộ, luận án đã thực hiện hai cải tiến đột phá: (1) Tối ưu hóa cơ chế hoán vị vị trí trống song song với đổi hướng spin, và (2) Tích hợp kỹ thuật vi chính tắc Wang-Landau để tính trực tiếp hàm mật độ trạng thái $g(E)$. Phương pháp này loại bỏ hoàn toàn hiện tượng trễ nhiệt động học và cho phép xác định chính xác các điểm chuyển tiếp năng lượng cực nhỏ.

3. Phát hiện bất ngờ nhất thu được từ dữ liệu mô phỏng là gì?
Trả lời: Phát hiện bất ngờ nhất là sự phân tách của chuyển pha đơn lẻ thành một chuỗi đa chuyển pha tuần tự khi giảm nồng độ phân tử ($c = 30%$). Cụ thể, mô hình Potts mở rộng ghi nhận tới 6 điểm chuyển tiếp nhiệt động vi mô riêng biệt ($T_1 = 0.531146$ đến $T_6 = 0.700268$), chứng minh rằng hiện tượng tan chảy từng lớp là một chuỗi các quá trình tái cấu trúc nhiệt động học liên tiếp của các nhóm phân tử bề mặt trước khi toàn bộ hệ chuyển sang pha isotropic.

4. Luận án có cung cấp quy trình tái lập (replication protocol) rõ ràng không?
Trả lời: Có. Luận án cung cấp đầy đủ và chi tiết toàn bộ thông số mô phỏng: cấu hình mạng lập phương ($N_s = 12^3, 18^3, 24^3$), nồng độ lấp đầy ($c = 30% - 100%$), các giá trị hằng số ghép cặp vi mô ($J_k = 1.0, 1.5$; $J_\perp = \pm 0.5$), hệ số thế Lennard-Jones ($V_0 = 0.2, 0.5, 1.0$), tiêu chí phẳng histogram ($H(E) \ge 80%$), và tiêu chuẩn dừng thuật toán ($\ln f \le 10^{-8}$), đảm bảo khả năng tái lập độc lập 100% trên mọi hệ thống tính toán khoa học.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao?
Trả lời: Lộ trình 10 năm tới bao gồm việc chuyển đổi mô hình từ mạng rời rạc sang không gian liên tục (off-lattice continuum), khảo sát hiệu ứng từ-quang và điện-quang phi tuyến dưới tác dụng của trường ngoài cực mạnh, và mở rộng ứng dụng mô hình sang các cấu trúc màng sinh học phức hợp trong công nghệ nano-y sinh học.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Thị Phương Thúy đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra, mang lại 6 đóng góp cốt lõi cho chuyên ngành Vật lý lý thuyết và Vật lý toán:

1. Thiết lập thành công mô hình vi mô hoàn chỉnh cho pha smectic dựa trên mô hình Potts linh động 6 trạng thái dị hướng kết hợp thế tương tác Lennard-Jones.
2. Cải tiến và tối ưu hóa vượt bậc phương pháp mô phỏng Monte Carlo thông qua sự kết hợp giữa thuật toán Metropolis hoán vị và kỹ thuật Wang-Landau vi chính tắc.
3. Chứng minh dứt khoát bản chất chuyển pha loại 1 của quá trình chuyển pha smectic - isotropic thông qua sự xuất hiện của cấu trúc phân bố xác suất năng lượng hai đỉnh ($P_L(E)$) tại nhiệt độ tới hạn.
4. Tái hiện và giải thích thành công cơ chế động học của hiện tượng tan chảy tuần tự từng lớp—phù hợp hoàn hảo với các quy luật thực nghiệm quốc tế kinh điển của Stoebe et al. (1994).
5. Làm sáng tỏ vai trò điều khiển vi mô của tỷ số tương tác nội lớp/liên lớp ($J_k/J_\perp$) và thế tương tác hút - đẩy tầm ngắn ($V_0$) trong việc quyết định cấu trúc pha và độ ổn định nhiệt.
6. Xác định ngưỡng nồng độ vật lý tới hạn ($c_{crit} \approx 20%$) ngăn cách giữa trạng thái có thể hình thành pha tinh thể lỏng và trạng thái bay hơi phân tán hoàn toàn.

Công trình đánh dấu bước tiến quan trọng trong việc ứng dụng vật lý thống kê tính toán hiện đại để giải quyết các bài toán cấu trúc vật chất ngưng tụ phức tạp, mở ra những hướng nghiên cứu liên ngành đầy triển vọng trong khoa học vật liệu và công nghệ quang điện tử tiên tiến.