计算化学公社

标题: 计算电离势IP时，高斯能否指定具体哪个位置失去一个电子？ [打印本页]

作者
Author: Jiangkunzhe 时间: 2025-6-26 20:50
标题: 计算电离势IP时，高斯能否指定具体哪个位置失去一个电子？
%mem=30GB
%nprocshared=22
%chk=MemAU+e-nbo.chk
# opt freq pop=nbo geom=connectivity def2tzvp em=gd3 m062x

Title Card Required

1 2
O                0.21582735 0.58273380 0.00000000
C                1.60607535 0.85614080 0.05981100
H                2.02821135 0.39979480 0.96879600
C                1.81679635 2.35840080 0.14674400
H                1.09825035 2.76544980 0.85847800
H                2.82563835 2.54109580 0.52500300
O                1.60670335 3.02502880 -1.07828300
H                2.37515335 2.86495480 -1.63582300
C                2.27715835 0.24340180 -1.15870000
H                1.80037935 0.66255880 -2.05248900
O                3.65122035 0.56314480 -1.12188400
H                4.09271035 -0.02342320 -1.74885000
C                2.05530635 -1.26102920 -1.16498000
H                2.52751835 -1.69172720 -0.27687000
O                2.68264235 -1.84657220 -2.28326500
H                2.07654035 -1.72664120 -3.02758900
C                0.55561635 -1.54431620 -1.15328400
H                0.37935235 -2.61478620 -1.00005900
O                0.08769035 -1.13660920 -2.41915500
H                -0.82399665 -1.43944620 -2.54642200
C                -0.07661665 -0.78989920 0.02595800
H                0.31622135 -1.23455620 0.95072600
O                -1.46365565 -0.93209720 0.12462100
C                -2.24057865 -0.30192420 -0.84916600
H                -1.64436965 0.46252180 -1.34487300
O                -2.57651965 -1.25320820 -1.92452700
C                -3.89213865 -1.62296320 -1.75409200
H                -4.56200865 -1.28118920 -2.53393900
C                -4.20371165 -1.84486620 -0.46427300
H                -3.35269065 -2.22659420 0.09420600
C                -4.71302665 -0.60587720 0.22384600
H                -4.78771365 -0.73527920 1.30117900
C                -3.48493465 0.35703780 -0.17256800
H                -3.06988065 0.85500380 0.70808100
C                -5.97074965 0.06941880 -0.35517300
H                -6.60764365 0.42996980 0.46073400
O                -6.69872465 -0.87803620 -1.12396600
H                -7.49133665 -0.45225620 -1.46389900
C                -5.45047265 1.21970880 -1.16936600
H                -6.09661365 1.86610180 -1.74951800
C                -4.13587265 1.36280880 -1.09140300
C                -3.32647465 2.39755580 -1.80329800
H                -2.71468465 1.90772880 -2.57235400
H                -2.63317065 2.86981180 -1.09440000
O                -4.19306865 3.34564780 -2.38728300
H                -3.67040465 3.96358280 -2.90419900
计算数据如下表格所示：
(, 下载次数 Times of downloads: 27)

请问各位大佬，请问高斯能否指定具体哪个位置失去一个电子呀？如果能，应该如何设置？非常感谢各位啦~

作者
Author: Stardust0831 时间: 2025-6-26 21:58
DFT计算一般不会指定某个原子的电荷，非要指定的话，试试看CDFT（约束性dft）

作者
Author: sobereva 时间: 2025-6-27 04:20
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “计算电离势IP” 改了，以后务必注意！

作者
Author: sobereva 时间: 2025-6-27 04:23
算IP时指定这个通常没有任何意义，电离过程中实际电荷分布该怎么变化就怎么变化。如果非要用下文的方法施加约束使得电离掉的一个电子完全来自于某个局部片段（通常很难解释得通此做法的意义、一般会被质疑引入莫名其妙的人为虚假因素），由于阳离子态的能量会变得更高，因此会得到偏大的电离能。

作者
Author: Jiangkunzhe 时间: 2025-6-27 15:02

sobereva 发表于 2025-6-27 04:20
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此 ...

sob老师，我已重新阅读新社员必读帖、公告栏内容，类似错误不会再犯，浪费您宝贵时间实属抱歉！

作者
Author: Jiangkunzhe 时间: 2025-6-27 15:02

sobereva 发表于 2025-6-27 04:23
算IP时指定这个通常没有任何意义，电离过程中实际电荷分布该怎么变化就怎么变化。如果非要用下文的方法施加 ...

我明白啦，感谢老师的答复~

欢迎光临计算化学公社 (http://bbs.keinsci.com/)