Skip to content

混合物GCMC：二氧化碳和甲烷在 MOF MATVAF 中的吸附

介绍

  • Bae 和共同作者发表了一篇使用混合配体MOF模拟 二氧化碳/甲烷分离的论文
  • 他们使用 MuSiC 代码进行了混合物 GCMC 模拟

✍ simulation.input

NumberOfInitializationCycles 100000
NumberOfEquilibrationCycles  0
NumberOfProductionCycles     100000

UseChargesFromCIFFile yes

UseMaxStep  yes
MaxStepPerCycle 1

RestartFile no
RandomSeed  0

NumberOfTrialPositions 10
NumberOfTrialOrientations 10

NumberOfBlocks 1
AdsorbateAllocateSpace 10240
NumberOfSimulations 1
SingleSimulation yes

InputFileType cif
FrameworkName MATVAF_clean-Bae-PartialCharge-fixed-Original-N2-Charge
UnitCells 0 3 3 2

ChargeMethod Ewald
Temperature  296
Pressure     500000

OverlapCriteria 1e5
CutOffVDW 12.8
CutOffCoulomb 12.0
EwaldPrecision 1e-6

Component 0 MoleculeName              CO2
            IdealGasRosenbluthWeight  1.0
            FugacityCoefficient       0.9730021872
            MolFraction               0.5
            TranslationProbability    1.0
            RotationProbability       1.0
            ReinsertionProbability    1.0
            IdentityChangeProbability 1.0
            SwapProbability           1.0
            CreateNumberOfMolecules   0

Component 1 MoleculeName              methane
            IdealGasRosenbluthWeight  1.0
            FugacityCoefficient       0.9894689168
            MolFraction               0.5
            TranslationProbability    1.0
            ReinsertionProbability    1.0
            IdentityChangeProbability 1.0
            SwapProbability           1.0
            CreateNumberOfMolecules   0

📝 特别注意事项！

  • gRASPA 支持使用 Peng-Robinson 状态方程计算逸度系数
    • 对每个成分使用关键字 FugacityCoefficient PR-EOS 调用 Peng-Robinson 系数计算
  • 使用 IdentityChangeProbability 1.0 改善采样
    • 不使用身份交换也可以，因为 GCMC 交换移动就足够了。
  • Bae 等人报告的值是 过量负载，而 gRASPA 目前报告的是 绝对负载
    • 💡 你可以尝试将 gRASPA 重启文件输入到 RASPA-2 并运行 零周期
  • 这将计算 gRASPA 重启文件中构型的过量负载。 查看 基础GCMC 基础以了解更多