Magnetization density distribution of Sr$_2$IrO$_4$: Deviation from a local $j_\text{eff}=1/2$ picture

$5d$ iridium oxides are of huge interest due to the potential for new quantum states driven by strong spin-orbit coupling. The strontium iridate Sr$_2$IrO$_4$ is particularly in the spotlight because of the novel $j_\text{eff}=1/2$ state consisting of a quantum superposition of the three $t_{2g}$ orbitals with nearly equal population, which stabilizes an unconventional Mott insulating state. Here, we report an anisotropic and aspherical magnetization density distribution measured by polarized neutron diffraction in a magnetic field up to 5~T at 4~K, which strongly deviates from a local \jeffHalf picture. Once reconstructed by the maximum entropy method and multipole expansion model refinement, the magnetization density shows cross-shaped positive four lobes along the crystallographic tetragonal axes with a large spatial extent, showing that the $xy$ orbital contribution is dominant. Theoretical considerations based on a momentum-dependent composition of the $j_\text{eff}=1/2$ orbital and an estimation of the different contributions to the magnetization density casts the applicability of an effective one-orbital $j_\text{eff}=1/2$ Hubbard model into doubt. The analogy to the superconducting copper oxide systems might thus be weaker than commonly thought.