A methodology is proposed for designing the stabilization process of polyacrylonitrile (PAN) fibers. In its core, this methodology is based on a model that describes the characteristic fiber length change during the treatment, through the de-convolution of the three main contributors (i.e. entropic shrinkage, creep, and chemical shrinkage). The model has the additional advantage of offering further insight into the physical and chemical phenomena taking place during the treatment. Validation of PAN-model prediction performance for different processing parameters was achieved as demonstrated by FTIR and DSC. Τensile testing revealed the effect of processing parameters on fiber quality, while model prediction demonstrated that ladder polymer formation is accelerated at temperatures over 200oC. Additionally, according the DSC and FTIR measurements predictions from the application of the model during stabilization seem to be more precise at high-temperature stabilization stages. It was shown that mechanical properties could be enhanced preferably by including a treatment step below 200oC, before the initiation of cyclization reactions. Further confirmation was provided via Raman spectroscopy, which demonstrated that graphitic like planes are formed upon stabilization above 200oC, and thus multistage stabilization is required to optimize synthesis of carbon fibers. Optical Microscopy proved that isothermal stabilization treatment did not severy alter the cross section geometry of PAN fiber monofilaments.