| 查看: 920 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有10人回复
-
2026国自然函评费到账
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
