I've been working for some time attempting to construct the CLAY force field (Cygan, 2004) using the potentials/Coulomb solvers available in ATK. My ultimate goal is to measure diffusion of water through a mesoporous silica nanopore (MSN). Per another question on the forum about creating gaps in the periodic images, the z lattice parameter is 30 Angstrom longer than the system (made up of water and the pore).
I'm running into several issues. First, I can't actually get a simulation to run without a few atoms breaking bonds and shooting off into the space beyond the defined lattice. Second, occasionally when it IS running, it stops and gives me the error or "Particle leaving simulation area."
I've checked the interatomic distances and charge balance several times.
I'm not going to describe the FF, etc, in too much detail here. It's all in the input file. I've included the input file and the tag files that go along with it. The FF is made up of the particle definitions, LJ spline potentials, and the SPME Coulomb solver. CLAY would typically involve several more particle types, but I only defined what I needed. The six particle types are the silicon in the bulk, the oxygen in the bulk(Ob), the hydroxyl oxygen attached to the pore surface(Ox), the hydroxyl hydrogen(Hx), water oxygen(Ow), and water hydrogen(Hw)
Can anyone reproduce what I'm seeing and have any advice on how to fix it?
