Трьохвимірна модель переносу заряду: легований резистор на основі арсеніду галію (GaAs).

У цьому навчальному посібнику показано, як перейти від простої 2D-моделі дрейфу-дифузії до повноцінної 3D-електричної симуляції для простого, але дуже наочного пристрою: легованого галій-арсеніду (GaAs) резистора. Він ідеально підходить для навчання та використання в практичних завданнях: фізичні процеси чітко відображені, залежність струму від напруги є майже лінійною, а результати наочно демонструють, як розподіляються потенціали та струми в просторі.

1. Створіть двовимірний резистор на основі арсеніду галію (GaAs).

1. Нова симуляція → Демонстрація GaAs → Легована дротова структура / резистор.
2. Геометрія: почніть із двовимірного простору, використовуючи точки сітки вздовж осей X і Y. (Наразі координату Z залиште вимкненою; увімкнемо її пізніше.)
3. Допування: визначте профіль акцепторів, який збільшується зліва направо (p-тип градієнта). «Дріт» – це просто прямокутний блок GaAs, до якого застосовано цей градієнт.
4. Контакти: протилежні поверхні діють як електроди (верхня/нижня або ліва/права, залежно від обраної вами системи координат). До однієї з них подаватиметься змінна напруга, а інша підтримуватиметься на рівні 0 В.
5. Рекомбінація: увімкніть просту модель SRH (аналітична форма). Для цієї демонстрації не потрібні динамічні пастки.

2. Запуск і перевірка на наявність помилок (у двовимірному просторі).

Натисніть «Виконати». У журналі:

• Тривалість кожного кроку має становити кілька мілісекунд (на звичайному ноутбуці – близько 6–8 мс).
• Контактні струми мають збігатися (у стаціонарному стані верхній і нижній струми мають бути однаковими). Якщо вони не збігаються, ймовірно, ви працюєте в часовій області або процес не зійшовся.
• Залишкові значення (загальна похибка) мають бути невеликими: значення, близьке до 10⁻⁹, є відмінним. Значення, такі як 10⁻⁶…10⁻³, свідчать про необхідність повторної перевірки налаштувань/сітки.

Відкрити вихідний файл:

• Вольт-амперна характеристика: спостерігається майже лінійна залежність, що відповідає очікуваній поведінці резистора. Незначні коливання можуть виникати через наявність просторового заряду та домішок.
• Графіки: відобразіть залежності енергії провідності (іноді позначається як «LUMO» в загальному графічному інтерфейсі), електростатичного потенціалу (ϕ) та концентрації носіїв заряду (n, p) від зміни прикладеної напруги.

3. Розширте можливості до тривимірного простору: сітка та контакти.

Увімкніть розмірність Z у редакторі сітки. Спочатку оберіть помірні розміри: Nx × Ny × Nz = 5 × 5 × 5. Обсяг обчислень і використання пам’яті приблизно пропорційні \( \mathcal{O}(N_x N_y N_z) \approx \mathcal{O}(N^3) \), тому 10×10×10 підійде для ноутбука, але 20×20×20 або більше можуть швидко призвести до значного навантаження.

У тривимірному просторі контакти представлені об’єктами з обмеженою площею (на відміну від двовимірного простору, де вони неявно охоплюють усю поверхню). Відкрийте редактор «Розміри/Контакти» та:

• Встановіть відповідну ширину/глибину для кожного контакту, щоб вона відповідала розміру пристрою (для електродів, що покривають всю поверхню обличчя), або оберіть часткові розміри (наприклад, 0,5 × 0,5 від розміру обличчя), щоб вивчити розподіл струму.
• Відрегулюйте зміщення, щоб розмістити провідні майданчики там, де вам потрібно (це демонструє асиметричне введення/виведення).

4. Запуск і аналіз (у 3D).

Запустіть 3D-моделювання (переконайтеся, що сітка не надто складна). Перевірте:

• Графік залежності напруги від струму (V–I характеристика): має залишатися близьким до лінійної форми (резистивної). Порівняйте зі значеннями для двовимірного випадку.
• Можливий вигляд ϕ у тривимірному просторі: спостерігайте за зміною градієнта та еволюцією тривимірного розподілу потенціалу. Використовуйте інструменти для створення зрізів, щоб «розрізати» пристрій і отримати чіткішу інтерпретацію.
• Потоки: визначте шість складових:
  • Jn,x, Jn,y, Jn,z
  • Jp,x, Jp,y, Jp,z
  Завдяки зміщеним контактам ви зможете побачити зони входу/виходу та розподіл струму. Повторюйте переміщення повзунка, щоб спостерігати за змінами напруги.
• Carriers: n(x,y,z,V), p

5. Практичні поради та можливі труднощі

• Спочатку створіть прототип у 1D/2D. Перевірте фізичні та чисельні моделі, а потім активуйте вісь Z лише за потреби.
• Використання більш дрібної сітки дає кращі результати, ніж застосування грубої сили. Збільшення кількості точок не обов’язково підвищує точність, якщо рівняння в часткових похідних розв’язуються з урахуванням невеликих похибок між вузлами. Збільшуйте розмір сітки лише для точного відображення реальних градієнтів.
• Перевірка збіжності: забезпечте малі значення похибок (≲10−8), відповідність контактних струмів і стабільність кроку обчислення. Якщо значення похибок стабілізуються на рівні 10−5…10−3, перегляньте параметри масштабування, зміщення або сітку.
• Вибір SRH: для демонстрації роботи цього резистора достатньо використовувати просту аналітичну модель SRH. Увімкніть динамічні пастки лише тоді, коли вам потрібні дані про кінетику пасток.
• Збільшення обсягу пам’яті: подвоєння сітки вздовж кожної осі призводить до збільшення кількості невідомих параметрів приблизно в 8 разів; слідкуйте за обсягом оперативної пам’яті для конфігурацій 20×20×20 і вище.