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

โครงการ:

ชื่อ

cscf - แก้สมการ Hartree-Fock

DESCRIPTION

โปรแกรม ซีเอสซีเอฟ ดำเนินขั้นตอนวนซ้ำเพื่อแก้สมการ Hartree-Fock

โปรแกรมนี้จำกัดความสมมาตรของ D2h และกลุ่มย่อยและอาชีพการโคจร
จะต้องเป็นจำนวนเต็ม ดังนั้น โมเมนตัมเชิงมุมบริสุทธิ์บางสถานะที่ได้มาจาก
ไม่สามารถรับอาชีพบางส่วนของ orbitals ที่เสื่อมโทรมด้วยรหัสปัจจุบัน สำหรับ
ตัวอย่าง สถานะ 2PIu (doublet PI u) ของ ONO เชิงเส้นที่ได้มาจากพลังงานต่ำสุด
การกำหนดค่าเชิงเส้น (pi u) 1 สามารถคำนวณได้เฉพาะเป็น 2B2u (doublet B2u) หรือ 2B3u
(doublet B 3u) ของสถานะ 2PIu (doublet PI u) และผลลัพธ์เชิงพื้นที่
ฟังก์ชันคลื่นจะไม่มีสมมาตร PI อย่างไรก็ตาม ในแง่หนึ่ง เรื่องนี้เป็นที่พึงปรารถนา
เนื่องจากพลังงานจะเป็นฟังก์ชันต่อเนื่องของมุมดัด การคำนวณพลังงาน
ของการกำหนดค่าที่โค้งงอเป็น 2B2u (doublet B 2u) หรือ 2B3u (doublet B 3u) และทำ 2PIu ล้วนๆ
(doublet PI u) ที่รูปทรงเชิงเส้นส่งผลให้เกิดความไม่ต่อเนื่องที่เด่นชัด

โดยส่วนใหญ่ สภาพแฝดที่เกิดจากการประกอบอาชีพสองครั้งของความเสื่อมทวีคูณ
ออร์บิทัล เช่น สถานะ 3A2 (ทริปเล็ต A 2) ที่เกิดจาก (e')2 หรือ (e")2
การกำหนดค่าในสมมาตร D3h หรือสถานะ 3SIGMAg (triplet SIGMA g) ของ a (pi g)2 หรือ (pi
u) การกำหนดค่า 2 ในสมมาตร Dinfh (D อินฟินิตี้ h) จะมีสมมาตรเชิงพื้นที่ที่เหมาะสม
สถานะเสื้อกล้ามที่เกิดจากการกำหนดค่าทางอิเล็กทรอนิกส์เดียวกันเหล่านี้มีอยู่โดยเนื้อแท้
การกำหนดค่าหลายรายการและด้วยเหตุนี้จึงไม่ได้แสดงด้วยการกำหนดค่าเดียว
ฟังก์ชั่นคลื่น

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

วิธีไฟล์ PK:

1. RC Raffenetti, เคม. สรีรวิทยา เลตต์. 20 (1973) 335.

สมมาตรโมเลกุลและการคำนวณ HF ของเปลือกปิด:

1. M.Dupuis และ HFKing, Int. เจ. ควอน. เคมี. 11 (1977) 613.

DIIS สำหรับเปลือกปิด:

1. ป. ปูเล, เคมี. สรีรวิทยา เลตต์. 73 (1980) 393.

2. ป. ปูเล, เจ. คอม เคมี. 3 (1982) 556.

ค่าสัมประสิทธิ์การมีเพศสัมพันธ์ (อัลฟาและเบต้า) สำหรับโอเพ่นเชลล์:

1. CCJ Roothaan รายได้ Mod สรีรวิทยา 32 (1960) 179.

การทำให้หมาด ๆ :

1. DR Hartree "การคำนวณโครงสร้างอะตอม" (Wiley: New York) 1957

2. MC Zerner และ M. Hehenberger, Chem. สรีรวิทยา เลตต์. 62 (1979) 550.

การเลื่อนระดับ:

1. VR Saunders และ IH Hillier, Int. เจ. ควอน. เคมี. 7 (1973) 699.

บรรจบกัน ก.พ

สำหรับกรณีเปลือกเปิดที่ยาก ขอแนะนำให้ใช้เปลือกปิดที่เหมาะสม
ให้รันการคำนวณก่อน (เพิ่มหรือลบอิเล็กตรอนส่วนเกิน) และเวกเตอร์ SCF นี้แล้ว
ใช้เป็นค่าเดาสำหรับฟังก์ชันคลื่นเปิดเปลือกที่ต้องการ สำหรับกรณี TCSCF จะเป็นเสมอ
ฉลาดที่จะเรียกใช้เชลล์ปิด (หรือบางทีอาจเป็นแฝดสามที่เหมาะสม) SCF ก่อนแล้วจึงใช้
นี่เป็นการเดาสำหรับ TCSCF

สำหรับระบบ open shell แนะนำให้เลื่อนระดับ 0.5 ถึง 3.0 เริ่มด้วย
มูลค่าสูง (2.0 - 3.0) สำหรับการคำนวณ SCF แรกแล้วลด (เป็น 0.5 - 1.0) สำหรับ
การวิ่งครั้งต่อๆ มาซึ่งใช้เวกเตอร์ SCF ที่บรรจบกันเป็นจุดเริ่มต้น

เป็นสิ่งสำคัญอย่างยิ่งที่จะต้องทราบว่าโค้ดเวอร์ชันนี้ไม่รองรับอีกต่อไป
แบบเปิด หนึ่ง ต้อง ใช้ ใหม่ คำหลัก อ้างอิง และ มัลติเพิล ไปยัง ระบุ ชนิด of SCF
จำเป็น

INPUT FORMAT

การขอ ซีเอสซีเอฟ โปรแกรมค้นหาผ่านเส้นทางคีย์เวิร์ดเริ่มต้น (first SCF แล้วก็ เริ่มต้น)
สำหรับคำหลักต่อไปนี้:

LABEL = เชือก
นี่คือสตริงอักขระที่จะรวมไว้ในเอาต์พุต สตริงนี้ไม่ได้ใช้
โดยโปรแกรม ไม่มีค่าเริ่มต้น

ดับเบิลยูเอฟเอ็น = เชือก
นี่คือประเภทของฟังก์ชันคลื่นที่ต้องการในท้ายที่สุด ค่าเริ่มต้นคือ SCF.

โอเพนไทป์ is ไม่ อีกต่อไป ได้รับการสนับสนุน

อ้างอิง = เชือก
ระบุประเภทของการคำนวณ SCF ที่ต้องการทำ สามารถเป็นหนึ่งใน อาร์เอชเอฟ
(สำหรับเสื้อกล้ามเปลือกปิด) โรห์ฟ (สำหรับการคำนวณ open shell ที่จำกัด) UHF
(สำหรับการคำนวณ open shell ที่ไม่จำกัด) ทูคอน (สำหรับการกำหนดค่าสองแบบ
เสื้อกล้าม) หรือ พิเศษ. ถ้า พิเศษ จะได้รับ จากนั้นคัปปลิ้งอัลฟ่าและเบตา
ต้องให้สัมประสิทธิ์กับ ALPHA และ เบต้า คีย์เวิร์ด ค่าเริ่มต้นคือ อาร์เอชเอฟ.

มัลติเพิล= จำนวนเต็ม
ระบุหลายหลากของโมเลกุล ค่าเริ่มต้นคือเสื้อกล้าม

ค่าใช้จ่าย= จำนวนเต็ม
ระบุประจุของโมเลกุล ผิดนัดคือ 0

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

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

เดอร์ไทป์ = เชือก
ระบุลำดับของอนุพันธ์ที่ต้องทำในที่สุด มันถูกใช้
โดย SCF โปรแกรมเพื่อกำหนดว่าต้องเขียนไฟล์บางไฟล์หรือไม่และยังเป็น
ใช้เพื่อกำหนดคอนเวอร์เจนซ์เริ่มต้นของฟังก์ชันคลื่น ค่าเริ่มต้นคือ
FIRST.

แม็กซิเตอร์ = จำนวนเต็ม
สิ่งนี้ทำให้จำนวนการวนซ้ำสูงสุด ค่าเริ่มต้นคือ 40

คอนเวอร์เจนซ์ = จำนวนเต็ม
สิ่งนี้ระบุว่าจะหลอมรวมฟังก์ชันคลื่นให้แน่นเพียงใด การบรรจบกันคือ
กำหนดโดยการเปรียบเทียบการเปลี่ยนแปลง RMS ในเมทริกซ์ความหนาแน่น ("เดลต้า P") กับ
ค่าที่กำหนด เกณฑ์การบรรจบกันคือ 10**(-จำนวนเต็ม). ค่าเริ่มต้นคือ 7 if
ทั้งสอง เดอร์ไทป์ = NONE และ ดับเบิลยูเอฟเอ็น = SCF จะได้รับและ 10 อย่างอื่น

การเปลี่ยนระดับ = จริง
สิ่งนี้ระบุการเลื่อนระดับ ค่าเริ่มต้นคือ 1

โดยตรง = บูล
ระบุว่าจะทำการคำนวณ SCF ด้วยเทคนิคอินทิกรัลตรงหรือไม่ NS
ค่าเริ่มต้นเป็นเท็จ

PRINT_MOS = บูล
ระบุว่าจะพิมพ์โมเลกุลออร์บิทัลหรือไม่ ค่าเริ่มต้นเป็นเท็จ

นอกจากนี้ยังมีพารามิเตอร์อินพุตที่ใช้กันน้อยกว่าจำนวนมาก ถ้าคุณไม่
ทำความเข้าใจความหมายของตัวเลือกต่อไปนี้ จากนั้นตรวจสอบให้แน่ใจว่าตัวเลือกเหล่านี้ไม่ปรากฏใน
ป้อนข้อมูล. ค่าเริ่มต้นจะใช้ได้ในกรณีส่วนใหญ่ เหล่านี้ระบุไว้
ด้วยคำสำคัญดังนี้

DELETE_INTS = บูล
ไฟล์อินทิกรัลจะถูกลบหาก ดับเบิลยูเอฟเอ็น = SCF และ เดอร์ไทป์ = FIRST or เดอร์ไทป์ = NONE.
หากคุณต้องการเก็บไฟล์อินทิกรัลไว้ ให้ตั้งค่า DELETE_INTS = เท็จ ค่าเริ่มต้นคือ
จริง

เรียงลำดับใหม่ = เชือก
พารามิเตอร์ควบคุมการเรียงลำดับใหม่ของออร์บิทัลระดับโมเลกุล หากตั้งค่าเป็น ก่อน แล้ว
เดา orbitals จากไฟล์ด่านจะถูกจัดลำดับใหม่ หากตั้งค่าเป็น หลัง, บรรจบกัน
ออร์บิทัลจะถูกจัดลำดับใหม่ก่อนที่จะเขียนไปยังไฟล์จุดตรวจ ไม่ว่าใน
กรณี มัวร์เดอร์ ต้องกำหนดพารามิเตอร์เพื่อระบุการเรียงลำดับแผนที่ใหม่ ค่าเริ่มต้นคือ
ไม่ต้องเรียงลำดับออร์บิทัลใหม่

มัวร์เดอร์ = จำนวนเต็ม_เวกเตอร์
นี่ระบุเวกเตอร์การจัดลำดับโมเลกุลการโคจรใหม่ จะใช้ก็ต่อเมื่อ
เรียงลำดับใหม่ ถูกตั้งค่า เวกเตอร์นี้แมปทุกออร์บิทัลกับดัชนีใหม่ของมัน เช่น มัวร์เดอร์ = (0
2 1) ระบุว่าหลังจากการเรียงลำดับ orbitals 1 และ 2 ใหม่จะถูกสลับ อันดับของ
เวกเตอร์นี้เท่ากับจำนวน MO ดัชนีอยู่ในลำดับพิทเซอร์
(เรียงลำดับตามความสมมาตร ตามด้วยพลังงานภายในแต่ละบล็อกสมมาตร) ฐาน-0 CSCF
มีแนวโน้มว่าจะล้มเหลวหาก MOORDER ที่กำหนดผสม orbitals จาก irreps ที่แตกต่างกัน ที่นั่น
ไม่มีค่าเริ่มต้น

ALPHA = จริง_เวกเตอร์
If โอเพนไทป์ = พิเศษจากนั้นพารามิเตอร์นี้จะให้ค่าสัมประสิทธิ์การมีเพศสัมพันธ์
จำนวนขององค์ประกอบในเวกเตอร์นี้คือ MM(MM+)/1 โดยที่ MM คือจำนวนของ
การแสดงที่ลดไม่ได้ซึ่งมีออร์บิทัลโมเลกุลที่ถูกครอบครองโดยลำพัง ที่นั่น
ไม่มีค่าเริ่มต้น

เบต้า = จริง_เวกเตอร์
If โอเพนไทป์ = พิเศษจากนั้นพารามิเตอร์นี้จะให้ค่าสัมประสิทธิ์การมีเพศสัมพันธ์แบบเบตา
จำนวนขององค์ประกอบในเวกเตอร์นี้คือ MM(MM+)/1 โดยที่ MM คือจำนวนของ
การแสดงที่ลดไม่ได้ซึ่งมีออร์บิทัลโมเลกุลที่ถูกครอบครองโดยลำพัง ที่นั่น
ไม่มีค่าเริ่มต้น

GUESS = เชือก
ตัวเลือกนี้กำหนดประเภทของการเดาเบื้องต้นที่ eigenvector CSCF จะใช้
ตัวเลือกเดียวที่ใช้ได้ในขณะนี้คือ : (1) GUESS = COREซึ่งทำให้การใช้งาน
core Hamiltonian eigenvector เพื่อเริ่มการคำนวณ (2) GUESS = AUTO ที่
ส่งผลให้เกิดความพยายามที่จะใช้เวกเตอร์ MO ในไฟล์จุดตรวจหรือหันไป
หลักเดาถ้าไม่มี eigenvector ในไฟล์นั้น ค่าเริ่มต้น if AUTO.

พิมพ์ = จำนวนเต็ม
นี่คือตัวเลือกการพิมพ์ ค่าเริ่มต้นคือ 0

MO_OUT = บูล
พิมพ์ออร์บิทัลด้วยความสมมาตรและอาชีพที่ส่วนท้ายของ
การคำนวณ ค่าเริ่มต้นเป็นจริง

หมุน = บูล
ออร์บิทัลของโมเลกุลจะไม่หมุนหากเป็นเท็จ การหมุนเท่านั้น
ส่งผลกระทบต่อออร์บิทัลเสมือนสำหรับระบบโอเพ่นเชลล์ พารามิเตอร์นี้จะต้องเป็นจริง
สำหรับการไล่ระดับสีที่มีความสัมพันธ์และต้องเป็นเท็จสำหรับอนุพันธ์อันดับสองและสูงกว่า
ค่าเริ่มต้นเป็นเท็จ if ดับเบิลยูเอฟเอ็น = SCF และจริงเป็นอย่างอื่น

ตรวจสอบ_ROT = บูล
ตรวจสอบการหมุนของโมเลกุลที่อธิบายข้างต้นเพื่อให้แน่ใจว่าไม่มีคอลัมน์ของ
เมทริกซ์ไอเกนเวคเตอร์ SCF ถูกสลับโดยการหมุน ไม่มีผลถ้า หมุน =
เท็จ. ค่าเริ่มต้นเป็นจริง

CHECK_MO_ORTHOGONALITY = บูล
ตรวจสอบว่าออร์บิทัลของโมเลกุลเป็นแบบออร์โธนอร์มอลหรือไม่ มีประโยชน์สำหรับการดีบักเท่านั้น NS
ค่าเริ่มต้นเป็นเท็จ

ดีไอเอส = บูล
สิ่งนี้กำหนดว่าจะใช้ diis หรือไม่ ค่าเริ่มต้นเป็นจริง

DISSTART = จำนวนเต็ม
นี่เป็นการวนซ้ำครั้งแรกที่จะใช้ DIIS ค่าเริ่มต้นคือ 0

เอ็นดีไอเอส = จำนวนเต็ม
ซึ่งจะทำให้จำนวนของเมทริกซ์ข้อผิดพลาดที่จะใช้ในขั้นตอน diis ค่าเริ่มต้น
คือ 6 สำหรับเชลล์แบบปิด 4 สำหรับเชลล์แบบเปิด และ 3 สำหรับ tcscf

ไดสดัม = จริง
สิ่งนี้ให้ปัจจัยการหน่วงสำหรับขั้นตอน diis ค่าเริ่มต้นคือ 0.0 สำหรับ
เชลล์ปิด 0.02 สำหรับเชลล์เปิด และ 0.01 สำหรับ tcscf

INCR = จริง
ค่านี้ใช้ใน tcscf เพื่อกำหนดความถี่ในการคำนวณค่าสัมประสิทธิ์ ci
จำนวนเล็กน้อย (~0.25) จะทำให้คำนวณใหม่เกือบทุก scf
การวนซ้ำ ค่าเริ่มต้นคือ 0.25

DYN_ACC = บูล
เมื่อดำเนินการ scf โดยตรง สิ่งนี้จะระบุว่าจะตัดความแม่นยำของอินทิกรัลไดนามิกหรือไม่
จะถูกนำไปใช้. ค่าดีฟอลต์เป็นจริง (ใช้ไดนามิกคัทออฟ) การวนซ้ำเบื้องต้นคือ
ดำเนินการด้วยอินทิกรัลที่ถูกต้องถึงหกหลัก หลังจากความหนาแน่นถูกหลอมรวมเป็น
ทำซ้ำ 10^-5 หรือ 30 ครั้ง ใช้ความถูกต้องครบถ้วนสมบูรณ์ ถ้าscf
ปัญหาการบรรจบกันมีประสบการณ์ในการปิดการใช้งานไดนามิกคัทออฟโดยการตั้งค่านี้
ตัวแปรเป็นเท็จอาจช่วยได้

ORTHOG_เท่านั้น = บูล
บางครั้งใน CASSCF หรือรูปแบบอื่นที่ไม่ใช่ HF/KS สำหรับการเพิ่มประสิทธิภาพการโคจร มันคือ
มีประโยชน์ในการจัดแนว MO ใหม่จากรูปทรงอื่นสำหรับเรขาคณิตปัจจุบัน so
สามารถใช้เป็นเดาเบื้องต้นสำหรับ MO ใหม่ได้ สามารถทำได้โดย
เรียกใช้ CSCF ด้วย ORTHOG_เท่านั้น = จริง. หลังจากที่ออร์บิทัลถูกปรับตั้งฉากแล้ว
โปรแกรมจะปิดโดยไม่ทำการคำนวณ SCF คีย์เวิร์ดนี้จะเป็น
ละเว้นหากไม่มีออร์บิทัลก่อนหน้าในไฟล์จุดตรวจ ค่าเริ่มต้นเป็น จริง
if ดับเบิลยูเอฟเอ็น = ดีแคส.

30 พฤษภาคม 1991 ซีเอสซีเอฟ(1)

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