Depending on how detailed you need the model to be we already have laser models in OptiSystem to model I-V curves and spatial modes (VCSEL). You can specify whether you want the laser to generate x – y polarized light.
Do you have information relating to cavity size, differential gain, cavity loss, etc.?