זוהי הפקודה detcas שניתן להריץ בספק האירוח החינמי של OnWorks באמצעות אחת מתחנות העבודה המקוונות המרובות שלנו, כגון Ubuntu Online, Fedora Online, אמולטור מקוון של Windows או אמולטור מקוון של MAC OS
תָכְנִית:
שֵׁם
detcas - אינטראקציה של תצורת קביעה / פעילה מלאה
תוכנית שדה עקבי עצמי בחלל
תיאור
התכנית detcas מבצע שדה שלם-פעיל-חלל עקבי-עצמי (CASSCF)
אופטימיזציה של אורביטלים מולקולריים באמצעות הליך דו-שלבי בשילוב עם
תוכנית האינטראקציה של תצורת הקובע DETCI. התוכנית פשוטה למדי ו
כיום משתמש בגישת ניוטון-ראפסון לעדכון האורביטלים, תוך שימוש בגישה פשוטה
הססיאן מסלולי משוער. ההתכנסות מואצת באמצעות היפוך ישיר של Pulay של
הליך המשנה האיטרטיבי (DIIS). הקוד נכתב כדי לאפשר יותר כללי
פונקציות גל שאינן כוללות בהכרח טיפול CI מלא של המרחב הפעיל. ב
במיוחד, כל פונקציית גל CI של מרחב פעיל מוגבל (RAS) הנתמכת על ידי DETCI
בשימוש, המאפשר פונקציות גל MCSCF מסוג RASSCF.
ביבליוגרפיה
הססיאן אורביטלי משוער:
1. G. Chaban, MW Schmidt, ו-MS Gordon, Theor. צ'ים. Acta 97, 88-95 (1997).
CI שטח פעיל מוגבל:
1. אלגוריתמי אינטראקציה של תצורה מבוססי דטרמיננטים עבור שלם ומוגבל
Configuration Interaction Spaces, J. Olsen, BO Roos, P. Jorgensen, and HJ Aa.
Jensen, J. Chem. פיזי. 89, 2185 (1988).
שטח פעיל מוגבל SCF:
1. P.-A. Malmqvist, A. Rendell, ו-BO Roos, J. Phys. Chem. 94, 5477 (1990).
תוכנית DETCI:
1. תקליטור שריל, אלגוריתמים חישוביים עבור רפרנסים מלאים ומרובים בקנה מידה גדול
תצורה אינטראקציה Wavefunctions, עבודת דוקטורט, אוניברסיטת ג'ורג'יה, אתונה,
GA, 1996.
2. CD Sherrill ו-HF Schaefer, The Configuration Interaction Method: Advances
בגישות מתואמות גבוהות, עו"ד. Quantum Chem. 34, 143-269 (1999).
קבצים נדרש
input.dat - קובץ קלט
file78 - אינטגרלים של אלקטרון אחד של MO (אופרטור fzc)
file72 - אינטגרלים של שני אלקטרונים MO
file73 - מטריצת צפיפות חלקיק אחת MO
file74 - מטריצת צפיפות שני חלקיקים MO
file75 - MO Lagrangian
קבצים מְעוּדכָּן
output.dat - קובץ פלט
file14.dat - תיעוד של אנרגיות ושיפועים מסלוליים
קלט פורמט
הארגומנטים הבאים של שורת הפקודה זמינים:
-שֶׁקֶט זה נותן את אותה תוצאה כמו PRINT=0.
-o fname
נותן את שם הקובץ עבור קובץ הפלט. ברירת המחדל היא output.dat.
קלט נוסף עבור תוכנית זו נקרא מהקובץ input.dat. הנפוץ יותר
מילות המפתח הן:
WFN = מחרוזת
סוגי פונקציית הגל הנתמכים הם CASSCF ו- RASSCF.
הִתכַּנְסוּת = מספר שלם
התכנסות רצויה על שיפוע המסלול. התכנסות מושגת כאשר RMS
של השגיאה בשיפוע המסלול קטן מ-10**(-n). ברירת המחדל היא 4 עבור
חישובי אנרגיה ו-7 עבור שיפועים.
DOCC = מערך_שלמות
וקטור זה נותן את מספר האורביטלים התפוסים כפול בכל irrep. יש
אין ברירת מחדל.
SOCC = מערך_שלמות
וקטור זה נותן את מספר האורביטלים תפוסים בודדים בכל irrep. יש
אין ברירת מחדל.
ENERGY_CONVERGENCE = מספר שלם
התכנסות רצויה על סך האנרגיה של MCSCF. ברירת המחדל היא 7.
FROZEN_DOCC = מערך_שלמות
מספר האורביטלים הנמוכים ביותר באנרגיה הנמוכה ביותר בכל אחד מהם בלתי ניתן לצמצום
ייצוג שממש יוקפא (לא מעודכן ב-MCSCF). ה
נעשה שימוש בסידור כותנה של הייצוגים הבלתי ניתנים לתיקון. ברירת המחדל היא ה
וקטור אפס.
FROZEN_UOCC = וקטור שלם
מספר האורביטלים הבלתי תפוסים באנרגיה הגבוהה ביותר בכל ייצוג בלתי ניתן לצמצום
שממש יוקפא (לא מעודכן ב-MCSCF). ברירת המחדל היא האפס
וֶקטוֹר.
RESTRICTED_DOCC = מערך_שלמות
מספר האורביטלים הנמוכים ביותר באנרגיה הנמוכה ביותר בכל אחד מהם בלתי ניתן לצמצום
ייצוג שיעבור אופטימיזציה אך יישמר כפול ב-MCSCF.
האורביטלים האלה מגיעים אחרי ה FROZEN_DOCC אורביטלים. ברירת המחדל היא האפס
וֶקטוֹר.
RESTRICTED_UOCC = מערך_שלמות
מספר האורביטלים הבלתי תפוסים באנרגיה הגבוהה ביותר בכל ייצוג בלתי ניתן לצמצום
אשר יעבור אופטימיזציה אך יישמר ללא תפוסה ב-MCSCF. האורביטלים האלה מגיעים
לפני FROZEN_UOCC אורביטלים. ברירת המחדל היא וקטור האפס.
NCASITER = מספר שלם
מספר מרבי של איטרציות כדי לייעל את האורביטלים. אפשרות זו צריך be
מפורט in מה היא תקלות סעיף of קלט, כי זה צריך להיות גלוי ל-
תוכנית בקרה PSI. ברירת המחדל היא 1.
הדפס = מספר שלם
אפשרות זו קובעת את מילוליות הפלט. ערך של 1 או 2 מציין
הדפסה מינימלית, ערך של 3 מציין הדפסה מילולית. ערכים של 4 או 5 הם
משמש לניפוי באגים. אין להשתמש ברמה 5 אלא אם כן מקרה הבדיקה קטן מאוד (למשל
STO H2O CISD).
מילות המפתח הנפוצות פחות הן:
DIIS_FREQ = מספר שלם
הליך האקסטרפולציה של DIIS יתבצע כל ניסיון n איטרציות. ברירת המחדל
הוא 1.
DIIS_MAX_VECS = מספר שלם
המספר המרבי של וקטורי תת-מרחב עבור הליך DIIS. אחרי המספר הזה
של וקטורים מגיעים, וקטורים ישנים יותר יושמטו מתת-המרחב כ
נחוץ. ברירת המחדל היא 8.
DIIS_MIN_VECS = מספר שלם
המספר המינימלי של וקטורי תת-מרחב לפני אינטרפולציית DIIS יכול להיות
מְבוּצָע. ברירת המחדל היא 2. לא הגיוני שיהיו ערכים פחות מ-2.
DIIS_START = מספר שלם
מספר האיטרציה בעת היפוך ישיר של Pulay של תת המרחב האיטרטיבי (DIIS)
יש להפעיל נוהל להאצת התכנסות. ברירת המחדל היא 3.
לפני איטרציה זו, וקטורים אינם מתווספים לתת-המרחב DIIS. הראשון
צעד DIIS לא יינקט עד DIIS_MIN_VECS וקטורים נמצאים בתת-מרחב DIIS.
LAG_FILE = מספר שלם
קובץ (מספר יחידה) לקריאת המטריצה הלגראנית. ערך ברירת המחדל הוא
כרגע 75.
OPDM_FILE = מספר שלם
קובץ (מספר יחידה) לקריאת מטריצת צפיפות החלקיקים האחד. ערך ברירת המחדל
כיום בן 73.
SCALE_GRAD = בוליאני
אומר אם להתאים את שיפוע המסלול לפי המסלול המשוער (אלכסוני).
היסיאן. ברירת המחדל היא TRUE.
SCALE_STEP = ממשי
גורם קנה מידה עבור שלב סיבוב המסלול. ברירת המחדל היא 1.0.
TPDM_FILE = מספר שלם
קובץ (מספר יחידה) לקריאת מטריצת צפיפות שני החלקיקים. ערך ברירת המחדל
כיום בן 74.
LEVEL_SHIFT = בוליאני
אומר אם לאפשר שינוי רמה של המטריצה ההסיאנית. זה יכול לשמש
כדי להבטיח שההסיאן יהיה חיובי מוגדר עבור האיטרציות ההתחלתיות. אם
ההסיאן אינו חיובי מובהק ליד התכנסות, אז זה עשוי להיות סימן של
אי יציבות מספרית ב-MCSCF. ברירת המחדל היא TRUE.
SHIFT = ממשי
אומר עד כמה יש להזיז את האלמנטים האלכסוניים של ההסיאן לרמה
הסטה. ברירת המחדל היא 0.01.
DETERM_MIN = ממשי
הערך המינימלי המותר עבור הקובע של המטריצה ההסיאנית אם הרמה משתנה
פועל. זה לא בשימוש אם LEVEL_SHIFT=FALSE. ברירת המחדל היא 0.00001.
STEP_MAX = ממשי
זהו הסיבוב המרבי המותר למסלול יחיד. ברירת המחדל היא 0.30.
FORCE_STEP = בוליאני
זה מאפשר למשתמש לעקוף את הצעד המחושב ולכפות צעד ב-a
כיוון מסוים. זה יכול להיות שימושי אם מנסים להרחיק את החישוב בכוח
מנקודת אוכף. ברירת המחדל היא FALSE.
FORCE_PAIR = מספר שלם
זהו המדד עבור הזוג העצמאי שיש לסובב אם
FORCE_STEP=TRUE. זה מתעלם אחרת. אין ברירת מחדל.
FORCE_VALUE = ממשי
זהו ערך סיבוב האורביטלים לסיבוב האורביטלים שצוין על ידי
FORCE_PAIR. זה משמש רק אם FORCE_STEP מוגדר ל-TRUE. ערך ברירת המחדל הוא
0.0.
CHECK_HESSIAN = בוליאני
חשב את הערכים העצמיים של ההסיאן האורביטלי. זה טוב לבדיקה
פתרונות נקודת אוכף. יש להשתמש ללא שינוי מפלס. ברירת המחדל היא
FALSE.
EIGEN_VECTORS = בוליאני
האם אתה רוצה לחשב גם את הווקטורים העצמיים של ההסיאן המסלולי. זה
התעלמו אם CHECK_HESSIAN=FALSE. ברירת המחדל היא FALSE.
8 במאי, 1998 detcas(1)
השתמש ב-detcas באינטרנט באמצעות שירותי onworks.net
