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 ... 6
Yes, the first component is then not changing it's number of molecules, while the second component is swapped in and out.

General / Re: error run simulate
« on: March 01, 2020, 09:40:13 AM »
Since I see some python stuff, so it is probably that you are using the numat-code (
You could ask them about their additional python routines.
Our version is strictly C, and is at:

Visualization / Re: Visualisation using Para View
« on: February 03, 2020, 09:32:40 PM »
I guess you need to figure out how Paraview scales these. For the density file, it is by default 150x150x150 points, and this needs to be converted to the real unit cell.

Visualization / Re: Imported PDB file
« on: February 03, 2020, 09:25:16 PM »
There are many pdb-files out there that are not according to the PDB specification. You could try to use "openbabel" to create a good pdb-file and see if that works.

Chemical potentials and fugacities are related by a trivial conversion and represent the same concept.
Fugacity is easier to use because it is always a positive number, while the chemical potential goes to minus infinity at infinite dilution.
Fugacity also has a more direct relation to pressure.

Simulation algorithms and theory / Re: about CFCMC potential
« on: January 16, 2020, 06:17:31 PM »
Energies are capped at high values to avoid numerical artifacts. In CFCMC, you can also use a small distance overlap criteria (as in the original CFCMC work of the Maginn-group).

You need a good force field to reproduce the correct occupancy. Also, the simulations are non-trivial and some people use parallel-tempering to compute the proper occupancy distributions.

General / Re: Flexible-MOFs forcefield Issue
« on: January 16, 2020, 06:13:10 PM »
The interactions are up to you. Bond-potentials and bend-potentials between the atoms are common, but force fields differ in the treatment of torsion-and improper torsion potentials, and VDW and electrostatics.

Output files / Re: Visualization of adsorption
« on: December 09, 2019, 08:54:35 AM »
You just import all the pdb- and cif-files you would like to visualize (so select multiple files during import).
Then you can set individual visualization properties, like origin, drawing-method, etc.

Either use a model that puts the charge-centers on the atomic centers (OPLS) or if the molecule is small you can also make the whole molecule rigid.

General / Re: Flexible-MOFs forcefield Issue
« on: November 08, 2019, 04:13:47 PM »
For a force field for a flexible MOF you need to define all the interactions similar to the flexible IRMOF-1 example.
The GenericMOFs defines only the interactions of molecules with the framework.

Simulation algorithms and theory / Re: MaxNumberOfCoulombicSites
« on: October 04, 2019, 01:18:16 PM »
You could do some basic debugging, like setting the PrintEvery to 1 and see if it is just slow or that you are really stuck at step 0.
The latter can happen if you try to put in more molecules than fit in the simulation volume.

Visualization / Re: Making the grid finer for adsorption surfaces
« on: October 03, 2019, 11:56:24 AM »
Currently it is using a grid of 128x128x128. Several opties are considered in future releases:
1) an option for 256x256x256
2) Use tessellation on newer videocards to create smoother surfaces.

Announcements / New review on force fields
« on: September 30, 2019, 06:41:57 PM »
An open access review that potentially could be of interest to some members on the forum, especially those that are interested in force fields, force field types and design, implementation (gradients, second derivatives, strain derivatives), Ewald summation, polarization, optimization, parameterization, Machine Learning, and General-Purpose GPU computing:

Simulation algorithms and theory / Re: MaxNumberOfCoulombicSites
« on: September 29, 2019, 11:45:36 AM »
In open ensembles (Gibbs, grand-canonical ensemble) the number of molecules fluctuates. So if more memory is needed the code needs to realloc the memory to make it larger.

Pages: [1] 2 3 ... 6