应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 应变计算不对称
51/1返回列表
查看: 1173  |  回复: 3
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

qinform

木虫 (著名写手)

[求助] 应变计算不对称已有1人参与

对一个二维体系进行DFT计算，a,b,c正交，但是两两不相等。c 包含了真空层20 Angstrom, 所以在施加应变的计算时忽略z方向的变化。
如果在x方向上施加应力，即改变x方向的晶格常数，比如+1%，让y方向晶格进行弛豫（通过脚本实现），可以发现y方向晶格常数有应变，比如 -0.9%.
然而，逆推回去，也就是从y方向 -0.9% 应变出发， 让x方向晶格自由弛豫，同样用脚本实现，会发现：此时优化得到的x方向应变不再是 +1%，优化值大概在+0.1%。
（1）请问，这种差别是什么原因？
（2）感觉需要消除这种差异，不知道怎么消除？
（3）我尝试了在添加应变的过程中保持体积不变，这样，由于z 方向晶格不变，相当于x和y的晶格常数乘积为一常数，这样就不需要利用脚本进行优化了。
这种办法自然保证了无论从那个方向算应变都是唯一的。不过，据说不是很合理。
      有什么建议？谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-10-31 23:57:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Amersly

新虫 (小有名气)
请问楼主如何通过脚本实现只优化某一个方向的晶格？
赞一下 回复此楼
4楼2017-10-12 11:30:16
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

stractor

金虫 (著名写手)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
qinform: 金币+3, ★★★很有帮助 2014-11-01 03:45:29
差别来源于应变的不对等。你在x方向施加应变，让y方向驰豫，也就是相当于x方向不平衡，而让y方向平衡。反之，y方向不平衡情况下而让x方向平衡。这两种情况是不同的，所以不对等。
但是，如果你的材料x,y方向是对称的，在x方向或y方向施加同样大小的应变，另一方向驰豫大小总是相同的。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
2楼2014-11-01 00:44:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qinform

木虫 (著名写手)
引用回帖:
2楼: Originally posted by stractor at 2014-11-01 00:44:30
差别来源于应变的不对等。你在x方向施加应变，让y方向驰豫，也就是相当于x方向不平衡，而让y方向平衡。反之，y方向不平衡情况下而让x方向平衡。这两种情况是不同的，所以不对等。
但是，如果你的材料x,y方向是对称 ...

是的，差别应该是体现了你说的应变不对等。
你认为“x方向不平衡，而让y方向平衡” 是不是 exx不等于0， 而eyy=0？
其中 exx为x 方向的应变（ax-ax0)/ax0, eyy 为y方向的应变。
我的目的是拟合gamma点能级随应变的变化，
考虑到一阶项应该有 a * exx +b * eyy， 其中a,b 为待定参数。
不知道“在x方向施加应变，让y方向驰豫”这样一个过程该是怎么拟合？
（1）认为exx,eyy都不为零，显然由于第一个帖子中所描述两种情形的不对等，拟合参数会显著不同；
（2）在x方向施加应变，让y方向驰豫，认为exx不为零，eyy=0，这样拟合可以得到一个参数a,
在y方向施加应变，让x方向驰豫，认为eyy不为零，exx=0，这样拟合可以得到一个参数b,
这时候也存在一个问题，显然一个方向施加应变，比如x方向，另外一个方向y 肯定会偏离平衡时的晶格常数，把y方向的应变取为0 似乎也不合适。
不知道问题表达清楚了没有？ 谢谢
赞一下 回复此楼
3楼2014-11-01 04:03:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭