- Welcome to iRASPA Community.
Recent posts
#51
Input files and parameters / Bulk Density Validation by MC
Last post by Bourne - August 07, 2024, 04:04:41 AMDear Dr. Dubbeldam and RASPA users,
I am a beginner with the RASPA software and am conducting a study on gas occurrence behaviors within porous rock. As part of my study, I am performing validation work on the fluid model by calculating the bulk density within an empty box at 333K using MC simulations. I have tested the TraPPE and EPM2 models. However, the results are unsatisfactory, as both bulk density values deviate from the data provided by NIST and the P-R equation, despite the systems being well equilibrated. This discrepancy is confusing, and I've tried various approaches without success. I have checked all the input files, and they appear to be correct. Could you please advise on how I might resolve this issue?
Thank you for your time.
Best regards,
Input files for TraPPE Model Validation:
SimulationType MonteCarlo
NumberOfCycles 20000
NumberOfInitializationCycles 6000
PrintEvery 1000
RestartFile no
ChargeMethod Ewald
EwaldPrecision 1e-5
Forcefield TraPPE
RemoveAtomNumberCodeFromLabel no
Box 0
BoxLengths 25 25 25
ExternalTemperature 333.0
ExternalPressure 5e6 10e6 15e6 20e6 25e6 30e6
Movies yes
WriteMoviesEvery 2000
Component 0 MoleculeName CO2
MoleculeDefinition TraPPE
TranslationProbability 0.5
RotationProbability 0.5
ReinsertionProbability 0.5
SwapProbability 1.0
CreateNumberOfMolecules 0
I am a beginner with the RASPA software and am conducting a study on gas occurrence behaviors within porous rock. As part of my study, I am performing validation work on the fluid model by calculating the bulk density within an empty box at 333K using MC simulations. I have tested the TraPPE and EPM2 models. However, the results are unsatisfactory, as both bulk density values deviate from the data provided by NIST and the P-R equation, despite the systems being well equilibrated. This discrepancy is confusing, and I've tried various approaches without success. I have checked all the input files, and they appear to be correct. Could you please advise on how I might resolve this issue?
Thank you for your time.
Best regards,
Input files for TraPPE Model Validation:
SimulationType MonteCarlo
NumberOfCycles 20000
NumberOfInitializationCycles 6000
PrintEvery 1000
RestartFile no
ChargeMethod Ewald
EwaldPrecision 1e-5
Forcefield TraPPE
RemoveAtomNumberCodeFromLabel no
Box 0
BoxLengths 25 25 25
ExternalTemperature 333.0
ExternalPressure 5e6 10e6 15e6 20e6 25e6 30e6
Movies yes
WriteMoviesEvery 2000
Component 0 MoleculeName CO2
MoleculeDefinition TraPPE
TranslationProbability 0.5
RotationProbability 0.5
ReinsertionProbability 0.5
SwapProbability 1.0
CreateNumberOfMolecules 0
#52
Output files / Calculation of MSD from GCMC a...
Last post by thor - August 03, 2024, 05:39:00 AMDear Professor David:
I want to use the final configuration of GCMC adsorption as the initial configuration for MD simulation to calculate MSD. So I renamed the Restart file obtained by adsorption to RestartInitial and used the following file as input. After the simulation is completed, I did not get the corresponding MSDOrderN. Why is this? Is my input file setting correct?
SimulationType MolecularDynamics
NumberOfCycles 1000000
NumberOfInitializationCycles 0
NumberOfEquilibrationCycles 1000000
PrintEvery 5000
PrintPropertiesEvery 1000
RestartFile yes
Ensemble NVT
TimeStep 0.0005
ComputeMSD yes
PrintMSDEvery 5000
Forcefield UFF-Direding
UseChargesFromCIFFile yes
ChargeMethod Ewald
Framework 0
FrameworkName COF2207-8
CutOffVDW 12.8
UnitCells 1 1 1
ExternalTemperature 298
ExternalPressure 1e5
Component 0 MoleculeName CO2
MoleculeDefinition TraPPE
TranslationProbability 0.5
RotationProbability 0.5
ReinsertionProbability 0.5
SwapProbability 1.0
CreateNumberOfMolecules 0
I want to use the final configuration of GCMC adsorption as the initial configuration for MD simulation to calculate MSD. So I renamed the Restart file obtained by adsorption to RestartInitial and used the following file as input. After the simulation is completed, I did not get the corresponding MSDOrderN. Why is this? Is my input file setting correct?
SimulationType MolecularDynamics
NumberOfCycles 1000000
NumberOfInitializationCycles 0
NumberOfEquilibrationCycles 1000000
PrintEvery 5000
PrintPropertiesEvery 1000
RestartFile yes
Ensemble NVT
TimeStep 0.0005
ComputeMSD yes
PrintMSDEvery 5000
Forcefield UFF-Direding
UseChargesFromCIFFile yes
ChargeMethod Ewald
Framework 0
FrameworkName COF2207-8
CutOffVDW 12.8
UnitCells 1 1 1
ExternalTemperature 298
ExternalPressure 1e5
Component 0 MoleculeName CO2
MoleculeDefinition TraPPE
TranslationProbability 0.5
RotationProbability 0.5
ReinsertionProbability 0.5
SwapProbability 1.0
CreateNumberOfMolecules 0
#53
General / iRASPA Potential Energy Surfac...
Last post by Cecilia - July 23, 2024, 11:46:59 AMHello,
I was wondering how I can use iRASPA to generate a potential energy surface a MOF. I had tried downloaded iRASPA from the Microsoft store and the Apple store but neither seems to have the "Adsorption Surface" section under the "Appearance" tab that I saw in the photos.
Furthermore, I was wondering if it is possible to get the energies of each molecule per MC step from RASPA simulations with the grand canonical ensemble?
Many thanks in advance!
I was wondering how I can use iRASPA to generate a potential energy surface a MOF. I had tried downloaded iRASPA from the Microsoft store and the Apple store but neither seems to have the "Adsorption Surface" section under the "Appearance" tab that I saw in the photos.
Furthermore, I was wondering if it is possible to get the energies of each molecule per MC step from RASPA simulations with the grand canonical ensemble?
Many thanks in advance!
#54
Visualization / Adsorption density distributio...
Last post by Z.GAO - July 03, 2024, 05:25:02 AMHi everyone,
In the post-processing of gas adsorption MOF, I want to make an adsorption density distribution map. MOF is a triclinic system, and the adsorption density map drawn by VisIt software is still a cubic structure. So I use iRASPA software to draw it. Although iRASPA can show the triclinic structure, there is no color scale, only one color, which makes it difficult to intuitively judge the priority adsorption site. How can I solve this problem?
In addition, when I use iRASPA to present the molecular structure, I find that some atoms are not bonded. How should I deal with this?
Best wishes!
Zheng
In the post-processing of gas adsorption MOF, I want to make an adsorption density distribution map. MOF is a triclinic system, and the adsorption density map drawn by VisIt software is still a cubic structure. So I use iRASPA software to draw it. Although iRASPA can show the triclinic structure, there is no color scale, only one color, which makes it difficult to intuitively judge the priority adsorption site. How can I solve this problem?
In addition, when I use iRASPA to present the molecular structure, I find that some atoms are not bonded. How should I deal with this?
Best wishes!
Zheng
#55
General / Difficult to install RASPA3 on...
Last post by Yu_Zhuochen - June 20, 2024, 12:10:33 PMI've run these following commands,and try to install raspa3 on ubuntu:
$sudo apt install build-essential git cmake
$sudo add-apt-repository 'deb http://apt.llvm.org/jammy/ llvm-toolchain-jammy-18 main'
$sudo add-apt-repository 'deb-src http://apt.llvm.org/jammy/ llvm-toolchain-jammy-18 main'
$sudo apt install libllvm-18-ocaml-dev libllvm18 llvm-18 llvm-18-dev llvm-18-runtime
$sudo apt install libc++-18-dev libc++abi-18-dev libomp-18-dev pybind11-dev python3-pybind11 liblapack-dev clang-tools-18
$cd src
$vi makefile-manual and set CXX=clang++-18
$make -f makefile-manual
However, there are several errors during installation:
In module 'std_atomic' imported from foundationkit/threadpool.ixx:31:
In module 'std_private_atomic_aliases' imported from /usr/lib/llvm-18/bin/../include/c++/v1/atomic:591:
In module 'std_cstdlib' imported from /usr/lib/llvm-18/bin/../include/c++/v1/__atomic/aliases.h:21:
/usr/lib/llvm-18/bin/../include/c++/v1/stdlib.h:123:59: error: 'abs' has different definitions in different modules; definition in module 'std_cstdlib' first difference is function body
123 | _LIBCPP_NODISCARD_EXT inline _LIBCPP_HIDE_FROM_ABI double abs(double __lcpp_x) _NOEXCEPT {
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
124 | return __builtin_fabs(__lcpp_x);
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
125 | }
| ~
/usr/lib/llvm-18/bin/../include/c++/v1/stdlib.h:123:59: note: but in '' found a different body
123 | _LIBCPP_NODISCARD_EXT inline _LIBCPP_HIDE_FROM_ABI double abs(double __lcpp_x) _NOEXCEPT {
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
124 | return __builtin_fabs(__lcpp_x);
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
125 | }
| ~
fatal error: too many errors emitted, stopping now [-ferror-limit=]
20 errors generated.
make: *** [makefile-manual:47: foundationkit/threadpool.pcm] Error 1
How can I solve this problem?
$sudo apt install build-essential git cmake
$sudo add-apt-repository 'deb http://apt.llvm.org/jammy/ llvm-toolchain-jammy-18 main'
$sudo add-apt-repository 'deb-src http://apt.llvm.org/jammy/ llvm-toolchain-jammy-18 main'
$sudo apt install libllvm-18-ocaml-dev libllvm18 llvm-18 llvm-18-dev llvm-18-runtime
$sudo apt install libc++-18-dev libc++abi-18-dev libomp-18-dev pybind11-dev python3-pybind11 liblapack-dev clang-tools-18
$cd src
$vi makefile-manual and set CXX=clang++-18
$make -f makefile-manual
However, there are several errors during installation:
In module 'std_atomic' imported from foundationkit/threadpool.ixx:31:
In module 'std_private_atomic_aliases' imported from /usr/lib/llvm-18/bin/../include/c++/v1/atomic:591:
In module 'std_cstdlib' imported from /usr/lib/llvm-18/bin/../include/c++/v1/__atomic/aliases.h:21:
/usr/lib/llvm-18/bin/../include/c++/v1/stdlib.h:123:59: error: 'abs' has different definitions in different modules; definition in module 'std_cstdlib' first difference is function body
123 | _LIBCPP_NODISCARD_EXT inline _LIBCPP_HIDE_FROM_ABI double abs(double __lcpp_x) _NOEXCEPT {
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
124 | return __builtin_fabs(__lcpp_x);
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
125 | }
| ~
/usr/lib/llvm-18/bin/../include/c++/v1/stdlib.h:123:59: note: but in '' found a different body
123 | _LIBCPP_NODISCARD_EXT inline _LIBCPP_HIDE_FROM_ABI double abs(double __lcpp_x) _NOEXCEPT {
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
124 | return __builtin_fabs(__lcpp_x);
| ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
125 | }
| ~
fatal error: too many errors emitted, stopping now [-ferror-limit=]
20 errors generated.
make: *** [makefile-manual:47: foundationkit/threadpool.pcm] Error 1
How can I solve this problem?
#56
General / Re: iRASPA cannot display the ...
Last post by Lingxiao Qin - June 10, 2024, 10:23:58 AMI have tried your input file, and both the 'DensityProfile_methane.vtk' and 'COMDensityProfile_methane.vtk' output files can be successfully imported into iRASPA. I suggest checking the versions of RASPA and iRASPA that you have installed. The latest version of RASPA is v2.0.47, so I doubt whether you really installed RASPA v2.0.4 (released in 2016) or it is just a typo. The VTK files generated by such an old version of RASPA might not be readable by iRASPA. I am not sure of the exact version of RASPA I am using, but it was at least released after 2022. I downloaded iRASPA from the Windows store, which was released in 2018. By the way, the lines
Code Select
Ensemble NVT
TimeStep 0.0005 in your input file seem unusual. From my understanding, these options are used for molecular dynamics simulations and should not be used in a Grand Canonical Monte Carlo simulation, so I commented these two lines when testing your input file. #57
I am trying to use the MM3_VDW but when I calculate the vDW energy for 2 atoms only using the same equation in the manual I am getting different value from the printed one in the output file
#58
General / Single point calculation for a...
Last post by Sara - June 04, 2024, 04:50:09 PMHi Dr David Dubbeldam,
Is there a way to compute the adsorbate-adsorbate energy for a defined XYZ coordinates for the both adsorbates? I mean I tried to put them in a box but raspa keeps changing the positions, I want it to read the exact positions in the .def file. Is that possible? it is like a single point calculation.
Is there a way to compute the adsorbate-adsorbate energy for a defined XYZ coordinates for the both adsorbates? I mean I tried to put them in a box but raspa keeps changing the positions, I want it to read the exact positions in the .def file. Is that possible? it is like a single point calculation.
#59
General / Re: How to make adsorption den...
Last post by John Walker - May 29, 2024, 10:14:56 AM大佬，你会了吗？
#60
General / iRASPA cannot display the dens...
Last post by John Walker - May 29, 2024, 09:16:10 AMDear:
I have recently been using MOF materials for methane adsorption, and have attempted to import VTK files generated by RASPA2.0.4 into the IRASPA software. I am only able to generate MOF structure diagrams, but not the methane diagrams. I imported using 'DensityProfile_methane' and 'COMDensityProfile_methane' from the VTK folder, and received the error 'Import failed, Cannot read VTK header (line: 910 file: iraspagui/mainwindow.cpp)'. Importing 'Framework_final.pdb' from the Movies folder does display the molecular structure, which confuses me. How can I produce a methane adsorption density map? I can do it with VisIt.
SimulationType MonteCarlo
NumberOfCycles 30000
NumberOfInitializationCycles 5000
PrintEvery 1000
ContinueAfterCrash no
#WriteBinaryRestartFileEvery 1000
ChargeMethod Ewald
EwaldPrecision 1e-6
Ensemble NVT
TimeStep 0.0005
Forcefield GenericMOFs
CutOffVDW 12
RemoveAtomNumberCodeFromLabel yes
Movies yes
WriteMoviesEvery 1000
RemoveAtomNumberCodeFromLabel yes
UseChargesFromCIFFile yes
Framework 0
FrameworkName HKUST-1
UnitCells 2 2 2
HeliumVoidFraction 0.39778
ExternalTemperature 298.15
ExternalPressure 1.0e+06
ComputeDensityProfile3DVTKGrid yes
WriteDensityProfile3DVTKGridEvery 1000
DensityProfile3DVTKGridPoints 150 150 150
AverageDensityOverUnitCellsVTK yes
DensityAveragingTypeVTK FullBox
Component 0 MoleculeName methane
MoleculeDefinition TraPPE
IdealGasRosenbluthWeight 1.0
TranslationProbability 1.0
RotationProbability 1.0
ReinsertionProbability 1.0
##CBMCProbability 1.0
SwapProbability 1.0
CreateNumberOfMolecules 0
I have recently been using MOF materials for methane adsorption, and have attempted to import VTK files generated by RASPA2.0.4 into the IRASPA software. I am only able to generate MOF structure diagrams, but not the methane diagrams. I imported using 'DensityProfile_methane' and 'COMDensityProfile_methane' from the VTK folder, and received the error 'Import failed, Cannot read VTK header (line: 910 file: iraspagui/mainwindow.cpp)'. Importing 'Framework_final.pdb' from the Movies folder does display the molecular structure, which confuses me. How can I produce a methane adsorption density map? I can do it with VisIt.
SimulationType MonteCarlo
NumberOfCycles 30000
NumberOfInitializationCycles 5000
PrintEvery 1000
ContinueAfterCrash no
#WriteBinaryRestartFileEvery 1000
ChargeMethod Ewald
EwaldPrecision 1e-6
Ensemble NVT
TimeStep 0.0005
Forcefield GenericMOFs
CutOffVDW 12
RemoveAtomNumberCodeFromLabel yes
Movies yes
WriteMoviesEvery 1000
RemoveAtomNumberCodeFromLabel yes
UseChargesFromCIFFile yes
Framework 0
FrameworkName HKUST-1
UnitCells 2 2 2
HeliumVoidFraction 0.39778
ExternalTemperature 298.15
ExternalPressure 1.0e+06
ComputeDensityProfile3DVTKGrid yes
WriteDensityProfile3DVTKGridEvery 1000
DensityProfile3DVTKGridPoints 150 150 150
AverageDensityOverUnitCellsVTK yes
DensityAveragingTypeVTK FullBox
Component 0 MoleculeName methane
MoleculeDefinition TraPPE
IdealGasRosenbluthWeight 1.0
TranslationProbability 1.0
RotationProbability 1.0
ReinsertionProbability 1.0
##CBMCProbability 1.0
SwapProbability 1.0
CreateNumberOfMolecules 0
SMF spam blocked by CleanTalk