Let's look at where the electromagnetic interaction comes from in hydrogen.
At first quantization you have a multiparticle system so the wavefunction is defined as $\psi=\psi(x_1,y_1,z_1,x_2,y_2,z_2,t)$ and the point is to write the Hamiltonian.
And the Hamiltonian comes from the Lagrangian. For a single particle of charge $q$ in an external electromagnetic potential the Lagrangian is $L(\vec x, \vec v)=-mc^2\sqrt{1-\frac{v^2}{c^2}}+q\vec A\cdot\vec v -q\Phi.$
Which leads to the Hamiltonian $H(\vec Q,\vec P,t)=\sqrt{(mc^2)^2+(\vec P-q\vec A)^2c^2}+q\Phi$ (which is approximately $mc^2 +(\vec P-q\vec A)^2/2m +q\Phi$ in the nonrelativistic limit).
If your goal is to find energy eigenstates you might want to consider the quantum versions of equilibrium and an isolated system so the center of mass of the system moves inertially and the two charged bodies exert roughly equal forces and so the $\approx$ 2000 times as massive proton moves pretty closely to inertially so in its frame the potential is all scalar and that goes inversely to $\sqrt{(x_1-x_2)^2+(y_1-y_2)^2+(z_1-z_2)^2}.$
This naturally suggests a change of variables where the two particle configuration is changed to a center of mass location and a relative position of the electron. The point of all this is that this is the kind of argument usually used to get a Hamiltonian that is a sum of a free particle (all kinetic) of the center of mass and a $P^2/2\mu -|e|/r$ for the relative position $\vec r$ and reduced mass $\mu.$
So the point is that you can see how things interact electromagnetically if you write the Lagrangian for all the charges, then do a Legendre transformation to get the Hamiltonian, then use that Hamiltonian or a nonrelativistic approximation (since electromagnetism is naturally already relativistic).
That is how it is done. And that's the first stab at an answer. Showing you how it is normally done, so you know the traditional approach for a Schrödinger equation level electromagnetic interaction.
So you generally need that Lagrangian for the electromagnetic interaction between the charges. But you want the Lagrangian to depend on the configuration and the time rate of change of position. But in general there isn't an obvious potential of one charge on the other that only depends on the current positions and velocities of the particles, there should be a time delay, you find the potential on the past light cone. But if you imagine the non relativistic limit as a high $c$ limit you could choose to use a wave equation for the potential with charge and current at the present time as the source. In the realm of nonrelativistic approximations the most common techniques for non experts is to either make something up and hope it is good or do an order by order (in v/c) expression for the system and truncate after some order.
An example system you can look at would be $H^+$, a proton with two electrons at which point you can ask yourself how the two electrons interact with each other electromagnetically. You can write a Lagrangian, get a Hamiltonian, approximate them to some nonrelativistic approximation, and then set up the quantum system. If you have two electrons you will have use the superselection rule, and will probably need to include spin to do the super selection (anti symmetrize the wavefunction under exchange of identical fermions).
And that's unfortunate, since it means you don't start with a known classical Lagrangian for a system with spin degrees of freedom. A most basic version could be to treat the spin as a magnetic dipole that is proportional to the spin, and include the terms in the Lagrangian corresponding to an interaction between a magnetic dipole and a vector potential.
OK, that's an attempt at an honest recap of how electromagnetic interactions are often dealt with in regular Schrödinger equation level first quantization. Sometimes people have external fields, do perturbation, etcetera. However you asked about a quantum source not an external source on a quantum feeler.
If you want to see an electromagnetic field outright, you want that to be an object in its own right. For that, at the classical level the potential itself is part of the configuration in addition to the locations of the charges. But when you write the Lagrangian for that it is a Lagrangian density and when it gets quantized you get quantum field theory which you said isn't what you want.
So that's another answer. To get a true electromagnetic field you need quantum field theory. If all you want is the electromagnetic effect of the quantum charged object then you can do it the same way you normally do, which I've outlined above.