You can try to add Zernike Standard Phase or similar Zernikes surfaces for this case.

Did you try phase surface?

Yeah, I tried with Zernike fringe phase and put numbers for zernike coefficients values z1, z2, z3 etc. Thus I can reduce the aberrations.  But, can you please say how the aberrations are changing ? As the values I am putting india the lens editor for phase surface are probably in lens unit (mm) but the Zernike standard coefficients are in wave units. So, is there any relation between them?

Thanks for the previous comments.

Please check this article:

https://support.zemax.com/hc/en-us/articles/1500005575422-How-to-model-a-black-box-optical-system-using-Zernike-coefficients

Okay, I will check. thanks