| 查看: 2324 | 回復(fù): 10 | |||
liuchenhan新蟲 (小有名氣)
|
[求助]
有關(guān)create_atoms相關(guān)命令 已有1人參與
|
|
本人腳本中先讀結(jié)構(gòu)文件nanonew.xyz,原子類型1-12,然后利用create_atoms命令生成一些其他原子,想安排這些原子的類型為13,但總是報錯,請各位幫忙指導(dǎo)一下,實現(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 報錯信息如下: 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ā)表 | |
|---|---|---|---|---|
|
[考研] 學(xué)碩274求調(diào)劑 +6 | Li李魚 2026-03-26 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +4 | 0857求調(diào)劑. 2026-03-24 | 4/200 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
[考研] 317求調(diào)劑 +5 | 十閑wx 2026-03-24 | 5/250 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
[考研] 一志愿陜師大生物學(xué)071000,298分,求調(diào)劑 +5 | SYA! 2026-03-23 | 5/250 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +3 | IceGuo 2026-03-26 | 4/200 |
|
|
[考研]
|
鐘llll 2026-03-26 | 4/200 |
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機 2026-03-25 | 4/200 |
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。,。,。,。i 2026-03-23 | 6/300 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
[考研] 318求調(diào)劑 +5 | plum李子 2026-03-21 | 8/400 |
|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +4 | 這么名字咋樣 2026-03-22 | 5/250 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚ly 2026-03-21 | 11/550 |
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運來來來 2026-03-22 | 5/250 |
|