A new mathematical model for monitoring the temporal evolution of the ice crystal size distribution during freezing in pharmaceutical solutions.

The freezing step plays a key role in the overall economy of the vacuum freeze-drying of pharmaceuticals, since the nucleation and crystal growth kinetics determine the number and size distribution of the crystals formed. In this work, a new mathematical model of the freezing step of a (bio)pharmaceutical solution is developed and validated. Both nucleation and crystal growth kinetics are modeled and included in a one-dimensional population balance (1D-PBM) that describes, given the product temperature measurement, the evolution of the pore size distribution during freezing. The developed model is coupled with the real-time measurements obtained from an infrared video camera. The ending time of the primary drying stage, and the maximum temperature inside the material, simulated through a simplified model of the process and the pore distribution forecast, resulted in good agreement with experimental values. The resulting Process Analytical Technology (PAT) has the potential to boost the development and optimization of a freeze-drying cycle and the implementation of a physically grounded Quality-by-Design approach in the manufacturing of pharmaceuticals. A more general mathematical model, including the aforementioned population balance, of a vial filled with a solution of sucrose was also developed and used to further validate the approach.