24小时热门版块排行榜    

查看: 4866  |  回复: 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 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见