Tổng quan luận án

Luận án tiến sĩ Vật lý với đề tài "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" được thực hiện bởi nghiên cứu sinh Hoàng Đình Hải 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ự thuộc Bộ Quốc phòng (Hà Nội, năm 2014). Công trình thuộc chuyên ngành Quang học, mã số 62.09.

Về tính cấp thiết và khoảng trống nghiên cứu, từ các nghiên cứu thực nghiệm đầu tiên của Arthur Ashkin và cộng sự vào các thập niên 1970 và 1986, bẫy quang học và kìm quang học đã trở thành công cụ quan trọng trong vật lý, hóa lý và lý sinh học nhằm giam giữ và điều khiển các đối tượng kích thước vi mô và nano như nguyên tử, chuỗi ADN, tế bào sống. Tuy nhiên, các mô hình bẫy quang học truyền thống chủ yếu hoạt động dựa trên chùm laser liên tục (CW). Khi áp dụng xung laser ngắn (như xung picosecond), nghiên cứu của Zhao và cộng sự đã chỉ ra rằng việc sử dụng một chùm xung Gauss đơn lẻ gặp hạn chế lớn: lực tán xạ biến thiên mạnh theo thời gian và phụ thuộc vào độ rộng xung, dẫn đến mất cân bằng lực tại vị trí bẫy trên trục truyền lan, khiến hạt không thể duy trì trạng thái ổn định. Để khắc phục điều này, cấu hình sử dụng hai chùm tia ngược chiều đã được M. Kawano và cộng sự đề xuất năm 2008, và sau đó được Mai Văn Lưu cùng Hồ Quang Quý phân tích lý thuyết vào năm 2009 nhằm triệt tiêu lực tán xạ thông qua tính đối ngẫu không gian. Mặc dù vậy, vấn đề phân tích định lượng sự ổn định của vi hạt điện môi trong chất lưu khi có sự cạnh tranh đồng thời giữa quang lực xung và lực ngẫu nhiên Brown, cũng như sự phụ thuộc của độ ổn định này vào các tham số quang - cơ của hệ thống, vẫn chưa được giải quyết một cách hoàn chỉnh. Luận án được triển khai nhằm giải quyết trực tiếp khoảng trống này.

Mục tiêu nghiên cứu của luận án được văn bản xác định cụ thể gồm hai nội dung chính:

  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ỡ nano 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 nghiên cứu của luận án là kìm quang học sử dụng hai chùm xung Gauss ngược chiều và đối tượng bẫy là vi hạt điện môi kích thước nano nhúng trong môi trường chất lưu.

Phạm vi nghiên cứu được giới hạn ở việc khảo sát lý thuyết tương tác trường bức xạ - vi hạt và mô phỏng số trên máy tính:

  • Về mặt đối tượng vi hạt: Xét các vi hạt điện môi hình cầu (vi cầu) đồng nhất có kích thước bán kính $a$ nằm trong khoảng $8\ \text{nm}$ đến $18\ \text{nm}$, thỏa mãn điều kiện gần đúng Rayleigh ($a \ll \lambda$), cụ thể là hạt thủy tinh lơ lửng trong môi trường nước có chiết suất $n_2 = 1{,}332$.
  • Về nguồn quang học: Hai chùm xung laser Gauss ngược chiều có bước sóng $\lambda = 1{,}064\ \mu\text{m}$, độ rộng xung ngắn cỡ picosecond ($\tau = 0{,}5 - 2\ \text{ps}$), được hội tụ bởi hệ quang có khẩu độ số trung bình $NA = 1$, bán kính thắt chùm $W_0$ từ $0{,}5\ \mu\text{m}$ đến $4\ \mu\text{m}$, năng lượng xung $U$ từ $0{,}1\ \mu\text{J}$ đến $5\ \mu\text{J}$.
  • Về mặt không gian và thời gian: Khảo sát quá trình động học và phân bố lực trong không gian ba chiều và trên các mặt phẳng pha $(\rho, z)$, $(\rho, t)$ và $(z, t)$.

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

Tổng quan tài liệu trong luận án điểm lại tiến trình phát triển của lý thuyết áp lực bức xạ và kỹ thuật bẫy quang học qua các mốc thời gian và công trình cụ thể:

  • Năm 1619, Johannes Kepler đưa ra nhận định rằng đuôi sao chổi luôn hướng ra xa Mặt Trời là do áp lực từ ánh sáng Mặt Trời.
  • Năm 1871, James Clerk Maxwell thiết lập cơ sở điện từ học chứng minh rằng sóng ánh sáng mang xung lượng và tạo ra áp lực trên bề mặt tiếp xúc khi chiếu vào vật thể.
  • Năm 1901, nhóm của Lebedev và nhóm của Nichols độc lập tiến hành các thí nghiệm đo đạc áp lực ánh sáng lên vật thể rắn, dù giá trị áp lực quan sát được rất nhỏ do cường độ dòng photon thấp.
  • Năm 1970 và 1971, sau khi nguồn laser ra đời, Arthur Ashkin đã thực hiện cân bằng áp suất ánh sáng với lực trọng trường để giữ hạt điện môi kích thước $20\ \mu\text{m}$. Đến năm 1986, Ashkin cùng Chu và cộng sự công bố bẫy quang học một chùm tia đầu tiên sử dụng thấu kính hội tụ có khẩu độ số lớn để giữ vi hạt kích thước từ $25\ \text{nm}$ đến $10\ \mu\text{m}$ trong môi trường nước.

Về phương pháp tính toán quang lực, tác giả tổng kết ba chế độ tính toán theo tỷ lệ giữa kích thước hạt $a$ và bước sóng ánh sáng $\lambda$:

  1. Chế độ Rayleigh áp dụng khi kích thước hạt nhỏ hơn rất nhiều so với bước sóng ($a \ll \lambda$).
  2. Chế độ quang hình áp dụng khi kích thước hạt lớn hơn nhiều so với bước sóng ($a \gg \lambda$).
  3. Chế độ Mie áp dụng khi kích thước hạt xấp xỉ bước sóng ($a \approx \lambda$).

Các công trình nghiên cứu mở rộng được luận án điểm lại bao gồm:

  • H. Kim khảo sát sự phụ thuộc của hiệu suất bẫy vào tỷ số giữa bán kính chùm tia và khẩu độ số $NA$ của hệ quang, qua đó đánh giá ảnh hưởng của độ nhớt chất lưu, kích thước vi hạt và công suất laser.
  • O. Stout phát triển phương pháp véc tơ để tính quang lực cho chùm tia có phân bố cường độ bất kỳ trên tiết diện ngang.
  • D. Bonessi và cộng sự kết hợp ba chế độ Rayleigh, Mie và quang hình để tính quang lực tác dụng lên các hạt có hình dạng bất kỳ.
  • Cui (2008) và Jian cùng cộng sự (2009) tập trung tính toán lực quang học tác dụng lên các vật thể kích thước nanomet.
  • Toàn bộ các hướng tiếp cận này được tổng hợp và hệ thống hóa trong công trình tổng quan của M.

Về hướng nghiên cứu sử dụng xung laser, từ năm 2005, nhiều tác giả chuyển từ nguồn laser sóng liên tục sang laser xung. Đáng chú ý, Zhao và cộng sự đã công bố kết quả tính toán quang lực cho kìm sử dụng một chùm xung Gauss, phân tích 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)$. Nghiên cứu của Zhao rút ra các nhận định: chuyển động Brown có thể bị triệt tiêu tạm thời khi hạt chịu tác động của xung ngắn; tuy nhiên giữa hai xung liên tiếp, nếu tần số 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$. Đặc biệt, Zhao chỉ ra rằng một chùm xung Gauss đơn lẻ không thể tạo ra bẫy quang học ổn định vì lực tán xạ biến thiên liên tục theo thời gian và độ rộng xung, dẫn đến tình trạng mất cân bằng lực đẩy tại điểm bẫy.

Để giải quyết nhược điểm này, M. Kawano và cộng sự (2008) đề xuất mô hình kìm quang học hai chùm tia ngược chiều. Năm 2009, Mai Văn Lưu và Hồ Quang Quý nghiên cứu quang lực của hệ hai chùm xung Gauss ngược chiều, chứng minh lực gradient và lực tán xạ có tính đối xứng qua tâm bẫy, giúp triệt tiêu hoàn toàn lực tán xạ dọc trục khi các tham số của hai chùm được cân bằng.

Khoảng trống nghiên cứu được luận án xác định: Dù các công trình trước đã khảo sát phân bố lực và các vùng bẫy hình thức trên mặt phẳng pha, nhưng việc phân tích định lượng sự ổn định của vi hạt dưới tác động đồng thời của quang lực xung và lực ngẫu nhiên Brown trong chất lưu vẫn chưa được thực hiện. Luận án lựa chọn giải quyết bài toán mô hình hóa động học vi hạt và xác định các điều kiện ổn định không gian - thời gian của hệ kìm hai chùm xung Gauss ngược chiều.

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

Luận án xây dựng cơ sở lý thuyết dựa trên mô hình tương tác trường điện từ với vi hạt điện môi hình cầu trong chế độ tán xạ Rayleigh ($a \ll \lambda$). Hạt được coi như một lưỡng cực điện cảm ứng đặt trong trường bức xạ, chịu tác dụng của lực Lorentz trung bình theo thời gian: $$\vec{F}P(\rho, z, t) = \vec{F}{grad} + \vec{F}t$$ Khi lấy trung bình theo thời gian trong điều kiện công suất laser ổn định, thành phần phụ thuộc đạo hàm từ trường triệt tiêu, lực tác dụng rút về lực gradient: $$\vec{F}{grad}(\rho, z, t) = \frac{1}{2} \sigma \nabla \vec{E}^2(\rho, z, t)$$ Trong đó, hệ số phân cực của hạt vi cầu điện môi được xác định theo biểu thức: $$\sigma = 4\pi n_2^2 \epsilon_0 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, $\epsilon_0$ là hằng số điện môi chân không, và $a$ là bán kính vi hạt.

Về cấu hình quang học, hệ thống gồm hai chùm xung laser Gauss lan truyền ngược chiều dọc theo trục $z$, được định hướng hội tụ đối xứng qua tâm kìm (gốc tọa độ $z = 0$). Để ngăn chặn hiện tượng giao thoa sóng tạo ra các vân nhiễu làm mất ổn định bẫy, hai chùm tia được phân cực vuông góc với nhau: chùm bên trái phân cực theo phương $\vec{x}$, chùm bên phải phân cực theo phương $\vec{y}$. Mặt thắt của chùm bên trái đặt tại tọa độ $z = -d/2$, mặt thắt của chùm bên phải đặt tại $z = d/2$, với $d$ là khoảng cách giữa hai mặt thắt.

Cường độ điện trường của chùm tia bên trái và bên phải được thiết lập theo các công thức: $$\vec{E}_l(\rho, z, t, d) = \vec{x} E_0 \frac{i k W_0^2}{i k W_0^2 + 2(z + d/2)} \exp\left{ -i\left[ k(z + d/2) - \omega_0 t \right] \right} \exp\left[ \frac{-2 i k (z + d/2) \rho^2}{(k W_0^2)^2 + 4(z + d/2)^2} \right] \exp\left[ -\frac{(k W_0^2)^2 \rho^2}{(k W_0^2)^2 + 4(z + d/2)^2} \right] \exp\left{ -\frac{[t - (z + d/2)/c]^2}{\tau^2} \right}$$ $$\vec{E}_r(\rho, z, t, d) = \vec{y} E_0 \frac{i k W_0^2}{i k W_0^2 + 2(z - d/2)} \exp\left{ -i\left[ k(z - d/2) - \omega_0 t \right] \right} \exp\left[ \frac{-2 i k (z - d/2) \rho^2}{(k W_0^2)^2 + 4(z - d/2)^2} \right] \exp\left[ -\frac{(k W_0^2)^2 \rho^2}{(k W_0^2)^2 + 4(z - d/2)^2} \right] \exp\left{ -\frac{[t + (z - d/2)/c]^2}{\tau^2} \right}$$ Biên độ điện trường $E_0$ tại trục liên hệ với năng lượng tổng $U$ của xung qua hệ thức: $$E_0^2 = \frac{4\sqrt{2} U}{n_2 \epsilon_0 c W_0^2 \pi^{3/2} \tau}$$ Cường độ chùm tia tương ứng với độ lớn véc tơ Poynting trung bình theo thời gian được tính toán độc lập cho từng chùm ($I_l$ và $I_r$), và cường độ tổng của hệ là $I = I_l + I_r$.

Từ trường cường độ này, các thành phần quang lực tác dụng lên hạt điện môi gồm:

  1. Lực tán xạ: $$\vec{F}{scat}(\rho, z, t, d) = \vec{z} \frac{n_2}{c} C{pr} I_l(\rho, z, t, d) - \vec{z} \frac{n_2}{c} C_{pr} I_r(\rho, z, t, d)$$ với tiết diện tán xạ $C_{pr} = \frac{128\pi^5 a^6}{3\lambda^4} \left( \frac{m^2 - 1}{m^2 + 2} \right)^2$.
  2. Lực gradient hướng tâm (quang lực ngang): $$\vec{F}_{grad,\rho}(\rho, z, t, d) = -\vec{\rho} \frac{2\sigma I_l \tilde{\rho}}{1 + 4(\tilde{z} + \tilde{d})^2} - \vec{\rho} \frac{2\sigma I_r \tilde{\rho}}{1 + 4(\tilde{z} - \tilde{d})^2}$$
  3. Lực gradient dọc trục kết hợp với lực tán xạ tạo thành quang lực dọc $F_z = F_{scat} + F_{grad,z}$.

Về mặt động học chất lưu, chuyển động của vi hạt được mô tả bằng phương trình vi phân ngẫu nhiên Langevin: $$m \frac{d\vec{\upsilon}}{dt} = \vec{f}{vis} + \vec{f}{total}$$ Trong đó:

  • $\vec{f}_{vis} = -\alpha \vec{\upsilon} = -6\pi\eta a \vec{\upsilon}$ là lực 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, $a$ là bán kính vi hạt).
  • $\vec{f}{total} = \vec{F}N + \vec{f}{Brown} + \vec{f}{gravity} + \vec{f}{Hydrate}$, bao gồm quang lực tổng hợp $\vec{F}N$, lực ngẫu nhiên Brown $\vec{f}{Brown}$, trọng lực $\vec{f}{gravity}$ và lực thủy động đàn hồi $\vec{f}_{Hydrate}$. Lực Brown biểu diễn dao động nhiệt ngẫu nhiên qua hàm phân bố ngẫu nhiên $h(t)$.

Phương pháp nghiên cứu thực tế là kết hợp giải tích lý thuyết và mô phỏng số trên phần mềm Matlab. Quy trình tính toán số tích phân phương trình Langevin theo các bước thời gian $\delta t$ từ $6\tau/2000$ đến $6\tau/10000$ để theo dõi tọa độ và vận tốc của hạt theo thời gian. Bộ thông số mô phỏng chuẩn bao gồm: bước sóng laser $\lambda = 1{,}064\ \mu\text{m}$, môi trường nước có $n_2 = 1{,}332$, hạt thủy tinh hình cầu bán kính $a = 10\ \text{nm}$, khẩu độ số $NA = 1$, bán kính thắt chùm $W_0 = 1\ \mu\text{m}$, độ rộng xung $\tau = 1\ \text{ps}$, năng lượng xung $U = 0{,}1\ \mu\text{J}$ hoặc $0{,}9\ \mu\text{J}$.

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 trình bày chi tiết cơ sở lý thuyết về quang lực, cấu hình hệ quang học hai chùm xung Gauss ngược chiều, và phân tích các yếu tố ảnh hưởng đến cường độ trường cũng như phân bố quang lực tác dụng lên hạt điện môi.

Nội dung khảo sát tập trung vào các kết quả mô phỏng không gian và thời gian:

  1. Phân bố cường độ tổng $I = I_l + I_r$ trong mặt phẳng pha $(\rho, z)$ tại thời điểm $t = \tau = 1\ \text{ps}$: Khi khảo sát với các khoảng cách $d$ khác nhau ($15\ \mu\text{m}, 10\ \mu\text{m}, 5\ \mu\text{m}, 0\ \mu\text{m}$), tác giả phát hiện tại $d = 15\ \mu\text{m}$ và $d = 10\ \mu\text{m}$, cường độ tổng hình thành hai đỉnh tách biệt dọc theo trục $z$. Khi giảm khoảng cách về $d \le 5\ \mu\text{m}$, hai đỉnh hợp nhất thành một đỉnh duy nhất có cường độ cực đại tại tâm bẫy ($z = 0, \rho = 0$). Khảo sát ảnh hưởng của bán kính thắt chùm $W_0$ ($2{,}0\ \mu\text{m}; 1{,}5\ \mu\text{m}; 1{,}0\ \mu\text{m}; 0{,}5\ \mu\text{m}$) tại $d = 10\ \mu\text{m}$ cho thấy khi $W_0$ giảm, cường độ cực đại tăng từ $3500\ \mu\text{W}/\mu\text{m}^2$ lên $9000\ \mu\text{W}/\mu\text{m}^2$ và phân bố Gauss bị phá vỡ cục bộ.
  2. Phân bố quang lực dọc $F_{grad,z}$ trong mặt phẳng pha $(z, t)$: Quang lực dọc hoàn toàn đối xứng qua gốc tọa độ $z = 0$, hình thành hai vùng bẫy ổn định dọc theo chiều truyền ánh sáng. Khi bán kính mặt thắt $W_0 = 0{,}5\ \mu\text{m}$, lực gradient dọc cực đại đạt giá trị $40 \times 10^3\ \text{pN}$; khi $W_0$ tăng lên $1{,}0\ \mu\text{m}$, lực cực đại giảm xuống $4{,}5 \times 10^3\ \text{pN}$ nhưng kích thước vùng bẫy mở rộng và hai vùng có xu hướng chồng lấn. Khảo sát theo độ rộng xung $\tau$ cho thấy với xung ngắn $\tau = 0{,}5\ \text{ps}$, lực cực đại là $8{,}5 \times 10^3\ \text{pN}$, trong khi với $\tau = 1\ \text{ps}$ lực giảm còn $4{,}5 \times 10^3\ \text{pN}$. Khi tăng khoảng cách $d$ ($5\ \mu\text{m}, 10\ \mu\text{m}, 15\ \mu\text{m}, 20\ \mu\text{m}$), hai vùng bẫy tách xa nhau, tạo ra khoảng ranh giới tự do ở giữa khiến hạt có thể dao động nhiệt không bị kiểm soát.
  3. Phân bố quang lực ngang $F_{grad,\rho}$ trong mặt phẳng pha $(\rho, t)$: Quang lực ngang tạo thành một hố thế thế năng đối xứng qua tâm $\rho = 0$, luôn có xu hướng kéo hạt về trục chùm tia. Khảo sát sự phụ thuộc của giá trị cực đại quang lực ngang vào bán kính mặt thắt tại $t = 1\tau, d = 10\ \mu\text{m}, z = 0$ cho thấy lực tăng dần và đạt giá trị cực đại tại $W_0 = 1{,}05\ \mu\text{m}$, sau đó giảm dần khi tiếp tục tăng $W_0$. Khi độ rộng xung $\tau$ tăng ($0{,}5\ \text{ps}; 1\ \text{ps}; 1{,}5\ \text{ps}$), giá trị lực ngang cực đại giảm nhưng thời gian tồn tại của hố thế trong chu kỳ xung kéo dài hơn.

Chương này cũng thiết lập cơ sở về chuyển động Brown và phương trình Langevin, phân tích lực ma sát nhớt $f_{vis} = -6\pi\eta a \upsilon$ theo định luật Stokes và tổng hợp các thành phần lực tác dụng lên hạt.

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

Chương 2 tập trung vào việc mô hình hóa toán học và thuật toán mô phỏng quá trình động học của vi hạt trong chất lưu dưới tác động kết hợp của lực ngẫu nhiên Brown và quang lực từ hai chùm xung Gauss ngược chiều.

Các nội dung và kết quả chính của chương gồm:

  • Thiết lập phương trình Langevin chi tiết cho hệ tọa độ cực và tọa độ hướng tâm, tích hợp biểu thức quang lực biến thiên theo thời gian của xung laser picosecond.
  • Xây dựng quy trình thuật toán mô phỏng trên Matlab. Tác giả khảo sát kiểm tra tính hội tụ và ổn định của phương pháp mô phỏng số đối với chuyển động Brown tự do của vi hạt thủy tinh trong nước tại các bước thời gian khác nhau: $\delta t = 6\tau/2000, 6\tau/4000, 6\tau/8000, 6\tau/10000$, ứng với tọa độ ban đầu $\rho_0 = 1\ \mu\text{m}$ và $\rho_0 = 0\ \mu\text{m}$.
  • Khảo sát sự cạnh tranh giữa quang lực và lực Brown trong thời gian một xung tác dụng: Khi xung laser xuất hiện và đạt đỉnh cường độ, quang lực gradient tăng vọt, chiếm ưu thế hoàn toàn so với lực dao động nhiệt ngẫu nhiên, kéo vi hạt nhanh chóng từ vị trí ban đầu dịch chuyển về tâm bẫy ($\rho = 0$). Sau khi xung laser đi qua, quang lực suy giảm, hạt chịu tác động chính của lực Brown và dao động ngẫu nhiên xung quanh vị trí cân bằng. Qua đó, chương xác lập khái niệm về thời gian bẫy hạt trong mỗi chu kỳ 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 phân tích định lượng ảnh hưởng của các tham số đầu vào của hệ quang và vi hạt lên tốc độ dịch chuyển về tâm kìm và biên độ dao động của hạt tại tâm.

Các tham số khảo sát gồm:

  1. Ảnh hưởng của vị trí ban đầu của vi hạt ($\rho_0$): Khi vi hạt ở các vị trí ban đầu khác nhau ($\rho_0 = 1\ \mu\text{m}, 2\ \mu\text{m}$), quỹ đạo động học cho thấy hạt nằm càng xa tâm thì gradient lực ngang càng lớn tại giai đoạn đầu của xung, kéo hạt chuyển động nhanh về tâm; tuy nhiên nếu hạt nằm ngoài phạm vi hố thế của chùm tia, lực quang học không đủ lớn để thắng lực Brown và hạt sẽ khuếch tán ra ngoài.
  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{,}1\ \mu\text{J}; 0{,}9\ \mu\text{J}; 5\ \mu\text{J}$. Khi tăng năng lượng xung, độ sâu hố thế quang học tăng, gia tốc kéo hạt về tâm tăng rõ rệt, đồng thời biên độ dao động nhiệt ngẫu nhiên của hạt tại tâm bẫy bị nén chặt hơn.
  3. Ảnh hưởng của bán kính mặt thắt chùm tia ($W_0$): Với các thông số $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}$, khảo sát các giá trị $W_0 = 1\ \mu\text{m}, 2\ \mu\text{m}, 3\ \mu\text{m}, 4\ \mu\text{m}$. Kết quả chỉ ra rằng tại $W_0 = 1\ \mu\text{m}$, quang lực gradient tập trung cao giúp hạt nhanh chóng bị hút về tâm bẫy. Khi tăng $W_0$ lên $2 - 4\ \mu\text{m}$, gradient cường độ giảm, hố thế thoải hơn làm giảm tốc độ dịch chuyển của hạt và thời gian bẫy ổn định bị rút ngắn.
  4. Ảnh hưởng của kích thước hạt ($a$): Khảo sát với các bán kính vi hạt $a = 8\ \text{nm}, 10\ \text{nm}, 14\ \text{nm}, 18\ \text{nm}$ tại $U = 0{,}9\ \mu\text{J}, W_0 = 1\ \mu\text{m}$. Vì lực gradient tỷ lệ với $a^3$ trong khi lực cản nhớt tỷ lệ với $a$, các hạt có kích thước lớn hơn ($14 - 18\ \text{nm}$) chịu quang lực kéo mạnh hơn nhưng cũng chịu lực cản nhớt lớn hơn, dẫn đến đặc tính dao động quanh tâm bẫy có tần số và biên độ khác biệt so với các hạt kích thước nhỏ ($8 - 10\ \text{nm}$).

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

Chương 4 là chương đóng góp trọng tâm về mặt định lượng, đề xuất khái niệm "vùng ổn định không gian - thời gian" của vi hạt trong kìm quang học. Vùng ổn định này được đặc trưng bởi hai thông số: khoảng thời gian ổn định của hạt trong bẫy ($\Delta t$) và đường kính vùng ổn định trên mặt phẳng vuông góc với trục truyền ($\Delta\rho$).

Chương khảo sát hệ thống ảnh hưởng của các nhóm tham số lên $\Delta t$ và $\Delta\rho$:

  1. Năng lượng xung laser ($U$): Tăng năng lượng xung làm tăng thời gian ổn định $\Delta t$ và mở rộng đường kính vùng ổn định $\Delta\rho$, giúp hạt được giam giữ chắc chắn hơn trong hố thế quang học.
  2. Bán kính thắt chùm ($W_0$): Sự phụ thuộc của $\Delta t$ và $\Delta\rho$ vào $W_0$ cho thấy một điểm cực trị; vùng ổn định đạt giá trị tối ưu quanh $W_0 \approx 1\ \mu\text{m}$. Khi $W_0$ quá nhỏ, phạm vi không gian của bẫy bị thu hẹp; khi $W_0$ quá lớn, độ sâu hố thế suy giảm làm giảm thời gian ổn định.
  3. Độ rộng xung laser ($\tau$): Xung có độ rộng $\tau$ lớn hơn cung cấp thời gian tương tác quang lực dài hơn trong mỗi chu kỳ, làm tăng thời gian bẫy hạt $\Delta t$.
  4. Tần số lặp xung laser ($f$): Khảo sát sự ổn định của vi hạt theo các tần số lặp xung khác nhau: $f = 1/6, 1/4, 1/3, 1/2$. Tần số lặp xung cao giúp giảm khoảng thời gian gián đoạn giữa các xung, không cho phép lực ngẫu nhiên Brown đẩy hạt trôi xa khỏi tâm bẫy, duy trì trạng thái giam giữ liên tục.
  5. Độ trễ thời gian giữa hai xung ($\delta T$): Khảo sát ảnh hưởng khi có độ trễ giữa hai chùm xung với các giá trị $\delta T = 0, 1\tau, 2\tau, 3\tau, 4\tau$. Khi $\delta T = 0$ (hai xung tới đồng thời tại tâm), bẫy đạt trạng thái ổn định cao nhất do lực tán xạ bị triệt tiêu hoàn toàn. Khi $\delta T$ tăng từ $1\tau$ đến $4\tau$, tính đồng thời bị phá vỡ, lực tán xạ của từng chùm không còn triệt tiêu lẫn nhau tại mỗi thời điểm, dẫn đến hố thế bị biến dạng và vùng ổn định bị thu hẹp nghiêm trọng.
  6. Bán kính vi hạt ($a$): Ảnh hưởng rõ rệt đến thời gian $\Delta t$ và đường kính $\Delta\rho$ do sự thay đổi đồng thời của quang lực và hệ số ma sát Stokes.
  7. Độ nhớt chất lưu ($\eta$): Khảo sát sự phụ thuộc của thời gian ổn định và đường kính vùng ổn định vào độ nhớt môi trường. Văn bản tổng hợp số liệu độ nhớt của một số chất lưu tại Bảng 4.1. Độ nhớt chất lưu càng lớn thì lực cản nhớt $f_{vis}$ càng mạnh, cản trở chuyển động khuếch tán ngẫu nhiên của lực Brown, giúp tăng thời gian ổn định $\Delta t$ của vi hạt trong bẫy.

Bảng tổng hợp dưới đây hệ thống hóa các thông số nghiên cứu chính được sử dụng và khảo sát trong toàn bộ luận án:

Tên thông số Ký hiệu Đơn vị Giá trị chuẩn / Khoảng khảo sát
Bước sóng chùm laser $\lambda$ $\mu\text{m}$ $1{,}064$
Khẩu độ số hệ vật kính $NA$ - $1{,}0$
Bán kính mặt thắt chùm tia $W_0$ $\mu\text{m}$ $0{,}5;\ 1{,}0;\ 1{,}05;\ 1{,}5;\ 2{,}0;\ 3{,}0;\ 4{,}0$
Khoảng cách giữa hai mặt thắt $d$ $\mu\text{m}$ $0;\ 5;\ 10;\ 15;\ 20$
Bán độ rộng xung laser $\tau$ $\text{ps}$ $0{,}5;\ 1{,}0;\ 1{,}5;\ 2{,}0$
Năng lượng xung laser $U$ $\mu\text{J}$ $0{,}1;\ 0{,}9;\ 5{,}0$
Tần số lặp xung (tương đối) $f$ - $1/6;\ 1/4;\ 1/3;\ 1/2$
Độ trễ giữa hai xung $\delta T$ $\tau$ $0;\ 1\tau;\ 2\tau;\ 3\tau;\ 4\tau$
Bán kính hạt điện môi (vi cầu) $a$ $\text{nm}$ $8;\ 10;\ 14;\ 18$
Chiết suất môi trường (nước) $n_2$ - $1{,}332$
Chiết suất hạt bẫy $n_1$ - Thủy tinh ($m = n_1/n_2 > 1$)
Tọa độ ban đầu của hạt $\rho_0$ $\mu\text{m}$ $0;\ 1;\ 2$
Bước thời gian mô phỏng $\delta t$ - $6\tau/2000;\ 6\tau/4000;\ 6\tau/8000;\ 6\tau/10000$

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

Các kết quả và đóng góp mới của luận án được thể hiện qua các luận điểm khoa học và thực tiễn cụ thể:

  • Về mặt lý luận:
    1. Luận án đã thiết lập thành công mô hình lý thuyết và 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 kìm quang học sử dụng hai chùm xung Gauss ngược chiều phân cực vuông góc, xét đến tác động đồng thời của quang lực biến thiên và lực ngẫu nhiên Brown trong môi trường chất lưu dựa trên phương trình vi phân ngẫu nhiên Langevin.
    2. Chứng minh một cách định lượng tính đối xứng của trường quang lực trong cấu hình hai chùm xung ngược chiều: lực tán xạ dọc trục bị triệt tiêu tại tâm bẫy, khắc phục nhược điểm mất ổn định của kìm quang học sử dụng một chùm xung Gauss.
    3. Đề xuất khái niệm khoa học mới về "vùng ổn định không gian - thời gian" của vi hạt trong bẫy quang học, xác định định lượng thông qua đường kính vùng ổn định $\Delta\rho$ và thời gian ổn định $\Delta t$, cho phép đánh giá mức độ giam giữ vi hạt một cách chính xác thay vì chỉ xem xét hố thế tĩnh.
  • Về mặt thực tiễn và ứng dụng:
    1. Xác định được miền tham số tối ưu cho hệ kìm quang học hai chùm xung Gauss ngược chiều: Bán kính mặt thắt chùm tia đạt hiệu quả bẫy quang lực ngang cực đại tại $W_0 \approx 1{,}05\ \mu\text{m}$; khoảng cách giữa hai mặt thắt chùm tia nên chọn $d \le 5\ \mu\text{m}$ để cường độ tổng tạo thành đỉnh đơn tập trung năng lượng tại tâm bẫy; độ trễ giữa hai xung phải đảm bảo $\delta T = 0$ để tránh làm lệch hố thế và phá vỡ sự triệt tiêu lực tán xạ.
    2. Xác lập mối liên hệ giữa các tham số laser (năng lượng xung $U$, độ rộng xung $\tau$, tần số lặp $f$), tham số hình học chùm tia ($W_0, d$), tính chất hạt ($a$) và môi trường chất lưu ($\eta$), cung cấp cơ sở dữ liệu tính toán định hướng cho việc thiết kế, chế tạo và vận hành các hệ kìm quang học thực nghiệm nhằm bẫy giữ các đối tượng nano, phân tử sinh học và tế bào sống mà không làm phá hủy mẫu do tích tụ nhiệt.

Các kết quả nghiên cứu trong luận án đã được công bố trong 07 công trình khoa học (gồm các bài báo trên tạp chí Communications in Physics, Chinese Optics Letters, Computational Methods for Science and Technology, Journal of Science and Arts, Journal of Physical Science and Application, và kỷ yếu Hội nghị Quang học - Quang phổ toàn quốc).

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

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

  • Toàn bộ các kết quả nghiên cứu trong luận án được thực hiện bằng phương pháp giải tích lý thuyết và mô phỏng số trên phần mềm máy tính (Matlab), chưa có điều kiện kiểm chứng trực tiếp bằng số liệu đo đạc thực nghiệm trên hệ kìm quang học vật lý.
  • Đối tượng hạt bẫy được giới hạn ở dạng vi cầu điện môi lý tưởng đồng nhất, đẳng hướng và kích thước nằm trong chế độ tán xạ Rayleigh ($a \ll \lambda$), chưa mở rộng sang các hạt có kích thước trung bình cỡ bước sóng (chế độ Mie) hoặc các vi hạt có hình dạng bất đối xứng.
  • Bỏ qua ảnh hưởng của hiện tượng hấp thụ ánh sáng của vi hạt và môi trường chất lưu, cũng như bỏ qua các hiệu ứng quang phi tuyến và biến đổi nhiệt độ cục bộ của chất lưu dưới tác động của chuỗi xung laser.

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

  • Xây dựng hệ thống thực nghiệm kìm quang học sử dụng hai chùm xung Gauss ngược chiều để kiểm chứng và hiệu chuẩn các kết quả mô phỏng số về vùng ổn định không gian - thời gian.
  • Mở rộng phạm vi tính toán cho các vi hạt điện môi có kích thước lớn hơn trong chế độ Mie và các hạt có hình học phức tạp (như phân tử sinh học hình que, chuỗi xoắn ADN).
  • Khảo sát bổ sung các hiệu ứng truyền nhiệt, vi dòng đối lưu chất lưu do laser gây ra và ảnh hưởng của sự bất đối xứng năng lượng giữa hai chùm tia trong điều kiện thực tế.

Giá trị tham khảo

Luận án có giá trị tham khảo chuyên môn đối với nhiều nhóm đối tượng học thuật và kỹ thuật:

  • Đối với nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Quang học, Quang điện tử, Vật lý kỹ thuật: Luận án là tài liệu tham khảo chi tiết về phương pháp thiết lập biểu thức quang lực xung picosecond và quy trình mô phỏng phương trình vi phân ngẫu nhiên Langevin trên Matlab để giải quyết bài toán động học vi hạt trong chất lưu.
  • Đối với các kỹ sư quang học và nhà thực nghiệm phát triển thiết bị bẫy quang học: Các đồ thị khảo sát và bộ tham số định lượng ở Chương 1 và Chương 4 (về bán kính thắt chùm tối ưu $W_0 \approx 1{,}05\ \mu\text{m}$, khoảng cách hai mặt thắt $d \le 5\ \mu\text{m}$, yêu cầu đồng bộ độ trễ xung $\delta T = 0$) là cơ sở thiết kế hệ quang hội tụ, lựa chọn khẩu độ số $NA$ của vật kính và tinh chỉnh bộ trễ xung quang học.
  • Đối với các nhà nghiên cứu trong lĩnh vực lý sinh và công nghệ nano: Công trình cung cấp cơ sở lý thuyết về việc sử dụng xung laser ngắn để giam giữ vi hạt nano, giúp giảm thiểu hư tổn do hiệu ứng nhiệt tích tụ của laser liên tục khi thao tác trên các mẫu sinh học nhạy cảm.

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

1. Vì sao kìm quang học sử dụng một chùm xung Gauss đơn lẻ không thể tạo ra bẫy quang học ổn định trên trục truyền lan?
Theo kết quả phân tích lý thuyết của Zhao và được luận án khẳng định lại, đối với một chùm xung Gauss đơn lẻ, lực tán xạ (đẩy hạt theo chiều truyền) biến thiên liên tục theo thời gian và phụ thuộc vào độ rộng xung. Tại vị trí bẫy, lực tán xạ này không được bù trừ hay cân bằng, khiến hạt liên tục bị đẩy dọc theo trục truyền lan của ánh sáng thay vì được giữ cố định.

2. Cấu hình hai chùm xung Gauss ngược chiều triệt tiêu lực tán xạ như thế nào?
Luận án bố trí hai chùm xung Gauss lan truyền ngược chiều nhau dọc theo trục $z$ với các thông số năng lượng và hình học đối xứng qua gốc tọa độ $z = 0$. Khi đó, lực tán xạ do chùm bên trái sinh ra hướng theo chiều $+z$ sẽ có độ lớn bằng và ngược chiều với lực tán xạ do chùm bên phải sinh ra hướng theo chiều $-z$. Lực tán xạ tổng hợp triệt tiêu tại tâm bẫy ($z = 0$), giúp hệ thống chỉ còn tác dụng của lực hồi phục gradient.

3. Hiện tượng giao thoa giữa hai chùm xung ngược chiều được khắc phục bằng cách nào?
Để tránh hiện tượng giao thoa ánh sáng tạo ra hệ vân sóng đứng làm biến dạng hố thế quang học, luận án sử dụng kính phân cực để thiết lập hướng phân cực của hai chùm tia vuông góc với nhau (chùm bên trái phân cực theo trục $\vec{x}$, chùm bên phải phân cực theo trục $\vec{y}$). Nhờ đó, hai chùm không giao thoa và cường độ tổng của hệ đơn giản là tổng độ lớn cường độ của từng chùm riêng biệt ($I = I_l + I_r$).

4. Bán kính mặt thắt chùm tia $W_0$ tối ưu để đạt quang lực ngang cực đại là bao nhiêu?
Kết quả khảo sát trong mặt phẳng pha $(\rho, t)$ tại vị trí $z = 0, t = 1\tau, d = 10\ \mu\text{m}$ cho thấy giá trị cực đại của quang lực ngang tăng dần theo $W_0$ và đạt đỉnh cao nhất tại $W_0 = 1{,}05\ \mu\text{m}$, sau đó giảm dần khi $W_0$ tiếp tục tăng. Do đó, $W_0 \approx 1{,}05\ \mu\text{m}$ là giá trị tối ưu để hệ kìm đạt hiệu quả giữ hạt theo phương ngang tốt nhất.

5. Khoảng cách giữa hai mặt thắt $d$ ảnh hưởng như thế nào đến cường độ tổng và sự ổn định của bẫy?
Khi khoảng cách $d$ lớn ($d = 10\ \mu\text{m}$ hoặc $15\ \mu\text{m}$), phân bố cường độ tổng hình thành hai đỉnh tách biệt và quang lực dọc tạo ra hai vùng bẫy rời rạc, để lại một khoảng ranh giới ở giữa khiến hạt có thể dao động Brown tự do. Khi giảm khoảng cách xuống $d \le 5\ \mu\text{m}$, hai đỉnh hợp nhất thành một đỉnh duy nhất tại tâm ($z = 0, \rho = 0$), tạo ra hố thế sâu và tập trung, giúp hạt được giữ ổn định tại gốc tọa độ.

Kết luận

Luận án tiến sĩ của tác giả Hoàng Đình Hải đã giải quyết bài toán mô hình hóa lý thuyết và mô phỏng độ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 phân cực vuông góc dưới tác động đồng thời của quang lực và lực ngẫu nhiên Brown trong chất lưu. Bằng việc giải phương trình vi phân ngẫu nhiên Langevin trên Matlab, tác giả đã chứng minh sự triệt tiêu hiệu dụng của lực tán xạ dọc trục, xác lập khái niệm và tiêu chí định lượng về vùng ổn định không gian - thời gian của kìm quang học. Công trình đã tìm ra bộ thông số tối ưu về bán kính thắt chùm ($W_0 \approx 1{,}05\ \mu\text{m}$), khoảng cách hai mặt thắt ($d \le 5\ \mu\text{m}$) và yêu cầu đồng bộ xung ($\delta T = 0$), đóng góp luận cứ khoa học quan trọng cho việc định hướng thiết kế và ứng dụng các hệ bẫy quang học phát xung trong vật lý và y sinh học.