Maybe there's a simple workaround. If you calculate two images.
One with ior 1.6584 of the interior of the crystal and another
with the other ior of 1.4864. And than combine this images in photoshop
or similar programs.
It's not fully simulation, but...
SciBorg