24小时热门版块排行榜    

查看: 4869  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zqlee

金虫 (小有名气)

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

使用三维周期边界条件,
npt控压,
稳定很快就稳定了,
但是压力变动很大,
咋么处理?
谢谢。
回复此楼
超精密加工技术
已阅   回复此楼   关注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的回帖
查看全部 11 个回答

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的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见