There is great interest in the adaptive optics (AO) science community to overcome the limitations imposed by incomplete knowledge of the point spread function (PSF). To address this limitation a program has been initiated at the W. M. Keck Observatory (WMKO) to demonstrate PSF determination for observations obtained with Keck AO science instruments. This paper aims to give a broad view of the progress achieved in this area. The concept and the implementation are briefly described. The results from on-sky on-axis NGS AO measurements using the NIRC2 science instrument are presented. On-sky performance of the technique is illustrated by comparing the reconstructed PSFs to NIRC2 PSFs. Accuracy of the reconstructed PSFs in terms of Strehl ratio and FWHM are discussed. Science cases for the first phase of science verification have been identified. More technical details of the program are presented elsewhere in the conference.