Density vs temperature profiles at different cooling rates and annealing simulations are presented in Fig. . At the highest cooling rate (5000 K/ns), significant deviations in densities are observed at higher temperatures during the heating simulation compared to the cooling simulations. This discrepancy may arise from the short simulation time, during which the polymer chains were not adequately equilibrated. As the cooling rates decrease, the observed deviation diminishes, suggesting that slower cooling rates or longer simulation times produce more accurate density profiles. Figure illustrates that density profiles and T g values from continuous cooling simulations at the lowest cooling rates (50 and 5 K/ns) closely match those from step-wise cooling simulations. The calculated glass transition temperatures differ by approximately 30 K between the various cooling rates.