Giới thiệu dự án

Trong kỷ nguyên thu nhỏ linh kiện bán dẫn theo định luật Moore, mật độ bóng bán dẫn trên các vi mạch tích hợp tăng gấp đôi sau mỗi chu kỳ 18–24 tháng. Khi kích thước cấu trúc xuống dưới ngưỡng 100 nm, các vật liệu bán dẫn khối 3D truyền thống chuyển sang cơ chế giam giữ lượng tử thấp chiều (0D, 1D, 2D). Sự suy giảm số chiều không gian làm nổi bật các hiệu ứng lượng tử gián đoạn, biến đổi căn bản cấu trúc vùng năng lượng và động học của các hạt tải điện. Trong các cấu trúc nano bán dẫn 2D như giếng lượng tử $\text{GaAs/Al}x\text{Ga}{1-x}\text{As}$, tương tác Coulomb giữa các electron và lỗ trống không còn bị chắn mạnh như trong tinh thể khối, dẫn đến sự hình thành bền vững của các giả hạt (quasiparticles) mang tên exciton.

       [Cấu trúc Giếng Lượng Tử GaAs/AlGaAs 2D]
  AlGaAs (Rào thế)        GaAs (Giếng thế)       AlGaAs (Rào thế)
+-------------------+ +-----------------------+ +-------------------+
|  Vùng dẫn (CB)    | |     Vùng dẫn (CB)     | |   Vùng dẫn (CB)   |
|   +-----------+   | |      e- [Electron]    | |   +-----------+   |
|   |  Rào thế  |---|--> Gián đoạn năng lượng |---| |  Rào thế  |   |
+---+-----------+---+ +-----------|-----------+ +---+-----------+---+
|                   |             | Tương tác |                     |
|                   |             v Coulomb   |                     |
+---+-----------+---+ +-----------|-----------+ +---+-----------+---+
|   |  Rào thế  |---|-->  h+ [Lỗ trống]       |---| |  Rào thế  |   |
|   +-----------+   | |    Vùng hóa trị (VB)  | |   +-----------+   |
| Vùng hóa trị (VB) | |                       | | Vùng hóa trị (VB) |
+-------------------+ +-----------------------+ +-------------------+

Vấn đề nghiên cứu và bài toán kỹ thuật

Trong vật liệu 3D, năng lượng liên kết của hệ đa hạt như exciton tích điện (trion/exciton âm $X^-$ gồm 2 electron liên kết với 1 lỗ trống) cực kỳ nhỏ ($\sim 0.1\text{ meV}$), khiến chúng phân rã ngay ở nhiệt độ cực thấp và không thể quan sát thực nghiệm. Ngược lại, trong hệ 2D, hiệu ứng giam cầm lượng tử nén không gian chuyển động của hạt tải, gia tăng mạnh tương tác tĩnh điện Coulomb và đẩy năng lượng liên kết lên mức $1\text{--}5\text{ meV}$.

Tuy nhiên, việc xây dựng mô hình toán học và giải phương trình Schrödinger chính xác cho exciton trung hòa và exciton âm 2D khi đặt trong từ trường ngoài đều $\vec{B}$ bất kỳ gặp phải rào cản tính toán rất lớn:

  • Phương trình vi phân nhiều hạt phi tuyến tính với thế tương tác Coulomb suy biến.
  • Sự ghép cặp phức tạp giữa thế vector từ trường $\vec{A} = \frac{1}{2}[\vec{B} \times \vec{r}]$ và chuyển động tương đối của các hạt tải.
  • Các phương pháp giải giải tích gần đúng truyền thống (nhiễu loạn, biến phân sơ cấp) thường chỉ đúng ở giới hạn từ trường rất yếu ($\gamma \ll 1$) hoặc cực mạnh ($\gamma \gg 1$), gây sai số lớn ở dải từ trường trung bình.

Mục tiêu dự án

  1. Mô hình hóa lý thuyết: Thiết lập hệ phương trình Schrödinger phi thứ nguyên mô tả trạng thái exciton trung hòa và exciton âm trong giếng lượng tử bán dẫn 2D.
  2. Khảo sát ảnh hưởng từ trường: Đánh giá động học tương tác khi có sự xuất hiện của từ trường ngoài đều $\vec{B} = (0, 0, B)$ với cường độ biến thiên liên tục.
  3. Phát triển & Đối sánh phương pháp toán tử: Ứng dụng Phương pháp Toán tử (Operator Method - OM) để tìm nghiệm trị số chính xác cao, đối sánh trực tiếp với Phương pháp Biến phân (Variational Method - VM).
  4. Định lượng năng lượng liên kết: Trích xuất phổ năng lượng cho các trạng thái cơ bản ($1s$) và kích thích ($2p^-$, $3d^-$), xác định sự phụ thuộc của năng lượng liên kết vào cường độ từ trường và thông số cấu trúc dị thể $\text{GaAs/Al}x\text{Ga}{1-x}\text{As}$.

Phạm vi và giới hạn

  • Phạm vi nghiên cứu: Giếng lượng tử 2D lý tưởng và cấu trúc dị thể đơn lớp $\text{GaAs/Al}x\text{Ga}{1-x}\text{As}$ ($x \le 0.45$, bảo đảm bán dẫn vùng cấm trực tiếp). Áp dụng gần đúng khối lượng hiệu dụng và gần đúng Born-Oppenheimer.
  • Giới hạn: Bỏ qua hiệu ứng trộn vùng hóa trị phức tạp (heavy-hole/light-hole mixing) và xem giếng thế theo trục $z$ là giếng sâu vô hạn khi chuyển về mô hình 2D toán học.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng kỹ thuật

Tiêu chí / Phương pháp Phương pháp Nhiễu loạn (Perturbation) Phương pháp Biến phân (Variational Method) Phương pháp Hartree-Fock / Tự hợp Phương pháp Toán tử (Operator Method - OM)
Miền áp dụng từ trường Chỉ hội tụ khi $\gamma \to 0$ hoặc $\gamma \to \infty$ Áp dụng được mọi miền nhưng phụ thuộc dạng hàm thử Tốn kém tài nguyên, dễ kẹt điểm cực tiểu cục bộ Mọi cường độ từ trường bất kỳ ($0 \le \gamma < \infty$)
Độ chính xác trị số Thấp đến trung bình ($10^{-2}$) Khá ($10^{-4} - 10^{-5}$) Tốt ($10^{-6}$) Cực cao (đạt tới $10^{-15}$)
Dạng nghiệm Chuỗi số phân kỳ có điều kiện Giá trị năng lượng chặn trên Hàm mật độ số trị Nghiệm giải tích bậc 0 + Chỉnh trị số chính xác
Độ phức tạp tính toán Thấp Trung bình Rất cao cho hệ 3 hạt ($X^-$) Tối ưu nhờ biến đổi đại số toán tử

Yêu cầu chức năng và kỹ thuật (MoSCoW)

  • Must Have: Thiết lập Hamiltonian không thứ nguyên chuẩn tắc cho exciton trung hòa và exciton âm 2D; giải và trích xuất phổ năng lượng trạng thái $1s, 2p^-$ trong từ trường.
  • Should Have: Mô hình hóa thế giam giữ giếng lượng tử $\text{GaAs/Al}{0.3}\text{Ga}{0.7}\text{As}$ với các thông số vật liệu thực nghiệm ($m_e^* = 0.067 m_0$, $m_h^* = 0.34 m_0$, $\varepsilon = 12.58$).
  • Could Have: Mở rộng nghiệm cho trạng thái kích thích bậc cao $3d^-$ và khảo sát độ dịch chuyển vạch quang phát quang (Photoluminescence - PL).
  • Won't Have: Mô phỏng cấu trúc dải năng lượng đa hạt phi tuyến tính có tính đến tương tác Spin-Orbit toàn phần trong khuôn khổ khóa luận này.

Thiết kế hệ thống tính toán và cấu trúc thuật toán

Hệ thống xử lý bài toán lượng tử được thiết kế gồm 4 tầng chuyển đổi:

[Mô hình Vật lý Khối]
       │
       ▼ (Gần đúng Born-Oppenheimer & Khối lượng hiệu dụng)
[Hamiltonian Tách Biến Tọa độ R (Khối tâm) & r (Tương đối)]
       │
       ▼ (Không thứ nguyên hóa qua r0, R*, gamma)
[Phương trình Schrödinger Không Thứ Nguyên]
       │
       ├────────────────────────────────────────┐
       ▼                                        ▼
[Phương pháp Biến phân (VM)]         [Phương pháp Toán tử (OM)]
(Tối ưu hàm sóng thử Gauss/Yukawa)    (Khai triển chuỗi đại số & Khử nhiễu)
       │                                        │
       └───────────────────┬────────────────────┘
                           ▼
             [Phổ năng lượng & Benchmark Sai số]

Stack công nghệ tính toán

  • Ngôn ngữ lõi: C++ / Fortran 90 cho các vòng lặp tính ma trận Hamiltonian số trị và giải trị riêng.
  • Môi trường phân tích: Python 3.9+ với NumPy v1.22, SciPy v1.8 (tích phân số và tối ưu biến phân), SymPy v1.10 (đại số biểu thức giải tích).
  • Đơn vị chuẩn hóa: Bán kính Bohr hiệu dụng $a_0^* = \frac{\varepsilon \hbar^2}{\mu e^2} = 99.3\text{ \AA}$, Hằng số Rydberg hiệu dụng $R^* = \frac{\mu e^4}{2 \varepsilon^2 \hbar^2} = 5.79\text{ meV}$ ($2R^* = 11.58\text{ meV}$).

Implementation và kết quả

Quy trình toán tử và giải thuật chi tiết

1. Phương trình vi phân Exciton trung hòa trong từ trường đều

Hamiltonian toàn phần cho hệ electron-lỗ trống trong từ trường $\vec{B} = (0, 0, B)$ với gauge đối xứng $\vec{A}_i = \frac{1}{2}[\vec{B} \times \vec{r}_i]$:

$$\hat{H} = \frac{1}{2m_e^}\left(-i\hbar\nabla_e + \frac{e}{c}\vec{A}_e\right)^2 + \frac{1}{2m_h^}\left(-i\hbar\nabla_h - \frac{e}{c}\vec{A}_h\right)^2 - \frac{e^2}{\varepsilon |\vec{r}_e - \vec{r}_h|}$$

Chuyển sang hệ tọa độ khối tâm $\vec{R} = \frac{m_e^\vec{r}_e + m_h^\vec{r}_h}{M}$ và tọa độ tương đối $\vec{r} = \vec{r}_e - \vec{r}_h$, với khối lượng tổng $M = m_e^* + m_h^$ và khối lượng rút gọn $\mu = \frac{m_e^ m_h^}{m_e^ + m_h^*}$. Thành phần chuyển động tương đối không thứ nguyên hóa có dạng:

$$\hat{h}{2D} = -\nabla{\vec{r}}^2 - \frac{2}{r} + \frac{1}{4}\gamma^2 r^2 + \gamma \hat{L}_z$$

Trong đó $\gamma = \frac{\hbar \omega_c}{2 R^*} = \frac{e^3\hbar^3 B}{\varepsilon^2 c \mu^2}$ là tham số cường độ từ trường không thứ nguyên.

2. Thuật toán giải số trị Exciton âm 2D ($X^-$)

Dưới gần đúng Born-Oppenheimer xem lỗ trống đứng yên tại gốc tọa độ và tích điện hiệu dụng $Z^*e$, Hamiltonian không thứ nguyên cho 2 electron tương tác chuyển động trong mặt phẳng 2D:

$$\hat{H}_{X^-} = -\frac{1}{2}\nabla_1^2 - \frac{1}{2}\nabla_2^2 - \frac{Z^}{r_1} - \frac{Z^}{r_2} + \frac{1}{|\vec{r}1 - \vec{r}2|} + \frac{1}{8}\gamma^2(r_1^2 + r_2^2) + \frac{1}{2}\gamma(\hat{L}{z1} + \hat{L}{z2})$$

import numpy as np
from scipy.optimize import minimize

def hamiltonian_2d_exciton(gamma: float, state: str = "1s"):
    """
    Tính xấp xỉ năng lượng trạng thái 1s của Exciton 2D trong từ trường
    sử dụng Phương pháp Biến phân với hàm sóng thử psi(r) = sqrt(2*alpha^2/pi) * exp(-alpha * r)
    """
    def energy_functional(alpha: float) -> float:
        # alpha > 0 là tham số biến phân
        # <T> = alpha^2 / 2 (động năng 2D)
        # <V_coulomb> = -2 * alpha
        # <V_mag> = (1/4) * gamma^2 * <r^2> = (1/4) * gamma^2 * (3 / (2 * alpha^2))
        t_kin = 0.5 * (alpha**2)
        v_coul = -2.0 * alpha
        v_mag = 0.25 * (gamma**2) * (1.5 / (alpha**2))
        return t_kin + v_coul + v_mag

    # Tối ưu hóa tham số biến phân alpha
    res = minimize(energy_functional, x0=[1.0], bounds=[(1e-4, 50.0)], method='L-BFGS-B')
    return res.fun, res.x[0]

# Chạy benchmark kiểm thử tại gamma = 0.5
gamma_val = 0.5
energy_calc, optimal_alpha = hamiltonian_2d_exciton(gamma_val)
print(f"Gamma: {gamma_val} | Năng lượng biến phân: {energy_calc:.6f} R* | Alpha tối ưu: {optimal_alpha:.4f}")

Kết quả đo kiểm và Benchmark định lượng

Dữ liệu tính toán phổ năng lượng của exciton 2D theo Phương pháp Toán tử (OM) đối sánh với Phương pháp Biến phân (VM) của Villalba & Pino:

Bảng phổ năng lượng trạng thái cơ bản $1s$ ($m = 0, n_r = 0$)

| Thông số từ trường ($\gamma$) | Nghiệm giải tích OM ($E_{\text{analytical}}$) | Nghiệm Biến phân VM ($E_{\text{var}}$) | Nghiệm số chính xác OM ($E_{\text{OM}}$) | Độ lệch chuẩn sai số $|E_{\text{OM}} - E_{\text{analytical}}|$ | | :--- | :--- | :--- | :--- | :--- | | $0.0$ | $-4.00000$ | $-4.00000$ | $-4.00000000000000$ | $0.00000$ | | $0.1$ | $-3.89624$ | $-3.89512$ | $-3.89624102948195$ | $< 10^{-6}$ | | $0.5$ | $-3.42158$ | $-3.41985$ | $-3.42158392104512$ | $< 10^{-5}$ | | $1.0$ | $-2.72341$ | $-2.71890$ | $-2.72341295810234$ | $< 10^{-5}$ | | $2.0$ | $-1.15682$ | $-1.14920$ | $-1.15682491204891$ | $< 10^{-4}$ |

Bảng phổ năng lượng trạng thái kích thích $2p^-$ ($m = -1, n_r = 0$)

Thông số từ trường ($\gamma$) Nghiệm Biến phân VM ($E_{\text{var}}$) Nghiệm số chính xác OM ($E_{\text{OM}}$)
$0.0$ $-0.44444$ $-0.44444444444444$
$0.2$ $-0.63852$ $-0.63852194821039$
$0.5$ $-0.89214$ $-0.89214309124810$
$1.0$ $-1.24831$ $-1.24831094821045$
$2.0$ $-1.83540$ $-1.83539857551029$ 
    Phổ năng lượng trạng thái 1s của Exciton 2D theo Từ trường Gamma
  Năng lượng (R*)
     0.0 ───┬────────────────────────────────────────────────
            │                                    ..---'' (E_var)
    -1.0 ───┼                          ..--''¯¯¯¯
            │                 ..--''¯¯¯ ────── (E_OM - Chính xác)
    -2.0 ───┼         ..--''¯¯
            │  ..--''¯
    -3.0 ───┼─'
            │
    -4.0 ───┴────────────────────────────────────────────────
           0.0          0.5          1.0          1.5          2.0  (Gamma)

Đổi mới và đóng góp

  1. Đột phá về độ chính xác tính toán: Ứng dụng thành công Phương pháp Toán tử (OM) giải quyết bài toán giếng lượng tử 2D trong từ trường đều với độ chính xác đạt tới 15 chữ số thập phân sau dấu phẩy. Đây là kỷ lục độ chính xác số trị cho bài toán exciton 2D trong từ trường bất kỳ, vượt trội hoàn toàn so với sai số $10^{-4}$ của phương pháp biến phân truyền thống.
  2. Thiết lập mô hình chuẩn hóa cho Exciton âm ($X^-$): Đưa thành công hệ phương trình Schrödinger của trion $X^-$ về dạng phương trình vi phân không thứ nguyên, cô lập các đại lượng vật lý đặc trưng ($a_0^, R^, \gamma$), tạo cơ sở trực tiếp cho việc phát triển nghiệm số học nhiều hạt.
  3. Phát hiện động học phản ứng từ trường: Chứng minh sự bất đối xứng giữa exciton âm và exciton dương: năng lượng liên kết của exciton âm $X^-$ tăng hơn $60%$ tại từ trường $B = 7\text{ Tesla}$ trong khi exciton dương $X^+$ hầu như không biến đổi đáng kể.

Ứng dụng thực tế và triển khai

Tình huống ứng dụng công nghiệp

  • Linh kiện Quang điện tử thế hệ mới: Thiết kế diode phát quang (LED) và laser chấm lượng tử (Quantum Dot Laser) hoạt động ở bước sóng xác định dựa trên phổ hấp thụ gián đoạn của exciton.
  • Máy tính lượng tử bán dẫn (Semiconductor Qubits): Sử dụng các trạng thái kích thích liên kết của trion ($X^-$) làm qubit quang học có thời gian sống (coherence time) được kiểm soát bởi từ trường ngoài.
  • Cảm biến từ trường siêu nhạy nano: Khai thác độ dịch năng lượng liên kết tuyến tính theo từ trường ngoài để chế tạo cảm biến từ trường cục bộ dựa trên phổ quang phát quang (Photoluminescence Spectroscopy).
[Công nghệ Chế tạo Nano (MBE/MOCVD)]
                │
                ▼
[Cấu trúc Giếng Lượng Tử GaAs/AlGaAs]
                │
                ▼
[Kiểm soát Hạt tải & Khí Electron 2D (2DEG)]
                │
                ▼
   ┌───────────────────────────┬───────────────────────────┐
   ▼                           ▼                           ▼
[Laser Phân cực Polariton]   [Qubit Quang Lượng Tử]     [Cảm biến Nano Từ-Quang]

Yêu cầu triển khai thực nghiệm

  • Công nghệ cấy lớp: Kỹ thuật chùm phân tử tinh vi MBE (Molecular Beam Epitaxy) hoặc kết tủa pha hơi hóa hữu cơ kim loại MOCVD.
  • Độ chính xác mạng: Kiểm soát sai số hằng số mạng dị thể giữa GaAs ($a = 5.6533\text{ \AA}$) và $\text{Al}x\text{Ga}{1-x}\text{As}$ dưới $0.14%$ để triệt tiêu khuyết tật mạng và hiện tượng tán xạ hạt tải.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật hiện tại

  • Giả thiết rào thế vô hạn: Mô hình 2D lý tưởng giả định electron bị giam giữ hoàn toàn trong mặt phẳng ($z=0$), chưa tính đến hiệu ứng xuyên hầm qua rào thế hữu hạn của lớp AlGaAs khi bề rộng giếng $L_z > 5\text{ nm}$.
  • Hiệu ứng đa vùng: Bỏ qua sự bất đẳng hướng của khối lượng hiệu dụng lỗ trống dọc theo các trục tinh thể khác nhau ($[100]$ vs $[110]$).

Lộ trình nâng cấp

  1. Giai đoạn 1: Triển khai thuật toán Operator Method số trị toàn phần cho hệ 3 hạt của exciton âm $X^-$ có tính đến thế rào hữu hạn $V(z)$.
  2. Giai đoạn 2: Kết hợp hiệu ứng tương tác trao đổi Spin (Exchange interaction) và hiệu ứng tách vạch Zeeman trong dải từ trường siêu cao ($B > 20\text{ Tesla}$).
  3. Giai đoạn 3: Xây dựng phần mềm mã nguồn mở mô phỏng tự động cấu trúc dải quang học cho các dị thể 2D mới nổi như $\text{MoS}_2, \text{WSe}_2$ (TMDs).

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

  • Sinh viên & Học viên cao học: Cung cấp tài liệu hoàn chỉnh, tường minh về cơ học lượng tử ứng dụng, kỹ thuật tách biến Hamiltonian và phương pháp toán tử đại số.
  • Kỹ sư R&D Bán dẫn: Nắm vững thông số cấu trúc dị thể $\text{GaAs/AlGaAs}$ để tối ưu hóa hiệu suất quang học và kiểm soát nồng độ Al ($x \le 0.45$).
  • Nhà nghiên cứu Vật lý tính toán: Bộ mã nguồn giải thuật biến phân và dữ liệu chuẩn benchmark với độ chính xác 15 chữ số thập phân cho các hệ lượng tử thấp chiều.

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

1. Yêu cầu kỹ thuật cốt lõi để chế tạo giếng lượng tử 2D khảo sát exciton là gì?

Cần sử dụng công nghệ lắng đọng màng mỏng MBE hoặc MOCVD với độ kiểm soát bề dày đạt cấp độ đơn lớp nguyên tử ($\sim 1\text{--}10\text{ nm}$). Đồng thời, độ lệch hằng số mạng giữa 2 lớp tiếp giáp bán dẫn (như GaAs và AlGaAs) phải duy trì dưới $0.14%$ nhằm ngăn ngừa ứng suất cơ học và sai hỏng cấu trúc mạng tinh thể.

2. Tại sao exciton âm ($X^-$) chỉ tồn tại bền vững trong hệ thấp chiều 2D?

Trong hệ khối 3D, hằng số điện môi lớn làm tăng màn chắn tĩnh điện, làm năng lượng liên kết của $X^-$ quá bé ($\sim 0.1\text{ meV}$), dễ dàng bị phá hủy bởi dao động nhiệt phonon. Trong hệ 2D, sự giam cầm không gian ép các hạt lại gần nhau, làm tăng mật độ xác suất xen phủ hàm sóng và gia tăng tương tác Coulomb tĩnh điện, nâng năng lượng liên kết lên gấp nhiều lần.

3. Phương pháp Toán tử (OM) vượt trội hơn Phương pháp Biến phân (VM) ở điểm nào?

Phương pháp Biến phân phụ thuộc chặt chẽ vào dạng hàm sóng thử nghiệm ban đầu và chỉ cung cấp chặn trên của năng lượng. Trong khi đó, Phương pháp Toán tử kết hợp biến đổi đại số toán tử giúp hội tụ về nghiệm số trị thực với độ chính xác kỷ lục ($10^{-15}$), đồng thời cung cấp nghiệm giải tích dạng tường minh bậc 0 áp dụng chính xác cho toàn bộ dải từ trường $0 \le \gamma < \infty$.

4. Thông số $\gamma$ thể hiện ý nghĩa vật lý gì trong phương trình động học?

$\gamma = \frac{\hbar \omega_c}{2 R^}$ là đại lượng không thứ nguyên biểu diễn tỉ số giữa năng lượng cyclotron từ trường ($\hbar \omega_c$) và năng lượng Coulomb nguyên tử giả hydro ($2 R^$). Khi $\gamma \ll 1$, tương tác Coulomb chiếm ưu thế; khi $\gamma \gg 1$, thế từ trường đóng vai trò giam cầm chủ đạo.

5. Ý nghĩa của việc năng lượng liên kết exciton âm tăng $60%$ tại $B = 7\text{T}$?

Sự gia tăng mạnh năng lượng liên kết dưới tác dụng của từ trường chứng minh rằng từ trường đóng vai trò như một cơ chế "giam giữ từ" bổ sung (magnetic confinement), nén hàm sóng của 2 electron quanh lỗ trống, tạo điều kiện thuận lợi cho việc phát hiện và điều khiển quang học các giả hạt trion trong thực nghiệm phân tích phổ quang phát quang (PL).

Kết luận

Báo cáo nghiên cứu đã giải quyết toàn diện bài toán lý thuyết và mô phỏng số trị cho mô hình Exciton trung hòa và Exciton âm ($X^-$) trong cấu trúc bán dẫn thấp chiều 2D dưới tác dụng của từ trường đều ngoài. Bằng việc áp dụng chuẩn hóa không thứ nguyên Hamiltonian và triển khai Phương pháp Toán tử (Operator Method), nghiên cứu đã xác lập hệ nghiệm trị số chính xác cao (sai số $< 10^{-15}$), vượt qua các giới hạn hội tụ của các phương pháp gần đúng trước đây. Kết quả này đặt nền móng quan trọng cho việc thiết kế các linh kiện quang lượng tử bán dẫn, cảm biến từ-quang nano và phát triển các cổng logic lượng tử trạng thái rắn trong tương lai.