在GCMC计算时出错，使用软件为RASPA

10yaoxue

_cell_length_a: 25.860400
_cell_length_b: 25.860400
_cell_length_c: 25.860400
_cell_length_alpha: 90.000000
_cell_length_beta: 90.000000
_cell_length_gamma: 90.000000
End reading cif-file
ReturnPseudoAtomNumber: Error!!!! :H_h2

附输入文件：SimulationType                MonteCarlo
NumberOfCycles                20000
NumberOfInitializationCycles  20000
PrintEvery                    500
RestartFile                   no

Forcefield                    TraPPE
ChargeMethod                  Ewald
EwaldPrecision                1e-6

Framework 0
FrameworkName MOF-5
UnitCells 1 1 1
HeliumVoidFraction 0.572
ExternalTemperature 273.15
ExternalPressure 1 10000 20000 30000 40000 50000 60000 70000 80000 90000 100000 110000 120000

Movies yes
WriteMoviesEvery 500

Component 0 MoleculeName                     hydrogen
            MoleculeDefinition               TraPPE
            TranslationProbability           0.5
            RotationProbability              0.5
            ReinsertionProbability           0.5
            SwapProbability                  1.0
            CreateNumberOfMolecules          0

一路向前

但是计算发现用这个软件计算出来的结果还是有很大偏差

一路向前

slxc920113 发表于 2023-11-8 14:33
https://autoff.readthedocs.io/en/latest/example.html#ce-ch-4-gcmc

谢谢

slxc920113

一路向前 发表于 2023-11-8 14:14
在用RASPA计算乙炔和CO2的时候，各位有没有好的计算参数，如怎么修改force_field_mixing_rules.def，C2H2.d ...

https://autoff.readthedocs.io/en/latest/example.html#ce-ch-4-gcmc

一路向前

在用RASPA计算乙炔和CO2的时候，各位有没有好的计算参数，如怎么修改force_field_mixing_rules.def，C2H2.def及pseduo_atom.def文件呢

一路向前

cdcdcd 发表于 2020-12-17 19:32
您好，我也遇到了类似的问题，请问您是怎么解决的呀

RASPA计算乙炔模拟需要修改力场参数

CJW

小瓶子 发表于 2017-1-22 22:32
楼主你好，我有个小问题想向你请教，在输入文件中，有
           TranslationProbability           0.5 ...

Swap 包含插入和删除两个操作，这里设为1这几个操作概率就一样了。当然没有必要保证这里的概率一样，大概只是为了写文章方便说。

NK-MOF

Zengwenmao 发表于 2022-1-22 00:56
您好，我也尝试把_后面的内容删除，但是请问这样会影响调用force_field_mixing_rules.def中对应的原子参 ...

可以在输入文件中加入关键词“RemoveAtomNumberCodeFromLabel yes”试一试。

Zengwenmao

NK-MOF 发表于 2021-9-14 15:34
把C2H2.def文件中的C_C2H2和H_C2H2改成C和H就可以了。

您好，我也尝试把_后面的内容删除，但是请问这样会影响调用force_field_mixing_rules.def中对应的原子参数吗？
比如：C2H2中的C原本应该调用C_C2H2的LJ参数，在去掉_C2H2后是否会调用C的参数？
是否有其他方法可以解决该问题？

NK-MOF

cdcdcd 发表于 2020-12-17 19:32
您好，我也遇到了类似的问题，请问您是怎么解决的呀

把C2H2.def文件中的C_C2H2和H_C2H2改成C和H就可以了。

cdcdcd

Yanglz 发表于 2019-11-3 18:35
您好，我在计算C2H2时也遇到了类似的问题，ReturnPseudoAtomNumber: Error!!!! :C_C2H2，请问您最后怎么解 ...

您好，我也遇到了类似的问题，请问您是怎么解决的呀

Yanglz

您好，我在计算C2H2时也遇到了类似的问题，ReturnPseudoAtomNumber: Error!!!! :C_C2H2，请问您最后怎么解决的？

senR_R

您好， 请问这个问题您最后是怎么解决的？

小瓶子

agent99 发表于 2017-1-23 06:51
我记得这些概率都是相对值吧

谢谢，学习了

agent99

小瓶子 发表于 2017-1-22 22:32
楼主你好，我有个小问题想向你请教，在输入文件中，有
           TranslationProbability           0.5 ...

我记得这些概率都是相对值吧