|
|
★ ★ ★ lvzhu2007(金币+3,VIP+0):3q!!!
256个电子,只有库仑力相互作用,在给定的温度下,求结构因子
你的库仑力怎么处理的?是否用周期性边界条件,还有你的MC每一步怎么取得。我处理过LJ势情况下的蒙特卡罗,60个粒子的,在保证每一步接受概率都为0.25~0.5的情况下,MCs的步数我在平衡态之前用了上亿步。
不知道你模拟的是什么体系,感觉只考虑电子,又要求结构因子的话,像是做X衍射之类的精修吧.你可以看看,算出来的结构因子合理否,这块就不熟悉了.
MC的达到平衡就是看你处理的量,是否涨落较小了. |
|