24小时热门版块排行榜    

查看: 1761  |  回复: 6

yiniangeqiu

新虫 (初入文坛)

[求助] 怎么用晶格动力学 计算phonon色散关系

小弟刚开始进入模拟领域。老板让我学习:
Write your own lattice dynamics code to calculate dispersion relations.
First, start with 1Dimensional, 2D and 3D cubic systems with springs.
搜素本论坛后也没有发现有用帖子.

对于1D chain with linear springs, 有公式,推理过程也可理解。
对于2D, 3D,不理解我需要求什么?
问题1: 我最终需要得到什么啊? 3-dimensional的色散关系图么?
问题2: 怎么得到阿? 是说由一个particle有波动,引起其他particle动了么?如果是这样,怎么得色散关系阿?
问题3: 哪些可以设置为已知量?哪些是未知量呢?已知量是spring constant, atom position么?
并求其他参考资料.
我是机械系刚开始进入分子动力模拟。所以,这方面理论不足。
多谢。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

longwei1221

铁虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
zh1987hs: 金币+2, 谢谢 2012-10-20 16:51:45
yiniangeqiu: 金币+2, 有帮助, 有帮助,虽然我不懂,但仍然奉上积分2 2012-10-21 01:00:57
我是看到phonon才点击的,我用castep计算过声子,可是晶格动力学没有专门去研究,计算得到的结果与一篇题目是“氟化锂晶格动力学”的结果是一样,但是我的输入文件是在CSD(http://www.crystallography.net/search.html )下载的cif文件然后倒入到MS中的,不知道是否对你有用
愿得一人心,白首不相离
2楼2012-10-20 16:46:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yiniangeqiu

新虫 (初入文坛)

那么你计算声子色散关系么?
什么是输入量?什么是输出量呢?
3楼2012-10-21 01:02:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

longwei1221

铁虫 (小有名气)

★ ★
chaizhm: 金币+2, 谢谢~ 2012-10-21 21:07:24
引用回帖:
3楼: Originally posted by yiniangeqiu at 2012-10-21 01:02:13
那么你计算声子色散关系么?
什么是输入量?什么是输出量呢?

在castep的task里面选择能量,然后在properties里面把phonon和density of state选项打上勾就可以了。输出的就是态密度和声子。你再点击得到的计算结果,castep-analysis里面选择你要分析的原子的声子谱就可以了,但愿对你有帮助
愿得一人心,白首不相离
4楼2012-10-21 18:41:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yiniangeqiu

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by longwei1221 at 2012-10-21 18:41:35
在castep的task里面选择能量,然后在properties里面把phonon和density of state选项打上勾就可以了。输出的就是态密度和声子。你再点击得到的计算结果,castep-analysis里面选择你要分析的原子的声子谱就可以了,但 ...

多谢。
5楼2012-10-23 23:06:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tangxian00

金虫 (著名写手)

老板要你自己编程算啊,比较麻烦。
一般是要求解矩阵的特征值,而矩阵要自己用code写出来。
要先学会计算机语言,比如fortran,C++等。
6楼2013-04-05 22:31:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

随便注册

银虫 (初入文坛)

参阅一般固体理论教材,对应1-3D体系,将晶格系统的动力学方程写成力常数和动力学矩阵的形式。构建好力常数,进行Fourier变换,写出动力学矩阵,对角化,就得出了一些列本征值。这些本征值是依赖于K的函数-----omega of K (K是晶格的波矢,或者是声子的波矢),它就是声子谱。大多数声子谱计算用的是第一性原理的大型软件。
根据导师给的任务,是让你做一些cross-filed Research ,相当开放的工作(虽然只能计算一些简单体系)。
所以,我目测楼主是在国外。
7楼2013-05-23 09:20:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yiniangeqiu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856材料与化工,270求调剂 +3 YXCT 2026-03-01 3/150 2026-03-01 19:51 by 无懈可击111
[考研] 材料学硕318求调剂 +9 February_Feb 2026-03-01 11/550 2026-03-01 19:47 by 无懈可击111
[考研] 265分求调剂不调专业和学校有行学上就 +5 礼堂丁真258 2026-02-28 7/350 2026-03-01 19:12 by Js512888
[考研] 298求调剂 +6 axyz3 2026-02-28 6/300 2026-03-01 19:00 by 18137688336
[考研] 一志愿中南大学理学化学 +4 15779376950 2026-03-01 5/250 2026-03-01 19:00 by Fff-1
[考研] 0856材料求调剂 +11 hyf hyf hyf 2026-02-28 12/600 2026-03-01 18:57 by 18137688336
[考研] 0857调剂 +3 一ll半 2026-02-28 3/150 2026-03-01 18:32 by 热情沙漠
[考研] 材料学调剂 +9 提神豆沙包 2026-02-28 11/550 2026-03-01 18:15 by ms629
[考研] 281求调剂 +4 2026计算机_诚心 2026-03-01 7/350 2026-03-01 17:20 by 2026计算机_诚心
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 化工专硕342,一志愿大连理工大学,求调剂 +3 kyf化工 2026-02-28 4/200 2026-03-01 16:49 by yywzz
[考研] 313求调剂 +3 水流年lc 2026-02-28 3/150 2026-03-01 16:01 by 新能源达人
[考研] 307求调剂 +5 wyyyqx 2026-03-01 5/250 2026-03-01 15:21 by Fff-1
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[考研] 寻找调剂 +4 LYidhsjabdj 2026-02-28 4/200 2026-03-01 10:56 by sunny81
[考研] 311求调剂 +9 南迦720 2026-02-28 10/500 2026-03-01 10:55 by sunny81
[硕博家园] 2025届双非化工硕士毕业,申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[考研] 292求调剂 +3 yhk_819 2026-02-28 3/150 2026-02-28 21:57 by gaoxiaoniuma
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见