24小时热门版块排行榜    

查看: 2106  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
请填处理意见