I'm surprised nobody mentioned the "electrostatic approximation" already. From what I understand, according to the electrostatic approximation, electric fields are curl-free. This means that the following derivative of the electrostatic field is equal to the zero vector:
\nabla\times\textbf{E}=\textbf{0}
In this equation, \nabla\times denotes the curl and \textbf{E} is the electrostatic field. Since \nabla\times\textbf{E}=\textbf{0}, that means there exists a scalar field such that
\textbf{E}\:=\:-\nabla\varphi
where \nabla\varphi denotes the gradient of \varphi. Therefore the vector field E has a vanishing curl because the curl of a gradient is always zero. The scalar field \varphi is called the electrostatic potential because it describes the potential energy that an electric charge will have when placed at a given point in the electrostatic field. So, this is why we define this quantity, because of the electrostatic approximation.
Also, the reason surface area of any surface is partitioned into sufficiently small area elements because the geometry of a general surface changes with position along the surface, so we partition the surface area and use a limiting process so that the largest area element will approach zero, and the smaller the area elements become, the more numerous they are, and the more accurate the value of the total surface area is. In other words, a double integral is used.
I believe the 'atomicity of charge' means that charge is discrete (quantized) because electric charges exist as individual subatomic particles, each one being a multiple of the elementary charge.