ນີ້ແມ່ນຄໍາສັ່ງ detcas ທີ່ສາມາດດໍາເນີນການໄດ້ໃນ OnWorks ຜູ້ໃຫ້ບໍລິການໂຮດຕິ້ງຟຣີໂດຍໃຊ້ຫນຶ່ງໃນຫຼາຍໆບ່ອນເຮັດວຽກອອນໄລນ໌ຂອງພວກເຮົາເຊັ່ນ Ubuntu Online, Fedora Online, Windows online emulator ຫຼື MAC OS online emulator
ໂຄງການ:
NAME
detcas - ການໂຕ້ຕອບການຕັ້ງຄ່າຕົວກໍານົດ / ສໍາເລັດການເຄື່ອນໄຫວ
ໂຄງການພາກສະຫນາມຕົນເອງສອດຄ່ອງກັບອາວະກາດ
ລາຍລະອຽດ
ໂຄງການ detcas ປະຕິບັດພາກສະຫນາມຕົນເອງທີ່ສອດຄ່ອງກັນໃນພື້ນທີ່ຄົບຖ້ວນສົມບູນ (CASSCF)
ການເພີ່ມປະສິດທິພາບຂອງວົງໂຄຈອນໂມເລກຸນໂດຍຜ່ານຂັ້ນຕອນສອງຂັ້ນຕອນໂດຍສົມທົບກັບ
ໂຄງການປະຕິສໍາພັນການຕັ້ງຄ່າຕົວກໍານົດ DETCI. ໂຄງການແມ່ນຂ້ອນຂ້າງງ່າຍດາຍແລະ
ປະຈຸບັນໃຊ້ວິທີການ Newton-Raphson ເພື່ອປັບປຸງວົງໂຄຈອນ, ນໍາໃຊ້ແບບງ່າຍດາຍ
ວົງໂຄຈອນປະມານ Hessian. Convergence ແມ່ນເລັ່ງໂດຍໃຊ້ການປີ້ນໂດຍກົງຂອງ Pulay ຂອງ
ຂັ້ນຕອນການຍ່ອຍສະຫຼາຍຊ້ຳ (DIIS). ລະຫັດໄດ້ຖືກຂຽນເພື່ອອະນຸຍາດໃຫ້ທົ່ວໄປຫຼາຍຂຶ້ນ
wavefunctions ທີ່ບໍ່ຈໍາເປັນຕ້ອງມີການປິ່ນປົວ CI ຢ່າງເຕັມທີ່ຂອງພື້ນທີ່ການເຄື່ອນໄຫວ. ໃນ
ໂດຍສະເພາະ, ການເຄື່ອນໄຫວທີ່ຈໍາກັດໃດໆ (RAS) wavefunction CI ສະຫນັບສະຫນູນໂດຍ DETCI ສາມາດເປັນ
ຖືກນໍາໃຊ້, ອະນຸຍາດໃຫ້ໃຊ້ wavefunctions MCSCF ຂອງປະເພດ RASSCF.
ຂໍ້ມູນອ້າງອິງ
ປະມານວົງໂຄຈອນ Hessian:
1. G. Chaban, MW Schmidt, ແລະ MS Gordon, Theor. ຊິມ. Acta 97, 88-95 (1997).
ພື້ນທີ່ເຄື່ອນໄຫວທີ່ຈຳກັດ CI:
1. ຂັ້ນຕອນການໂຕ້ຕອບການໂຕ້ຕອບຂອງການຕັ້ງຄ່າທີ່ອີງໃສ່ຕົວກໍານົດສໍາລັບການຄົບຖ້ວນສົມບູນແລະຖືກຈໍາກັດ
Configuration Interaction Spaces, J. Olsen, BO Roos, P. Jorgensen, ແລະ HJ Aa.
Jensen, J. Chem. ຟີຊິກ. 89, 2185 (1988).
ພື້ນທີ່ເຄື່ອນໄຫວທີ່ຈຳກັດ SCF:
1. ປ-ກ. Malmqvist, A. Rendell, ແລະ BO Roos, J. Phys. ເຄັມ. 94, 5477 (1990).
ໂຄງການ DETCI:
1. CD Sherrill, ສູດການຄິດໄລ່ການຄິດໄລ່ສຳລັບຂະໜາດໃຫຍ່ເຕັມ ແລະຫຼາຍອ້າງອີງ
ການຕັ້ງຄ່າການໂຕ້ຕອບ Wavefunctions, PhD thesis, University of Georgia, Athens,
GA, 1996.
2. CD Sherrill ແລະ HF Schaefer, ວິທີການປະຕິສໍາພັນການຕັ້ງຄ່າ: ກ້າວຫນ້າ
ໃນວິທີການທີ່ມີຄວາມສໍາພັນສູງ, Adv. Quantum Chem. 34, 143-269 (1999).
ເອກະສານ ຕ້ອງການ
input.dat - ໄຟລ໌ປ້ອນຂໍ້ມູນ
file78 - MO one-electron integratals (fzc operator)
file72 - MO ສອງ-ເອເລັກໂຕຣນິກ
file73 - MO one-particle density matrix
file74 - ມາຕຣິກເບື້ອງຄວາມໜາແໜ້ນສອງອະນຸພາກ MO
file75 - MO Lagrangian
ເອກະສານ ປັບປຸງ
output.dat - ໄຟລ໌ຜົນຜະລິດ
file14.dat - ບັນທຶກຂອງພະລັງງານ ແລະ gradients ຂອງວົງໂຄຈອນ
ປັດໄຈນໍາເຂົ້າ ຮູບແບບ
ການໂຕ້ຖຽງແຖວຄໍາສັ່ງຕໍ່ໄປນີ້ແມ່ນມີຢູ່:
- ງຽບ ນີ້ເຮັດໃຫ້ຜົນໄດ້ຮັບດຽວກັນກັບ ພິມ=0.
-o ຊື່ສຽງ
ໃຫ້ຊື່ໄຟລ໌ສໍາລັບໄຟລ໌ຜົນຜະລິດ. ຄ່າເລີ່ມຕົ້ນຂອງ output.dat.
ການປ້ອນຂໍ້ມູນເພີ່ມເຕີມສໍາລັບໂຄງການນີ້ແມ່ນອ່ານຈາກໄຟລ໌ input.dat. ການນໍາໃຊ້ທົ່ວໄປຫຼາຍ
ຄໍາສໍາຄັນແມ່ນ:
WFN = string
ປະເພດຟັງຊັນຄື້ນທີ່ຮອງຮັບແມ່ນ CASSCF ແລະ RASSCF.
ຄວາມເຊື່ອ ໝັ້ນ = integer
ຕ້ອງການການລວມກັນຢູ່ໃນເສັ້ນໂຄ້ງ gradient. Convergence ແມ່ນບັນລຸໄດ້ໃນເວລາທີ່ RMS
ຄວາມຜິດພາດໃນການໄລ່ສີຂອງວົງໂຄຈອນແມ່ນຫນ້ອຍກ່ວາ 10**(-n). ຄ່າເລີ່ມຕົ້ນແມ່ນ 4 ສໍາລັບ
ການຄິດໄລ່ພະລັງງານ ແລະ 7 ສໍາລັບ gradients.
DOCC = integer_array
vector ນີ້ໃຫ້ຈໍານວນຂອງວົງໂຄຈອນທີ່ຄອບຄອງສອງເທົ່າໃນແຕ່ລະ irrep. ມີ
ບໍ່ມີຄ່າເລີ່ມຕົ້ນ.
SOCC = integer_array
vector ນີ້ໃຫ້ຈໍານວນຂອງ orbitals ຄອບຄອງດຽວໃນແຕ່ລະ irrep. ມີ
ບໍ່ມີຄ່າເລີ່ມຕົ້ນ.
ENERGY_CONVERGENCE = integer
Convergence ທີ່ຕ້ອງການກ່ຽວກັບພະລັງງານ MCSCF ທັງຫມົດ. ຄ່າເລີ່ມຕົ້ນແມ່ນ 7.
FROZEN_DOCC = integer_array
ຈໍານວນຂອງພະລັງງານຕ່ໍາສຸດສອງເທົ່າຄອບຄອງວົງໂຄຈອນໃນແຕ່ລະ irreducible
ການເປັນຕົວແທນທີ່ຈະຖືກແຊ່ແຂງ (ບໍ່ໄດ້ປັບປຸງໃນ MCSCF). ໄດ້
ຄໍາສັ່ງຝ້າຍຂອງຕົວແທນ irredicible ໄດ້ຖືກນໍາໃຊ້. ຄ່າເລີ່ມຕົ້ນແມ່ນ
ສູນ vector.
FROZEN_UOCC = integer_vector
ຈໍານວນຂອງພະລັງງານສູງສຸດວົງໂຄຈອນ unoccupied ໃນແຕ່ລະຕົວແທນ irreducible
ເຊິ່ງຈະຖືກແຊ່ແຂງ (ບໍ່ໄດ້ປັບປຸງໃນ MCSCF). ຄ່າເລີ່ມຕົ້ນແມ່ນສູນ
ເວັກເຕີ
RESTRICTED_DOCC = integer_array
ຈໍານວນຂອງພະລັງງານຕ່ໍາສຸດສອງເທົ່າຄອບຄອງວົງໂຄຈອນໃນແຕ່ລະ irreducible
ການເປັນຕົວແທນທີ່ຈະຖືກປັບປຸງໃຫ້ດີທີ່ສຸດແຕ່ຖືກເກັບຮັກສາສອງເທົ່າຢູ່ໃນ MCSCF.
ວົງໂຄຈອນເຫຼົ່ານີ້ມາຫຼັງຈາກ FROZEN_DOCC ວົງໂຄຈອນ. ຄ່າເລີ່ມຕົ້ນແມ່ນສູນ
ເວັກເຕີ
RESTRICTED_UOCC = integer_array
ຈໍານວນຂອງພະລັງງານສູງສຸດວົງໂຄຈອນ unoccupied ໃນແຕ່ລະຕົວແທນ irreducible
ເຊິ່ງຈະຖືກປັບປຸງໃຫ້ເໝາະສົມ ແຕ່ບໍ່ໄດ້ຢູ່ໃນ MCSCF. ວົງໂຄຈອນເຫຼົ່ານີ້ມາ
ກ່ອນທີ່ຈະໄດ້ FROZEN_UOCC ວົງໂຄຈອນ. ຄ່າເລີ່ມຕົ້ນແມ່ນສູນ vector.
NCASITER = integer
ຈຳນວນສູງສຸດຂອງການເຮັດຊ້ຳເພື່ອປັບແຕ່ງວົງໂຄຈອນ. ທາງເລືອກນີ້ ຄວນ be
ລະບຸ in ໄດ້ ຄ່າເລີ່ມຕົ້ນ ສ່ວນ of ການປ້ອນຂໍ້ມູນ, ເນື່ອງຈາກວ່າມັນຈໍາເປັນຕ້ອງໄດ້ຮັບການສັງເກດເຫັນ
ໂຄງການຄວບຄຸມ PSI. ຄ່າເລີ່ມຕົ້ນເປັນ 1.
PRINT = integer
ທາງເລືອກນີ້ກໍານົດ verbosity ຂອງຜົນຜະລິດໄດ້. ຄ່າຂອງ 1 ຫຼື 2 ລະບຸ
ການພິມຫນ້ອຍທີ່ສຸດ, ມູນຄ່າຂອງ 3 ກໍານົດການພິມ verbose. ຄ່າຂອງ 4 ຫຼື 5 ແມ່ນ
ໃຊ້ສໍາລັບການດີບັກ. ຢ່າໃຊ້ລະດັບ 5 ເວັ້ນເສຍແຕ່ກໍລະນີທົດສອບມີຂະຫນາດນ້ອຍຫຼາຍ (ຕົວຢ່າງ
STO H2O CISD).
ຄໍາທີ່ໃຊ້ຫນ້ອຍທີ່ສຸດແມ່ນ:
DIIS_FREQ = integer
ຂັ້ນຕອນການຂຸດຄົ້ນ DIIS ຈະຖືກພະຍາຍາມທຸກໆຄັ້ງ n ຊ້ຳ. ຄ່າເລີ່ມຕົ້ນ
ແມ່ນ 1.
DIIS_MAX_VECS = integer
ຈຳນວນສູງສຸດຂອງ vectors subspace ສໍາລັບຂັ້ນຕອນ DIIS. ຫຼັງຈາກຕົວເລກນີ້
ຂອງ vectors ແມ່ນບັນລຸໄດ້, vectors ເກົ່າຈະຖືກຫຼຸດລົງຈາກ subspace ເປັນ
ຈໍາເປັນ. ຄ່າເລີ່ມຕົ້ນແມ່ນ 8.
DIIS_MIN_VECS = integer
ຈໍາ ນວນ ຕໍາ ່ ສຸດ ຂອງ vectors subspace ກ່ອນ ທີ່ ຈະ interpolation DIIS ເປັນ
ປະຕິບັດ. ຄ່າເລີ່ມຕົ້ນແມ່ນ 2. ມັນບໍ່ມີຄວາມຫມາຍທີ່ຈະມີຄ່າຫນ້ອຍກວ່າ 2.
DIIS_START = integer
ຕົວເລກຊໍ້າຄືນເມື່ອ Pulay's Direct Inversion of the Iterative Subspace (DIIS)
ຂັ້ນຕອນການສໍາລັບການເລັ່ງຂອງ convergence ຄວນໄດ້ຮັບການເປີດ. ຄ່າເລີ່ມຕົ້ນແມ່ນ 3.
ກ່ອນທີ່ຈະເຮັດຊ້ຳນີ້, vectors ບໍ່ໄດ້ຖືກເພີ່ມໃສ່ພື້ນທີ່ຍ່ອຍ DIIS. ທໍາອິດ
ຂັ້ນຕອນ DIIS ຈະບໍ່ຖືກປະຕິບັດຈົນກ່ວາ DIIS_MIN_VECS vectors ຢູ່ໃນ subspace DIIS.
LAG_FILE = integer
ໄຟລ໌ (ຈໍານວນຫນ່ວຍ) ສໍາລັບການອ່ານ matrix lagrangian. ຄ່າເລີ່ມຕົ້ນແມ່ນ
ປະຈຸບັນ 75.
OPDM_FILE = integer
ໄຟລ໌ (ຈໍານວນຫົວໜ່ວຍ) ສໍາລັບການອ່ານມາຕຣິກເບື້ອງຄວາມຫນາແຫນ້ນຂອງອະນຸພາກ. ຄ່າເລີ່ມຕົ້ນ
ປະຈຸບັນແມ່ນ 73.
SCALE_GRAD = ປຸ້ຍ
ບອກວ່າຈະປັບຂະໜາດຂອງສີໂຄ້ງຕາມວົງໂຄຈອນໂດຍປະມານ (ເສັ້ນຂວາງ) ຫຼືບໍ່
ເຮສຊຽນ. ຄ່າເລີ່ມຕົ້ນແມ່ນ TRUE.
SCALE_STEP = ທີ່ແທ້ຈິງ
ປັດໄຈຂະຫນາດສໍາລັບຂັ້ນຕອນການຫມຸນຂອງວົງໂຄຈອນ. ຄ່າເລີ່ມຕົ້ນແມ່ນ 1.0.
TPDM_FILE = integer
ໄຟລ໌ (ຈໍານວນຫນ່ວຍ) ສໍາລັບການອ່ານມາຕຣິກເບື້ອງຄວາມຫນາແຫນ້ນຂອງສອງອະນຸພາກ. ຄ່າເລີ່ມຕົ້ນ
ປະຈຸບັນແມ່ນ 74.
LEVEL_SHIFT = ປຸ້ຍ
ບອກວ່າຈະອະນຸຍາດໃຫ້ມີການປ່ຽນລະດັບຂອງ hessian matrix ຫຼືບໍ່. ນີ້ສາມາດຖືກນໍາໃຊ້
ເພື່ອຮັບປະກັນວ່າ hessian ມີຄວາມແນ່ນອນໃນທາງບວກສໍາລັບການເລີ່ມຕົ້ນ iterations. ຖ້າ
hessian ແມ່ນບໍ່ແນ່ນອນໃນທາງບວກຢູ່ໃກ້ກັບ convergence, ຫຼັງຈາກນັ້ນນີ້ອາດຈະເປັນສັນຍານຂອງ
ຄວາມບໍ່ສະຖຽນລະພາບທາງດ້ານຕົວເລກໃນ MCSCF. ຄ່າເລີ່ມຕົ້ນແມ່ນ TRUE.
SHIFT = ທີ່ແທ້ຈິງ
ບອກວ່າອົງປະກອບທາງຂວາງຂອງ hessian ຄວນຖືກປ່ຽນເປັນລະດັບເທົ່າໃດ
ການເຄື່ອນຍ້າຍ. ຄ່າເລີ່ມຕົ້ນແມ່ນ 0.01.
DETERM_MIN = ທີ່ແທ້ຈິງ
ຄ່າຕໍ່າສຸດທີ່ອະນຸຍາດໃຫ້ມີຄ່າສໍາລັບຕົວກໍານົດຂອງ hessian matrix ຖ້າລະດັບການປ່ຽນ
ເປີດຢູ່. ອັນນີ້ບໍ່ໄດ້ໃຊ້ຖ້າ LEVEL_SHIFT=FALSE. ຄ່າເລີ່ມຕົ້ນແມ່ນ 0.00001.
STEP_MAX = ທີ່ແທ້ຈິງ
ນີ້ແມ່ນການຫມຸນວົງໂຄຈອນດຽວສູງສຸດທີ່ອະນຸຍາດ. ຄ່າເລີ່ມຕົ້ນແມ່ນ 0.30.
FORCE_STEP = ປຸ້ຍ
ນີ້ອະນຸຍາດໃຫ້ຜູ້ໃຊ້ສາມາດ overide ຂັ້ນຕອນການຄິດໄລ່ແລະບັງຄັບຂັ້ນຕອນໃນ a
ທິດທາງໂດຍສະເພາະ. ນີ້ສາມາດເປັນປະໂຫຍດຖ້າພະຍາຍາມບັງຄັບການຄິດໄລ່ອອກໄປ
ຈາກຈຸດ saddle. ຄ່າເລີ່ມຕົ້ນແມ່ນ FALSE.
FORCE_PAIR = integer
ນີ້ແມ່ນດັດຊະນີສໍາລັບຄູ່ເອກະລາດທີ່ຈະຖືກຫມຸນຖ້າ
FORCE_STEP=TRUE. ອັນນີ້ຖືກລະເລີຍເປັນຢ່າງອື່ນ. ບໍ່ມີຄ່າເລີ່ມຕົ້ນ.
FORCE_VALUE = ທີ່ແທ້ຈິງ
ນີ້ແມ່ນຄ່າການຫມຸນຂອງວົງໂຄຈອນສໍາລັບການຫມຸນວົງໂຄຈອນທີ່ລະບຸໄວ້ໂດຍ
FORCE_PAIR. ອັນນີ້ຖືກໃຊ້ພຽງແຕ່ຖ້າ FORCE_STEP ຖືກຕັ້ງເປັນ TRUE. ຄ່າເລີ່ມຕົ້ນແມ່ນ
0.0.
CHECK_HESSIAN = ປຸ້ຍ
ຄິດໄລ່ຄ່າ eigenvalues ຂອງ hessian ວົງໂຄຈອນ. ນີ້ແມ່ນດີສໍາລັບການກວດສອບ
ການແກ້ໄຂຈຸດ saddle. ຄວນໃຊ້ໂດຍບໍ່ມີການປ່ຽນລະດັບ. ຄ່າເລີ່ມຕົ້ນແມ່ນ
ບໍ່ຖືກຕ້ອງ.
EIGEN_VECTORS = ປຸ້ຍ
ທ່ານຕ້ອງການຄິດໄລ່ eigenvectors ຂອງ orbital hessian ເຊັ່ນກັນ. ນີ້ແມ່ນ
ບໍ່ສົນໃຈຖ້າ CHECK_HESSIAN=FALSE. ຄ່າເລີ່ມຕົ້ນແມ່ນ FALSE.
8 ພຶດສະພາ, 1998 detcas(1)
ໃຊ້ detcas ອອນໄລນ໌ໂດຍໃຊ້ບໍລິການ onworks.net
