Μετρήσεις CELIV σε ηλιακά κύτταρα περοβσκίτη: Μέρος α
1. Επισκόπηση
Σε αυτό το εκπαιδευτικό υλικό, θα ξεκινήσουμε παρουσιάζοντας τις βασικές έννοιες της τεχνικής «Εξαγωγής Φορτίου με Γραμμική Αύξηση της Τάσης» (CELIV), περιγράφοντας πώς χρησιμοποιείται για τη διερεύνηση της κινητικότητας και της πυκνότητας των φορέων σε ημιαγωγούς λεπτών υμενίων. Αφού εδραιώσουμε το θεωρητικό υπόβαθρο, θα προχωρήσουμε γρήγορα στα πρακτικά βήματα, δείχνοντας πώς να διαμορφώσετε και να εκτελέσετε μια προσομοίωση CELIV στο OghmaNano. Αυτός ο συνδυασμός θεωρίας και πρακτικής προσομοίωσης θα παρέχει μια σαφή κατανόηση της μεθόδου και τα εργαλεία για να την εφαρμόσετε στις δικές σας συσκευές.
2. Ιστορικό του CELIV.
Η τεχνική Εξαγωγής Φορτίου μέσω Γραμμικής Αύξησης της Τάσης (CELIV) είναι μια πειραματική μέθοδος που χρησιμοποιείται για τη μελέτη της κινητικότητας και της πυκνότητας των φορέων φορτίου σε οργανικά και περοβσκιτικά ηλιακά κύτταρα. Σε μια τυπική διάταξη, ένα ηλιακό κύτταρο συνδέεται με μια γεννήτρια σήματος που εφαρμόζει μια γραμμικά αυξανόμενη αντίστροφη τάση, ενώ η προκύπτουσα μεταβατική τιμή του ρεύματος καταγράφεται σε ένα οσκιλλόσκοπιο. Προαιρετικά, μπορεί να προστεθεί συνεχής φωτισμός ή ένας σύντομος παλμός λέιζερ για τη δημιουργία φορέων φορτίου. Ένα σχηματικό διάγραμμα αυτής της διαμόρφωσης εμφανίζεται στο ??.
Όταν εφαρμόζεται η ράμπα τάσης, το μετρούμενο ρεύμα αποτελείται από δύο διακριτές συνιστώσες: μια ορθογώνια βασική γραμμή που προέρχεται από τη γεωμετρική χωρητικότητα της συσκευής και μια επιπλέον κορυφή που προέρχεται από την εξαγωγή κινητών φορέων φορτίου στον ημιαγωγό. Ο όρος της γεωμετρικής χωρητικότητας είναι πάντα παρών, ενώ η κορυφή που σχετίζεται με τους φορείς παρέχει άμεσες πληροφορίες σχετικά με την κινητικότητα και την πυκνότητα των φορέων. Ένα παράδειγμα αυτών των μεταβατικών τάσεων και ρευμάτων φαίνεται στο ??.
Η κινητικότητα των φορέων μπορεί να προσδιοριστεί άμεσα από το παροδικό σήμα CELIV, μετρώντας τον χρόνο στον οποίο το ρεύμα εξαγωγής φτάνει στο μέγιστο. Στην απλούστερη περίπτωση ομοιόμορφης κατανομής των φορέων και αμελητέας ανασυνδυασμού, η κινητικότητα δίνεται από την εξίσωση CELIV: \[ \mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}} \], όπου \(d\) είναι το πάχος του ενεργού στρώματος, \(A\) είναι ο ρυθμός αύξησης της τάσης (\(A = \mathrm{d}V/\mathrm{d}t\)) και \(t_{\text{max}}\) είναι ο χρόνος στον οποίο εμφανίζεται η μέγιστη τιμή του ρεύματος. Αυτή η σχέση συνδέει το πειραματικά παρατηρούμενο παροδικό σήμα με τη θεμελιώδη παράμετρο μεταφοράς που μας ενδιαφέρει. ( Αρχική εξαγωγή από τον Juska)
3. Περιορισμοί του CELIV.
Η εξαγωγή της εξίσωσης κινητικότητας CELIV βασίζεται σε διάφορες απλοποιητικές υποθέσεις. Πρώτον, υποτίθεται ότι μόνο ένας τύπος φορέα φορτίου κυριαρχεί στη μεταφορά, έτσι ώστε μόνο ένα είδος φορέα να κινείται όταν εφαρμόζεται η κλίση τάσης. Δεύτερον, υποτίθεται ότι οι φορείς είναι ομοιόμορφα κατανεμημένοι και γλιστρούν ομαλά έξω από τη συσκευή με τρόπο παρόμοιο με την ιδανική περίπτωση που εμφανίζεται στην αριστερή πλευρά του ??. Τρίτον, οι διαδικασίες ανασύνθεσης θεωρούνται αμελητέες κατά την εξαγωγή, έτσι ώστε να μην υπάρχει σημαντική απώλεια φορέων. Στην πραγματικότητα, αυτές οι συνθήκες σπάνια ικανοποιούνται. Όπως απεικονίζεται στη δεξιά πλευρά του ??, οι φορείς φορτίου συχνά καταλαμβάνουν ρηχές, μέτριες ή βαθιές παγίδες και η καθυστερημένη απελευθέρωσή τους παράγει διευρυμένα ή παραμορφωμένα μεταβατικά φαινόμενα εξαγωγής. Όπως και με πολλές πειραματικές τεχνικές, η κινητικότητα που εξάγεται με τη χρήση της CELIV θα πρέπει επομένως να θεωρείται ως μια προσέγγιση και όχι ως η πραγματική μικροσκοπική τιμή, και στην πραγματικότητα η φαινομενική κινητικότητα μπορεί να εξελιχθεί κατά τη διάρκεια του ίδιου του μεταβατικού φαινομένου ( Μια κριτική επισκόπηση της CELIV που συζητά τους περιορισμούς της). Αυτός ο περιορισμός είναι ιδιαίτερα σημαντικός σε διαταραγμένους οργανικούς ημιαγωγούς, όπου η ενεργειακή διαταραχή και η παγίδευση επηρεάζουν σημαντικά τη μεταφορά. Αντίθετα, τα υλικά περοβσκίτη συχνά παρουσιάζουν υψηλότερες κινητικότητες και χαμηλότερες πυκνότητες παγίδων, καθιστώντας την τυπική ανάλυση CELIV πιο αξιόπιστη και ευκολότερη στην ερμηνεία.
💡 Δοκιμάστε να αξιολογήσετε την κατανόησή σας: Προσπαθήστε να εφαρμόσετε τη θεωρία CELIV σε μια απλή περίπτωση χρήσης.
Ας υποθέσουμε ότι έχετε μια συσκευή με πάχος d = 200 nm, ρυθμό αύξησης της τάσης A = 2 × 106 V/s και μέγιστη τιμή ρεύματος που παρατηρείται στο tmax = 5 µs. Χρησιμοποιώντας την εξίσωση κινητικότητας CELIV, εκτιμήστε την κινητικότητα των φορέων.
Εμφάνισε την απάντηση.
Η κινητικότητα του CELIV δίνεται από την εξίσωση \[
\mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}}
\]. Αντικαθιστώντας τις τιμές: d = 200 × 10-9 m, A = 2 × 106 V/s, tmax = 5 × 10-6 s.
\[ \mu = \frac{2 (200 × 10^{-9})^{2}}{3 (2 × 10^{6}) (5 × 10^{-6})^{2}} \approx 1.1 × 10^{-4} \; \text{cm}^2/\text{Vs} \]
Αυτό το απλό παράδειγμα δείχνει πώς η χρονική στιγμή της μέγιστης τιμής σε ένα παροδικό φαινόμενο CELIV μπορεί να μετατραπεί σε μια εκτίμηση της κινητικότητας.
4: Εκτέλεση μιας προσομοίωσης CELIV στο OghmaNano.
Κάντε κλικ στο "Νέα προσομοίωση". Αυτό θα ανοίξει τη βιβλιοθήκη με τους διαθέσιμους τύπους συσκευών, όπως φαίνεται στην εικόνα ??. Κάντε διπλό κλικ στο "Κελιά περοβσκίτη" (που είναι επισημασμένο με κόκκινο) για να ανοίξετε τον φάκελο με τα παραδείγματα περοβσκίτη. Θα δείτε μια λίστα με τις προκαθορισμένες προσομοιώσεις, συμπεριλαμβανομένων των "Συσκευή MAPbI₃", "Ηλιακή κυψέλη περοβσκίτη" και ειδικών επιδείξεων CELIV, όπως φαίνεται στην εικόνα ??. Για αυτό το σεμινάριο, επιλέξτε το παράδειγμα "Ηλιακή κυψέλη περοβσκίτη – CELIV". Όταν σας ζητηθεί, αποθηκεύστε την προσομοίωση σε έναν φάκελο στον οποίο έχετε δικαιώματα εγγραφής.
Αφού επιλέξετε το παράδειγμα ηλιακού στοιχείου περοβσκίτη – CELIV, ανοίγει το κύριο παράθυρο προσομοίωσης (δείτε ??). Για να ξεκινήσετε τον υπολογισμό, κάντε κλικ στην επιλογή "Εκτέλεση προσομοίωσης" (μπλε εικονίδιο αναπαραγωγής) ή πατήστε το πλήκτρο F9. Στη συνέχεια, το OghmaNano θα επιλύσει τις εξισώσεις διάχυσης-μεταφοράς που εξαρτώνται από το χρόνο και θα δημιουργήσει την παροδική καμπύλη CELIV.
jv.csv (δεδομένα καμπύλης JV), time_j.csv (ρεύμα εξαγωγής έναντι χρόνου), time_v.csv (εφαρμοσμένη τάση έναντι χρόνου) και optical_output (κατανομές πεδίου). Κάντε διπλό κλικ σε οποιοδήποτε αρχείο για να δείτε τα αντίστοιχα αποτελέσματα στα ενσωματωμένα εργαλεία γραφικών παραστάσεων.Βήμα 5: Δείτε τα αποτελέσματα του CELIV.
Άνοιξε το Έξοδος ταμπλέτα??) και κάντε διπλό κλικ. time_v.csvΑυτό απεικονίζει το πρόγραμμα της εφαρμοζόμενης τάσης – την κλίση CELIV – όπως φαίνεται στο
??Στη συνέχεια, κάντε διπλό κλικ. time_j.csv για να εμφανιστεί η μεταβατική τιμή του ρεύματος εξαγωγής, όπως στο ??Θα παρατηρήσετε ότι το παροδικό στοιχείο έχει το οικείο σχήμα CELIV, αλλά φαίνεται να είναι **ανεστραμμένο**: σύμφωνα με τη σύμβαση σημάτων στο OghmaNano, το ρεύμα που εξέρχεται από τη συσκευή είναι αρνητικό. Όταν η κλίση στο
?? εφαρμόζεται,
η φόρτιση εξάγεται από τη συσκευή, έτσι η συσκευή παρέχει αρνητικό ρεύμα (μέγιστη τιμή εξαγωγής).
Η μεταγενέστερη μέγιστη τιμή με αντίθετο πρόσημο (θετική στην
??) συμβαίνει όταν η τάση πόλωσης επανέρχεται και το φορτίο επιστρέφει στη συσκευή, αφού η παροδική κατάσταση έχει επιστρέψει στην αρχική της κατάσταση.
Το διάγραμμα στο ?? απεικονίζει το παροδικό φαινόμενο CELIV, αφού πολλαπλασιάσουμε το προσομοιωμένο ρεύμα με το -1, ώστε να ταιριάζει με την συμβατική μορφή παρουσίασης που χρησιμοποιείται στα πειράματα. Σε αυτή τη μορφή, η βασική γραμμή αντιστοιχεί στη γεωμετρική χωρητικότητα της συσκευής, ενώ η απότομη θετική κορυφή υποδεικνύει την εξαγωγή κινητών φορέων φορτίου κατά τη διάρκεια της εφαρμοζόμενης ράμπας τάσης. Το μέγεθος και ο χρονισμός αυτής της κορυφής είναι οι βασικές ποσότητες που χρησιμοποιούνται για τον προσδιορισμό της κινητικότητας και της πυκνότητας των φορέων. Η παρουσίαση του παροδικού φαινομένου σε αυτή τη συμβατική μορφή διευκολύνει τη σύγκριση των προσομοιώσεων OghmaNano απευθείας με τις δημοσιευμένες μετρήσεις CELIV. Για να εξαχθεί η κινητικότητα, προσδιορίστε τον χρόνο στον οποίο το ρεύμα φτάνει στο μέγιστο (\(t_\text{max}\)) και αντικαταστήστε αυτή την τιμή στην εξίσωση κινητικότητας CELIV \(\mu = \tfrac{2d^2}{3At_\text{max}^2}\), όπου \(d\) είναι το πάχος του ενεργού στρώματος και \(A\) είναι ο ρυθμός της εφαρμοζόμενης ράμπας τάσης.
div style="clear: both;">🧪 Παράδειγμα: Διαβάστε τον χρόνο κορυφής του σήματος CELIV από το ?? και υπολογίστε την κινητικότητα. Υποθέστε ένα πάχος της ενεργού στρώσης d = 600 nm, μια ταχύτητα αύξησης της τάσης A = 5 × 105 V/s και έναν μετρούμενο χρόνο κορυφής tmax = 3 µs. Ποια τιμή κινητικότητας προκύπτει;
Εμφάνισε την απάντηση.
Χρησιμοποιήστε τη σχέση CELIV μεταξύ των \[ \mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}} \], \(d = 600\times10^{-9}\,\mathrm{m}\), \(A = 5\times10^{5}\,\mathrm{V\,s^{-1}}\) και \(t_{\text{max}} = 3\times10^{-6}\,\mathrm{s}\).
\[ \mu = \frac{2(600\times10^{-9})^{2}}{3\,(5\times10^{5})\,(3\times10^{-6})^{2}} \approx 5.3\times10^{-8}\ \mathrm{m^{2}\,V^{-1}\,s^{-1}} \] Μετατροπή σε \(\mathrm{cm^{2}\,V^{-1}\,s^{-1}}\) (\(1\ \mathrm{m^{2}} = 10^{4}\ \mathrm{cm^{2}}\)): \[ \mu \approx 5.3\times10^{-4}\ \mathrm{cm^{2}\,V^{-1}\,s^{-1}}. \]
Αυτή η τιμή βρίσκεται εντός του τυπικού εύρους που αναφέρεται για διαταραγμένους οργανικούς ημιαγωγούς. (Εάν η τιμή \({t_\text{max}}\) αλλάζει ανάλογα με την ταχύτητα αύξησης ή τον φωτισμό, θα δείτε ότι η τιμή \(\mu\) μετατοπίζεται ανάλογα.)
👉 Επόμενα βήματα:
- Προχωρήστε στην επόμενη ενότητα της προσομοίωσης CELIV — Μέρος Β — όπου θα ρυθμίσουμε την ταχύτητα αύξησης της τάσης, θα πειραματιστούμε με την ένταση του φωτός και θα εξερευνήσουμε πώς η αλλαγή της κινητικότητας των ηλεκτρονίων/οπών επηρεάζει το μεταβατικό φαινόμενο του CELIV.
- Επιστρέψτε στην οδηγία προσομοίωσης ηλιακών κυψελών περοβσκίτη για να δείτε την πλήρη επισκόπηση του εκπαιδευτικού προγράμματος.