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.

Messages - David Dubbeldam

Pages: [1] 2 3 ... 15
Visualization / Re: RASPA visualization (Density Plot)
« on: September 25, 2022, 06:01:19 PM »
That depends on the VTK-version. But VTK-file plotting is now included in iRASPA, provided you have written the VTK with the latest version of RASPA.

General / Re: Phase state in MC simulation of grand canonical ensemble
« on: September 25, 2022, 05:58:14 PM »
Based on a Peng-Robinson equation of state it calculates the phase and density of the pure fluid (without framework). This density is used to calculate the the excess adsorption.

Input files and parameters / Re: desorption-error
« on: September 25, 2022, 05:55:10 PM »
You can use a debugger to find out where it crashes and then what input error caused it.

Input files and parameters / Re: desorption-example
« on: September 25, 2022, 05:51:53 PM »
Desorption is done by going from a high (saturation) pressure back to low pressure.

Visualization / Re: TIP4P visualization and bond lengths
« on: September 25, 2022, 05:47:07 PM »
You can select whether an atom is visible or not (checkbox).

The bond-distance criteria is a standard one. This is not yet easily changeable.

Input files and parameters / Re: Problem in adsorption
« on: September 25, 2022, 05:42:33 PM »
I would suggest to talk to experimental people and ask why their data goes up.
I expect this is due to adsorption at the external surface. You would expect that the results for adsorption should level off to a saturation value.

General / Re: _atom_site_disorder_group in .cif file
« on: September 25, 2022, 05:36:46 PM »
There are no standard methods. By making choices you need to create a structure with no symmetry (P1) so that all atom positions are defined.

General / Re: simulation of reaction
« on: September 25, 2022, 05:33:15 PM »
Instead of a Box use a Framework.

Input files and parameters / Re: Idealgasrosenbluchweight
« on: July 14, 2022, 07:43:34 PM »
The ideal gas Rosenbluth weight is only depending on temperature. For each temperature you can pre-compute it (see examples).

The 'CutOffVDW' should not be different for the gas and liquid phase, otherwise you're using different force fields for them.
The 'CutOffCoulomb' should be set to something slightly smaller than half the box size. This significantly speeds up the computation for the gas-phase since the gas-phase box can be large.
Importantly, both can be set per system (just list them after the 'Box' statement).

General / Re: list of configurational moves in .def file
« on: June 22, 2022, 10:13:44 PM »
This for the PartialReinsertionProbability move, where only a few atoms are regrown and the others are left fixed.
So, in this case there are 4 config move defined (and you see 5 lines). For each line you have the number of fixed atoms, and the the list of the fixed atom ids.
The other atoms are regrown using CBMC. This move has high efficiency, because the molecule is already present and is only partially regrown.
Note that in CBMC all branches need to be grown at the same time, so you can not keep more than one branch fixed (in a branched molecule).

General / Re: How to get movie from MD
« on: June 14, 2022, 09:33:53 PM »
An "MD path" is the path or trajectory that a particle follows in MD. if you ran too short, then the molecules have not moved enough.  Diffusion is defined as the limiting behavior of the Mean-Squared Displacement (MSD) at long-times. For diffusional processes, the slope must become linear at long times. If you did not reach that regime yet, then you have ran too short.
Slow diffusion might be outside the reach of MD, e.g. very slow diffusion can not be measured using MD. The maximum is probability up to 10^-14 m^2/s.

General / Re: How to get movie from MD
« on: June 09, 2022, 01:11:09 PM »
In general no, because diffusion coefficients are by definition the long-time behavior of molecules.
In any short MD-path you might even see opposite behavior. Of course for illustrative purposes you could just select an MD path that shows what you want to show.

These radii are only needed for flexible frameworks, and used as bond-criteria to determine bonding.
For other cases these values do not matter.

General / Re:negative MSD slope
« on: June 04, 2022, 07:06:06 PM »
See literature on how to compute diffusion coefficients from MSD. The slops of the MSD at long times can not be negative.
You're probably looking at a situation where MD is not capable of computing diffusion because it is too slow.

Pages: [1] 2 3 ... 15