News:

SMF - Just Installed!

Main Menu

Recent posts

#61
Visualization / Density Plot
Last post by saif - August 21, 2023, 06:49:01 PM
Hi,

I'm seeking advice.

I'm attempting to plot the DensityProfile.vtk file that I obtained from my simulation run. However, I'm encountering an issue where the plot appears completely white without any visible content.
All the data inside the file consists of very small values.

-----------------------------------------------------------
Part of my file:

# vtk DataFile Version 1.0
CELL_PARAMETERS 19.152000 19.152000 19.152000 90.000000 90.000000 90.000000
ASCII
DATASET STRUCTURED_POINTS
DIMENSIONS 150 150 150
ORIGIN 0.000000 0.000000 0.000000
SPACING 0.127680 0.127680 0.127680
POINT_DATA 3375000
SCALARS scalars double
LOOKUP_TABLE default
0
0.000584217
0.00084387
0.00183199
0.00266144
0.00421935
0.00569792
0.00809249
0.00786169
0.0121387

-----------------------------------------------------------
Thank you
#62
General / Density Histogram from GCMC si...
Last post by Cigdem - August 16, 2023, 08:32:36 AM
Dear Prof. Dubbeldam and dear RASPA Community,

With adding these lines to the simulation input, we can obtain density profiles as a result of the GCMC simulation:

ComputeDensityProfile3DVTKGrid yes
WriteDensityProfile3DVTKGridEvery 1000
DensityProfile3DVTKGridPoints 150 150 150
AverageDensityOverUnitCellsVTK yes
DensityAveragingTypeVTK FullBox
FreeEnergyAveragingTypeVTK FullCell


However, we encounter a problem when we want to run GCMC simulations to obtain density profiles from four-component gas mixture adsorption in a MOF.
The error message is long but starts and ends as given below:

" *** Error in `simulate': free(): invalid next size (normal): 0x000056087ce427e0 ***

......
....
..
.

/var/spool/slurm/d/job5763373/slurm_script: line 67: 14511 Aborted                 $COMMAND $INPUT

RASPA finished. Return code is 134".

I tried with different cif files and different gas definitions, but, as far as I note, it only happens when I run the simulation for a mixture of more than three components.
I would be glad to learn if there is something that I am missing and avoids the production of density profiles for a mixture of more than three components.

Thank you,
#63
Input files and parameters / RDF for adsorption
Last post by mhaeri - July 22, 2023, 06:31:25 PM
Hi Professor Dubeldam,
I have a question about the implementation of rdf:
The RASPA documents contain a clear description of rdf calculation for pure fluids.
How can I get rdf for systems involving the adsorbent-adsorbate?
I know that rdf for adsorbate-adsorbent, as well as adsorbate-adsorbate, are routinely reported in the literature using RASPA.
But I am not sure how it is implemented. It looks like it is more than just replacing the box with the framework. Could you kindly provide me with a sample simulation.input script ?
Best regards,
Mas
#64
General / Re: Can the cavity of a MOF be...
Last post by Lingxiao Qin - June 27, 2023, 01:40:37 PM
Quote from: David Dubbeldam on June 13, 2023, 01:45:24 PM
In the gallery, then primitives, then IRMOF-1, you can see such an example.

Thank you for your reply. I do find the example. However, I am struggling to reproduce such a picture. The illustrations in the gallery are very compelling, but I think it may be difficult for a beginner like me to figure out how to draw those pictures without a tutorial. I was wondering if it would be possible for you to write a tutorial about how to draw the pictures in the gallery in the future. I understand that creating tutorials requires time and effort, but I believe this tutorial would be invaluable to the community.
#65
General / Re: Inquiry regarding reconven...
Last post by kokojiang - June 24, 2023, 08:46:18 AM
Hi, actually I met with the same problem. I still dont know how to solve this yet. I tried this to make it continue, but I DONT know if it is correct:

1.Copy the whole simulation files, including input,output...
2.Rename  "Restart" folder to "RestartInitial"
3.Add "RestartFile  yes" to the simulation.input file and run

This will start new simulation cycles after the last results.

In the Raspa manual, it says "ContinueAfterCrash yes" will work, however, I tried it many times and it won't. My senior said he used to run it successfully, but now he can't either. Maybe this is becuase we uploaded Raspa to the latest version?

I also tried to interrupt some still running simulations, change the input file, increase the number of cycles, add "ContinueAfterCrash yes", and continue. After this, I checked the output file, it still continue after the original inputfile, not the one I changed. For example, the origin inputfile:"NumberOfCycles   100000", changed to "NumberOfCycles   200000", it still finished at 100000 sycles.

I am still finding a way that can continue the simulation in the same outputfile directory.
#66
Input files and parameters / Re: using grids for energy bia...
Last post by wei123 - June 22, 2023, 09:05:35 AM
Yes, the adsorption using the grid was much higher than without using the grid (about 10 times more already within a few thousand MC steps until I stopped the simulation, and the adsorption increases continuously as the simulation runs). I assume it is because the framework was charged, because the energy grid was created for the framework without the presence of cations.
#67
General / Re: RxMC ensemble
Last post by Mohammad Ahmadi - June 18, 2023, 08:47:52 AM
I mean the sentence "how can the values of TranslationProbability, RotationProbability, ReinsertionProbability be calculated?"
  That is, is there a formula to calculate these values or do they have to be guessed؟
#68
General / Re: RxMC ensemble
Last post by Mohammad Ahmadi - June 18, 2023, 08:24:14 AM
The meaning of fluctuating behavior is that, for example, in a temperature range, the results obtained from the conversion behave like a sine function and do not have a specific trend, sometimes they increase and sometimes decrease.
#69
General / Re: Atomization energies in pa...
Last post by David Dubbeldam - June 13, 2023, 02:11:06 PM
In RASPA2 you input the final result: the partition-factor. In more recent RASPA versions the Log of the partition function to avoid over/underflow.
It is an user input, and up to the user to compute it.
#70
Input files and parameters / Re: pseudo_atom radii
Last post by David Dubbeldam - June 13, 2023, 02:09:19 PM
You can look these radii up in handbooks.
SMF spam blocked by CleanTalk