We study hairy black holes in Einstein-Maxwell-complex scalar theory in four-dimensional asymptotically global anti-de Sitter (AdS) spacetime when the Robin boundary conditions are imposed on the scalar field. This setup is dual to the double trace deformation of strongly interacting field theory on R × S ² by charged scalar operators. We identify the instability of the Reissner-Nordström-AdS (RNAdS) black holes under the Robin boundary conditions and construct backreacted geometries branching at the onset of the instability. Also considering associated horizonless geometries called boson stars, we obtain phase diagrams with fairly rich structure in the grand canonical ensemble depending on the boundary condition parameter or the deformation parameter, where phase transition occurs between thermal AdS, RNAdS, charged boson stars, and hairy black holes.