| 查看: 2035 | 回复: 18 | ||||
| 本帖产生 1 个 模拟EPI ,点击这里进行查看 | ||||
xiaowu787木虫 (正式写手)
|
[交流]
【讨论】高温高压的体系需要注意的环节已有4人参与
|
|||
|
希望LAMMPS模拟的模型中心产生一个高温高压(有时间寿命)的局部区,然后观测温度和压力的向周围扩展的情况。可是温度高压强高时会出现lost atoms,现在 第一个疑问是:周期边界条件下的image atom怎么不起作用了--周期性的边界条件不是可以保证粒子数不变? 第二: 局部升温的选择: 1. fix all nve+fix 局部 temp/rescale,可是压强无法控制 2. fix all nve+fix 局部 nvt or npt? 符合实际情况吗? 3. fix all npt+fix 局部 nvt or npt? 符合实际情况吗? 4. fix all nvt+fix 局部 nvt or npt? 符合实际情况吗? 那种选择更好一些? |
» 收录本帖的淘帖专辑推荐
动力学或热力学计算 |
» 猜你喜欢
请问有评职称,把科研教学业绩算分排序的高校吗
已经有4人回复
孩子确诊有中度注意力缺陷
已经有12人回复
2025冷门绝学什么时候出结果
已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
AI论文写作工具:是科研加速器还是学术作弊器?
已经有3人回复
2026博士申请-功能高分子,水凝胶方向
已经有6人回复
论文投稿,期刊推荐
已经有4人回复
硕士和导师闹得不愉快
已经有13人回复
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
» 本主题相关价值贴推荐,对您同样有帮助:
高温高压反应釜防火措施
已经有5人回复
高温高压反应釜防火措施
已经有10人回复
求《安全血液和血液制品》原本或相关题库
已经有4人回复
【求助】一个高温高压催化试验台的若干问题
已经有10人回复
【转帖】大力推进煤矿瓦斯抽采利用
已经有7人回复
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
2楼2010-09-29 18:02:49
xiaowu787
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1411.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.1小时
- 虫号: 326816
- 注册: 2007-03-18
- 性别: GG
3楼2010-09-29 20:37:00
manutdhk
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 5 (幼儿园)
- 金币: 5939.4
- 散金: 70
- 红花: 5
- 帖子: 94
- 在线: 97.3小时
- 虫号: 983544
- 注册: 2010-03-27
- 专业: 核技术及其应用
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+4):谢谢 2010-09-29 22:31:27
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+4):谢谢 2010-09-29 22:31:27
|
记得先前有位版友也提到过 高温的时候 lost atoms。 我不是用lammps的,所以具体不清楚什么原因。但是据经验分析,lost atoms,是不是因为粒子速度过大,在步长之内的位移过大,以至于越出了边界? 温度的扩散容易模拟和控制,但是压强的扩散就麻烦一点了。一般来说,压强可以通过外加应力,或者改变晶格参数来实现。如果不是恒压,应该选取NPT的系综,而且体积大小、形状都允许发生变化。 建议,先做温度扩散,再做压力扩散,然后再综合考虑。 |
4楼2010-09-29 21:08:59
xiaowu787
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1411.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.1小时
- 虫号: 326816
- 注册: 2007-03-18
- 性别: GG
5楼2010-09-29 21:15:44
manutdhk
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 5 (幼儿园)
- 金币: 5939.4
- 散金: 70
- 红花: 5
- 帖子: 94
- 在线: 97.3小时
- 虫号: 983544
- 注册: 2010-03-27
- 专业: 核技术及其应用
6楼2010-09-29 21:33:56
manutdhk
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 5 (幼儿园)
- 金币: 5939.4
- 散金: 70
- 红花: 5
- 帖子: 94
- 在线: 97.3小时
- 虫号: 983544
- 注册: 2010-03-27
- 专业: 核技术及其应用
7楼2010-09-29 21:35:40
xiaowu787
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1411.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.1小时
- 虫号: 326816
- 注册: 2007-03-18
- 性别: GG
8楼2010-09-29 22:28:35
manutdhk
木虫 (小有名气)
- 模拟EPI: 1
- 应助: 5 (幼儿园)
- 金币: 5939.4
- 散金: 70
- 红花: 5
- 帖子: 94
- 在线: 97.3小时
- 虫号: 983544
- 注册: 2010-03-27
- 专业: 核技术及其应用
9楼2010-09-29 23:07:27
xiaowu787
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1411.9
- 散金: 95
- 红花: 1
- 帖子: 550
- 在线: 430.1小时
- 虫号: 326816
- 注册: 2007-03-18
- 性别: GG
★ ★ ★ ★
zh1987hs(金币+4):谢谢,鼓励交流 2010-10-01 15:44:50
zh1987hs(金币+4):谢谢,鼓励交流 2010-10-01 15:44:50
|
时间步长设置 0.01fs 300K初始弛豫,时间步长改为 1fs 后开始升温2000K(整体NVE+部分temp/rescale(速度标度))也可以进行。但是有个问题--由于整体NVE,导致在部分进行速度标度是,温度向周围传递的很快(形成连续的温度梯度),希望的温度梯度最好不是很连续--中间快速升温时,两端温度上升的越慢越好,当随着时间的推移,温度真正传递到这里时才出现真正的升温过程。 [ Last edited by xiaowu787 on 2010-9-30 at 09:54 ] |
10楼2010-09-30 09:51:16














回复此楼
