Summary

רגש החמיריות לחישוב ספקטרום הקליטה אופטי ואלקטרואופטיקה תכונות של אגרגטים מולקולריים ומוצקים

Published: May 27, 2020
doi:

Summary

כאן, אנו מציגים פרוטוקול עבור parametrizing מילטוניאן הדוקה מחייב מחייבת לחישוב ספקטרום קליטה אופטי ותכונות אלקטרואופטיקה של חומרים מולקולריים מתוך העקרונות הראשונים של הקוונטים כימית חישובים.

Abstract

תכנון רציונלי של אגרגטים מולקולריים ומוצקים לאפליקציות אלקטרואופטיקה מסתמך על יכולתנו לחזות את תכונות החומרים הללו בשיטות תאורטיות וחישוביות. עם זאת, מערכות מולקולריות גדולות שבהן הפרעה היא משמעותית מכדי להיחשב במגבלה פרעתיות לא ניתן לתאר באמצעות העקרונות הראשונים כימיה קוונטית או תיאורית הלהקה. מידול רב-גוני הוא גישה מבטיחה להבנת ואופטימיזציה של המאפיינים האופאלקטרוניים של מערכות אלה. היא משתמשת העקרונות הראשונים שיטות קוונטיות כדי לחשב את המאפיינים של מולקולות בודדות, ולאחר מכן בונה מודל החמיריות של אגרגטים מולקולריים או חומרים בתפזורת על בסיס חישובים אלה. במאמר זה אנו מציגים פרוטוקול לבניית מילטוניאן מחייב המייצג את המצבים הנלהבים של חומר מולקולרי בבסיס של הפרקאל הפרקיות: זוגות אלקטרונים-חורים המותאמים למולקולות בודדות העושות את החומר. The מילטוניאן parametrization הציע כאן לחשבון עבור זיווגים מרגש בין מולקולות, כמו גם עבור הקיטוב האלקטרוסטטי של צפיפות האלקטרונים על מולקולה על ידי התפלגות הטעינה על מולקולות המקיפים. כגון מודל הלטונוריאנים ניתן להשתמש כדי לחשב ספקטרום קליטה אופטי ותכונות אלקטרואופטיקה אחרות של אגרגטים מולקולריים ומוצקים.

Introduction

בשני העשורים האחרונים, מוצקים וסרטים העשויים ממולקולות אורגניות צבורים מצאו יישומים מרובים במכשירים אלקטרואופטיקה. להתקנים המבוססים על חומרים מסוג זה יש תכונות אטרקטיביות רבות, כולל משקל קטן, גמישות, צריכת חשמל נמוכה ופוטנציאל לייצור זול באמצעות הדפסת הזרקת דיו. מציג המבוסס על דיודות פולטות אור אורגניים (oleds) מחליפים מציג גבישי נוזלי כמדינה של האמנות עבור טלפונים ניידים, מחשבים נישאים, טלויזיות, ומכשירים אלקטרוניים אחרים1,2,3,4. החשיבות של OLEDs עבור התאורה יישומים צפוי להגדיל בשנים הקרובים4. הביצועים של התקנים אורגניים וולטאית הוא שיפור בהתמדה, עם יעילות המרת החשמל מעל 16% דיווחו לאחרונה עבור הצומת יחיד אורגני תאים סולריים5. חומרים אורגניים יש גם את הפוטנציאל לשבש טכנולוגיות אחרות, כגון תקשורת סיב אופטי, שבו השימוש בהם מאפשר פיתוח של מודולטורים אלקטרו אופטיים עם רוחב מאוד גבוה של 15 thz ומעלה6,7.

האתגר העיקרי במיטוב חומרים מולקולריים של מצב מוצק עבור יישומים ב אלקטרואופטיקה היא שבדרך כלל המאפיינים שלהם תלויים מאוד במבנה הננו של החומר. תהליך הייצור מאפשר להגדיר את ננו מבנה של חומר במידה מסוימת באמצעות טכניקות צמיחה מבוקרת, כגון התצהיר אדים כימיים,8 בתבנית של מולקולות פעילות אופטית על חומר אחר (כלומר, מטריצה פולימר9,10), תרמי ריפוי11,12, וכו ‘. עם זאת, הפרעה ננו-סקאלה הוא מהותי לחומרים המולקולריים ביותר, בדרך כלל לא ניתן לחסל לחלוטין. לכן, ההבנה כיצד ההפרעה משפיעה על המאפיינים של חומר ומציאת דרכים להנדס אותו להשגת ביצועים מיטביים היא חיונית לעיצוב הרציונלי של חומרים אלקטרואופטיקה אורגניים.

מידת ההפרעה בחומרים מולקולריים היא בדרך כלל גדולה מכדי להתייחס אליו כאל מבנה גבישי תקופתי עם מבנה אלקטרוני שניתן לתארו על ידי תיאורית הלהקה. מצד שני, מספר המולקולות שיש לכלול בסימולציה כדי לשכפל את המאפיינים של חומר בתפזורת או סרט הוא גדול מדי כדי להשתמש בעקרונות הראשונים שיטות כימיות קוואנטיות כמו תיאוריה פונקציונלית צפיפות (dft)13,14 ו-הזמן התלוי בצפיפות תיאוריה פונקציונלית (TD-dft)15,16. מולקולות אורגניות עם יישומים אלקטרואופטיקה בדרך כלל יש π גדול יחסית מערכות מצומדות; לרבים מהם יש גם קבוצות תורמים וקבלה. לכידת ההתנהגות הנכונה להעברת טעינה במולקולות כאלה הוא חיוני כדי לחשב את המאפיינים אלקטרואופטיקה שלהם, אבל זה יכול להתבצע רק באמצעות פונקציונליות היברידית לטווח ארוך מתוקן TD-dft17,18,19,20. חישובים המשתמשים פונקציונלים כאלה בקנה מידה סופר באופן ליניארי עם גודל המערכת, כיום, הם מעשיים רק למידול את המאפיינים אלקטרואופטיקה של מולקולות אורגניות בודדות או אגרגטים מולקולריים קטנים שניתן לתאר באמצעות לא יותר מ ~ 104 פונקציות בסיס אטומי. שיטת סימולציה שיכולה לתאר חומרים מופסדר שמורכבים ממספר גדול של ראשון יהיה מאוד שימושי עבור דוגמנות מערכות אלה.

היקף האינטראקציות הבינמולקולריות בחומרים מולקולריים הוא לעתים קרובות דומה או קטן יותר מסדר הווריאציה בפרמטרים האנרגטיים (כגון אנרגיות המדינה או אנרגיות עירור) בין מולקולות בודדות העושות את החומר. במקרים כאלה, מידול רב-גוני הוא הגישה המבטיחה ביותר להבנת ואופטימיזציה של המאפיינים האופאלקטרוניים של מערכות מולקולריות גדולות בסדר גודל21,22,23. גישה זו משתמשת העקרונות הראשונים שיטות קוונטיות (בדרך כלל DFT ו TD-DFT) כדי לחשב במדויק את המאפיינים של מולקולות בודדות שמרכיבים את החומר. המילטוניאן של מדגם חומר גדול מספיק כדי לייצג את החומר המולקולרי בצובר (אולי, על ידי העסקת גבולות תקופתיים) הוא נבנה אז באמצעות הפרמטרים שחושבו עבור מולקולות בודדות. מילטוניאן זה יכול לשמש כדי לחשב את הפרמטרים אלקטרואופטיקה של צבירה מולקולרית גדול, סרט דק, או חומר מולקולרי בצובר.

מודלים מדגם הינם מחלקת משנה של דגמים מרובי-סקאלה, שבהם המדינות המתרגשות של חומר מולקולרי מיוצגות בבסיס של מרגש: צמדי אלקטרון-חור, המאוגדים באטרקציה של קולון24,25. למידול תהליכי מדינה נרגשים רבים, מספיקה לכלול רק את פרנקל מרגש26, שם האלקטרון והחור מותאמים לאותה מולקולה. העברת מטען, היכן שהאלקטרון והחור מותאמים למולקולות שונות, ייתכן שיהיה צורך לכלול אותם במקרים מסוימים (למשל, כאשר מידול ההפרדה של התביעות במערכות התורמות)27,28. למרות ההתרגשות מודלים הם מודלים multiscale שניתן לparametrized באמצעות חישובי העיקרון הראשון בלבד על מולקולות בודדות, הם עדיין בחשבון עבור אינטראקציות בין-מולקולריות. שני סוגי האינטראקציה העיקריים שהם יכולים להסביר הם (a) מצמדים מרגש בין מולקולות המאפיינות את היכולת של הרגש להתאפיין על פני או להעביר בין מולקולות ו-(ב) הקיטוב האלקטרוסטטי של צפיפות האלקטרונים על מולקולה על ידי התפלגות האישום על מולקולות המקיפים. בעבר הצגנו כי שני הגורמים הללו חשובים למידול התכונות האופטיות והאלקטרו-אופטיות של אגרגטים מולקולריים, כגון ספקטרום הקליטה האופטי29 והראשונה היפרפולביות30.

במאמר זה, אנו מציגים פרוטוקול עבור parametrizing, מודלים שניתן להשתמש בהם כדי לחשב את ספקטרום אופטי ותכונות אלקטרואופטיקה אחרות של אגרגטים מולקולריים גדולים וחומרים מולקולריים בצובר. מילטוניאן ההתרגשות הוא הניח להיות מחייב מילטוניאן24,25,

Equation 1

כאשר εi הוא האנרגיה עירור שלהמולקולה האניבחומר, bij הוא זיווג מרגש בין האיי והמולקולות jth , âi ו- âאני הם מפעילי הבריאהוהשמדה , בהתאמה, עבור המדינה נרגש על המולקולההiבחומר. הפרמטרים המילטוניאן ההנדוניים נמצאים באמצעות חישובי TD-DFT המבוצעים על מולקולות בודדות העושות את החומר. בחישובי TD-DFT אלה, התפלגות החיוב על כל המולקולות האחרות בחומר מיוצגת על-ידי הטבעה אלקטרוסטטית של האשמות בנקודות האטום לחשבון בנוגע לקיטוב אלקטרוסטטי של דחיסות אלקטרונית של מולקולה. האנרגיות עירור, εi, עבור מולקולות בודדות נלקחים ישירות הפלט TD-dft החישוב. מצמדים מרגש, bij, בין מולקולות מחושבים באמצעות שיטת הקוביה בצפיפות שיטה31, עם הקרקע-להתרגש צפיפות מעבר המדינה עבור מולקולות אינטראקציה שנלקחו הפלט של חישוב TD-dft ב- גאוס32 ו לאחר עיבוד באמצעות מנתח הפונקציה multiwfn משולבים-גל33. להדמיית המאפיינים של מוצקים מולקולריים בצובר, ניתן להחיל תנאי גבול תקופתיים על מילטוניאן.

הפרוטוקול הנוכחי מחייב שלמשתמש תהיה גישה ל-32 הגאוסו- multiwfn33 תוכניות. הפרוטוקול נבדק באמצעות 16 , מהדורה B1 ו- multiwfn גירסה 3.3.8, אך יש לעבוד גם עבור גירסאות עדכניות אחרות של תוכניות אלה. בנוסף, הפרוטוקול משתמש בכלי שירות מותאם אישית ב-C++ ובמספר של פיתון מותאם אישית 2.7 וקבצי script של Bash, קוד המקור שעבורו מסופק תחת הרישיון הציבורי הכללי של GNU (גירסה 3) ב-https://github.com/kocherzhenko/ExcitonicHamiltonian. החישובים נועדו להתבצע במחשב שמפעיל מערכת הפעלה ממשפחת Unix/Linux.

Protocol

1. פיצול מערכת multi-מולקולרית למולקולות בודדות ליצור את המבנה של המערכת שעבורה המילטוניאן הפוניים צריך להיות בנוי בפורמט הקובץ המולקולרי של Tripos MOL2. מבנה זה יכול להיות תמונה של דינמיקה מולקולרית או סימולציה מונטה קרלו של המערכת. אם כל המולקולות במערכת מורכבות מאותו מספר של אטומים,…

Representative Results

בסעיף זה אנו מציגים תוצאות הנציג למחשוב את ספקטרום הקליטה האופטית של צבירה של שישה מולקולות YLD 124, המוצג באיור 3a, שבו מבנה הצבירה הושגה מסימולציה מונטה קרלו גסה. YLD 124 הוא כרומטופור-העברת-מטען טיפוסי המורכבת מקבוצה של אלקטרון-תרומת דיאתיל אמין עם הגנה מפני הקבוצות של tert-בוטיל…

Discussion

השיטה המוצגת כאן מאפשרת התאמות אישיות מרובות. לדוגמה, ניתן לשנות את הפרמטרים של חישובי DFT ו-TD-DFT, כולל הצפיפות הפונקציונלית, ערכת הבסיס וההגדרה הספציפית של חיובי הנקודה האטומית.

שימוש בפונקציונליות מתוקנת של טווח ארוך, כגון ωB97X, ωB97XD או ωPBE, מומלץ על מנת להשיג צפיפויות מעבר סבי…

Disclosures

The authors have nothing to disclose.

Acknowledgements

אנו מודים לד ר אנדריאס טילמון (אוק רידג ‘ המעבדה הלאומית), ד ר לואיס ג’ונסון (אוניברסיטת וושינגטון), וד ר ברוס רובינסון (אוניברסיטת וושינגטון) לפיתוח התוכנית לסימולציות גסות מונטה קרלו ששימשו להפקת מבנה המערכת המולקולרית שהוצגה בסעיף תוצאות הנציג. A.A.K. ו-P.F.G. נתמכים על ידי פרס מחקר משותף מהמכללה למדעים, מפרץ מזרח לאטלסיב. M.H. נתמך על ידי מלגת חלוץ לנצח מהמרכז לחקר הסטודנטים, מפרץ מזרח לאטלסיב. C.M.I. ו-אס. פי. או נתמכים על ידי משרד ההגנה של ארה ב (הצעה 67310-CH-נציג) תחת משרד חיל האוויר של המחלקה למחקר מדעי חומרים אורגניים.

Materials

Gaussian 16, revision B1
Multiwfn version 3.3.8
GNU compiler collection version 9.2
python 2.7.0

References

  1. Tsujimura, T. . OLED Display Fundamentals and Applications, 2nd Ed. , (2017).
  2. Barnes, D. LCD or OLED: Who Wins. SID Symposium Digest of Technical Papers. 44 (1), 26-27 (2013).
  3. Mizukami, M., et al. Flexible Organic Light-Emitting Diode Displays Driven by Inkjet-Printed High-Mobility Organic Thin-Film Transistors. IEEE Electron Device Letters. 39 (1), 39-42 (2018).
  4. Koden, M. . OLED Displays and Lighting. , (2017).
  5. Fan, B., et al. Achieving Over 16% Efficiency for Single-Junction Organic Solar Cells. Science China Chemistry. 62 (6), 746-752 (2018).
  6. Dalton, L. R., Gunter, P., Jazbinsek, M., Kwon, O. P., Sullivan, P. A. . Organic Electro-Optics and Photonics: Molecules, Polymers, and Crystals. , (2015).
  7. Robinson, B. H., et al. Optimization of Plasmonic-Organic Hybrid Electro-Optics. Journal of Lightwave Technology. 36 (21), 5036-5047 (2018).
  8. Yu, D., Yang, Y. Q., Chen, Z., Tao, Y., Liu, Y. F. Recent Progress on Thin-Film Encapsulation Technologies for Organic Electronic Devices. Optics Communications. 362 (1), 43-49 (2016).
  9. Wanapun, D., Hall, V. J., Begue, N. J., Grote, J. G., Simpson, G. J. DNA-Based Polymers as Chiral Templates for Second-Order Nonlinear Optical Materials. ChemPhysChem. 10 (15), 2674-2678 (2009).
  10. Siao, Y. Y., et al. Orderly Arranged NLO Materials on Exfoliated Layered Templates Based on Dendrons with Alternating Moieties at the Periphery. Polymer Chemistry. 4 (9), 2747-2759 (2013).
  11. Sepeai, S., Salleh, M. M., Yahaya, M., Umar, A. A. Improvement of White Organic Light Emitting Diode Performances by an Annealing Process. Thin Solid Films. 517 (16), 4679-4683 (2009).
  12. Mao, G., et al. Considerable Improvement in the Stability of Solution Processed Small Molecule OLED by Annealing. Applied Surface Science. 257 (17), 7394-7398 (2011).
  13. Parr, R. G., Yang, W. . Density Functional Theory of Atoms and Molecules. , (1989).
  14. Dreizlerm, R. M., Gross, E. K. U. . Density Functional Theory: An Approach to the Quantum Many-body Problem. , (1990).
  15. Burke, K., Werschnik, J., Gross, E. K. U. Time-Dependent Density Functional Theory: Past, Present, and Future. Journal of Chemical Physics. 123 (6), 062206 (2005).
  16. Ullrich, C. . Time-Dependent Density-Functional Theory: Concepts and Applications. , (2011).
  17. Vydrov, O. A., Scuseria, G. E. Assessment of a Long-Range Corrected Hybrid Functional. Journal of Chemical Physics. 125 (23), 234109 (2006).
  18. Tawada, Y., Tsuneda, T., Yanagisawa, S. A Long-Range-Corrected Time-Dependent Density Functional Theory. Journal of Chemical Physics. 120 (18), 5425 (2004).
  19. Rohrdanz, M. A., Herbert, J. M. Simultaneous Benchmarking of Ground- and Excited-State Properties with Long-Range-Corrected Density Functional Theory. Journal of Chemical Physics. 129 (3), 034107 (2008).
  20. Autschbach, J. Charge-Transfer Excitations and Time-Dependent Density Functional Theory: Problems and Some Proposed Solutions. ChemPhysChem. 10 (11), 1757-1760 (2008).
  21. Nelson, J., Kwiatkowski, J. J., Kirkpatrick, J., Frost, J. M. Modeling Charge Transport in Organic Photovoltaic Materials. Accounts of Chemical Research. 42 (11), 1768-1778 (2009).
  22. Walker, A. B. Multiscale Modeling of Charge and Energy Transport in Organic Light-Emitting Diodes and Photovoltaics. Proceedings of the IEEE. 97 (9), 1587-1596 (2009).
  23. Wang, L., Li, Q., Shuai, Z., Chenc, L., Shic, Q. Multiscale Study of Charge Mobility of Organic Semiconductor with Dynamic Disorders. Physical Chemistry Chemical Physics. 12 (13), 3309-3314 (2010).
  24. Davydov, A. S. . Theory of Molecular Excitons. , (1971).
  25. Agranovich, V. M. Excitations in Organic Solids. International Series of Monographs on Physics. 142, (2008).
  26. Frenkel, J. On the Transformation of Light into Heat in Solids. I. Physical Review. 37 (1), 17-44 (1931).
  27. Kocherzhenko, A. A., Lee, D., Forsuelo, M. A., Whaley, K. B. Coherent and Incoherent Contributions to Charge Separation in Multichromophore Systems. Journal of Physical Chemistry C. 119 (14), 7590-7603 (2015).
  28. Lee, D., Forsuelo, M. A., Kocherzhenko, A. A., Whaley, K. B. Higher-Energy Charge Transfer States Facilitate Charge Separation in Donor-Acceptor Molecular Dyads. Journal of Physical Chemistry C. 121 (24), 13043-13051 (2017).
  29. Kocherzhenko, A. A., Sosa Vazquez, X. A., Milanese, J. M., Isborn, C. M. Absorption Spectra for Disordered Aggregates of Chromophores Using the Exciton Model. Journal of Chemical Theory and Computation. 13 (8), 3787-3801 (2017).
  30. Kocherzhenko, A. A., et al. Unraveling Excitonic Effects for the First Hyperpolarizabilities of Chromophore Aggregates. Journal of Physical Chemistry C. 123 (22), 13818-13836 (2019).
  31. Krueger, B., Scholes, G., Fleming, G. Calculation of Couplings and Energy-Transfer Pathways between the Pigments of LH2 by the ab Initio Transition Density Cube Method. Journal of Physical Chemistry B. 102 (27), 5378-5386 (1998).
  32. Frisch, M. J., et al. Gaussian 16, Revision B.01. Gaussian, Inc. , (2016).
  33. Lu, T., Chen, F. Multiwfn: A Multifunctional Wavefunction Analyzer. Journal of Computational Chemistry. 33 (5), 580-592 (2012).
  34. Chai, J. D., Head-Gordon, M. Systematic Optimization of Long-Range Corrected Hybrid Density Functionals. Journal of Chemical Physics. 128 (8), 084106 (2008).
  35. Hehre, W., Ditchfield, R., Pople, J. Self-Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian-Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules. Journal of Chemical Physics. 56 (5), 2257-2261 (1972).
  36. Hariharan, P., Pople, J. The Influence of Polarization Functions on Molecular Orbital Hydrogenation Energies. Theoretica chimica acta. 28 (3), 213-222 (1973).
  37. Breneman, C. M., Wiberg, K. B. Determining Atom-Centered Monopoles from Molecular Electrostatic Potentials. The Need for High Sampling Density in Formamide Conformational Analysis. Journal of Computational Chemistry. 11 (3), 361-373 (1990).
  38. Mulliken, R. S. Electronic Population Analysis on LCAO-MO MolecularWave Functions. I. Journal of Chemical Physics. 23 (10), 1833-1840 (1955).
  39. Jen, A., et al. Exceptional Electro-Optic Properties through Molecular Design and Controlled Self-Assembly. Proceedings of SPIE. 5935, 593506 (2005).
  40. Hirata, S., Head-Gordon, M. Time-Dependent Density Functional Theory Within the Tamm-Dancoff Approximation. Chemical Physics Letters. 314 (3-4), 291-299 (1999).
  41. Randolph, K. A., Myers, L. L. . Basic Statistics in Multivariate Analysis. , 11-34 (2013).
  42. Garrett, K., et al. Optimum Exchange for Calculation of Excitation Energies and Hyperpolarizabilities of Organic Electro-optic Chromophores. Journal of Chemical Theory and Computation. 10 (9), 3821-3831 (2014).
  43. Sekino, H., Maeda, Y. Polarizability and Second Hyperpolarizability Evaluation of Long Molecules by the Density Functional Theory with Long-Range Correction. Journal of Chemical Physics. 126 (1), 014107 (2007).
  44. Johnson, L. E., Dalton, L. R., Robinson, B. H. Optimizing Calculations of Electronic Excitations and Relative Hyperpolarizabilities of Electrooptic Chromophores. Accounts of Chemical Research. 47 (11), 3258-3265 (2014).
  45. Lee, J., et al. Molecular Mechanics Simulations and Improved Tight-Binding Hamiltonians for Artificial Light Harvesting Systems: Predicting Geometric Distributions, Disorder, and Spectroscopy of Chromophores in a Protein Environment. Journal of Physical Chemistry B. 122 (51), 12292-12301 (2018).
  46. Bellinger, D., Pflaum, J., Brüning, C., Engel, V., Engels, B. The Electronic Character of PTCDA Thin Films in Comparison to Other Perylene-Based Organic Semi-conductors: Ab Initio-, TD-DFT and Semi-Empirical Computations of the Opto-Electronic Properties of Large Aggregates. Physical Chemistry Chemical Physics. 19 (3), 2434 (2017).
  47. Zuehlsdorff, T. J., Isborn, C. M. Combining the Ensemble and Franck-Condon Approaches for Calculating Spectral Shapes of Molecules in Solution. The Journal of Chemical Physics. 148 (2), 024110 (2018).
  48. Zuehlsdorff, T. J., Isborn, C. M. Modeling Absorption Spectra of Molecules in Solution. International Journal of Quantum Chemistry. 119 (1), 25719 (2019).
  49. Plötz, P. A., Megow, J., Niehaus, T., Kühn, O. All-DFTB Approach to the Parametrization of the System-Bath Hamiltonian Describing Exciton-Vibrational Dynamics of Molecular Assemblies. Journal of Chemical Theory and Computation. 14 (10), 5001-5010 (2018).
  50. Tillack, A., Johnson, L., Eichinger, B., Robinson, B. H. Systematic Generation of Anisotropic Coarse-Grained Lennard-Jones Potentials and Their Application to Ordered Soft Matter. Journal of Chemical Theory and Computation. 12 (9), 4362-4374 (2016).

Play Video

Cite This Article
Kocherzhenko, A. A., Shedge, S. V., Germaux, P. F., Heidarian, M., Isborn, C. M. Excitonic Hamiltonians for Calculating Optical Absorption Spectra and Optoelectronic Properties of Molecular Aggregates and Solids. J. Vis. Exp. (159), e60598, doi:10.3791/60598 (2020).

View Video