Let $Y=3k-1,X=n+1$, then, as @GerryMyerson said in the comments, you have an elliptic curve $Y^2+Y=X^3-X$, and as such it can only have finitely many integer solutions.
You can use Sage to find them:
sage: E=EllipticCurve([0,0,1,-1,0])
sage: E
Elliptic Curve defined by y^2 + y = x^3 - x over Rational Field
sage: E.integral_points()
[(-1 : 0 : 1), (0 : 0 : 1), (1 : 0 : 1), (2 : 2 : 1), (6 : 14 : 1)]
Or you can look it up in the elliptic curve database at LMFDB.
Three of the solutions have $Y=k=0$, the other two give $k=1,P_k=1$ and $k=5,P_k=35$, so the solutions @MarioCarneiro found are all of them.