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

#91
You are correct (and the manual has an error). You can put a list of pressures in the input, but only one temperature.
In practice, using a list of pressures is a bad idea, unless you really want to compute adsorption+desorption and examine hysteresis. The downside is that all simulations have to be run sequentially.

In practice, people develop scripts to submit batch-scripts to submit large number of jobs (for example: for various MOFs, various temperatures, various molecules, compute the isotherm over a pressure range).

An an inspiration, see the example scrips in the 'scripts' directory of the RASPA source.
You fill in the top-part, run the script and it will generate the input-files. A second script will read the data out of the output-files and put in a more useful format that can be used for plotting.
#92
FIXED_BOND is for fixing a single bond, while the full molecule can still be flexible via bend and torsion potentials. An example is the TraPPE model for united-atom alkanes. In Monte Carlo it is easy to fix bonds because the bonds are _generated_. In Molecular Dynamics algorithms like shake needs to be used.
I would say that is better not to use fixed bonds, but rather use a harmonic bond-potential.
#93
I would think they should be identical "in principle", if the force field (i.e. bond, bend, torsion parameters) is the same for both. So the molecule has the same type of atoms, same structure except R vs S.
Like you say, in practice they differ slightly, but the difference should vanish in the limit of running extremely long.
#94
Your molecule has not diffused at all. I would first do a test with a small molecule so that you see what a graph looks like for a molecule that diffuses.
The y-axis is the MSD, in Angstrom squared. You can get diffusion results when the molecules diffuse at least through the repeating unit (so unit-cell-width squared).
In this case, MD is probably not capable of computing diffusion because it is too slow in your system.
#95
See the many articles of Martin and Siepmann on this topic.

To compute pressure:
ComputeMolecularPressure yes
#96
You can use a tool like "openbabel" to convert between different file formats.
#97
See the examples for molecules. RIGID_BOND is used for molecules that are fully rigid, for example water or CO2. The bonding information is used for the CBMC-biasing.
RIGID_BOND is a type of bond-potential that has only one bond-distance.
#98
MIL-101 is very large indeed. RASPA is not optimized for speed.  You can use a debugger to see where it goes wrong.
#99
General / Re: High-Pressure Simulations
July 21, 2021, 08:51:45 AM
9% at those high pressure is still very, very good I would say. The validity of such a simple methane model is not "endless". Any classical model will show weaknesses at extreme conditions. Perhaps a full-atom model might do better, but these are harder to paremeterize.
#100
1) As many cycles and molecules as you need. The more the better, but the slower the simulations. That's always a trade-off.
In equilibrium, the macroscopic properties no longer change. The amount of cycles needed to reach equilibrium depends on the system.
To limit finite size effects you need a large system. In practice you can simulate various system sizes as see when your properties converge to a plateau value as a function of system size.

2) in a NVT or muVT or Gibbs NVT simulation the energy fluctuates. This energy fluctuation is related to a physical property called the heat capacity.

3) There are no rules of thumbs, but rather it will take experience to get a feel for this.
#101
Yes, all Qeq are differently implemented. The code in RASPA is based on the method of Wilmer and Snurr.
https://pubs.acs.org/doi/10.1021/jz3008485
#102
General / Re: parameter input
June 08, 2021, 02:36:34 PM
The 'pseudo_atoms.def' contains predefined atoms, and this is often the way to go for molecules.
This especially important for united atom types etc.

For frameworks, read from a CIF file, there is always an element defined (_atom_site_type_symbol), and the type (_atom_site_label) is automatically added when found. The mass is a lookup from a database in 'src/scattering_factors.c'.
#103
Input files and parameters / Re: Density histograms
June 08, 2021, 02:28:59 PM
This function writes out a VTK-file, a 3D histograms of the density. See the examples on Visualization.
You could also read this file yourself.
#104
The DREIDING force field uses specific functional forms and parameters. You can look these up and use (most of) these in RASPA.
#105
Input looks, at first sight, okay to me. Remove as much as possible, i.e. component 1-3 etc, and check if it works. Then from add the components one by one to see where it goes wrong.
In general, change stuff in small steps, otherwise it is harder to debug.