| 查看: 810 | 回复: 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 |
回复此楼» 猜你喜欢
材料调剂
已经有10人回复
-
一志愿哈尔滨工业大学085600英一数二337分求调剂
已经有5人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有11人回复
-
求调剂，一志愿南京航空航天大学 ，080500材料科学与工程学硕
已经有9人回复
-
本科211，293分请求调剂
已经有4人回复
-
求调剂
已经有3人回复
-
085400电子信息319求调剂（接受跨专业调剂）
已经有5人回复
-
268求调剂
已经有9人回复
-
一志愿双非085502，267分，过四级求调剂
已经有3人回复
-
一志愿085404，总分291，四级已过，求调剂
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
