Messages

thank you, sir

Hi Prof. Dubbeldam,

When using Monte Carlo method to carry out a simulation, i can gaint 'average loading absolute' and 'average loading excess' of absorbate. However i don't know the theory of obtaining the adsorption capacity by Monte Carlo method. What is the theory of this process, or what can I learn to understand the theory?

Thank you,
Jiang

Dear Dubbeldam,

I'm using RASPA to simulate the alkane adsorption in Montmorillonite. But in the output , the absolute adsorption is zero, the excess loading is negative, I don't know why and how to get the correct amount of adsorption.

Kind regards,
Yadi Jiang

*************
This is the component's critical parameter：792K, 0.98MPa, 0.8.
(In the simulation, temperature is 358K, pressure is 30MPa,40MPa.)

*************
This is the simulation.input file:

SimulationType                MonteCarlo
NumberOfCycles                100000
NumberOfInitializationCycles  50000
PrintEvery                    1000
RestartFile                   no

ContinueAfterCrash no
WriteBinaryRestartFileEvery 5000

ChargeMethod                    Ewald
Forcefield                      clayff  //Custom force field
CutOffVDW                       12.5
RemoveAtomNumberCodeFromLabel   no

Framework 0
FrameworkName 10-MMT
UseChargesFromCIFFile no
UnitCells 3 3 3
HeliumVoidFraction  1
ExternalTemperature 358
ExternalPressure 3e7 4e7
Movies yes
WriteMoviesEvery  5000
ComputeDensityProfile3DVTKGrid           yes
WriteDensityProfile3DVTKGridEvery        1000
DensityProfile3DVTKGridPoints            150 150 150
AverageDensityOverUniteCellsVTK          yes
DensityAverageTypeVTK                    FullBox

Component 0 MoleculeName              C22
            MoleculeDefinition        TraPPE
            TranslationProbability    1.0
            RotationProbability       1.0
            ReinsertionProbability    1.0
            IdentityChangeProbability 0.0
            SwapProbability           1.0
            CreateNumberOfMolecules   0


***************
The force field parameter:

# general rule for shifted vs truncated
shifted
# general rule tailcorrections
no
# number of defined interactions
4
# type interaction, parameters.    IMPORTANT: define shortest matches first, so that more specific ones overwrites these
O_             lennard-jones    78.263        3.5532
Al_            lennard-Jones    0.00067       4.7943
Si_            lennard-Jones    0.000926      3.7064
Ca_            lennard-Jones    50.362        3.2237

# general mixing rule for Lennard-Jones
Lorentz-Berthelot