Semicrystalline polymers are ubiquitous, yet despite their fundamental and industrial importance, the theory of homogeneous nucleation from a melt remains a subject of debate. A key component of the controversy is that polymer crystallization is a non-equilibrium process, making it difficult to distinguish between effects that are purely kinetic and those that arise from the underlying thermodynamics. Due to computational cost constraints, simulations of polymer crystallization typically employ non-equilibrium molecular dynamics techniques with large degrees of undercooling that further exacerbate the coupling between thermodynamics and kinetics. In a departure from this approach, in this study, we isolate the near-equilibrium nucleation behavior of a simple model of a melt of short, semiflexible oligomers. We employ several Monte Carlo methods and compute a phase diagram in the temperature–density plane along with two-dimensional free energy landscapes (FELs) that characterize the nucleation behavior. The phase diagram shows the existence of ordered nematic and crystalline phases in addition to the disordered melt phase. The minimum free energy path in the FEL for the melt–crystal transition shows a cooperative transition, where nematic order and monomer positional order move in tandem as the system crystallizes. This near-equilibrium phase transition mechanism broadly agrees with recent evidence that polymer stiffness plays an important role in crystallization but differs in the specifics of the mechanism from several recent theories. We conclude that the computation of multidimensional FELs for models that are larger and more fine-grained will be important for evaluating and refining theories of homogeneous nucleation for polymer crystallization.

Nonsolvent-induced phase separation (NIPS) is a popular method for creating polymeric particles with internal microstructure, but many fundamental questions remain surrounding the kinetics of the complex coupled mass transfer and phase separation processes. In this work, we use simulations of a phase-field model to examine how (i) finite domain boundaries of a polymer droplet and (ii) solvent/nonsolvent miscibility affect the NIPS process. To isolate the effects of phase separation kinetics and solvent/nonsolvent mass transfer on the NIPS process, we study two different cases. First, we investigate droplet concentrations that originate inside the two-phase region, where phase separation kinetics alone governs the microstructure. Second, we investigate the effects of solvent/nonsolvent mass transfer by studying droplet concentrations that begin outside the two-phase region, where both phase separation kinetics and mass transfer play a role. In both cases, we find that qualitative NIPS behavior is a strong function of the relative location of the initial droplet composition with respect to the phase diagram. We also find that polymer/nonsolvent miscibility competes with solvent/nonsolvent miscibility in driving NIPS kinetic behavior. Finally, we examine polymer droplets undergoing solvent/nonsolvent exchange and find that the model predicts droplets that shrink with nearly Fickian diffusion kinetics. We conclude with a brief perspective on the state of simulations of NIPS processes and some recommendations for future work.