24小时热门版块排行榜    

查看: 2061  |  回复: 8

meizhixin

新虫 (小有名气)

[交流] 能量应变法求弹性常数问题

想用能量应变法求弹性常数,其中由defvector.f得到defvector.x这一步,介绍中说是使用g77 -o defector.x defector.f命令得到,为什么我的总不成功,得不到defvctor.x文件,有没有计算成功的,求指教,谢谢

发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
meizhixin: 金币+10, 谢谢回复,很有帮助 2017-03-06 10:09:59
你有该defector.f的源代码吗?
如果有的话,在你的linux 终端下,运行命令 gfortran   -o defector.x defector.f也是可以的。
或者运行gfortran defector.f命令得到a.out,然后将a.out改名为defector.x  (mv   a.out defector.x)。
如果你用g77的话,那你的装g77的编译器。
2楼2017-03-05 21:51:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meizhixin

新虫 (小有名气)

引用回帖:
2楼: Originally posted by valenhou001 at 2017-03-05 21:51:50
你有该defector.f的源代码吗?
如果有的话,在你的linux 终端下,运行命令 gfortran   -o defector.x defector.f也是可以的。
或者运行gfortran defector.f命令得到a.out,然后将a.out改名为defector.x  (mv   a. ...

我用的defvector.f就是老师您写的采用vasp如何计算晶体弹性常数里面的,我用您说的方法还是不行,出现下面图示的情况,我百度了下说是fortran编译错误,建议用gfortran -ffree-form 加defvctor.f最后得到a.out,但是得到的a.out里面是乱码,这是对的吗?
能量应变法求弹性常数问题



发自小木虫Android客户端
3楼2017-03-06 10:07:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
3楼: Originally posted by meizhixin at 2017-03-06 11:07:21
我用的defvector.f就是老师您写的采用vasp如何计算晶体弹性常数里面的,我用您说的方法还是不行,出现下面图示的情况,我百度了下说是fortran编译错误,建议用gfortran -ffree-form 加defvctor.f最后得到a.out,但 ...

你是直接从pdf文件上拷贝下来的吗?

“但是得到的a.out里面是乱码”
a.out是二进制码,你直接编辑它,当然看到的是一堆乱码。这个就是编译得到的可执行文件。

» 本帖已获得的红花(最新10朵)

4楼2017-03-06 10:27:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
祝福
5楼2017-03-06 11:10:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meizhixin

新虫 (小有名气)

送红花一朵
引用回帖:
4楼: Originally posted by valenhou001 at 2017-03-06 10:27:05
你是直接从pdf文件上拷贝下来的吗?

“但是得到的a.out里面是乱码”
a.out是二进制码,你直接编辑它,当然看到的是一堆乱码。这个就是编译得到的可执行文件。...

别人分享的,我直接在网上下载的脚本,已经可以计算了,谢谢老师!
6楼2017-03-06 12:27:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小虫虫LG

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
4楼: Originally posted by valenhou001 at 2017-03-06 10:27:05
你是直接从pdf文件上拷贝下来的吗?

“但是得到的a.out里面是乱码”
a.out是二进制码,你直接编辑它,当然看到的是一堆乱码。这个就是编译得到的可执行文件。...

侯博士你好,现在用vasp自己可以直接算弹性常数了,那么vasp直接算和这里手动算,是一样的吗
7楼2017-03-06 14:56:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
7楼: Originally posted by 小虫虫LG at 2017-03-06 15:56:03
侯博士你好,现在用vasp自己可以直接算弹性常数了,那么vasp直接算和这里手动算,是一样的吗...

不一样。这个defvector.f手动计算,是基于能量更应变的关系。 vasp自带的是基于应力和应变,基于线性响应的方法来计算的。
两种计算出来的弹性常数在数值上略有差别,但是一般上它们应该是很接近。

» 本帖已获得的红花(最新10朵)

8楼2017-03-06 15:02:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小虫虫LG

木虫 (正式写手)

送红花一朵
引用回帖:
8楼: Originally posted by valenhou001 at 2017-03-06 15:02:20
不一样。这个defvector.f手动计算,是基于能量更应变的关系。 vasp自带的是基于应力和应变,基于线性响应的方法来计算的。
两种计算出来的弹性常数在数值上略有差别,但是一般上它们应该是很接近。...

知道了,非常感谢!
9楼2017-03-06 15:06:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 meizhixin 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见