24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2124  |  回复: 17

潇潇雨歇99

新虫 (小有名气)


[交流] 用LAMMPS模拟Ni3AL降温,体积和能量没有突变

我想在2500K到300K降温,每100K降温10W步,弛豫10W步,一共450W步。
但是第一次计算的时候发现体系在2500到2400K之间体积突变,明显不对,又换了个势函数,这下直接没有突变了。
附上in文件
# LAMMPS Ni3Al

units                metal
boundary        p p p
atom_style        atomic

read_data Ni3Al1000.dat


pair_style eam/alloy
pair_coeff * * Al-Ni2009.eam.alloy Ni Al


timestep        0.001
thermo                50
neighbor        0.6 bin
neigh_modify        every 20 delay 0 check yes


variable        N equal step

compute                3        all pe/atom
compute                4        all ke/atom
compute                5        all coord/atom 4.0


velocity        all create 2500 825577 dist gaussian

dump                1 all cfg 20000 a.*.cfg id type xs ys zs c_3 c_4 c_5
dump                2 all atom 20000 dump.2500

fix                        1 all npt temp 2500 2500 0.02 iso 0 0 0.2 drag 1
run                        100000
unfix                1
write_restart        restart.$N

undump                1
undump                2

label                loop

variable        x index 2500 2400 2300 2200 2100 2000 1900 1800 1700 1600 1500 1400 1300 1200 1100 1000 900 800 700 600 500 400
variable         y index 2400 2300 2200 2100 2000 1900 1800 1700 1600 1500 1400 1300 1200 1100 1000  900 800 700 600 500 400 300
dump                1 all cfg 5000 a$y.*.cfg id type xs ys zs c_3 c_4 c_5
dump                2 all atom 5000 dump.$y


# Cooling steps:
fix                        1 all npt temp $x $y 0.01 iso 0 0 0.2 drag 1
run                        100000
unfix                1
write_restart        restart.$y.$N

# Equilibrium steps:
fix                 2 all npt temp $y $y 0.02 iso 0 0 0.2 drag 1
run                        100000
unfix                2


undump                1
undump                2

write_restart        restart.$N.$y

next                x
next                y
jump         Ni3Al.in         loop
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

tableman

木虫之王 (文学泰斗)



潇潇雨歇99(金币+1): 谢谢参与
祝福楼主了。
3楼2016-04-30 00:03:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

完全是火

铜虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
为什么要逐步降温而不直接从2500降温到300K呢?
13楼2016-05-03 14:14:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫: 金币+0.5, 给个红包,谢谢回帖
你的Ni3Al1000.dat 是怎么得到的?如果是1000K下的Ni3Al,肯定没有变化。
如果初始模型没有问题,建议跑一下升温过程
14楼2016-05-03 15:54:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

潇潇雨歇99

新虫 (小有名气)


引用回帖:
14楼: Originally posted by abinitio at 2016-05-03 15:54:02
你的Ni3Al1000.dat 是怎么得到的?如果是1000K下的Ni3Al,肯定没有变化。
如果初始模型没有问题,建议跑一下升温过程

是用一个叫persons crystal data文件里找到Ni3Al的单晶胞,放进MS里超晶胞10*10*10得到一个4631个原子的文件,然后用VESTA换成XYZ格式再放进matlab导出成data。中间需要自己补一下masses信息和坐标范围。1000的意思是10*10*10的超晶胞,不是1000K
15楼2016-05-04 10:01:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

潇潇雨歇99

新虫 (小有名气)


引用回帖:
13楼: Originally posted by 完全是火 at 2016-05-03 14:14:18
为什么要逐步降温而不直接从2500降温到300K呢?

因为要弛豫过程中的数据呀,直接降下去的话没有平衡,没法处理数据
16楼2016-05-04 10:02:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

潇潇雨歇99

新虫 (小有名气)


突然发现variable  命令少复制了几行,不过这个不是重点。按理说降温过程会过冷,可是我第一次算的迅速形成非晶,双体分布函数第二峰都裂了。密度是7点多。
第二次没有凝固,始终是液态,密度一直是6点多。
建模方式是师兄告诉我的。他算的另一种金属,并没有什么问题
17楼2016-05-04 10:04:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

影歌happy

新虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
请问我想降温行成晶体怎么搞?

发自小木虫IOS客户端
18楼2017-05-15 19:44:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2016-04-29 23:44   回复  
潇潇雨歇99(金币+1): 谢谢参与
syhorchid4楼
2016-04-30 00:09   回复  
潇潇雨歇99(金币+1): 谢谢参与
xhmaohan5楼
2016-04-30 02:04   回复  
潇潇雨歇99(金币+1): 谢谢参与
2016-04-30 05:44   回复  
潇潇雨歇99(金币+1): 谢谢参与
Anita_1237楼
2016-04-30 07:23   回复  
潇潇雨歇99(金币+1): 谢谢参与
2016-04-30 07:28   回复  
潇潇雨歇99(金币+1): 谢谢参与
2016-04-30 07:33   回复  
潇潇雨歇99(金币+1): 谢谢参与
2016-04-30 07:40   回复  
潇潇雨歇99(金币+1): 谢谢参与
2016-04-30 07:45   回复  
潇潇雨歇99(金币+1): 谢谢参与
2016-04-30 07:50   回复  
相关版块跳转 我要订阅楼主 潇潇雨歇99 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 评审有感 +3 popular289 2026-05-18 6/300 2026-05-18 15:17 by gy116024
[基金申请] 青C资助名额大幅增加! +12 西葫芦炒鸡蛋 2026-05-13 16/800 2026-05-18 10:02 by Equinoxhua
[基金申请] 重磅!青年科学基金项目(C类)资助增幅预计超过50% +7 水和泥不是水泥 2026-05-13 10/500 2026-05-18 07:50 by 水和泥不是水泥
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 xx7gd5zq4e 2026-05-15 6/300 2026-05-17 19:36 by Equinoxhua
[考博] 找博士生导师 +4 小代想上岸 2026-05-15 5/250 2026-05-17 17:44 by 脱颖而出
[考博] 2026博士还有哪些学校有名额 +6 小王求读研 2026-05-15 7/350 2026-05-17 16:54 by 知音湖畔
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +13 dd921ww 2026-05-12 16/800 2026-05-17 12:38 by 于轩
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:11 by 11n4dfd8yn
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 7hlccowb3h 2026-05-15 4/200 2026-05-17 07:46 by 11n4dfd8yn
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:31 by 11n4dfd8yn
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 7/350 2026-05-17 06:55 by 11n4dfd8yn
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:25 by ue3ir18jc3
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:57 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:50 by ue3ir18jc3
[基金申请] 请问大佬b0816评完了吗 +4 市民华南虎 2026-05-12 8/400 2026-05-16 19:54 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-16 17:45 by j6b2pdz07o
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 5/250 2026-05-16 16:20 by czyzsu
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
请填处理意见