小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
pwscf的原子坐标ATOMIC_POSITIONS设置
pwscf中输入原子坐标部分，我能直接用ms的*.cell文件中的原子坐标吗？能的话ATOMIC_POSITIONS后面{}我要填什么啊？
之前填的是crystal，但是我用xcrysden导入input文件，出来的结构跟ms里差很多
xcrysden中出现的
ms中建的模型
ecutwfc我大概测试了下，50-60左右，这是第一步加密K点算的，所以比较密。我确定ibrav的时候从materials studio里看行吗？
我没用过materials studio，不晓得是否能从那里面直接看（能不能直接用肯定也取决于你在QE的input里的具体设置），不过至少可以作为一个参考的。 ibrav具体怎么取还要看具体情况，而且可能性往往不止一种。比如这里你也可以设置ibrav=0，然后指定crystal lattice vector，或者还可以把你的结构当作三斜的，指定a, b, c, alpha, beta, gamma，看个人爱好，以及如何设置方便。不过不管怎么说，楼主都要对程序里的结构参数以及要做的体系的实际结构有一个了解，才方便进行设置的。
关于ecutwfc，我不知道楼主为什么要加密K点，但我不认为加密K点就需要提高ecutwfc，貌似没有太大关系。
我不是说提高ecutwfc是为的加密K点。我的ecutwfc经过收敛测试选取的；我是计算电声耦合，这里截取的input是第一步加密K点的计算
再请问，pwscf中的输出文件中，能查看所建结构的空间群吗？
麻烦你了，谢谢~
不能的，要用其他的软件分析查看。
比如？我现在只知道xcrysden可以导入，但是它貌似也不能看空间群吧
一般的找对称性的软件，比如findsym。需要你利用原子坐标设置一个findsym的输入文件。
嗯，谢谢你啦，我试试看
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
PWSCF中设置原胞坐标
在PWSCF中如何设置下图的原胞坐标？
还有这个呢？
好的，谢谢。
当然还可以根据原胞结构来设置的。
是啊，谢谢。
是的，这也是办法吧。
你好，请问下，在xcrysden1.4.1下看不到crystal95/96/03/06 input 这个的，就是无法打开；
想装1.5.24的xcrysden，就是装不上，不知道问题在哪里？
我的系统是Red Hat Enterprise 6.1。
你好，为什么我配置的包是这样到结果呢？
这可是花了好长时间才整出来的啊。
我的Xcrysden这一项是灰色的，无语了。。。
没***crystal，用xcConfigure设置XCrySDen时也没指定crystal软件包位置，自然那一项是灰色的了～
说实话，很不喜欢这种回复
还是希望多有souledge这样的人在，顶一下二楼的帖子！
是哦，其实是不知道，才问的啊。
这方面的资料有好多啊，注意去看看吧。
谢谢.请问 用XCrySDen可以建CNT的模型么?
请问Souledge专家，当ibrav=0利用惯用晶胞建立晶体结构时所得的结果找不到全部的对称性，怎么解决呢？
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研[转载]PWscf计算IR和Raman谱的步骤和注意事项
1、建立你所需要的晶体结构
在这步中，最重要的是测试
2）pw.x是否能找到晶体的全部对称性；
2、vc-relax&
&relax等，这取决于你是否想让你计算的结构加压力，或者应力什么的。
前面两步是前提，但这里就不多说了，我主要谈谈如何计算红外和拉曼谱
3、用优化好的晶体结构进行scf计算
&&&title =
&& & prefix =
&&&calculation =
restart_mode = 'from_scratch' ,
&& &wf_collect =
&& & outdir =
&& &pseudo_dir =
&etot_conv_thr = 1.0d-5 ,
&forc_conv_thr = 1.0d-4 ,
&& &tstress =
&& &tprnfor =
&&&nstep = 100
&&&ibrav =
&& & celldm(1) =
&& & celldm(3) =
&& & nat =
&& &ntyp =
&& &ecutwfc =
&&&nosym =
&ELECTRONS
&electron_maxstep = 100,
&&&conv_thr =
ATOMIC_SPECIES
&16-S.GGA.fhi.UPF
Zn&&65.409&
&30-Zn.GGA.fhi.UPF
ATOMIC_POSITIONS crystal
&-0.&&-0.&
K_POINTS automatic
&&16 16 16 0 0 0
结果文件中，有
& & 12 Sym.Ops. (no
inversion)
这说明pw.x找到了12个对称操作。如何判断pw.x找到了全部对称操作呢？我是根据materials
studio判断的，这个问题有空再讨论。
总之，pw.x一定要找到全部的对称操作才行。
除了要找到全部对称性外，另外就是使用较密的k-mesh，为后续的声子性质计算做准备；
还有一个需要注意的，我这里用了模守恒赝势，是从abinit网站上下的，然后用fhi2upf.x工具转换的，在pwtools下有。
4、用ph.x计算gamma点的声子性质，即一阶红外和拉曼性质
ph.x & ph.in & ph.out
ph.in 文件如下
phonon calculation for ZnS
&&tr2_ph=1.0d-12,
&&prefix='ph',
&&epsil=.true.,
&&trans=.true.,
&&lraman=.true.,
&&amass(1)=32.066,
&&amass(2)=65.409,
&&outdir='./',
&&fildyn='ZnS_IR.dynG',
&&fildrho='ZnS_IR.drho',
0.0 0.0 0.0
这里我要说明下
&&tr2_ph=1.0d-12,
这好像是默认值，可以调高些，计算量增加的不多，但是能够提高声子的计算精确度；
&&epsil=.true.,
&&trans=.true.,
&&lraman=.true.,
&&这三个设置成true，就能计算红外和拉曼了
&&fildyn='ZnS_IR.dynG',
这个是输出的动力学矩阵文件，在后处理中需要对该文件处理，大家可以看看里面是什么内容；
下面贴出ph.out 文件的最后输出信息
**************************************************************************
&&&omega( 1)
& 1.007969 [THz] =&
&& &33.622437
&&&omega( 2)
& 1.274035 [THz] =&
&& &42.497508
&&&omega( 3)
& 1.274035 [THz] =&
&& &42.497508
&&&omega( 4)
& 2.280533 [THz] =&
&& &76.070895
&&&omega( 5)
& 2.280533 [THz] =&
&& &76.070895
&&&omega( 6)
& 6.220740 [THz] =&
&&&207.502947
&&&omega( 7)
& 8.545154 [THz] =&
&&&285.037559
&&&omega( 8)
& 8.616679 [THz] =&
&&&287.423368
&&&omega( 9)
& 8.616679 [THz] =&
&&&287.423368
&&&omega(10)
& 8.870194 [THz] =&
&&&295.879779
&&&omega(11)
& 8.870194 [THz] =&
&&&295.879779
&&&omega(12)
&10.242682 [THz] =&
&&&341.661368
**************************************************************************
symmetry, C_6v (6mm)&&point
&&&omega(&&1
&33.6&&[cm-1]&
&&&omega(&&2
&42.5&&[cm-1]&
&&&omega(&&4
&76.1&&[cm-1]&
&&&omega(&&6
&&&207.5&&[cm-1]&
&&&omega(&&7
&&&285.0&&[cm-1]&
&&&omega(&&8
&&&287.4&&[cm-1]&
&&&omega( 10 -
&&&295.9&&[cm-1]&
&&&omega( 12 -
&&&341.7&&[cm-1]&
**************************************************************************
频率大小，对称性都给出了，不可约表示也给出来了
但是红外，拉曼的强度没有给出
下面进行下一步！
5，用dynmat.x 计算红外和拉曼强度
dynmat.x & ZnS001.in &
ZnS001.out
ZnS001.in文件
fildyn='ZnS_IR.dynG',
asr='crystal',
q(1)=0,q(2)=0,q(3)=1,
filout='ZnS_001.out'
filmol='ZnS_001.mold'
filxsf='ZnS_001.axsf'/
计算得到的
ZnS001.out文件
&&&IR cross
sections are in (D/A)^2/amu units
&&&Raman cross
sections are in A^4/amu units
&&&multiply
Raman by 0.140845 for Clausius-Mossotti correction
&-0.0000& &
0.0000& &&
&0.0155& &
0.0000& &&
&0.0372& &
0.0000& &&
&0.0132& &
& 2.2805& &
0.0000& &&
&0.1087& &
& 2.2805& &
0.0000& &&
&0.1087& &
& 207.50& &
0.0000& &&
&0.0000& &
& 288.35& &
8.1445& &&
&1.6782& &
& 288.35& &
8.1445& &&
&1.6782& &
& 295.88& &
0.0000& &&
&2.3796& &
& 295.88& &
0.0000& &&
&2.3796& &
&10.2427& &
0.0000& &&
&0.0000& &
&10.5003& &
9.0605& &&
&3.2269& &
IR 和 Raman 栏里就是表示计算的强度
可以自己写个小的进行高斯展宽的程序，那样就能画个拉曼谱了
如果有需要这个小程序的，可以站内消息联系。这里就不多谈这个小程序了
我的感觉：
1，计算红外和拉曼对系统的要求比较高，不是什么系统都是适合用来计算红外拉曼性质的。
比如强关联的体系，因为不好找赝势，而且需要很多测试。
2，红外和拉曼的计算量还是挺大的，ph.x的运行时间比较久，所以大的系统，比如原子数比较多，或者对称性很低的晶体结构，
可能需要的时间更多。
3，小分子的红外和拉曼我没有计算过，但我觉得还是可以计算的。不过有些特殊的结构，比如有5重对称性的分子，C60，估计不好
计算。因为晶体群对称性中没有5重对称性，这样的话，由于对称性的减小，计算红外和拉曼将要花太多的时间。如果有计算过小分子
红外和拉曼的虫友，可以互相交流下。
4，提供下我发表的几篇简单的工作和我经常看的和红外拉曼计算相关的文献
& Y. C. Cheng, H. T. Chen, X. X. Li, X. L. Wu, J.
Zhu, S. H. Li, and P. K. Chu, J. Appl. Phys. 105, 083511
& Y. C. Cheng, C. Q. Jin, F. Gao, X. L. Wu, W.
Zhong, S. H. Li, and P. K. Chu, J. Appl. Phys. 106, 123505
& J. Sun, X. F. Zhou, J. Chen, Y. X. Fan, H. T.
Wang, X. J. Guo, J. L. He, and Y. J. Tian, Phys. Rev. B 74, 193101
& M. Lazzeri, and F. Mauri, Phys. Rev. Lett. 90,
& S. Baroni, S. de Gironcoli, A. Dal Corso, and P.
Giannozzi, Rev. Mod. Phys. 73, 515 (2001).
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。