| 查看: 2323 | 回復(fù): 10 | |||
liuchenhan新蟲 (小有名氣)
|
[求助]
有關(guān)create_atoms相關(guān)命令 已有1人參與
|
|
本人腳本中先讀結(jié)構(gòu)文件nanonew.xyz,原子類型1-12,然后利用create_atoms命令生成一些其他原子,想安排這些原子的類型為13,但總是報(bào)錯,請各位幫忙指導(dǎo)一下,實(shí)現(xiàn)我想要的功能,非常感謝。 units metal dimension 3 boundary s s s neighbor 2.0 nsq #force field ******************************************************* atom_style atomic read_data nanonew.xyz lattice custom 4.32 a1 1.0 0.0 0.0 a2 0.0 0.577350269 0.0 a3 0.0 0.0 1.0 basis 0.0 0.5 0.0 basis 0.3333333333333 0.5 0.0 basis 0.5 0.0 0.0 basis 0.83333333333 0.0 0.0 #lattice custom 3.52 a1 1.0 0.0 0.0 a2 0.0 0.577350269 0.0 a3 0.0 0.0 1.0 basis 0.0 0.0 0.0 region boxx block 0 3 0 3 0 2 #create_box 1 box create_atoms 13 region boxx mass 13 28.0 pair_style lj/cut 3.0 pair_coeff * * 0.024 3.4 velocity all create 300 58127640 mom yes rot yes dist gaussian units box #300K fix 3 all temp/berendsen 300 300 0.05 #set the T of mobile is 300K while the initial value is300K thermo_style custom step temp etotal vol #thermo_modify lost warn thermo 100 timestep 0.0005 dump id1 all custom 1000 dump.lammpstrj id type x y z compute mytempall all temp/com fix Tall all ave/time 1 100 100 c_mytempall file Tall.txt dump last all custom 40000 dump.last type x y z run 40000 報(bào)錯信息如下: LAMMPS (14 May 2013) Reading data file ... orthogonal box = (-0.1 -0.1 -0.1) to (28.7 40.7 100) 1 by 1 by 1 MPI processor grid 240 atoms Lattice spacing in x,y,z = 4.32 2.49415 4.32 ERROR: Illegal create_atoms command (../create_atoms.cpp:109) |
分子動力學(xué) |
|
應(yīng)該先設(shè)定盒子,然后再在盒子里加region,在region里加原子。參考命令: region boxx block 0 5 0 5 0 4 create_box 1 box region boxx block 1 4 1 4 1 3 create_atoms 13 region boxx 另外,一般是先建立box,然后才能read_data導(dǎo)入數(shù)據(jù),所以你需要調(diào)整輸入命令的順序! |
新蟲 (小有名氣)
新蟲 (小有名氣)
木蟲 (正式寫手)

新蟲 (初入文壇)
新蟲 (小有名氣)
新蟲 (小有名氣)
新蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
孅華 2026-03-22 | 17/850 |
|
|---|---|---|---|---|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +10 | 路癡小琪 2026-03-27 | 10/500 |
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩,求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 303求調(diào)劑 +7 | 安憶靈 2026-03-22 | 8/400 |
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研]
|
鐘llll 2026-03-26 | 4/200 |
|
|
[考研] 324求調(diào)劑 +4 | wysyjs25 2026-03-21 | 4/200 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +9 | holy J 2026-03-21 | 9/450 |
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。,。,。,。i 2026-03-23 | 6/300 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|