Trójwymiarowa symulacja dryfu i dyfuzji: rezystor z domieszkowanego arsenku galu (GaAs).

W tym samouczku pokazano, jak przejść od prostego, dwuwymiarowego modelu dryfu i dyfuzji do pełnej, trójwymiarowej symulacji elektrycznej dla prostego, ale bardzo pouczającego urządzenia: rezystora z domieszkowanego arsenku galu (GaAs). Jest to idealne rozwiązanie do celów dydaktycznych i jako przykład: fizyka jest zrozumiała, charakterystyka prądowo-napięciowa jest prawie liniowa, a wyniki wyraźnie ilustrują, jak potencjały i prądy rozkładają się w przestrzeni.

1. Utwórz dwuwymiarowy rezystor z arsenku galu (GaAs).

1. Nowa symulacja → Prezentacja GaAs → Drut/rezystor domieszkowany.
2. Geometria: zacznij od pracy w przestrzeni dwuwymiarowej, definiując punkty siatki wzdłuż osi X i Y. (Na razie pozostaw oś Z wyłączoną; włączymy ją później).
3. Domieszkowanie: należy zdefiniować profil akceptora, który wzrasta od lewej do prawej strony (gradient typu p). „Przewód” to po prostu prostokątny blok GaAs, w którym zastosowano ten gradient.
4. Elektrody: przeciwległe powierzchnie pełnią funkcję elektrod (górna/dolna lub lewa/prawa, w zależności od przyjętej konwencji osi). Na jednej z nich przykładane jest zmienne napięcie, a druga utrzymywana jest na poziomie 0 V.
5. Rekombinacja: umożliwia prostą analizę SRH (w formie analitycznej). Do celów demonstracyjnych nie są wymagane dynamiczne pułapki.

2. Uruchom i sprawdź poprawność działania (w 2D).

Kliknij „Uruchom”. W dzienniku:

• Czas trwania pojedynczego kroku powinien wynosić kilka milisekund (na przeciętnym laptopie, około 6–8 ms).
• Wartości prądów styków są identyczne (w stanie ustalonym wartość górna równa się wartości dolnej). Jeśli tak nie jest, prawdopodobnie pracujesz w domenie czasowej lub algorytm nie zbiegł się.
• Wartości resztkowe (całkowity błąd) powinny być niewielkie: wartość rzędu 10−9 lub mniejsza jest bardzo dobra. Wartości takie jak 10−6…10−3 wskazują na konieczność ponownej weryfikacji ustawień/siatki.

Otwórz wynik.

• Krzywa charakterystyki prądowo-napięciowej (JV): wykazuje prawie liniowe zachowanie, co jest typowe dla rezystora. Niewielkie odchylenia mogą wystąpić z powodu obecności ładunków przestrzennych i domieszek.
• Przykładowe wykresy: przedstawiają zależność pasma przewodnictwa (czasami oznaczone jako „LUMO” w ogólnych schematach), potencjału elektrostatycznego (ϕ) oraz gęstości nośników (n, p) od zmieniającego się napięcia polaryzacyjnego.

3. Rozszerzenie do trzech wymiarów: siatka i punkty styku.

Włącz wymiar Z w edytorze siatki. Na początek wybierz rozsądne wartości: Nx × Ny × Nz = 5 × 5 × 5. Obciążenie procesora i zużycie pamięci rosną mniej więcej w tempie \( \mathcal{O}(N_x N_y N_z) \approx \mathcal{O}(N^3) \), więc 10×10×10 będzie odpowiednie dla laptopa, ale 20×20×20 lub więcej może szybko spowodować znaczne obciążenie.

W modelu 3D kontakty są obiektami o ograniczonej powierzchni (w przeciwieństwie do modelu 2D, gdzie domyślnie obejmują całą powierzchnię). Otwórz edytor wymiarów/kontaktów i:

• Ustaw szerokość/głębokość każdego z elektrod tak, aby odpowiadała rozmiarowi urządzenia (w przypadku elektrod pokrywających całą powierzchnię) lub wybierz częściowe wymiary (np. 0,5 × 0,5 powierzchni twarzy), aby zbadać rozkład prądu.
• Dostosuj przesunięcia, aby umieścić elektrody w żądanym miejscu (demonstruje asymetryczne wprowadzanie/pobieranie).

4. Uruchom i przeanalizuj (w 3D)

Uruchom symulację 3D (zadbaj o to, aby siatka była stosunkowo prosta). Sprawdź:

• Krzywa charakterystyki prądowo-napięciowej (JV): powinna zachowywać się w sposób zbliżony do liniowego (rezystywnego). Porównać z wartościami uzyskanymi w symulacji dwuwymiarowej.
• Potencjal ϕ w przestrzeni trójwymiarowej: obserwuj, jak zmienia się składowa stała i śledź ewolucję rozkładu potencjału w przestrzeni 3D. Użyj narzędzi do tworzenia przekrojów, aby „przeciąć” urządzenie i uzyskać lepszą interpretację.
• Prądy: zaznacz na wykresie sześć elementów:
  • Jn,x, Jn,y, Jn,z
  • Jp,x, Jp,y, Jp,z
  Dzięki zastosowaniu elektrod pomiarowych z przesuniętym kontaktem można zaobserwować obszary wejścia/wyjścia oraz rozkład prądu. Powtarzaj sekwencję zdjęć za pomocą suwaka, aby monitorować zmiany napięcia.
• Carriers: n(x,y,z,V), p

5. Praktyczne wskazówki i potencjalne problemy

• Najpierw należy stworzyć prototyp w 1D/2D. Następnie należy sprawdzić poprawność obliczeń i symulacji fizycznych, a dopiero potem włączyć komponent Z, gdy będzie to konieczne.
• Siatka obliczeniowa daje lepsze rezultaty niż metody oparte na sile. Większa liczba punktów niekoniecznie oznacza większą dokładność, jeśli równania różniczkowe cząstkowe są rozwiązywane z zachowaniem ściśle określonych wartości resztkowych między węzłami. Zwiększaj gęstość siatki tylko w celu dokładnego odwzorowania rzeczywistych gradientów.
• Sprawdzanie zbieżności: należy zapewnić małe wartości reszt (≲10−8), dopasowanie prądów w punktach kontaktu oraz stabilne kroki czasowe. Jeśli wartości reszt stabilizują się na poziomie 10−5…10−3, należy ponownie przeanalizować skalowanie, przesunięcia lub siatkę.
• Wybór metody SRH: prosta, analityczna metoda SRH jest wystarczająca do demonstracji działania tego rezystora. Włącz dynamiczne pułapki tylko wtedy, gdy potrzebujesz analizy kinetyki pułapek.
• Skalowanie pamięci: podwojenie liczby elementów siatki wzdłuż każdej osi powoduje około ośmiokrotny wzrost liczby niewiadomych; należy monitorować zużycie pamięci RAM, szczególnie w przypadku konfiguracji 20×20×20 i wyższych.