I'm trying to run a silica-water system using tip4p water model for water. The issue is that the force field that I use for silica also has a coulombic contribution. Since tip4p is a 4-site charge model, how to overlay these two pair styles together to get an accurate description of silica-water systems.