24小时热门版块排行榜    

查看: 5052  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 此成果不能导入原因:元数据必填信息不完整,可 进行补充。 +3 Kittylucky 2026-03-02 4/200 2026-03-02 09:48 by kingkocxr
[考研] 材料类求调剂 +11 wana_kiko 2026-02-28 14/700 2026-03-02 08:46 by 聪明的大松鼠
[考研] 322求调剂 +3 熊境喆 2026-03-01 3/150 2026-03-02 08:44 by houyaoxu
[考研] 材料复试调剂 +4 学材料的点 2026-03-01 5/250 2026-03-02 08:26 by houyaoxu
[考研] 欢迎采矿、地质、岩土、计算机、人工智能等专业的同学报考 +4 pin8023 2026-02-28 6/300 2026-03-02 06:35 by 汪!?!
[考研] 求调剂 +5 yunziaaaaa 2026-03-01 6/300 2026-03-01 23:57 by ccp273206157
[考研] 材料学硕318求调剂 +11 February_Feb 2026-03-01 13/650 2026-03-01 23:53 by ccp273206157
[考研] 一志愿郑大材料学硕298分,求调剂 +5 wsl111 2026-03-01 5/250 2026-03-01 23:45 by 暮雨星晴
[基金申请] 成果系统访问量大,请一小时后再尝试。---NSFC啥时候好哦,已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 26考研报考西工大材料308分求调剂 +3 weizhong123 2026-03-01 3/150 2026-03-01 21:42 by 公瑾逍遥
[考研] 0805总分292,求调剂 +7 幻想之殇 2026-03-01 7/350 2026-03-01 21:22 by 公瑾逍遥
[考研] 化工299分求调剂 一志愿985落榜 +5 嘻嘻(*^ω^*) 2026-03-01 5/250 2026-03-01 19:47 by 无际的草原
[考研] 295复试调剂 +3 简木ChuFront 2026-03-01 3/150 2026-03-01 14:27 by zzxw520th
[考研] 298求调剂 +9 人间唯你是清欢 2026-02-28 12/600 2026-03-01 14:23 by Ducount.Y
[考研] 材料284求调剂,一志愿郑州大学英一数二专硕 +10 想上岸的土拨鼠 2026-02-28 10/500 2026-03-01 14:12 by yc258
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
信息提示
请填处理意见