מודלי הפרוטוקול, תופעות תרמיות וקוונטיות בקטליזה הטרוגנית בשלב נוזלי. זוהי הראשונה, לדעתנו, לשלב מכניקת קוונטים עם דגימה מלאה של סביבה נוזלית מפורשת. תצורות של מולקולות נוזליות הנוצרות מפרוטוקול זה מייצגות את אלה הצפויות בתנאי תגובה ממשיים ו ניתן להשתמש בהן כדי לחקור תופעות ברמה מולקולרית התלויות בסידורים המרחביים של מולקולות.
התצורות של מולקולות נוזליות הנוצרות על ידי פרוטוקול זה מספקות תובנה לגבי התפקידים שממס ממלא על קטליזה הטרוגנית בשלב נוזלי. אם אתה מנסה פרוטוקול זה בפעם הראשונה, אני ממליץ לך קודם לוודא שיש לך גישה ובאפשרותך להפעיל את תוכנת VASP, MCPliQ, VMD ו- LAMMPS. הדגימה של ההליך תהיה טיאנג'ון שי, סטודנטית לתואר שני מהמעבדה שלי.
לאחר יצירת מבנה הפילוח, צור את קבצי הקלט של LAMMPS עבור הדמיית NPT וצייד את נפח התא באמצעות FFMD. העתק את קובץ הקלט של LAMMPS לספריה העובדת. ערכו את משתנה הקבוצה בשורה 34 כדי לציין את אינדקסי סוגי האטום עבור אטומי החמצן והמים של מימן המים.
במשתנה הקבוצה בשורה 35, כדי לציין את אינדקסי סוג האטום עבור אטומי הפלטינה והספוג. בשורה 17 של קובץ הקלט, ערוך את משתנה שלב הריצה כדי להגדיר את משך הדמיית NPT מספיק זמן כדי לכלול הפעלת שיווי משקל והפעלת ייצור. בצע את תוכנת LAMMPS, על-ידי הקלדת הפקודה בממשק שורת הפקודה, הכולל את המידע של מספר ליבות ה- CPU לשימוש ואת שם קובץ ההפעלה LAMMPS.
מזעור האנרגיה מחדד את תצורת מולקולת המים ואחריו סימולציית FFMD המבוצעת במספר קבוע של מולקולות מים, נפח וטמפרטורה כדי להביא את המים לטמפרטורת הסימולציה. סימולציה נוספת של FFMD פועלת לאחר מכן במספר קבוע של מולקולות מים, לחץ וטמפרטורה כדי לקבוע את הגובה הנכון פיזית של תיבת הסימולציה. קבצי הפלט ישמשו מאוחר יותר.
לאחר סימולציית NPT, להתוות את הגובה של supercell נגד הזמן. הנקודה שבה הוא יורד לערך מצב יציב היא הנקודה בסימולציית NPT כאשר הייצור יכול להתחיל. ודא שיווי משקל של סימולציית NPT על ידי הבטחת כי התנודות בגובה של supercell הם מינימליים או התכנסו לערך יציב.
אם מתרחשות תנודות גדולות, פתח את הקלט. לצייד את הקובץ ולהקטין את timestep על קו 92 כדי לחדש את תצורת מולקולת המים ולבצע את תוכנת LAMMPS כמו קודם. כדי להתחיל, הקלד בממשק שורת הפקודה כדי להפעיל את קובץ ה- Script.
קובץ Script זה מפיק את גובה supercell הממוצע מחלק הפעלת הייצור של הדמיית NPT לקובץ TXT. קובץ ה- Script קורא את אורך ממד z של התא במרווחים של 1, 000 femtosecond, שהוא מרווח הזמן המוגדר כברירת מחדל להדפסת מידע ב- LAMMPS. אם רצוי מרווח הדפסה אחר, ניתן לשנות אותו על-ידי עריכת קו 20 של get_npt_lz.
py סקריפט ו קו 16 של קלט LAMMPS. לצייד את הקובץ. התסריט מזהה ומשליך את שני הננו-שניות הראשונות בשווי ערכי lz, מכיוון שהן מהוות את חלק האקוויוויאליציה של הסימולציה.
ניתן לשנות את משך הפעלת האקוויוויאליות על-ידי עריכת שורה 19 של הקובץ. שלושת הננו-שניות הנותרות מהוות את חלק הייצור, ולכן משמשות לתפוקת אורך ממד z הממוצע. בנוסף, קובץ ה- Script מפיק פלט של קובץ TXT אחר המספק ערכים של lz כפונקציה של timestep וכן קובץ PNG, המתווה את אותם נתונים.
ניתן להשתמש בעלילה כדי לאמת את שיווי המשקל של סימולציית NPT. כדי לשחזר את תא העל באמצעות הגובה הממוצע שנקבע ב- NPT, העתק את קובץ הנתונים שנוצר קודם לכן לספריה עובדת חדשה ושנה את שמו כ- data.myadsorbate. לאחר מכן ערכו את קובץ הנתונים החדש כדי לשנות zlo ל- 0.0 ו- zhi לערך lz מפלט הערך הממוצע בקובץ TXT.
העתק את קובץ הקלט של LAMMPS לספריה העובדת החדשה. ערכו את משתנה הקבוצה בשורה 32 כדי לציין את אינדקסי סוגי האטום עבור אטומי מימן חמצן ומים ואת משתנה הקבוצה בשורה 33 כדי לציין את אינדקסי סוג האטום עבור אטומי הפלטינה והספוג. לאחר מכן, בקו 16, ערוך את משתנה runStep כך שהוא ארוך מספיק כדי לכלול הפעלת שיווי משקל והפעלה של ייצור.
הקלד את הפקודה לביצוע LAMMPS בממשק שורת הפקודה כדי להפעיל את תוכנת LAMMPS. פעולה זו תפעיל הדמיית NVT קבועה על מולקולות המים, וקובץ פלט המפתח ייווצר. סימולציית NVT כוללת חלק שיווי משקל וחלק ייצור.
חלק הייצור מתחיל כאשר האנרגיה של המערכת התווה נגד רמות הזמן כבוי. כדי להתחיל, לפתוח את התסריט לכל החיים קשר מימן. שנה את משתנה ההתחלה בפועל בשורה 22 כדי להגדיר את פרק הזמן של מסגרת הזמן הראשונה.
שנה את משתנה timestep בשורה 23 כדי להגדיר באיזו תדירות מסגרות נכתבות לקובץ המסלול של LAMMPS. שנה N_first ו- N-last בשורה 24 ו- 25 כדי להגדיר את הteps הראשון ו האחרון, ושנה את המשתנה לעולם לא בשורה 26 כדי לקבוע אם מסגרות עוקבות ייחשבו או ידלגו עליהן. הגדר את מספר השורות לכל מקטע מסגרת של קובץ המסלול על-ידי שינוי משתנה קו המסגרת בשורה 27.
בנוסף, ערוך שורות 31 עד 35, כדי לציין אילו סוגי אטום בתוך הנתונים. קובץ myadsorbate שייך לסופוירבאט, וסוגי האטומים שייכים למולקולות המים. התסריט מנתח את תצורות המים בהפעלת הייצור וקובע אם מולקולות מים כלשהן נקשרות למימן לסוחת.
לאחר מכן הוא סופר את זמן הסימולציה שכל קשר מימן נשאר שלם ומדווח על מידע זה כהפצה של אורך חיים של קשרי מימן ביחידות של פיקו-שניות. LAMMPS כותב את התצורה של מולקולות מים לקובץ כל 1, 000 femtoseconds, המהווה את ברירת המחדל בקובץ הקלט LAMMPS שסופק. הוא מזהה ומשליך את שתי תצורות nanoseconds הראשון בשווי של הקובץ כפי שהם מהווים את החלק שיווי משקל של הסימולציה ומשתמש בשלושת ננו שניות הנותרים כדי לחשב אורך חיים של קשר מימן.
כדי להפעיל את קובץ ה- Script, הקלד בממשק שורת הפקודה. לאחר מכן הוא יוצר קובץ DAT. התווה את הנתונים בקובץ כדי להציג את ההתפלגות של משכי חיים של קשרי מימן שהתרחשו במהלך הדמיית NVT.
כדי לקבוע את תוספת הזמן לשימוש עבור מרווח הזמן לדגימת זמן, השתמש בהפרש קבוע של זמן הגדול או שווה לאורך החיים המרבי של קשר מימן. קבע את מספר התצורות מהפעלת הייצור של מסלול NVT FFMD כדי לדגום כך שהזמן המינימלי בין תצורות יהיה שווה או גדול ממרווח דגימת הזמן שזוהה קודם לכן. בקובץ ה- Script של חילוץ המסגרות שנכתב קודם לכן, ערוך את ערך ברירת המחדל עבור משתנה מספר המסגרות בשורה 21 כדי לציין את מספר התצורות שיש לחלץ.
כדי להפעיל את קובץ ה- Script, הקלד את שם קובץ ה- Script בממשק שורת הפקודה. פעולה זו תפרט רשימה של זמני הדמיה המתאימים לתצורות שיש לחלץ מקובץ הסימולציה של NVT. תצורות אלה יכולות לשמש כמבנים התחלתיים בסימולציות AIMD או QM.
בהליך זה, FFMD שימש ליצירת התצורה הראשונית של מולקולות מים. סימולציית AIMD מראה כי מולקולת מים כי הוא במקור מימן מלוכדים לספוג אלכוהול סוכר על משטח פלטינה-111 מפשט את המימן מן adsorbate אלכוהול ומפקיד מימן שני על פני הפלטינה-111. המבנים של מולקולות מים נוזליים תלויים בהגדרות קלט.
לקביעתם בצורה לא נכונה יכולות להיות השפעות לא מכוונות על מבני המים. בנתון זה, הצד השמאלי הוא המבנה ההתחלתי להפעלת FFMD. והיד ימין נמצא בתוך פיקו-שניה אחת מהתחלת הסימולציה.
סימולציית FFMD מתפוצצת עקב הגדרות כוח גדולות מבחינה פיזית, מה שגורם למולקולות מים להתרחק מפני השטח. ניתן להשתמש בתצורות במכניקת הקוונטים, או בסימולציית QMM, או להשתמש בהן כדי לנתח סטטיסטיקות הקשורות למיקומים המרחביים של המולקולות. טכניקה זו סוללת את הדרך לחוקרים לחקור את התפקידים שיש לסביבות תגובה נוזלית על קטליזה על ידי יצירת תצורות בפועל של מולקולות נוזליות בממשקים קטליטיים.
