Best eigensolver for the first 500 eigenvalues on periodic mesh

Hi all,

I am looking to solve Poisson-Schrodinger on a 3D domain with periodic BC on x-perp and y-perp domains, which requires the first several hundred eigenpairs.

I tried to use in-house eigensolver PINVT. However, it became very slow for such a large amound of eigenpairs. I then turned to ngsPETSc. However, I then found that the periodic meshing is not yet supported in the MPI-parallel version.

I think I might have overlooked something. Is there any in-house eigensolver that fits this problem?

Best,
Chuning