We propose a systematic and simulation-free strategy for coarse graining of homopolymer melts, where each chain of Nm monomers is uniformly divided into N segments, with the spatial position of each segment corresponding to the center-of-mass of its monomers. We use integral-equation theories suitable for the study of equilibrium properties of polymers, instead of many-chain molecular simulations, to obtain the structural and thermodynamic properties of both original and coarse-grained (CG) systems, and quantitatively examine how the effective pair potentials between CG segments and the thermodynamic properties of CG systems vary with N. Our systematic and simulation-free strategy is much faster than those using many-chain simulations, thus effectively solving the transferability problem in coarse graining, and provides the quantitative basis for choosing the appropriate N-values. It also avoids the problems caused by finite-size effects and statistical uncertainties in many-chain simulations. Taking the simple hard-core Gaussian thread model [K. S. Schweizer and J. G. Curro, Chem. Phys. 149, 105 (1990)] as the original system, we demonstrate our strategy applied to structure-based coarse graining, which is quite general and versatile, and compare in detail the various integral-equation theories and closures for coarse graining. Our numerical results show that the effective CG potentials for various N and closures can be collapsed approximately onto the same curve, and that structure-based coarse graining cannot give thermodynamic consistency between original and CG systems at any N< Nm.

To highlight the importance of quantitative and parameter-fitting-free comparisons among different models/methods, we revisited the comparisons made by Groot and Madden [J. Chem. Phys. 108, 8713 (1998)] and Chen et al. [J. Chem. Phys. 122, 104907 (2005)] between their dissipative particle dynamics (DPD) simulations of the DPD model and the self-consistent field (SCF) calculations of the "standard" model done by Matsen and Bates [Macromolecules 29, 1091 (1996)] for diblock copolymer (DBC) A-B melts. The small values of the invariant degree of polymerization used in the DPD simulations do not justify the use of the fluctuation theory of Fredrickson and Helfand [J. Chem. Phys. 87, 697 (1987)] by Groot and Madden, and their fitting between the DPD interaction parameters and the Flory-Huggins χ parameter in the "standard" model also has no rigorous basis. Even with their use of the fluctuation theory and the parameter-fitting, we do not find the "quantitative match" for the order-disorder transition of symmetric DBC claimed by Groot and Madden. For lamellar and cylindrical structures, we find that the system fluctuations/correlations decrease the bulk period and greatly suppress the large depletion of the total segmental density at the A-B interfaces as well as its oscillations in A- and B-domains predicted by our SCF calculations of the DPD model. At all values of the A-block volume fractions in the copolymer f (which are integer multiples of 0.1), our SCF calculations give the same sequence of phase transitions with varying χN as the "standard" model, where N denotes the number of segments on each DBC chain. All phase boundaries, however, are shifted to higher χN due to the finite interaction range in the DPD model, except at f = 0.1 (and 0.9), where χN at the transition between the disordered phase and the spheres arranged on a body-centered cubic lattice is lower due to N = 10 in the DPD model. Finally, in 11 of the total 20 cases (f-χN combinations) studied in the DPD simulations, a morphology different from the SCF prediction was obtained due to the differences between these two methods.

It is concluded that the mean-field coil–globule transition of a polymer chain of finite length N immersed in a small-molecule solvent exhibits the type-I behavior; that is, it is either a first-order phase transition, a critical point, or a crossover depending on the location of the critical point. It becomes a second-order phase transition with respect to the solvent equality characterized by the Flory–Huggins parameter χ (or equivalently the second virial coefficient υ or the temperature T) only in the limit of N → ∞. Even in this limit, it still has the type-I behavior with respect to υN 1/2 (or equivalently (1 – 2χ)N 1/2).

Using fast lattice Monte Carlo (FLMC) simulations (Wang, Q. Sof t Matter 2009, 5, 4564) and the corresponding polymer lattice field theories, including the lattice self-consistent field and Gaussian-fluctuation (LGF) theories, we studied a model system of incompressible homopolymer melts on a hexagonal lattice, where each lattice site is occupied by a total of ρ 0 ≥ 1 polymer segments. We generalized the cooperative motion algorithm (Pakula, T. Macromolecules 1987, 20, 679), as well as the related vacancy diffusion algorithm (Reiter, J.; Edling, T.; Pakula, T. J. Chem. Phys. 1990, 93, 837), originally proposed for the self-and mutual-avoiding walk (where ρ 0 = 1) to the case of ρ 0 > 1, where our generalized algorithm is highly efficient (i.e., nearly rejection-free). On the other hand, we extended the method of Wang (Wang, Z.-G. Macromolecules 1995, 28, 570) to calculate various single-chain properties in LGF theory. Direct comparisons between FLMC and LGF results, both of which are based on the same Hamiltonian (thus without any parameter-fitting between them), unambiguously and quantitatively reveal the effects of non-Gaussian fluctuations neglected by the latter. We found that FLMC results approach LGF predictions with increasing ρ 0 , and that the leading order of non-Gaussian fluctuation effects on the singlechain properties is inversely proportional to ρ 0 2 . Our work suggests that theories capturing the first-order non-Gaussian fluctuation effects may give quantitative agreement with FLMC simulations of incompressible homopolymer melts at ρ 0 ≥ 2 in two and three dimensions.

No abstract