You would need OptiFDTD for this application. Not necessarily the OptiFDTD simulation engine but the additional software that comes with it, OptiMode. In OptiMode you can construct photonic crystal fibers and simulate the modes that exist. You can also run a parameter sweep to find which wavelength corresponds to cutoff.
Regards