The knowledge of the exact structure of the optical system point spread function (PSF) enables a high-quality image reconstruction in fluorescence microscopy. Accurate PSF models account for the vector nature of light and the phase and amplitude modifications. Most existing real-space-based PSF models fall into a sampling pitfall near the center position, yielding to the violation of energy conservation. In this work, we present a novel, to the best of our knowledge, Fourier-based techniques for computing vector PSF and compare them to the state-of-the-art. Our methods are shown to satisfy the physical condition of the imaging process. They are reproducible, computationally efficient, easy to implement, and easy to modify to represent various imaging modalities.