如何把实验中的恒电位转化为Gaussian里的电场

www1

各位老师同学们好，我是还没入门的小白，请教各位问题，希望获得帮助。
这个课题需要计算团簇在一个恒电位下的一系列情况。之前查到是用field=read这样的加外场的指令来实现。但是目前已知的信息只有-500mv的恒电位，这样的如何量化为电场（a.u.）的大小呢？
用简单的公式去进行计算不太现实，团簇会直接接触到电极表面，感觉如果直接按照电极之间是匀强电场来处理的话误差会很大。

这个问题可能对于有经验的老师或同学比较简单，但是我确实不知道该如何处理。

谢谢

wzkchem5

你这个计算应该是电解过程，和外加电场计算完全是不一样的。
高斯里field=read指的是在绝缘介质里，加一个电场，让你的分子被这个电场极化，但是不得失电子；而你的体系应该是允许分子得失电子的，与此同时这个电位导致的电场非常小，几乎没有化学效应，所以实验上观察到的外加电位的效应应该都是因为得失电子导致的。这种情况下应该这样做：（1）用https://pubs.acs.org/doi/10.1021/acs.jctc.6b00252的方法计算这个体系的各级氧化还原电势；（2）比较计算的氧化还原电势和实验电位，判断分子在实验条件下所处的氧化态；（3）按照那个氧化态做常规计算就行了，不需要加电场。

www1

wzkchem5 发表于 2022-6-21 21:28
你这个计算应该是电解过程，和外加电场计算完全是不一样的。
高斯里field=read指的是在绝缘介质里，加一个 ...

谢谢老师，这个团簇是作为催化剂使用的，确实是存在这得失电子的情况，甚至配体有可能脱落。
这是一个比较大的Au团簇，有450个原子，请问这样的体系是不是不太适合用Gaussian呢，谢谢。

wzkchem5

www1 发表于 2022-6-22 02:46
谢谢老师，这个团簇是作为催化剂使用的，确实是存在这得失电子的情况，甚至配体有可能脱落。
这是一个比 ...

这个太大了，不仅高斯很难算动，其他软件也都算不太动，除非用半经验方法。建议把体系截得小一些再算，或者用ONIOM算，内层用DFT，外层用xTB。

www1

wzkchem5 发表于 2022-6-22 14:44
这个太大了，不仅高斯很难算动，其他软件也都算不太动，除非用半经验方法。建议把体系截得小一些再算，或 ...

好的，谢谢老师。