You can use OptiFDTD to model this crystal.
The only thing is that you have to define the layers by hand (without using the photonic crystal tool) or using our script language. You also cannot use the PWE band solver to find bandgaps on this type of structure since it relies on the photonic crystal tool.