What is a "charged system" in this specific context?

I'm doing a polarization calculation for the first time. To specify one of the input parameters (the "reference point") correctly, I have to figure out if my system is "charged." If it is not charged, I don't have to worry about the reference point.

On one hand, I understand that a system made up of water molecules (for example) is not charged. But if some of the water molecules dissociate or we add salt that ends up as Cl- and a corresponding +1 ion, I'm not sure. I.e., does "charged system" refer to the entire system or does it mean that we need to account for what's going on locally inside the system?

• Could you add some information about what type of system you are trying to model, the software, type of calculation etc.? In molecular dynamics for example, the charged system would mean the charge of the main molecule (e.g. protein) is non-zero, even if the total system (i.e. the simulation system with solvents) may have zero charge. Apr 6, 2021 at 17:02
• +1 but I agree with NTS. You are asking what "charged means in this context" but what precisely is the context? Who said that you have to figure out if the system is "charged" in order to specify one of the input parameters correctly? Which software, and where did it say you have to do that? Apr 6, 2021 at 17:08
• @ShoubhikRMaiti, thank you. I would like to do MD in CP2K on systems that contain water with and without metal and dissolved salt. The documentation only shows how to switch on calculation of dipole moments. There's a subsection that mentions "reference point." A forum post I found elsewhere says that all choices of reference point should lead to the same polarization (for a Berry phase calculation) for a non-charged system. I have no idea if my system is considered to be charged or not. It is charge-neutral overall but like I said, I don't know if it's the overall charge that matters.
– NTS
Apr 6, 2021 at 17:17
• @SRMaiti was the response by NTS helpful? Oct 30, 2021 at 20:24
• @NikeDattani Yeah it was helpful, but I have only experience with classical MD, not much on material science, so I left it to others to answer the question. Nov 1, 2021 at 10:42

As a simple example of the origin dependence for charged systems, consider the dipole of $$\ce{H2+}$$, which should have a bond length of roughly $$\pu{2bohr}$$ and the charge split evenly among the two hydrogens. If we put the origin at one of the atoms: $$|\mu_\text{atom_origin}|=(\pu{0.5 e} \times \pu{0 bohr})+(\pu{0.5 e} \times \pu{2 bohr})=\pu{1 e*bohr}$$
However, if we put the origin at the center of mass (COM) (which is also the center of charge): $$|\mu_\text{COM_origin}|=(\pu{0.5 e} \times \pu{-1 bohr})+(\pu{0.5 e} \times \pu{1 bohr})=\pu{0 e*bohr}$$