Tổng quan luận án

Luận án Tiến sĩ Vật lý chuyên ngành Quang học (mã số: 62.09) của nghiên cứu sinh Hoàng Đình Hải, được thực hiện dưới sự hướng dẫn khoa học của PGS. Hồ Quang Quý tại Viện Khoa học và Công nghệ Quân sự - Bộ Quốc phòng (Hà Nội, 2014), mang tên: "Nghiên cứu ảnh hưởng của một số thông số lên kìm quang học sử dụng hai chùm xung Gauss ngược chiều".

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

Bẫy quang học (optical trap) hay kìm quang học (optical tweezers) từ khi được Arthur Ashkin và cộng sự phát hiện và thử nghiệm thành công (1970, 1986) đã trở thành một công cụ vi thao tác quan trọng trong vật lý, hóa lý và y sinh học, cho phép bẫy giữ và điều khiển các vi hạt, nguyên tử, chuỗi ADN hoặc tế bào sống mà không cần tiếp xúc cơ học.

Trong thực tế kỹ thuật laser hiện đại, hầu hết các hệ laser phát ra từ buồng cộng hưởng gương cầu đều có phân bố không gian dạng chùm Gauss và thường hoạt động ở chế độ điều biến xung ngắn (xung nano giây, pico giây hoặc femto giây). Tuy nhiên, các nghiên cứu lý thuyết kinh điển về kìm quang học phần lớn chỉ khảo sát trong chế độ sóng phẳng hoặc chùm laser phát liên tục (CW). Một số nghiên cứu tính toán cho kìm quang học sử dụng một chùm xung Gauss (tiêu biểu như nhóm nghiên cứu của Zhao) đã chỉ ra rằng: việc dùng một chùm xung laser đơn lẻ rất khó tạo ra bẫy ổn định trên trục truyền lan, do thành phần lực tán xạ (scattering force) biến thiên phụ thuộc vào thời gian và độ rộng xung, dẫn đến mất cân bằng lực dọc tại vị trí bẫy.

Để khắc phục nhược điểm của kìm đơn chùm, mô hình kìm quang học sử dụng hai chùm xung Gauss truyền ngược chiều đã được đề xuất nhằm triệt tiêu lực tán xạ dọc trục. Mặc dù một số công trình ban đầu đã bước đầu tính toán quang lực và vẽ các mặt phẳng pha, việc phân tích định lượng chi tiết về sự chuyển động động học của vi hạt trong chất lưu dưới tác động đồng thời của quang lực xung và lực chuyển động Brown (Brownian motion), cũng như ảnh hưởng của các tham số quang - cơ - nhiệt đến độ ổn định của vi hạt trong bẫy vẫn còn bị bỏ ngỏ. Trong thực nghiệm, vi hạt hay tế bào sống trong bẫy không đứng yên tuyệt đối mà dao động xung quanh tâm bẫy do độ cứng (stiffness) hữu hạn của bẫy và sự tán xạ nhiệt từ môi trường. Do đó, việc nghiên cứu có hệ thống động học của vi hạt và xác lập các điều kiện ổn định không gian - thời gian là yêu cầu cấp thiết để định hướng cho thiết kế thực nghiệm.

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

Luận án xác định hai mục tiêu nghiên cứu cụ thể:

  1. Đưa ra được các luận cứ có tính khoa học, xây dựng các điều kiện để có thể ổn định được các vi hạt có kích thước cỡ nanomet bằng bẫy quang học sử dụng hai chùm xung Gauss ngược chiều.
  2. Nghiên cứu được các mối quan hệ tương quan giữa các tham số quang, tham số cơ, tham số nhiệt tham gia trong quá trình ổn định hạt điện môi trong kìm quang học sử dụng hai chùm xung Gauss ngược chiều.

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

  • Đối tượng nghiên cứu: Kìm quang học sử dụng hai chùm xung Gauss ngược chiều. Đối tượng bị bẫy là vi hạt điện môi có kích thước nanomet (vi cầu thủy tinh) lơ lửng trong môi trường chất lưu (nước và một số chất lưu khác).
  • Phạm vi nghiên cứu:
    • Chế độ tán xạ Rayleigh: xét kích thước bán kính vi hạt $a$ nhỏ hơn nhiều so với bước sóng ánh sáng ($\lambda = 1{,}064\ \mu\text{m}$, $a \ll \lambda$, khảo sát các bán kính hạt từ $8\text{ nm}$ đến $18\text{ nm}$).
    • Cấu hình chùm tia: hai chùm laser Gauss đối xứng truyền ngược chiều dọc theo trục $z$, có phân cực vuông góc với nhau để triệt tiêu hiện tượng giao thoa.
    • Môi trường nghiên cứu: vi hạt nhúng trong chất lưu, chịu lực cản nhớt Stokes và lực thăng giáng nhiệt ngẫu nhiên (lực Brown).
    • Phương pháp giải quyết: mô hình hóa lý thuyết và mô phỏng số trên môi trường Matlab.

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 được điểm luận

Lịch sử nghiên cứu áp lực ánh sáng và bẫy quang học được tác giả hệ thống hóa qua các mốc thời gian và công trình tiêu biểu:

  • Khởi nguồn cơ sở lý thuyết áp lực bức xạ: Johannes Kepler (1619) giải thích hiện tượng đuôi sao chổi hướng ra xa Mặt Trời do áp lực ánh sáng. James Clerk Maxwell (1871) chứng minh bằng lý thuyết điện từ rằng sóng ánh sáng mang xung lượng và tạo ra áp lực lên bề mặt tiếp xúc. Peter Lebedev (1901) và Ernest Fox Nichols (1901) độc lập thực hiện các thí nghiệm đo đạc áp lực ánh sáng (light pressure) lên vật thể.
  • Sự ra đời của bẫy quang học: Arthur Ashkin (1970, 1971) sử dụng chùm laser hội tụ để cân bằng áp lực bức xạ với trọng lực, giữ hạt điện môi kích thước $20\ \mu\text{m}$ trong chất lỏng và chất khí; đồng thời đề xuất cấu hình bẫy 3 chiều dùng hai chùm tia truyền ngược chiều. Năm 1986, Ashkin, Steven Chu và cộng sự chế tạo thành công bẫy quang học đơn chùm sử dụng vật kính có khẩu độ số lớn ($NA$) để giữ các vi hạt có đường kính từ $25\text{ nm}$ đến $10\ \mu\text{m}$ trong nước.
  • Phương pháp tính quang lực theo các chế độ kích thước:
    • Chế độ Rayleigh ($a \ll \lambda$): xem vi hạt như một lưỡng cực điện cảm ứng chịu lực Lorentz trung bình thời gian, bao gồm lực gradient hướng về cực đại cường độ và lực tán xạ đẩy theo chiều truyền sóng.
    • Chế độ quang hình ($a \gg \lambda$): tính toán dựa trên định luật khúc xạ, phản xạ và bảo toàn động lượng tia sáng.
    • Chế độ Mie ($a \approx \lambda$): áp dụng lý thuyết tán xạ điện từ trường đầy đủ của Mie.
    • Các mở rộng khác: O. Stout sử dụng phương pháp véctơ tính quang lực cho chùm có phân bố cường độ bất kỳ; D. Bonessi và cộng sự kết hợp ba chế độ tính cho hạt có hình dạng bất kỳ; Cui (2008) và Jian (2009) tính toán quang lực tác dụng lên các cấu trúc kích thước nano. H. Kim nghiên cứu sự phụ thuộc của hiệu suất bẫy vào tỉ số bán kính chùm trên khẩu độ số của hệ quang, khảo sát ảnh hưởng của độ nhớt môi trường, bán kính hạt và công suất laser.
  • Nghiên cứu bẫy quang học sử dụng xung laser:
    • Zhao và cộng sự công bố các tính toán quang lực cho kìm dùng một chùm xung Gauss, khảo sát phân bố lực gradient ngang và lực dọc trên các mặt phẳng pha $(\rho, z)$, $(\rho, t)$ và $(z, t)$. Nhóm nghiên cứu chỉ ra rằng: chuyển động Brown bị triệt tiêu tạm thời khi bẫy bằng xung ngắn; nhưng giữa hai xung nếu tần số lặp nhỏ hơn $100\text{ Hz}$, hạt sẽ chuyển động Brown tự do trong thể tích cỡ $10^{-12}\text{ m}^3$. Quan trọng nhất, lực tán xạ biến đổi theo thời gian và độ rộng xung nên không thể giữ hạt ổn định trên trục truyền tia.
    • M. Kawano và cộng sự (2008) đề xuất sử dụng hai chùm laser ngược chiều để bẫy quang học vi hạt.
    • Mai Văn Lưu, Hồ Quang Quý / H. Lou (2009) tính toán quang lực của hai chùm xung Gauss ngược chiều, chứng minh lực gradient và lực tán xạ có tính chất đối ngẫu qua tâm bẫy, loại bỏ được sự phụ thuộc đơn hướng của lực tán xạ vào độ rộng xung.

Khoảng trống nghiên cứu luận án lựa chọn giải quyết

Mặc dù Kawano và Lou đã chỉ ra tính đối xứng của quang lực trong hệ hai chùm ngược chiều, các công trình này mới chỉ dừng lại ở việc khảo sát phân bố lực trong các mặt phẳng pha lý thuyết. Khoảng trống then chốt mà luận án tập trung giải quyết là:

  • Chưa có mô hình động học hoàn chỉnh mô tả chuyển động của vi hạt nano trong chất lưu dưới tác động đồng thời của hệ lực quang học điều biến theo thời gian và lực ngẫu nhiên nhiệt Brown.
  • Chưa lượng hóa được ảnh hưởng của các tham số quang (năng lượng xung $U$, độ rộng xung $\tau$, bán kính thắt chùm $W_0$, khoảng cách giữa hai mặt thắt $d$, độ trễ xung $\delta T$, tần số lặp xung $f$) và tham số cơ - nhiệt (bán kính hạt $a$, độ nhớt môi trường $\eta$) lên quỹ đạo và tốc độ dịch chuyển của vi hạt.
  • Thiếu một định nghĩa định lượng rõ ràng về "vùng ổn định không gian - thời gian" (thời gian ổn định và kích thước không gian dao động) của hạt bị bẫy trong hệ xung laser ngược chiều.

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

Khung lý thuyết và khái niệm sử dụng

Luận án xây dựng trên nền tảng lý thuyết tương tác giữa trường điện từ laser với hạt điện môi và lý thuyết chuyển động Brown trong chất lưu:

  1. Trường điện từ của chùm xung Gauss: Cường độ điện trường của chùm xung Gauss truyền dọc trục $z$ được mô tả trong gần đúng cận trục: $$E(\rho, z, t) = \hat{x} E_0 \left{ \frac{W_0}{\sqrt{W_0^2 + 4z^2/k^2}} \right} \exp\left{ -i[kz - \omega t] \right} \exp\left[ -\frac{\rho^2}{W^2(z)} \right] \exp\left[ -2\left(\frac{t - z/c}{\tau}\right)^2 \right]$$ Trong cấu hình hai chùm ngược chiều, chùm trái có tâm thắt tại $z = -d/2$, chùm phải có tâm thắt tại $z = d/2$ ($d$ là khoảng cách giữa hai mặt thắt). Hai chùm có phương phân cực vuông góc nhau, do đó cường độ tổng $I(\rho, z, t, d)$ là tổng đại số cường độ của từng chùm riêng biệt: $$I(\rho, z, t, d) = I_L(\rho, z, t, d) + I_R(\rho, z, t, d)$$

  2. Quang lực trong chế độ tán xạ Rayleigh ($a \ll \lambda$): Hạt điện môi được xem như một lưỡng cực cảm ứng có độ phân cực: $$\alpha = 4\pi \varepsilon_0 n_2^2 a^3 \left( \frac{m^2 - 1}{m^2 + 2} \right)$$ với $n_2$ là chiết suất của môi trường chất lưu, $n_1$ là chiết suất hạt điện môi, $m = n_1 / n_2$ là tỉ số chiết suất tương đối, $a$ là bán kính hạt.

    Quang lực trung bình tác dụng lên hạt gồm:

    • Lực gradient (kéo hạt về vị trí có mật độ năng lượng cực đại): $$\vec{F}{grad} = \frac{1}{2} \alpha \nabla \langle E^2 \rangle = \frac{2\pi n_2 a^3}{c} \left( \frac{m^2 - 1}{m^2 + 2} \right) \nabla I$$ Lực gradient được phân tích thành thành phần hướng tâm $F{grad,\rho}$ theo tọa độ xuyên tâm $\rho$ và thành phần dọc trục $F_{grad,z}$.
    • Lực tán xạ (đẩy hạt theo chiều truyền sóng): $$\vec{F}_{scat} = \frac{8\pi n_2 k^4 a^6}{3c} \left( \frac{m^2 - 1}{m^2 + 2} \right)^2 I \vec{z} = \frac{128\pi^5 n_2 a^6}{3 c \lambda^4} \left( \frac{m^2 - 1}{m^2 + 2} \right)^2 I \vec{z}$$

  3. Phương trình động học Langevin: Chuyển động của vi hạt trong chất lưu chịu tác dụng đồng thời của lực quang học, lực cản ma sát nhớt của môi trường và lực dao động nhiệt ngẫu nhiên: $$m \frac{d^2\vec{r}}{dt^2} + \gamma \frac{d\vec{r}}{dt} = \vec{F}{opt}(\vec{r}, t) + \vec{F}{Brown}(t)$$ Trong đó:

    • $m$ là khối lượng hạt vi cầu.
    • $\gamma = 6\pi \eta a$ là hệ số ma sát nhớt theo định luật Stokes ($\eta$ là độ nhớt động lực học của chất lưu).
    • $\vec{F}{opt} = \vec{F}{grad} + \vec{F}_{scat}$ là tổng quang lực tác dụng lên hạt.
    • $\vec{F}{Brown}(t)$ là lực ngẫu nhiên Brown gây ra bởi va chạm nhiệt của các phân tử môi trường, thỏa mãn điều kiện trung bình thống kê $\langle \vec{F}{Brown}(t) \rangle = 0$ và tương quan thời gian $\langle \vec{F}{Brown}(t) \vec{F}{Brown}(t') \rangle = 2 k_B T \gamma \delta(t - t')$.

Phương pháp nghiên cứu và quy trình tính toán

Luận án kết hợp giải tích lý thuyết và mô phỏng số:

  • Phương pháp mô phỏng số: Thiết lập thuật toán giải phương trình vi phân ngẫu nhiên Langevin trên phần mềm tính toán khoa học Matlab.
  • Quy trình mô phỏng:
    • Khảo sát sự hội tụ của nghiệm số bằng cách chia nhỏ bước thời gian $\delta t = 6\tau/2000$, $6\tau/4000$, $6\tau/8000$, $6\tau/10000$.
    • Tính toán phân bố không gian và thời gian của cường độ sáng và quang lực trong các mặt phẳng pha $(\rho, z)$, $(\rho, t)$, $(z, t)$.
    • Mô phỏng quỹ đạo chuyển động ngẫu nhiên của vi hạt từ các vị trí xuất phát ban đầu khác nhau trong mặt phẳng tiêu bản, theo dõi quá trình hạt bị hút về tâm bẫy và dao động xung quanh tâm.
  • Tham số chuẩn sử dụng trong mô phỏng:
Đại lượng / Tham số Ký hiệu Giá trị danh định trong mô phỏng Đơn vị
Bước sóng chùm laser $\lambda$ $1{,}064$ $\mu\text{m}$
Chiết suất chất lưu (nước) $n_2$ $1{,}332$ Không thứ nguyên
Khẩu độ số của hệ vật kính $NA$ $1{,}0$ Không thứ nguyên
Bán kính thắt chùm laser $W_0$ $1{,}0$ (khảo sát $0{,}5 - 4{,}0$) $\mu\text{m}$
Bán độ rộng xung laser $\tau$ $1{,}0$ (khảo sát $0{,}5 - 2{,}0$) $\text{ps}$
Bán kính hạt điện môi (thủy tinh) $a$ $10$ (khảo sát $8 - 18$) $\text{nm}$
Năng lượng xung laser $U$ $0{,}9$ (khảo sát đến $5{,}0$) $\mu\text{J}$
Khoảng cách hai mặt thắt $d$ $10$ (khảo sát $0 - 20$) $\mu\text{m}$

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

Chương 1. Tổng quan về kìm quang học sử dụng hai chùm xung Gauss ngược chiều

Chương 1 thiết lập mô hình cấu hình quang học, biểu thức giải tích của cường độ điện trường tổng và các thành phần quang lực, đồng thời khảo sát sự biến thiên của chúng theo các tham số hình học và quang học:

  • Cấu hình quang hai chùm ngược chiều: Chùm laser được phân tách làm hai nhánh phân cực vuông góc, qua các gương lái tia và hội tụ đối xứng vào buồng mẫu. Mặt phẳng tiêu điểm của hai chùm được dịch chuyển cách nhau khoảng cách $d$ ($z = -d/2$ và $z = d/2$).
  • Phân bố cường độ sáng tổng:
    • Khảo sát biến thiên theo khoảng cách $d$ ($d = 15\ \mu\text{m}$, $10\ \mu\text{m}$, $5\ \mu\text{m}$, $0\ \mu\text{m}$): Khi $d = 15\ \mu\text{m}$ và $10\ \mu\text{m}$, phân bố cường độ tổng xuất hiện hai đỉnh rõ rệt đối xứng qua gốc tọa độ. Khi khoảng cách giảm xuống $d \le 5\ \mu\text{m}$, hai đỉnh này chập lại thành một đỉnh duy nhất có cường độ cực đại nằm ngay tại tâm bẫy ($z = 0, \rho = 0$).
    • Khảo sát theo bán kính thắt chùm $W_0$ ($W_0 = 2{,}0;\ 1{,}5;\ 1{,}0;\ 0{,}5\ \mu\text{m}$ với $d = 10\ \mu\text{m}$): Khi $W_0$ lớn, cường độ tổng nhỏ và phân bố tương đối đều. Khi $W_0$ giảm dần, cường độ cực đại tăng vọt từ $3500\ \mu\text{W}/\mu\text{m}^2$ (với $W_0 = 2\ \mu\text{m}$) lên đến $\sim 9000\ \mu\text{W}/\mu\text{m}^2$ (với $W_0 = 1\ \mu\text{m}$); tại $W_0 < 1\ \mu\text{m}$, dạng phân bố Gauss bị phá vỡ và tập trung cường độ rất cao tại hai đỉnh hẹp gần gốc tọa độ.
  • Phân bố quang lực dọc và ngang:
    • Quang lực dọc trục ($F_z$): Trên mặt phẳng pha $(z, t)$, xuất hiện hai vùng bẫy ổn định đối xứng nhau qua gốc tọa độ $z = 0$. Khi tăng bán kính thắt chùm $W_0$ từ $0{,}5\ \mu\text{m}$ lên $2\ \mu\text{m}$, kích thước vùng bẫy tăng lên nhưng giá trị cực đại của lực giảm mạnh (từ $F_{z,\max} \approx 1{,}032\ \text{pN}$ tại $W_0 = 0{,}5\ \mu\text{m}$ giảm xuống mức thấp hơn nhiều ở $W_0 = 2\ \mu\text{m}$).
    • Ảnh hưởng của độ rộng xung $\tau$: Khảo sát $\tau$ từ $0{,}5\text{ ps}$ đến $2\text{ ps}$ cho thấy kích thước vùng bẫy tăng theo độ rộng xung, nhưng độ lớn lực cực đại giảm tỉ lệ nghịch: tại $\tau = 0{,}5\text{ ps}$ lực cực đại đạt $8{,}5 \cdot 10^3\ \text{pN}$, trong khi tại $\tau = 1\text{ ps}$ lực cực đại giảm xuống còn $4{,}5 \cdot 10^3\ \text{pN}$.
    • Ảnh hưởng của khoảng cách $d$: Khi tăng $d$ ($5, 10, 15, 20\ \mu\text{m}$), hai vùng bẫy tách xa nhau, tạo ra một vùng ranh giới ở giữa mà tại đó quang lực triệt tiêu, khiến hạt chuyển động nhiệt tự do.
    • Quang lực ngang ($F_{grad,\rho}$): Tạo ra hố thế thế năng giam giữ hạt theo phương xuyên tâm $\rho$. Lực ngang cực đại đạt được tại thời điểm đỉnh xung ($t = 1\tau$) và giảm nhanh khi hạt đi ra xa trục quang học.

Chương 2. Quá trình động học của vi hạt

Chương 2 giải quyết bài toán động lực học của vi hạt dưới tác động đồng thời của quang lực và lực Brown trong môi trường chất lưu:

  • Kiểm tra tính hội tụ của mô hình số: Thực hiện mô phỏng chuyển động Brown thuần túy của vi hạt thủy tinh trong nước ở vị trí ban đầu $\rho_0 = 1\ \mu\text{m}$ và $\rho_0 = 0\ \mu\text{m}$ với các bước thời gian mô phỏng khác nhau: $\delta t = 6\tau/2000$, $6\tau/4000$, $6\tau/8000$, $6\tau/10000$. Kết quả cho thấy khi $\delta t \le 6\tau/8000$, quỹ đạo thống kê của hạt đạt trạng thái hội tụ ổn định, đảm bảo độ chính xác cho việc tích phân phương trình Langevin.
  • Sự cạnh tranh giữa quang lực và lực Brown:
    • Trong khoảng thời gian xung laser xuất hiện (thời gian cỡ vài pico giây), quang lực gradient chiếm ưu thế hoàn toàn so với lực ngẫu nhiên Brown. Lực này gia tốc và kéo hạt điện môi dịch chuyển rất nhanh về phía tâm bẫy.
    • Khi xung laser đi qua (khoảng thời gian nghỉ giữa các xung), quang lực biến mất, vi hạt chịu sự chi phối hoàn toàn của lực cản nhớt môi trường và lực Brown, bắt đầu khuếch tán ngẫu nhiên tự do.
  • Quá trình động học khi có quang lực: Mô phỏng vị trí của vi hạt trên mặt phẳng tiêu bản với các bộ tham số năng lượng khác nhau ($U = 0{,}9\ \mu\text{J}$, $U = 5\ \mu\text{J}$) cho thấy: nếu năng lượng xung đủ lớn, hạt sẽ bị giam giữ hiệu quả trong một bán kính dao động giới hạn quanh gốc tọa độ. Thời gian kéo hạt về tâm và quỹ đạo chuyển động phụ thuộc trực tiếp vào biên độ quang lực xung.

Chương 3. Ảnh hưởng của các tham số lên quá trình động học của vi hạt

Chương 3 đi sâu khảo sát định lượng sự phụ thuộc của quá trình chuyển dịch của vi hạt về tâm bẫy vào 4 tham số chính:

  1. Ảnh hưởng của vị trí ban đầu ($\rho_0$): Khảo sát các vị trí xuất phát ban đầu $\rho_0 = 1\ \mu\text{m}$, $2\ \mu\text{m}$. Hạt càng ở xa trục chùm tia thì lực gradient tác dụng ban đầu càng nhỏ (do phân bố Gauss suy giảm ở phần đuôi), dẫn đến thời gian cần thiết để hạt bị hút về tâm càng dài. Nếu hạt ở ngoài phạm vi tác dụng hữu hiệu của hố thế quang học, lực Brown sẽ làm hạt khuếch tán ra xa và thoát khỏi bẫy.
  2. Ảnh hưởng của năng lượng tổng ($U$): Khảo sát với các mức năng lượng $U = 0{,}9\ \mu\text{J}$ và $5\ \mu\text{J}$. Năng lượng xung tăng làm tăng trực tiếp độ dốc thế năng và độ lớn của lực gradient. Do đó, tốc độ kéo hạt về tâm bẫy tăng lên đáng kể, đồng thời biên độ dao động tản mạn do lực Brown quanh tâm bẫy bị nén lại rõ rệt.
  3. Ảnh hưởng của bán kính thắt chùm ($W_0$): Khảo sát các giá trị $W_0 = 1, 2, 3, 4\ \mu\text{m}$ (trong điều kiện $U = 0{,}9\ \mu\text{J}$, $\tau = 1\text{ ps}$, $\rho_0 = 1\ \mu\text{m}$, $\lambda = 1{,}064\ \mu\text{m}$, $a = 10\text{ nm}$, khoảng thời gian $t = 0 - 6\text{ ps}$):
    • Khi $W_0 = 1\ \mu\text{m}$, mật độ quang thông tại tâm rất lớn, hạt bị kéo mạnh về tâm và định vị chặt chẽ.
    • Khi tăng $W_0$ lên $2, 3, 4\ \mu\text{m}$, mật độ năng lượng giảm theo tỉ lệ nghịch với $W_0^2$, lực gradient giảm làm tốc độ chuyển dịch của hạt về tâm chậm đi; hạt dao động với biên độ rộng hơn trên mặt phẳng tiêu bản.
  4. Ảnh hưởng của kích thước hạt ($a$): Khảo sát các bán kính hạt $a = 8\text{ nm}$, $10\text{ nm}$, $14\text{ nm}$, $18\text{ nm}$ (với $U = 0{,}9\ \mu\text{J}$, $\tau = 1\text{ ps}$, $\rho_0 = 1\ \mu\text{m}$, $W_0 = 1\ \mu\text{m}$):
    • Do lực gradient tỷ lệ thuận với thể tích hạt ($a^3$) trong khi hệ số cản nhớt chỉ tỷ lệ với $a$, các hạt có kích thước lớn hơn ($14\text{ nm}$, $18\text{ nm}$) chịu quang lực vượt trội hơn hẳn so với lực ma sát và lực Brown, bị kéo về tâm nhanh hơn và dao động quanh tâm với biên độ nhỏ hơn hạt $a = 8\text{ nm}$.

Chương 4. Ảnh hưởng của các tham số lên vùng ổn định

Chương 4 là phần trọng tâm định lượng của luận án, đề xuất khái niệm mới và thiết lập các mối quan hệ chi phối vùng giam giữ ổn định:

  • Khái niệm "Vùng ổn định không gian - thời gian": Trong điều kiện chiếu xung lặp lại, vi hạt không đứng yên tại một điểm toán học cố định. Luận án định nghĩa vùng ổn định là một miền không gian hình trụ ("pillar") giới hạn trên mặt phẳng tiêu bản, đặc trưng bởi hai thông số:
    • Đường kính vùng ổn định ($\Delta \rho$): phạm vi không gian lớn nhất mà vi hạt dao động quanh tâm bẫy.
    • Thời gian ổn định ($\Delta t$): khoảng thời gian liên tục mà hạt được giam giữ an toàn bên trong đường kính $\Delta \rho$ mà không bị thoát ra ngoài dưới tác động của lực Brown.
       Mô hình định tính vùng ổn định không gian - thời gian
             +-----------------------------------------+
             |         Đường kính vùng ổn định:        |
             |               Δρ (hoặc ds)              |
             |    |<----------------------------->|    |
             |                                         |
Thời gian    |    . - . - . - . - . - . - . - . - .    |
bẫy giữ      |    |       Quỹ đạo hạt dao động    |    |
ổn định:     |    |       xung quanh tâm bẫy:     |    |
Δt           |    |       (Cạnh tranh giữa        |    |
             |    |        Quang lực xung và      |    |
             |    |        chuyển động Brown)     |    |
             |    . - . - . - . - . - . - . - . - .    |
             +-----------------------------------------+
  • Kết quả khảo sát ảnh hưởng của các tham số lên vùng ổn định:
    • Năng lượng xung laser ($U$): Khi tăng năng lượng xung, thời gian ổn định $\Delta t$ tăng tuyến tính, trong khi đường kính vùng ổn định $\Delta \rho$ thu hẹp lại, hạt được định vị chính xác hơn.
    • Bán kính thắt chùm ($W_0$): Khi $W_0$ tăng, gradient cường độ giảm, làm giảm thế năng bẫy; do đó thời gian ổn định $\Delta t$ suy giảm và đường kính $\Delta \rho$ tăng rộng ra.
    • Bán độ rộng xung ($\tau$): Độ rộng xung tăng từ $0{,}5\text{ ps}$ lên các giá trị lớn hơn làm kéo dài thời gian tác dụng của quang lực trong mỗi chu kỳ xung, làm thay đổi kích thước vùng ổn định tương ứng.
    • Tần số lặp xung ($f$): Khảo sát sự ổn định của vi hạt phụ thuộc vào các tần số lặp xung khác nhau. Tần số xung càng cao, khoảng thời gian trễ tự do giữa hai xung càng ngắn, hạt chưa kịp khuếch tán xa bởi lực Brown thì đã chịu tác động của xung tiếp theo, giúp bẫy đạt độ ổn định cao.
    • Độ trễ thời gian giữa hai xung ngược chiều ($\delta T$): Khảo sát các mức trễ $\delta T = 0$, $1\tau$, $2\tau$, $3\tau$, $4\tau$. Khi $\delta T = 0$ (hai xung đến tâm bẫy đồng thời), bẫy đạt tính đối xứng và ổn định tốt nhất. Khi độ trễ $\delta T$ tăng lên ($1\tau - 4\tau$), tính đồng thời bị phá vỡ, làm suy giảm thời gian ổn định của hạt.
    • Bán kính vi hạt ($a$): Hạt có bán kính lớn hơn ($a = 18\text{ nm}$) có thời gian ổn định $\Delta t$ dài hơn và đường kính vùng ổn định $\Delta \rho$ hẹp hơn so với hạt nhỏ ($a = 8\text{ nm}$).
    • Độ nhớt của môi trường chất lưu ($\eta$): Luận án tổng hợp bảng số liệu độ nhớt của một số chất lưu (Bảng 4.1). Môi trường có độ nhớt càng cao thì hệ số ma sát nhớt càng lớn, cản trở chuyển động khuếch tán ngẫu nhiên của lực Brown, làm giảm đường kính vùng ổn định $\Delta \rho$ và tăng thời gian giữ hạt $\Delta t$.

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

Đóng góp mới về mặt lý luận và khoa học

  1. Thiết lập mô hình động học vi hạt hoàn chỉnh cho bẫy xung ngược chiều: Luận án đã xây dựng thành công quy trình mô phỏng quá trình động học của vi hạt điện môi kích thước nanomet trong chất lưu dựa trên phương trình vi phân ngẫu nhiên Langevin, tính toán đồng thời cả quang lực xung biến thiên theo thời gian (chế độ Rayleigh) và lực thăng giáng nhiệt Brown.
  2. Đề xuất và lượng hóa khái niệm "Vùng ổn định không gian - thời gian": Luận án đã đưa ra một định nghĩa khoa học chuẩn xác để đánh giá chất lượng của kìm quang học điều khiển bằng xung laser: không xét trạng thái đứng yên tuyệt đối mà xác lập vùng ổn định qua hai đại lượng đo được là đường kính dao động $\Delta \rho$ và khoảng thời gian lưu giữ $\Delta t$.
  3. Làm rõ cơ chế cạnh tranh lực: Chứng minh định lượng quy luật chuyển dịch của hạt về tâm bẫy trong thời gian có xung và quy luật khuếch tán do dao động nhiệt trong thời gian ngắt xung; khẳng định tính đối ngẫu và sự triệt tiêu hiệu quả của lực tán xạ dọc trục trong cấu hình hai chùm xung ngược chiều.

Đóng góp mới về mặt thực tiễn và ứng dụng

  1. Xác lập bộ thông số tối ưu cho thiết kế kìm quang học: Cung cấp các đường cong định lượng chỉ rõ mối quan hệ giữa các tham số quang học của chùm tia ($U, \tau, W_0, d, \delta T, f$) và tham số môi trường vi hạt ($a, \eta$), giúp xác định các ngưỡng tham số cần thiết để giữ hạt nano ổn định.
  2. Định hướng ứng dụng vi thao tác: Kết quả tính toán cung cấp cơ sở dữ liệu khoa học tin cậy phục vụ việc định hướng thiết kế các thiết bị kìm quang học laser xung dùng trong thao tác tế bào sống, phân tử sinh học (chuỗi ADN), hoặc vi hạt kim loại/điện môi mà không gây phá hủy mẫu do nhiệt độ tích lũy cao như laser liên tục công suất lớn.

Danh mục các công trình khoa học đã công bố của tác giả

Kết quả nghiên cứu của luận án được công bố trong 07 bài báo khoa học:

  1. CT1: Hoang Dinh Hai, "Influence of Energy and Duration of Laser Pulses on Stability of Dielectric Nanoparticles in Optical Trap", Communications in Physics.
  2. CT2: Mai Van Luu, Hoang Dinh Hai and Donan Zhuang, "The Simulation of the Stabilizing Process of Dielectric Nanoparticle in Optical Trap using Counter-propagating Pulsed Laser beams", Chinese Optics Letters.
  3. CT3: Mai Van Luu, Hoang Dinh Hai, "The Influence of Parameters on Space-Time 'Pillar' in Optical Tweezers using Counter-propagating Pulsed Laser beams", Computational Methods for Science and Technology, Special Issue (2) (Ba Lan), 2010.
  4. CT4: Hoang Dinh Hai, "The simulated influence of optical parameters on stable space-time pillar of nano-particle in optical tweezers using pulsed laser beams", Journal of Science / Communications in Physics.
  5. CT5: Hoang Dinh Hai, "The simulation of the stabilizing process of glass nanoparticle in optical tweezers using series of laser pulses", Communications in Physics.
  6. CT6: Hoang Dinh Hai, "Dynamics of the dielectric nanoparticle in temporal-incoherent optical tweezers", Báo cáo Hội nghị Quang học - Quang phổ toàn quốc lần thứ VII, ISSN 1859-4271, 2012.
  7. CT7: Hoang Dinh Hai, "Dynamic of the dielectric nanoparticle in optical tweezers using counter-propagating pulsed laser beams", Journal of Physical Science and Application (USA).

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

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

  • Các kết luận khoa học chủ yếu dựa trên nền tảng tính toán lý thuyết và mô phỏng số thông qua phần mềm Matlab, chưa được đối chứng trực tiếp bằng các dữ liệu đo đạc thực nghiệm từ một hệ kìm quang học xung hoàn chỉnh trong phòng thí nghiệm.
  • Giới hạn phạm vi tính toán cho vi hạt điện môi có dạng hình cầu lý tưởng và áp dụng trong gần đúng Rayleigh ($a \ll \lambda$); chưa giải quyết cho các hạt có kích thước tương đương bước sóng ($a \approx \lambda$, chế độ Mie) hoặc các hạt có cấu trúc hình học bất đối xứng phức tạp.
  • Giả thiết hai chùm laser xung ngược chiều hoàn toàn đồng nhất về cấu trúc chùm tia, năng lượng và phân bố cường độ; trong thực tế thực nghiệm luôn tồn tại sự sai lệch nhất định về căn chỉnh quang học và chất lượng chùm tia giữa hai nhánh.

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

  • Mở rộng mô hình động học cho các hạt có hình dạng bất kỳ hoặc hạt có cấu trúc dị tính (như cấu trúc lõi-vỏ core-shell, hạt nano kim loại có hiệu ứng plasmon).
  • Mở rộng bài toán sang chế độ Mie để áp dụng cho các vi hạt và tế bào sinh học có kích thước cỡ micrômét.
  • Xây dựng hệ bố trí thí nghiệm thực tế để kiểm chứng các thông số tính toán về vùng ổn định không gian - thời gian đối với vi hạt thủy tinh và tế bào sống trong dung dịch.

Giá trị tham khảo

Luận án là tài liệu tham khảo hữu ích cho các nhóm đối tượng chuyên môn sau:

  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Vật lý quang học, Quang tử học, Vật lý kỹ thuật: Phần lý thuyết (Chương 1 và Chương 2) cung cấp phương pháp dẫn giải chi tiết từ phương trình Maxwell, quang lực Lorentz trong gần đúng Rayleigh đến phương trình động học vi phân ngẫu nhiên Langevin, phục vụ tốt cho công tác giảng dạy và nghiên cứu chuyên sâu về bẫy quang học.
  • Các nhà nghiên cứu trong lĩnh vực Lý sinh và Công nghệ nano: Cung cấp phương pháp luận và các thông số mô phỏng chính xác về hành vi chuyển động Brown của hạt kích thước nano trong chất lưu, làm cơ sở để nghiên cứu các tương tác cơ - sinh học ở cấp độ đơn phân tử.
  • Kỹ sư thiết kế thiết bị quang - cơ điện tử và quang học y sinh: Chương 3 và Chương 4 cung cấp các bảng tham số, biểu đồ tương quan về ảnh hưởng của độ rộng xung, độ trễ, tần số lặp và khẩu độ số, trực tiếp hỗ trợ việc tính toán, lựa chọn linh kiện quang học và nguồn laser phù hợp khi chế tạo kìm quang học xung ngược chiều.

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

1. Tại sao kìm quang học sử dụng một chùm xung Gauss lại khó bẫy giữ vi hạt ổn định trên trục truyền tia?

Trong bẫy đơn chùm xung Gauss, quang lực tác dụng lên hạt gồm hai thành phần: lực gradient kéo hạt vào tâm hội tụ và lực tán xạ đẩy hạt dọc theo chiều truyền tia sáng. Khi hoạt động ở chế độ xung, cường độ sáng biến thiên theo thời gian làm cho lực tán xạ biến đổi liên tục và phụ thuộc trực tiếp vào độ rộng xung. Do không có lực đối kháng dọc trục, lực tán xạ không thể tự cân bằng tại tiêu điểm, khiến hạt liên tục bị đẩy lệch khỏi vị trí bẫy trên trục truyền tia sáng.

2. Cấu hình hai chùm xung Gauss ngược chiều khắc phục nhược điểm mất ổn định dọc trục bằng cơ chế nào?

Cấu hình này sử dụng hai chùm tia giống nhau truyền ngược chiều dọc theo trục $z$, được thiết lập có phương phân cực vuông góc nhau để tránh hiện tượng giao thoa sóng. Khi đó, lực tán xạ do chùm bên trái sinh ra sẽ triệt tiêu với lực tán xạ do chùm bên phải sinh ra nhờ tính đối ngẫu qua tâm bẫy. Kết quả là lực tán xạ dọc trục bị triệt tiêu, chỉ còn lại lực gradient kéo hạt về vị trí đối xứng giữa hai chùm tia, giúp ổn định vị trí hạt dọc theo trục $z$.

3. Khái niệm "vùng ổn định không gian - thời gian" của vi hạt trong bẫy xung được định nghĩa như thế nào?

Do vi hạt chịu tác động của chuyển động nhiệt ngẫu nhiên Brown từ môi trường chất lưu, hạt không đứng yên tuyệt đối mà dao động xung quanh tâm bẫy. Luận án định nghĩa vùng ổn định là một miền không gian hình trụ ("pillar") giới hạn trên mặt phẳng tiêu bản, được xác định bởi hai thông số định lượng: đường kính vùng ổn định $\Delta \rho$ (độ dịch chuyển không gian lớn nhất của hạt quanh tâm bẫy) và thời gian ổn định $\Delta t$ (khoảng thời gian liên tục mà hạt được giam giữ trong phạm vi $\Delta \rho$ trước khi bị lực Brown đẩy ra ngoài).

4. Khi bán kính mặt thắt chùm tia $W_0$ và độ rộng xung $\tau$ tăng lên, quang lực cực đại và vùng bẫy thay đổi ra sao?

Theo kết quả mô phỏng của luận án, khi bán kính mặt thắt $W_0$ hoặc độ rộng xung $\tau$ tăng lên, phạm vi không gian của vùng bẫy trên các mặt phẳng pha có xu hướng mở rộng ra. Tuy nhiên, giá trị cực đại của quang lực ($F_{\max}$) lại suy giảm đáng kể. Cụ thể, khi $W_0$ giảm từ $2\ \mu\text{m}$ xuống $0{,}5\ \mu\text{m}$, quang lực dọc cực đại tăng vọt; tương tự, khi giảm bán độ rộng xung $\tau$ từ $1\text{ ps}$ xuống $0{,}5\text{ ps}$, giá trị quang lực dọc cực đại tăng từ $4{,}5 \cdot 10^3\ \text{pN}$ lên $8{,}5 \cdot 10^3\ \text{pN}$.

Kết luận

Luận án Tiến sĩ của Hoàng Đình Hải đã giải quyết trọn vẹn bài toán mô hình hóa động học của vi hạt điện môi kích thước nanomet trong kìm quang học sử dụng hai chùm xung Gauss ngược chiều nhúng trong chất lưu. Bằng việc giải số phương trình Langevin tích hợp đồng thời quang lực Rayleigh và lực ngẫu nhiên Brown, công trình đã làm sáng tỏ cơ chế cạnh tranh giữa bẫy quang học và thăng giáng nhiệt trong từng chu kỳ xung laser. Khái niệm "vùng ổn định không gian - thời gian" cùng hệ thống số liệu khảo sát toàn diện về các tham số quang, cơ, nhiệt là đóng góp khoa học có giá trị lý thuyết và thực tiễn cao, định hướng hiệu quả cho việc chế tạo và ứng dụng kìm quang học laser xung trong vật lý và y sinh học hiện đại.