1
Output files / Re: Energy landscape
« on: April 18, 2022, 07:56:22 PM »
Recently, these capabilities are implemented in iRASPA, where you play around with these energy-surfaces on the fly (computed very fast on the GPU).
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.
2
Simulation algorithms and theory / Re: Reduction of calculation effort by usage of block-pockets
« on: April 18, 2022, 07:47:20 PM »
Many of the moves use biasing (Multiple-first beads method). As soon as overlap is detected the move will be rejected.
Do you really need all the atoms outside the first cylindrical ring? If not, then you could
RestrictMovesToCylinder yes
and define a cylinder. Everything outside the cylinder will be considered an overlap.
Do you really need all the atoms outside the first cylindrical ring? If not, then you could
RestrictMovesToCylinder yes
and define a cylinder. Everything outside the cylinder will be considered an overlap.
3
General / Re: CIF file
« on: April 18, 2022, 07:43:09 PM »
CIF files are published by experimental groups to convert the structural data in a standardized way. The are uploaded to databases like CSD database.
For simulation codes, the atomic information is necessary _input_ for the simulations, because the positions of all the atoms need to be specified.
Programs like 'Openbabel' can convert between many formats.
For simulation codes, the atomic information is necessary _input_ for the simulations, because the positions of all the atoms need to be specified.
Programs like 'Openbabel' can convert between many formats.
4
Input files and parameters / Re: lennard jones parameters
« on: April 18, 2022, 07:39:16 PM »
Input units for the Lennard-Jones are epsilon/k_B in units of Kelvin, and sigma in units of Angstrom (see manual, also for units and unit conversions).
5
Input files and parameters / Re: UNIT CELL and affect on adsorption output ---- important
« on: April 18, 2022, 07:34:27 PM »
The smallest perpendicular width of the simulation box needs to be at least twice the cutoff (so as to satisfy the minimum image convention).
But you probably need larger cells to reduce finite size effects. A too large cell is computationally expensive.
But you probably need larger cells to reduce finite size effects. A too large cell is computationally expensive.
6
General / Re: Enthalpy of adsorption for binary mixture
« on: March 14, 2022, 07:13:39 PM »
The fluctuation formula has to be used in the grand-conical ensemble, i.e. the molecules must fluctuate.
Numerical problems might also occur at high loadings where the lack of accepted insertions/deletions can become a problem.
Numerical problems might also occur at high loadings where the lack of accepted insertions/deletions can become a problem.
7
Input files and parameters / Re: HYDRATED_ION_WATER potential
« on: March 14, 2022, 07:09:14 PM »
You do a: grep HYDRATED_ION_WATER *.c
in the source code and get the files back where this is used.
In the file 'potentials.c' you can check whether it is read okay, or printed out wrongly, or whether it is something with the mixing rule.
in the source code and get the files back where this is used.
In the file 'potentials.c' you can check whether it is read okay, or printed out wrongly, or whether it is something with the mixing rule.
8
Input files and parameters / Re: ExternalTemperature does not work like ExternalPressure
« on: March 14, 2022, 07:05:28 PM »
Yes, that is the standard way of simulation desorption. Note that is best to start from an equilibrated restart-file at the highest pressure (usually obtained from the adsorption branch).
9
Simulation algorithms and theory / Re: can I simulate the adsorption isothermal of NaCl solution in MOF
« on: March 14, 2022, 07:01:26 PM »
That can not really be done without implementing that. The system must stay charge-neutral, so you cannot individually swap na and cl in and out. You will need to do that in pairs.
10
Input files and parameters / Re: amorphous polymer can be used as adsorbent
« on: March 14, 2022, 06:58:12 PM »
The positions of the atoms of the framework are 'input'. There might be tools available (builders) that create these initial structures for you.
11
Input files and parameters / Re: Segmentation fault (core dumped )
« on: March 14, 2022, 06:56:34 PM »
See the tips in section 5.5 of the manual on debugging to find the cause of this.
Basically, compile with the flag '-g' and run the executable via a debugger.
Basically, compile with the flag '-g' and run the executable via a debugger.
12
General / Re: Error loading VTK output and lack of H on methane
« on: March 14, 2022, 06:53:50 PM »
1) Methane is often modeled as a united-atom (a single bead).
2) The VTK output that can be read by iRASPA is a VTK density-file (see 'ComputeDensityProfile3DVTKGrid').
2) The VTK output that can be read by iRASPA is a VTK density-file (see 'ComputeDensityProfile3DVTKGrid').
13
Pictures and movies / Re: creating movies in iRASPA
« on: December 02, 2021, 10:58:18 AM »
Select both files when loading them into iRASPA. The first tab in the project-pane shows the project. The second tab shows the two "movies" that you have loaded: the cif file with 1 frame, and the pdb with many frames. The third tab shows the frames of the movie.
Note: there is help in iRASPA (see menu item 'Help').
Note: there is help in iRASPA (see menu item 'Help').
14
Simulation algorithms and theory / Re: Acceptance rate of swap (insertion and deletion moves)
« on: December 02, 2021, 10:54:24 AM »
I would say > 0.1%. Lower then that and you will really have to simulate a long time. With biasing you can increases the acceptance at the cost of computing more.
15
General / Re: NPT Simulations
« on: December 02, 2021, 10:48:39 AM »
I have never used it, not sure if it works. You will have to check the VolumeChange move and modify to only change in the z-direction leaving the framework intact.