Electronic nematicity plays an important role in iron-based superconductors. These materials have a layered structure and the theoretical description of their magnetic and nematic transitions has been well established in the two-dimensional approximation, i.e., when the layers can be treated independently. However, the interaction between iron layers mediated by electron tunneling may cause nontrivial three-dimensional behavior. Starting from the simplest model for orbital nematic in a single layer, we investigate the influence of interlayer tunneling on the bulk nematic order and a possible preemptive state where this order is only formed near the surface. We found that the interlayer tunneling suppresses the bulk nematicity, which makes favorable the formation of a surface nematic order above the bulk transition temperature. The purely electronic tunneling Hamiltonian, however, favors a nematic order parameter that alternates from layer to layer. The uniform bulk state typically observed experimentally may be stabilized by the coupling with the elastic lattice deformation. Depending on the strength of this coupling, we found three regimes: (i) surface nematic and alternating bulk order, (ii) surface nematic and uniform bulk order, and (iii) uniform bulk order without the intermediate surface phase. The intermediate surface-nematic state may resolve the current controversy about the existence of a weak nematic transition in the compound ${\mathrm{BaFe}}_2{\mathrm{As}}_{2-x}{\mathrm{P}}_x$.

• Received 27 April 2016
• Revised 26 August 2016

DOI:https://doi.org/10.1103/PhysRevB.94.094509