Οδηγός για το GaAs (Μέρος Β): Από 3D σε 2D σε 1D (Αφαίρεση του ελαττώματος)

1. Εισαγωγή: επιλογή μεταξύ των μεθόδων διάχυσης-μεταφοράς σε 3D, 2D και 1D.

Μια σημαντική απόφαση κατά τη δημιουργία μοντέλων στις προσομοιώσεις ημιαγωγικών συσκευών είναι η επιλογή της κατάλληλης διάστασης: 3D, 2D ή 1D. Τα μοντέλα με μεγαλύτερη διάσταση απαιτούν περισσότερους υπολογιστικούς πόρους, αλλά τα μοντέλα με μικρότερη διάσταση ενδέχεται να παραλείψουν σημαντικές φυσικές ιδιότητες εάν υπάρχουν πλευρικές επιδράσεις. Αυτή η ενότητα επικεντρώνεται στην κατανόηση του πότε απαιτείται πραγματικά μια μεγαλύτερη διάσταση. Στο Μέρος Α, η δίοδος GaAs περιείχε μια κάθετη ατέλεια που ανάγκαζε την πλευρική ροή ρεύματος, καθιστώντας το πρόβλημα εγγενώς τρισδιάστατο. Εδώ, αφαιρούμε αυτή την ασυμμετρία απενεργοποιώντας την ατέλεια, μετατρέποντας έτσι την ίδια συσκευή σε μια ομοιόμορφη δομή που μπορεί να επιλυθεί σε 3D, 2D και 1D χρησιμοποιώντας τις ίδιες φυσικές παραμέτρους.

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

2. Αφαιρέστε την πλευρική ασυμμετρία (απενεργοποιήστε τα ελαττωματικά αντικείμενα).

Στο Μέρος Α, η κάθετη ατέλεια προκάλεσε πλευρική ροή ρεύματος και κατέστησε τη συσκευή εγγενώς τρισδιάστατη. Για να μελετήσουμε πώς η διάσταση από μόνη της επηρεάζει τον χρόνο λειτουργίας και τα αποτελέσματα, αφαιρούμε τώρα αυτή την πηγή πλευρικής διακύμανσης, διατηρώντας παράλληλα ανέπαφη την υποκείμενη γεωμετρία. Αντί να διαγράψουμε την ατέλεια, θα την απενεργοποιήσουμε, έτσι ώστε η συσκευή να γίνει πλευρικά ομοιογενής στις κατευθύνσεις x και z.

Το ελάττωμα αποτελείται από δύο ξεχωριστά αντικείμενα. Στην τρισδιάστατη προβολή της συσκευής, κάντε δεξί κλικ σε ένα από τα μπλοκ του ελαττώματος και επιλέξτε "Επεξεργασία αντικειμένου" (όχι "Διαγραφή"), όπως φαίνεται στην εικόνα ??. Αυτό ανοίγει τον επεξεργαστή αντικειμένων (??), όπου μπορείτε να ενεργοποιήσετε ή να απενεργοποιήσετε την πράσινη επιλογή "Αντικείμενο ενεργοποιημένο". Μόλις απενεργοποιηθεί, το αντικείμενο επισημαίνεται με ένα κόκκινο "Χ" και τα πεδία του γίνονται γκρι (??). Επαναλάβετε αυτή τη διαδικασία για το δεύτερο μπλοκ του ελαττώματος.

Αφού απενεργοποιηθούν και τα δύο μπλοκ ελαττωμάτων, η κύρια προβολή της προσομοίωσης θα πρέπει να μοιάζει με το ??. Οι όγκοι των ελαττωμάτων αντικαθίστανται από μικρές σφαιρικές ενδείξεις, υποδεικνύοντας ότι τα αντικείμενα εξακολουθούν να υπάρχουν στο έργο, αλλά δεν συμμετέχουν πλέον στην προσομοίωση. Αυτό αποκαθιστά την πλευρική ομοιομορφία, διατηρώντας παράλληλα την αρχική γεωμετρία, επιτρέποντάς μας να αλλάξουμε τις διαστάσεις στο επόμενο βήμα χωρίς να εισαγάγουμε νέα φυσικά χαρακτηριστικά.

3. Σύγκριση προσομοιώσεων 3Δ, 2Δ και 1Δ: διαστάσεις, χρόνος εκτέλεσης και αποτελέσματα.

Αφού απενεργοποιηθεί η ανωμαλία, η συσκευή είναι πλέον ομοιόμορφη στις κατευθύνσεις x και z. Αυτό μας επιτρέπει να μελετήσουμε την επίδραση της διάστασης από μόνη της, χωρίς να εισάγουμε νέα φυσικά χαρακτηριστικά. Σε αυτήν την ενότητα, θα εκτελέσετε την ίδια δίοδο GaAs ως πρόβλημα διάχυσης-μεταφοράς σε 3D, 2D και 1D, θα συγκρίνετε πόσο χρόνο διαρκεί κάθε προσομοίωση και θα επαληθεύσετε ότι τα ηλεκτρικά αποτελέσματα παραμένουν αμετάβλητα. Ξεκινήστε ανοίγοντας την καρτέλα "Electrical" και κάνοντας κλικ στο "Electrical mesh" (δείτε ??). Πρώτα, βεβαιωθείτε ότι οι x, y και z είναι όλα ενεργοποιημένα, ώστε η προσομοίωση να είναι πλήρως τρισδιάστατη (??). Με το πλέγμα σε αυτήν την κατάσταση, εκτελέστε την προσομοίωση και παρακολουθήστε την έξοδο στο τερματικό καθώς εξελίσσεται.

Κατά τη διάρκεια της εκτέλεσης, δώστε προσοχή στις πληροφορίες χρονισμού που εμφανίζονται στην οθόνη. Η δεξιά στήλη αναφέρει τον αριθμό των χιλιοστών του δευτερολέπτου που διαρκεί κάθε βήμα της προσομοίωσης (??). Σε έναν τυπικό φορητό υπολογιστή, η τρισδιάστατη περίπτωση θα πρέπει να διαρκεί περίπου δεκάδες έως εκατοντάδες χιλιοστά του δευτερολέπτου ανά βήμα, γεγονός που αντιστοιχεί σε συνολικό χρόνο εκτέλεσης περίπου 30 δευτερολέπτων. Όταν ολοκληρωθεί η προσομοίωση, σημειώστε την τελική τιμή ρεύματος που εμφανίζεται στο τέλος της εκτέλεσης (για παράδειγμα, η γραμμή που αντιστοιχεί στο top = 1,00 V).

Στη συνέχεια, επιστρέψτε στον επεξεργαστή πλέγματος ηλεκτρικού πεδίου και απενεργοποιήστε την κατεύθυνση z. Το πλέγμα θα πρέπει πλέον να αντιπροσωπεύει μια προσομοίωση 2D (y–x) (??). Εκτελέστε ξανά την προσομοίωση και συγκρίνετε την έξοδο της κονσόλας (??). Θα πρέπει να παρατηρήσετε μια σημαντική μείωση στον χρόνο εκτέλεσης, με τους χρόνους ανά βήμα να μειώνονται σε μόλις λίγα χιλιοστά του δευτερολέπτου. Σημαντικό είναι ότι η τελική τιμή ρεύματος που αναφέρεται στο τέλος της εκτέλεσης θα πρέπει να είναι αριθμητικά πανομοιότυπη με την περίπτωση 3D.

Τέλος, απενεργοποιήστε και την κατεύθυνση x, αφήνοντας μόνο την κατεύθυνση μεταφοράς y. Αυτό απλοποιεί το πρόβλημα σε μια καθαρά μονοδιάστατη προσομοίωση διάχυσης-μεταφοράς (??). Εκτελέστε τον επιλυτή για τελευταία φορά και ελέγξτε την έξοδο στην κονσόλα (??). Σε αυτή την περίπτωση, ο χρόνος ανά βήμα είναι ουσιαστικά μηδενικός, ωστόσο το τελικό ρεύμα στην κορυφή, όπου η τάση είναι 1,00 V, είναι το ίδιο με αυτό που παρατηρήθηκε στις δισδιάστατες και τρισδιάστατες προσομοιώσεις.

Αυτή η σύγκριση απεικονίζει μια θεμελιώδη αρχή της μοντελοποίησης συσκευών. Καθώς μειώνεται η διαστατικότητα, η προσομοίωση γίνεται σημαντικά ταχύτερη επειδή ο αριθμός των άγνωστων παραμέτρων μειώνεται ραγδαία (συχνά κατά προσέγγιση, με ρυθμό που σχετίζεται με την τρίτη δύναμη της πυκνότητας του πλέγματος). Ωστόσο, επειδή η συσκευή είναι πλέον ομοιογενής στην πλευρική διάσταση, η μείωση της διαστατικότητας δεν αλλάζει το φυσικό αποτέλεσμα. Λύνετε τις ίδιες εξισώσεις για την ίδια δομή, απλώς εκφρασμένες σε μια πιο αποδοτική αναπαράσταση συντεταγμένων. Το βασικό συμπέρασμα είναι ότι η υψηλότερη διαστατικότητα δεν παράγει αυτόματα διαφορετική ή καλύτερη φυσική συμπεριφορά. Προσθέτει αξία μόνο όταν η ίδια η συσκευή περιέχει πραγματική πλευρική μεταβολή. Μόλις αυτή η μεταβολή αφαιρεθεί, οι τρισδιάστατες, δισδιάστατες και μονοδιάστατες προσομοιώσεις συγκλίνουν στην ίδια ηλεκτρική λύση – διαφέρουν μόνο ως προς το υπολογιστικό κόστος.

Βήμα 5: Παρατηρήστε ότι τα αριθμητικά αποτελέσματα είναι πανομοιότυπα.

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

Αυτό ήταν για το Μέρος Β. Στο Μέρος Γ, θα εξετάσουμε θέματα αριθμητικής σταθερότητας που μπορεί να προκύψουν σε προσομοιώσεις 2D και 3D, και θα αρχίσουμε να μετατρέπουμε αυτή τη δίοδο σε ένα πρόβλημα παρόμοιο με αυτό μιας ηλιακής κυψέλης, εισάγοντας φως και παρατηρώντας πώς η φωτεινότητα αλλάζει την απόκριση της συσκευής.