| 1 | 1/1 | 返回列表 |
| 查看: 793 | 回復(fù): 0 | |||
| 【懸賞金幣】回答本帖問題,作者yxchen22將贈送您 15 個金幣 | |||
yxchen22新蟲 (初入文壇)
|
[求助]
raspa output Qst出現(xiàn)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,力場,小分子,和骨架結(jié)構(gòu)都是用的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 不知道這是為什么,做了很多測試和檢查,一直都算不出來數(shù)據(jù)。換成乙烯小分子的時候,算出來所有的block都是nan |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 081200計算機科學(xué)與技術(shù)數(shù)一英一306求調(diào)劑 +3 | intankt 2026-03-08 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 282分材料專業(yè)求調(diào)劑院校 +16 | 楓橋ZL 2026-03-09 | 23/1150 |
|
|
[考研] 279求調(diào)劑 +3 | 莫xiao 2026-03-10 | 4/200 |
|
|
[考研] 327分求調(diào)劑086 +4 | 西紅柿?小帥 2026-03-09 | 7/350 |
|
|
[考研] 0817化學(xué)工程319求調(diào)劑 +7 | lv945 2026-03-08 | 9/450 |
|
|
[考研]
|
likeihood 2026-03-06 | 16/800 |
|
|
[考研] 物理學(xué)求調(diào)劑 +4 | sx宋 2026-03-05 | 4/200 |
|
|
[考研] 294 英二數(shù)二物化 求調(diào)劑 +6 | 米飯團不好吃 2026-03-09 | 6/300 |
|
|
[考研] 310 070300化學(xué)求調(diào)劑 +4 | 撲風(fēng)鈴的貓 2026-03-08 | 5/250 |
|
|
[考研] 0832食品科學(xué)與工程293調(diào)劑 +3 | 東東不東 2026-03-07 | 3/150 |
|
|
[考研] 334求調(diào)劑 +5 | Trying] 2026-03-07 | 5/250 |
|
|
[考研] 考研求調(diào)劑,工科,總分304 +4 | 矮子不想長高了 2026-03-08 | 7/350 |
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-08 | 3/150 |
|
|
[考研] 070300化學(xué)求調(diào)劑292分 +3 | 打烊eee 2026-03-07 | 3/150 |
|
|
[考研] 一志愿武漢理工085601,初試301分,請問能調(diào)劑到湖北嗎 +3 | 肖yang 2026-03-06 | 3/150 |
|
|
[考研] 276求調(diào)劑 本科天津工業(yè)大學(xué) 一志愿北京化工大學(xué)0855 +4 | 路lyh123 2026-03-06 | 6/300 |
|
|
[考研] 求調(diào)劑推薦 +4 | 微辣不吃 2026-03-06 | 4/200 |
|
|
[考研] 求調(diào)劑 +4 | 呼呼?~+123456 2026-03-06 | 4/200 |
|
|
[考研] 復(fù)試調(diào)劑 +5 | 呼呼?~+123456 2026-03-05 | 5/250 |
|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|