Modularity functions maximization with nonnegative relaxation facilitates community detection in networks

► Community detection equals to trace maximization of modularity Laplacian. ► An efficient method to optimize modularity function with nonnegative relaxation. ► Our solutions are very close to the ideal community indicator matrix. ► Fuzzy communities can be found by the near-orthogonal columns of the solution.