Навчальний посібник з GaAs (Частина C): Числова стабільність і перетворення діода на фотодетектор.
1. Вступ
У цьому розділі ви дізнаєтеся, чому числова стабільність стає ключовим питанням у дво- та тривимірних моделюваннях дрейфу-дифузії, і чому артефакти часто виникають саме в режимі низьких струмів. Це не є помилкою в програмному забезпеченні, а скоріше фундаментальним наслідком розв’язання тісно пов’язаних рівнянь переносу заряду в напівпровідниках у надзвичайно широкому динамічному діапазоні. Розуміння цієї поведінки є важливим під час інтерпретації кривих залежності струму від напруги, карт щільності струму та знімків реалістичних багатовимірних моделей пристроїв.
2. Чому числова стабільність має більше значення у двовимірних і тривимірних моделях.
Моделювання дрейфу та дифузії напівпровідників зазвичай охоплює дуже широкий динамічний діапазон. Концентрації носіїв заряду можуть відрізнятися на 10–20 і більше порядків величини між основними та другорядними областями (наприклад, сильно легована n-типу речовина порівняно з другорядними дірками). Коли ці масштаби потрапляють у пов’язану нелінійну систему (рівняння Пуассона плюс рівняння неперервності для електронів і дірок), отриманий якобіан може стати погано обумовленим. На практиці це означає, що невеликі числові похибки округлення можуть посилюватися і перетворюватися на помітні артефакти, особливо коли справжній фізичний струм близький до нуля.
Більшість чисельних методів розв’язання рівнянь переносу-дифузії використовують обчислення з подвійною точністю. Подвійне число містить приблизно 53 біти мантиси, що відповідає близько 15–16 десятковим знакам точності. Хоча цього достатньо для переважної більшості задач, це не є нескінченною точністю. Коли рівняння містять величини, які відрізняються на багато порядків, основною чисельною проблемою стають взаємне скорочення, масштабування та приведення матриць до належного вигляду, а не просто швидкість процесора. В одномірному випадку ці проблеми часто менш виражені, оскільки існує менше ступенів свободи та менше шляхів взаємодії. У двомірному та тримірному випадках додаткові просторові напрямки призводять до більших матриць і сильнішої взаємодії, тому числова чутливість, як правило, проявляється швидше, особливо при низьких значеннях струму. У решті цього розділу ви дізнаєтеся, як ці числові ефекти виникають на практиці, як вони проявляються в багатовимірних моделюваннях і як правильно їх інтерпретувати під час аналізу двомірних і тримірних результатів.
3. Забезпечення світловипромінювання: перетворення діода на фотодіод.
До цього моменту пристрій на основі арсеніду галію моделювався як діод, який не пропускає світло. У цьому розділі ми активуємо освітлення та спостерігаємо, як пристрій реагує на оптичне генерування носіїв заряду. Це дає нам змогу пов’язати обговорення чисельної стабільності з попереднього розділу з фізично значущою кривою залежності струму та напруги при освітленні та визначити, де на практиці найчастіше виникають артефакти при низьких значеннях струму.
Щоб увімкнути підсвічування, перейдіть на вкладку «Оптика» (??) і встановіть значення параметра «Інтенсивність світла (сонячні одиниці)» на 1,0. Після ввімкнення в основному вікні пристрою над стовпчиком з’являться зелені стрілки, що вказуватимуть на падаючі фотони (??).
🧭 Для цього розділу ми свідомо використовуємо двовимірну симуляцію: відкрийте редактор електричної сітки та переконайтеся, що параметри x і y увімкнено, а параметр z – вимкнено. Це дозволяє скоротити час обчислень, водночас демонструючи чисельну поведінку, яку ми хочемо дослідити.
Запустіть симуляцію та відкрийте jv_contact0.csv (верхній контакт). Отримана крива має нагадувати фізично правдоподібну характеристику фотоелемента, але якщо ви збільшите масштаб у зоні низьких струмів, ви часто помітите невеликі нерівності або «згустки» (??).
На цьому етапі два ефекти накладаються один на одного та підсилюють один одного. По-перше, коли реальний фізичний струм близький до нуля, числова чутливість, про яку йшлося раніше, стає помітною на отриманій кривій залежності струму від напруги. По-друге, абсолютний фотострум є відносно невеликим, оскільки структура все ще сконфігурована як діод, а не як оптимізована сонячна батарея. Цей невеликий фотострум зміщує робочу точку глибше в область близьких до нуля значень струму, що ще більше посилює видимість числових артефактів.
jv_contact0.csv). Нерівномірність при низьких значеннях струму відображає числову чутливість, коли реальний струм близький до нуля.4. Діагностика числових шумів за допомогою серії знімків просторової густини струму
Щоб зрозуміти походження невеликих відхилень, які спостерігаються на кривій залежності фотоструму від напруги, корисно проаналізувати внутрішні процеси, які виконує обчислювальна програма. Це можна зробити за допомогою вихідних даних snapshots/, які зберігають значення величин, розв’язаних у просторі, як функцію від кроку моделювання (зазвичай це прикладена напруга). У цьому розділі ми зосереджуємося на Jn.csv, тобто на вертикальній густині струму електронів.
Відкрийте папку snapshots/ і створіть графік Jn.csv. У результаті ви повинні отримати графіки, подібні до ?? і ??, які відповідають низькому та вищому значенням прикладеної напруги відповідно.
За низької полярності (при зміщенні повзунка вліво), абсолютна величина струму надзвичайно мала. У цьому режимі графік щільності струму може виглядати просторово нерівномірним: ви фактично бачите, наскільки складно обчислити величину, яка фізично близька до нуля. У міру збільшення полярності та значного відхилення струму від нуля, просторовий розподіл Jn.csv стає більш рівномірним і зрозумілим з фізичної точки зору. Ці невеликі просторові нерівномірності при дуже низькому струмі безпосередньо впливають на отриману криву залежності струму від напруги (JV), тому область низьких значень струму на кривій JV може виглядати дещо «нерівномірно».
Один із практичних способів зменшити видимість цих артефактів – змінити місце, де обчислюється струм. Струми, які обчислюються дуже близько до контактів, як правило, більш чутливі до числових шумів, ніж струми, які обчислюються в середині пристрою. У цьому посібнику струм уже обчислюється в середній точці пристрою, і це налаштування можна знайти в розділі «Електричні параметри» в меню конфігурації «Дрейфова дифузія» → «Розв’язувач». Обчислення струму в середній точці зменшує чутливість до граничних ефектів і зазвичай забезпечує більш чітку криву залежності струму від напруги (JV) при низьких значеннях струму.
Jn.csv при низькому рівні зміщення у двовимірному просторі. Фактичний струм дуже малий, тому стає помітною числова чутливість.
Jn.csv при вищій напрузі зміщення. Як тільки струм стає помітно відмінним від нуля, просторовий профіль стає більш плавним і реалістичним.
5. Усунення проблеми низького фотоструму: не направляйте промінь через товстий металевий шар.
Ви, можливо, помітили, що отриманий досі фотострум є нижчим, ніж зазвичай можна було б очікувати для пристрою на основі арсеніду галію (GaAs). У цьому прикладі така поведінка є переважно запланованою. Спочатку структура була розроблена як діод, а не як оптимізований фотодіод або сонячна батарея: верхній алюмінієвий контакт достатньо товстий, щоб сильно відбивати та поглинати падаюче світло. Іншими словами, моделювання фактично намагається «пропустити світло крізь метал».
Щоб виправити це, відкрийте редактор шарів (??) і зменште товщину верхнього алюмінієвого (Al) контакту з 100 нм до 10 нм. Повторно запустіть моделювання та знову відкрийте графік залежності струму від напруги при освітленні.
Тепер ви повинні помітити значно більший фотострум і набагато більш плавний діапазон низьких струмів. Це покращення досягається не шляхом зміни чисельного методу, а шляхом збільшення реального фізичного сигналу. Віддаляючи пристрій від режиму, близького до нульового струму, чисельна межа фактично виводиться за межі видимості, і раніше помітні артефакти зникають.
6. Підсумок та основні висновки
У цьому навчальному посібнику, що складається з трьох частин, ви розглянули, як розмірність, числова стабільність і фізична конфігурація взаємодіють у моделюванні процесів дрейфу та дифузії:
- У частині А ви провели тестування повністю тривимірного діода на основі арсеніду галію (GaAs), який містив штучно створений вертикальний дефект, візуалізували зони концентрації струму та проаналізували внутрішні параметри за допомогою
snapshots/. - У розділі B ви усунули дефект і показали, що після того, як структура стає однорідною в бічній площині, тривимірне, двовимірне та одновимірне моделювання дає однаковий електричний результат, хоча час обчислення змінюється в рази.
- У розділі C ви використали освітлення для створення фотопристрою, застосували просторові графіки щільності струму для діагностики числової чутливості та продемонстрували, як низький фотострум підсилює видимі артефакти.
Основна ідея: чисельні артефакти найбільш помітні, коли реальний фізичний струм близький до нуля. Режими низького струму — незалежно від того, чи викликані вони зміщенням, рекомбінацією або оптичним блокуванням — за своєю суттю є найскладнішими випадками для будь-якого чисельного розв’язувача рівнянь переносу-дифузії. Як тільки струм стає значно відмінним від нуля, той самий чисельний метод зазвичай демонструє набагато кращу збіжність.
Куди рухатися далі: дослідіть, як змінюється числова стабільність зі збільшенням динамічного діапазону. Корисними напрямками є зміна параметрів рекомбінації за механізмом SRH або Оже, варіація контрасту легування або рухливості, збільшення інтенсивності освітлення або введення більш реалістичних параметрів селективності контакту. Кожен із цих підходів впливає на розв’язувач по-різному та допомагає сформувати інтуїтивне розуміння того, коли потрібне багатовимірне моделювання, яке водночас є чисельно надійним.