Heaviside–Lorentz units

Heaviside–Lorentz units (or Lorentz–Heaviside units) constitute a system of units and quantities that extends the CGS with a particular set of equations that defines electromagnetic quantities, named for Oliver Heaviside and Hendrik Antoon Lorentz. They share with the CGS-Gaussian system that the electric constant ε₀ and magnetic constant µ₀ do not appear in the defining equations for electromagnetism, having been incorporated implicitly into the electromagnetic quantities. Heaviside–Lorentz units may be thought of as normalizing ε₀ = 1 and µ₀ = 1, while at the same time revising Maxwell's equations to use the speed of light c instead.

The Heaviside–Lorentz unit system, like the International System of Quantities upon which the SI system is based, but unlike the CGS-Gaussian system, is rationalized, with the result that there are no factors of 4π appearing explicitly in Maxwell's equations. That this system is rationalized partly explains its appeal in quantum field theory: the Lagrangian underlying the theory does not have any factors of 4π when this system is used. Consequently, electromagnetic quantities in the Heaviside–Lorentz system differ by factors of √4π in the definitions of the electric and magnetic fields and of electric charge. It is often used in relativistic calculations, and are used in particle physics. They are particularly convenient when performing calculations in spatial dimensions greater than three such as in string theory.

Motivation

In the mid-late 19th century, electromagnetic measurements were frequently made in either the so-named electrostatic (ESU) or electromagnetic (EMU) systems of units. These were based respectively on Coulomb's and Ampere's Law. Use of these systems, as with to the subsequently developed Gaussian CGS units, resulted in many factors of 4π appearing in formulas for electromagnetic results, including those without any circular or spherical symmetry.

For example, in the CGS-Gaussian system, the capacitance of sphere of radius r is r while that of a parallel plate capacitor is ⁠A/4πd⁠, where A is the area of the smaller plate and d is their separation.

Heaviside, who was an important, though somewhat isolated, early theorist of electromagnetism, suggested in 1882 that the irrational appearance of 4π in these sorts of relations could be removed by redefining the units for charges and fields. In his 1893 book Electromagnetic Theory, Heaviside wrote in the introduction:

It is not long since it was taken for granted that the common electrical units were correct. That curious and obtrusive constant 4π was considered by some to be a sort of blessed dispensation, without which all electrical theory would fall to pieces. I believe that this view is now nearly extinct, and that it is well recognised that the 4π was an unfortunate and mischievous mistake, the source of many evils.

In plain English, the common system of electrical units involves an irrationality of the same kind as would be brought into the metric system of weights and measures, were we to define the unit area to be the area, not of a square with unit side, but of a circle of unit diameter. The constant π would then obtrude itself into the area of a rectangle, and everywhere it should not be, and be a source of great confusion and inconvenience. So it is in the common electrical units, which are truly irrational.

Now, to make a mistake is easy and natural to man. But that is not enough. The next thing is to correct it: When a mistake has once been started, it is not necessary to go on repeating it for ever and ever with cumulative inconvenience. — Oliver Heaviside (1893)

Length–mass–time framework

As in the Gaussian system (G), the Heaviside–Lorentz system (HL) uses the length–mass–time dimensions. This means that all of the units of electric and magnetic quantities are expressible in terms of the units of the base quantities length, time and mass.

Coulomb's equation, used to define charge in these systems, is F = q^G
₁q^G
₂ / r² in the Gaussian system, and F = q^HL
₁q^HL
₂ / (4πr²) in the HL system. The unit of charge then connects to 1 dyn⋅cm² = 1 statC² = 4π HLC², where 'HLC' is the HL unit of charge. The HL quantity q^HL describing a charge is then √4π times larger than the corresponding Gaussian quantity. There are comparable relationships for the other electromagnetic quantities (see below).

The commonly used set of units is the called the SI, which defines two constants, the vacuum permittivity (ε₀) and the vacuum permeability (μ₀). These can be used to convert SI units to their corresponding Heaviside–Lorentz values, as detailed below. For example, SI charge is √ε₀L³M / T². When one puts ε₀ = 8.854 pF/m, L = 1 cm, M = 1 g, and T = 1 s, this evaluates to 9.409669×10⁻¹¹ C, the SI-equivalent of the Heaviside–Lorentz unit of charge.

Comparison of Heaviside–Lorentz with other systems of units

This section has a list of the basic formulas of electromagnetism, given in the SI, Heaviside–Lorentz, and Gaussian systems. Here E and D are the electric field and displacement field, respectively, B and H are the magnetic fields, P is the polarization density, M is the magnetization, ρ is charge density, J is current density, c is the speed of light in vacuum, ϕ is the electric potential, A is the magnetic vector potential, F is the Lorentz force acting on a body of charge q and velocity v, ε is the permittivity, χ_e is the electric susceptibility, μ is the magnetic permeability, and χ_m is the magnetic susceptibility.

Maxwell's equations

Main article: Maxwell's equations

Maxwell's equations in SI, Heaviside–Lorentz, and Gaussian quantities

Name	SI quantities	Heaviside–Lorentz quantities	Gaussian quantities
Gauss's law (macroscopic)	$\nabla \cdot \mathbf {D} ^{\textsf {SI}}=\rho _{\text{f}}^{\textsf {SI}}$	$\nabla \cdot \mathbf {D} ^{\textsf {HL}}=\rho _{\text{f}}^{\textsf {HL}}$	$\nabla \cdot \mathbf {D} ^{\textsf {G}}=4\pi \rho _{\text{f}}^{\textsf {G}}$
Gauss's law (microscopic)	$\nabla \cdot \mathbf {E} ^{\textsf {SI}}=\rho ^{\textsf {SI}}/\varepsilon _{0}$	$\nabla \cdot \mathbf {E} ^{\textsf {HL}}=\rho ^{\textsf {HL}}$	$\nabla \cdot \mathbf {E} ^{\textsf {G}}=4\pi \rho ^{\textsf {G}}$
Gauss's law for magnetism	$\nabla \cdot \mathbf {B} ^{\textsf {SI}}=0$	$\nabla \cdot \mathbf {B} ^{\textsf {HL}}=0$	$\nabla \cdot \mathbf {B} ^{\textsf {G}}=0$
Maxwell–Faraday equation (Faraday's law of induction)	$\nabla \times \mathbf {E} ^{\textsf {SI}}=-{\frac {\partial \mathbf {B} ^{\textsf {SI}}}{\partial t}}$	$\nabla \times \mathbf {E} ^{\textsf {HL}}=-{\frac {1}{c}}{\frac {\partial \mathbf {B} ^{\textsf {HL}}}{\partial t}}$	$\nabla \times \mathbf {E} ^{\textsf {G}}=-{\frac {1}{c}}{\frac {\partial \mathbf {B} ^{\textsf {G}}}{\partial t}}$
Ampère–Maxwell equation (macroscopic)	$\nabla \times \mathbf {H} ^{\textsf {SI}}=\mathbf {J} _{\text{f}}^{\textsf {SI}}+{\frac {\partial \mathbf {D} ^{\textsf {SI}}}{\partial t}}$	$\nabla \times \mathbf {H} ^{\textsf {HL}}={\frac {1}{c}}\mathbf {J} _{\text{f}}^{\textsf {HL}}+{\frac {1}{c}}{\frac {\partial \mathbf {D} ^{\textsf {HL}}}{\partial t}}$	$\nabla \times \mathbf {H} ^{\textsf {G}}={\frac {4\pi }{c}}\mathbf {J} _{\text{f}}^{\textsf {G}}+{\frac {1}{c}}{\frac {\partial \mathbf {D} ^{\textsf {G}}}{\partial t}}$
Ampère–Maxwell equation (microscopic)	$\nabla \times \mathbf {B} ^{\textsf {SI}}=\mu _{0}\left(\mathbf {J} ^{\textsf {SI}}+\varepsilon _{0}{\frac {\partial \mathbf {E} ^{\textsf {SI}}}{\partial t}}\right)$	$\nabla \times \mathbf {B} ^{\textsf {HL}}={\frac {1}{c}}\mathbf {J} ^{\textsf {HL}}+{\frac {1}{c}}{\frac {\partial \mathbf {E} ^{\textsf {HL}}}{\partial t}}$	$\nabla \times \mathbf {B} ^{\textsf {G}}={\frac {4\pi }{c}}\mathbf {J} ^{\textsf {G}}+{\frac {1}{c}}{\frac {\partial \mathbf {E} ^{\textsf {G}}}{\partial t}}$

The electric and magnetic fields can be written in terms of the potentials A and ϕ. The definition of the magnetic field in terms of A, B = ∇ × A, is the same in all systems of units, but the electric field is ${\textstyle \mathbf {E} =-\nabla \phi -{\frac {\partial \mathbf {A} }{\partial t}}}$ in the SI system, but ${\textstyle \mathbf {E} =-\nabla \phi -{\frac {1}{c}}{\frac {\partial \mathbf {A} }{\partial t}}}$ in the HL or Gaussian systems.

Other basic laws

Other electrostatic laws in SI, Heaviside–Lorentz, and Gaussian quantities

Name	SI quantities	Heaviside–Lorentz quantities	Gaussian quantities
Lorentz force	$\mathbf {F} =q^{\textsf {SI}}\left(\mathbf {E} ^{\textsf {SI}}+\mathbf {v} \times \mathbf {B} ^{\textsf {SI}}\right)$	$\mathbf {F} =q^{\textsf {HL}}\left(\mathbf {E} ^{\textsf {HL}}+{\frac {1}{c}}\mathbf {v} \times \mathbf {B} ^{\textsf {HL}}\right)$	$\mathbf {F} =q^{\textsf {G}}\left(\mathbf {E} ^{\textsf {G}}+{\frac {1}{c}}\mathbf {v} \times \mathbf {B} ^{\textsf {G}}\right)$
Coulomb's law	$\mathbf {F} ={\frac {1}{4\pi \varepsilon _{0}}}{\frac {q_{1}^{\textsf {SI}}q_{2}^{\textsf {SI}}}{r^{2}}}\mathbf {\hat {r}}$	$\mathbf {F} ={\frac {1}{4\pi }}{\frac {q_{1}^{\textsf {HL}}q_{2}^{\textsf {HL}}}{r^{2}}}\mathbf {\hat {r}}$	$\mathbf {F} ={\frac {q_{1}^{\textsf {G}}q_{2}^{\textsf {G}}}{r^{2}}}\mathbf {\hat {r}}$
Electric field of stationary point charge	$\mathbf {E} ^{\textsf {SI}}={\frac {1}{4\pi \varepsilon _{0}}}{\frac {q^{\textsf {SI}}}{r^{2}}}\mathbf {\hat {r}}$	$\mathbf {E} ^{\textsf {HL}}={\frac {1}{4\pi }}{\frac {q^{\textsf {HL}}}{r^{2}}}\mathbf {\hat {r}}$	$\mathbf {E} ^{\textsf {G}}={\frac {q^{\textsf {G}}}{r^{2}}}\mathbf {\hat {r}}$
Biot–Savart law	$\mathbf {B} ^{\textsf {SI}}={\frac {\mu _{0}}{4\pi }}\oint {\frac {I^{\textsf {SI}}d\mathbf {l} \times \mathbf {\hat {r}} }{r^{2}}}$	$\mathbf {B} ^{\textsf {HL}}={\frac {1}{4\pi c}}\oint {\frac {I^{\textsf {HL}}d\mathbf {l} \times \mathbf {\hat {r}} }{r^{2}}}$	$\mathbf {B} ^{\textsf {G}}={\frac {1}{c}}\oint {\frac {I^{\textsf {G}}d\mathbf {l} \times \mathbf {\hat {r}} }{r^{2}}}$

Dielectric and magnetic materials

Below are the expressions for the macroscopic fields $\mathbf {D}$, $\mathbf {P}$, $\mathbf {H}$ and $\mathbf {M}$ in a material medium. It is assumed here for simplicity that the medium is homogeneous, linear, isotropic, and nondispersive, so that the susceptibilities are constants.

SI quantities		Heaviside–Lorentz quantities		Gaussian quantities
Dielectric	Magnetic	Dielectric	Magnetic	Dielectric	Magnetic
${\begin{aligned}\mathbf {D} ^{\textsf {SI}}&=\varepsilon _{0}\mathbf {E} ^{\textsf {SI}}+\mathbf {P} ^{\textsf {SI}}\\&=\varepsilon \mathbf {E} ^{\textsf {SI}}\end{aligned}}$	${\begin{aligned}\mathbf {B} ^{\textsf {SI}}&=\mu _{0}(\mathbf {H} ^{\textsf {SI}}+\mathbf {M} ^{\textsf {SI}})\\&=\mu ^{\textsf {SI}}\mathbf {H} ^{\textsf {SI}}\end{aligned}}$	${\begin{aligned}\mathbf {D} ^{\textsf {HL}}&=\mathbf {E} ^{\textsf {HL}}+\mathbf {P} ^{\textsf {HL}}\\&=\varepsilon \mathbf {E} ^{\textsf {HL}}\end{aligned}}$	${\begin{aligned}\mathbf {B} ^{\textsf {HL}}&=\mathbf {H} ^{\textsf {HL}}+\mathbf {M} ^{\textsf {HL}}\\&=\mu ^{\textsf {HL}}\mathbf {H} ^{\textsf {HL}}\end{aligned}}$	${\begin{aligned}\mathbf {D} ^{\textsf {G}}&=\mathbf {E} ^{\textsf {G}}+4\pi \mathbf {P} ^{\textsf {G}}\\&=\varepsilon \mathbf {E} ^{\textsf {G}}\end{aligned}}$	${\begin{aligned}\mathbf {B} ^{\textsf {G}}&=\mathbf {H} ^{\textsf {G}}+4\pi \mathbf {M} ^{\textsf {G}}\\&=\mu ^{\textsf {G}}\mathbf {H} ^{\textsf {G}}\end{aligned}}$
$\mathbf {P} ^{\textsf {SI}}=\chi _{\text{e}}^{\textsf {SI}}\varepsilon _{0}\mathbf {E} ^{\textsf {SI}}$	$\mathbf {M} ^{\textsf {SI}}=\chi _{\text{m}}^{\textsf {SI}}\mathbf {H} ^{\textsf {SI}}$	$\mathbf {P} ^{\textsf {HL}}=\chi _{\text{e}}^{\textsf {HL}}\mathbf {E} ^{\textsf {HL}}$	$\mathbf {M} ^{\textsf {HL}}=\chi _{\text{m}}^{\textsf {HL}}\mathbf {H} ^{\textsf {HL}}$	$\mathbf {P} ^{\textsf {G}}=\chi _{\text{e}}^{\textsf {G}}\mathbf {E} ^{\textsf {G}}$	$\mathbf {M} ^{\textsf {G}}=\chi _{\text{m}}^{\textsf {G}}\mathbf {H} ^{\textsf {G}}$
$\varepsilon ^{\textsf {SI}}/\varepsilon _{0}=1+\chi _{\text{e}}^{\textsf {SI}}$	$\mu ^{\textsf {SI}}/\mu _{0}=1+\chi _{\text{m}}^{\textsf {SI}}$	$\varepsilon ^{\textsf {HL}}=1+\chi _{\text{e}}^{\textsf {HL}}$	$\mu ^{\textsf {HL}}=1+\chi _{\text{m}}^{\textsf {HL}}$	$\varepsilon ^{\textsf {G}}=1+4\pi \chi _{\text{e}}^{\textsf {G}}$	$\mu ^{\textsf {G}}=1+4\pi \chi _{\text{m}}^{\textsf {G}}$

Note that The quantities $\varepsilon ^{\textsf {SI}}/\varepsilon _{0}$, $\varepsilon ^{\textsf {HL}}$ and $\varepsilon ^{\textsf {G}}$ are dimensionless, and they have the same numeric value. By contrast, the electric susceptibility $\chi _{\text{e}}$ is dimensionless in all the systems, but has different numeric values for the same material: $$\chi _{\text{e}}^{\textsf {SI}}=\chi _{\text{e}}^{\textsf {HL}}=4\pi \chi _{\text{e}}^{\textsf {G}}$$ The same statements apply for the corresponding magnetic quantities.

Advantages and disadvantages of Heaviside–Lorentz units

Advantages

• The formulas above are clearly simpler in HL units compared to either SI or Gaussian units. As Heaviside proposed, removing the 4π from the Gauss law and putting it in the Force law considerably reduces the number of places the π appears compared to Gaussian CGS units.
• Removing the explicit 4π from the Gauss law makes it clear that the inverse-square force law arises by the E field spreading out over the surface of a sphere. This allows a straightforward extension to other dimensions. For example, the case of long, parallel wires extending straight in the z direction can be considered a two-dimensional system. Another example is in string theory, where more than three spatial dimensions often need to be considered.
• The equations are free of the constants ε₀ and μ₀ that are present in the SI system. (In addition ε₀ and μ₀ are overdetermined, because ε₀ μ₀ = 1 / c².)

The below points are true in both Heaviside–Lorentz and Gaussian systems, but not SI.

• The electric and magnetic fields E and B have the same dimensions in the Heaviside–Lorentz system, meaning it is easy to recall where factors of c go in the Maxwell equation. Every time derivative comes with a 1 / c, which makes it dimensionally the same as a space derivative. In contrast, in SI units [E] / [B] is [c].
• Giving the E and B fields the same dimension makes the assembly into the electromagnetic tensor more transparent. There are no factors of c that need to be inserted when assembling the tensor out of the three-dimensional fields. Similarly, ϕ and A have the same dimensions and are the four components of the 4-potential.
• The fields D, H, P, and M also have the same dimensions as E and B. For vacuum, any expression involving D can simply be recast as the same expression with E. In SI units, D and P have the same units, as do H and M, but they have different units from each other and from E and B.

Disadvantages

• Despite Heaviside's urgings, it proved difficult to persuade people to switch from the established units. He believed that if the units were changed, "[o]ld style instruments would very soon be in a minority, and then disappear ...". Persuading people to switch was already difficult in 1893, and in the meanwhile there have been more than a century's worth of additional textbooks printed and voltmeters built.
• Heaviside–Lorentz units, like the Gaussian CGS units by which they generally differ by a factor of about 3.5, are frequently of rather inconvenient sizes. The ampere (coulomb/second) is reasonable unit for measuring currents commonly encountered, but the ESU/s, as demonstrated above, is far too small. The Gaussian CGS unit of electric potential is named a statvolt. It is about 300 V, a value which is larger than most commonly encountered potentials. The henry, the SI unit for inductance is already on the large side compared to most inductors; the Gaussian unit is 12 orders of magnitude larger.
• A few of the Gaussian CGS units have names; none of the Heaviside–Lorentz units do.

Textbooks in theoretical physics use Heaviside–Lorentz units nearly exclusively, frequently in their natural form (see below), HL system's conceptual simplicity and compactness significantly clarify the discussions, and it is possible if necessary to convert the resulting answers to appropriate units after the fact by inserting appropriate factors of c and ε₀. Some textbooks on classical electricity and magnetism have been written using Gaussian CGS units, but recently some of them have been rewritten to use SI units. Outside of these contexts, including for example magazine articles on electric circuits, Heaviside–Lorentz and Gaussian CGS units are rarely encountered.

Translating expressions and formulas between systems

To convert any expression or formula between the SI, Heaviside–Lorentz or Gaussian systems, the corresponding quantities shown in the table below can be directly equated and hence substituted. This will reproduce any of the specific formulas given in the list above.

Equivalence expressions between SI, Heaviside–Lorentz, and Gaussian unit systems

Name	SI units	Heaviside–Lorentz units	Gaussian units
electric field, electric potential	${\sqrt {\varepsilon _{0}}}\left(\mathbf {E} ^{\textsf {SI}},\varphi ^{\textsf {SI}}\right)$	$\left(\mathbf {E} ^{\textsf {HL}},\varphi ^{\textsf {HL}}\right)$	${\frac {1}{\sqrt {4\pi }}}\left(\mathbf {E} ^{\textsf {G}},\varphi ^{\textsf {G}}\right)$
displacement field	${\frac {1}{\sqrt {\varepsilon _{0}}}}\mathbf {D} ^{\textsf {SI}}$	$\mathbf {D} ^{\textsf {HL}}$	${\frac {1}{\sqrt {4\pi }}}\mathbf {D} ^{\textsf {G}}$
charge, charge density, current, current density, polarization density, electric dipole moment	${\frac {1}{\sqrt {\varepsilon _{0}}}}\left(q^{\textsf {SI}},\rho ^{\textsf {SI}},I^{\textsf {SI}},\mathbf {J} ^{\textsf {SI}},\mathbf {P} ^{\textsf {SI}},\mathbf {p} ^{\textsf {SI}}\right)$	$\left(q^{\textsf {HL}},\rho ^{\textsf {HL}},I^{\textsf {HL}},\mathbf {J} ^{\textsf {HL}},\mathbf {P} ^{\textsf {HL}},\mathbf {p} ^{\textsf {HL}}\right)$	${\sqrt {4\pi }}\left(q^{\textsf {G}},\rho ^{\textsf {G}},I^{\textsf {G}},\mathbf {J} ^{\textsf {G}},\mathbf {P} ^{\textsf {G}},\mathbf {p} ^{\textsf {G}}\right)$
magnetic B field, magnetic flux, magnetic vector potential	${\frac {1}{\sqrt {\mu _{0}}}}\left(\mathbf {B} ^{\textsf {SI}},\Phi _{\text{m}}^{\textsf {SI}},\mathbf {A} ^{\textsf {SI}}\right)$	$\left(\mathbf {B} ^{\textsf {HL}},\Phi _{\text{m}}^{\textsf {HL}},\mathbf {A} ^{\textsf {HL}}\right)$	${\frac {1}{\sqrt {4\pi }}}\left(\mathbf {B} ^{\textsf {G}},\Phi _{\text{m}}^{\textsf {G}},\mathbf {A} ^{\textsf {G}}\right)$
magnetic H field	${\sqrt {\mu _{0}}}\ \mathbf {H} ^{\textsf {SI}}$	$\mathbf {H} ^{\textsf {HL}}$	${\frac {1}{\sqrt {4\pi }}}\mathbf {H} ^{\textsf {G}}$
magnetic moment, magnetization	${\sqrt {\mu _{0}}}\left(\mathbf {m} ^{\textsf {SI}},\mathbf {M} ^{\textsf {SI}}\right)$	$\left(\mathbf {m} ^{\textsf {HL}},\mathbf {M} ^{\textsf {HL}}\right)$	${\sqrt {4\pi }}\left(\mathbf {m} ^{\textsf {G}},\mathbf {M} ^{\textsf {G}}\right)$
relative permittivity, relative permeability	$\left({\frac {\varepsilon ^{\textsf {SI}}}{\varepsilon _{0}}},{\frac {\mu ^{\textsf {SI}}}{\mu _{0}}}\right)$	$\left(\varepsilon ^{\textsf {HL}},\mu ^{\textsf {HL}}\right)$	$\left(\varepsilon ^{\textsf {G}},\mu ^{\textsf {G}}\right)$
electric susceptibility, magnetic susceptibility	$\left(\chi _{\text{e}}^{\textsf {SI}},\chi _{\text{m}}^{\textsf {SI}}\right)$	$\left(\chi _{\text{e}}^{\textsf {HL}},\chi _{\text{m}}^{\textsf {HL}}\right)$	$4\pi \left(\chi _{\text{e}}^{\textsf {G}},\chi _{\text{m}}^{\textsf {G}}\right)$
conductivity, conductance, capacitance	${\frac {1}{\varepsilon _{0}}}\left(\sigma ^{\textsf {SI}},S^{\textsf {SI}},C^{\textsf {SI}}\right)$	$\left(\sigma ^{\textsf {HL}},S^{\textsf {HL}},C^{\textsf {HL}}\right)$	$4\pi \left(\sigma ^{\textsf {G}},S^{\textsf {G}},C^{\textsf {G}}\right)$
resistivity, resistance, inductance	$\varepsilon _{0}\left(\rho ^{\textsf {SI}},R^{\textsf {SI}},L^{\textsf {SI}}\right)$	$\left(\rho ^{\textsf {HL}},R^{\textsf {HL}},L^{\textsf {HL}}\right)$	${\frac {1}{4\pi }}\left(\rho ^{\textsf {G}},R^{\textsf {G}},L^{\textsf {G}}\right)$

As an example, starting with the equation $$\nabla \cdot \mathbf {E} ^{\textsf {SI}}=\rho ^{\textsf {SI}}/\varepsilon _{0},$$ and the equations from the table $${\begin{aligned}{\sqrt {\varepsilon _{0}}}\ \mathbf {E} ^{\textsf {SI}}&=\mathbf {E} ^{\textsf {HL}}\\{\frac {1}{\sqrt {\varepsilon _{0}}}}\rho ^{\textsf {SI}}&=\rho ^{\textsf {HL}}\,.\end{aligned}}$$

Moving the factor across in the latter identities and substituting, the result is $$\nabla \cdot \left({\frac {1}{\sqrt {\varepsilon _{0}}}}\mathbf {E} ^{\textsf {HL}}\right)=\left({\sqrt {\varepsilon _{0}}}\rho ^{\textsf {HL}}\right)/\varepsilon _{0},$$ which then simplifies to $$\nabla \cdot \mathbf {E} ^{\textsf {HL}}=\rho ^{\textsf {HL}}.$$