Reduced Density Matrix Optimization for the Two-Dimensional Hubbard Model