24小时热门版块排行榜    

查看: 1778  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 318求调剂 +4 plum李子 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[考研] 一志愿武理材料工程348求调剂 +3  ̄^ ̄゜汗 2026-03-19 5/250 2026-03-22 14:03 by  ̄^ ̄゜汗
[考研] 305分求调剂(食品工程) +3 Sxy112 2026-03-21 5/250 2026-03-22 11:35 by 杨杨杨紫
[考研] 一志愿华中科技大学071000,求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 269专硕求调剂 +5 金恩贝 2026-03-21 5/250 2026-03-21 22:37 by zhyzzh
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
信息提示
请填处理意见