计算化学公社

 找回密码 Forget password
 注册 Register
Views: 14337|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[Lammps] lammps中一维周期性边界条件如何添加

[复制链接 Copy URL]

155

帖子

0

威望

2876

eV
积分
3031

Level 5 (御坂)

lammps中,我了解到只能通过slab去除z方向的周期性,像这种模型能不能实现?
另外有人建议我在CNT四周加很大的真空层以去除周期性,但是关于库仑力使用Ewald或者PPPm方法,这样的误差是无法排除的,有文献提到极化力矩矫正不知道软件里面能不能实现。

CNT.jpg (118.82 KB, 下载次数 Times of downloads: 61)

该CNT沿轴向具有周期性,其他两个方向飞周期性

该CNT沿轴向具有周期性,其他两个方向飞周期性

1093

帖子

6

威望

6269

eV
积分
7482

Level 6 (一方通行)

2#
发表于 Post on 2016-5-28 19:57:07 | 只看该作者 Only view this author
lammps没用过,如果你需要的话我可以告诉你NAMD里面怎么做

294

帖子

0

威望

3528

eV
积分
3822

Level 5 (御坂)

3#
发表于 Post on 2016-5-29 10:36:41 | 只看该作者 Only view this author
fhh2626 发表于 2016-5-28 19:57
lammps没用过,如果你需要的话我可以告诉你NAMD里面怎么做

我觉的可以讲讲,在NAMD中实现,学习一下

155

帖子

0

威望

2876

eV
积分
3031

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2016-5-29 13:21:10 | 只看该作者 Only view this author
fhh2626 发表于 2016-5-28 19:57
lammps没用过,如果你需要的话我可以告诉你NAMD里面怎么做

谢谢你的回复,那麻烦你讲一下NAMD的做法,大家都学习一下!

2

帖子

0

威望

214

eV
积分
216

Level 3 能力者

5#
发表于 Post on 2016-5-29 15:51:45 | 只看该作者 Only view this author
本帖最后由 onigirinofruit 于 2016-5-29 15:53 编辑

刚学lammps,用boundary命令只将x设为p,yz不设置为周期性不行吗?

155

帖子

0

威望

2876

eV
积分
3031

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2016-5-29 16:24:02 | 只看该作者 Only view this author
onigirinofruit 发表于 2016-5-29 15:51
刚学lammps,用boundary命令只将x设为p,yz不设置为周期性不行吗?

这样的话库伦相互作用就无法使用Ewald等远程力方法了

1093

帖子

6

威望

6269

eV
积分
7482

Level 6 (一方通行)

7#
发表于 Post on 2016-5-30 12:33:09 | 只看该作者 Only view this author
ruanyang 发表于 2016-5-29 10:36
我觉的可以讲讲,在NAMD中实现,学习一下

首先是建立无限长的CNT,一般用VMD先画有限长的,然后再PSF里面手动添加几个跨越周期性边界的键长和键角
然后NAMD可以直接设定一维或者二维的PBC。PME肯定是不能用了,可以用新的MSM方法,或者直接把cutoff改大点儿
Z轴两边的压力差可以用tclforce或者tclBC实现

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1

查看全部评分 View all ratings

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-23 06:40 , Processed in 0.206235 second(s), 31 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list