| 查看: 952 | 回复: 0 | |||
gswylq木虫 (正式写手)
|
[求助]
已知动力学矩阵求声子振动频率
|
|
QE中有一个例子计算Si的声子谱，下面是si.dyn_G文件，请问如何手工实现动力学矩阵对角化得到声子振动频率呢？（声子频率在下方） Dynamical matrix file 1 2 2 10.2000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 1 'Si ' 25598.3672624059 1 1 0.0000000000 0.0000000000 0.0000000000 2 1 0.2500000000 0.2500000000 0.2500000000 Dynamical Matrix in cartesian axes q = ( 0.000000000 0.000000000 0.000000000 ) 1 1 0.27605479 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.27605479 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.27605479 0.00000000 1 2 -0.27549986 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 -0.27549986 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 -0.27549986 0.00000000 2 1 -0.27549986 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 -0.27549986 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 -0.27549986 0.00000000 2 2 0.27605479 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.27605479 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.00000000 0.27605479 0.00000000 Diagonalizing the dynamical matrix q = ( 0.000000000 0.000000000 0.000000000 ) ************************************************************************** omega( 1) = 0.484378 [THz] = 16.157120 [cm-1] ( 0.254905 0.000000 0.128027 0.000000 -0.647018 0.000000 ) ( 0.254905 0.000000 0.128027 0.000000 -0.647018 0.000000 ) omega( 2) = 0.484378 [THz] = 16.157120 [cm-1] ( 0.659031 0.000000 -0.021585 0.000000 0.255367 0.000000 ) ( 0.659031 0.000000 -0.021585 0.000000 0.255367 0.000000 ) omega( 3) = 0.484378 [THz] = 16.157120 [cm-1] ( 0.026485 0.000000 -0.695085 0.000000 -0.127104 0.000000 ) ( 0.026485 0.000000 -0.695085 0.000000 -0.127104 0.000000 ) omega( 4) = 15.270849 [THz] = 509.380687 [cm-1] ( -0.489933 0.000000 -0.509869 0.000000 0.000000 0.000000 ) ( 0.489933 0.000000 0.509869 0.000000 0.000000 0.000000 ) omega( 5) = 15.270849 [THz] = 509.380687 [cm-1] ( 0.214690 0.000000 -0.206295 0.000000 0.641366 0.000000 ) ( -0.214690 0.000000 0.206295 0.000000 -0.641366 0.000000 ) omega( 6) = 15.270849 [THz] = 509.380687 [cm-1] ( -0.462465 0.000000 0.444383 0.000000 0.297740 0.000000 ) ( 0.462465 0.000000 -0.444383 0.000000 -0.297740 0.000000 ) ************************************************************************** |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有69人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有11人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
找到一些相关的精华帖子，希望有用哦~
pwscf中的动力学矩阵的单位
已经有9人回复
- 关于力常数矩阵与声子谱计算时的q网格
已经有10人回复
- Abinit计算声子问题
已经有10人回复
- 怎么消除虚频，细节
已经有10人回复
- 非gamma点处的虚频意味着什么？
已经有14人回复
- 请问VASP5.2中如何设置INCAR才可以计算（gamma点的）phonon frequency？
已经有17人回复
- 【求助】怎样看出体系Gamma点的每个频率对应的原子振动方向呢？
已经有11人回复
- 【求助】怎样区分LO模和TO模
已经有11人回复
- 【求助】用vasp怎么计算声子谱和电声相互作用矩阵元
已经有6人回复
- 【求助】考虑温度的情况下，已知结构如何计算应力矩阵呢？
已经有5人回复
- 【求助】IR spectrum分析
已经有9人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
