Hollow lattice tensor gauge theories with bosonic matter

Higher rank gauge theories are generalizations of electromagnetism where, in addition to overall charge conservation, there is also conservation of higher rank multipoles such as the total dipole moment. In this work we study a four-dimensional lattice tensor gauge theory coupled to bosonic matter which has second rank tensor electric and magnetic fields and charge conservation on individual planes. Starting from the Hamiltonian, we derive the lattice action for the gauge fields coupled to <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"> <a:mi>q</a:mi> <a:mo>=</a:mo> <a:mn>1</a:mn> </a:math> , 2 charged scalars. We use the action formulation to carry out Monte Carlo simulations to map the phase diagram as a function of the gauge ( <c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"> <c:mi>β</c:mi> </c:math> ) and matter ( <e:math xmlns:e="http://www.w3.org/1998/Math/MathML" display="inline"> <e:mi>κ</e:mi> </e:math> ) couplings. We compute the nature of correlators at strong and weak coupling in the pure gauge theory and compare the results to numerical simulations. Simulations show that the naive weak coupling regime (small <g:math xmlns:g="http://www.w3.org/1998/Math/MathML" display="inline"> <g:mi>κ</g:mi> </g:math> , large <i:math xmlns:i="http://www.w3.org/1998/Math/MathML" display="inline"> <i:mi>β</i:mi> </i:math> ) does not survive in the thermodynamic limit. Instead, the strong coupling confined phase spans the whole phase diagram. It is a proliferation of instantons that destroys the weak coupling phase and we show, via a duality transformation, that the expected strong confinement is present in the analog of Wilson line correlators. For finite matter coupling at <k:math xmlns:k="http://www.w3.org/1998/Math/MathML" display="inline"> <k:mi>q</k:mi> <k:mo>=</k:mo> <k:mn>1</k:mn> </k:math> we find a single thermodynamic phase albeit with a first-order phase transition terminating in a critical end point. For <m:math xmlns:m="http://www.w3.org/1998/Math/MathML" display="inline"> <m:mi>q</m:mi> <m:mo>=</m:mo> <m:mn>2</m:mn> </m:math> it is known that the X-cube model with <o:math xmlns:o="http://www.w3.org/1998/Math/MathML" display="inline"> <o:msub> <o:mi mathvariant="double-struck">Z</o:mi> <o:mn>2</o:mn> </o:msub> </o:math> fractonic topological order is recovered deep in the Higgs regime. The simulations indeed reveal a distinct Higgs phase in this case.