计算化学公社

标题: 编程计算小分子能量以及其他参数 [打印本页]

作者
Author: Markmahao 时间: 2022-4-12 12:28
标题: 编程计算小分子能量以及其他参数
各位老师好，最近突然萌生手撕量化计算的想法（主要想计算能量、偶极矩、极化率等），感觉这样更好理解计算的原理和流程。因为不是纯量化专业，水平比较差，请各位老师推荐几个教程和参考书，我可以先跟着做一做学一学。

作者
Author: niobium 时间: 2022-4-12 12:36
建议放弃

作者
Author: sobereva 时间: 2022-4-12 12:51
把能量、偶极矩、极化率的计算公式搞明白就够了。若自己从头写代码，包括电子积分、SCF迭代那些，完全是无用功

作者
Author: zjxitcc 时间: 2022-4-12 13:04
本帖最后由 zjxitcc 于 2022-4-12 13:07 编辑

想法很好，就是延期指日可待。如果你有耐心把Modern Quantum Chemistry前3章看完、做上部分习题，再考虑编程的问题。前3章就有HF能量、偶极矩的公式。

作者
Author: Markmahao 时间: 2022-4-12 13:32

sobereva 发表于 2022-4-12 12:51
把能量、偶极矩、极化率的计算公式搞明白就够了。若自己从头写代码，包括电子积分、SCF迭代那些，完全是无 ...

谢谢sob老师

作者
Author: Markmahao 时间: 2022-4-12 13:32

niobium 发表于 2022-4-12 12:36
建议放弃

哈哈哈

作者
Author: Markmahao 时间: 2022-4-12 13:37

zjxitcc 发表于 2022-4-12 13:04
想法很好，就是延期指日可待。如果你有耐心把Modern Quantum Chemistry前3章看完、做上部分习题，再考虑编 ...

哈哈哈，感谢！已经毕业了

单纯出于兴趣，后面遇到不懂的，还请帮忙解答

作者
Author: wzkchem5 时间: 2022-4-12 15:36
建议不要写从头算程序，而是写半经验程序。第一写起来快，第二你优化一下半经验参数、能提高偶极矩和极化率的精度以后，结合MD，用autocorrelation的方法算一些红外和拉曼光谱，可以发文章。尤其是拉曼，正常做从头算计算，拉曼活性需要用弥散基组，半经验方法虽然是极小基，但是有希望通过优化参数把弥散函数的贡献等效地考虑进来。

作者
Author: flyingchow 时间: 2022-4-12 21:07
简单写一个应该不算难，github上一搜一大把。因为好多美国高校的本科生计算化学课project就是自己写一段简单的SCF。
其实说是自己写，大多数都是在github上“借鉴”的。
举个最简单的例子，链接里面这个是只算1s轨道的SCF：
https://github.com/nickelandcopper/HartreeFockPythonProgram

作者
Author: Markmahao 时间: 2022-4-12 23:41

wzkchem5 发表于 2022-4-12 15:36
建议不要写从头算程序，而是写半经验程序。第一写起来快，第二你优化一下半经验参数、能提高偶极矩和极化率 ...

感谢您的建议！

作者
Author: Markmahao 时间: 2022-4-12 23:41

flyingchow 发表于 2022-4-12 21:07
简单写一个应该不算难，github上一搜一大把。因为好多美国高校的本科生计算化学课project就是自己写一段简 ...

谢谢您，我先去瞅瞅

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