detcas - ออนไลน์ในคลาวด์

โครงการ:

ชื่อ

detcas - การโต้ตอบการกำหนดค่าตัวกำหนด / เปิดใช้งานให้เสร็จสมบูรณ์
โปรแกรมภาคสนามที่สม่ำเสมอในอวกาศ

DESCRIPTION

โปรแกรม ดีแคส ดำเนินการฟิลด์ที่สอดคล้องกันในตัวเองของพื้นที่แอ็คทีฟสมบูรณ์ (CASSCF)
การเพิ่มประสิทธิภาพของออร์บิทัลโมเลกุลผ่านขั้นตอนสองขั้นตอนร่วมกับ
โปรแกรมโต้ตอบการกำหนดค่าดีเทอร์มิแนนต์ DETCI โปรแกรมค่อนข้างง่ายและ
ปัจจุบันใช้แนวทางของ Newton-Raphson เพื่ออัปเดต orbitals โดยใช้ simple
ประมาณวงโคจร Hessian การบรรจบกันถูกเร่งโดยใช้การผกผันของ .โดยตรงของ Pulay
ขั้นตอนย่อยซ้ำ (DIIS) รหัสถูกเขียนขึ้นเพื่อให้กว้างขึ้น
ฟังก์ชันคลื่นที่ไม่จำเป็นต้องมีการรักษา CI เต็มรูปแบบของพื้นที่แอ็คทีฟ ใน
โดยเฉพาะอย่างยิ่ง ฟังก์ชันคลื่น CI แบบจำกัดพื้นที่ใช้งาน (RAS) ที่รองรับโดย DETCI สามารถทำได้
ใช้แล้ว อนุญาตให้ใช้ฟังก์ชันคลื่น MCSCF ของประเภท RASSCF

ข้อมูลอ้างอิง

Hessian วงโคจรโดยประมาณ:

1. G. Chaban, MW Schmidt และ MS Gordon, Theor ชิม. แอคตา 97, 88-95 (1997).

พื้นที่ใช้งานที่ถูกจำกัด CI:

1. อัลกอริธึมการโต้ตอบการกำหนดค่าตามตัวกำหนดเพื่อความสมบูรณ์และถูกจำกัด
Configuration Interaction Spaces, J. Olsen, BO Roos, P. Jorgensen และ HJ Aa
เจนเซ่น, เจ. เคม. สรีรวิทยา 89, 2185 (1988).

SCF พื้นที่ใช้งานที่ถูกจำกัด:

1. ป.-อ. Malmqvist, A. Rendell และ BO Roos, J. Phys. เคมี. 94, 5477 (1990).

โปรแกรม DETCI:

1. CD Sherrill อัลกอริธึมการคำนวณสำหรับการอ้างอิงแบบเต็มขนาดใหญ่และแบบหลายตัว
Configuration Interaction Wavefunctions, วิทยานิพนธ์ระดับปริญญาเอก, มหาวิทยาลัยจอร์เจีย, เอเธนส์,
จอร์เจีย, 1996.

2. CD Sherrill และ HF Schaefer วิธีการโต้ตอบการกำหนดค่า: ความก้าวหน้า
ในแนวทางที่มีความสัมพันธ์สูง, Adv. ควอนตัมเคม. 34, 143-269 (1999).

ไฟล์ ที่จำเป็น

input.dat - ไฟล์อินพุต
file78 - MO อินทิกรัลอิเล็กตรอนหนึ่งตัว MO (ตัวดำเนินการ fzc)
file72 - MO อินทิกรัลสองอิเล็กตรอน
file73 - MO เมทริกซ์ความหนาแน่นหนึ่งอนุภาค
file74 - เมทริกซ์ความหนาแน่นสองอนุภาค MO
ไฟล์ 75 - MO ลากรองจ์

ไฟล์ ให้กับคุณ

output.dat - ไฟล์เอาต์พุต
file14.dat - บันทึกพลังงานและการไล่ระดับของวงโคจร

INPUT FORMAT

มีอาร์กิวเมนต์บรรทัดคำสั่งต่อไปนี้:

-เงียบ สิ่งนี้ให้ผลลัพธ์เช่นเดียวกับ พิมพ์=0.

-o ชื่อ
ให้ชื่อไฟล์สำหรับไฟล์ที่ส่งออก ค่าเริ่มต้นเป็น output.dat

อินพุตเพิ่มเติมสำหรับโปรแกรมนี้สามารถอ่านได้จากไฟล์ input.dat ที่ใช้กันมากขึ้น
คำหลักคือ:

ดับเบิลยูเอฟเอ็น = เชือก
ประเภทฟังก์ชันคลื่นที่รองรับคือ CASSCF และ RASSCF

คอนเวอร์เจนซ์ = จำนวนเต็ม
คอนเวอร์เจนซ์ที่ต้องการบนเกรเดียนต์ของวงโคจร การบรรจบกันจะเกิดขึ้นเมื่อ RMS
ของข้อผิดพลาดในการไล่ระดับวงโคจรน้อยกว่า 10**(-n) ค่าเริ่มต้นคือ 4 สำหรับ
การคำนวณพลังงานและ 7 สำหรับการไล่ระดับสี

อคส = จำนวนเต็ม_array
เวกเตอร์นี้ให้จำนวนของออร์บิทัลที่ถูกครอบครองสองเท่าในแต่ละ irrep มี
ไม่มีค่าเริ่มต้น

สกสค = จำนวนเต็ม_array
เวกเตอร์นี้ให้จำนวนของออร์บิทัลที่ถูกครอบครองโดยลำพังในแต่ละ irrep มี
ไม่มีค่าเริ่มต้น

พลังงาน_คอนเวอร์เจนซ์ = จำนวนเต็ม
การบรรจบกันที่ต้องการในพลังงาน MCSCF ทั้งหมด ค่าเริ่มต้นคือ 7

แช่แข็ง_DOCC = จำนวนเต็ม_array
จำนวนออร์บิทัลที่มีพลังงานต่ำสุดเป็นสองเท่าในแต่ละออร์บิทัลที่ลดไม่ได้
การเป็นตัวแทนซึ่งจะถูกระงับอย่างแท้จริง (ไม่ได้อัปเดตใน MCSCF) NS
ใช้การเรียงลำดับฝ้ายของการแสดงที่ปฏิเสธไม่ได้ ค่าเริ่มต้นคือ
เวกเตอร์ศูนย์

แช่แข็ง_UOCC = จำนวนเต็ม_เวกเตอร์
จำนวนออร์บิทัลว่างพลังงานสูงสุดในแต่ละการแสดงที่ลดไม่ได้
ซึ่งจะถูกระงับอย่างแท้จริง (ไม่ได้อัปเดตใน MCSCF) ค่าเริ่มต้นคือศูนย์
เวกเตอร์.

จำกัด_DOCC = จำนวนเต็ม_array
จำนวนออร์บิทัลที่มีพลังงานต่ำสุดเป็นสองเท่าในแต่ละออร์บิทัลที่ลดไม่ได้
การเป็นตัวแทนซึ่งจะได้รับการปรับให้เหมาะสม แต่ยังคงครอบครอง MCSCF สองเท่า
ออร์บิทัลเหล่านี้มาหลัง แช่แข็ง_DOCC ออร์บิทัล ค่าเริ่มต้นคือศูนย์
เวกเตอร์.

จำกัด_UOCC = จำนวนเต็ม_array
จำนวนออร์บิทัลว่างพลังงานสูงสุดในแต่ละการแสดงที่ลดไม่ได้
ซึ่งจะได้รับการปรับให้เหมาะสมแต่ยังคงว่างอยู่ใน MCSCF ออร์บิทัลเหล่านี้มา
ก่อนที่จะ แช่แข็ง_UOCC ออร์บิทัล ค่าเริ่มต้นคือเวกเตอร์ศูนย์

เอ็นซีซิเตอร์ = จำนวนเต็ม
จำนวนการวนซ้ำสูงสุดเพื่อปรับออร์บิทัลให้เหมาะสม ตัวเลือกนี้ น่า be
ที่ระบุไว้ in เริ่มต้น ส่วน of อินพุตเพราะมันจะต้องมองเห็นได้กับ
โปรแกรมควบคุม PSI ค่าเริ่มต้นเป็น 1

พิมพ์ = จำนวนเต็ม
ตัวเลือกนี้กำหนดความละเอียดของเอาต์พุต ค่า 1 หรือ 2 ระบุ
การพิมพ์ขั้นต่ำ ค่า 3 ระบุการพิมพ์แบบละเอียด ค่าของ 4 หรือ 5 are
ใช้สำหรับแก้จุดบกพร่อง อย่าใช้ระดับ 5 เว้นแต่กรณีทดสอบมีขนาดเล็กมาก (เช่น
STO H2O CISD)

คำหลักที่ใช้น้อยกว่าคือ:

DIIS_FREQ = จำนวนเต็ม
ขั้นตอนการอนุมาน DIIS จะพยายามทุก n การทำซ้ำ ค่าเริ่มต้น
คือ 1

DIIS_MAX_VECS = จำนวนเต็ม
จำนวนเวกเตอร์ซับสเปซสูงสุดสำหรับโพรซีเดอร์ DIIS หลังเลขนี้
ของเวกเตอร์ถึงแล้ว เวกเตอร์ที่เก่ากว่าจะถูกทิ้งจากสเปซย่อยเป็น
จำเป็น. ค่าเริ่มต้นคือ 8

DIIS_MIN_VECS = จำนวนเต็ม
จำนวนต่ำสุดของเวกเตอร์ซับสเปซก่อนการประมาณค่า DIIS จะเป็น
ดำเนินการ ค่าเริ่มต้นคือ 2 ไม่ควรมีค่าน้อยกว่า 2

DIIS_START = จำนวนเต็ม
จำนวนการวนซ้ำเมื่อ Direct Inversion of the Iterative Subspace (DIIS) ของ Pulay
ควรเปิดขั้นตอนการเร่งการบรรจบกัน ค่าเริ่มต้นคือ 3
ก่อนการวนซ้ำนี้ เวกเตอร์จะไม่ถูกเพิ่มลงในสเปซย่อย DIIS ครั้งแรก
ขั้นตอน DIIS จะไม่ถูกดำเนินการจนกว่า DIIS_MIN_VECS เวกเตอร์อยู่ในซับสเปซ DIIS

LAG_FILE = จำนวนเต็ม
ไฟล์ (หมายเลขหน่วย) สำหรับอ่านเมทริกซ์ลากรองจ์ ค่าเริ่มต้นคือ
ปัจจุบัน 75.

OPDM_FILE = จำนวนเต็ม
ไฟล์ (หมายเลขหน่วย) สำหรับอ่านเมทริกซ์ความหนาแน่นหนึ่งอนุภาค ค่าเริ่มต้น
ปัจจุบันคือ 73

SCALE_GRAD = บูล
บอกว่าจะปรับขนาดการไล่ระดับของวงโคจรโดยประมาณการ (แนวทแยง) orbital
เฮสเซียน ค่าเริ่มต้นคือ TRUE

SCALE_STEP = จริง
ตัวประกอบมาตราส่วนสำหรับขั้นตอนการหมุนของวงโคจร ค่าเริ่มต้นคือ 1.0

TPDM_FILE = จำนวนเต็ม
ไฟล์ (หมายเลขหน่วย) สำหรับอ่านเมทริกซ์ความหนาแน่นสองอนุภาค ค่าเริ่มต้น
ปัจจุบันคือ 74

LEVEL_SHIFT = บูล
บอกว่าจะอนุญาตให้เลื่อนระดับของเมทริกซ์เฮสเซียนหรือไม่ ใช้ได้เลย
เพื่อให้แน่ใจว่า hessian นั้นแน่นอนในเชิงบวกสำหรับการทำซ้ำในตอนเริ่มต้น ถ้า
hessian ไม่เป็นบวกแน่นอนใกล้ convergence แล้วนี่อาจเป็นสัญญาณของ
ความไม่แน่นอนเชิงตัวเลขใน MCSCF ค่าเริ่มต้นคือ TRUE

SHIFT = จริง
บอกว่าองค์ประกอบในแนวทแยงของ hessian ควรเลื่อนระดับเท่าใด
ขยับ ค่าเริ่มต้นคือ 0.01

DETERM_MIN = จริง
ค่าต่ำสุดที่อนุญาตสำหรับดีเทอร์มีแนนต์ของเมทริกซ์เฮสเซียนหากระดับขยับ
เปิดอยู่ ค่านี้จะไม่ถูกใช้หาก LEVEL_SHIFT=FALSE ค่าเริ่มต้นคือ 0.00001

STEP_MAX = จริง
นี่คือการหมุนวงล้อเดี่ยวสูงสุดที่อนุญาต ค่าเริ่มต้นคือ 0.30

บังคับ_ขั้นตอน = บูล
ซึ่งช่วยให้ผู้ใช้สามารถแทนที่ขั้นตอนที่คำนวณได้และบังคับขั้นตอนใน a
ทิศทางเฉพาะ สิ่งนี้มีประโยชน์หากพยายามบังคับการคำนวณออกไป
จากจุดอาน ค่าเริ่มต้นคือ FALSE

FORCE_PAIR = จำนวนเต็ม
นี่คือดัชนีสำหรับคู่อิสระที่จะหมุนถ้า
FORCE_STEP=จริง สิ่งนี้ถูกละเว้นเป็นอย่างอื่น ไม่มีค่าเริ่มต้น

FORCE_VALUE = จริง
นี่คือค่าการหมุนของออร์บิทัลสำหรับการหมุนออร์บิทัลที่ระบุโดย
FORCE_PAIR. ใช้เฉพาะเมื่อตั้งค่า FORCE_STEP เป็น TRUE ค่าเริ่มต้นคือ
0.0.

ตรวจสอบ_HESSIAN = บูล
คำนวณค่าลักษณะเฉพาะของเฮสเซียนออร์บิทัล นี้เป็นสิ่งที่ดีสำหรับการตรวจสอบสำหรับ
โซลูชั่นอานม้า ควรใช้โดยไม่มีการเลื่อนระดับ ค่าเริ่มต้นคือ
FALSE

EIGEN_VECTORS = บูล
คุณต้องการคำนวณเวกเตอร์ลักษณะเฉพาะของเฮสเซียนออร์บิทัลด้วยหรือไม่ นี่คือ
ละเว้นถ้า CHECK_HESSIAN=FALSE ค่าเริ่มต้นคือ FALSE

8 พฤษภาคม 1998 ดีแคส(1)

ใช้ detcas ออนไลน์โดยใช้บริการ onworks.net