I'm relatively new to MD simulation, and am studying a highly viscous system of very concentrated aqueous biphasic systems of two different electrolytes. The higher-concentration systems are notoriously slow to equilibrate, it can take over 500 ns to equilibrate (from literature). When I run npt equilibrations with the barostat set to 1 bar, the average pressure through the simulation increases with the concentration - for example, for 10 molal in both phases, av. pressure is 2 bar, whereas for 20 molal it's 7 bar. However, the density of the systems doesn't change through the sim. Would this suggest that the systems are not equilibrated, or is this just to do with how the pressure is calculated?
I'm aware that the pressure can fluctuate a lot with RMSD in the hundreds. I just thought it weird that the average pressure would increase with increasing concentration?
I'm using a modified AMBER force field for one phase with scaled charges of 0.85 for ions, the Madrid force field for the other, and TIP4P/2005 for the water. Some boxes are 5 x 5 x 10 nm and others 10 x 10 x 30 nm