Possible three-dimensional nematic odd-parity superconductivity in Sr$_2$RuO$_4$

The superconducting pairing in Sr$_2$RuO$_4$ is widely considered to be chiral $p$-wave with $\vec{d}_{\boldsymbol k} \sim (k_x + ik_y)\hat{z}$, which belongs to the $E_u$ representation of the crystalline $D_{4h}$ group. However, this superconducting order appears hard to reconcile with a number of key experiments. In this paper, based on symmetry analysis we discuss the possibility of odd-parity pairing with inherent three-dimensional (3D) character enforced by the inter-orbital interlayer coupling and the sizable spin-orbit coupling (SOC) in the material. We focus on a yet unexplored $E_u$ pairing, which contains finite $(k_z \hat x$, $k_z \hat y)$-component in the gap function. Under appropriate circumstances a novel time-reversal invariant nematic pairing can be realized. This nematic superconducting state could make contact with some puzzling observations on Sr$_2$RuO$_4$, such as the absence of spontaneous edge current and no evidences of split transitions under uniaxial strains.