ABSTRACT We use the polaris radiative transfer code to produce simulated circular polarization Zeeman emission maps of the cyanide (CN) J = 1–0 molecular line transition for two types of protostellar envelope magnetohydrodynamic simulations. Our first model is a low-mass disc envelope system (box length L = 200 au), and our second model is the envelope of a massive protostar (L = 104 au) with a protostellar wind and a CN-enhanced outflow shell. We compute the velocity-integrated Stokes I and V, as well as the implied V/I polarization percentage, for each detector pixel location in our simulated emission maps. Our results show that both types of protostellar environments are in principle accessible with current circular polarization instruments, with each containing swaths of envelope area that yield percentage polarizations that exceed the 1.8 per cent nominal sensitivity limit for circular polarization experiments with the Atacama Large Millimeter/submillimeter Array. In both systems, high-polarization (≳1.8 per cent) pixels tend to lie at an intermediate distance away from the central star and where the line-centre opacity of the CN emission is moderately optically thin (τLC ∼ 0.1–1). Furthermore, our computed V/I values scale roughly with the density-weighted mean line-of-sight magnetic field strength, indicating that Zeeman observations can effectively diagnose the strength of envelope-scale magnetic fields. We also find that pixels with large V/I are preferentially co-located where the absolute value of the velocity-integrated V is also large, suggesting that locations with favourable percentage polarization are also favourable in terms of raw signal.