detcas - 云端在线

程序：

您的姓名

detcas - 行列式配置交互/完全活动
空间自洽现场计划

商品描述

该程序 德卡斯 执行完全活动空间自洽场 (CASSCF)
通过两步程序优化分子轨道
行列式配置交互程序 DETCI。 该程序相当简单
目前使用 Newton-Raphson 方法来更新轨道，采用简单的
近似轨道 Hessian。 使用 Pulay 的直接反演加速收敛
迭代子空间 (DIIS) 过程。 已编写代码以允许更一般的
不一定具有对活动空间的完整 CI 处理的波函数。 在
特别是，DETCI 支持的任何受限活动空间 (RAS) CI 波函数都可以是
使用，允许 RASSCF 类型的 MCSCF 波函数。

参考文献：

近似轨道 Hessian：

1. G. Chaban、MW Schmidt 和 MS Gordon、Theor。 哼。 学报 97, 88-95 (1997)。

受限活动空间 CI：

1. 基于行列式的完全和受限配置交互算法
配置交互空间，J. Olsen、BO Roos、P. Jorgensen 和 HJ Aa。
詹森，J. 化学。 物理89, 2185 (1988)。

受限活动空间 SCF：

1. P.-A. Malmqvist、A. Rendell 和 BO Roos、J. Phys。 化学94, 5477 (1990)。

DETCI 计划：

1. CD Sherrill，用于大规模全参考和多参考的计算算法
配置相互作用波函数，博士论文，乔治亚大学，雅典，
GA，1996 年。

2. CD Sherrill 和 HF Schaefer，配置交互方法：进展
在高度相关的方法中，Adv。 量子化学34, 143-269 (1999)。

FILES 所需

input.dat - 输入文件
file78 - MO 单电子积分（fzc 算子）
file72 - MO 二电子积分
file73 - MO 单粒子密度矩阵
file74 - MO 二粒子密度矩阵
file75 - MO 拉格朗日

FILES 已更新

output.dat - 输出文件
file14.dat - 能量和轨道梯度的记录

INPUT FORMAT

以下命令行参数可用：

-安静的 这给出了相同的结果 打印=0.

-o 名
给出输出文件的文件名。 默认为 output.dat。

该程序的附加输入是从文件 input.dat 中读取的。 比较常用的
关键词是：

WFN更多 = 绳子
支持的波函数类型为 CASSCF 和 RASSCF。

收敛 = 整数
轨道梯度上所需的收敛。 当 RMS 达到收敛
的轨道梯度误差小于 10**(-n)。 默认值为 4
能量计算和 7 梯度。

DOCC = 整数数组
这个向量给出了每个不可逆的双占用轨道的数量。 有
没有默认。

SOCC = 整数数组
这个向量给出了每个不可逆的单占轨道的数量。 有
没有默认。

能源汇聚 = 整数
总 MCSCF 能量所需的收敛。 默认值为 7。

FROZEN_DOCC 文件 = 整数数组
每个不可约中最低能量双占据轨道的数量
将被冻结的表示（未在 MCSCF 中更新）。 这
使用不可分割的表示的棉花排序。 默认是
零向量。

冰雪奇缘_UOCC = 整数向量
每个不可约表示中最高能量未占据轨道的数量
这将被冻结（未在 MCSCF 中更新）。 默认为零
向量。

RESTRICTED_DOCC 文件 = 整数数组
每个不可约中最低能量双占据轨道的数量
表示将被优化但在 MCSCF 中保持双重占用。
这些轨道出现在 FROZEN_DOCC 文件 轨道。 默认为零
向量。

RESTRICTED_UOCC = 整数数组
每个不可约表示中最高能量未占据轨道的数量
这将被优化但在 MCSCF 中保持空闲。 这些轨道来
前 冰雪奇缘_UOCC 轨道。 默认值是零向量。

卡西特 = 整数
优化轨道的最大迭代次数。 这个选项 应该 be
指定 in 此 DEFAULT 部分 of 输入，因为它需要对
控制程序 PSI。 默认为 1。

打印 = 整数
此选项确定输出的详细程度。 1 或 2 的值指定
最小打印，值为 3 指定详细打印。 4 或 5 的值是
用于调试。 不要使用级别 5，除非测试用例非常小（例如
STO H2O CISD）。

不太常用的关键字是：

DIIS_频率 = 整数
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 子空间中。

滞后文件 = 整数
用于读取拉格朗日矩阵的文件（单元号）。 默认值为
目前 75。

OPDM_文件 = 整数
用于读取单粒子密度矩阵的文件（单位编号）。 默认值
目前是 73。

SCALE_GRAD 等级 = 布尔
告诉是否按近似（对角线）轨道缩放轨道梯度
黑森州。 默认值为真。

比例尺 = 真实
轨道旋转步骤的比例因子。 默认值为 1.0。

TPDM_文件 = 整数
用于读取两粒子密度矩阵的文件（单位编号）。 默认值
目前是 74。

水平移动 = 布尔
告诉是否允许 hessian 矩阵的水平移动。 这个可以用
以确保 hessian 对于开始迭代是正定的。 如果
Hessian 不是正定接近收敛，那么这可能是
MCSCF 中的数值不稳定性。 默认值为真。

SHIFT = 真实
告诉粗麻布的对角线元素应该移动多少水平
转移。 默认值为 0.01。

DETERM_MIN = 真实
如果水平移动，hessian 矩阵的行列式的最小允许值
正在。 如果 LEVEL_SHIFT=FALSE，则不使用此选项。 默认值为 0.00001。

STEP_MAX = 真实
这是允许的最大单轨道旋转。 默认值为 0.30。

力步 = 布尔
这允许用户覆盖计算的步骤并强制执行一个步骤
特定的方向。 如果试图强制计算，这可能很有用
从鞍点。 默认值为 FALSE。

力对 = 整数
这是要旋转的独立对的索引，如果
FORCE_STEP=真。 否则将被忽略。 没有默认值。

力值 = 真实
这是旋转由指定的轨道的轨道旋转值
FORCE_PAIR。 这仅在 FORCE_STEP 设置为 TRUE 时使用。 默认值为
0.0.

检查_黑森 = 布尔
计算轨道 hessian 的特征值。 这有利于检查
鞍点解。 应在没有电平转换的情况下使用。 默认是
假。

特征向量 = 布尔
您是否也想计算轨道 hessian 的特征向量。 这是
如果 CHECK_HESSIAN=FALSE，则忽略。 默认值为 FALSE。

<span>2017-05-22</span> 德卡斯(1)

使用 onworks.net 服务在线使用 detcas