 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (4735)
- >導(dǎo)師招生 (692)
- >考研 (675)
- >基金申請(qǐng) (385)
- >文獻(xiàn)求助 (292)
- >考博 (178)
- >蟲(chóng)友互識(shí) (174)
- >休閑灌水 (104)
- >碩博家園 (82)
- >論文投稿 (68)
- >招聘信息布告欄 (56)
- >博后之家 (52)
- >論文道賀祈福 (35)
- >綠色求助(高懸賞) (24)
- >找工作 (22)
- >SciFinder/Reaxys (21)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4295 | 回復(fù): 4 | ||
Materialman新蟲(chóng) (初入文壇)
|
[求助]
Lammps運(yùn)行結(jié)果問(wèn)題！無(wú)數(shù)據(jù)，出現(xiàn)nan！ 已有1人參與
|
|
本人在進(jìn)行能量?jī)?yōu)化的時(shí)候，有輸出的結(jié)果但是輸出的數(shù)據(jù)都為nan,不知道是什么回事，請(qǐng)各位大神指點(diǎn)迷津！ 以下為相應(yīng)的in文件 #------模擬環(huán)境初始化------ boundary p p p #邊界周期性設(shè)定，自由非周期、自由非周期、自由非周期 units real #模型單元模式設(shè)定為real dimension 3 #幾何維度設(shè)定為3D atom_style full #原子模式設(shè)定為full neighbor 2.0 bin #設(shè)定neighbor截?cái)喟霃?力場(chǎng)截?cái)喟霃剑愃芃S中Buffer width） #------微觀模型的建立------ # 定義原子結(jié)構(gòu) pair_style lj/class2/coul/cut 10.0 10.0 #設(shè)定PCFF勢(shì)函數(shù)的對(duì)勢(shì)（類似MS中的Nonbond potential）信息：對(duì)勢(shì)包含相互作用為范德華力和庫(kù)侖力，范德華力和庫(kù)侖力截?cái)喟霃剑愃芃S中Cutoff distance） bond_style class2 angle_style class2 dihedral_style class2 improper_style class2 read_data ../bench/WJ/PP.lammps05 #------設(shè)置原子相關(guān)參數(shù)------ velocity all create 500 102486 rot yes dist & gaussian units box #------動(dòng)力學(xué)設(shè)置------ # 能量?jī)?yōu)化 fix NPT all npt temp 500.0 500.0 100.0 iso 1.0 1.0 1000.0 compute pe all pe variable en equal c_ke+c_pe fix 1 all ave/time 1 10000 10000 c_temp c_ke c_pe v_en file ../bench/WJ/data/relax_PE20.dat1 minimize 1.0e-7 0 1000000 1000000 thermo 100000 thermo_modify lost warn timestep 1 run 1000000 運(yùn)行結(jié)果： LAMMPS (8 Jul 2013) Scanning data file ... 1 = max bonds/atom 6 = max angles/atom 9 = max dihedrals/atom 4 = max impropers/atom Reading data file ... orthogonal box = (0 0 0) to (39.0755 39.0755 39.0755) 1 by 2 by 2 MPI processor grid 7320 atoms 7260 bonds 14400 angles 21060 dihedrals 9600 impropers Finding 1-2 1-3 1-4 neighbors ... 4 = max # of 1-2 neighbors 6 = max # of 1-3 neighbors 12 = max # of 1-4 neighbors 16 = max # of special neighbors WARNING: Resetting reneighboring criteria during minimization (../min.cpp:173) Setting up minimization ... Memory usage per processor = 20.985 Mbytes Step Temp E_pair E_mol TotEng Press Volume 0 500 43041.617 45689.834 99639.726 200249.35 59664.174 236 500 8978.4944 9795.6924 29682.461 64908.573 59664.174 Loop time of 40.9642 on 4 procs for 236 steps with 7320 atoms Minimization stats: Stopping criterion = energy tolerance Energy initial, next-to-last, final = 88731.4515056 18774.1886865 18774.1868369 Force two-norm initial, final = 3071.24 323.327 Force max component initial, final = 870.624 81.8178 Final line search alpha, max atom move = 1.5559e-08 1.273e-06 Iterations, force evaluations = 236 548 Pair time (%) = 22.6427 (55.2744) Bond time (%) = 7.31013 (17.8451) Neigh time (%) = 0.33576 (0.819642) Comm time (%) = 10.531 (25.7078) Outpt time (%) = 0 (0) Other time (%) = 0.144612 (0.353019) Nlocal: 1830 ave 2187 max 1564 min Histogram: 1 0 1 0 1 0 0 0 0 1 Nghost: 13004 ave 13494 max 12423 min Histogram: 1 0 0 0 0 1 1 0 0 1 Neighs: 819009 ave 1.10617e+06 max 669830 min Histogram: 2 0 0 1 0 0 0 0 0 1 Total # of neighbors = 3276035 Ave neighs/atom = 447.546 Ave special neighs/atom = 11.6721 Neighbor list builds = 7 Dangerous builds = 0 Setting up run ... Memory usage per processor = 19.4735 Mbytes Step Temp E_pair E_mol TotEng Press Volume 236 500 8978.4944 9795.6924 29682.461 64908.573 59664.174 100000 -nan 0 -nan -nan -nan -nan 200000 -nan 0 -nan -nan -nan -nan 300000 -nan 0 -nan -nan -nan -nan 400000 -nan 0 -nan -nan -nan -nan 500000 -nan 0 -nan -nan -nan -nan 600000 -nan 0 -nan -nan -nan -nan 700000 -nan 0 -nan -nan -nan -nan 800000 -nan 0 -nan -nan -nan -nan 900000 -nan 0 -nan -nan -nan -nan 1000000 -nan 0 -nan -nan -nan -nan 1000236 -nan 0 -nan -nan -nan -nan Loop time of 22495.5 on 4 procs for 1000000 steps with 7320 atoms 2.png |
回復(fù)此樓» 猜你喜歡
材料與化工求調(diào)劑
已經(jīng)有7人回復(fù)
-
304求調(diào)劑
已經(jīng)有3人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有12人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有12人回復(fù)
-
286求調(diào)劑
已經(jīng)有10人回復(fù)
-
324分 085600材料化工求調(diào)劑
已經(jīng)有3人回復(fù)
-
317求調(diào)劑
已經(jīng)有6人回復(fù)
-
307求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助一個(gè)關(guān)于Linux下的VASP故障的問(wèn)題
已經(jīng)有6人回復(fù)
- 用lammps做模擬時(shí)如將新修改的勢(shì)函數(shù)讓lammps讀取運(yùn)行
已經(jīng)有9人回復(fù)
- Findit運(yùn)行問(wèn)題
已經(jīng)有15人回復(fù)
- lammps運(yùn)行出錯(cuò)：ERROR: Illegal fix nvt/npt/nph command (fix_nh.cpp:317)
已經(jīng)有7人回復(fù)
- 修改lammps程序的neighbor命令
已經(jīng)有4人回復(fù)
- 急！lammps安裝修改makefile.g++文件后運(yùn)行出錯(cuò)。
已經(jīng)有12人回復(fù)
- 經(jīng)過(guò)MS模型導(dǎo)入到LAMMPS，在運(yùn)算時(shí)出現(xiàn)的錯(cuò)誤求助
已經(jīng)有22人回復(fù)
- ms構(gòu)建好的模型轉(zhuǎn)化為lammps data文件
已經(jīng)有5人回復(fù)
- 關(guān)于lammps中的in文件運(yùn)行時(shí)的錯(cuò)誤問(wèn)題查找
已經(jīng)有10人回復(fù)
- LAMMPS計(jì)算無(wú)反應(yīng)，nan
已經(jīng)有4人回復(fù)
- LAMMPS運(yùn)行結(jié)果出現(xiàn)nan
已經(jīng)有20人回復(fù)
- 怎樣輸出應(yīng)力值
已經(jīng)有6人回復(fù)
- 我這樣運(yùn)行LAMMPS軟件中自帶的例子有錯(cuò)誤嗎？
已經(jīng)有6人回復(fù)
- 【求助】LAMMPS-運(yùn)行報(bào)錯(cuò)
已經(jīng)有9人回復(fù)
- 【求助】LAMMPS運(yùn)行錯(cuò)誤：Out of range atoms - cannot compute PPPM
已經(jīng)有7人回復(fù)
- 【求助】lammps中多元化合物的晶格創(chuàng)建問(wèn)題
已經(jīng)有14人回復(fù)
- 【討論】lammps運(yùn)行eam/fs出錯(cuò)
已經(jīng)有38人回復(fù)
2樓2014-09-25 09:43:25
Materialman
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1.5
- 帖子: 21
- 在線: 31小時(shí)
- 蟲(chóng)號(hào): 2765757
- 注冊(cè): 2013-10-30
- 專業(yè): 高分子材料與環(huán)境
3樓2014-09-28 09:00:41
wangluwei
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2.5
- 帖子: 9
- 在線: 15.1小時(shí)
- 蟲(chóng)號(hào): 5784496
- 注冊(cè): 2017-02-28
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
|
請(qǐng)問(wèn)lj勢(shì)用的時(shí)候可以不涉及到具體的原子么？?jī)煞N類型的原子該怎么寫(xiě)呢？我先用了具體的原子這樣寫(xiě)的 # LJ 12-6 system with npt ensemble units lj atom_style atomic boundary p p p lattice fcc 1.073 lattice fcc 0.8 region box block 0 10 0 5 0 5 create_box 2 box create_atoms 1 box create_atoms 2 box mass 1 1.0 mass 2 1.0 velocity all create 0.01 872877 timestep 0.01 dump 1 all xyz 1000 lj.xyz pair_style lj/cut 2.5 pair_coeff 1 1 1.0 1.0 2.5 pair_coeff 2 2 1.6 1.1765 2.5 pair_coeff 1 2 0.2 1.0882 2.5 #pair_modify shift yes mix geometric neighbor 0.3 bin neigh_modify every 10 delay 0 check yes thermo 1000 fix 1 all npt temp 0.01 0.01 1.00 iso 0.0024 0.0024 100.0 run 50000 unfix 1 thermo 1000 運(yùn)行結(jié)果總出現(xiàn)nan怎么回事？不用具體的原子怎么寫(xiě)代碼呢？ 新手跪謝�。。�！ |
4樓2017-04-15 21:43:32
Stoudemire29
金蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 862.5
- 帖子: 235
- 在線: 55.6小時(shí)
- 蟲(chóng)號(hào): 5963450
- 注冊(cè): 2017-03-13
- 專業(yè): 結(jié)構(gòu)工程
5樓2018-11-10 11:11:56
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 Materialman 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|