Hey Katya,

If you're using hybrid mode NSC, make sure the stop is before the NSC group and use an 'object-specific' aperture definition like entrance pupil diameter or float-by-stop. You'll probably also need to use a plane phase reference and have pretty high sampling. Check out the sample file Documents\Zemax\Samples\Physical Optics\Lenslet Array.ZMX for an example of a similar system.

BUT

There's a little-known capability in Non-sequential that might do what you need better. Check out the sample file Samples\Non-sequential\Coherence Interference and Diffraction\Interference Example 4- Diffraction Limited Imaging.ZMX. Note the use of the Huygens Mode flag on detector object 6. This adds the rays landing on the detector as a coherent sum of plane waves to give the production of a full PSF rather than a geometric point image:

This may be the easier way of getting what you want. I recommend raeding the documentation on the Huygens Mode of the detector

Quick follow up, check out https://my.zemax.com/en-US/Knowledge-Base/kb-article/?ka=KA-01402 for a full discussion of the POP and lenslet approach. The same thing should work for POP and DMD.