Hi,
I wanted to test the Comb potential by simulating a simple bulk structure of cuprous oxide. Unfortunately the structure isn't stable in the simulation. I attached my input file and the restart file for the bulk configuration, because I don't know where the problem might be. I'm pretty sure that the configuration file for cuprous oxide is correct (It's based on a cubic cell with a lattice parameter of 4.27 A), although it looks like I started with an off equilibrium configuration.
Be that as it may, the problem is, that I'm not familiar with the comb potential, especially the charge equilibration is totally new to me. So maybe the error is located in a wrong setting and therefore it would be nice if the comb experts can take a loot at my input file for crucial mistakes. I should add, that wether or not defining new start values of the charges, the result keeps the same. The structure isn't stable and the temperature rises up to ca. 4000 K.
Best regards,
Andreas Rühl
#------------- inputfile ---------------------------
---------- Initialize Simulation