detci - אונליין בענן

תָכְנִית:

שֵׁם

detci - תוכנית אינטראקציה של תצורת דטרמיננט

תיאור

התכנית דטסי מבצע אלכסון של האופרטור ההמילטוני האלקטרוני הלא-רלטיביסטי ב-
על בסיס דטרמיננטים של סלייטר. ניתן לבחור את קבוצת הדטרמיננטים שבהם נעשה שימוש (מרחב CI) ב
מגוון דרכים. התוכנית יכולה להתמודד עם כל מרחב CI שניתן לנסח כ-
CI מרחב פעיל מוגבל. זה כולל CISD, CISDT, CISDTQ וכו', עד CI מלא, כ
כמו גם CI רב-הפניות שבהן ההפניות נבחרות ככל הגורמים הדטרמיננטיים ב
אשר עד n אלקטרונים מעוררים במרחב פעיל של MO. זה כולל CISD[T],
CISD[TQ] ו-CI מסדר שני (SOCI).

ביבליוגרפיה

CI של מרחב פעיל מוגבל:

1. אלגוריתמים של אינטראקציה תצורתית מבוססת דטרמיננטים עבור ניתוחים מלאים ומוגבלים
מרחבי אינטראקציה של תצורה, ג'יי אולסן, ב.ו. רוס, פ. יורגנסן, וה.י. א.
ג'נסן, J. Chem. Phys. 89, 2185 (1988).

2. מעבר למגבלה של מיליארד בחישובי אינטראקציה מלאה (FCI),
ג'. אולסן, פ. יורגנסן וג'. סימונס, Chem. Phys. Lett. 169, 463 (1990).

תת-מרחבים וירטואליים שלישוניים (RAS IV):

1. פונקציות גל וריאציוניות קומפקטיות הכוללות פונקציות משולשות וארבע-פונקציות מוגבלות
עירורים, CD Sherrill ו-HF Schaefer, J. Phys. Chem. 100, 6069-6075
(1996).

תוכנית DETCI:

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

קבצים נדרש

input.dat - קובץ קלט
קובץ 71 - אינטגרלים של אלקטרון אחד שעברו טרנספורמציה
קובץ 72 - אינטגרלים של שני אלקטרונים שעברו טרנספורמציה

זמני קבצים מְשׁוּמָשׁ

קובץ 50 - אלכסון של המילטוניאן
קובץ 51 - וקטורי CI
קובץ 52 - וקטורי סיגמא
קובץ53 - קובץ D (וקטורי תיקון)

קבצים מְעוּדכָּן

output.dat - קובץ פלט

קלט פורמט

הארגומנטים הבאים של שורת הפקודה זמינים:

-שֶׁקֶט זה נותן את אותה התוצאה כמו הדפס=0.

-o fname
נותן את שם הקובץ עבור קובץ הפלט. ברירת המחדל היא output.dat.

-e אפשרות זו גורמת לכתיבת סך אנרגיית ה-CI או האנרגיות לקובץ בשם
detci_energies.dat.

-c ערך
נותן התכנסות רופפת יותר על וקטור ה-CI, שימושי בחישובי DETCAS.
ערך הוא מספר ממשי, לא מספר שלם כמו ב הִתכַּנְסוּתההתכנסות שבה נעשה שימוש
יהיה המפסיד יותר של ערך ו הִתכַּנְסוּת.

קלט נוסף עבור תוכנית זו נקרא מהקובץ. מילות המפתח הבאות הן
תָקֵף:

הִתכַּנְסוּת = מספר שלם
התכנסות רצויה על וקטור CI. התכנסות מושגת כאשר ה-RMS של ה-
השגיאה בוקטור CI קטנה מ-10**(-n). ברירת המחדל היא 4 לאנרגיות ו-7
עבור גרדיאנטים. זהו אינו אותו קריטריון התכנסות וקטור CI כפי שנמצא ב
גוגאצי.

DOCC = מערך_שלם
וקטור זה נותן את מספר האורביטלים המאוכלסים כפליים בכל אי-רפ. יש
אין ברירת מחדל.

SOCC = מערך_שלם
וקטור זה נותן את מספר האורביטלים המאוכלסים באופן יחיד בכל אי-רפ. יש
אין ברירת מחדל.

שיטת_דיאגנוסטיקה = מחרוזת
זה מציין באיזו שיטה יש להשתמש באלכסון ההמילטוניאן.
אפשרויות תקפות הן: RSP, כדי ליצור את כל מטריצת H ולבצע אלכסון באמצעות libciomr
כדי לקבל את כל הערכים העצמיים (nb דורש זיכרון עצום); אולסן, להשתמש במונח של אולסן
שיטת תת-מרחב הפוך מותנית (1990); מיטרושנקוב, כדי להשתמש ב-2x2
שיטת אולסן/דיווידסון; ו DAVIDSON (אוֹ SEM) להשתמש בהתרחבות הסימולטנית של ליו
שיטה, הזהה לשיטת דייווידסון אם נמצא רק שורש אחד.
קיים גם מצב ניפוי שגיאות SEM, מבחן סמסטר. ה SEM השיטה היא הכי
חזק, אך הוא דורש גם 2(N*M)+1 וקטורי CI בדיסק, כאשר N הוא המקסימום
מספר האיטרציות ו-M הוא מספר השורשים.

מרכך מקדים = מחרוזת
זה מציין את סוג המוקדם שיש להשתמש בו באלכסון שנבחר
שיטה. האפשרויות התקפות הן: DAVIDSON אשר מקרוב את מטריצת ההמילטונית
על ידי האלמנטים האלכסוניים; H0BLOCK_INV אשר משתמש בהמילטוניאן מדויק של
H0_BLOCKSIZE והופך אותו במפורש; GEN_DAVIDSON אשר עושה ספקטרלי
פירוק של H0BLOCK; ITER_INV באמצעות גישה איטרטיבית כדי להשיג את
וקטור תיקון של H0BLOCK. ה-H0BLOCK_INV, GEN_DAVIDSON ו-ITER_INV
כל הגישות שוות ערך פורמלית, אך ה-ITER_INV פחות חישובי
יקר. ברירת המחדל היא DAVIDSON.

הפניה = מחרוזת
זה מציין את סוג פונקציית הייחוס. זהו RHF או ROHF. UHF ו-
TWOCON אינם נתמכים. עבור ROHF, ריבוי של 1 מרמז על מעטפת פתוחה.
סינגלט. התוכנית תרוץ עבור סינגלטים פתוחים, אך היא טרם
מותאם כראוי לשימוש בהתייחסות נכונה של שני דטרמיננטים במקרה זה, כך
הפעלה עם הפניות סינגלט פתוחות אינה מומלצת למעט עבור CI מלאים.

עדכון = מחרוזת
DAVIDSON משתמש בנוסחת וקטור העדכון או התיקון הסטנדרטית של DAVIDSON, בעוד
אולסן משתמש בווקטור תיקון OLSEN. ברירת המחדל היא DAVIDSON.

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

HD_AVE = מחרוזת
HD_EXACT משתמש באנרגיות האלכסון המדויקות שמובילות לוקטורי התפשטות אשר
לשבור סימטריית ספין. HD_KAVE ממוצע האנרגיות האלכסוניות מעל צימוד ספין
קבוצה המניבה וקטורי התפשטות טהורים של ספין. ORB_ENER משתמש בסכום של מסלולים
קירוב אנרגיה נותן ספין וקטורי התפשטות טהורים אך בדרך כלל מכפיל את
מספר איטרציות של דיווידסון. אוונגליסטים משתמש בסכומים ובהפרשים של מסלולים
אנרגיות עם אנרגיית הייחוס SCF כדי לייצר וקטורי התפשטות טהורים של ספין.
ליינינגר קירוב אשר מחסיר את תרומת האלקטרון האחד מה-
אנרגיות מסלוליות, מכפילות ב-0.5, ומוסיפות בחזרה את תרומת האלקטרון האחד
בתוך, לייצר וקטורי התפשטות טהורים של ספין ופותחו על ידי שלי ועובדים כ
וכן אוונגליסטים.

NODFILE = בוליאני
אפשרי רק אם NUM_ROOTS = 1. משתמש במרחב הווקטור האחרון בקובץ BVEC כדי
לכתוב DVEC מאפס במקום להשתמש בקובץ DVEC נפרד.

התכנסות אנרגיה = מספר שלם
התכנסות רצויה על אנרגיית ה-CI. ברירת המחדל היא 6 עבור אנרגיות נקודה בודדת
ו-8 עבור גרדיאנטים או CASSCF.

EX_LVL = מספר שלם
רמת עירור עבור עירורים לתוך אורביטלים וירטואליים (ברירת מחדל 2, כלומר CISD).

VAL_EX_LVL = מספר שלם
רמת עירור עבור הפניות באורביטלים של RAS II. ברירת המחדל היא אפס.

FCI = בוליאני
אם דגל זה מוגדר ל נכון, אז אחסון המחרוזות פשוט יותר עבור Full
CI והחישוב דורש פחות תקורה. עם זאת, התוצאות הסופיות צריכות להיות
להיות זהים לאלה כאשר FCI = שֶׁקֶרעלול לגרום לתוצאות בלתי צפויות אם FCI =
נכון אבל EX_LVL אינו עולה בקנה אחד עם CI מלא.

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

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

FREEZE_CORE = בוליאני
אפשרות זו קובעת האם יש לכלול את אורביטלי הליבה הקפואים
באופן מרומז (אמת) או במפורש (שקר). במקרה הראשון, האנרגטי
תרומות מהאורביטלים הקפואים של הליבה מקופלות לתוך האלקטרון החד-אלקטרוני
אינטגרלים ולתוך "אנרגיית הליבה הקפואה" המחושבת על ידי תוכנית הטרנספורמציה.
ברירת המחדל היא אמת.

ייצוא_וקטור = בוליאני
זה מציין האם לאחסן וקטורים מתכנסים בסוף הריצה.
וקטור/ים מאוחסנים בפורמט שקוף כך שתוכניות אחרות יוכלו להשתמש בו
בקלות. הפורמט מצוין ב src/lib/libqt/slaterdset.h. ברירת המחדל היא
שקר.

NUM_EXPORT = מספר שלם
If ייצוא_וקטור מוגדר כ-true, פעולה זו קובעת את מספר הווקטורים ש
יש לייצא בסוף הריצה. ברירת המחדל היא 1.

GUESS_VECTOR = מחרוזת
זה מציין איזה סוג של וקטור ניחוש יש להשתמש באיטרציה של CI. נכון לעכשיו
משמש רק בשיטת איטרציה של SEM. ערכים מקובלים הם יחידה עבור וקטור יחידה
ניחוש (NUM_ROOTS ו NUM_INIT_VECS שניהם חייבים להיות 1); H0_BLOCK להשתמש בווקטורים עצמיים
מתת-מטריצת בלוק H0 (ברירת מחדל); DFILE כדי להשתמש NUM_ROOTS התכנסו בעבר
וקטורים בקובץ D; ו MP2 להשתמש בפונקציית הגל של MP2 כניחוש (לא עובד)
כרגע).

גודל בלוק H0 = מספר שלם
פרמטר זה מציין את גודל בלוק ה-"H0" של ההמילטוניאן, שהוא
נפתר במדויק. נבחרו n דטרמיננטים בעלי אנרגיית SCF הנמוכה ביותר, ו-a
תת-מטריצת ההמילטוניאן נוצרת באמצעות דטרמיננטים אלה. תת-מטריצת זו היא
משמש להאצת התכנסות איטרציות ה-CI ב- אולסן ו מיטרושנקוב
סכמות איטרציה, וגם למצוא ניחוש התחלתי טוב עבור ה- SEM שיטה אם
GUESS_VECTOR = H0_BLOCKברירת המחדל היא 40. שים לב שהתוכנית עשויה לשנות את
גודל נתון עבור Ms=0 מקרים (Ms0 = נכון) אם הוא קובע שבלוק H0 כולל
רק איבר אחד מתוך זוג דטרמיננטים הקשורים על ידי סימטריית היפוך זמן. עבור
גדלי בלוקים קטנים מאוד, זה יכול להעלות על הדעת לבטל את כל בלוק H0;
התוכנה אמורה להדפיס אזהרות אם זה קורה.

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

MAX_DET = מספר שלם
קובע את המספר המרבי של דטרמיננטים; אם מרחב ה-CI גדול מזה, ה-
התוכנית בוטלת. אפשרות זו קיימת כדי להבטיח שחישובים גדולים מאוד לא יבוצעו
מופעל בטעות. במהלך שלב הפיתוח הנוכחי, ברירת המחדל היא 10000, אבל
זה יועלה בקרוב.

מקסיטר = מספר שלם
מספר איטרציות מקסימלי לאלכסון ההמילטוניאן. ברירת המחדל היא 12.

Ms0 = בוליאני
If נכון, השתמש ברכיב Ms=0 של המצב. ברירת המחדל היא נכון אם קליפה סגורה ו
ל שֶׁקֶר אחרת. קשור ל- S פָּרָמֶטֶר.

הדפס = מספר שלם
ערך זה מציין את מספר הגורמים שיודפסו, יחד עם
המקדמים שלהם, ברשימת הגורמים החשובים ביותר ב-CI הסופי
וקטור. ברירת המחדל היא 20.

NUM_ROOTS = מספר שלם
ערך זה נותן את מספר השורשים שיש לקבל מהחילוני
משוואות. ברירת המחדל היא אחת. אם נדרש יותר משורש אחד, הגדר שיטת_דיאגנוסטיקה
ל SEM (או, במקרים קטנים מאוד, RSP or מבחן סמסטר).

NUM_INIT_VECS = מספר שלם
מספר הווקטורים ההתחלתיים לשימוש בהליך האיטרטיבי של CI. ברירת המחדל היא
מספר השורשים.

OPDM = בוליאני
If נכון חשב את מטריצת צפיפות החלקיקים האחד ועשה OPDM_WRITE ברירת מחדל ל
נכון. ערך ברירת המחדל של OPDM is שֶׁקֶר.

OPDM_FILE = מספר שלם
קובץ (מספר יחידה) לכתיבת מטריצת צפיפות החלקיק הבודד אם OPDM_WRITE =
נכוןערך ברירת המחדל הוא כעת 73.

OPDM_WRITE = בוליאני
סמן אם לכתוב את מטריצת צפיפות החלקיקים האחד לדיסק או לא.

OPDM_PRINT = בוליאני
סמן אם להדפיס את מטריצת צפיפות החלקיקים האחד או לא.

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

WRTNOS = בוליאני
דגל מציין האם לכתוב את האורביטלים הטבעיים של CI ל-PSIF_CHKPT.

ORBSFILE = מספר שלם
קובץ (מספר יחידה) לכתיבת אורביטלים טבעיים שונים של CI. ערך ברירת המחדל הוא
76.

OPDM_AVE = בוליאני
סמן האם לממוצע את ה-OPDM על פני מספר שורשים על מנת לקבל
מטריצת צפיפות של חלקיק אחד בממוצע מצב. ניתן לאלכסון מטריצת צפיפות זו
כדי לקבל את האורביטלים הטבעיים של CI.

ORBS_ROOT = מספר שלם
דגל להגדרת מספר השורש שעבורו ייכתבו אורביטלים טבעיים של CI
PSIF_CHKPT. ערך ברירת המחדל הוא 1 (השורש הנמוך ביותר).

הדפס = מספר שלם
אפשרות זו קובעת את רמת הפירוט של הפלט. ערך של 1 או 2 מציין
הדפסה מינימלית, ערך של 3 מציין הדפסה מפורטת. ערכים של 4 או 5 הם
משמש לניפוי שגיאות. אין להשתמש ברמה 5 אלא אם כן מקרה הבדיקה קטן מאוד (למשל
אחסון H2O CISD).

שורש = מספר שלם
השורש שיש לרשום עבורו את מטריצת צפיפות שני החלקיקים (צפיפות החלקיק האחד
מטריצות נכתבות עבור כל השורשים). שימושי עבור CASSCF או CI ספציפיים למצב
אופטימיזציה במצב מעורר.

S = מספר שלם
הערך של מספר הקוונטים של הספין S ניתן על ידי אפשרות זו. ברירת המחדל היא 0
(סינגלט). הדבר היחיד שמשתמשים בו בפועל הוא קביעת הפאזה של
החצי העודף של וקטור CI כאשר משתמשים ברכיב Ms=0 (כלומר, Ms0 =
נכון). במקרים בהם S אינו מספר שלם, אין צורך להזין פרמטר זה
מכיוון שמצב כזה לא יכול לכלול רכיב Ms=0.

TPDM = בוליאני
If נכון חשב את מטריצת צפיפות שני החלקיקים ועשה TPDM_WRITE ברירת מחדל ל
נכון. ערך ברירת המחדל של TPDM is שֶׁקֶר.

קובץ_TPDM = מספר שלם
קובץ (מספר יחידה) לכתיבת מטריצת צפיפות שני החלקיקים אם TPDM_WRITE =
נכוןערך ברירת המחדל הוא כעת 74.

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

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

יש גם קלט פחות נפוץ, שמשתמשים מתחילים ב-PSI לא יצטרכו להשתמש בו.
שימוש.

בנדזולי = בוליאני
השתמשו בכמה שגרות לחישוב סיגמא בהתבסס על המאמרים של בנדזולי ואחרים.
נראה איטי יותר ולא משתלם; עלול להיעלם בסופו של דבר. עובד רק עבור
CI מלא ואני לא זוכר אם יכולתי לראות איך התוכנית החכמה שלהם יכולה להיות
הורחב ל-RAS באופן כללי.

CALC_SSQ = בוליאני
אם TRUE, חשב את ערך התוחלת של האופרטור S^2 עבור CI הסופי
פונקציית גל עבור כל שורש. באופן עקרוני, DETCI אמור להניב S^2 פונקציות עצמיות.
ברירת המחדל היא FALSE.

כווץ בגודל מספר שלם
נותן את מספר הווקטורים שיש לשמור כאשר תת-מרחב דייווידסון קורס (ראה
מקסנווקט אם גדול מאחד, תת-המרחב המכווץ שומר על הטוב ביותר
אומדן וקטור CI עבור n האיטרציות הקודמות. ברירת המחדל היא 1.

יחידת_TMP_FIRST = מספר שלם
נותן את מספר הקובץ (יחידה) המשויך לקובץ הגירסה הראשון בו השתמש DETCI.
קבצי גירוד אחרים ממוספרים ברצף מנקודה זו, בסדר
H(diag), C, S, D. כל אחד מהקבצים הלוגיים הללו תופס מספר קבצים פיזיים
שצוינו על ידי פרמטרי הקלט המעורפלים עוד יותר NUM_HD_TMP_UNITS,
NUM_C_TMP_UNITS, NUM_S_TMP_UNITS, NUM_D_TMP_UNITS. המשתמש יכול גם לציין
נקודות התחלה שונות עבור כל אחת מהקבוצות הללו באמצעות הפרמטרים
יחידת_TMP_HD_FIRST_ וכן הלאה. פיצול קובץ על פני מספר יחידות עשוי לסייע
הימנעו מבעיית גודל המספר השלם בטיפול בקבצים גדולים הקיימת ב
DETCI ובספריות PSI I/O; אבל שוב, לא בדקתי את זה כדי לראות מה
קורה. היחידה הראשונה של כל מקטע מודפסת תחת הכותרת קבצים ב
פלט הפרמטר שמתחיל את ריצת ה-DETCI.

פדרציית חוץ = בוליאני
קובע האם האורביטלים הקפואים של הליבה נחשבים כקפואים באמת (כלומר,
נעדר לחלוטין מהחישוב, פדרציית חוץ = נכון) או האם הם נוכחים אבל
מוגבל לתפוסה כפולה (פדרציית חוץ = שֶׁקֶרבתוכנית GUGA CI, זה
ההבחנה בין מה שהוא מכנה אורביטלים FZC ו-COR. באופן כללי, ה-
אינטגרלים עבור אורביטלים קפואים של ליבה אינם נחוצים על ידי DETCI, אך ייתכן שיהיה צורך בהם.
עבור MCSCF או גרדיאנטים.

אייקור = מספר שלם
מציין כיצד לטפל באחסון במאגר של וקטורי CI. ערך של 0 הופך את התוכנית ל...
לבצע קלט/פלט תת-בלוק RAS אחד בכל פעם; 1 משתמש בווקטורי CI שלמים בכל פעם; ו-2
משתמש בבלוק irrep אחד בכל פעם. ערכים של 0 או 2 גורמים לחוסר יעילות מסוים ב-
קלט/פלט (הדורש קריאות מרובות של וקטור C בעת בניית H באיטרטיבי)
תת-מרחב אם DIAG_METHOD = SEM), אך דורשים פחות זיכרון ליבה.

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

מקסנווקט = מספר שלם
נותן את המספר המקסימלי של וקטורי תת-מרחב דייווידסון שניתן להחזיק בדיסק למשך
מקדם ה-CI ווקטורי הסיגמא. (יש וקטור H(diag) אחד ומספר
(מספר הווקטורים D שווה למספר השורשים). כאשר מספר הווקטורים בדיסק
מגיע לערך של מקסנווקט, תת-המרחב של דייווידסון יקרוס ל
כווץ בגודל וקטורים עבור כל שורש. זה מאוד מועיל לחיסכון בשטח דיסק.
ברירת המחדל היא מקסיטר * NUM_ROOTS + NUM_INIT_VECS.

מעורב = בוליאני
זה קובע האם עירורים "מעורבים" של RAS II/RAS III מותרים ב-CI
שטח. זה שימושי להצבת אילוצים נוספים על RAS CI.

MIXED4 = בוליאני
זה דומה למילת המפתח MIXED, אך מתייחס לעירורים לתוך RAS IV.

נזירות = מספר שלם
מספר קבצי הגרסה הקודמת לשימוש באחסון וקטורי C (וגם עבור הסיגמא)
וקטורים).

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

קובץ OEI = מספר שלם
מילת מפתח זו מאפשרת למשתמש לציין את קובץ האינטגרל החד-אלקטרוני שעבר טרנספורמציה.
ברירת המחדל היא 71.

הדפס_CIBLKS = בוליאני
מציין האם התוכנית צריכה להדפיס סיכום של כל הבלוקים ב-
וקטור CI (אשר ניתן ליצוק לצורת מטריצה, ראה הפניות).

R4S = בוליאני
מגביל את מחרוזות RAS IV לקבוצה המינימלית, וחוסך זיכרון. אם אתה
מודאג לגבי אפשרות זו, עליך לפנות לדוד לקבלת ייעוץ אלא אם כן אתה
מומחה DETCI.

REF_SYM = מספר שלם
אפשרות זו מאפשרת למשתמש לחפש וקטורי CI בעלי אי-רפ שונה מזה של ה-
הפניה. זה כנראה הגיוני רק עבור CI מלא, וסביר להניח שזה לא יהיה
עבודה עם ניחושים של וקטור יחידה. המספור מתחיל מאפס עבור ה-total-
אי-רפציה סימטרית

REPL_OTF = בוליאני
אומר ל-DETCI האם לבצע החלפות מחרוזות תוך כדי תנועה או לא. יכול לחסוך
כמות עצומה של זיכרון (במיוחד עבור CI קטועים) אבל קצת רעוע ו
לא נבדק כבר זמן מה. כפי שאני זוכר, זה עובד רק עבור סוגים מסוימים של
חישובי RAS. צרו קשר עם דיוויד לקבלת סיוע. בסופו של דבר, תוך כדי תנועה
יש לחדש את הדברים החלופיים בצורה הרבה יותר חכמה כדי שלא ייקח
עידנים של זמן מעבד. עבודה בכיוון זה החלה פעם אחת וייתכן שתסתיים
בסופו של דבר.

איתחול = בוליאני
אפשרות זו מאפשרת למשתמש לחדש איטרציה של DETCI שהסתיימה
בטרם עת. זה מניח שוקטורי ה-CI וה-sigma נמצאים בדיסק; מספר ה-
וקטורים שצוינו על ידי RESTART_VECS מכווץ לווקטור אחד לכל שורש.

RESTART_VECS = מספר שלם
If איתחול = נכון זה מציין את מספר וקטורי ה-CI (והסיגמא) שיש לקרוא מהם
דיסק. בדרך כלל זהו מספר האיטרציות שהושלמו בהצלחה מ-
הריצה הקודמת כפול מספר השורשים עבור אותה ריצה.

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

TEI_FILE = מספר שלם
מילת מפתח זו מאפשרת למשתמש לציין את קובץ האינטגרל בעל שני האלקטרונים שעבר טרנספורמציה.
ברירת המחדל היא 72.

MPN = בוליאני
כאשר אפשרות זו היא TRUE, DETCI יחשב את סדרת MPn לסדר k כאשר k
נקבע על ידי maxnvect. עבור מערכות פתוחות (REF=ROHF, WFN = ZAPTN), DETCI
יחשב את סדרת ZAPTn. GUESS_VECTOR חייב להיות מוגדר כ-TRUE. HD_OTF חייב להיות
יש להגדיר ל-TRUE. יש להגדיר את HD_AVE ל-orb_ener.

שמור_MPN2 = מספר שלם
כאשר MPN הוא TRUE ו-WIGNER הוא TRUE, אפשרות זו הופכת לתקפה. אם מוגדרת ל-1
אזי אנרגיה MP(2n-1) נשמרת. אם מוגדר ל-2 אזי אנרגיה MP(2n-2) נשמרת. אם בכלל
ערך אחר MPn אנרגיה נשמרת. ברירת המחדל היא 0.

9 פברואר, 1996 דטסי(1)

השתמש ב-detci באינטרנט באמצעות שירותי onworks.net