简单Multiwfn使用问题和波函数分析问题答疑专贴

sobereva

snljty 发表于 2020-4-30 17:07
请问老师，Multiwfn输出各种图里面文字的字体是什么？是dislin库决定的么，这个编程时候能改么？谢谢。

是的
dislin库支持的字体很有限，Multiwfn用的已经是我尝试过的最好的了。没法改

snljty

sobereva 发表于 2020-5-1 02:59
是的
dislin库支持的字体很有限，Multiwfn用的已经是我尝试过的最好的了。没法改

谢谢老师~

函数与激情

sob老师，我看Multiwfn 手册4.A.9部分最后处理激子耦合如下所述：
(5) Load the dimer .chg file into Multiwfn, enter main function 7 and select option -2, then input atom list of the two monomers in turn, the excitonic coupling energy will be printed immedately.
可是我将两个单体chg文件以及总的chg文件用Multiwnf打开都不行啊，会闪退。下面是我其中两个单体的chg文件，麻烦老师有空的时候看一下。

1. C    -7.259304    1.429563    4.429785-164.8392771359
   
2. C    -7.673013    2.422267    3.545682 -12.1407241425
   
3. O    -7.317362    1.662322    5.763943 101.2490187813
   
4. C    -6.786764    0.220981    3.899750 170.6254460528
   
5. C    -7.590898    2.187221    2.176190 -21.9269769743
   
6. H    -8.050678    3.359025    3.939351  28.9904426781
   
7. H    -7.016531    0.841431    6.210532 -38.1070779712
   
8. C    -6.708365   -0.001074    2.526141  54.2790870187
   
9. C    -6.309024   -0.942128    4.640448-162.8329425439
   
10. H    -7.912167    2.971434    1.496650  43.8632720243
    
11. C    -7.100940    0.993620    1.643793 -93.7223650431
    
12. C    -6.179103   -1.395328    2.240383-271.6340381477
    
13. C    -5.893594   -1.999809    3.633072 -71.2182803083
    
14. O    -6.243346   -1.028800    5.861082  80.2227052749
    
15. H    -7.037168    0.854390    0.570011  47.2277687395
    
16. H    -5.281213   -1.356048    1.616415  93.6207674926
    
17. H    -6.920555   -1.980130    1.686780  92.7937250323
    
18. H    -4.840739   -2.247693    3.786284  61.1114332046
    
19. H    -6.462400   -2.915279    3.823633  62.4380159678
    

复制代码

1. C    -2.985891    0.001073    4.458905  58.9466905675
   
2. C    -2.907493   -0.220981    3.085296 161.4882072807
   
3. C    -2.593317   -0.993620    5.341253 -95.5386202835
   
4. C    -3.515153    1.395328    4.744663-277.0532018255
   
5. C    -2.434953   -1.429563    2.555261-165.3475721891
   
6. C    -3.385232    0.942128    2.344598-154.5911183277
   
7. C    -2.103359   -2.187221    4.808856 -19.5761752543
   
8. H    -2.657089   -0.854390    6.415034  46.6479987709
   
9. C    -3.800662    1.999809    3.351973 -68.0534460399
   
10. H    -4.413043    1.356048    5.368630  94.4067117927
    
11. H    -2.773702    1.980129    5.298266  94.6418953181
    
12. O    -2.376895   -1.662323    1.221103 101.8633189791
    
13. C    -2.021244   -2.422267    3.439364  -9.5358154421
    
14. O    -3.450910    1.028800    1.123964  78.6112101365
    
15. H    -1.782090   -2.971434    5.488395  43.9430086234
    
16. H    -4.853517    2.247693    3.198762  60.2597302499
    
17. H    -3.231857    2.915279    3.161413  59.7733907183
    
18. H    -2.677726   -0.841431    0.774514 -37.7270529954
    
19. H    -1.643579   -3.359026    3.045694  26.8408399204
    

复制代码

snljty

函数与激情 发表于 2020-5-6 21:12
sob老师，我看Multiwfn 手册4.A.9部分最后处理激子耦合如下所述：
(5) Load the dimer .chg file into Mul ...

不知道你这个chg是什么软件的chg格式，Multiwfn指的是内置的格式，最后一列是原子电荷，不可能每个原子带一百多个电荷。

函数与激情

snljty 发表于 2020-5-6 22:18
不知道你这个chg是什么软件的chg格式，Multiwfn指的是内置的格式，最后一列是原子电荷，不可能每个原子带 ...

我是用cubegen产生的 也许我弄错了 我检查下

函数与激情

sob老师，之前的操作有误，我重新弄了下现在chg文件里的charge也正常了，为什么算出来的耦合能是负值呢 这有着什么含义呢

1. Input atom list for fragment 1, e.g. 1,4,6-9
   
2. Input 0 can exit
   
3. 1-19
   
4. Input atom list for fragment 2, e.g. 1,4,6-9
   
5. 20-38
   
6. Electrostatic interaction energy:       -0.59 kJ/mol       -0.14 kcal/mol
   

复制代码

刘小英雄steven

老师好，关于静电势ESP和电子定域化函数ELF的问题需要请教您。我发现二维碳材料的缺陷位置ESP值较小，这样似乎说明该区域静电势由电子主导，且更加容易吸引亲电气体CO2（电荷分析发现CO2吸附后呈负电荷）。但是某篇Advanced Materials 发表的文章提到（见附图）：文章ELF分析提到："缺陷区域ELF值小，说明缺陷位置较周围而言电子密度较小,以至于缺陷附近的碳原子“patically positive charge” "。这与ESP分析结论是否矛盾？而且这样似乎不利于向亲电的CO2传递电子。我想，可否在文章这么解释：缺陷区域电子束缚力小，虽然电子密度小，但是向气体输出电子的能力高（相比无缺陷位置）？

lianghuashou

sob老师，关于电子-空穴的分析，我一直有个疑问。对于电子激发到高能级轨道的过程来说，通常是电子的来处是空穴， 电子的去处是电子。但是对于去激发的过程来说，电子从高能级跃迁回到低能级，这时候怎么区别电子和空穴呢。感谢sob老师。

sobereva

lianghuashou 发表于 2020-5-13 11:44
sob老师，关于电子-空穴的分析，我一直有个疑问。对于电子激发到高能级轨道的过程来说，通常是电子的来处是 ...

“去激发”指代不明。如果你是指发射过程，把空穴和电子反过来解释就行了。如果你说的去激发是指去激发组态函数，这在空穴-电子分析中已经自然而然考虑了，用户不需要顾虑。

sobereva

刘小英雄steven 发表于 2020-5-13 09:07
老师好，关于静电势ESP和电子定域化函数ELF的问题需要请教您。我发现二维碳材料的缺陷位置ESP值较小，这样 ...

ELF和电子密度没有任何直接相关性，完全没有正比关系

那个缺陷的特征应当用通过静电势来考察，算出来计算为明显负值说明带正电的物质倾向于插入这里

另外，要搞清楚静电势的表达式。某处的静电势不是看此处的电子密度的大小，而是由附近一片区域内的电子密度和原子核电荷共同决定的。只有搞清楚物理意义、计算表达式，才能正确进行分析。

13163031955

计算单点能时出现l502错误，请问该怎么解决的，谢谢大家。以下是输入文件
%mem=8000MW
%nprocshared=36
# mp2/genecp geom=connectivity sp

Title Card Required

0 2
Cl                 4.51486600    1.41871400    0.00107500
Cl                -4.51469400   -1.41898600    0.00106800
Cl                -1.44995000    1.97217700   -0.00026900
Cl                 1.44963200   -1.97200800   -0.00026300
Au                -2.99717900    0.25458000    0.00006100
Ag                 0.00004400    0.00024500   -0.00078300
Au                 2.99718500   -0.25470400    0.00005800

snljty

《Gaussian 常见报错及解决方法；新手求助报错时的注意事项》
http://bbs.keinsci.com/thread-4829-1-1.html
顺便，这种体系没什么道理用MP2

lianghuashou

谢谢sob老师的回复，对于发射过程来说，比如电子的跃迁是 从119 MO回到118 MO， 所以对发射过程来说， 119 MO 是电子，118MO是空穴，这样是对的吗。谢谢。

sobereva

13163031955 发表于 2020-5-14 17:15
计算单点能时出现l502错误，请问该怎么解决的，谢谢大家。以下是输入文件
%mem=8000MW
%nprocshared=36

合并帖子合并错了...本来是要合并到http://bbs.keinsci.com/thread-806-1-1.html

这种老生常谈的问题发帖前先搜论坛
仔细看
解决SCF不收敛问题的方法
http://sobereva.com/61

用MP2会被内行审稿人批得很惨，仔细看
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272（http://bbs.keinsci.com/thread-536-1-1.html）

sobereva

lianghuashou 发表于 2020-5-15 10:42
谢谢sob老师的回复，对于发射过程来说，比如电子的跃迁是 从119 MO回到118 MO， 所以对发射过程来说， 119  ...

显然119是空穴，118是电子。激发态时候119是占据的