24小时热门版块排行榜    

查看: 8048  |  回复: 15
【奖励】 本帖被评价12次,作者chemvasp增加金币 9
本帖产生 1 个 1ST强帖 ,点击这里进行查看

chemvasp

新虫 (初入文坛)


[资源] VASP并行可执行软件包,可对晶胞参数进行部分优化

下载地址:http://ishare.iask.sina.com.cn/f/17331535.html

在优化表面,二维或一维体系时,必须固定某个或某几个晶胞方向不优化。
VASP本身没有这个功能,但是提供了constr_cell_relax.F文件来设置。
根据文件中的提示,添加几行代码后,就可以任意固定晶胞方向了。

使用方法:
在目录下建立文件OPTCELL,第一行输入三个数字(比如110),
数字1表示优化,数字0表示不优化;三个数字依次对应于
xyz方向。比如110表示z轴不优化;001表示x和y轴都不优化。
注1. 无OPTCELL文件则完全是原版软件的功能;
注2. 只适用与正交晶胞,所有角都是直角。

附:修改后的constr_cell_relax.F文件
      SUBROUTINE CONSTR_CELL_RELAX(FCELL)
      USE PREC
      REAL(Q) FCELL(3,3), SAVE(3)
      LOGICAL FILFLG
      INTEGER ICELL(3)

      INQUIRE(FILE='OPTCELL',EXIST=FILFLG)
      IF (FILFLG) THEN
        OPEN(67,FILE='OPTCELL',FORM='FORMATTED',STATUS='OLD')
          READ(67,"(3I1)" (ICELL(I),I=1,3)
        CLOSE(67)
        DO I=1,3
          SAVE(I)=FCELL(I,I)
        ENDDO
        FCELL=0.0d0
        DO I=1,3
          IF (ICELL(I)==1) FCELL(I,I)=SAVE(I)
        ENDDO
      ENDIF

!     just one simple example
!     relaxation in x directions only
!      SAVE=FCELL(1,1)
!      FCELL=0   ! F90 style: set the whole array to zero
!      FCELL(1,1)=SAVE

      RETURN
      END SUBROUTINE
回复此楼

» 收录本帖的淘贴专辑推荐

第一性原理计算-tg vasp计算资料 第一性原理计算经验 材料计算模拟软件
VASP资料 量化软件 有效质量 网贴收集
vasp soft study vasp

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cavediger

金虫 (正式写手)


★★★★★ 五星级,优秀推荐

挺好,有创新,支持
2楼2011-07-22 17:08:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雁儿霏霏

木虫 (正式写手)


★★★★★ 五星级,优秀推荐

强烈支持
3楼2011-07-23 08:08:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wutianhuan

铁杆木虫 (著名写手)


★★★★★ 五星级,优秀推荐

谢谢,相信会很有用的,呵呵
4楼2011-07-23 12:21:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dycding

铁杆木虫 (初入文坛)


★★★ 三星级,支持鼓励

请问我在CentOS5.5上按说明方法安装时出现了无法执行二进制文件vasp52mpi.x(cannot execute binary file)这是什么原因?有什么方法可以解决吗?谢谢!
5楼2011-07-30 15:47:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dragon5061

铁杆木虫 (正式写手)


引用回帖:
572389楼: Originally posted by chemvasp at 2011-07-22 15:11:45:
下载地址:http://ishare.iask.sina.com.cn/f/17331535.html

在优化表面,二维或一维体系时,必须固定某个或某几个晶胞方向不优化。
VASP本身没有这个功能,但是提供了constr_cell_relax.F文件来设置。
根据文

如果不是非正交晶胞的话,该如何实现呢?
6楼2012-05-22 11:17:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangll99

木虫 (职业作家)


★★★★★ 五星级,优秀推荐

请问楼主,如果三个数都选1,那么晶胞的边的夹角会不会变化?
7楼2012-06-28 09:13:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljzhou86

金虫 (正式写手)


引用回帖:
572389楼: Originally posted by chemvasp at 2011-07-22 15:11:45
下载地址:http://ishare.iask.sina.com.cn/f/17331535.html

在优化表面,二维或一维体系时,必须固定某个或某几个晶胞方向不优化。
VASP本身没有这个功能,但是提供了constr_cell_relax.F文件来设置。
根据文 ...

你好,我优化一个2D构型,(011),x轴不优化,放宽y,z轴优化,但最终x,y,z都在优化,似乎这于这种情况不适合,你检查下代码是哪有问题?当然如果是(110)这样的话,程序就没有问题。
9楼2012-08-01 00:07:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljzhou86

金虫 (正式写手)


★★★★★ 五星级,优秀推荐

你好,我优化一个2D构型,(011),x轴不优化,放宽y,z轴优化,但最终x,y,z都在优化,似乎这于这种情况不适合,你检查下代码是哪有问题?当然如果是(110)这样的话,程序就没有问题。
10楼2012-08-01 00:34:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hbl626

铁虫 (正式写手)


★★★★★ 五星级,优秀推荐

太牛了,谢谢
11楼2012-08-18 16:58:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liancsh

铜虫 (小有名气)


★★★★★ 五星级,优秀推荐

非常感谢!
12楼2012-11-07 08:27:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duxiu727

金虫 (小有名气)


★★★★★ 五星级,优秀推荐

引用回帖:
6楼: Originally posted by dragon5061 at 2012-05-22 11:17:08
如果不是非正交晶胞的话,该如何实现呢?...

同问,如果非正交呢
13楼2013-10-21 22:12:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyqh

木虫 (著名写手)


★ 一星级,一般

顶顶更健康
14楼2015-01-01 22:55:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangjl580

新虫 (初入文坛)


引用回帖:
10楼: Originally posted by ljzhou86 at 2012-08-01 00:34:09
你好,我优化一个2D构型,(011),x轴不优化,放宽y,z轴优化,但最终x,y,z都在优化,似乎这于这种情况不适合,你检查下代码是哪有问题?当然如果是(110)这样的话,程序就没有问题。

我也发现了 011肯定是有问题的  代码有问题 暂时还没学程序,看不懂哪里有问题
15楼2015-01-31 19:44:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

蓝凌小未

铜虫 (正式写手)


★★★★★ 五星级,优秀推荐

引用回帖:
10楼: Originally posted by ljzhou86 at 2012-08-01 00:34:09
你好,我优化一个2D构型,(011),x轴不优化,放宽y,z轴优化,但最终x,y,z都在优化,似乎这于这种情况不适合,你检查下代码是哪有问题?当然如果是(110)这样的话,程序就没有问题。

你好~我也遇到了这个问题,请问你是如何解决的,谢谢大神虫友~
16楼2019-05-14 16:01:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2012-07-25 16:41   回复  
五星好评  
相关版块跳转 我要订阅楼主 chemvasp 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 入职高校3年发表10+SCI,尽人事听天命 +31 kaoyan250 2024-05-27 42/2100 2024-06-02 12:09 by ECOblance
[考博] 导师不让硕转博,让我去国外读博,能理解吗? +12 萧山幽谷 2024-05-29 20/1000 2024-06-02 12:01 by yuan0806
[基金申请] 2024杰青和万人领军什么时候会评 +3 墨香琴韵 2024-06-02 3/150 2024-06-02 10:31 by 摩天思瑞
[硕博家园] 讨厌酒桌文化 +7 hahamyid 2024-05-31 7/350 2024-06-02 10:01 by 御龙老大
[硕博家园] 每到中夜,情难自抑 +34 sioc-sunj 2024-05-28 58/2900 2024-06-02 02:53 by csyky2007
[硕博家园] 每天学术时间不能保证,能保证的只有: +10 hahamyid 2024-05-27 10/500 2024-06-01 21:11 by 小小芝麻影
[考博] 24or25材料专业申博 +4 农夫三拳有点痛 2024-05-30 11/550 2024-06-01 14:45 by Napoleonsky
[基金申请] B口人才项目 +9 WOWO159357 2024-05-29 19/950 2024-06-01 14:24 by linxuhuizj
[论文投稿] 求Sci期刊推荐 10+4 甄小鱼 2024-05-30 7/350 2024-06-01 10:41 by bobvan
[基金申请] 博后特别资助状态变化 +24 随梦而飞2017 2024-05-30 35/1750 2024-06-01 10:10 by 青岛阳仔
[论文投稿] 纠结选哪一个期刊,电化学领域 50+8 Freya163 2024-05-28 10/500 2024-05-31 15:09 by wzykobe
[文学芳草园] 对对子啊 +5 天若孤独 2024-05-29 7/350 2024-05-31 09:00 by wjykycg
[文学芳草园] 物是人非 +4 myrtle 2024-05-30 4/200 2024-05-30 15:05 by mapenggao
[论文投稿] 《控制理论与应用》期刊收版面费吗? +4 ygj2015 2024-05-29 4/200 2024-05-30 06:59 by nono2009
[硕博家园] 又想换工作 +15 brightmj 2024-05-27 20/1000 2024-05-29 23:25 by zyqchem
[论文投稿] 真急着毕业,CPB主编终审17天了,邮件催稿了两次,就是一点动静没有 5+3 kkkk夏 2024-05-28 6/300 2024-05-29 11:18 by hitsdu
[基金申请] E05青基有几个评审 +4 KYXY123 2024-05-28 4/200 2024-05-28 19:25 by popt2t
[有机交流] 奇怪的物质 100+4 桃桃PXS 2024-05-27 7/350 2024-05-28 10:22 by 091602
[硕博家园] 课题做到一半担心做不出来怎么办 +4 mkpooooo 2024-05-27 7/350 2024-05-27 11:56 by kanmiaolucky
[硕博家园] 周日 +6 1加油哦棒 2024-05-26 9/450 2024-05-27 10:30 by hahamyid
信息提示
请填处理意见