24小时热门版块排行榜    

查看: 749  |  回复: 0

seven23

新虫 (初入文坛)

[求助] 跑NPT出现如下错误,请大神们指点一下,谢谢

一个含1条高分子的盒子,先NVT,正常,再NPT出现如下错误,请大家看看如何解决,非常感谢啊!!!


DISCOVER Molecular Simulation Program

                Version:   2010.1
                Build:     Sep 22 2010
                Date:      Tue Mar 12 17:38:14 2013
                User Name: WINDOWS-IO9TQ04$
                Host Name: WINDOWS-IO9TQ04
                Host Type: Windows
                Threads:   1


---------------------------------------------------------------
Checked out license feature: MS_discover_site [for msi] (1 copy)
---------------------------------------------------------------

randomSeed is set to 81094

Line   7:BTCL> autoEcho off

---------------------------------------------------------------
Checked out license feature: MS_compass_site [for msi] (1 copy)
---------------------------------------------------------------

INPUT FILES
___________

File Type                Name
---------                ----

Forcefield                D:\ms\MATERI~1.5X6\share\Discover\res/compass.frc
Molecular data                thelast.mdf
Coordinate                thelast.car

Periodic Boundary Conditions
____________________________

                Length (A)                Angle (degrees)
                ----------                ---------------

                a 20.32590                alpha  90.00000
                b 20.32590                beta   90.00000
                c 20.32590                gamma  90.00000


MOLECULAR TOPOLOGY
__________________

Number of molecules:                     1
Number of residues:                     1
Number of atoms:                   916 (asymmetric unit: 916)
Number of atom types:                    10
Number of bonds:                   915
Number of consolidated angles:          1726
Number of consolidated torsions:  2383
Number of bond_bond_1_3s:          2383
Number of angle-angles:                  3198
Number of out-of-planes:            50

FORCEFIELD OPTIONS
__________________

Filename                  : compass.frc
Definition name                  : cff91
Version                          : 2.8
Last modification date          : 3/6/2007
# of automatic parameters : 0

NONBOND ENERGY CUTOFFS
______________________

                Cutoff (A)        Spline Width (A)        Buffer Width (A)

vdW                  9.50                     1.00                     0.50
Coulomb                  9.50                     1.00                     0.50

Summation method for vdW        : Group based
Summation method for Coulomb        : Group based
Dielectric                        : 1.00
NOTE: Pressure and energy tail corrections will be added to all subsequent energy evaluations


MOLECULAR DYNAMICS
__________________

Ensemble           :       NPT

Temperature        :    298.00 K
Control Method     :   Andersen, Collision frequency is every 265 steps

Pressure           :    0.0001 GPa
Stress             : Sxx =    0.0000 GPa Syy =    0.0000 GPa Szz =    0.0000 GPa
                   : Syz =    0.0000 GPa Sxz =    0.0000 GPa Sxy =    0.0000 GPa
Control Method : Berendsen , tau =   0.1443 psec , compressibility =   0.5000

Timestep           :      1.00 fs
Duration           : 100000.00 fs
Integration Method : Velocity Verlet

Initial Velocities :   Random Velocities from Boltzmann distribution
Initial Temp.      :    298.00 K


Error: Energy expression: syntax error in expression "(2.0 * 3.14159) / -1.#IND"
Error: Dynamics: Integrator in dynamics failed
    while executing
"expr (2.0 * $pi) / $vol"
    invoked from within
"set energyFactor [expr (2.0 * $pi) / $vol]..."
    (procedure "fullTailCorrection" line 15)
    invoked from within
"fullTailCorrection 9.50 0"
    invoked from within
"dynamics \
    time = 100000.00 \
    timestep = 1.00 \
    initial_temperature = 298.00 \
    ensemble = NPT \
    temperature_control_method = ander ..."

Total time used by DISCOVER: 6 mins 34.45 secs (394.45 secs)

Completion date: Tue Mar 12 17:44:49 2013
Exiting Discover: An ERROR has occurred!
回复此楼
祝福大家都好运····
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 seven23 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 7lpolszZVXgi 2026-08-14 4/200 2026-08-15 01:08 by 4wMiSEwB6436
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7lpolszZVXgi 2026-08-14 5/250 2026-08-15 00:52 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 7/350 2026-08-15 00:40 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] filecode +6 cratir 2026-08-14 10/500 2026-08-14 18:29 by 笑叹辞穷
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +5 工厂打螺丝 2026-08-08 5/250 2026-08-14 11:03 by 玖戈弋
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
信息提示
请填处理意见