| 查看: 1203 | 回复: 3 | ||
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的晶格常数乘积为一常数，这样就不需要利用脚本进行优化了。 这种办法自然保证了无论从那个方向算应变都是唯一的。不过，据说不是很合理。 有什么建议？谢谢 |
回复此楼» 猜你喜欢
346求调剂[0856]
已经有7人回复
-
一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！
已经有3人回复
-
310求调剂
已经有3人回复
-
机械专硕299求调剂至材料
已经有4人回复
-
070300化学319求调剂
已经有7人回复
-
08工科 320总分 求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有12人回复
-
307求调剂
已经有3人回复
-
265求调剂
已经有3人回复
-
301求调剂
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有应变就一定有应力吗？
已经有16人回复
- 关于BCC-Fe的弹性常数的能量-应变计算
已经有47人回复
- 电阻-应力（应变）曲线拟合
已经有5人回复
- 压缩应力应变曲线问题
已经有14人回复
- vasp计算应力--应变曲线时应力值发生突变的问题
已经有8人回复
- 拉伸试件不同，计算应力应变的方法相同吗
已经有12人回复
- 拉伸应变的计算
已经有4人回复
- 请教表面吸附DOS计算中Ef很大于0后仍然不回0
已经有7人回复
- 求助VASP计算弹性常数的问题
已经有5人回复
- 请教如何计算应力-应变曲线下的面积
已经有19人回复
- 对称分子在ATK软件模拟计算下I-V曲线不对称是什么原因？可以提供文献出处最好
已经有5人回复
- 镁合金压缩应力应变曲线
已经有17人回复
- 【求助】如何计算应变对应的应力值
已经有5人回复
- 【求助】CASTEP计算弹性常数得到应力应变曲线
已经有17人回复
- 【求助】算VASP弹性常数时的能量-应变出错
已经有18人回复
- 【讨论】单向拉伸过程中如何计算应变速率呢
已经有13人回复
- 【求助】castep计算应力应变曲线结果文件中的两个stress tensor
已经有5人回复
- 【求助】MS计算应力应变曲线
已经有31人回复
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
2楼2014-11-01 00:44:30
qinform
木虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 17458.5
- 红花: 4
- 帖子: 1482
- 在线: 1157.8小时
- 虫号: 97142
- 注册: 2005-11-08
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
是的，差别应该是体现了你说的应变不对等。 你认为“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
4楼2017-10-12 11:30:16
