As often discussed in this context, I disagree with this very unusual treatment of the electromagnetic potentials and the claim that a split of the fields in this sense makes much sense from the point of view of physics. It also doesn't take into account the very important feature of gauge invariance, which is among the most important subjects to introduce students to modern physics already in the context of classical electromagnetism.
The standard derivation (using the most simple case of "microscopic classical electrodynamics", i.e., the vacuum equations taking into account all charges and currents explicitly) is as follows. One starts from the homogeneous Maxwell equations (using Heaviside Lorentz units)
$$\vec{\nabla} \times \vec{E}+\frac{1}{c} \partial_t \vec{B}=0, \quad \vec{\nabla} \cdot \vec{B}=0.$$
From Helmholtz's decomposition theorem from the 2nd equation it follows that there is a vector potential such that
$$\vec{B}=\vec{\nabla} \times \vec{A}.$$
It's immediately clear that for the given magnetic field, ##\vec{A}## is only determined up to a gradient field of an arbitrary scalar field, i.e., the vector potential
$$\vec{A}'=\vec{A}-\vec{\nabla} \chi$$
for any scalar field ##\chi## discribes the same physics (gauge invariance).
Plugging this into the first equation, leads to
$$\vec{\nabla} \times (\vec{E}+1/c \partial_t \vec{A})=0,$$
and from Helmholtz's theorem this implies the existence of a scalar potential for the quantitiy in parantheses:
$$\vec{E}+\frac{1}{c} \partial_t \vec{A} = -\vec{\nabla} \Phi$$
or
$$\vec{E}=-\vec{\nabla} \Phi - \frac{1}{c} \partial_t \vec{A}.$$
Under the gauge transformation you have to also adapt this scalar potential
$$\vec{A}'=\vec{A}-\vec{\nabla} \chi \; \Rightarrow\; \Phi'=\Phi+\frac{1}{c} \partial_t \chi.$$
Indeed then
$$-\vec{\nabla} \Phi'-\frac{1}{c} \partial_t \vec{A}'=-\vec{\nabla} \Phi-\frac{1}{c} \partial_t \vec{A}=0.$$
The potentials then obey equations, derived from the inhomogeneous Maxwell equations, which however are gauge invariant and thus determine the potentials only up to a gauge transformation, and thus to make the solutions unique one has to impose a gauge constraint, fixing (or at least partially fixing) the gauge. Usual choices are the Lorenz gauge,
$$\vec{\nabla} \cdot \vec{A}+\frac{1}{c} \partial_t \Phi=0$$
and the Coulomb gauge
$$\vec{\nabla} \cdot \vec{A}=0.$$
This demonstrates that any split of the electromagnetic field components based on pieces of the expressions in terms of the potentials is not gauge invariant and thus physically not easy to interpret.
It is utmost important to stress that the direct application of Helmholtz's decomposition theorem to the potentials is not very convenient but one should rather use the retarded propgator of the d'Alembert operator (Lorenz gauge) and the Green's function of the Laplace operator (for the scalar potential in the Coulomb gauge). The use of the Helmholtz's decomposition theorem in the way shown in this Insights article is overly complicated leading to nonlocal descriptions of the fields. The final solution in terms of the retarded propagator for the fields (the socalled Jefimenko equations) which is the unique solution of the Maxwell equations in terms of the physical sources of the electromagnetic fields, which are the charge-current distributions and not some non-local integrals shown in the article, though they may be formally correct from a purely mathematical point of view.
The introduction of a vector potential for the electric field or a scalar potential for the magnetic field can be sometimes for calculations in special cases (though I'm not aware of the use of a vector potential for the electric field in the literature).