Let's say we have a system of $n$ stationary charges interacting via Coulomb potential. Let's ignore possible external electromagnetic fields. Moreover the system is quantum, and its wave function is $\psi(R_1, R_2, ...,R_n)$ where $R_k$ denotes the coordinate of charge $k$-th. The question is the following: given a gauge transformation with a phase $\alpha(R_1,R_2, ...,R_n)$, is it correct to write the wave function transformation as \begin{equation} \psi'(R_1, R_2, ...,R_n) = e^{iQ\alpha(R_1,R_2, ...,R_n)} \psi(R_1, R_2, ...,R_n) \, , \end{equation} where $Q$ is the total charge, i.e $Q=\sum_{i} Q_i$. How do you derive the total charge conservation then?
What about the following relation: \begin{equation} \psi'(R_1, R_2, ...,R_n) = e^{iQ_1g_1(R_1)} e^{iQ_2g_2(R_2)} \cdots \psi(R_1, R_2, ...,R_n) \, , \end{equation} where $g_l$ are some functions. I am thinking of phases which change at each charge position independently from each other. Can it be still considered a gauge transformation?
This case probably holds only insofar $\psi$ may be written as a product, i.e., $\psi_1(R_1)\psi_2(R_2)\cdots \psi_n(R_n)$. Does this imply that the charges are non-interacting at all? Please let me know what do you think.