请教关于弯曲石墨烯模型如何构建的问题

白皛白

             老师们好，最近想要构建一个弯曲的石墨烯类似的模型，也看到一篇做了相关结构的文献，在这篇文献里是通过更改晶格常数来诱导平面的石墨烯出现弯曲的，大概意思是将晶格常数减小至平坦结构的95%或90%来使石墨烯出现不同程度的弯曲这样……我也想要效仿这种改变晶格常数的方法，先ISIF=3将晶格完全弛豫，然后再ISIF=2，ab两个方向上的晶格常数全部缩小至原本的90%，再做结构优化。但是却发现石墨烯并没有发生弯曲之类的形变……
             所以想要请教下各位老师，直接缩小晶格常数是否真的能让石墨烯发生弯曲呢？如果这个方法行不通的话应该怎样去构建不同弯曲程度的石墨烯呢？
             感谢老师们阅读我的提问，还望您能不吝赐教。

superbaba

是否可以先做一个GO的结构优化 然后把含氧基团删了

白皛白

superbaba 发表于 2023-3-10 13:03
是否可以先做一个GO的结构优化 然后把含氧基团删了

非常感谢您的回复，今天尝试去调整了下初始结构，然后再做结构优化就优化出弯曲的结构了。
原本的结构没有产生弯曲的原因似乎是因为初始结构在Z轴上太过于规整了，今天试着稍微调整了下原子的位置，把一部分原子抬高了0.1A~0.2A左右之后，再缩小晶格常数优化就会出现和文献里类似的弯曲了。
我对这个现象的粗浅的理解是，或许缩小晶格常数就相当于从ab两个方向对二维材料施加压力，对于太过于平坦的结构来说只在平面的两个方向上施加压力就比较难产生弯曲。但是如果结构原本就有些弯曲的话，再施加压力就会容易产生弯曲了，尽管最初的弯曲是非常非常小的，大概就像是沿着折痕折纸一样。不过这只是我的猜测，有理解错误的地方还请各位老师指正。

superbaba

白皛白 发表于 2023-3-10 20:39
非常感谢您的回复，今天尝试去调整了下初始结构，然后再做结构优化就优化出弯曲的结构了。
原本的结构没 ...

看了一下CONTCAR 原子确实几乎都在一个平面上😂