| 查看: 10759 | 回复: 17 | |||||||||||||
| 【奖励】 本帖被评价14次,作者chemvasp增加金币 10.4 个 | |||||||||||||
| 本帖产生 1 个 1ST强帖 ,点击这里进行查看 | |||||||||||||
[资源]
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 |
» 猜你喜欢
垃圾破二本职称评审标准
已经有19人回复
职称评审没过,求安慰
已经有53人回复
毕业后当辅导员了,天天各种学生超烦
已经有5人回复
26申博自荐
已经有3人回复
A期刊撤稿
已经有4人回复
» 本主题相关价值贴推荐,对您同样有帮助:
关于vasp收敛和优化
已经有24人回复
ansys中接触单元的参数如何定义
已经有10人回复
vasp.5.2.12 说明书
已经有4人回复
正交设计试验后,是否以及怎么样进行后续的进一步优化?
已经有6人回复
在vasp计算里,有没有好的只优化原子位置和晶胞体积的方法?
已经有8人回复
vasp 那些参数设置可以影响到磁性计算的磁矩
已经有12人回复
dmol在对晶胞参数进行优化时是否应该固定晶胞角度(fix angles)
已经有12人回复
利用CASTEP对建立的晶胞进行几何优化的问题 (急!!!)
已经有13人回复
【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
【求助】vasp优化晶格参数a
已经有13人回复
【求助】侯老师的vasp入门的晶格优化部分没看懂,问一下
已经有25人回复
【求助】能否发一份用vasp计算超晶胞某个表面的输入文件给我
已经有7人回复
【求助】Castep中结构优化后晶胞参数差别较大 vasp高手也帮我看一下
已经有22人回复
【求助】vasp编译时,报错:cannot find -lmkl_em64t
已经有11人回复
7楼2012-06-28 09:13:48
2楼2011-07-22 17:08:16
3楼2011-07-23 08:08:35
4楼2011-07-23 12:21:24
5楼2011-07-30 15:47:27
6楼2012-05-22 11:17:08
9楼2012-08-01 00:07:12
10楼2012-08-01 00:34:09
11楼2012-08-18 16:58:31
12楼2012-11-07 08:27:16
13楼2013-10-21 22:12:53
14楼2015-01-01 22:55:12
15楼2015-01-31 19:44:25
16楼2019-05-14 16:01:47
17楼2024-07-22 16:30:53
|
本帖内容被屏蔽 |
18楼2024-07-23 00:35:50
简单回复
2012-07-25 16:41
回复
五星好评













(ICELL(I),I=1,3)
回复此楼