| 查看: 4646 | 回復(fù): 22 | |||
[交流]
求助:lammps非周期邊界 已有6人參與
|
|
向大家請教一下,lammps里的非周期邊界f,s,m是什么意思?物理意義是什么?手冊上只是形象的描述,不算解釋吧。 lammps里其他的命令是怎么利用非周期邊界條件的?只是用于原子丟失的判定嗎? 假如我設(shè)置一個很大的盒子,原子只占里面的很小部分,肯定不會丟失,這時f,s,m下是一樣的結(jié)果嗎? |




木蟲 (初入文壇)

|
in文件: dimension 3 units metal boundary p p f atom_style full neighbor 3.0 bin neigh_modify delay 5 every 1 check yes read_data gufinal2.data # force field pair_style hybrid airebo 2.5 0 0 lj/cut 10.0 pair_coeff * * airebo CH.airebo C C C C pair_coeff 1 2*4 lj/cut 0.00284 3.40 pair_coeff 2 3*4 lj/cut 0.00284 3.40 pair_coeff 3 4 lj/cut 0.00284 3.40 # define region before define group region rfix block INF INF INF INF INF -8 units box region r1 block INF INF INF INF -8 -4.5 units box region r2 block INF INF INF INF -4.5 -1 units box region r3 block INF INF INF INF -1 INF units box region th1 block INF INF INF 2 -8 INF units box region th2 block INF INF 31 INF -8 INF units box group gfix region rfix group g1 region r1 group g2 region r2 group g3 region r3 group gth1 region th1 group gth2 region th2 group gth union gth1 gth2 group gthg12 subtract gth g3 group gthg3 subtract gth gthg12 group g12 union g1 g2 group gf12 union gfix g12 # compute and velocity compute temp6 gthg3 temp/com variable vg3x equal vcm(g3,x)#x component of velocity of the center of mass of group g3 variable pg3x equal xcm(g3,x)#x component of position of the center of mass of group g3 variable pg3y equal xcm(g3,y)#y component of position of the center of mass of group g3 variable pg3z equal xcm(g3,z)#z component of position of the center of mass of group g3 velocity gfix set 0.0 0.0 0.0 units box velocity g12 create 10 4928459 rot yes dist gaussian units box velocity g3 create 10 4928459 rot yes dist gaussian units box # fixes fix 1 all nve fix 2 gthg12 langevin 10 10 0.3 699483 zero yes fix 3 gthg3 langevin 10 10 0.1 699483 zero yes fix_modify 3 temp temp6 fix 4 gfix setforce 0.0 0.0 0.0 #relax min_modify dmax 0.01 minimize 1.0e-6 1.0e-8 1000 100000 neigh_modify exclude group gfix gfix timestep 0.001 thermo 100 thermo_style custom step v_pg3x v_pg3y v_pg3z v_vg3x dump 1 all xyz 200 gufinal2.xyz dump 2 all atom 200 gufinal2.lammpstrj run 50000 velocity g3 set 1.00 NULL NULL units box sum yes run 100000 盒子尺寸: 0.000000000 85.2168 xlo xhi 0.000000000 34.4332 ylo yhi -8.500000000 18.257102924 zlo zhi 模擬的是碳納米管在石墨烯上運(yùn)動,z方向原子不會丟失。 只改 boundary p p f 這句, v_pg3x v_pg3y v_pg3z v_vg3x 結(jié)果不一樣。 修改data文件里的 -8.500000000 18.257102924 zlo zhi 使盒子尺寸稍微變化一點,結(jié)果也不一樣。 |




| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +9 | 申子申申 2026-03-19 | 15/750 |
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
[考研] 310求調(diào)劑 +3 | baibai1314 2026-03-16 | 3/150 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +7 | 錦鯉0909 2026-03-17 | 7/350 |
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 323求調(diào)劑 +3 | 洼小桶 2026-03-18 | 3/150 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 生物學(xué)調(diào)劑招人!。 +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|