cp2k inp帮助
- EXT_RESTART
-
- RESTART_FILE_NAME
- RESTART_VEL
- GLOBAL
-
-
- PHYSCON
- PRINT_LEVEL
-
- LOW
- PROGRAM
- PROGRAM_NAME
-
- ATOM
- CP2K
- FARMING
- MC_ANALYSIS
- OPTIMIZE_BASIS
- OPTIMIZE_INPUT
- SWARM
- TEST
- TMC
- PROJECT_NAME
- RUN_TYPE
-
- BAND
- ENERGY
- GEO_OPT
- MD
- NONE
- VIBRATIONAL_ANALYSIS
- SEED
- TRACE
- FORCE_EVAL
-
- DFT
-
- BASIS_SET_FILE_NAME
- CHARGE
- LSD
- MGRID
-
- CUTOFF
- NGRIDS
- REL_CUTOFF
- RELATIVE_CUTOFF
- POISSON
- POTENTIAL_FILE_NAME
-
- E_DENSITY_CUBE
- LOWDIN
- MO_CUBES
-
- NHOMO
- NLUMO
- MULLIKEN
- QS
-
- CDFT
-
- ATOM_GROUP
-
- ATOMS
- COEFF
- CONSTRAINT_TYPE
- ATOMIC_CHARGES
- BECKE_CONSTRAINT
-
- ADJUST_SIZE
- ATOMIC_RADII
- CAVITY_CONFINE
- CAVITY_SHAPE
- CUTOFF_TYPE
- ELEMENT_CUTOFF
- EPS_CAVITY
- IN_MEMORY
- SHOULD_SKIP
- DUMMY_ATOMS
-
- ATOMS
- OUTER_SCF
- PROGRAM_RUN_INFO
- STRENGTH
- TARGET
- TYPE_OF_CONSTRAINT
-
- BECKE
- HIRSHFELD
- NONE
- EPS_DEFAULT
- EXTRAPOLATION
-
- ASPC
- PS
- EXTRAPOLATION_ORDER
- METHOD
-
- AM1
- DFTB
- GAPW
- GAPW_XC
- GPW
- LRIGPW
- MNDO
- MNDOD
- OFGPW
- PDG
- PM3
- PM6
- PM6-FM
- PNNL
- RIGPW
- RM1
- XTB
- SE
-
- COULOMB
-
- CUTOFF
- EXCHANGE
-
- CUTOFF
- RESTART_FILE_NAME
- SCF
-
- CHOLESKY
-
- INVERSE
- INVERSE_DBCSR
- OFF
- REDUCE
- RESTORE
- EPS_LUMO
- EPS_SCF
- MAX_ITER_LUMOS
- MAX_SCF
- OT
-
- ALGORITHM
-
- IRAC
- STRICT
- MINIMIZER
-
- CG
- DIIS
- N_BROYDEN
- N_DIIS
- N_HISTORY_VEC
- NDIIS
- PRECONDITIONER
-
- FULL_ALL
- FULL_KINETIC
- FULL_SINGLE
- FULL_SINGLE_INVERSE
- FULL_S_INVERSE
- NONE
- OUTER_SCF
-
- CDFT_OPT
-
- CONTINUE_LS
- FACTOR_LS
- JACOBIAN_FREQ
- JACOBIAN_RESTART
- JACOBIAN_STEP
- JACOBIAN_TYPE
- MAX_LS
- EPS_SCF
- EXTRAPOLATION_ORDER
- MAX_SCF
- OPTIMIZER
-
- BISECT
- BROYDEN
- DIIS
- NEWTON
- NEWTON_LS
- NONE
- SD
- SECANT
- STEP_SIZE
- TYPE
-
- BASIS_CENTER_OPT
- CDFT_CONSTRAINT
- DDAPC_CONSTRAINT
- S2_CONSTRAINT
- NONE
-
- RESTART
- SCF_GUESS
-
- ATOMIC
- RESTART
- SPIN_POLARIZED
- UKS
- UNRESTRICTED_KOHN_SHAM
- WFN_RESTART_FILE_NAME
- XC
-
- HF
-
- FRACTION
- INTERACTION_POTENTIAL
-
- CUTOFF_RADIUS
- POTENTIAL_TYPE
- T_C_G_DATA
- MEMORY
- SCREENING
-
- EPS_SCHWARZ
- SCREEN_ON_INITIAL_P
- VDW_POTENTIAL
-
- POTENTIAL_TYPE
-
- NONE
- NON_LOCAL
- PAIR_POTENTIAL
- PAIR_POTENTIAL
-
- CALCULATE_C9_TERM
- LONG_RANGE_CORRECTION
- PARAMETER_FILE_NAME
- REFERENCE_C9_TERM
- REFERENCE_FUNCTIONAL
- TYPE
- XC_FUNCTIONAL
-
- 参数
-
- B3LYP
- BEEFVDW
- BLYP
- BP
- HCTH120
- LDA
- NONE
- NO_SHORTCUT
- OLYP
- PADE
- PBE
- PBE0
- TPSS
- PBE
-
- SCALE_C
- SCALE_X
- PBE_HOLE_T_C_LR
-
- CUTOFF_RADIUS
- SCALE_X
- XC_GRID
-
- XC_DERIV
-
- SPLINE2_SMOOTH
- XC_SMOOTH_RHO
-
- NN10
- METHOD
-
- FIST
- MIXED
- QuickStep
- QS
- MIXED
-
- GENERIC
-
- ERROR_LIMIT
- MIXING_FUNCTION
- VARIABLES
- GROUP_PARTITION
- MAPPING
-
- FORCE_EVAL
-
- DEFINE_FRAGMENTS
- FORCE_EVAL_MIXED
-
- FRAGMENT
- MIXING_TYPE
-
- GENMIX
- LINEAR_COMBINATION
- MM
-
- FORCEFIELD
-
- BEND
-
- ATOMS
- K
- KIND
- THETA0
- BOND
-
- ATOMS
- K
- KIND
-
- CHARMM
- HARMONIC
- R0
- CHARGE
-
- ATOM
- CHARGE
- EI_SCALE14
- PARM_FILE_NAME
- parmtype
-
- AMBER
- CHM
- NONBONDED
-
- EAM
-
- ATOMS
- PARM_FILE_NAME
- GENPOT
-
- ATOMS
- FUNCTION
- PARAMETERS
- RCUT
- VALUES
- VARIABLES
- LENNARD-JONES
-
- ATOMS
- EPSILON
- SIGMA
- RCUT
- SPLINE
-
- EMAX_SPLINE
- EPS_SPLINE
- VDW_SCALE14
-
- DIPOLE
-
- EACH
- FILENAME
- PERIODIC
- FF_INFO
- FF_PARAMETER_FILE
- POISSON
-
- EWALD
-
- EWALD_TYPE
-
- NONE
- PME
- SPME
- ALPHA
- GMAX
- O_SPLINE
- MT
-
- ALPHA
- REL_CUTOFF
- PERIODIC
-
- NONE
- POISSON
- POISSON_SOLVER
-
- PERIODIC
- WAVELET
- PSOLVER
-
- FORCES
- GRID_INFORMATION
- PROPERTIES
-
- RESP
-
- CONSTRAINT
-
- ATOM_LIST
- EQUAL_CHARGES
-
- V_RESP_CUBE
- SPHERE_SAMPLING
-
- RMAX_KIND
- RMIN_KIND
- STRIDE
- STRESS_TENSOR
-
- ANALYTICAL
- SUBSYS
-
- CELL
-
- ABC
- MULTIPLE_UNIT_CELL
- PERIODIC
- COLVAR
-
- TORSION
-
- ATOMS
- COORD
-
- UNIT
-
- angstrom
- KIND
-
- BASIS_SET
- POTENTIAL
-
- TOPOLOGY_INFO
-
- AMBER_INFO
- TOPOLOGY
-
- CENTER_COORDINATES
- CHARGE_BETA
- CONN_FILE_FORMAT
-
- AMBER
- GENERATE
- OFF
- PSF
- CONN_FILE_NAME
- CONNECTIVITY
- COORD_FILE_NAME
- COORD_FILE_FORMAT
-
- CIF
- CP2K
- CRD
- G96
- OFF
- PDB
- XTL
- XYZ
- COORDINATE
- DUMP_PDB
- DUMP_PSF
- GENERATE
-
- BONDLENGTH_MAX
- MULTIPLE_UNIT_CELL
- MOTION
-
- BAND
-
- BAND_TYPE
-
- IT-NEB
- NPROC_REP
- NUMBER_OF_REPLICA
- OPTIMIZE_BAND
-
- DIIS
-
- MAX_STEPS
- N_DIIS
- OPT_TYPE
-
- DIIS
- REPLICA
-
- COORD
- CONSTRAINT
-
- COLLECTIVE
-
- COLVAR
- INTERMOLECULAR
-
- T
- F
- RESTRAINT
-
- K
- TARGET
- FIXED_ATOMS
-
- COMPONENTS_TO_FIX
-
- Z
- LIST
- GEO_OPT
-
- OPTIMIZER
-
- BFGS
- CG
- LBFGS
- MAX_FORCE
- MAX_ITER
- RMS_FORCE
- MD
-
- ANGVEL_TOL
- ANGVEL_ZERO
- BAROSTAT
-
- PRESSURE
- TIMECON
- COMVEL_TOL
- ENSEMBLE
-
- NPE_I
- NVE
- NVT
- RESPA
- STEPS
- STEP_START_VAL
- TIMESTEP
- TEMPERATURE
- THERMOSTAT
-
- NOSE
-
- LENGTH
- MTS
- TIMECON
- YOSHIDA
-
- ENERGY
-
- EACH
- FORCES
- PROGRAM_RUN_INFO
-
- ADD_LAST
-
- NO
- NUMERIC
- SYMBOLIC
- COMMON_ITERATION_LEVELS
- EACH
- RESTART
- RESTART_HISTORY
- TRAJECTORY
-
- ADD_LAST
-
- NUMERIC
- EACH
- FILENAME
- FORMAT
-
- ATOMIC
- DCD
- DCD_ALIGNED_CELL
- PDB
- XMOL
- XYZ
- VELOCITIES
- MULTIPLE_FORCE_EVALS
-
- FORCE_EVAL_ORDER
- MULTIPLE_SUBSYS
- OPTIMIZE_INPUT
-
- ACCURACY
- FORCE_MATCHING
-
- COMPARE_ENERGIES
-
- EACH
-
- GEO_OPT
- MD
- POWELL_OPT
- QS_SCF
- ENERGY_WEIGHT
- FRAME_START
- GROUP_SIZE
- OPTIMIZE_FILE_NAME
- REF_FORCE_FILE_NAME
- REF_TRAJ_FILE_NAME
- SHIFT_QM
- SHIFT_MM
- HISTORY
-
- EACH
- MAX_FUN
- METHOD
-
- FORCE_MATCHING
- RESTART
-
- BACKUP_COPIES
- COMMON_ITERATION_LEVELS
- EACH
- FILENAME
- STEP_SIZE
- VARIABLE
-
- FIXED
- LABEL
- VALUE
EXT_RESTART
RESTART_FILE_NAME
指定要读取的重新启动文件(或任何其他输入文件)的名称。仅使用与重新启动相关的字段(除非使用本节中的关键字关闭)
RESTART_VEL
从外部文件获取速度 默认值:true
GLOBAL
全局信息
控制物理常数和数学常数的打印
PHYSCON
如果printkey处于活动状态,则打印物理常数 默认值:true
PRINT_LEVEL
打印级别,默认MEDIUM
LOW
打印很少信息
PROGRAM
同PROGRAM_NAME
PROGRAM_NAME
指定运行哪个程序 默认值:CP2K
ATOM
运行单原子计算
CP2K
运行一个CP2K包
FARMING
在一次运行中运行N个独立作业
MC_ANALYSIS
运行(树)蒙特卡罗(Monte Carlo )轨迹文件分析
OPTIMIZE_BASIS
为给定的一组训练结构创建MOLOPT或ADMM基础的工具
OPTIMIZE_INPUT
优化CP2K输入中参数的工具
SWARM
运行基于群集的计算
TEST
做一些基准测试和测试
TMC
使用其他输入文件运行树蒙特卡罗算法
PROJECT_NAME
项目名(中间计算文件用这个命名) 别名:PROJECT
RUN_TYPE
仿真类型
BAND
Band方法 nudged elastic band (NEB)在已知反应物和产物之间寻找鞍点和最小能量路径的方法 http://blog.sina.com.cn/s/blog_6ac7d7930102uyw9.html http://theory.cm.utexas.edu/henkelman/research/saddle/neb/ http://theory.cm.utexas.edu/henkelman/pubs/sheppard08_134106.pdf
ENERGY
计算能量
GEO_OPT
几何优化 (键距离,角度等)
MD
分子动力学
NONE
不执行任何任务
VIBRATIONAL_ANALYSIS
振动分析
SEED
全局(伪)随机数生成器的初始种子,用于创建正态高斯分布随机数流。应正好为1或6个正整数值。复制单个值以用6个数字填充整个种子数组 默认值:2000
TRACE
是否应该写入程序执行的调试跟踪 默认值:否
FORCE_EVAL
计算力或能量
DFT
密度泛函理论参数 https://www.zhihu.com/question/30510921
BASIS_SET_FILE_NAME
basis set 文件名
CHARGE
系统的总电荷
LSD
同UKS
MGRID
多晶格信息
CUTOFF
最细网格截止,对SE或DFTB默认值为1
NGRIDS
要使用的多重网格数 默认值:4
REL_CUTOFF
确定高斯映射的网格,给出用于alpha=1的高斯的截止值。对于高精度的结果或具有可变单元格的模拟,可能需要50±10Ry的值。2.3之前的版本使用默认值30Ry。 默认值:4.00000000E+001 默认单位:[Ry]
RELATIVE_CUTOFF
同REL_CUTOFF
POISSON
设置泊松分解器。 同MM->POISSON
POTENTIAL_FILE_NAME
赝势文件名 赝势: 在近似中,将原子的核电子(即非价电子)以及原子核共同产生的一个复杂的势置换成一个“有效势”
DFT代码中可能的打印选项部分
E_DENSITY_CUBE
使用电子密度和(对于LSD计算)自旋密度控制立方体文件的打印
LOWDIN
控制Lowdin(自旋)总体分析的打印
MO_CUBES
控制分子轨道立方体的打印
NHOMO
如果printkey被激活,这个值控制作为HOMO(最高占用轨道)的数量,-1则特征值全部打印 默认值:1
NLUMO
如果printkey被激活,这个值控制作为NLUMO(最低未占用轨道)的数量,-1则特征值全部打印 默认值:0
MULLIKEN
控制Mulliken(旋转)总体分析的打印
QS
设置Quickstep框架所需的参数
CDFT
设置受约束DFT计算所需的参数。ATOM_GROUP部分的每次重复都定义了一个新的约束。约束在单独的外部SCF循环中聚合,设置从外部SCF部分读取。支持的约束:Becke和Gaussian Hirshfeld(部分) https://www.cp2k.org/howto:cdft 对 Kohn-Sham 能量方程 E K S E_\mathrm{KS} EKS加约束 E C D F T [ ρ , λ ⃗ ] = max λ ⃗ min ρ ( E K S [ ρ ] + ∑ c λ c [ ∑ i = ↑ , ↓ ∫ w c i ( r ) ρ i ( r ) d r − N c ] ) E_\mathrm{CDFT}[\rho, \vec{\lambda}] =\max_{\vec{\lambda}}\min_\rho \left( E_\mathrm{KS}[\rho] + \sum_c \lambda_c \left[\sum_{i = \uparrow, \downarrow} \int w_c^i(\mathbf{r})\rho^i(\mathbf{r})d\mathbf{r} - N_c \right] \right) ECDFT[ρ,λ ]=maxλ minρ(EKS[ρ]+∑cλc[∑i=↑,↓∫wci(r)ρi(r)dr−Nc]) 其中 λ ⃗ = [ λ 1 , λ 2 , ⋯ ] T \vec\lambda = [\lambda_1, \lambda_2, \cdots]^T λ =[λ1,λ2,⋯]T是(“strength(s) of the constraint potential(s)”), N c N_c Nc是target value of the constraint w i ( r ) w^i(\mathbf{r}) wi(r)各原子权重 系数COEFF
ATOM_GROUP
定义一组原子以用于CDFT约束。此部分的每次重复都会创建一个新约束。
ATOMS
指定包含在约束组中的原子列表
COEFF
定义原子列表中原子的系数 w i ( r ) w^i(\mathbf{r}) wi(r)。接受+/-1.0的值
CONSTRAINT_TYPE
确定要应用的约束类型。 默认值:CHARGE ####### ALPHA alpha自旋密度约束 ####### BETA beta自旋密度约束 ####### CHARGE 总电荷密度约束(rho_alpha + rho_beta)。 ####### MAGNETIZATION 磁化密度约束(rho_alpha - rho_beta)
ATOMIC_CHARGES
使用选定的权重函数计算原子CDFT电荷(Z=Z_core-Z_CDFT)。对于基于碎片的约束,电荷与碎片参考状态有关,即Z=Z_CDFT-Z_frag_reference。注意:如果原子数大于默认pw_pool max cache,则在解除分配原子网格期间,计算原子CDFT费用将提示无害警告。 默认值:FALSE
BECKE_CONSTRAINT
定义影响BECKE权重函数构造的设置
ADJUST_SIZE
用原子半径调整Becke细胞边界,以生成异核截止轮廓。这些半径是用关键字“原子半径”定义的 默认值:FALSE
ATOMIC_RADII
定义原子半径以生成异核截止剖面。按每个元素在输入坐标中出现的顺序为每个元素指定一个值 默认单位:[angstrom]
CAVITY_CONFINE
激活高斯腔约束。仅在型腔内部计算约束。空腔是由以约束原子为中心的球形高斯加和形成的 默认值:FALSE
CAVITY_SHAPE
指定使用的高斯腔的类型 默认值:DEFAULT ####### COVALENT 用共价半径(covalent radii )构造高斯函数 ####### DEFAULT 使用共价半径(单位为埃)来构造高斯,但对于半径大于此值的元素,固定为1.0_dp半径 ####### SINGLE 所有原子类型的单高斯分布,半径由CAVITY_RADIUS给定 ####### USER 使用用户定义的半径(关键字原子半径)构造高斯。 ####### VDW 用范德华半径构造高斯函数
CUTOFF_TYPE
指定在构建BECKE权重函数时使用的截断类型 默认值:GLOBAL ####### ELEMENT 对所有元素使用不同的值。从ELEMENT_CUTOFF读取的值 ####### GLOBAL 对所有元素使用单个值。从GLOBAL_CUTOFF读取
ELEMENT_CUTOFF
定义特定于元素的截止值,以确定每个实空间栅格点上哪些原子对权重函数有贡献。按元素在坐标中出现的顺序为每个元素指定一个值 默认单位:[angstrom]
EPS_CAVITY
创建空腔的密度阈值。忽略高斯密度低于阈值的网格点 默认值:1.00000000E-006
IN_MEMORY
在约束的初始形成过程中,预计算Becke约束产生的梯度,并将其存储在内存中。与限制、内存密集型结合使用。如果未计算力,则不会执行任何操作。 默认值:FALSE
SHOULD_SKIP
若网格点距离所有约束原子的GLOBAL_CUTOFF更远,则直接移动到下一个网格点,从而节省计算资源 默认值:FALSE
DUMMY_ATOMS
定义一组额外的原子,仅计算原子CDFT电荷。不能包含在节ATOM_GROUP中的任何约束原子
ATOMS
指定DUMMY_ATOMS组中包含的原子列表
OUTER_SCF
同SCF->OUTER_SCF
PROGRAM_RUN_INFO
同PRINT->PROGRAM_RUN_INFO
STRENGTH
约束力常数(拉格朗日乘数) λ ⃗ \vec\lambda λ 。为每个约束组指定一个值 默认值:0.00000000E+000
TARGET
约束目标值 N c N_c Nc。为每个约束组指定一个值。目标值是定义约束的原子上所需的价电子数、自旋矩或α或β电子数,如果采用两组原子之间的相对约束,则适当乘以原子系数。请注意,核心电荷不会从目标值中减去 默认值:0.00000000E+000
TYPE_OF_CONSTRAINT
指定所用约束的类型。 默认值:NONE
BECKE
BECKE约束。
HIRSHFELD
高斯Hirshfeld约束。部分实施:不使用力。没有经过广泛测试
NONE
没有约束
EPS_DEFAULT
尝试将所有EPS_xxx设置为可使能量校正达到EPS_DEFAULT的值 默认值:1.00000000E-010
EXTRAPOLATION
波函数的外推策略,例如MD。并非所有选项都适用于所有模拟方法。建议使用PS和ASPC,另请参见EXTRAPOLATION_ORDER。 默认值:ASPC
ASPC
始终稳定的预测-校正,类似于PS,但倾向于MD稳定性,而不是初始猜测精度(不适用于K点)
PS
EXTRAPOLATION_ORDER
PS或ASPC外推顺序(通常为2-4)。更高的阶数可能带来更高的精度,但对于大型系统来说,也需要付出一定的代价。在某些情况下,高阶外推不稳定,需要降低阶数 默认值:3
METHOD
指定应采用的电子结构方法 默认值:GPW
AM1
AM1半经验
DFTB
基于密度泛函的紧约束
GAPW
高斯和增广平面波方法
GAPW_XC
仅适用于XC的高斯和增广平面波方法
GPW
高斯和平面波方法
LRIGPW
恒等式方法的局部分解
MNDO
MNDO半经验
MNDOD
MNDO-d半经验
OFGPW
无轨道GPW方法
PDG
PDG半经验
PM3
PM3半经验
PM6
PM6半经验
PM6-FM
PM6-FM半经验
PNNL
PNNL半经验
RIGPW
HXC项的恒等式分解方法
RM1
RM1 半经验
XTB
GFN xTB扩展紧密绑定
SE
建立半经验方法所需的参数
COULOMB
SE计算中库仑项评估的设置参数
CUTOFF
计算库仑积分的原子截止半径截止。对于非周期计算,默认值正好是完整单元维度,以便评估所有对交互。相反,对于定期计算,使用默认数值。 默认值:1.20000000E+001 默认单位:[angstrom]
EXCHANGE
SE计算中交换项和核心哈密顿项评估的设置参数
CUTOFF
用于计算交换积分的原子截止半径截止。对于非周期计算,默认值正好是完整单元维度,以便评估所有对交互。相反,对于定期计算,默认值是单元格尺寸的1/4和输入中指定的值(明确定义或默认数值)之间的最小值。 默认值:1.20000000E+001 默认单位:[angstrom]
RESTART_FILE_NAME
同WFN_RESTART_FILE_NAME
SCF
自洽场运行参数 https://baike.baidu.com/item/%E9%87%8F%E5%AD%90%E5%8A%9B%E5%AD%A6%E7%9A%84%E8%87%AA%E6%B4%BD%E5%9C%BA%E8%BF%91%E4%BC%BC%E6%B3%95 自洽场方法(self-consistent field method, SCF),即自洽场近似法,是量子力学中迭代求解多粒子系统薛定谔方程的基本方法。其基本思想是首先按照某种方法给出波函数的一个估计,然后利用这个估计来计算电子密度,再通过电子密度来得到哈密顿量中与粒子间相互作用有关的项,再进行薛定谔