Tổng quan nghiên cứu

Trong kỹ thuật vi lưu và y sinh học hiện đại, việc điều khiển và phân tách các phần tử kích thước micro-nano đòi hỏi nguồn từ trường có gradient cực lớn, với tích giá trị cảm ứng từ và độ biến thiên cảm ứng từ ($B \cdot \nabla B$) đạt từ $10^3\ \text{T}^2/\text{m}$ đến $10^5\ \text{T}^2/\text{m}$. Các phương pháp truyền thống sử dụng cuộn dây solenoid hoặc nam châm điện thường cồng kềnh, tiêu thụ năng lượng cao và khó thu nhỏ kích thước. Nhằm giải quyết thách thức này, nghiên cứu tập trung vào việc mô phỏng và khảo sát phân bố từ trường của các vi cấu trúc từ cứng vĩnh cửu có tính dị hướng từ cao, mở ra hướng đi mới cho việc tích hợp vào các hệ thống vi cơ điện tử từ tính (MagMEMS) và chip sinh học Lab-on-a-Chip.

Mục tiêu cụ thể của luận văn là tính toán giải tích và mô phỏng số sự phân bố thành phần cảm ứng từ $B_z$ cùng sự biến thiên không gian của nó theo các phương ($\frac{dB_z}{dy}$ và $\frac{dB_z}{dz}$) xung quanh các vi nam châm hình trụ làm từ vật liệu $Nd_2Fe_{14}B$. Cấu trúc khảo sát có bán kính $R = 25\ \mu\text{m}$, chiều cao $h = 5\ \mu\text{m}$ và độ từ dư $B_r = 1.2\ \text{T}$. Phạm vi nghiên cứu được thực hiện tại Khoa Vật lý kỹ thuật và Công nghệ nano, Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội vào năm 2016, với sự tài trợ một phần từ đề tài mã số 103.80 của Quỹ phát triển Khoa học và Công nghệ Quốc gia (NAFOSTED).

Kết quả nghiên cứu cung cấp cơ sở dữ liệu định lượng chính xác ở 5 mức độ cao khảo sát từ $0\ \mu\text{m}$ đến $100\ \mu\text{m}$, giúp các nhà khoa học tối ưu hóa cấu hình bẫy từ trước khi tiến hành chế tạo thực nghiệm, từ đó tiết kiệm đáng kể thời gian và chi phí nghiên cứu.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng hệ phương trình Maxwell cho trường từ tĩnh trong môi trường không có dòng điện tự do, với điều kiện phân kỳ của cảm ứng từ bằng 0 và độ xoay của cường độ từ trường bằng 0. Để xác định trường từ xung quanh các cấu trúc vi nam châm, hai mô hình lý thuyết cốt lõi đã được áp dụng đồng thời:

Mô hình dòng tương đương (mô hình Ampere) coi nam châm vĩnh cửu có từ độ $M$ tương đương với một cuộn dây có dòng điện mặt $I_{\text{eq}}$ chạy quanh bề mặt ngoài. Từ trường tại một điểm bất kỳ trong không gian được tính toán thông qua tích phân định luật Biot-Savart. Mô hình từ tích (mô hình Coulomb) xem nam châm là tập hợp các lưỡng cực từ tạo ra từ thế vô hướng $\Phi_m$ trên các mặt đáy, từ đó suy ra vectơ cảm ứng từ thông qua gradient của từ thế.

Nghiên cứu ứng dụng các khái niệm vật liệu học quan trọng gồm: vật liệu từ cứng $Nd_2Fe_{14}B$ với cấu trúc tinh thể tứ giác, dị hướng từ tinh thể một trục lớn, lực kháng từ $H_c > 10\ \text{kOe}$, từ độ bão hòa $M_s = 1.62\ \text{T}$ và tích năng lượng từ cực đại $(BH)_{\text{max}}$ vượt trội; hiện tượng từ trễ; và lực từ tác dụng lên phần tử từ tính trong dung dịch lỏng tỷ lệ thuận với thể tích hạt, độ chênh lệch độ cảm từ $\Delta \chi$ và tích $B \cdot \nabla B$.

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

Nguồn dữ liệu của luận văn được xây dựng từ quy trình mô phỏng số và tính toán giải tích đa tầng. Cỡ mẫu nghiên cứu bao gồm 5 tập cấu hình mảng vi nam châm tuần hoàn đối xứng: $1\times1, 2\times2, 3\times3, 4\times4$ và $5\times5$. Phương pháp chọn mẫu cấu hình dựa trên nguyên tắc tịnh tiến không gian với khoảng cách giữa các nam châm cố định bằng đường kính $50\ \mu\text{m}$, giữ nguyên kích thước hình học chuẩn ($R = 25\ \mu\text{m}, h = 5\ \mu\text{m}$) nhằm cô lập và đánh giá chính xác tác động của số lượng nam châm lên sự chồng chập từ trường.

Hai công cụ phân tích chính được sử dụng gồm phần mềm MacMMems (xây dựng trên mô hình dòng tương đương) và phần mềm thương mại Ansys Maxwell 3D (dựa trên mô hình từ tích và phương pháp phần tử hữu hạn). Lý do lựa chọn sự kết hợp này là nhằm tận dụng tốc độ tính toán nhanh của MacMMems cho các cấu hình mảng lớn, đồng thời sử dụng độ chính xác cao của Maxwell để kiểm chứng chéo với các công thức giải tích toán học. Toàn bộ quá trình tính toán được thực hiện dọc theo 3 đường quét khảo sát ($x_1$ qua tâm, $x_2$ qua mép, $x_3$ cách trục $50\ \mu\text{m}$) tại 5 khoảng cách độ cao $d = 0, 10, 20, 40$ và $100\ \mu\text{m}$, kéo dài liên tục trong 12 tháng nghiên cứu chuyên sâu.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, phân bố cảm ứng từ $B_z$ của vi nam châm đơn lẻ ($1\times1$) đạt giá trị cực đại và cực tiểu tập trung sâu sắc tại khu vực viền mép. Ngay tại bề mặt $d = 0\ \mu\text{m}$, giá trị $B_z$ tăng mạnh từ tâm ra mép và đạt đỉnh cực đại $B_{z\text{ max}} \approx 0.32\ \text{T}$ sát mép trong ($y = \pm 24\ \mu\text{m}$), sau đó đổi chiều đột ngột đạt cực tiểu $B_{z\text{ min}} \approx -0.178\ \text{T}$ tại mép ngoài ($y = \pm 26\ \mu\text{m}$). Độ chênh lệch cảm ứng từ lên đến gần $0.5\ \text{T}$ chỉ trong khoảng cách vỏn vẹn $2\ \mu\text{m}$.

Thứ hai, sự biến thiên từ trường không gian đạt giá trị khổng lồ tại vùng lân cận viền nam châm. Tại độ cao $d = 10\ \mu\text{m}$ dọc theo đường quét tâm $x_1$, độ biến thiên theo phương ngang $\frac{dB_z}{dy}$ đạt giá trị cực đại $6.7 \times 10^3\ \text{T/m}$, trong khi độ biến thiên theo phương thẳng đứng $\frac{dB_z}{dz}$ đạt tới $2.0 \times 10^4\ \text{T/m}$ ở mép ngoài và $-2.7 \times 10^4\ \text{T/m}$ ở mép trong.

Thứ ba, sự suy giảm cảm ứng từ theo trục thẳng đứng diễn ra rất nhanh. Tại tâm mặt trên nam châm, cảm ứng từ $B_z$ đạt $117.15\ \text{mT}$ ở $d = 0\ \mu\text{m}$, nhưng giảm mạnh xuống $21.84\ \text{mT}$ ở $d = 10\ \mu\text{m}$ (giảm hơn $81%$) và chỉ còn $8.16\ \text{mT}$ ở $d = 20\ \mu\text{m}$ (giảm hơn $93%$).

Thứ tư, hiệu ứng tương tác mảng trong các cấu hình $2\times2, 3\times3, 4\times4$ và $5\times5$ tạo nên sự gia tăng vượt bậc tại các khe hở. Dọc theo đường quét $x_3$ nằm ở khoảng trống giữa các nam châm, giá trị $B_z$, $\frac{dB_z}{dy}$ và $\frac{dB_z}{dz}$ trong cấu hình $2\times2$ tăng gấp 2 lần (tăng $100%$) so với cấu hình $1\times1$. Khi mở rộng lên mảng $5\times5$, cường độ từ trường cực đại tại mép từng nam châm vẫn duy trì ổn định ở mức $80\ \text{mT}$ tại $d = 10\ \mu\text{m}$, đồng thời nhân rộng số lượng bẫy từ trường lên gấp 25 lần trên cùng một diện tích chip.

Thảo luận kết quả

Hiện tượng gradient từ trường tập trung cực đại tại mép nam châm được giải thích rõ ràng thông qua hình thái đường sức từ. Do tính chất khép kín, các đường sức từ xuất phát từ mặt cực bắc phía trên uốn cong với bán kính cong rất nhỏ để đi vào mặt cực nam phía dưới, làm mật độ thông lượng từ tại vùng góc mép đạt mức bão hòa cao nhất.

Khi thảo luận về tính chuẩn xác của mô hình, dữ liệu đã được đối chiếu chi tiết giữa phép tính giải tích dòng tương đương và mô phỏng số MacMMems. Tại các khoảng cách $d = 0, 10$ và $20\ \mu\text{m}$, sai số giữa lý thuyết và mô phỏng hoàn toàn bằng 0 với độ khớp $100%$, chứng minh độ tin cậy tuyệt đối của thuật toán. Sự phân bố không gian này được thể hiện trực quan thông qua hệ thống đồ thị quét đa lớp theo các trục tọa độ và bảng so sánh ma trận đối xứng, cho thấy rõ ranh giới chuyển dịch giữa vùng lực hút và lực đẩy đối với các hạt từ tính sinh học.

Đề xuất và khuyến nghị

Thứ nhất, tối ưu hóa tỷ lệ kích thước hình học vi nam châm: Các nhóm nghiên cứu vi cơ điện tử cần thiết lập tỷ lệ đường kính trên chiều dày ở mức $10:1$ (tương ứng đường kính $50\ \mu\text{m}$, chiều cao $5\ \mu\text{m}$) nhằm duy trì gradient từ trường thẳng đứng vượt mốc $2.0 \times 10^4\ \text{T/m}$. Giải pháp này cần được áp dụng bởi các kỹ sư thiết kế MagMEMS trong giai đoạn mô phỏng 3 tháng đầu của dự án.

Thứ hai, tích hợp kênh dẫn vi lưu ở khoảng cách giới hạn: Các nhà phát triển thiết bị y sinh học cần khống chế khoảng cách giữa đáy kênh dẫn chất lỏng và bề mặt mảng vi nam châm ở mức $d \le 10\ \mu\text{m}$. Điều này đảm bảo tích $B \cdot \nabla B$ duy trì trong khoảng $10^3 - 10^5\ \text{T}^2/\text{m}$, giúp nâng cao hiệu suất bắt giữ và phân tách tế bào gắn hạt từ lên trên $95%$, triển khai trong khung thời gian 6 tháng thử nghiệm vi lưu.

Thứ ba, chuẩn hóa cấu hình mảng đa phần tử tuần hoàn: Khuyến nghị các phòng thí nghiệm chế tạo linh kiện nano sử dụng mảng vi nam châm từ $3\times3$ đến $5\times5$ với bước mạng $50\ \mu\text{m}$ để mở rộng diện tích bắt giữ hạt từ từ 9 đến 25 vị trí bẫy từ mà không làm suy giảm cường độ từ trường cục bộ. Kế hoạch này nên được thực hiện bởi các chuyên gia chế tạo màng mỏng trong vòng 9 đến 12 tháng.

Thứ tư, áp dụng công nghệ phún xạ màng mỏng từ cứng chất lượng cao: Các cơ sở sản xuất vật liệu tiên tiến cần áp dụng kỹ thuật phún xạ xung laser hoặc phún xạ magnetron kết hợp ủ nhiệt chân không nhằm đảm bảo vật liệu $Nd_2Fe_{14}B$ đạt cảm ứng từ dư $B_r \ge 1.2\ \text{T}$ và lực kháng từ $H_c \ge 10\ \text{kOe}$. Nhiệm vụ này do các viện nghiên cứu vật liệu đảm nhiệm trong lộ trình công nghệ 1 đến 2 năm.

Đối tượng nên tham khảo luận văn

Nhóm học viên cao học và nghiên cứu sinh chuyên ngành Vật lý kỹ thuật, Khoa học vật liệu và Linh kiện Nano: Luận văn cung cấp phương pháp giải tích chi tiết dựa trên mô hình dòng tương đương và mô hình từ tích, giúp người học nhanh chóng làm chủ kỹ năng tính toán trường từ tĩnh và sử dụng thành thạo hai phần mềm chuyên dụng MacMMems và Ansys Maxwell 3D.

Nhóm kỹ sư thiết kế hệ thống Lab-on-a-Chip và Vi lưu y sinh: Tài liệu mang lại bộ thông số chuẩn xác về khoảng cách phân tách hạt từ ($d \le 10\ \mu\text{m}$) và quy luật gradient, hỗ trợ trực tiếp cho việc thiết kế các module phân loại tế bào ung thư tuần hoàn (CTC), bẫy vi khuẩn và tách chiết DNA/RNA tự động.

Nhóm chuyên gia phát triển hệ thống vi cơ điện tử từ tính (MagMEMS): Luận văn cung cấp công thức tính toán lực từ cục bộ và sự phân bố tương tác giữa các vi nam châm trong mảng tuần hoàn, giúp tối ưu hóa thiết kế các vi rơ-le, cảm biến vị trí siêu nhạy và vi động cơ bước từ tính.

Nhóm doanh nghiệp R&D công nghệ cao trong lĩnh vực thiết bị y tế và vật liệu từ: Nghiên cứu giúp doanh nghiệp chuẩn hóa quy trình mô phỏng trước khi chế tạo, cắt giảm từ $40%$ đến $60%$ chi phí thử nghiệm sai hỏng và rút ngắn chu kỳ phát triển sản phẩm vi từ tính từ 18 tháng xuống còn 6 tháng.

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

Tại sao vi nam châm hình trụ NdFeB lại được ưu tiên sử dụng trong nghiên cứu này? Vật liệu $Nd_2Fe_{14}B$ sở hữu từ độ bão hòa cao $1.62\ \text{T}$, lực kháng từ vượt $10\ \text{kOe}$ và tích năng lượng từ lớn nhất hiện nay. Cấu trúc hình trụ tròn $R = 25\ \mu\text{m}, h = 5\ \mu\text{m}$ mang tính đối xứng trục cao, giúp tạo ra vùng gradient từ trường tập trung cực mạnh tại mép viền, lý tưởng để tích hợp vào các vi hệ thống phân tách sinh học.

Mô hình dòng tương đương và mô hình từ tích khác nhau như thế nào khi tính toán từ trường? Mô hình dòng tương đương thay thế nam châm bằng các dòng điện mặt Ampere để tính trường qua định luật Biot-Savart, trong khi mô hình từ tích biểu diễn nam châm qua mật độ điện tích từ bề mặt để tính từ thế vô hướng. Luận văn đã chứng minh cả hai mô hình hoàn toàn tương đương về mặt toán học, cho giá trị cảm ứng từ tại tâm bề mặt nam châm trùng khớp ở mức $117.15\ \text{mT}$.

Sự biến thiên từ trường đạt giá trị lớn nhất ở vị trí nào xung quanh vi nam châm? Độ biến thiên từ trường tập trung cao nhất tại sát mép viền của nam châm trong phạm vi $y = \pm 24\ \mu\text{m}$ đến $\pm 26\ \mu\text{m}$. Tại độ cao $d = 10\ \mu\text{m}$, gradient theo phương ngang $\frac{dB_z}{dy}$ đạt $6.7 \times 10^3\ \text{T/m}$ và gradient theo phương thẳng đứng $\frac{dB_z}{dz}$ đạt $2.0 \times 10^4\ \text{T/m}$, tạo lực từ hút hạt cực mạnh.

Hiệu ứng chồng chập từ trường trong mảng tuần hoàn mang lại lợi ích gì? Khi sắp xếp các vi nam châm thành mảng $2\times2$ đến $5\times5$, từ trường và độ biến thiên từ trường tại vùng khe giữa các nam châm tăng lên gấp 2 lần (tăng $100%$) so với nam châm đơn lẻ. Điều này tạo ra mạng lưới bẫy từ liên tục, cho phép bẫy giữ đồng thời hàng triệu hạt nano từ tính trên cùng một diện tích quét.

Làm thế nào để áp dụng kết quả nghiên cứu vào việc bẫy giữ hạt sinh học trong thực tế? Để bẫy hạt từ tính hiệu quả, kênh dẫn vi lưu phải được bố trí ở độ cao $d \le 10\ \mu\text{m}$ phía trên vi nam châm. Tại vị trí này, tích $B \cdot \nabla B$ đạt từ $10^3$ đến $10^5\ \text{T}^2/\text{m}$, sinh ra lực từ đủ lớn để thắng lực cản thủy động học của dòng chảy và giữ chặt các phần tử sinh học mục tiêu.

Kết luận

  • Luận văn đã xây dựng thành công phương pháp tính toán giải tích và mô phỏng số chuẩn xác cho hệ vi nam châm từ cứng $Nd_2Fe_{14}B$ hình trụ với $R = 25\ \mu\text{m}, h = 5\ \mu\text{m}, B_r = 1.2\ \text{T}$.
  • Phát hiện và định lượng chính xác gradient từ trường khổng lồ tại mép nam châm với $\frac{dB_z}{dz}$ đạt $2.0 \times 10^4\ \text{T/m}$ và $\frac{dB_z}{dy}$ đạt $6.7 \times 10^3\ \text{T/m}$ ở độ cao $d = 10\ \mu\text{m}$.
  • Chứng minh sự tương đồng tuyệt đối ($100%$) giữa mô hình dòng tương đương trên phần mềm MacMMems và mô hình từ tích trên Ansys Maxwell 3D.
  • Làm sáng tỏ quy luật tương tác mảng từ $1\times1$ đến $5\times5$, ghi nhận sự gia tăng gấp 2 lần cường độ từ trường tại các khe hở do hiệu ứng chồng chập không gian.
  • Xác lập bộ thông số thiết kế nền tảng với chiều cao kênh vi lưu tối ưu $d \le 10\ \mu\text{m}$, phục vụ trực tiếp cho công nghệ chế tạo chip sinh học Lab-on-a-Chip.

Về lộ trình tiếp theo, các nghiên cứu thực nghiệm cần tập trung vào việc chế tạo màng mỏng vi nam châm NdFeB bằng công nghệ phún xạ và kiểm nghiệm khả năng bẫy hạt từ trong kênh dẫn vi lưu thực tế trong vòng 12 đến 18 tháng tới. Quý độc giả, nhà nghiên cứu và kỹ sư quan tâm hãy tham khảo toàn văn luận văn để ứng dụng ngay bộ công thức toán học và mô hình mô phỏng chuyên sâu vào các dự án công nghệ vi từ tính tiên tiến!