24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 724  |  回复: 8
当前主题已经存档。

doudouding

木虫 (小有名气)

[交流] 【求助】pwscf可以计算弹性常数吗

需要计算弹性常数,在manual里找不到弹性常数的计算方法,是不是pwscf不可以计算弹性常数
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zxzj05(金币+1,VIP+0):3Q! 鼓励讨论! 11-4 20:57
yjmaxpayne(金币+1,VIP+0):多谢分享! 11-7 11:52
方法一:计算不同体积下的体系能量
然后利用ev.x进行拟合,可以拟合出来的体弹性模量
http://blog.sina.com.cn/s/blog_5f15ead20100f5jx.html

方法二:计算声子谱,得到力常数
然后自己写程序得到弹性模量(具体方法参考科恩和黄昆的名著)
2楼2009-11-04 12:15:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

doudouding

木虫 (小有名气)


fegg7502(金币+1,VIP+0):thank you very much! 11-8 03:46
非常感谢,用ev.x拟合得到的只是体弹性模量,如果是想计算得到Cij,有什么办法吗?
3楼2009-11-04 13:52:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zxzj05(金币+1,VIP+0):3Q! 鼓励讨论! 11-4 20:57
第二种方法好像可以得到Cij
只是可能,具体过程我忘记了:)
4楼2009-11-04 20:13:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wuchenwf(金币+1,VIP+0):谢谢回帖交流 11-5 23:26
yjmaxpayne(金币+1,VIP+0):多谢分享 11-7 11:51
yjmaxpayne(金币+0,VIP+0):能否将理论方面的一些细节过程给大伙儿稍微解释一下? 毕竟力常数是一个体系的重要指标。 11-7 11:52
理论上,什么软件都可以算弹性常数,
真理是一点点接近的!
5楼2009-11-05 21:51:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watermall

金虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjmaxpayne(金币+1,VIP+0):多谢分享 11-7 11:51
可以计算,但是因为pwscf里本身没有这个功能,但是可以通过分步计算,就能得到所需要的Cij的。我曾用自己的方法计算六角结构的弹性系数的
用声子谱计算也可以,但可能工作量比较大
用自己的方法也可以,但可能比较烦
6楼2009-11-07 11:26:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
fegg7502(金币+1,VIP+0):thank you very much! 11-8 03:46
引用回帖:
Originally posted by watermall at 2009-11-7 11:26:
可以计算,但是因为pwscf里本身没有这个功能,但是可以通过分步计算,就能得到所需要的Cij的。我曾用自己的方法计算六角结构的弹性系数的
用声子谱计算也可以,但可能工作量比较大
用自己的方法也可以,但可能比 ...

您好;)
如何通过分步计算,就能得到所需要的Cij的?
能说的详细些么?
谢谢;P
7楼2009-11-07 16:25:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watermall

金虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
fegg7502(金币+1,VIP+0):thank you very much! 11-8 03:47
引用回帖:
Originally posted by xirainbow at 2009-11-7 16:25:


您好;)
如何通过分步计算,就能得到所需要的Cij的?
能说的详细些么?
谢谢;P

1. 优化好晶体结构;
2. 对优化好的晶体结构进行形变,再relax下原子位置,得到应力;
通过应力应变关系,得到Cij的关系;
3. 当然,第二步要做多次,这取决于晶体的结构对称性;
注:
1. control卡里有 tstress 这个关键词,是计算应力的
8楼2009-11-07 19:09:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
fegg7502(金币+1,VIP+0):thank you very much! 11-8 03:47
引用回帖:
Originally posted by watermall at 2009-11-7 19:09:

1. 优化好晶体结构;
2. 对优化好的晶体结构进行形变,再relax下原子位置,得到应力;
通过应力应变关系,得到Cij的关系;
3. 当然,第二步要做多次,这取决于晶体的结构对称性;
注:
1. control卡里有 t ...

哦;)
原来您是从Cij最原始的定义出发
进行数值计算的:)
非常感谢您的指点;P
9楼2009-11-07 22:53:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 doudouding 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见