Thermodynamics and partitioning of homopolymers into a slit-A grand canonical Monte Carlo simulation study

Grand canonical ensemble Monte Carlo simulation (GCMC) combined with the histogram reweighting technique was used to study the thermodynamic equilibrium of a homopolymer solution between a bulk and a slit pore. GCMC gives the partition coefficients that agree with those from canonical ensemble Monte Carlo simulations in a twin box, and it also gives results that are not accessible through the regular canonical ensemble simulation such as the osmotic pressure of the solution. In a bulk polymer solution, the calculated osmotic pressure agrees very well with the scaling theory predictions both for the athermal polymer solution and the theta solution. However, one cannot obtain the osmotic pressure of the confined solution in the same way since the osmotic pressure of the confined solution is anisotropic. The chemical potentials in GCMC simulations were found to differ by a translational term from the chemical potentials obtained from canonical ensemble Monte Carlo simulations with the chain insertion method. This confirms the equilibrium condition of a polymer solution partition between the bulk and a slit pore: the chemical potentials of the polymer chain including the translational term are equal at equilibrium. The histogram reweighting method enables us to obtain the partition coefficients in the whole range of concentrations based on a limited set of simulations. Those predicted bulk-pore partition coefficient data enable us to perform further theoretical analysis. Scaling predictions of the partition coefficient at different regimes were given and were confirmed by the simulation data.