基底的厚度是波长的5-10倍时，该如何仿真？

最近在模拟一篇PRL文章的结果，Dynamical Electric and Magnetic Metamaterial Response at Terahertz Frequencies
样品是3微米厚的铜样品做在670微米厚的砷化镓基底上，根据文章响应频率在0.5THz，这个基底对THz波来说太厚了，对S曲线影响很大，使S曲线上下震荡。我选了一个单元尺寸为50*50*(670+3)，单位都为微米，利用Master/slave边界，上下设PML（其中一个FSS），平面波激励，仿真结果差异很大，S曲线上下震荡（见附图）。
若只模拟铜样品，结果差别也很大。
不知碰到这种基底相对波长较厚的情况，应该怎样建立单元模型（50*50*?），怎样模拟?

也在方面的但没遇到过，希望高手能够解答以下

有没有高手给指点一下呢?
我可不可以在建单元模型时,只包括基底的一部分,例如只有20微米厚,这时单元50*50*(20+3).我试了试,震动趋势确实不明显了,但S参数的曲线还是和文献差别挺大的.
我一是搞不清引起S参数值上下振动的原因?
二是怀疑PML边界的设置不太正确.我一个设的是FSS,另一个就是简单的rad.不知当单元模型取的薄时,PML在基底内的那个面边界如何设?