| 5 | 1/1 | 返回列表 |
| 查看: 2327 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
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ā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 302分初試 0854 +3 | 伶可樂 2026-04-02 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 348求調(diào)劑 +6 | 吳彥祖24k 2026-04-02 | 6/300 |
|
|
[考研] 求調(diào)劑!生物與醫(yī)藥專碩 +4 | 逆轉(zhuǎn)陸先生 2026-04-01 | 4/200 |
|
|
[考研] 265求調(diào)劑11408 +4 | 劉小鹿lu 2026-03-27 | 4/200 |
|
|
[考研] 材料考研調(diào)劑 +9 | Gs大王 2026-04-02 | 9/450 |
|
|
[考研] 生物學(xué)308分求調(diào)劑(一志愿華東師大) +5 | 相信必會光芒萬?/a> 2026-03-31 | 6/300 |
|
|
[考研] 377求調(diào)劑 +3 | RASKIN 2026-04-02 | 3/150 |
|
|
[考研] 化學(xué)工程專碩324分,一志愿中國礦業(yè)大學(xué)求調(diào)劑 +7 | 耿耿1314 2026-04-01 | 7/350 |
|
|
[考研] 085602化學(xué)工程268分蹲調(diào)劑 +8 | 月照花林。 2026-04-01 | 8/400 |
|
|
[考研] 298求調(diào)劑 +4 | 什么是胖頭魚 2026-03-30 | 6/300 |
|
|
[考研] 353求調(diào)劑 +4 | 拉鉤不許變 2026-04-01 | 4/200 |
|
|
[考研] 求調(diào)劑 +4 | DADA怪 2026-03-31 | 4/200 |
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級,六級521,求調(diào)劑 +10 | 美味可樂雞翅 2026-03-26 | 12/600 |
|
|
[考研] 358求調(diào)劑 +3 | 王向陽花 2026-03-31 | 3/150 |
|
|
[考研] 085601英二數(shù)二求調(diào)劑 總分325 +4 | 余航航 2026-03-31 | 4/200 |
|
|
[考研] 求收留 +8 | 1943443204 2026-03-28 | 8/400 |
|
|
[考研] 085600材料與化工調(diào)劑 +16 | kikiki7 2026-03-30 | 16/800 |
|
|
[考研]
|
Gymno 2026-03-30 | 6/300 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚 2026-03-27 | 12/600 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|