News:

SMF - Just Installed!

Main Menu
Menu

Show posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.

Show posts Menu

Messages - David Dubbeldam

#181
Output files / Re: Histograms produced by raspa
August 31, 2019, 07:42:21 PM
Histograms do not really have units on the y-axis. The x-axis, for energy, would be Kelvin.
You can easily figure these things out by checking the source, and check what is exactly counted and how it is written out.
See source file: sample.c, and the function 'void SampleEnergyHistogram(int Switch)'.
#182
Saving of snapshots is more easily done by adding them to a movie. But you could run short simulations and copy the restart-file saved at the end to another directory if the goal is to create many restart-files.
#183
That is not possible. The atoms for the molecule file needs to be defined and the charge is taken from the pseudo_atoms.def file.
So, currently, the same type of atom of a molecule must have the same properties.
#184
You could add your way of computing the pressure in the source code, or modify the implemented calculation of the pressure.
#185
To reliably model adsorption/diffusion, the mobility of the cations needs to be taken into account. For example, in FAU there are many different cation-sites with certain "occupancy" ratios. Your cation model needs to be able to reproduce that.
That is why it is hard to reach an equilibrated system and it is a good idea to, when you have generated that, to store these positions in a restart-file.
Since cations are charged, the usual CBMC biasing can lead to numerical issues, but for single-atom cations you can achieve the same with a random-translation that tries to place the cation randomly in the box.
Look in the literature for papers on modeling of zeolites with cations for more pointers.
#186
What do you mean with "there were no more than 2 number of adsorbent in each cycle", do you have multiple frameworks? do you mean "adsorbent" which is the framework, or "adsorbate" which means the molecules inside the framework.
"So I think something wrong about the water, maybe the increased charge of it in the MC simulation.", have you tried making it charge-neutral?
#187
Are you comparing your simulation to other simulations or your simulations to experiments?
For the latter you need a bigger picture view and compare the simulated isotherm to the experimental one. Usually at some pressure the isotherm shoots up due to water-water clustering, but the details of where that happens depends on many factors (especially defects).
Also, the water-model is never perfect, nor is the force field in general.
#188
General / Re: Using RASPA for Nanocarbons
July 06, 2019, 04:36:59 PM
Just use a cif-file describing your carbon-nanotube.
#189
Visualization / Re: VTK visualization
July 06, 2019, 04:34:44 PM
For combining VTK-densities and other VTK-files use the appropriate number of unit cells.
#190
The simulation-length is set by specifying the NumberOfCycles (see the manual).
#191
The FH-LJ is a potential between a pair of atoms. There are no mixing rules for atoms that defined using different potential forms. But you can straightforwardly define the interactions of all pairs of atoms that you would like (see the manual).
#192
Use the restart-file to start from the equilibrated systems (see manual).
You can take the results from different simulations and average them.
#193
Input files and parameters / Re: CIF files help
July 06, 2019, 04:24:19 PM
P1 is safer. Specifying the spacegroup is not unique, unless the Hall-symbols are used. Materials Studio and other codes do not do that. So, yes you can use symmetry, but if it does not work, it is usually because the space group is not set properly.
Using P1 avoids these issues.
#194
You could change the source code to achieve that.
Look at the function 'int VolumeMove(void)' in mc_moves.c and change the isotropic scaling into only scaling in the z-direction.
#195
In RASPA, the ensemble follows from the MC-moves.
So for NPT Gibbs you would use:
1) VolumeChangeProbability
2) GibbsSwapProbability
3) GibbsIdentityChangeProbability (assuming you have more than one component, otherwise you would use NVT-Gibbs).
4) thermalization moves like TranslationProbability, RotationProbability, ReinsertionProbability etc.

So the difference between NVT-Gibbs and NPT-Gibbs is that NVT-Gibbs uses a special volume move 'GibbsVolumeChangeProbability' so that the total volume is fixed. NPT-Gibbs uses the normal volume 'VolumeChangeProbability' that adjust each simulation box according to the specified pressure.