24小时热门版块排行榜    

查看: 5075  |  回复: 10

zqlee

金虫 (小有名气)

[交流] 【讨论】npt压力控制问题 已有7人参与

使用三维周期边界条件,
npt控压,
稳定很快就稳定了,
但是压力变动很大,
咋么处理?
谢谢。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

超精密加工技术
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ghcacj(金币+1):谢谢 2010-12-21 14:36:28
能不能给出更详细的数据呢?比如使用什么控压方法,什么参数?压力变化究竟有多大?
2楼2010-12-21 11:01:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zqlee

金虫 (小有名气)

★ ★ ★ ★ ★
zh1987hs(金币+5):thank you 2010-12-21 21:20:51
这是我的程序,温度可以稳定住,但是压力变化范围在1000以上。
# 3D tersoff potential cutting

units          metal
boundary       p p p

atom_style     atomic
lattice        diamond 3.5667

region         box block 0 10 0 10 0 10
create_box     1 box
create_atoms   1 box

#set up the pair style
pair_style     tersoff
pair_coeff     * * c.tersoff C

neighbor       0.3 bin
neigh_modify   delay 5 check yes

mass           * 12


# init^M
velocity       all create 297 8000

#energy minimization
minimize       1.0e-10 1.0e-10 10000 10000

#relaxation
fix            1 all npt 297 297 0.1 xyz 1 1 0.1 drag 0.2


timestep       0.0005
dump           1 all atom 1000 init.lammpstrj

thermo         100
run            1000000
超精密加工技术
3楼2010-12-21 13:30:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wally8962

木虫 (著名写手)

★ ★
zh1987hs(金币+1):thank you 2010-12-21 21:21:03
zh1987hs(金币+1):thank you 2010-12-21 21:21:11
fluctuation,

毕竟桠强是宏观概念,在这个小尺度上这样的涨落是正常的。
4楼2010-12-21 15:54:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)


zh1987hs(金币+1):thank you 2010-12-21 21:21:08
固体压力涨落尤其明显
5楼2010-12-21 16:00:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zqlee

金虫 (小有名气)

那也就是说控制不住的,
有没有别的好办法,
谢谢。
超精密加工技术
6楼2010-12-23 07:42:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ghcacj(金币+3):谢谢 2010-12-23 15:54:32
压力的涨落本来就是很大的。你的情况是正常的。1000个bar还不到一个GPa。通常NVE模拟,压力涨落在GPa以内都是很正常的。
7楼2010-12-23 09:14:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

andong1988

金虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖
引用回帖:
3楼: Originally posted by zqlee at 2010-12-21 13:30:28:
这是我的程序,温度可以稳定住,但是压力变化范围在1000以上。
# 3D tersoff potential cutting

units          metal
boundary       p p p

atom_style     atomic
lattice        diamond 3.5667

r ...

楼主,你好啊。呵呵。请教一下您的势函数文件“C.tersoff”是从哪里找到的呢?我也想用c的势函数。可只有SiC.tersoff啊。
非常感谢。
哈哈哈哈
8楼2012-02-12 12:58:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rush7

金虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我到现在也没搞明白到底波动多少算稳定……
9楼2012-04-25 18:41:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
1248151楼: Originally posted by 老虎大王 at 2010-12-23 09:14:12:
压力的涨落本来就是很大的。你的情况是正常的。1000个bar还不到一个GPa。通常NVE模拟,压力涨落在GPa以内都是很正常的。

请教老虎大王,NVE系综是不控制压强的,所以压强涨落是正常的,但npt系综是控制压强的,从理论上来说,压强值的波动不应过大吧,在NPT上,波动的幅度多少算是合适呢?
10楼2012-04-27 16:24:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zqlee 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学一志愿985,分数349求调剂 +4 zxts12 2026-03-21 7/350 2026-03-22 09:57 by zxts12
[考研] 化学调剂 +5 yzysaa 2026-03-21 5/250 2026-03-21 22:12 by peike
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 北科281学硕材料求调剂 +5 tcxiaoxx 2026-03-20 5/250 2026-03-20 21:35 by laoshidan
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
信息提示
请填处理意见