简单量子化学问题答疑专帖（旧版，已关帖）

明园 · 发表于 Post on 2019-12-16 15:24:22

社长，这个中间体整体不带电荷，gjf文件中的电荷和自选多重度，就直接按默认的就可以吗？

elvisng · 发表于 Post on 2019-12-16 15:27:14

如果想用 ORCA 做单点, 但想用 unrestricted，请问该如何修改以下的关键词？

例如：
1. ! DLPNO-CCSD(T) normalPNO RIJK cc-pVTZ cc-pVTZ/JK cc-pVTZ/C
2. ! wB97M-V def2-TZVP def2/J RIJCOSX grid4 gridx4
3. ! PWPB95 D3 def2-QZVPP def2/J def2-QZVPP/C RIJCOSX grid4 gridx4

谢谢

lianghuashou · 发表于 Post on 2019-12-16 17:32:51

Sob老师，我在计算电子振动的发射光谱，首先报错信息是结构变化大，因此我的输入改为
#P B3LYP/genecp  opt freq=(FC,emission,ReadFCHT) scf(nosymm)

Title Card Required

0 1
   6       -3.270488 4.135032 0.291606
   6       -1.887551 4.363492 0.339629
   6       -1.022957 3.262388 0.238430
。。。。。

basis set

ForceFCCalc

c:/XXX/pt_t1_n2_ph_01.chk

但是依然报错，报错信息是 Initial State
-------------
Data taken from checkpoint file "ForceFCCalc"
Error termination in NtrErr:
ntran open failure returned to fopen.

请问sob老师，这个怎么解决，谢谢

ssfm · 发表于 Post on 2019-12-16 18:05:46

sobereva 发表于 2019-12-13 19:27
看最新版Multiwfn手册4.A.13节的第9部分

PS：本帖只回复量化问题，Multiwfn相关问题请在本论坛的波函 ...

谢谢老师的耐心回复

zj123_45 · 发表于 Post on 2019-12-16 19:53:18

sobereva 发表于 2019-12-13 19:26
解决SCF不收敛问题的方法
http://sobereva.com/61

好的老师，问题已解决，感谢老师

675543414 · 发表于 Post on 2019-12-16 20:18:33

本帖最后由 675543414 于 2019-12-16 21:39 编辑

sob老师您好，我有一个体系，反应物是两个自由基分子（0 2,0 2），产物为闭壳层（0 1）。
请问我在寻找过渡态的时候，只需要改成ub3lyp 嘛？电荷和自旋多重态用（0 1）？

lastzealot · 发表于 Post on 2019-12-16 22:11:27

本帖最后由 lastzealot 于 2019-12-16 22:17 编辑

Sob老师您好，在看文献时看到了元素铀的SO pseudopotential（如图所示），我上网查了一下只看到ECP60MWB等一系列标量相对论的基组和pp，没有找到所谓的spin orbit赝势。
请问老师，这种赝势应该去哪里下载？

（molpro manual中I的SOspin orbit赝势如下所示）
! Iodine-ECP (Dirac-Fock) with SO-coupling
!
ecp,I,46,4,3;
1; 2, 1.00000000,    0.00000000;                                     ! lokal term = 0
2; 2, 3.50642001,    83.09814545; 2, 1.74736492,    5.06370919; ! s-terme
4; 2, 2.99860773, 1/3* 81.88444526; 2, 3.01690894, 2/3* 83.41280402; ! p-terms with weights
2, 1.59415934, 1/3*  2.32392477; 2, 1.19802939, 2/3*  2.72079843;
4; 2, 1.03813792, 2/5*  6.40131754; 2, 1.01158599, 3/5*  6.21328827; ! d-terms with weights
2, 2.04193864, 2/5* 19.11604172; 2, 1.99631017, 3/5* 19.08465909;
4; 2, 2.64971585,-3/7* 24.79106489; 2, 2.75335574,-4/7* 24.98147319; ! f-terms with weights
2, 0.49970082,-3/7*  0.27936581; 2, 0.79638982,-4/7*  0.70184261;
4; 2, 2.99860773,-2/3* 81.88444526; 2, 3.01690894, 2/3* 83.41280402; ! ECP-SO for p-terms
2, 1.59415934,-2/3*  2.32392477; 2, 1.19802939, 2/3*  2.72079843;
4; 2, 1.03813792,-2/5*  6.40131754; 2, 1.01158599, 2/5*  6.21328827; ! ECP-SO for d-terms
2, 2.04193864,-2/5* 19.11604172; 2, 1.99631017, 2/5* 19.08465909;
4; 2, 2.64971585, 2/7* 24.79106489; 2, 2.75335574,-2/7* 24.98147319; ! ECP-SO for f-terms
2, 0.49970082, 2/7*  0.27936581; 2, 0.79638982,-2/7*  0.70184261;
!
! Iodine-basis
!
s,I,0.2027624,0.4080619,0.8212297,1.6527350,3.3261500;
c,1.5,-0.4782372,-0.5811680,0.2617769,0.4444120,-0.1596560;
s,I,0.05,0.1007509;
p,I,0.2027624,0.4080619,0.8212297,1.6527350,3.3261500;
c,1.5,0.4251859,0.2995618,0.0303167,-0.2064228,0.0450858;
p,I,0.05,0.1007509,0.01;  ! diffuse p-Funktion wegen evt. neg. Part.Ldg
d,I,0.2,0.4;
f,I,0.3;

谢谢老师

lastzealot · 发表于 Post on 2019-12-16 22:15:13

didi_dudu 发表于 2015-3-18 15:43
请问一下常用ECP中如ECPXXMDF的MDF是什么意思？是multiconfiguration Dirac±Hartree±
Fock (MCDHF) 的进 ...

Library keywords in case b) are of the form ECPnXY ; n is the number of core electrons which are replaced by the pseudopotential, X denotes the reference system used for generating the pseudopotential (X = S: states of the single-valence-electron ion; X = M: states of the neutral atom and of ions with low charge), and Y stands for the theoretical level of the reference data (Y = HF: Hartree-Fock, Y = WB: quasi-relativistic; Y = DF: relativistic). For one- or two-valence electron atoms X = S, Y = DF is a good choice, while otherwise X = M, Y = WB (or Y = DF) is recommended. (For light atoms, or for the discussion of relativistic effects, the corresponding Y = HF pseudopotentials may be useful.) Additionally, spin-orbit (SO) potentials and core-polarization potentials (CPP) are available, to be used in connection with case b)
ECPs, but these are not currently contained in the library, so explicit input is necessary here (cf. below).

sobereva · 发表于 Post on 2019-12-16 23:48:25

lastzealot 发表于 2019-12-16 22:11
Sob老师您好，在看文献时看到了元素铀的SO pseudopotential（如图所示），我上网查了一下只看到ECP60MWB等 ...

Stuttgart赝势基组官网上，选择molpro之类的支持旋轨势的程序，在定义页面上就能看到旋轨势

参考此文
用NWChem做SODFT在DFT计算中考虑旋轨耦合效应
http://sobereva.com/368（http://bbs.keinsci.com/thread-5573-1-1.html）

sobereva · 发表于 Post on 2019-12-16 23:50:01

675543414 发表于 2019-12-16 20:18
sob老师您好，我有一个体系，反应物是两个自由基分子（0 2,0 2），产物为闭壳层（0 1）。
请问我在寻找过 ...

这种情况，通常过渡态是闭壳层的。就算真是开壳层，光写个U什么用也没有，该怎么考虑看
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

你先当闭壳层找过渡态，如果找到了，做个波函数稳定性测试，如果稳定，说明就是闭壳层状态

sobereva · 发表于 Post on 2019-12-16 23:51:14

lianghuashou 发表于 2019-12-16 17:32
Sob老师，我在计算电子振动的发射光谱，首先报错信息是结构变化大，因此我的输入改为
#P B3LYP/genecp op ...

明显是把ForceFCCalc当文件名读取了。位置调换一下

sobereva · 发表于 Post on 2019-12-16 23:53:07

elvisng 发表于 2019-12-16 15:27
如果想用 ORCA 做单点, 但想用 unrestricted，请问该如何修改以下的关键词？

例如：

DFT的加上UKS，后HF的加上UHF

sobereva · 发表于 Post on 2019-12-16 23:54:19

shalene 发表于 2019-12-16 13:25
The authors should clarify in which cases an implicit solvation scheme was used. 这话是啥意思？用了g ...

说明白哪些任务用了溶剂模型，用了什么溶剂模型

lastzealot · 发表于 Post on 2019-12-16 23:59:48

sobereva 发表于 2019-12-16 23:48
Stuttgart赝势基组官网上，选择molpro之类的支持旋轨势的程序，在定义页面上就能看到旋轨势

参考此文
...

感谢Sob老师帮助

sobereva · 发表于 Post on 2019-12-17 00:02:59

明园发表于 2019-12-16 15:24
社长，这个中间体整体不带电荷，gjf文件中的电荷和自选多重度，就直接按默认的就可以吗？

电荷用0，不同自旋多重度都算算，哪个能量低用哪个

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 简单量子化学问题答疑专帖（旧版，已关帖）

求助，中间体整体不带电荷的结构。

浏览过的版块