Tổng quan luận án

Bản luận án tiến sĩ vật lý với tiêu đề "Điều khiển độ căng của phân tử ADN trong dung môi phi tuyến bằng kìm quang học" được thực hiện bởi nghiên cứu sinh Thái Đình Trung, chuyên ngành Quang học (mã số: 62.09), dưới sự hướng dẫn khoa học của PGS. Hồ Quang Quý và TS. Đoàn Hoài Sơn tại Trường Đại học Vinh vào năm 2017.

Tính cấp thiết và khoảng trống nghiên cứu

Kìm quang học (bẫy quang học) kể từ khi được Arthur Ashkin phát minh vào năm 1970 đã trở thành công cụ quan trọng trong vật lý thực nghiệm và lý sinh phân tử để giam giữ, thao tác vi hạt điện môi và các đại phân tử sinh học. Trong nghiên cứu phân tử Acid DeoxyriboNucleic (ADN), việc kéo căng và gỡ xoắn chuỗi đơn phân tử đóng vai trò trung tâm nhằm tìm hiểu các cơ chế sao chép (replication), phiên mã (transcription) cũng như tương tác giữa protein và ADN. Vì phân tử ADN có đường kính rất nhỏ (khoảng 2 nm), các thí nghiệm không thể chiếu bẫy trực tiếp mà phải liên kết một đầu ADN với mặt đế thủy tinh và đầu kia gắn với vi hạt điện môi polystyrene có vai trò như một móc kéo.

Vấn đề tồn tại trong các kỹ thuật thực nghiệm hiện hành là việc điều khiển vị trí vi hạt và duy trì độ căng của chuỗi ADN chủ yếu dựa trên các cơ cấu điện - cơ, cụ thể là dịch chuyển vị trí vết hội tụ của chùm laser hoặc dịch chuyển bàn mẫu bằng gốm áp điện (piezoelectric stage). Cơ cấu cơ học này gặp hạn chế lớn về độ chính xác khi theo dõi và điều khiển ở thang đo nanomet, đồng thời có độ trễ nhất định khi cố gắng ổn định độ căng mong muốn của phân tử. Mặt khác, các hàm lực đàn hồi bán thực nghiệm theo mô hình chuỗi liên tục (Worm-Like Chain - WLC) cổ điển khi đưa vào phương trình động lực học Langevin bộc lộ các điểm kỳ dị toán học không phù hợp với thực tế vật lý (lực đàn hồi tiến tới vô cùng khi độ giãn bằng chiều dài cực đại, lực khác không khi độ giãn nhỏ hơn chiều dài bền). Do đó, yêu cầu đặt ra là phải hiệu chỉnh hàm lực đàn hồi cho phù hợp với bản chất vật lý của ADN và đề xuất một phương pháp điều khiển vị trí hạt cũng như độ căng ADN bằng hiệu ứng quang học thuần túy (toàn quang) thông qua môi trường phi tuyến Kerr thay cho chuyển động cơ học.

Mục tiêu nghiên cứu

Trên cơ sở nguyên lý hoạt động của kìm quang học và tính chất cơ học đàn hồi của chuỗi phân tử ADN, luận án giải quyết phương trình Langevin tổng quát dưới sự cạnh tranh của các lực (quang lực, lực đàn hồi, lực ma sát nhớt Stokes và lực dao động nhiệt Brown), khảo sát quá trình động học và các điều kiện ổn định vi hạt liên kết với ADN, từ đó đề xuất mô hình kìm quang học điều khiển độ căng của chuỗi ADN trong không gian ba chiều bằng cách thay đổi công suất laser trong môi trường dung môi phi tuyến Kerr.

Nội dung nghiên cứu

Để thực hiện mục tiêu trên, luận án xác định bốn nội dung nghiên cứu cụ thể:

  1. Hiệu chỉnh biến số của biểu thức lực đàn hồi của phân tử ADN dạng chuỗi con sâu (WLC) trong phương trình Langevin tổng quát nhằm mô tả chính xác động lực học của hạt liên kết với phân tử ADN dưới tác động của quang lực.
  2. Thiết lập biểu thức giải tích của quang lực tác dụng lên vi hạt điện môi trong không gian ba chiều đối với trường hợp vi hạt nhúng trong môi trường tuyến tính và môi trường phi tuyến Kerr ở chế độ Rayleigh.
  3. Sử dụng phương pháp mô phỏng số khảo sát quá trình động lực học của vi hạt trong bẫy quang học; phân tích sự phụ thuộc của động học và trạng thái ổn định của hạt vào các tham số của chùm laser (cường độ đỉnh, bán kính thắt chùm).
  4. Đề xuất mô hình kìm quang học phi tuyến và khảo sát quá trình điều khiển vị trí cân bằng của vi hạt cùng độ căng của phân tử ADN trong không gian ba chiều bằng việc thay đổi công suất laser.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu:
    1. Kìm quang học sử dụng một chùm tia Gaussian trong môi trường tuyến tính và kìm quang học phi tuyến (vi hạt nhúng trong môi trường chất lưu phi tuyến Kerr).
    2. Lực đàn hồi của các chuỗi phân tử ADN (dạng plasmid và thực khuẩn thể phage $\lambda$) và độ căng tương ứng.
    3. Động lực học của vi hạt điện môi gắn với phân tử ADN chịu tác dụng đồng thời của quang lực, lực đàn hồi và các lực môi trường.
  • Phạm vi nghiên cứu: Luận án giới hạn trong phạm vi xây dựng mô hình lý thuyết, phân tích giải tích và mô phỏng số quá trình động lực học; không bao gồm chế tạo thực nghiệm vật lý trong phòng thí nghiệm.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước đây

Luận án điểm lại hệ thống các công trình nghiên cứu nền tảng về kìm quang học và cơ học phân tử ADN:

  • Nguyên lý bẫy quang học: A. Ashkin và cộng sự (1970, 1986) đặt nền móng cho kỹ thuật kìm quang học thông qua việc khảo sát áp suất bức xạ, chứng minh chùm laser hội tụ mạnh có thể bẫy và điều khiển vi hạt điện môi trong không gian ba chiều nhờ lực gradient vượt qua lực tán xạ và áp suất bức xạ nhiệt. Các cơ sở tính toán quang lực trong chế độ Rayleigh ($\lambda \gg 2a$) và chế độ Mie/quang hình ($\lambda \le 2a$) được xác lập bởi Harada, Askin (1996) và Rohrbach (2002).
  • Ứng dụng kìm quang học trong cơ học sinh học: Kỹ thuật bẫy hạt nhanh chóng được triển khai để nghiên cứu hệ thống sinh học: xác định liên kết đuôi vi khuẩn, đo lực sinh ra do sự vận động của đơn phân tử protein (Finer và cs., 1994), đo độ biến dạng màng tế bào hồng cầu (Hénon và cs., 1999; Sleep và cs., 1999).
  • Kéo căng và xác định đàn hồi ADN:
    • Wang và cộng sự (1997) thực hiện thí nghiệm kéo căng đơn phân tử ADN bằng kìm quang học kết hợp liên kết vi hạt polystyrene và neo RNA polymerase trên đế kính.
    • Bustamante và cộng sự (2000) tái lập thí nghiệm, khẳng định khả năng dùng kìm quang học khảo sát quá trình kéo căng, gỡ xoắn ADN nhằm phục vụ phân tích quá trình tái bản và phiên mã.
    • Strick và cộng sự (2000) tổng quát hóa quá trình xoắn và kéo căng của đơn phân tử ADN, xác định suất căng đạt khoảng 1000 pN và phân tích ảnh hưởng của chuyển động Brown.
    • Singer và cộng sự (2001) khảo sát động lực học của vi cầu bằng cách thay đổi quang lực ba chiều qua dao động tắt dần để suy ra hệ số cản.
    • Soni và cộng sự (2002) phát triển hệ thống hai bẫy laser kết hợp bàn gốm áp điện ba chiều để điều khiển vị trí hai đầu phân tử ADN với độ chính xác cỡ nanomet.
    • Hegner và cộng sự (2002) thiết lập đường cong đặc trưng giữa lực đàn hồi và độ biến dạng cho các mẫu ADN khác nhau.
    • Huang và cộng sự (2015) đề xuất bẫy quang học plasmonic vượt giới hạn nhiễu xạ; Jadhav và cộng sự (2016) sử dụng bẫy quang học định lượng tương tác giữa protein và ADN.
  • Mô hình đàn hồi của polyme sinh học:
    • Marko và Siggia (1995) thiết lập công thức nội suy cho mô hình chuỗi con sâu (WLC) mô tả đường cong lực - độ giãn đối với chuỗi ADN sợi đôi xoắn kép.
    • Baumann và cs. (1997, 2000), Bustamante và cs. (1994, 2000), Mack (2002), Sharma, Hamdi phát triển các biểu thức gần đúng và mô hình WLC mở rộng có tính đến biến dạng đàn hồi của khung phân tử ở dải lực cao (10 - 50 pN).

Khoảng trống luận án lựa chọn giải quyết

Tổng quan tài liệu chỉ ra hai khoảng trống quan trọng:

  1. Về mặt cơ cấu điều khiển: Tất cả các hệ kìm quang học đo đạc trước đây đều dùng cơ chế dịch chuyển cơ - điện thụ động (thay đổi vị trí bẫy bằng cách di chuyển cụm thấu kính hoặc đế mẫu bằng gốm áp điện). Cơ chế này tạo ra rung động cơ học, độ trễ và khó duy trì độ ổn định liên tục ở thang đo dưới nanomet trong quá trình đo đạc động lực học.
  2. Về mặt giải tích mô hình: Biểu thức lực đàn hồi WLC của Marko - Siggia và Bustamante xem biến số độ biến dạng $x$ là chiều dài tổng tính từ gốc tọa độ ($0$ đến $L$). Khi $x \to L$, hàm số cho lực đàn hồi $F \to \infty$, điều này mâu thuẫn với thực tế là ADN có suất căng giới hạn $S \approx 1000\text{ pN}$ và sẽ đứt tại độ giãn xác định. Đồng thời khi $x \le L_b$ (chiều dài bền), biểu thức cũ vẫn tính ra lực dương dù phân tử chưa chịu biến dạng kéo. Sự không tương thích này làm cho các hàm giải tích cũ không thể tích hợp trực tiếp vào phương trình Langevin để giải số động lực học vi hạt trong kìm quang học. Luận án giải quyết triệt để hai vấn đề này bằng việc hiệu chỉnh biến số hàm WLC và đưa ra mô hình bẫy quang học trong dung môi phi tuyến Kerr.

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở lý thuyết và khung phân tích

Luận án kết hợp lý thuyết đàn hồi cơ học đại phân tử, lý thuyết điện từ cổ điển và quang học phi tuyến:

  • Cấu trúc phân tử ADN: ADN sợi đôi dạng xoắn kép (mô hình Watson - Crick dạng B) có đường kính 2 nm, chiều cao mỗi chu kỳ xoắn gồm 10 cặp nucleotide là 3,4 nm (khoảng cách giữa hai bậc thang là 0,34 nm). Các liên kết hóa trị photphodieste dọc theo mạch đơn đảm bảo tính ổn định và các liên kết hydro nối các cặp base bổ sung ($A-T$ bằng 2 liên kết, $G-C$ bằng 3 liên kết). Luận án khảo sát hai dạng cấu trúc hình học: phân tử dạng vòng (plasmid DNA) và phân tử chuỗi kép thẳng (phage $\lambda$ DNA).
  • Lý thuyết biến dạng đàn hồi thanh mỏng: Dựa trên định luật nén/kéo tuyến tính, ứng suất $\sigma = Y \frac{\Delta L}{L_0}$, trong đó $Y$ là suất Young, $S = \pi r^2 Y$ là suất căng (mô đun kéo). Đối với ADN sợi đôi đường kính 2 nm, $Y \approx 3 \cdot 10^8\text{ Pa}$, suy ra $S \approx 1000\text{ pN}$. Độ cứng kéo của thanh được tính qua năng lượng tích trữ nhiệt động: $k_{rod} = \frac{S}{L_0}$. Năng lượng biến dạng uốn được xác định bởi: $$E_{bend} = \frac{B}{2} \int \left(\frac{1}{R}\right)^2 ds$$ với hằng số uốn $B = YI = \frac{\pi}{4} Y r^4$. Chiều dài bền $L_p$ (hoặc ký hiệu $L_b$) liên hệ với độ cứng uốn qua biểu thức: $$L_p = \frac{B}{k_B T}$$ Đối với ADN sợi đôi xoắn kép trong điều kiện sinh lý, $L_p \approx 50\text{ nm}$.
  • Mô hình chuỗi liên tục (WLC): Mô hình WLC xem đại phân tử như một thanh đàn hồi liên tục chịu dao động nhiệt. Ở vùng lực kéo nhỏ ($F < k_B T / L_p \approx 0,1\text{ pN}$), lực bị chi phối bởi entropy: $$F \approx \frac{3 k_B T}{2 L_p} \frac{x}{L_0}$$ Ở vùng lực lớn ($1 - 10\text{ pN}$), lực tỷ lệ nghịch với bình phương khoảng cách tới giới hạn: $$F \approx \frac{k_B T}{4 L_p (1 - x/L_0)^2}$$ Công thức nội suy Marko - Siggia liên kết hai vùng: $$F(x) = \frac{k_B T}{L_p} \left[ \frac{x}{L_0} + \frac{1}{4(1 - x/L_0)^2} - \frac{1}{4} \right]$$
  • Lý thuyết quang lực trong chế độ Rayleigh: Kích thước vi hạt điện môi polystyrene thỏa mãn điều kiện Rayleigh khi bán kính hạt $a$ nhỏ hơn nhiều so với bước sóng laser ($\lambda \gg 2a$). Khi đó hạt được coi như một lưỡng cực điện điểm phân cực tuyến tính chịu tác dụng của lực Lorentz: $$\vec{F} = q (\vec{E} + \vec{v} \times \vec{B})$$ Sử dụng hệ phương trình Maxwell và lấy trung bình theo chu kỳ dao động quang học ($10^{14}\text{ Hz}$), quang lực phân tách thành hai thành phần chính:
    1. Lực tán xạ ($F_{scat}$): Tác động đẩy hạt dọc theo trục truyền sóng ($z$), tỷ lệ thuận với cường độ chùm tia $I(r)$ và tiết diện tán xạ $\sigma_s = \frac{8\pi}{3} k^4 a^6 \left(\frac{m^2-1}{m^2+2}\right)^2$: $$\vec{F}_{scat}(r) = \frac{n_b}{c} \sigma_s I(r) \hat{z}$$
    2. Lực gradient ($F_{grad}$): Kéo vi hạt về phía có mật độ năng lượng điện trường cực đại (vùng thắt chùm laser), tỷ lệ với độ phân cực $\alpha = 4\pi \epsilon_0 n_f^2 a^3 \left(\frac{m^2-1}{m^2+2}\right)$: $$F_{grad,z}(z) = \frac{2\pi n_b a^3}{c} \left(\frac{m^2-1}{m^2+2}\right) \frac{\partial I}{\partial z}$$ $$F_{grad,\rho}(\rho) = \frac{2\pi n_f a^3}{c} \left(\frac{m^2-1}{m^2+2}\right) \frac{\partial I}{\partial \rho}$$ trong đó $m = n_b / n_f$ là tỷ số chiết suất giữa hạt điện môi ($n_b$) và chất lưu ($n_f$).
  • Môi trường phi tuyến Kerr: Khi chất lưu xung quanh hạt có tính phi tuyến quang học Kerr, chiết suất phụ thuộc cường độ laser theo hệ thức $n = n_f + n_{fnl} I$. Sự tự hội tụ hoặc tự tán xạ làm thay đổi cấu trúc chùm tia Gaussian, làm dịch chuyển vị trí thắt chùm $w_{0,1}$ và tiêu cự phi tuyến $f_{nl}$ theo công suất laser truyền vào.

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

Luận án áp dụng phương pháp nghiên cứu lý thuyết kết hợp mô hình hóa giải tích và mô phỏng số:

  • Thiết lập phương trình vi phân ngẫu nhiên Langevin mô tả chuyển động của vi hạt có khối lượng $M$: $$M \frac{dv}{dt} = \vec{f}{vis} + \vec{f}{Brown} + \vec{F}{optic} + \vec{F}{el}$$ trong đó $\vec{f}{vis} = -\gamma \vec{v} = -6\pi\eta a \vec{v}$ là lực ma sát nhớt theo định luật Stokes ($\eta$ là độ nhớt động học); $\vec{f}{Brown} = \sqrt{2D} h(t)$ là lực va đập nhiệt ngẫu nhiên với hệ số khuếch tán Einstein - Smoluchowski $D = \frac{k_B T}{\gamma}$ và $h(t) \in [-1, 1]$ là biến ngẫu nhiên Gauss; $\vec{F}{optic} = \vec{F}{grad} + \vec{F}{scat}$ là tổng quang lực; $\vec{F}{el}$ là lực đàn hồi của chuỗi ADN.
  • Giải số phương trình Langevin bằng thuật toán lặp sai phân bước thời gian để xác định quỹ đạo vị trí $\vec{r}(t)$ và vận tốc $\vec{v}(t)$ của hạt từ vị trí ban đầu đến khi rơi vào tâm bẫy quang học và đạt trạng thái cân bằng.

Nguồn tư liệu và dữ liệu

Dữ liệu thông số vật lý của các mẫu ADN được tác giả kế thừa từ các công trình đo đạc thực nghiệm của Bustamante và cs. (1994, 2000), Smith và cs. (1996), Baumann và cs. (1997, 2000) đối với ADN plasmid và phage $\lambda$ trong các dung dịch điện ly có nồng độ ion muối khác nhau ($\text{Na}^+$, $\text{Spd}^{3+}$, $\text{CoHex}$).

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

Chương 1: Kìm quang học cho phân tử ADN

Chương 1 hệ thống hóa cơ sở sinh học và cơ học của phân tử ADN, chỉ ra các khiếm khuyết trong các biểu thức WLC hiện hữu và đề xuất công thức hiệu chỉnh.

Tác giả phân tích mô hình WLC cổ điển dạng Marko - Siggia và dạng mở rộng: $$F(x) = \frac{k_B T}{L_b} \left[ \frac{x}{L} + \frac{1}{4(1 - x/L)^2} - \frac{1}{4} \right]$$ và chỉ ra ba mâu thuẫn vật lý:

  1. Khi độ biến dạng $x \to L$, mẫu số tiến về $0$ dẫn đến $F \to \infty$. Thực tế, khi bị kéo căng tới chiều dài cực đại, lực đàn hồi của ADN chỉ đạt giá trị hữu hạn cỡ vài chục pN trước khi đạt tới suất căng đứt gãy ($S \approx 1000\text{ pN}$).
  2. Khi $0 < x \le L_b$, phân tử chưa bị kéo giãn khỏi chiều dài bền cân bằng nhưng biểu thức cũ vẫn cho giá trị $F > 0$. Lực đàn hồi đúng bản chất phải triệt tiêu ($F = 0$) trong toàn bộ miền $x \le L_b$.
  3. Biến số $x$ trong các biểu thức cũ biểu diễn khoảng cách tuyệt đối từ gốc tọa độ thay vì độ giãn biến dạng thực tế $\Delta x$, do đó không thể gán tương thích với tọa độ vị trí của hạt trong phương trình vi phân Langevin.

Để khắc phục, tác giả đề xuất thay thế biến số $x/L$ bằng đại lượng dịch chuyển tương đối $(x - L_b)/L$. Biểu thức lực đàn hồi hiệu chỉnh có dạng: $$F(x) = \frac{k_B T}{L_b} \left[ \frac{x - L_b}{L} + \frac{1}{4\left(1 - \frac{x - L_b}{L}\right)^2} - \frac{1}{4} \right]$$ Và biểu thức WLC mở rộng có xét đến suất căng $S$: $$\frac{F(x) L_b}{k_B T} = \frac{1}{4\left(1 - \frac{x - L_b}{L} + \frac{F(x)}{S}\right)^2} - \frac{1}{4} + \frac{x - L_b}{L} - \frac{F(x)}{S}$$

Hàm hiệu chỉnh thỏa mãn các điều kiện biên vật lý:

  • $F(x) = 0$ khi $x \le L_b$.
  • $F(x) > 0$ khi $x > L_b$.
  • Tại độ giãn bằng chiều dài cực đại ($x = L$), giá trị lực đạt hữu hạn: $$F(L) = \frac{k_B T}{L_b} \left[ 1 - \frac{L_b}{L} + \frac{1}{4(L_b / L)^2} - \frac{1}{4} \right]$$

Luận án dẫn giải các thông số thực nghiệm của ADN plasmid trong các dung môi đệm ion khác nhau nhằm phục vụ kiểm chứng mô phỏng:

TT Thành phần đệm dung dịch Chiều dài bền $L_b$ (nm) Hằng số lực $S$ (pN) Chiều dài tổng $L$ (nm)
1 10 mM $\text{Na}^+$ 47,0 1008 1343
2 10 mM $\text{Na}^+$ và 10 $\mu\text{M}$ $\text{Spd}^{3+}$ 40,4 1166 1353
3 10 mM $\text{Na}^+$ và 100 $\mu\text{M}$ $\text{Spd}^{3+}$ 48,0 1202 1313
4 10 mM $\text{Na}^+$ và 200 $\mu\text{M}$ $\text{Spd}^{3+}$ 40,0 1253 1318
5 10 mM $\text{Na}^+$ và 300 $\mu\text{M}$ $\text{Spd}^{3+}$ 38,3 - -
6 10 mM $\text{Na}^+$ và 100 $\mu\text{M}$ $\text{CoHex}$ 44,6 - -
7 186 mM $\text{Na}^+$ và 100 $\mu\text{M}$ $\text{Spd}^{3+}$ 48,8 - -
8 186 mM $\text{Na}^+$ và 100 $\mu\text{M}$ $\text{CoHex}$ 32,0 - -

Kết quả so sánh trên đồ thị giữa công thức cũ và công thức hiệu chỉnh cho thấy: với mẫu số 1, khi độ biến dạng tiến sát $L = 1343\text{ nm}$, công thức cũ dự báo lực tăng vọt lên tới $10^5\text{ pN}$ (vượt xa suất căng đứt), trong khi công thức hiệu chỉnh cho lực đàn hồi đạt khoảng 16 pN. Đối với toàn bộ các mẫu plasmid, lực đàn hồi cực đại dao động trong khoảng 16 - 40 pN, hoàn toàn phù hợp với dữ liệu đo đạc thực nghiệm kéo căng của Bustamante.

Chương 2: Quá trình động lực học của vi hạt liên kết với phân tử ADN trong kìm quang học

Chương 2 tập trung khảo sát động học bẫy hạt trong bẫy quang học một chùm tia tuyến tính bằng phương trình Langevin.

Tác giả phân tích chi tiết cơ chế truyền xung lượng giữa photon và hạt điện môi:

  • Hiện tượng phản xạ bề mặt tạo nên biến thiên xung lượng $\Delta \vec{p} = 2 \hbar \vec{k} \cos\theta_i$, hình thành lực tán xạ hướng dọc theo chiều truyền ánh sáng.
  • Hiện tượng khúc xạ qua hạt làm lệch phương truyền tia sáng, tạo ra biến thiên xung lượng ngược hướng làm xuất hiện lực gradient hướng về vị trí có mật độ photon cao nhất (tâm bẫy).
  • Trong chế độ Rayleigh ($\lambda \gg 2a$), lực gradient ngang, gradient dọc và lực tán xạ được tính toán chính xác theo phân bố chùm tia Gaussian có bán kính thắt chùm $w_0$ và cường độ đỉnh $I_0$.

Mô hình hệ vật lý gồm: một vi hạt polystyrene có bán kính $a$, khối lượng $M$, liên kết với một đầu chuỗi ADN; đầu còn lại của ADN neo chặt tại bề mặt tấm thủy tinh. Toàn bộ hệ nhúng trong dung dịch chất lưu. Chùm laser Gaussian chiếu hội tụ tạo bẫy quang học giữ hạt.

Phương trình động lực học Langevin được thiết lập chi tiết dưới dạng hệ phương trình vi phân từng tọa độ: $$M \frac{d^2 \rho}{dt^2} + 6\pi\eta a \frac{d\rho}{dt} + F_{optic,\rho}(\rho, z) + F_{el,\rho}(\rho, z) = \sqrt{2D} h_\rho(t)$$ $$M \frac{d^2 z}{dt^2} + 6\pi\eta a \frac{dz}{dt} + F_{optic,z}(\rho, z) + F_{el,z}(\rho, z) = \sqrt{2D} h_z(t)$$

Thông qua giải số phương trình Langevin, chương 2 rút ra các kết quả động học quan trọng:

  1. Đặc trưng quỹ đạo vị trí - thời gian: Vi hạt dao động tắt dần dưới tác dụng của lực cản nhớt và nhanh chóng rơi vào vị trí cân bằng ổn định trong bẫy sau khoảng thời gian cỡ mili giây.
  2. Ảnh hưởng của cường độ đỉnh chùm laser ($I_0$): Khi tăng $I_0$, quang lực gradient tăng tuyến tính, làm tăng độ cứng của bẫy ($k_{trap}$), kéo hạt về gần tâm thắt chùm hơn, đồng thời làm tăng độ giãn và lực căng của chuỗi ADN.
  3. Ảnh hưởng của bán kính thắt chùm ($w_0$): Chùm tia hội tụ càng hẹp ($w_0$ nhỏ) thì gradient điện trường càng lớn, dẫn đến lực bẫy càng mạnh và vị trí cân bằng của vi hạt càng ổn định trước các dao động nhiệt ngẫu nhiên của lực Brown.

Chương 3: Điều khiển độ căng của phân tử ADN bằng kìm quang học Kerr

Chương 3 giải quyết bài toán cốt lõi của luận án: thiết kế cơ chế điều khiển toàn quang học vị trí cân bằng của vi hạt và độ căng của phân tử ADN bằng cách sử dụng môi trường phi tuyến Kerr.

Nguyên lý hoạt động dựa trên hiệu ứng chiết suất phi tuyến: khi chùm laser có công suất cao truyền qua môi trường chất lưu có hệ số phi tuyến $n_{fnl}$, môi trường đóng vai trò như một thấu kính mỏng phi tuyến có tiêu cự $f_{nl}$ phụ thuộc vào cường độ ánh sáng: $$f_{nl} \approx \frac{w_0^2}{4 n_{fnl} d I_0}$$ trong đó $d$ là bề dày của môi trường chất lưu Kerr.

Sự xuất hiện của thấu kính phi tuyến làm tái phân bố trường laser Gaussian:

  • Vị trí thắt chùm dịch chuyển dọc trục $z$ một đoạn $\Delta z(P)$ tỷ lệ với công suất laser $P$.
  • Bán kính thắt chùm tái tạo $w_{0,1}$ biến đổi tương ứng, làm thay đổi gradient cường độ ánh sáng trong không gian ba chiều.

Tác giả đề xuất mô hình kìm quang học sử dụng hai nguồn laser kết hợp:

  1. Nguồn laser kiểm soát (laser mạnh): Điều chỉnh công suất để làm biến đổi chiết suất môi trường chất lưu Kerr, từ đó điều khiển vị trí thắt chùm (tâm bẫy) dịch chuyển dọc trục hoặc mở rộng/thu hẹp phân bố trường bẫy.
  2. Nguồn laser bẫy (laser yếu): Tạo quang lực bẫy giữ vi hạt tại vị trí cân bằng được thiết lập bởi trường laser tổng.

Dưới sự cạnh tranh giữa quang lực phi tuyến điều khiển và lực đàn hồi WLC hiệu chỉnh của chuỗi ADN, vị trí cân bằng của vi hạt dịch chuyển có quy luật:

  • Điều khiển hướng tâm ($\rho$): Thay đổi công suất laser làm thay đổi gradient ngang, cho phép điều chỉnh tọa độ bán kính của hạt vi mô với độ chính xác dưới nanomet.
  • Điều khiển dọc trục ($z$): Bằng cách tăng hoặc giảm công suất laser mạnh, tâm bẫy dịch chuyển dọc theo trục quang học, kéo dãn hoặc chùng chuỗi ADN neo trên đế kính.

Kết quả mô phỏng khẳng định: độ căng của chuỗi phân tử ADN (từ trạng thái chùng ban đầu $F = 0$ khi $x \le L_b$ đến trạng thái căng cực đại $16 - 40\text{ pN}$) hoàn toàn có thể khống chế và điều chỉnh liên tục thông qua việc điều biến công suất laser mà không cần sử dụng bất kỳ chuyển động cơ học nào của bàn mẫu.

Kết quả và những đóng góp mới

Đóng góp mới về mặt lý luận

  1. Hiệu chỉnh giải tích hàm lực đàn hồi WLC: Đã phát hiện và khắc phục thành công các mâu thuẫn toán học và vật lý của mô hình Marko - Siggia và Bustamante. Việc đưa biến số $(x - L_b)/L$ vào biểu thức giúp hàm lực đàn hồi phản ánh đúng bản chất vật lý: lực bằng không khi độ dài chưa vượt qua chiều dài bền ($x \le L_b$), triệt tiêu điểm dị thường vô hạn tại chiều dài cực đại, cho giá trị lực hữu hạn phù hợp với thực nghiệm kéo đứt, và đồng nhất tọa độ biến số với phương trình động lực học Langevin.
  2. Thiết lập mô hình động lực học hoàn chỉnh cho vi hạt liên kết ADN: Xây dựng thành công phương trình vi phân ngẫu nhiên Langevin trong không gian ba chiều tích hợp đồng thời lực ma sát nhớt Stokes, lực ngẫu nhiên Brown, quang lực Rayleigh đầy đủ (lực tán xạ, lực gradient ngang và gradient dọc) và lực đàn hồi ADN hiệu chỉnh; phân tích định lượng ảnh hưởng của các thông số quang học ($I_0, w_0$) đến độ ổn định của bẫy.
  3. Lý thuyết kìm quang học phi tuyến Kerr: Xác lập cơ sở lý thuyết về sự cạnh tranh giữa lực quang học phi tuyến và lực đàn hồi đại phân tử sinh học; mô hình hóa sự tái phân bố chùm laser Gaussian qua thấu kính phi tuyến Kerr và cơ chế dịch chuyển tâm bẫy quang học bằng công suất laser.

Đóng góp mới về mặt thực tiễn và giải pháp đề xuất

  1. Đề xuất cấu hình kìm quang học toàn quang: Đưa ra mô hình kìm quang học hai chùm tia trong dung môi phi tuyến Kerr cho phép điều khiển vị trí vi hạt và độ căng của phân tử ADN hoàn toàn bằng quang học.
  2. Khắc phục nhược điểm của hệ cơ - điện: Giải pháp điều khiển bằng công suất laser loại bỏ hoàn toàn nhu cầu dịch chuyển cơ học của hệ thấu kính hoặc bàn áp điện, loại trừ rung động cơ học và độ trễ, nâng cao độ ổn định và độ chính xác của phép đo cơ học phân tử ở thang đo nanomet và dải lực piconewton.

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

Hạn chế của luận án

  • Toàn bộ kết quả trong luận án được xây dựng trên cơ sở tính toán lý thuyết, mô hình hóa giải tích và mô phỏng số; chưa có điều kiện bố trí hệ thực nghiệm vật lý để đo đạc kiểm chứng trực tiếp trên mẫu ADN thực tế trong dung môi Kerr.
  • Mô hình quang lực chỉ giới hạn khảo sát trong chế độ tán xạ Rayleigh ($\lambda \gg 2a$), áp dụng cho các vi hạt điện môi có kích thước rất nhỏ so với bước sóng; chưa mở rộng cho các hạt kích thước micro thuộc chế độ Mie hoặc quang hình.
  • Mô hình dung môi giả định chất lưu phi tuyến Kerr đồng nhất, đẳng hướng và không xét đến các hiệu ứng thứ cấp như hấp thụ phi tuyến (hấp thụ hai photon) hay gradient nhiệt độ do laser công suất cao gây ra trong chất lưu.

Hướng nghiên cứu tiếp

  • Mở rộng lý thuyết tính toán quang lực và phương trình Langevin sang chế độ Mie ($a \approx \lambda$) và chế độ quang học sóng đầy đủ nhằm áp dụng cho các loại hạt bẫy kích thước lớn hơn.
  • Thiết kế và triển khai mô hình thực nghiệm thực tế sử dụng các chất lưu có hệ số phi tuyến Kerr cao và tương thích sinh học để kiểm chứng các đường cong dịch chuyển vị trí hạt theo công suất laser.
  • Khảo sát các cấu hình chùm tia laser phức tạp hơn như chùm tia Bessel không nhiễu xạ hoặc chùm tia xoáy quang học mang mômen động lượng quỹ đạo để điều khiển đồng thời cả độ căng và độ xoắn của phân tử ADN.

Giá trị tham khảo

Bản luận án mang lại giá trị tham khảo chuyên môn rõ rệt đối với các nhóm đối tượng học thuật và ứng dụng:

  • Nghiên cứu sinh và học viên cao học ngành Quang học, Vật lý laser và Vật lý lý thuyết: Luận án cung cấp hệ thống công thức giải tích chuẩn xác về quang lực trong chế độ Rayleigh, phương pháp mô hình hóa chùm laser Gaussian qua môi trường phi tuyến Kerr và thuật toán giải số phương trình vi phân ngẫu nhiên Langevin.
  • Các nhà nghiên cứu trong lĩnh vực Lý sinh học phân tử và Kỹ thuật y sinh:
    • Phần hiệu chỉnh hàm lực đàn hồi WLC ở Chương 1 là tài liệu tham khảo trực tiếp cho các nghiên cứu mô phỏng cơ học biến dạng của ADN, ARN và các chuỗi polypeptide/protein.
    • Số liệu tổng hợp về chiều dài bền, chiều dài cực đại và suất căng của các mẫu plasmid ADN trong các điều kiện đệm ion muối khác nhau là dữ liệu chuẩn hữu ích cho việc so sánh thực nghiệm.
  • Kỹ sư thiết kế thiết bị quang - cơ điện tử và bẫy vi hạt: Đề xuất về cơ chế dịch chuyển tâm bẫy toàn quang bằng hiệu ứng Kerr ở Chương 3 gợi mở hướng phát triển các thế hệ bẫy quang học không chuyển động cơ học (non-mechanical optical manipulation), ứng dụng trong vi lưu (microfluidics) và phòng thí nghiệm trên chip (lab-on-a-chip).

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

1. Tại sao không thể sử dụng trực tiếp công thức WLC của Marko - Siggia trong phương trình Langevin mô phỏng kìm quang học?

Công thức WLC của Marko - Siggia sử dụng biến số $x$ là tọa độ tuyệt đối tính cả chiều dài bền $L_b$, dẫn đến lực đàn hồi khác không ngay cả khi phân tử chưa bị kéo giãn ($x \le L_b$). Đặc biệt, khi phân tử bị kéo giãn tới chiều dài cực đại ($x \to L$), lực đàn hồi tiến tới vô cùng ($\infty$). Trong phương trình vi phân ngẫu nhiên Langevin, điểm kỳ dị vô hạn này khiến các thuật toán giải số bị phân kỳ và không thể phản ánh trạng thái cân bằng lực thực tế giữa quang lực và lực đàn hồi khi hạt điện môi dao động quanh tâm bẫy.

2. Biến đổi cốt lõi trong công thức hiệu chỉnh lực đàn hồi của tác giả là gì và có ý nghĩa vật lý như thế nào?

Biến đổi cốt lõi là thay thế tỷ số độ dài tuyệt đối $x/L$ bằng tỷ số biến dạng tương đối $(x - L_b)/L$. Về mặt vật lý, biến đổi này đảm bảo rằng lực đàn hồi chỉ xuất hiện khi chiều dài phân tử vượt qua chiều dài bền ($x > L_b$); khi $x \le L_b$, lực đàn hồi triệt tiêu hoàn toàn ($F = 0$). Ngoài ra, tại độ dài cực đại ($x = L$), lực đàn hồi đạt giá trị xác định (khoảng 16 - 40 pN đối với plasmid và phage $\lambda$), nhỏ hơn nhiều so với suất căng đứt gãy ($S \approx 1000\text{ pN}$), hoàn toàn phù hợp với thực nghiệm kéo căng phân tử.

3. Tại sao thí nghiệm kìm quang học cần liên kết phân tử ADN với vi hạt điện môi polystyrene thay vì bẫy trực tiếp ADN?

Phân tử ADN có kích thước mặt cắt ngang rất bé (đường kính chỉ khoảng 2 nm), quá nhỏ so với bước sóng ánh sáng laser ($0,5 - 1\text{ }\mu\text{m}$) nên độ phân cực quang học không đủ lớn để tạo ra quang lực gradient thắng được chuyển động nhiệt Brown. Hạt vi cầu polystyrene có đường kính lớn hơn nhiều, chiết suất thích hợp ($n_b \approx 1,57 - 1,59$) và không hấp thụ ánh sáng trong vùng laser sử dụng, đóng vai trò như một "tay cầm" hoặc "móc kéo" nhận xung lượng từ trường laser để truyền lực căng đến chuỗi ADN.

4. Cơ chế điều khiển độ căng của phân tử ADN trong dung môi phi tuyến Kerr hoạt động như thế nào mà không cần hệ cơ học?

Khi chiếu chùm laser có công suất biến thiên vào chất lưu có tính phi tuyến Kerr ($n = n_f + n_{fnl} I$), chùm tia tự làm biến đổi chiết suất môi trường và biến khối chất lưu thành một thấu kính mỏng phi tuyến có tiêu cự $f_{nl}$ thay đổi theo công suất ánh sáng. Thấu kính này làm dịch chuyển vị trí thắt chùm $w_{0,1}$ (tâm bẫy quang học) dọc trục truyền sóng và thay đổi gradient cường độ. Bằng cách điều chỉnh công suất của nguồn laser kiểm soát, tâm bẫy quang học dịch chuyển trong không gian, làm thay đổi khoảng cách giữa hạt và điểm neo của ADN trên đế kính, qua đó kiểm soát chính xác độ căng của phân tử hoàn toàn bằng ánh sáng.

5. Những tham số quang học nào của chùm laser ảnh hưởng trực tiếp đến trạng thái cân bằng của vi hạt trong bẫy?

Các tham số quang học chính gồm cường độ đỉnh của chùm tia ($I_0$) và bán kính thắt chùm ($w_0$). Cường độ đỉnh $I_0$ quyết định độ lớn của quang lực gradient và độ cứng bẫy ($k_{trap}$); khi $I_0$ tăng, lực bẫy tăng giúp giữ hạt ổn định hơn trước dao động nhiệt Brown. Bán kính thắt chùm $w_0$ quyết định độ dốc phân bố năng lượng; chùm tia hội tụ càng hẹp ($w_0$ nhỏ) thì gradient điện trường càng lớn, tạo quang lực gradient ngang và dọc càng mạnh, giúp định vị vi hạt chính xác tại tâm bẫy.

Kết luận

Luận án tiến sĩ của tác giả Thái Đình Trung đã giải quyết hoàn chỉnh bài toán lý thuyết về điều khiển độ căng của đơn phân tử ADN bằng kìm quang học trong dung môi phi tuyến Kerr. Bằng việc hiệu chỉnh thành công biến số trong mô hình WLC, tác giả đã khắc phục triệt để điểm kỳ dị toán học của các công thức đàn hồi cổ điển, cho phép tích hợp chính xác lực đàn hồi vào phương trình ngẫu nhiên Langevin để mô phỏng động lực học hạt bẫy trong không gian ba chiều. Đóng góp nổi bật của công trình là đề xuất giải pháp điều khiển vị trí vi hạt và độ căng chuỗi ADN bằng biến thiên công suất laser trong môi trường phi tuyến, mở ra hướng tiếp cận toàn quang học giúp nâng cao độ chính xác trong các thao tác vi mô trên đại phân tử sinh học.