Tổng quan luận án

Luận án tiến sĩ vật lý của tác giả Nguyễn Thị Lâm Quỳnh với đề tài "Ảnh hưởng của sự lượng tử hoá do giảm kích thước lên hiệu ứng Ettingshausen và hiệu ứng Peltier trong siêu mạng và hố lượng tử" thuộc chuyên ngành Vật lí lí thuyết và vật lí toán (mã số: 9440130.01), được thực hiện tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội dưới sự hướng dẫn khoa học của PGS. Nguyễn Bá Đức và GS. Nguyễn Quang Báu (năm 2023).

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

Kể từ khi kỹ thuật epitaxy chùm phân tử (Molecular Beam Epitaxy - MBE) ra đời vào những năm 60 của thế kỷ XX, khả năng chế tạo các màng mỏng bán dẫn chất lượng cao đã tạo tiền đề cho sự hình thành của các cấu trúc bán dẫn thấp chiều có kích thước nano. Trong các hệ này, chuyển động của hạt tải bị giới hạn trong các vùng kích thước cỡ bước sóng De Broglie, làm gián đoạn phổ năng lượng theo các phương giam giữ. Khi hệ thấp chiều chịu tác động đồng thời của các trường ngoài như điện trường không đổi, từ trường và sóng điện từ mạnh, các định luật bảo toàn năng xung lượng bị thay đổi, dẫn đến sự xuất hiện của các hiệu ứng quang phi tuyến, hiệu ứng cao tần và các hiệu ứng động nhiệt điện dị thường.

Nghiên cứu các hiệu ứng từ - nhiệt - điện, tiêu biểu là hiệu ứng Ettingshausen (EE) và hiệu ứng Peltier (PE), có ý nghĩa khoa học và ứng dụng thực tiễn trong công nghệ làm lạnh trạng thái rắn không cần chi tiết chuyển động cơ học hay chất lỏng tuần hoàn, cũng như phát triển vật liệu "nam châm vĩnh cửu nhiệt điện" từ các hợp kim loại $\text{SmCo}_5$. Tuy nhiên, phương pháp phương trình động Boltzmann cổ điển chỉ áp dụng được trong vùng nhiệt độ cao. Phương pháp phương trình động lượng tử cho phép mở rộng nghiên cứu trên toàn dải nhiệt độ. Trước thời điểm luận án thực hiện, các công trình về EE và PE trong hệ bán dẫn thấp chiều hai chiều (2D) và một chiều (1D) chỉ mới tính đến sự giam cầm của điện tử mà chưa tính đến sự giam cầm của phonon. Ảnh hưởng của sự giam cầm phonon lên EE và PE trong các hệ 2D dưới tác dụng của trường ngoài là khoảng trống nghiên cứu mà luận án hướng tới giải quyết.

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

Mục tiêu tổng quát của luận án là nghiên cứu ảnh hưởng của sự giam cầm phonon lên hai hiệu ứng từ - nhiệt - điện điển hình là hiệu ứng Ettingshausen (EE) và hiệu ứng Peltier (PE) trong các hệ bán dẫn hai chiều (2D), bao gồm hố lượng tử (QW), siêu mạng bán dẫn pha tạp (DSS) và siêu mạng bán dẫn hợp phần (CSS) khi đặt trong điện trường không đổi, từ trường và sóng điện từ mạnh.

Các nhiệm vụ nghiên cứu cụ thể được xác định theo 4 nội dung:

  1. Thiết lập phương trình động lượng tử cho hàm phân bố của điện tử trong các hệ bán dẫn 2D (QW, DSS, CSS) khi có mặt điện trường không đổi, từ trường, sóng điện từ mạnh và có xét đến sự giam cầm của cả điện tử và phonon.
  2. Giải phương trình động lượng tử, tìm hàm phân bố không cân bằng của điện tử, tính toán mật độ dòng toàn phần và mật độ thông lượng nhiệt; từ đó xác lập biểu thức giải tích cho các ten-xơ động, hệ số Ettingshausen (EC) và hệ số Peltier (PC).
  3. Tính toán số và vẽ đồ thị sự phụ thuộc của EC và PC vào các đại lượng đặc trưng của trường ngoài cùng các tham số cấu trúc vật liệu đối với các hệ 2D tiêu biểu.
  4. Đánh giá các kết quả thu được, so sánh với kết quả trong bán dẫn khối và trong các hệ 2D khi chưa kể tới sự giam cầm của phonon nhằm làm rõ vai trò của phonon giam cầm.

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

  • Đối tượng nghiên cứu: Hiệu ứng Ettingshausen và hiệu ứng Peltier thông qua các đại lượng đặc trưng là hệ số Ettingshausen (EC, ký hiệu $P$) và hệ số Peltier (PC, ký hiệu $\Pi$).
  • Phạm vi cấu trúc: Các cấu trúc bán dẫn thấp chiều hai chiều (2D) gồm hố lượng tử parabol (PQW), siêu mạng bán dẫn pha tạp (DSS) và siêu mạng bán dẫn hợp phần (CSS). Cấu trúc tính toán số đại diện là vật liệu bán dẫn $\text{GaAs/AlGaAs}$.
  • Phạm vi tương tác và gần đúng:
    • Xét các quá trình tán xạ hoặc hấp thụ không quá một photon ($\ell = 0, \pm 1$).
    • Tương tác điện tử giam cầm - phonon giam cầm được giả định là cơ chế tương tác trội.
    • Xét hai cơ chế tán xạ chính: tán xạ điện tử giam cầm - phonon âm giam cầm (CAP) và tán xạ điện tử giam cầm - phonon quang giam cầm (COP).

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

Lịch sử nghiên cứu các hiệu ứng và công nghệ liên quan

Luận án điểm lại các mốc phát triển lý thuyết và thực nghiệm:

  • Kỹ thuật màng mỏng: John R. Arthur và Alfred Y. Cho phát triển kỹ thuật epitaxy chùm phân tử (MBE) vào những năm 60 của thế kỷ XX, mở ra khả năng chế tạo các màng dị cấu trúc siêu mỏng chất lượng cao.
  • Hiệu ứng Ettingshausen: Albert von Ettingshausen (1850–1932) và Walther Hermann Nernst (1864–1941) phát hiện khi nghiên cứu hiệu ứng Hall trong Bismuth; hạt tải chịu tác dụng của lực từ làm xuất hiện gradient nhiệt độ ngang $\nabla T$ vuông góc với dòng điện $\vec{j}$ và từ trường $\vec{B}$.
  • Hiệu ứng Peltier: Jean Charles Athanase Peltier (1785–1845) phát hiện năm 1834 khi cho dòng điện chạy qua lớp tiếp xúc Đồng - Bismuth, làm xuất hiện sự hấp thụ hoặc tỏa nhiệt tại mặt tiếp xúc.
  • Cấu trúc hố lượng tử: Herbert Kroemer đề xuất khái niệm hố lượng tử (QW) vào năm 1963 dựa trên sự gián đoạn cực tiểu vùng dẫn giữa các lớp bán dẫn khác nhau.

Các phương pháp nghiên cứu lý thuyết trong tài liệu

  • Phương pháp phương trình động Boltzmann cổ điển được sử dụng trong các nghiên cứu truyền thống về hệ số nhiệt điện, nhưng bị giới hạn nghiêm ngặt ở nhiệt độ cao.
  • Phương pháp hàm Green và phương pháp tích phân phiếm hàm được phát triển để xử lý các bài toán tương tác lượng tử phức tạp.
  • Phương pháp phương trình động lượng tử (Quantum Kinetic Equation) được áp dụng hiệu quả để thu được kết quả giải tích đúng trên toàn bộ dải nhiệt độ, đặc biệt ở vùng nhiệt độ thấp và có trường ngoài mạnh.

Vị trí và đóng góp của luận án so với y văn

Các công trình khoa học công bố trước đó nghiên cứu về EE và PE trong hệ hai chiều (2D) và một chiều (1D) mới chỉ xét đến sự lượng tử hóa chuyển động của điện tử (electron confinement) và xem phonon như trong môi trường khối (un-confined phonon). Đồng thời, một số nghiên cứu khác chỉ ra rằng sự giam cầm phonon (phonon confinement) có ảnh hưởng đáng kể đến hiệu ứng radio điện và hiệu ứng Hall lượng tử. Luận án của Nguyễn Thị Lâm Quỳnh định vị tại giao điểm này: xây dựng hoàn chỉnh bức tranh giải tích và số trị cho EE và PE khi tính đồng thời cả sự giam cầm của điện tử và phonon (CAP và COP) dưới tác dụng đồng thời của điện trường không đổi, từ trường và sóng điện từ mạnh.


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

Cơ sở lý thuyết và hình thức luận

Luận án sử dụng hình thức luận lượng tử hóa lần thứ hai và phương trình chuyển động Heisenberg để mô tả hệ điện tử - phonon trong trường ngoài.

Hamiltonian toàn phần của hệ trong trường ngoài có dạng: $$H = H_e + H_{ph} + H_{e-ph} + H_{field}$$

Trong đó:

  • $a_{\vec{k}}^+, a_{\vec{k}}$ là toán tử sinh và hủy điện tử ở trạng thái sóng $\vec{k}$.
  • $b_{\vec{q}}^+, b_{\vec{q}}$ là toán tử sinh và hủy phonon với véc-tơ sóng $\vec{q}$.
  • Thế vô hướng của điện trường và từ trường: $\phi(\vec{q}) = (2\pi i)^3 \left[ e\vec{E} + \omega_c [\vec{q}, \vec{h}] \right] \frac{\partial}{\partial \vec{q}} \delta(\vec{q})$, với $\omega_c = \frac{eB}{m_e c}$ là tần số cyclotron và $\vec{h} = \frac{\vec{B}}{B}$.
  • Thế véc-tơ $\vec{A}(t)$ thỏa mãn $-\frac{1}{c}\frac{\partial \vec{A}(t)}{\partial t} = \left(e\vec{E} - \frac{\varepsilon - \varepsilon_F}{T}\nabla T\right)\sin\omega t$.

Hằng số tương tác điện tử - phonon trong bán dẫn khối:

  • Tán xạ điện tử - phonon âm (Acoustic phonon): $$|C_{\vec{q}}^A|^2 = \frac{\hbar \varsigma^2 q^2}{2\rho v_s V_0}$$
  • Tán xạ điện tử - phonon quang (Optical phonon): $$|C_{\vec{q}}^O|^2 = \frac{2\pi e^2 \hbar\omega_0}{\varepsilon_0 V_0 q^2}\left(\frac{1}{\chi_\infty} - \frac{1}{\chi_0}\right)$$
Ký hiệu đại lượng Tên gọi đại lượng vật lý Đơn vị / Thứ nguyên trong luận án
$m_e$ Khối lượng hiệu dụng của điện tử $\text{kg}$ (tính theo $m_0$)
$e$ Điện tích hiệu dụng của điện tử $\text{C}$ (tính theo $e_0$)
$\vec{E}$ / $\vec{E}_0, \omega$ Điện trường không đổi / Biên độ, tần số sóng điện từ $\text{V/m}$ / $\text{rad/s}$ (Hz)
$\vec{B}$ / $\omega_c$ Từ trường ngoài / Tần số cyclotron $\text{T}$ / $\text{rad/s}$
$\varepsilon_F$ Năng lượng Fermi của điện tử $\text{meV}$
$\tau(\varepsilon)$ Thời gian phục hồi năng lượng của điện tử $\text{s}$
$v_s$ / $\rho$ Vận tốc sóng âm / Mật độ khối lượng $\text{m/s}$ / $\text{g/cm}^3$
$\varsigma$ Hằng số thế biến dạng $\text{eV}$
$\chi_0, \chi_\infty$ Độ điện thẩm tĩnh và độ điện thẩm cao tần Không thứ nguyên
$K_L$ Hệ số dẫn nhiệt của mạng tinh thể $\text{W/(m}\cdot\text{K)}$
$N, n$ Chỉ số mức Landau và mức con giam cầm điện tử Số nguyên
$m$ Chỉ số lượng tử giam cầm của phonon Số nguyên

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

  1. Phương pháp giải tích lượng tử:

    • Sử dụng phương trình chuyển động Heisenberg đối với hàm phân bố hạt tải $n_{\vec{k}}(t) = \langle a_{\vec{k}}^+ a_{\vec{k}} \rangle_t$ và hàm tương quan điện tử - phonon $F_{\vec{k}_1, \vec{k}2, \vec{q}}(t) = \langle a{\vec{k}1}^+ a{\vec{k}2} b{\vec{q}} \rangle_t$.
    • Giải phương trình vi phân tuyến tính cấp một không thuần nhất bằng phương pháp biến thiên hằng số Lagrange.
    • Áp dụng khai triển hàm sóng Bessel $e^{\pm i\alpha \sin\theta} = \sum_{\mu} J_\mu(\alpha) e^{\pm i\mu\theta}$ và giữ lại gần đúng bậc một qua cường độ trường ngoài ($\ell = 0, \pm 1$).
    • Sử dụng gần đúng khí điện tử suy biến hoàn toàn: $f_0(\varepsilon_{\vec{k}}) = n_0 \theta(\varepsilon_F - \varepsilon_{\vec{k}})$.
    • Tính toán mật độ dòng toàn phần $\vec{J} = \int \vec{R}(\varepsilon) d\varepsilon = \hat{\sigma}\vec{E} + \hat{\beta}\nabla T$ và mật độ thông lượng nhiệt $\vec{Q} = \frac{1}{e}\int (\varepsilon - \varepsilon_F)\vec{R}(\varepsilon) d\varepsilon - K_L \nabla T = \hat{\gamma}\vec{E} + \hat{\xi}\nabla T$.
    • Thiết lập biểu thức hệ số Ettingshausen $P$ và hệ số Peltier $\Pi$: $$P = \frac{1}{4\pi B}\frac{\sigma_{xx}\gamma_{xy} - \sigma_{xy}\gamma_{xx}}{\sigma_{xx}[\beta_{xx}\gamma_{xx} - \sigma_{xx}(\xi_{xx} - K_L)]}$$ $$\Pi = \gamma_{xx}\sigma_{xx}^{-1}$$
  2. Phương pháp tính số và công cụ mô phỏng:

    • Sử dụng phần mềm Matlab để lập trình tính toán số trị các biểu thức tích phân và giải tích phức tạp của ten-xơ động.
    • Dữ liệu tham số cho cấu trúc $\text{GaAs/AlGaAs}$: $m_e = 0{,}067 m_0$, $e = 2{,}07 e_0$, $\varepsilon_F = 50\text{ meV}$, $\tau(\varepsilon_F) = 10^{-12}\text{ s}$, $v_s = 5378\text{ m/s}$, $\rho = 5{,}32\text{ g/cm}^3$, $\omega_z = 2\cdot 10^{12}\text{ Hz}$, $L_y = 2\text{ nm}$.

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

Chương 1: Tổng quan về hiệu ứng Ettingshausen và hiệu ứng Peltier trong bán dẫn khối

  • Cơ chế vật lý: Khi đặt mẫu bán dẫn trong từ trường $\vec{B} = (0, 0, B)$ và điện trường $\vec{E} = (E, 0, 0)$ có mật độ dòng $j_x$, lực Lorentz làm lệch các điện tích theo trục $y$, tạo nên gradient nhiệt độ $\nabla_y T$. Hệ số Ettingshausen được xác định theo công thức $P = -\frac{\nabla T}{j B}$. Hiệu ứng Peltier biểu diễn mật độ thông lượng nhiệt sinh ra tỉ lệ với mật độ dòng điện: $\vec{Q} = \Pi \vec{j}$.
  • Thiết lập lý thuyết lượng tử: Chương này thiết lập hệ phương trình vi phân Heisenberg cho hàm tương quan $F_{\vec{k}_1, \vec{k}2, \vec{q}}(t)$, dẫn ra nghiệm giải tích chứa các hàm Bessel $J_s(\vec{a}\vec{q}) J\ell(\vec{a}\vec{q})$ với $\vec{a} = \frac{e\vec{E}_0}{m_e \omega^2}$.
  • Kết quả giải tích ten-xơ: Tìm được biểu thức các thành phần ten-xơ độ dẫn $\sigma_{xx}, \sigma_{xy}$, ten-xơ nhiệt điện $\beta_{xx}$, ten-xơ hệ số Peltier $\gamma_{xx}, \gamma_{xy}$ và ten-xơ dẫn nhiệt $\xi_{xx}$ trong bán dẫn khối. Đây là hệ quy chiếu lý thuyết chuẩn để đối sánh với các cấu trúc thấp chiều trong các chương sau.

Chương 2: Ảnh hưởng của sự lượng tử hóa do giảm kích thước lên hiệu ứng Ettingshausen và hiệu ứng Peltier trong hố lượng tử với thế giam giữ parabol

  • Mô hình thế giam giữ: Hố lượng tử parabol (PQW) có thế $V(z) = \frac{1}{2}m_e \omega_z^2 z^2$ với bề rộng hố $L$. Điện trường đặt dọc trục $x$, sóng điện từ mạnh dọc trục $y$ ($\vec{E}_0(t) = (0, E_0\sin\omega t, 0)$) và từ trường dọc trục $z$ ($\vec{B} = (0, 0, B)$).
  • Phổ năng lượng lượng tử hóa:
    • Điện tử bị giam cầm theo phương $z$: Phổ năng lượng nhận các giá trị gián đoạn $\varepsilon_{N,n}^{\text{PQW}}(k_y) = \left(N + \frac{1}{2}\right)\hbar\omega_c + \left(n + \frac{1}{2}\right)\hbar\omega_z + \frac{m_e v_d^2}{2} - \hbar v_d k_y$, trong đó hàm sóng phụ thuộc đa thức Hermite $H_n(z/\ell_z)$.
    • Tần số phonon âm giam cầm (CAP): $\omega_{m,\vec{q}\perp}^{\text{PQW}-A} = v_s \sqrt{q\perp^2 + (m\pi/L)^2}$.
    • Tần số phonon quang giam cầm (COP): $\omega_{m,\vec{q}\perp}^{\text{PQW}-O} = \sqrt{\omega_0^2 - \upsilon^2 [q\perp^2 + (m\pi/L)^2]}$.
  • Điều kiện cộng hưởng từ - phonon - photon (MPPRC): Phương trình động lượng tử cho thấy xuất hiện các đỉnh cộng hưởng khi thỏa mãn điều kiện: $$(N - N')\hbar\omega_c + (n - n')\hbar\omega_z + eE\bar{r} \pm \hbar\omega_m \pm \hbar\omega = 0$$ Thành phần $\hbar\omega_m$ ($= \hbar v_s m\pi/L$ với CAP hoặc $\hbar\omega_m^O$ với COP) đóng góp trực tiếp vào vị trí các đỉnh cộng hưởng, tạo nên sự dịch chuyển và tách vạch phổ so với trường hợp phonon không giam cầm (un-CAP, un-COP).
  • Phát hiện qua tính số trên $\text{GaAs/AlGaAs}$:
    • Hệ số EC dao động mạnh với biên độ lớn ở dải tần số sóng điện từ thấp $\omega < 10^{12}\text{ Hz}$. Khi $\omega > 10^{12}\text{ Hz}$, EC dao động điều hòa với biên độ nhỏ dần.
    • Trường hợp phonon giam cầm (CAP, COP) cho giá trị cực đại của EC và PC khác biệt rõ rệt so với trường hợp un-CAP và un-COP. Khi cho chỉ số lượng tử $m \to 0$, các biểu thức giải tích của luận án tiệm cận chính xác về kết quả của trường hợp phonon không giam cầm đã công bố trước đây.
Thông số khảo sát trong PQW Giá trị vật lý thiết lập Hành vi của hệ số EC ($P$) và PC ($\Pi$)
Tần số sóng điện từ $\omega$ $0 - 6\cdot 10^{12}\text{ Hz}$ ($E_0 = 2\cdot 10^4\text{ V/m}$) Xuất hiện các đỉnh cộng hưởng MPPRC; dao động mạnh ở $\omega < 10^{12}\text{ Hz}$
Nhiệt độ môi trường $T$ $0 - 20\text{ K}$ ($E = 1{,}3\cdot 10^3\text{ V/m}$) EC và PC thay đổi phi tuyến, CAP cho độ lớn vượt trội so với un-CAP
Từ trường ngoài $B$ $0 - 3\text{ T}$ ($T = 100\text{ K}$) Phổ xuất hiện cấu trúc cộng hưởng Landau chuyển mức giữa các chỉ số $N, N'$
Biên độ sóng điện từ $E_0$ $10^4 - 10^7\text{ V/m}$ ($T = 300\text{ K}$) Cường độ trường ngoài làm biến điệu độ cao của các đỉnh hấp thụ photon

Chương 3: Ảnh hưởng của sự lượng tử hóa do giảm kích thước lên hiệu ứng Ettingshausen và hiệu ứng Peltier trong siêu mạng pha tạp

  • Cấu trúc siêu mạng pha tạp (DSS): Khảo sát cấu trúc gồm các lớp bán dẫn cùng loại được pha tạp xen kẽ các lớp donor ($n$-type) và acceptor ($p$-type), tạo nên thế chu kỳ theo phương giam giữ.
  • Biểu thức giải tích: Thiết lập phương trình động lượng tử cho DSS có tính đến sự giam cầm đồng thời của điện tử và phonon. Biểu thức giải tích của ten-xơ độ dẫn $\sigma_{i\nu}$ và ten-xơ nhiệt điện $\beta_{i\nu}$ phụ thuộc vào chu kỳ siêu mạng và chỉ số lượng tử $m$.
  • Kết quả tính số:
    • Khảo sát sự biến thiên của các thành phần ten-xơ độ dẫn và ten-xơ nhiệt điện theo từ trường $B$ với $E_0 = 4{,}8\cdot 10^4\text{ V/m}$.
    • Hệ số PC và EC phụ thuộc cộng hưởng vào tần số cyclotron $\omega_c$.
    • Sự phụ thuộc vào nồng độ pha tạp được khảo sát tại điều kiện $T = 5\text{ K}$ và $B = 1{,}5\text{ T}$, chỉ ra ngưỡng nồng độ tối ưu để đạt hệ số chuyển đổi nhiệt điện cao.

Chương 4: Ảnh hưởng của sự lượng tử hóa do giảm kích thước lên hiệu ứng Ettingshausen và hiệu ứng Peltier trong siêu mạng hợp phần

  • Cấu trúc siêu mạng hợp phần (CSS): Khảo sát hệ cấu trúc tuần hoàn tạo bởi hai vật liệu bán dẫn khác nhau với bề rộng lớp $d$ (khảo sát cụ thể tại $d = 15\text{ nm}$).
  • Cơ chế tán xạ: Khảo sát đầy đủ cả hai kênh tán xạ điện tử giam cầm - CAP và điện tử giam cầm - COP trong điều kiện có sóng điện từ mạnh.
  • Kết quả tính số:
    • Khảo sát sự phụ thuộc của EC và PC vào từ trường với các mức biên độ $E_0 = 10^6\text{ V/m}$ và $E_0 = 2\cdot 10^6\text{ V/m}$ tại tần số $\omega = 3\cdot 10^{12}\text{ Hz}$.
    • Đồ thị phụ thuộc nhiệt độ $T$ trong dải nhiệt độ thấp ($T = 3\text{ K}$) và dải nhiệt độ cao cho thấy sự chuyển đổi vai trò chủ đạo từ tán xạ CAP sang tán xạ COP.

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

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

  1. Thiết lập khung lý thuyết hoàn chỉnh: Xây dựng thành công hệ phương trình động lượng tử cho hàm phân bố điện tử trong các hệ bán dẫn hai chiều (QW, DSS, CSS) dưới tác dụng đồng thời của điện trường không đổi, từ trường và sóng điện từ mạnh, có tính đến sự giam cầm của cả điện tử và phonon (CAP và COP).
  2. Xác lập biểu thức giải tích mới: Tìm ra các công thức giải tích tường minh cho các ten-xơ động (độ dẫn điện, hệ số nhiệt điện, hệ số dẫn nhiệt điện tử), hệ số Ettingshausen (EC) và hệ số Peltier (PC). Các biểu thức này chứa tường minh các số lượng tử giam cầm điện tử ($N, n$) và số lượng tử giam cầm phonon ($m$).
  3. Phát hiện điều kiện cộng hưởng mới: Xác lập điều kiện cộng hưởng từ - phonon - photon (MPPRC) trong hố lượng tử parabol và siêu mạng bán dẫn, chỉ ra sự đóng góp của năng lượng phonon giam cầm $\hbar\omega_m$ vào các dịch chuyển vạch phổ cộng hưởng.
  4. Tính nhất quán và bao hàm: Chứng minh giải tích rằng khi chỉ số giam cầm phonon $m \to 0$, các biểu thức giải tích của luận án suy biến chính xác về kết quả của bài toán phonon không giam cầm (un-confined) và tiệm cận về trường hợp bán dẫn khối khi kích thước hệ tăng vô hạn.

Đóng góp về mặt thực tiễn và công bố khoa học

  • Cơ sở dữ liệu vật liệu: Cung cấp tập dữ liệu tính số định lượng chi tiết trên vật liệu bán dẫn $\text{GaAs/AlGaAs}$, làm cơ sở định hướng tính toán thiết kế cho các linh kiện nhiệt điện vi mô và thiết bị làm lạnh trạng thái rắn.
  • Công bố khoa học: Các kết quả chính của luận án đã được công bố trong 07 công trình khoa học, bao gồm:
    • 02 bài báo trên các tạp chí quốc tế thuộc danh mục ISI.
    • 02 bài báo trên các tạp chí quốc tế thuộc danh mục Scopus.
    • 03 bài báo trên các tạp chí khoa học chuyên ngành trong nước.

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

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

  • Luận án hoàn toàn là công trình nghiên cứu lý thuyết giải tích và tính toán mô phỏng số, chưa có điều kiện kết hợp đo đạc thực nghiệm trực tiếp trên các mẫu màng mỏng nano chế tạo bằng MBE.
  • Các biểu thức giải tích thu được dựa trên phép gần đúng tuyến tính đối với cường độ trường ngoài và chỉ xét quá trình hấp thụ/phát xạ không quá một photon ($\ell = 0, \pm 1$).
  • Chỉ khảo sát hai cơ chế tán xạ chính là phonon âm (thế biến dạng) và phonon quang (phân cực), chưa xét đến tán xạ tạp chất ion hóa, tán xạ áp điện hay tán xạ giữa các phân vùng năng lượng.

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

  • Mở rộng lý thuyết phương trình động lượng tử cho các quá trình đa photon ($\ell > 1$) dưới tác động của sóng điện từ có cường độ siêu mạnh.
  • Nghiên cứu ảnh hưởng của sự giam cầm phonon lên EE và PE trong các cấu trúc thấp chiều khác như dây lượng tử (1D), chấm lượng tử (0D) và siêu mạng dị thể cấu trúc phức tạp.
  • Khảo sát bổ sung các cơ chế tán xạ khác (tán xạ tạp chất, tán xạ nhám bề mặt) trong điều kiện nhiệt độ cực thấp hoặc nồng độ pha tạp cao.
  • Phối hợp với các nhóm nghiên cứu thực nghiệm để tiến hành đo đạc kiểm chứng các đỉnh cộng hưởng MPPRC trên các cấu trúc dị thể $\text{GaAs/AlGaAs}$.

Giá trị tham khảo

Đối tượng thụ hưởng

  • Nghiên cứu sinh và học viên cao học: Chuyên ngành Vật lý lý thuyết, Vật lý toán, Vật lý chất rắn, Vật lý bán dẫn và Khoa học vật liệu nano.
  • Giảng viên và nhà nghiên cứu: Dùng làm tài liệu tham khảo chuyên sâu cho các học phần Vật lý bán dẫn thấp chiều, Lý thuyết lượng tử chất rắn và Các hiệu ứng động nhiệt điện.
  • Kỹ sư công nghệ bán dẫn: Tham khảo các tham số vật liệu và điều kiện cộng hưởng để định hướng tối ưu hóa cấu trúc vi linh kiện nhiệt điện và cảm biến từ - nhiệt.

Các nội dung có giá trị tham khảo cao

  • Quy trình toán học từng bước thiết lập và giải phương trình vi phân Heisenberg cho hàm tương quan điện tử - phonon bằng phương pháp biến thiên hằng số.
  • Kỹ thuật tính toán các ten-xơ động học lượng tử dưới tác dụng đồng thời của ba trường ngoài (điện trường, từ trường, sóng điện từ).
  • Thuật toán và mô hình tính số trên phần mềm Matlab mô phỏng hệ số Ettingshausen và hệ số Peltier trong hố lượng tử và siêu mạng.

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

1. Hiệu ứng Ettingshausen và hiệu ứng Peltier khác nhau như thế nào về bản chất vật lý trong luận án?

Hiệu ứng Ettingshausen (EE) là hiệu ứng từ - nhiệt - điện xuất hiện khi đặt vật liệu trong từ trường vuông góc với dòng điện, dẫn đến sự tích tụ hạt tải ở một phía do lực Lorentz và tạo ra gradient nhiệt độ ngang $\nabla T$ (đặc trưng bởi hệ số $P = -\nabla T / (j B)$). Hiệu ứng Peltier (PE) là hiện tượng sinh ra hoặc hấp thụ nhiệt tại mặt tiếp xúc giữa hai vật liệu khác nhau khi có dòng điện chạy qua, đặc trưng bởi hệ số $\Pi = \vec{Q} / \vec{j}$.

2. Sự giam cầm của phonon (confined phonon) làm thay đổi tần số dao động mạng như thế nào so với phonon không giam cầm?

Trong hố lượng tử parabol, do hiệu ứng kích thước, thành phần véc-tơ sóng theo phương giam giữ bị lượng tử hóa thành $q_m = \frac{m\pi}{L}$ (với $m$ là số nguyên). Khi đó, tần số của phonon âm giam cầm trở thành $\omega_{m,\vec{q}\perp}^A = v_s \sqrt{q\perp^2 + (m\pi/L)^2}$ và tần số phonon quang giam cầm trở thành $\omega_{m,\vec{q}\perp}^O = \sqrt{\omega_0^2 - \upsilon^2 [q\perp^2 + (m\pi/L)^2]}$, tạo nên sự phụ thuộc tường minh vào chỉ số $m$.

3. Điều kiện cộng hưởng từ - phonon - photon (MPPRC) trong hố lượng tử parabol có dạng toán học như thế nào?

Điều kiện cộng hưởng MPPRC được xác lập dưới dạng: $$(N - N')\hbar\omega_c + (n - n')\hbar\omega_z + eE\bar{r} \pm \hbar\omega_m \pm \hbar\omega = 0$$ Trong đó $N, N'$ là chỉ số mức Landau, $n, n'$ là chỉ số mức con điện tử, $\omega_z$ là tần số giam giữ thế parabol, $\omega_m$ là tần số phonon giam cầm (CAP hoặc COP) và $\omega$ là tần số sóng điện từ mạnh.

4. Tại sao luận án lại sử dụng phương pháp phương trình động lượng tử thay vì phương trình động Boltzmann cổ điển?

Phương trình động Boltzmann cổ điển chỉ áp dụng chính xác ở vùng nhiệt độ cao và không mô tả được đầy đủ các hiệu ứng lượng tử hóa năng lượng gián đoạn cũng như các chuyển mức lượng tử khi có trường ngoài mạnh. Phương pháp phương trình động lượng tử (dựa trên Hamiltonian lượng tử hóa lần thứ hai và phương trình chuyển động Heisenberg) cho phép thu được nghiệm giải tích đúng trên toàn bộ dải nhiệt độ, bao gồm cả vùng nhiệt độ thấp và điều kiện trường ngoài biến thiên nhanh.

5. Luận án đã sử dụng bộ thông số vật liệu nào để tính toán số cho cấu trúc $\text{GaAs/AlGaAs}$?

Bộ thông số được sử dụng gồm: khối lượng hiệu dụng điện tử $m_e = 0{,}067 m_0$, điện tích hiệu dụng $e = 2{,}07 e_0$, năng lượng Fermi $\varepsilon_F = 50\text{ meV}$, thời gian phục hồi năng lượng $\tau(\varepsilon_F) = 10^{-12}\text{ s}$, vận tốc sóng âm $v_s = 5378\text{ m/s}$, mật độ khối lượng $\rho = 5{,}32\text{ g/cm}^3$, tần số giam giữ $\omega_z = 2\cdot 10^{12}\text{ Hz}$ và bề rộng hố $L_y = 2\text{ nm}$.


Kết luận

Luận án tiến sĩ của tác giả Nguyễn Thị Lâm Quỳnh đã xây dựng hoàn chỉnh lý thuyết lượng tử về hiệu ứng Ettingshausen và hiệu ứng Peltier trong các hệ bán dẫn hai chiều gồm hố lượng tử parabol, siêu mạng pha tạp và siêu mạng hợp phần dưới tác động đồng thời của điện trường không đổi, từ trường và sóng điện từ mạnh. Bằng việc áp dụng phương pháp phương trình động lượng tử, công trình đã xác lập hệ biểu thức giải tích của các ten-xơ động và các hệ số nhiệt điện có tính đến sự giam cầm đồng thời của cả điện tử và phonon (CAP và COP). Các kết quả tính số trên cấu trúc $\text{GaAs/AlGaAs}$ đã làm sáng tỏ điều kiện cộng hưởng từ - phonon - photon (MPPRC) và chứng minh vai trò rõ rệt của phonon giam cầm so với trường hợp phonon không giam cầm và bán dẫn khối. Toàn bộ kết quả nghiên cứu đã được công bố trên 07 bài báo khoa học uy tín, cung cấp cơ sở lý thuyết có giá trị cho lĩnh vực vật lý bán dẫn thấp chiều và công nghệ vật liệu nhiệt điện nano.