24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1809  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +4 拾光任染 2026-04-06 4/200 2026-04-11 12:01 by Delta2012
[考研] 305求调剂 +5 77Qi 2026-04-07 5/250 2026-04-11 11:45 by zhq0425
[考研] 267求调剂 +3 一个叶子zsw 2026-04-05 4/200 2026-04-11 11:23 by zhq0425
[考研] 调剂求收留 +18 果然有我 2026-04-10 19/950 2026-04-11 11:06 by may_新宇
[考研] 电子信息279求调剂,有书读就行 +5 wwwooden 2026-04-08 8/400 2026-04-11 09:30 by zhq0425
[考研] 346,工科0854求调剂,专硕 +7 moser233 2026-04-10 8/400 2026-04-11 08:52 by 猪会飞
[考研] 化学工程与技术324调剂 +23 孙常华 2026-04-09 25/1250 2026-04-11 00:07 by 骑牛渡寒江
[考研] 0856专硕求调剂 希望是a区院校 +21 好好休息好不好 2026-04-09 24/1200 2026-04-10 16:58 by luoyongfeng
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
[考研] 070300化学学硕311分求调剂 +18 梁富贵险中求 2026-04-04 20/1000 2026-04-09 11:18 by 哒哒哒呱呱呱
[考研] 求调剂 +3 猪肉墩粉条cc 2026-04-08 4/200 2026-04-09 10:05 by 猪肉墩粉条cc
[考研] 一志愿北京2,材料与化工308求调剂 +17 熊二想上岸 2026-04-04 18/900 2026-04-09 09:13 by leyan1127
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 277、学硕,求调剂 数一104, +11 瓶子PZ 2026-04-07 12/600 2026-04-07 23:30 by 一只好果子?
[考研] 计算机11408 287 求调剂 +3 LiLe5 2026-04-07 3/150 2026-04-07 23:15 by shanqishi
[考研] 0854电子信息319求调剂(接受跨专业调剂) +5 星星不眨眼喽 2026-04-05 6/300 2026-04-07 22:16 by hemengdong
[考研] 0854求调剂 +9 亨氏番茄沙司 2026-04-06 10/500 2026-04-07 14:37 by shdgaomin
[考研] 信工所11408 340分 本科西安交大自动化 +3 moontrek 2026-04-06 3/150 2026-04-07 09:56 by chongya
[考研] 考研调剂 +11 小sun要好运 2026-04-04 11/550 2026-04-05 08:02 by qlm5820
信息提示
请填处理意见