Τρισδιάστατη προσομοίωση διάχυσης-μεταφοράς φορτίου: αντιστάτης GaAs με προσθήκη προσμίξεων.

Αυτό το εκπαιδευτικό υλικό δείχνει πώς να μεταβείτε από μια απλή 2D προσομοίωση διάχυσης-μεταφοράς σε μια πλήρη 3D ηλεκτρική προσομοίωση για μια απλή –αλλά πολύ κατατοπιστική– συσκευή: έναν αντιστάτη γαλλίου-αρσενικού (GaAs) με προσθήκη προσμίξεων. Είναι ιδανικό για διδασκαλία και παραδείγματα: η φυσική είναι σαφής, η καμπύλη τάσης-ρεύματος είναι σχεδόν γραμμική και τα αποτελέσματα απεικονίζουν με ωραίο τρόπο πώς κατανέμονται οι δυναμικές και τα ρεύματα στον χώρο.

1. Δημιουργήστε την δισδιάστατη αντίσταση GaAs.

1. Νέα προσομοίωση → Επίδειξη με GaAs → Ενσωματωμένο σύρμα/αντιστάτη.
2. Γεωμετρία: ξεκινήστε σε 2D με σημεία πλέγματος κατά μήκος των αξόνων X και Y. (Αφήστε τον άξονα Z απενεργοποιημένο προς το παρόν· θα τον ενεργοποιήσουμε αργότερα.)
3. Ντόπινγκ: ορίστε ένα προφίλ αποδεκτών που αυξάνεται από αριστερά προς τα δεξιά (κλίση τύπου p). Το «σύρμα» είναι απλώς ένα ορθογώνιο μπλοκ GaAs με αυτή την κλίση.
4. Οι επαφές: οι αντίθετες επιφάνειες λειτουργούν ως ηλεκτρόδια (πάνω/κάτω ή αριστερά/δεξιά, ανάλογα με τη σύμβαση του άξονα που χρησιμοποιείτε). Η μία θα υποστεί μεταβολή τάσης, ενώ η άλλη θα διατηρηθεί σε 0 V.
5. Ανασυνδυασμός: ενεργοποιήστε την απλή SRH (αναλυτική μορφή). Δεν απαιτούνται δυναμικές παγίδες για αυτή την επίδειξη.

2. Εκτέλεση και έλεγχος λειτουργικότητας (2Δ)

Κάντε κλικ στο "Εκτέλεση". Στο αρχείο καταγραφής:

• Οι χρόνοι εκτέλεσης κάθε βήματος θα πρέπει να είναι λίγα χιλιοστά του δευτερολέπτου (σε έναν μέτριας απόδοσης φορητό υπολογιστή, της τάξης των 6–8 χιλιοστών του δευτερολέπτου).
• Τα ρεύματα επαφής πρέπει να ταιριάζουν (η τιμή στην κορυφή πρέπει να είναι ίση με την τιμή στο κάτω μέρος σε σταθερή κατάσταση). Εάν δεν ταιριάζουν, πιθανότατα βρίσκεστε σε ανάλυση στον τομέα του χρόνου ή η προσομοίωση δεν έχει συγκλίνει.
• Τα υπολείμματα (συνολικό σφάλμα) θα πρέπει να είναι μικρά: ≲10−9 είναι εξαιρετικό. Τιμές όπως 10−6…10−3 υποδεικνύουν ότι θα πρέπει να ελεγχθούν ξανά οι ρυθμίσεις/το πλέγμα.

Άνοιγμα εξόδου:

• Καμπύλη τάσης-ρεύματος (JV): παρουσιάζει σχεδόν γραμμική συμπεριφορά, όπως αναμένεται για έναν αντιστάτη. Μικρές διακυμάνσεις μπορεί να εμφανιστούν λόγω του φαινομένου της χωρικής φόρτισης και της προσθήκης προσμίξεων.
• Διαγράμματα: απεικονίστε την EC (ζώνη αγωγιμότητας· μερικές φορές επισημαίνεται ως «LUMO» σε γενικές διεπαφές χρήστη), το ϕ (ηλεκτροστατικό δυναμικό) και τις πυκνότητες φορέων (n, p) καθώς αυξάνεται η τάση πόλωσης.

3. Επέκταση σε 3D: πλέγμα και επαφές.

Ενεργοποιήστε τη διάσταση 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)

Εκτελέστε την τρισδιάστατη προσομοίωση (διατηρήστε το πλέγμα σε λογικό μέγεθος). Ελέγξτε:

• Καμπύλη JV: θα πρέπει να παραμένει περίπου γραμμική (αντιστάθμιση). Συγκρίνετε με τις τιμές 2D.
• Δυνατότητα απεικόνισης του ϕ σε 3D: παρατηρήστε την εξέλιξη της κλίσης και της τρισδιάστατης κατανομής του δυναμικού. Χρησιμοποιήστε εργαλεία τομής/διατομής για να «κόψετε» τη συσκευή και να διευκολύνετε την ερμηνεία.
• Ρεύματα: σχεδιάστε τα έξι συστατικά στοιχεία:
  • Jn,x, Jn,y, Jn,z
  • Jp,x, Jp,y, Jp,z
  Με τις μετατοπισμένες επαφές, θα δείτε τις περιοχές εισόδου/εξόδου και τη διάχυση του ρεύματος. Χρησιμοποιήστε την επαναλαμβανόμενη λειτουργία του ρυθμιστικού για να παρατηρήσετε τις τάσεις της τάσης.
• Carriers: n(x,y,z,V), p

5. Πρακτικές συμβουλές και πιθανές παγίδες.

• Δημιουργήστε πρώτα ένα πρωτότυπο σε 1Δ/2Δ. Ελέγξτε την ακρίβεια της φυσικής προσομοίωσης και των αριθμητικών υπολογισμών και, στη συνέχεια, ενεργοποιήστε τη διάσταση Ζ μόνο όταν είναι απαραίτητο.
• Η χρήση ενός πιο πυκνού πλέγματος είναι πιο αποτελεσματική από την απλή αύξηση της υπολογιστικής ισχύος. Η αύξηση του αριθμού των σημείων δεν συνεπάγεται απαραίτητα μεγαλύτερη ακρίβεια, εάν οι μερικές διαφορικές εξισώσεις επιλύονται με μικρά σφάλματα μεταξύ των κόμβων. Αυξήστε το πλέγμα μόνο για να αποτυπώσετε τις πραγματικές κλίσεις.
• Έλεγχοι σύγκλισης: επιβολή μικρών υπολοίπων (≲10−8), αντιστοίχιση ρευμάτων επαφής και σταθεροποίηση των χρονικών βημάτων. Εάν τα υπόλοιπα σταθεροποιηθούν στο 10−5…10−3, επανεξετάστε την κλίμακα, τις αποκλίσεις ή το πλέγμα.
• Επιλογή SRH: η απλή αναλυτική μέθοδος SRH είναι κατάλληλη για αυτή την επίδειξη με αντιστάσεις. Ενεργοποιήστε τις δυναμικές παγίδες μόνο όταν χρειάζεστε δεδομένα για την κινητική των παγίδων.
• Κλιμάκωση μνήμης: διπλασιασμός του πλέγματος κατά μήκος κάθε άξονα, περίπου 8 φορές περισσότερες άγνωστες μεταβλητές· παρακολουθήστε τη διαθέσιμη μνήμη RAM για επίπεδα 20×20×20 και άνω.