Methodical evaluation of Boyle temperatures using Mayer sampling Monte Carlo with application to polymers in implicit solvent

The Boyle temperature, TB, for an n-segment polymer in solution is the temperature where the second osmotic virial coefficient, A2, is zero. This characteristic is of interest for its connection to the polymer condensation critical temperature, particularly for n → ∞. TB can be measured experimentally or computed for a given model macromolecule. For the latter, we present and examine two approaches, both based on the Mayer-sampling Monte Carlo (MSMC) method, to calculate Boyle temperatures as a function of model parameters. In one approach, we use MSMC calculations to search for TB, as guided by the evaluation of temperature derivatives of A2. The second approach involves numerical integration of an ordinary differential equation describing how TB varies with a model parameter, starting from a known TB. Unlike general MSMC calculations, these adaptations are appealing because they neither invoke a reference for the calculation nor use special averages needed to avoid bias when computing A2 directly. We demonstrate these methods by computing TB lines for off-lattice linear Lennard-Jones polymers as a function of chain stiffness, considering chains of length n ranging from 2 to 512 monomers. We additionally perform calculations of single-molecule radius of gyration Rg and determine the temperatures Tθ, where linear scaling of Rg2 with n is observed, as if the polymers were long random-walk chains. We find that Tθ and TB seem to differ by 6% in the n → ∞ limit, which is beyond the statistical uncertainties of our computational methodology. However, we cannot rule out systematic error relating to our extrapolation procedure as being the source of this discrepancy.