  With the onset of the era of gravitational-wave (GW) astronomy, the search
for continuous gravitational waves (CGWs), which remain undetected to date, has
intensified in more ways than one. Rapidly rotating neutron stars with
non-axisymmetrical deformations are the main targets for CGW searches. The
extent of this quadrupolar deformation is measured by the maximum ellipticity
that can be sustained by the crust of a neutron star and it places an upper
limit on the CGW amplitudes emitted by such systems. In this paper, following
previous works on this subject, we calculate the maximum ellipticity of a
neutron star generated by the Lorentz force exerted on it by the internal
magnetic fields. We show that the ellipticity of stars deformed by such a
Lorentz force is of the same order of magnitude as previous theoretical and
astrophysical constraints. We also consider if this ellipticity can be further
enhanced by crustal surface currents. We discover that this is indeed true;
surface currents at crustal boundaries are instrumental towards enhancing the
ellipticity of magnetized neutron stars.
