| 查看: 879 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者yxchen22将赠送您 15 个金币 | ||
[求助]
raspa output Qst出现nan
|
||
|
我用raspa测试手册案例：LTA4A+96Na+CO2进行MC模拟 input文件如下： SimulationType MonteCarlo NumberOfCycles 25000 NumberOfInitializationCycles 10000 PrintEvery 1000 RestartFile no Forcefield GarciaPerez2006 RemoveAtomNumberCodeFromLabel yes ModifyOxgensConnectedToAluminium yes Framework 0 FrameworkName LTA4A UnitCells 1 1 1 ExternalTemperature 300.0 ExternalPressure 100000.0 component 0 MoleculeName sodium MoleculeDefinition TraPPE TranslationProbability 1.0 RandomTranslationProbability 1.0 ExtraFrameworkMolecule yes CreateNumberOfMolecules 96 Component 1 MoleculeName CO2 MoleculeDefinition TraPPE BlockPockets yes BlockPocketsFilename LTA_SI TranslationProbability 1.0 ReinsertionProbability 1.0 SwapProbability 1.0 ExtraFrameworkMolecule no CreateNumberOfMolecules 0 然后算出来的Qst是nan，力场，小分子，和骨架结构都是用的raspa自带的。并且LTA_SI用这个力场吸附CO2可以算出Qst； Enthalpy of adsorption: ======================= Enthalpy of adsorption component 0 [sodium] ------------------------------------------------------------- Block[ 0] 5163279119389143040.00000 [-] Block[ 1] -nan [-] Block[ 2] -3631306914516404736.00000 [-] Block[ 3] -nan [-] Block[ 4] -812704293754330112.00000 [-] ------------------------------------------------------------------------------ Average -nan +/- nan [K] -nan +/- nan [KJ/MOL] Note: Ug should be subtracted to this value Note: The heat of adsorption Q=-H Enthalpy of adsorption component 1 [CO2] ------------------------------------------------------------- Block[ 0] 2206708.68421 [-] Block[ 1] -nan [-] Block[ 2] -1261447.81955 [-] Block[ 3] -nan [-] Block[ 4] -315586.95489 [-] ------------------------------------------------------------------------------ Average -nan +/- nan [K] -nan +/- nan [KJ/MOL] Note: Ug should be subtracted to this value Note: The heat of adsorption Q=-H Total enthalpy of adsorption from components and measured mol-fraction ---------------------------------------------------------------------- Block[ 0] 3542102898865602048.00000 [-] Block[ 1] -nan [-] Block[ 2] -2485242509054886400.00000 [-] Block[ 3] -nan [-] Block[ 4] -554462315282042368.00000 [-] ------------------------------------------------------------------------------ Average -nan +/- nan [K] -nan +/- nan [KJ/MOL] Note: Ug should be subtracted to this value Note: The heat of adsorption Q=-H 不知道这是为什么，做了很多测试和检查，一直都算不出来数据。换成乙烯小分子的时候，算出来所有的block都是nan |
回复此楼» 猜你喜欢
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
有谁可曾问过你过的还好吗？
已经有21人回复
-
E0414, 我的本子有没有希望？
已经有7人回复
-
一篇论文同时出现在两个期刊，一模一样，这算不算学术不端，请各位老师斧正。
已经有12人回复
-
希望面上有个好结果
已经有7人回复
-
三区计算机方向期刊推荐
已经有5人回复
-
sci论文二审求助
已经有5人回复
-
函评
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
