小白请教问题，如何模拟非晶态的光激发效果？

luanshixiao

纯小白，刚知道量子化学是模拟非周期性结构，第一性原理是模拟周期性结构。
目前的课题是 模拟非晶态Ge-Sb-S玻璃受光照后的结构变化，但作为一个从没接触过计算化学的小白不知道该选择什么程序进行学习？

根据目前的思考，这个问题我觉得可以分成两步：

• 生成非晶态玻璃的结构模型
• 模拟光照
  

那么就产生了两个问题：

• 如何生成非晶态玻璃的结构模型？
  

玻璃是通过将物质升温至熔化之后再淬冷形成的，分子动力学一般只能模拟100ps的量级，相比实际淬冷速度快很多个数量级，所以不太能模拟玻璃；
还有一种方法是蒙特卡洛法，已有文献应用RMC来实现玻璃的形成，但似乎有点复杂。

• 如何模拟光照？
  

统计了以往的文献，主要是三种程序：Gaussian、VASP、WIEN2k。Gaussian的做法是激发至LUMO模拟光激发；WIEN2k的做法是将一个电子移除价带来模拟光激发。
但Gaussian属于量子化学，而WIEN2k属于第一性原理，所以哪一种更合适呢？按理说玻璃这种非周期性结构应该用Gaussian？

还请各位不吝赐教！

sobereva

“Gaussian的做法是激发至LUMO模拟光激发”这种说法是严重错误的，见
正确地认识分子的能隙(gap)、HOMO和LUMO
http://sobereva.com/543（http://bbs.keinsci.com/thread-16758-1-1.html）

你研究的这个问题绝对不是什么简单的东西。先去搜文献，如果有人算过类似问题的，去效仿；如果找不到类似文献，强烈建议以你现在的知识水平和经验不要去搞这个，否则很难腾出什么有意义的数据来。“光照”的考虑绝对不是一个简单的事情。做计算应当先从相对简单、研究得比较成熟的问题上手做。

luanshixiao

sobereva 发表于 2020-6-28 19:43
“Gaussian的做法是激发至LUMO模拟光激发”这种说法是严重错误的，见
正确地认识分子的能隙(gap)、HOMO和L ...

谢谢你的答复，类似的文献是有的，虽然目前的知识水平还看不懂。。
其次这是老板安排的死任务，没办法