应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » 其他 » 怎么用晶格动力学 计算phonon色散关系
71/1返回列表
查看: 1760  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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么？
并求其他参考资料.
我是机械系刚开始进入分子动力模拟。所以，这方面理论不足。
多谢。
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-10-20 03:50:08
已阅   回复此楼   关注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 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856求调剂285 +8 吕仔龙 2026-02-28 8/400 2026-03-01 17:25 by 刘兵
[考研] 290求调剂 +9 材料专硕调剂； 2026-02-28 11/550 2026-03-01 17:21 by sunny81
[考研] 0856化工专硕求调剂 +10 董boxing 2026-03-01 10/500 2026-03-01 17:19 by 刘兵
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +4 虫虫虫虫虫7 2026-03-01 7/350 2026-03-01 16:52 by caszguilin
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 材料学硕318求调剂 +8 February_Feb 2026-03-01 10/500 2026-03-01 16:48 by sunny81
[考研] 0805总分292，求调剂 +4 幻想之殇 2026-03-01 4/200 2026-03-01 16:47 by ms629
[考研] 材料化工调剂 +9 今夏不夏 2026-03-01 10/500 2026-03-01 16:01 by hmn_wj
[考研] 304求调剂 +6 曼殊2266 2026-02-28 7/350 2026-03-01 15:14 by wjLi2017
[考研] 材料工程274求调剂 +3 Lilithan 2026-03-01 3/150 2026-03-01 14:58 by ms629
[考研] 303求调剂 +4 今夏不夏 2026-03-01 4/200 2026-03-01 14:46 by 嘟嘟小浣熊
[考研] 课题组接收材料类调剂研究生 +3 gaoxiaoniuma 2026-02-28 4/200 2026-03-01 14:30 by jjj三跨
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
[考研] 306分材料调剂 +3 chuanzhu川烛 2026-03-01 4/200 2026-03-01 12:32 by houyaoxu
[考研] 298求调剂 +5 axyz3 2026-02-28 5/250 2026-03-01 06:45 by 刘兵
[考研] 272求调剂 +4 田智友 2026-02-28 4/200 2026-03-01 06:43 by 刘兵
[考研] 材料调剂 +4 爱擦汗的可乐冰 2026-02-28 4/200 2026-03-01 00:38 by 猫猫球alter
[基金申请] 面上模板改不了页边距吧？ +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
信息提示
关闭
请填处理意见
关闭