We investigate in underdoped cuprates possible coexistence of the superconducting (SC) order at zero momentum and pair density wave (PDW) at momentum ${bf Q}=(pi, pi)$ in the presence of a Neel order. By symmetry, the $d$-wave uniform singlet pairing $dS_0$ can coexist with the $d$-wave triplet PDW $dT_{bf Q}$, and the $p$-wave singlet PDW $pS_{bf Q}$ can coexist with the $p$-wave uniform triplet $pT_0$. At half filling, we find the novel $pS_{bf Q}+pT_0$ state is energetically more favorable than the $dS_0+dT_{bf Q}$ state. At finite doping, however, the $dS_0+dT_{bf Q}$ state is more favorable. In both types of states, the variational triplet parameters, $dT_{bf Q}$ and $pT_0$, are of secondary significance. Our results point to a fully symmetric $mathrm{Z_2}$ quantum spin liquid with spinon Fermi surface in proximity to the Neel order at zero doping, and to intertwined $d$-wave triplet PDW fluctuations and spin moment fluctuations along with the dominant $d$-wave singlet SC at finite doping. The results are obtained by variational quantum Monte Carlo simulations.
Download