Tổng quan về luận án

Nghiên cứu cấu trúc và động học chuyển pha của tinh thể lỏng (Liquid Crystals - LCs) đóng vai trò nền tảng trong vật lý chất rắn, vật lý thống kê hiện đại và công nghệ quang điện tử. Trong số các pha trung gian (mesophases), pha smectic thermotropic sở hữu mức độ trật tự cao khi các phân tử hình khuẩn que (calamitic molecules) vừa duy trì trật tự định hướng tầm xa (orientational order) theo trục ưu tiên $\vec{n}$ (director), vừa tổ chức thành các mặt phẳng lớp song song với trật tự vị trí định kỳ (positional order). Quá trình chuyển pha smectic - isotropic (Sm-I) – chuyển đổi trực tiếp từ trạng thái tinh thể lỏng phân lớp sang trạng thái lỏng đẳng hướng vô trật tự – là một bài toán phức tạp bậc nhất của cơ học thống kê.

Mặc dù thực nghiệm của T. Stoebe và cộng sự (1994) trên tạp chí Physical Review Letters đã phát hiện hiệu ứng tan chảy từng lớp (layer-by-layer melting) độc đáo trên màng tự do $H10F5MOPP$ và $H8F7MOPP$ tuân theo quy luật hàm số mũ $N \propto t^{-\nu}$ ($\nu \approx 0.4$), các mô hình lý thuyết trường trung bình liên tục (continuum mean-field theories) như lý thuyết Landau-de Gennes vẫn chưa thể mô tả trọn vẹn bản chất vi mô của động học tan chảy này. Hạn chế cốt lõi nằm ở việc thiếu vắng một mô hình mạng vi mô rời rạc kết hợp đồng thời tương tác spin định hướng, bậc tự do tịnh tiến và tương tác đẩy - hút thực tế giữa các tâm phân tử.

Luận án tiến sĩ của Nghiên cứu sinh Nguyễn Thị Phương Thúy 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" (Chuyên ngành Vật lý lý thuyết và Vật lý toán, Mã số 9440103; 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; hướng dẫn bởi PGS. Ngô Văn Thanh và TS. Nguyễn Thị Lâm Hoài, với sự cố vấn của GS. H. Diep - Đại học Cergy, Pháp) đã giải quyết triệt để khoảng trống nghiên cứu này.

Luận án đặt ra ba câu hỏi nghiên cứu trọng tâm:

  1. RQ1: Làm thế nào để mô hình hóa đồng thời tính linh động tịnh tiến và tương tác dị hướng vi mô của phân tử tinh thể lỏng trong pha smectic trên một mạng rời rạc?
  2. RQ2: Cơ chế vi mô nào chi phối quá trình tan chảy cấu trúc lớp và bản chất nhiệt động học của chuyển pha Sm-I tại các dải nồng độ phân tử khác nhau?
  3. RQ3: Sự kết hợp giữa tương tác trao đổi spin tầm ngắn và thế năng tương tác cặp Lennard-Jones làm thay đổi bức tranh chuyển pha Sm-I và độ bền nhiệt của hệ như thế nào?

Ba giả thuyết khoa học tương ứng được xác lập:

  • H1: Mô hình Potts linh động $6$ trạng thái ($q=6$ Mobile Potts model) với tương tác dị hướng giữa các lớp ($J_k, J_\perp$) có thể tái lập chính xác trật tự lớp smectic và trật tự định hướng không gian.
  • H2: Chuyển pha Sm-I mang bản chất chuyển pha loại 1 (first-order phase transition), đặc trưng bởi sự tồn tại đồng thời của hai pha và xuất hiện rào cản năng lượng tự do tại nhiệt độ tới hạn $T_c$.
  • H3: Việc bổ sung thế năng Lennard-Jones ($V_0$) vào Hamiltonian sẽ ngăn chặn hiện tượng bay hơi sớm của các phân tử ở nồng độ thấp ($c \le 30%$), tái hiện cơ chế phân rã từng lớp phù hợp với thực nghiệm.

Khung lý thuyết của luận án tích hợp mô hình Potts linh động, cơ học thống kê Boltzmann, thuật toán Metropolis cải tiến và kỹ thuật lấy mẫu mật độ trạng thái Wang-Landau (Wang-Landau density of states sampling). Nghiên cứu được thực thi trên các kích thước mạng tinh thể lập phương $N_s = 12^3, 18^3, 24^3$ với dải nồng độ phân tử $c = 30% - 100%$, tạo bước đột phá trong việc làm sáng tỏ cơ chế vi mô của hiện tượng chuyển pha Sm-I.

Literature Review và Positioning

Nghiên cứu chuyển pha tinh thể lỏng khởi đầu từ quan sát thực nghiệm năm 1888 của Friedrich Reinitzer trên hợp chất cholesteryl benzoate và các khảo sát quang học phân cực của Otto Lehmann (1889). Năm 1922, Georges Friedel phân loại chính thức các pha mesophase thành nematic, smectic và cholesteric dựa trên trật tự phân tử. Pierre-Gilles de Gennes (Nobel Vật lý 1991) đã mở rộng lý thuyết chuyển pha Ginzburg-Landau sang hệ tinh thể lỏng, định nghĩa các tham số trật tự tensor cho pha nematic và tham số trật tự phức cho pha smectic.

Trong dòng nghiên cứu thực nghiệm về chuyển pha Sm-I, công trình bước ngoặt của Stoebe et al. (1994) trên màng mỏng 25 lớp phân tử đo hệ số phản xạ quang học $R = cN^2$ đã chỉ ra quá trình tan chảy tuần tự: tại $85^\circ\text{C}$, 10 lớp ngoài cùng tan rã; khi nhiệt độ tăng dần, số lớp giảm tiếp tục còn 11, 9, 8, 7, 6, 5, 4, 3, 2 trước khi phá vỡ hoàn toàn ở $112^\circ\text{C}$. Các thực nghiệm tiếp theo của Chao et al. (1996) bằng nhiễu xạ electron, Johnson et al. (1997) và Mirantsev (1995) trên hợp chất không chứa fluor ($54COOBC$) xác nhận tính phổ quát của hiện tượng tan chảy theo lớp.

                      ┌────────────────────────────────────────┐
                      │    LÝ THUYẾT TRƯỜNG LIÊN TỤC           │
                      │  Landau-de Gennes / Frank Elastic /    │
                      │  Maier-Saupe (1958 - 1991)             │
                      └──────────────────┬─────────────────────┘
                                         │ (Chưa giải thích vi mô)
                                         ▼
┌────────────────────────────────────────┴────────────────────────────────────────┐
│                        TRANH LUẬN HỌC THUẬT VỀ CƠ CHẾ                           │
├────────────────────────────────────────┬────────────────────────────────────────┤
│ Trường phái Tạo mầm Vòng lệch mạng     │ Trường phái Biến dạng Bề mặt Mặt khum  │
│ Pankratz (1999) & Géminard (1997)      │ Dolganov et al. (2022)                 │
└────────────────────────────────────────┬────────────────────────────────────────┘
                                         │ (Cần mô hình vi mô rời rạc)
                                         ▼
                      ┌────────────────────────────────────────┐
                      │    MÔ HÌNH MẠNG VI MÔ LINH ĐỘNG        │
                      │  GS. H. Diep et al. (2020)             │
                      └──────────────────┬─────────────────────┘
                                         │ (Hạn chế: Chưa rõ động học lân cận Tc)
                                         ▼
                      ┌────────────────────────────────────────┐
                      │         ĐÓNG GÓP CỦA LUẬN ÁN           │
                      │  Mô hình Potts Mở rộng + Wang-Landau   │
                      │  Tái hiện chính xác Tan chảy Từng lớp  │
                      └────────────────────────────────────────┘

Về mặt cơ chế, y văn quốc tế tồn tại hai luồng quan điểm đối lập:

  • Trường phái 1 (Pankratz 1999; Géminard 1997; de Jeu 2015): Cho rằng sự tan chảy theo lớp xuất phát từ quá trình tạo mầm nhiệt động (nucleation) và sự phát triển của các vòng lệch mạng rối loạn (dislocation loops) do mất cân bằng giữa lực căng bề mặt và lực đàn hồi.
  • Trường phái 2 (Dolganov et al., 2022): Khẳng định biến dạng bề mặt mặt khum (meniscus deformation) và gradient thế hóa học ở biên mới là nguyên nhân kích hoạt tan rã cấu trúc màng.

Về phương pháp mô phỏng, các mô hình spin truyền thống (Lebwohl-Lasher 1972) chỉ gán spin cố định tại các nút mạng, bỏ qua bậc tự do tịnh tiến. Mô hình hạt cứng (spherocylinder) hay Gay-Berne (1981) tiêu tốn tài nguyên tính toán cực lớn, khó khảo sát các kích thước hệ thống đủ lớn để phân tích hiệu ứng kích thước hữu hạn. Năm 2020, GS. H. Diep và cộng sự đề xuất mô hình Potts linh động 6 trạng thái sử dụng thuật toán Metropolis; tuy nhiên, mô phỏng ban đầu gặp hiện tượng thăng giáng giả tạo và chưa khảo sát được quá trình tan chảy từng lớp sắc nét. Luận án định vị chính xác vào giao điểm này: kết hợp mô hình Potts vi mô với thế tương tác Lennard-Jones và phương pháp lấy mẫu mật độ trạng thái Wang-Landau, tái hiện trực tiếp bức tranh tan rã theo lớp từ các tương tác vi mô cơ bản.

Đó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 (Julius Ashkin & Edward Teller 1943; Renfrey Potts 1951; F. Y. Wu 1982) từ mạng cố định sang mạng linh động có chỗ trống (dilute mobile lattice), gán $q=6$ trạng thái định hướng tương ứng với 6 hướng không gian 3D ($\pm Ox, \pm Oy, \pm Oz$).

Hamiltonian của hệ được thiết lập dưới dạng vi mô dị hướng:

$$H_1 = - \sum_{\langle i,j \rangle} J_{ij} \delta_{\sigma_i, \sigma_j}$$

Trong đó, tương tác trao đổi $J_{ij}$ phân tách thành $J_k$ (tương tác giữa hai phân tử cùng mặt phẳng lớp) và $J_\perp$ (tương tác giữa hai phân tử ở hai mặt phẳng lân cận).

Để bổ sung tương tác tầm xa và năng lượng liên kết chất lỏng, luận án phát triển Hamiltonian mở rộng:

$$H_2 = - \sum_{\langle i,j \rangle} J_{ij} \delta_{\sigma_i, \sigma_j} + \sum_{\langle i,j \rangle} V_{LJ}(r_{ij})$$

với thế Lennard-Jones rút gọn $V_{LJ}(r) = 4\varepsilon [(\sigma/r)^{12} - (\sigma/r)^6]$ được chuẩn hóa thông qua tham số cường độ $V_0$.

Ba mệnh đề lý thuyết cốt lõi được luận chứng:

  • Mệnh đề P1: Tỷ số bất đối xứng tương tác $\kappa = J_\perp / J_k$ chi phối trực tiếp độ cứng vững trật tự smectic. Khi $J_\perp > 0$ (sắt từ), các lớp liên kết pha đồng nhất; khi $J_\perp < 0$ (phản sắt từ), hệ xuất hiện cấu trúc phản trật tự lớp xen kẽ.
  • Mệnh đề P2: Độ linh động phân tử tỉ lệ nghịch với nồng độ chiếm chỗ mạng $c = N_s / N_L$. Sự giảm mật độ cục bộ tại bề mặt tạo ra gradient năng lượng tự do kích hoạt sự mất trật tự phân lớp từ ngoài vào trong.
  • Mệnh đề P3: Thế năng tương tác cặp $V_0$ đóng vai trò áp suất liên kết nội tại. Khi $V_0 > 0$, hệ duy trì trạng thái ngưng tụ ổn định, chuyển hóa quá trình bay hơi phân tử thành quá trình tái sắp xếp nội tại, làm xuất hiện chuỗi chuyển pha nhiều bước sắc nét.

Luận án khẳng định chuyển pha Sm-I là chuyển pha loại 1 thông qua việc chứng minh hàm phân bố xác suất năng lượng $P_L(E)$ tại nhiệt độ chuyển pha có dạng hai đỉnh cực đại rõ rệt ngăn cách bởi một cực tiểu entropy, phản ánh sự tồn tại của nhiệt ẩn (latent heat) và sức căng bề mặt giữa pha smectic và pha isotropic.

                    HÀM PHÂN BỐ NĂNG LƯỢNG TẠI NHIỆT ĐỘ CHUYỂN PHA Tc
               PL(E)
                 ▲
                 │        Đỉnh Pha Smectic            Đỉnh Pha Isotropic
                 │           (Trật tự)                   (Hỗn loạn)
                 │             ┌───┐                        ┌───┐
                 │            │     │                      │     │
                 │           │       │                    │       │
                 │          │         │                  │         │
                 │         │           │    Rào cản ΔF   │           │
                 │        │             │   (Nhiệt ẩn)  │             │
                 │       │               └───┐        ┌──┘               │
                 │      │                     └──────┘                    │
                 │     │                                                   │
                 └─────┴───────────────────────────────────────────────────┴───►
                                             Năng lượng E

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

Khung phân tích của luận án tích hợp ba trụ cột:

  1. Mô hình spin Potts 6 trạng thái: Đại diện cho tính dị hướng quang học và trật tự định hướng của các phân tử calamitic.
  2. Cơ chế khuếch tán nút mạng trống: Đại diện cho tính chảy và độ nhớt của trạng thái lỏng.
  3. Thế tương tác đẩy - hút Lennard-Jones: Đại diện cho thể tích loại trừ của lõi cứng và lực hút van der Waals giữa các mạch alkyl.

Điều kiện biên của khung phân tích được xác định rõ: áp dụng cho tinh thể lỏng thermotropic đơn thành phần, cấu trúc phân tử đối xứng hai đầu (non-polar), bỏ qua tương tác chiral và điện trường ngoài.

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

Thiết kế nghiên cứu

Luận án theo đuổi lập trường bản thể luận thực tại khoa học (scientific realism) và phương pháp luận thực chứng tính toán (computational positivism) thuộc vật lý thống kê. Nghiên cứu sử dụng thiết kế đa quy mô: từ tham số vi mô ở cấp độ nút mạng quy đổi sang các đại lượng nhiệt động vĩ mô thông qua trung bình tập hợp canonical.

Mạng mô phỏng là mạng lập phương 3D kích thước $L \times L \times L$ ($L \in {12, 18, 24}$), tổng số nút mạng $N_L = L^3$ (từ 1.728 đến 13.824 vị trí). Số lượng phân tử $N_s$ được kiểm soát chính xác theo nồng độ $c \in {30%, 50%, 60%, 80%, 100%}$.

┌────────────────────────────────────────────────────────────────────────┐
│                        THIẾT KẾ MÔ PHỎNG HỆ THỐNG                      │
├───────────────────────────────────┬────────────────────────────────────┤
│ Kích thước mạng lập phương        │ $N_s = 12^3, 18^3, 24^3$           │
├───────────────────────────────────┼────────────────────────────────────┤
│ Dải nồng độ phân tử ($c$)         │ $30\%, 50\%, 60\%, 80\%, 100\%$   │
├───────────────────────────────────┼────────────────────────────────────┤
│ Tương tác nội lớp ($J_k$)         │ $1.0, 1.5, 3.0$                    │
├───────────────────────────────────┼────────────────────────────────────┤
│ Tương tác liên lớp ($J_\perp$)    │ $+0.5$ (Ferro), $-0.5$ (Anti-ferro)│
├───────────────────────────────────┼────────────────────────────────────┤
│ Cường độ thế Lennard-Jones ($V_0$)│ $0.2, 0.5, 1.0$                    │
└───────────────────────────────────┴────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình mô phỏng thực hiện qua hai thuật toán bổ trợ:

  1. Thuật toán Metropolis cải tiến:

    • Tại mỗi bước Monte Carlo Step (MCS), một phân tử được chọn ngẫu nhiên.
    • Phân tử thực hiện thử nghiệm hai bước động học độc lập: (i) Thay đổi trạng thái spin $\sigma_i \to \sigma_i'$ với xác suất $W = \min(1, \exp(-\Delta E / k_B T))$; (ii) Nhảy sang một trong các vị trí nút mạng trống lân cận gần nhất (NN) với xác suất tương tự.
    • Cải tiến then chốt: Tách biệt và cân bằng tỷ lệ giữa bước quay định hướng và bước dịch chuyển không gian, loại bỏ hoàn toàn hiện tượng bẫy cấu hình cục bộ (metastable states) làm lệch điểm chuyển pha.
    • Thời gian hồi phục cân bằng: $10^5$ MCS, thời gian lấy mẫu thống kê: $10^6$ MCS.

  2. Kỹ thuật Wang-Landau mật độ trạng thái:

    • Tính toán trực tiếp hàm mật độ trạng thái $g(E)$ không phụ thuộc vào nhiệt độ.
    • Lược đồ cập nhật: Mỗi khi trạng thái có năng lượng $E$ được thăm dò, mật độ trạng thái được biến đổi $g(E) \leftarrow g(E) \times f$ và biểu đồ năng lượng $H(E) \leftarrow H(E) + 1$.
    • Hệ số hiệu chỉnh ban đầu $f_0 = e^1 \approx 2.71828$, giảm dần theo quy tắc $f_{i+1} = \sqrt{f_i}$ khi biểu đồ $H(E)$ đạt độ phẳng tiêu chuẩn ($> 80%$).
    • Tiêu chuẩn hội tụ nghiêm ngặt: Quá trình dừng khi $f_{\text{final}} = \exp(10^{-8}) \approx 1.00000001$.
                        SƠ ĐỒ THUẬT TOÁN WANG-LANDAU
                                     │
                                     ▼
                      ┌─────────────────────────────┐
                      │ Khởi tạo g(E)=1, H(E)=0,    │
                      │ f = exp(1)                  │
                      └──────────────┬──────────────┘
                                     │
                                     ▼
                      ┌─────────────────────────────┐
                      │ Thử nghiệm bước nhảy MC     │
                      │ E_old ──► E_new             │
                      └──────────────┬──────────────┘
                                     │
                                     ▼
                      ┌─────────────────────────────┐
                      │ Xác suất chuyển:            │
                      │ P = min(1, g(E_old)/g(E_new)│
                      └──────────────┬──────────────┘
                                     │
                                     ▼
                      ┌─────────────────────────────┐
                      │ Cập nhật:                   │
                      │ g(E) = g(E)*f               │
                      │ H(E) = H(E) + 1             │
                      └──────────────┬──────────────┘
                                     │
                                     ▼
                      ┌─────────────────────────────┐
             Không    │ Kiểm tra độ phẳng biểu đồ:   │
         ┌────────────┤ H(E) >= 80% <H(E)>          │
         │            └──────────────┬──────────────┘
         │                           │ Đúng
         │                           ▼
         │            ┌─────────────────────────────┐
         │            │ Giảm f = sqrt(f), H(E) = 0  │
         │            └──────────────┬──────────────┘
         │                           │
         │                           ▼
         │    Không   ┌─────────────────────────────┐
         └────────────┤ f < exp(10^-8) ?            │
                      └──────────────┬──────────────┘
                                     │ Đúng
                                     ▼
                      ┌─────────────────────────────┐
                      │ Kết thúc: Xuất g(E) chuẩn   │
                      └─────────────────────────────┘

Data và phân tích

Từ mật độ trạng thái $g(E)$, các hàm nhiệt động lực học vĩ mô được tái tạo liên tục ở mọi nhiệt độ thông qua phân bố canonical:

$$Z(T) = \sum_E g(E) e^{-E/k_B T}$$

$$U(T) = \langle E \rangle = \frac{1}{Z(T)} \sum_E E g(E) e^{-E/k_B T}$$

$$C_v(T) = \frac{\langle E^2 \rangle - \langle E \rangle^2}{k_B T^2}$$

Tham số trật tự định hướng $M$ và độ cảm thăng giáng trật tự $\chi$ được xác định theo:

$$M = \frac{q \max(N_1, N_2, \dots, N_q) - N_s}{(q-1)N_s}$$

$$\chi = \frac{N_s}{k_B T} \left( \langle M^2 \rangle - \langle M \rangle^2 \right)$$

Hệ số khuếch tán $D$ được trích xuất từ độ dịch chuyển bình phương trung bình (Mean Squared Displacement - MSD) theo hệ thức Einstein:

$$D = \lim_{t \to \infty} \frac{\langle |\vec{r}(t) - \vec{r}(0)|^2 \rangle}{6t}$$

Toàn bộ hệ thống thuật toán được lập trình song song hóa và thực thi trên Hệ thống tính toán hiệu năng cao HPC01 của Viện Hàn lâm Khoa học và Công nghệ Việt Nam.

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

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

  1. Xác lập bản chất chuyển pha loại 1 của Sm-I: Kết quả hàm phân bố năng lượng $P_L(E)$ thu được từ kỹ thuật Wang-Landau tại nhiệt độ chuyển pha của hệ nồng độ $c=100%$ ($T_c = 0.8790$), $c=80%$ ($T_1 = 0.8320$), và $c=60%$ ($T_1 = 0.6010$) đều hiển thị cấu trúc hai đỉnh cực đại rõ nét. Sự gián đoạn đột ngột của tham số trật tự $M$ từ giá trị xấp xỉ $1.0$ tụt về $0$ kèm theo đỉnh phân kỳ nhọn của nhiệt dung riêng $C_v$ là bằng chứng nhiệt động không thể bác bỏ về chuyển pha loại 1.
                      BẢNG THÔNG SỐ CHUYỂN PHA THEO NỒNG ĐỘ VÀ MÔ HÌNH
┌────────────────┬──────────────┬──────────────┬────────────────────────────────┐
│ Nồng độ ($c$)  │ Mô hình      │ $T_c$ / $T_i$│ Đặc điểm chuyển pha            │
├────────────────┼──────────────┼──────────────┼────────────────────────────────┤
│ $100\%$        │ Potts chuẩn  │ $0.8790$     │ 1 bước nhọn, loại 1 điển hình  │
├────────────────┼──────────────┼──────────────┼────────────────────────────────┤
│ $80\%$         │ Potts chuẩn  │ $0.8320$     │ 1 bước, dịch chuyển nhiệt độ   │
├────────────────┼──────────────┼──────────────┼────────────────────────────────┤
│ $60\%$         │ Potts chuẩn  │ $0.6010$     │ Bắt đầu xuất hiện thăng giáng  │
├────────────────┼──────────────┼──────────────┼────────────────────────────────┤
│ $50\%$         │ Potts chuẩn  │ $0.747, 0.949$│ Phân rã cấu trúc theo nhiệt độ │
├────────────────┼──────────────┼──────────────┼────────────────────────────────┤
│ $50\%$         │ Mở rộng ($V_0=1$) │ $0.7641, 0.8796$ │ 2 bước chuyển pha rõ rệt   │
├────────────────┼──────────────┼──────────────┼────────────────────────────────┤
│ $30\%$         │ Mở rộng ($V_0=1$) │ $T_1 \dots T_3$ │ 3 đỉnh $C_v$: $0.7253, 0.7916, 0.8791$ │
├────────────────┼──────────────┼──────────────┼────────────────────────────────┤
│ $30\%$         │ Mở rộng ($V_0=0.2$)│ $T_1 \dots T_6$ │ 6 bước: $0.5311 \to 0.7002$    │
└────────────────┴──────────────┴──────────────┴────────────────────────────────┘
  1. Mô phỏng thành công hiện tượng tan chảy theo từng lớp (Layer-by-layer Melting): Thông qua ảnh chụp trạng thái cấu hình spin ở các mức nhiệt độ khác nhau (tại $T=0.747$ ngay sau khi bắt đầu tan chảy và $T=0.949 > T_c$), luận án chứng minh: các lớp bề mặt bị mất trật tự định hướng và khuếch tán trước, trong khi các lớp lõi trung tâm vẫn duy trì trật tự smectic hoàn hảo. Đây là công trình mô phỏng vi mô đầu tiên tái hiện thành công hiện tượng Stoebe quan sát được trong thực nghiệm năm 1994.
                               CƠ CHẾ TAN CHẢY TỪNG LỚP
               ┌──────────────────────────────────────────────────┐
               │ [Lớp 1] ░ ░ ░ ░ ░ ░ ░ (Tan chảy trước ở bề mặt)  │
               ├──────────────────────────────────────────────────┤
               │ [Lớp 2] | | | | | | | (Bắt đầu mất trật tự)      │
               ├──────────────────────────────────────────────────┤
               │ [Lớp 3] | | | | | | | (Lõi duy trì trật tự cao)  │
               ├──────────────────────────────────────────────────┤
               │ [Lớp 4] | | | | | | | (Lõi duy trì trật tự cao)  │
               ├──────────────────────────────────────────────────┤
               │ [Lớp 5] | | | | | | | (Bắt đầu mất trật tự)      │
               ├──────────────────────────────────────────────────┤
               │ [Lớp 6] ░ ░ ░ ░ ░ ░ ░ (Tan chảy trước ở bề mặt)  │
               └──────────────────────────────────────────────────┘

  1. Khám phá chuỗi chuyển pha đa bước ở nồng độ thấp ($c=30%-50%$): Trong mô hình Potts mở rộng với $V_0 = 1.0$, hệ tại $c=50%$ xuất hiện 2 đỉnh chuyển pha riêng biệt tại $T_1 = 0.7641$ và $T_2 = 0.8796$. Đặc biệt, tại $c=30%$, hệ xuất hiện 3 đỉnh chuyển pha tại $T_1 = 0.7253$, $T_2 = 0.7916$, $T_3 = 0.8791$. Khi giảm tương tác $V_0 = 0.2$, đường cong nhiệt dung $C_v$ ghi nhận tới 6 đỉnh phân kỳ cực đại tại $T_1 = 0.531146, T_2 = 0.604178, T_3 = 0.635270, T_4 = 0.656481, T_5 = 0.681227, T_6 = 0.700268$. Điều này phản ánh sự tan rã độc lập của từng cặp lớp phân tử.

  2. Tác động định hình cấu trúc của tính phân cực tương tác liên lớp ($J_\perp$): Khi tương tác liên lớp là sắt từ ($J_\perp = +0.5$), tham số trật tự $M$ suy giảm theo các bước chuyển pha rõ rệt. Khi chuyển sang tương tác phản sắt từ ($J_\perp = -0.5$), hệ hình thành cấu trúc smectic phản song song, làm thay đổi hoàn toàn đường cong cảm từ $\chi$ và động học khuếch tán phân tử $D(T)$.

  3. Vai trò ổn định hóa của thế Lennard-Jones ($V_0$): Trong mô hình Potts thuần túy, ở nồng độ loãng ($c \le 30%$), các phân tử bị bay hơi hoàn toàn thành pha khí trước khi kịp quan sát chuyển pha lỏng. Việc tích hợp $V_0 \in [0.2, 1.0]$ tạo ra một giếng thế liên kết, duy trì trạng thái lỏng cô đặc cục bộ và cho phép quan sát đầy đủ động học chuyển pha Sm-I.

Implications đa chiều

  • Về mặt lý thuyết: Kết quả luận án giải quyết tranh cãi kéo dài giữa mô hình mầm lệch mạng (dislocation loop) và biến dạng mặt khum (meniscus), chứng minh rằng sự thăng giáng mật độ phân tử và mất trật tự định hướng cục bộ tại lớp bề mặt là động lực nội tại kích hoạt quá trình bóc tách từng lớp.
  • Về mặt phương pháp: Khẳng định tính ưu việt tuyệt đối của kỹ thuật Wang-Landau so với thuật toán Metropolis truyền thống trong việc xử lý các rào cản năng lượng tự do sâu của chuyển pha loại 1 trong vật liệu mềm.
  • Về mặt ứng dụng thực tiễn: Các thông số nhiệt độ chuyển pha và điều kiện bền vững của cấu trúc lớp smectic cung cấp cơ sở dữ liệu định lượng cho việc thiết kế màn hình tinh thể lỏng sắt điện (Ferroelectric LCDs - FLCDs) với tốc độ đáp ứng micro-giây, cảm biến nhiệt độ màng mỏng chính xác cao ($18^\circ\text{C} - 34^\circ\text{C}$), và các bộ điều biến quang học (optical switches) dựa trên hiệu ứng lưỡng chiết.

Limitations và Future Research

Luận án thừa nhận một số giới hạn nghiên cứu:

  1. Không gian định hướng rời rạc: Việc giới hạn vector định hướng trong 6 trạng thái rời rạc ($q=6$) của mô hình Potts chưa phản ánh hoàn toàn tính liên tục trong không gian $O(3)$ của vector director $\vec{n}$.
  2. Mạng tinh thể cứng bức thiết kế: Các phân tử bị ràng buộc chuyển động trên các nút mạng lập phương cố định, chưa xét đến biến dạng co giãn liên tục của mạng (off-lattice continuum).
  3. Bỏ qua tương tác điện trường ngoài: Chưa tích hợp số hạng tương tác giữa mômen lưỡng cực phân tử với điện trường ngoài $\vec{E}$, vốn là yếu tố vận hành cốt lõi trong các linh kiện hiển thị thực tế.

Chương trình nghiên cứu tiếp theo được hoạch định:

  • Phát triển mô hình liên tục ngoài nút mạng (Off-lattice Monte Carlo) kết hợp thế Gay-Berne và thuật toán Wang-Landau song song GPU.
  • Mở rộng nghiên cứu sang các biến thể smectic nghiêng (Smectic C và Chiral Smectic C*), khảo sát hiện tượng phân cực điện tự phát (spontaneous ferroelectric polarization).
  • Mô phỏng động học chuyển pha trong không gian giam giữ nano (nanopores) và bề mặt có xử lý định hướng (anchoring substrates).

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

Luận án tạo ra tác động học thuật sâu sắc trong lĩnh vực vật lý chất ngưng tụ và vật lý tính toán, mở ra phương pháp tiếp cận mới để khảo sát các hệ vật chất mềm phức tạp. Công trình cung cấp giải pháp mô hình hóa có chi phí tính toán tối ưu hơn hàng chục lần so với phương pháp động lực học phân tử (Molecular Dynamics) toàn nguyên tử mà vẫn nắm bắt trọn vẹn bản chất nhiệt động học của chuyển pha.

Về mặt công nghiệp, kết quả nghiên cứu hỗ trợ trực tiếp cho các phòng R&D vật liệu quang điện tử trong việc tối ưu hóa tỷ lệ pha trộn mesogen nhằm mở rộng dải nhiệt độ hoạt động của thiết bị hiển thị và cảm biến quang nhiệt.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận khung phương pháp luận mô phỏng Monte Carlo nâng cao (Metropolis cải tiến và Wang-Landau) kèm mã nguồn giải thuật tối ưu trên siêu máy tính.
  • Các nhà vật lý lý thuyết và cơ học thống kê: Sở hữu bộ dữ liệu kiểm chứng chuẩn xác về phân bố năng lượng $P_L(E)$ và tính chất chuyển pha loại 1 trên mạng Potts linh động.
  • Kỹ sư R&D công nghệ hiển thị và quang học: Ứng dụng các quy luật ổn định nhiệt của cấu trúc phân lớp smectic để tổng hợp các hợp chất tinh thể lỏng có độ bền pha vượt trội.

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ì?

Đóng góp độc đáo nhất là việc mở rộng thành công mô hình Potts 6 trạng thái từ mạng cố định sang mạng linh động có chỗ trống kết hợp thế năng Lennard-Jones ($V_0$). Mô hình này dung hòa được lý thuyết spin thống kê với tính linh động khuếch tán của chất lỏng, giải thích chính xác hiện tượng chuyển pha loại 1 và sự sụp đổ trật tự từng lớp mà lý thuyết Landau-de Gennes liên tục không thể mô tả chi tiết ở cấp độ vi mô.

2. Cải tiến phương pháp luận có gì vượt trội so với các công trình quốc tế tiền nhiệm?

So với nghiên cứu của GS. H. Diep et al. (2020) chỉ dùng Metropolis đơn thuần dễ bị bẫy ở trạng thái giả cân bằng, luận án đã: (i) Cải tiến quy trình cập nhật hai bước động học độc lập (quay spin và nhảy vị trí); (ii) Ứng dụng kỹ thuật Wang-Landau để tính trực tiếp hàm mật độ trạng thái $g(E)$ với chuẩn hội tụ $f = \exp(10^{-8})$, cho phép tính toán chính xác hàm nhiệt dung $C_v$ và đường phân bố năng lượng $P_L(E)$ mà không gặp phải hiện tượng chậm tới hạn (critical slowing down).

3. Phát hiện nào gây bất ngờ nhất về mặt số liệu?

Phát hiện bất ngờ nhất là sự xuất hiện của cấu trúc đa đỉnh phân kỳ trên đường nhiệt dung riêng $C_v$ tại các nồng độ loãng $c=30%$ và $c=50%$ trong mô hình Potts mở rộng ($V_0=0.2$ ghi nhận tới 6 đỉnh chuyển pha từ $T_1=0.531146$ đến $T_6=0.700268$). Điều này chứng minh rằng ở mật độ thấp, quá trình tan chảy không diễn ra đồng loạt mà phân rã thành một chuỗi vi chuyển pha tuần tự của từng lớp riêng biệt.

4. Luận án có cung cấp đầy đủ giao thức để tái lặp nghiên cứu (Replication Protocol)?

Có. Luận án mô tả chi tiết:

  • Cấu hình mạng tinh thể lập phương $L = 12, 18, 24$;
  • Bộ tham số tương tác $J_k = 1.0, 1.5, 3.0$; $J_\perp = \pm 0.5$; $V_0 = 0.2, 0.5, 1.0$;
  • Tỷ lệ nồng độ $c = 30% - 100%$;
  • Tiêu chuẩn phẳng biểu đồ $H(E) \ge 80%$ và ngưỡng dừng hệ số $f_{\text{final}} = \exp(10^{-8})$;
  • Quy chuẩn tính toán các đại lượng $M, \chi, C_v, D$ trên hệ thống HPC01 - VAST.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình hướng tới: (1) Tích hợp trường ngoài (điện trường, từ trường) vào mô hình Potts linh động; (2) Chuyển đổi mô hình từ mạng rời rạc sang mô hình liên tục không gian thực sử dụng thế Gay-Berne sửa đổi; (3) Mở rộng nghiên cứu sang màng tự do kích thước lớn để mô phỏng chính xác số mũ thực nghiệm $N \propto t^{-0.4}$ của Stoebe; (4) Khảo sát các pha lyotropic và tương tác màng sinh học phospholipid.

Kết luận

  1. Luận án đã phát triển thành công mô hình vi mô toàn diện – mô hình Potts linh động 6 trạng thái kết hợp tương tác dị hướng và thế năng Lennard-Jones – mô tả chính xác cấu trúc phân lớp và trật tự định hướng của tinh thể lỏng smectic.
  2. Thiết lập thành công quy trình mô phỏng chính xác cao nhờ kết hợp thuật toán Metropolis cải tiến và kỹ thuật Wang-Landau mật độ trạng thái, vượt qua rào cản năng lượng tự do và hiện tượng chậm tới hạn.
  3. Chứng minh dứt khoát bản chất chuyển pha smectic - isotropic là chuyển pha loại 1 thông qua cấu trúc phân bố xác suất hai đỉnh $P_L(E)$, bước nhảy gián đoạn của tham số trật tự $M$ và đỉnh nhọn phân kỳ của nhiệt dung $C_v$.
  4. Lần đầu tiên tái hiện bằng mô phỏng số cơ chế động học tan chảy từng lớp (layer-by-layer melting) từ bề mặt vào lõi, cung cấp bằng chứng vi mô giải thích các số liệu thực nghiệm kinh điển của Stoebe et al. (1994).
  5. Khám phá hiện tượng phân rã đa bước (lên đến 6 đỉnh nhiệt độ chuyển pha tại $c=30%, V_0=0.2$) và làm rõ vai trò quyết định của thế năng $V_0$ trong việc ổn định hóa trạng thái lỏng cô đặc.
  6. Mở ra một hướng nghiên cứu mới trong vật lý tính toán chất mềm, cung cấp công cụ mô hình hóa hiệu năng cao cho các bài toán chuyển pha phức tạp trong vật lý và công nghệ vật liệu tương lai.