Consistent and transferable coarse-grained model for semidilute polymer solutions in good solvent

Consistent and transferable coarse-grained model for semidilute polymer solutions in good solvent

Giuseppe D’Adamo Physics Department, University of L’Aquila, Via Vetoio, 67100 L’Aquila, Italy    Andrea Pelissetto Physics Department, “Sapienza” University of Rome, and INFN, Sezione di Roma I, P.le A.Moro 2, 00185 Rome, Italy    Carlo Pierleoni Physics Department, University of L’Aquila and CNISM, UdR L’Aquila, Via Vetoio, 67100 L’Aquila, Italy
July 5, 2019

We present a coarse-grained model for linear polymers with a tunable number of effective atoms (blobs) per chain interacting by intra- and inter-molecular potentials obtained at zero density. We show how this model is able to accurately reproduce the universal properties of the underlying solution of athermal linear chains at various levels of coarse-graining and in a range of chain densities which can be widened by increasing the spatial resolution of the multiblob representation, i.e., the number of blobs per chain. The present model is unique in its ability to quantitatively predict thermodynamic and large scale structural properties of polymer solutions deep in the semidilute regime with a very limited computational effort, overcoming most of the problems related to the simulations of semidilute polymer solutions in good solvent conditions.


I Introduction

The last decade has witnessed a considerable effort in developing coarse-grained (CG) models to bridge the length-scale gap between the microscopic scale and the meso(macro)-scopic scale typical of soft-matter and biological systems. Two general strategies have been attempted: the structure-based route in which CG models are tuned to reproduce some structural properties at the local scale and the thermodynamic-based route in which the models are required to reproduce solvation free energies. Also mixed strategies have been developed and applied to a large variety of physical systems. An overview of methods and recent applications can be found in Refs. Voth2009 (); PCCP2009 (); SM2009 (); FarDisc2010 ().

The structure-based route groups a number of atoms into effective “particles” and assumes state-dependent pair interactions between them. These potentials are derived from the local structure of the atomic-level system, using the Iterative Boltzmann Inversion (IBI) Carbone:2008p2254 (); Fritz:2009p1721 (); Peter:2009p1734 (); Tschop98 (), the Inverse Monte Carlo (IMC) Lyubartsev:2009p2497 () or by liquid-state theory techniques Hansen:2005p292 (); Clark:2010p2483 (). The state dependence of the effective pair potentials hides the underlying many-body character and entropic content of the effective interactions and poses the question of their transferability. Indeed, if the effective potentials are not transferable, setting up the CG model requires deriving the effective interactions for all thermodynamic states of interest, seriously limiting the benefits of the coarse-graining strategy. Generally speaking, the problem of transferability remains unsolved.

In the specific case of linear polymers in solutions, structure-based CG models, where a single coil is mapped onto a single soft blob, have been introduced quite some time ago Flory:1950p2484 (); Grosberg:1982p2265 (); Dautenhahn:1994p2250 (); Murat:1998p1980 (); Bolhuis:2001p268 (); D'Adamo2012 (); Pelissetto:2005p296 (). For two isolated coils in athermal solvent the effective interaction between the centers of mass of the coils is roughly Gaussian, of the order of at overlap, and with a range of the order of the coil radius of gyration. This single-blob model is only accurate in the dilute regime, in which chain overlaps can be neglected. Extension of the single-blob CG model to the semidilute regime can be obtained by allowing the use of density dependent pair potentialsBolhuis:2001p268 () which can be determined from the radial distribution function between the centers of mass of two coils computed from a full monomer simulation at finite density. This strategy however suffers for the transferability problem since setting up the coarse-grained model requires simulations of the underlying full monomer system at all densities of interest. Another limitation of the single-blob model with density dependent potential is in representing non-homogeneous systems since the interaction should depend on the local density which is not known beforehand and some kind of self-consistent procedure should be developed. Furthermore, representing polymers as soft spherically symmetric particles is not always appropriate. For instance, in studying polymers adsorbed on surfaces, like polymer brushes or polymer-coated colloids, it is clear that the anchorage to the surface breaks the rotational symmetry of the chains, an effect that must be taken into account in any accurate coarse-grained modelColuzza2008 (). A further example is in modelling solutions of A-B block copolymers which cannot be represented as soft particles interacting by a spherically symmetric pair potential Pierleoni:2006p159 (); HansenJP:2006p2248 (); Sambriski:2007p781 (); Capone2009 (); Gross:2010p2528 (); Capone2010 (); Capone2011 ().

In a semidilute solution of linear chains of monomers in a volume , chain density, , is larger than the overlap density ( is the isolated coil radius of gyration, the monomer size and the scaling exponent) while the monomer density is still very small, a condition that can always be satisfied if chains are long enough Doi (). An accurate description of the thermodynamic and large-scale structural behavior of polymers in these conditions can in principle be obtained by using CG chains, in which a number of the original monomers are grouped in one effective monomer (blob) to map the original chain of monomers onto a chain of effective blobs (multiblob model). If the level of coarse graining, i.e. the number of blobs per chain, is such that the blob density is below the overlap blob density ( is the zero-density radius of gyration of the blob), then it is expected that zero-density potentials between blobs or different chains can be safely used. These ideas were discussed in Refs. Pierleoni:2007p193 (); Pelissetto:2009p287 (); Vettorel:2010p1733 () and recently reviewed in Ref. DAdamoG:2011p2477 (). Although in principle very appealing, the problem with this approach is how to obtain the intramolecular potential. Indeed, this inherently many-body potential is in principle of increasing complexity when increasing the number of blobs per chain, and since blobs are tethered together in some fixed topology, a zero-density expansion cannot be invoked to decompose it into a sum of two-body, three-body, etc. terms. In Ref. DAdamoG:2011p2477 () a coarse-grained model for good-solvent conditions with four effective monomers per chain (tetramer) was developed and throughly studied. Such a model, set up at zero density, was found to be accurate up to reduced polymer density supporting the multiblob ideas.

In this paper we present the extension of the tetramer model to chains with an arbitrary number of blobs, an extension needed to explore the semidilute regime at high polymer concentration. We show that a careful parametrization of the -body interactions up to is enough to obtain a very accurate CG potential field fully transferable both with the number of blobs per chain and with the density of chains, provided that the length of the effective chains is chosen to ensure that the reduced blob density remains always low. The paper is organized as follows. In section II we will provide a theoretical framework on which any accurate coarse-grained model should be based. In the following section III we present our multiblob model and report results for systems of chains of varying length both at zero density (section III.a) and at finite density in the semidilute regime (section III.b). In section IV we collect our conclusion and perspectives. Finally an appendix reports the tetramer potentials as obtained by the IBI and their explicit parametrization.

Ii Theoretical framework

A crucial requirement of any accurate coarse-grained representation (CGR) representation of the underlying full-monomer (FM) chain model is to preserve the value of the radius of gyration for all values of the number of blobs . Let us consider a chain of monomers mapped onto a chain of blobs, each representing the center of mass of a subchain of monomers. The fundamental relation among the chain radius of gyration , the radius of gyration of the chain of blobs, and the mean blob radius of gyration , holds for any single configuration and therefore on average. Here is the average of the square radius of gyration of all blobs in the chain




where is the position of the -th monomer in the chain. If is the center of mass of the -th blob in the chain,


is defined by


In the scaling limit, it has been found that DAdamoG:2011p2477 ()


where Clisby:2010p2249 (). Here and in the following we will use a hat to indicate zero-density averages. For , and .

Any consistent multiblob (MB) model must obey these relations when varying . A proper definition of , independent on , including the prefactor, is crucial since chain structural properties are universal only if distances are expressed in terms of rescaled distances . Moreover, at finite density results from different values of should be compared at the same value of the polymer volume fraction . Changing the definition of with , changes the definition of both and and the comparison among models with different becomes meaningless.

In Ref. DAdamoG:2011p2477 () a MB model with four blobs per chain (tetramer) was developed and throughly studied. Central potentials between first, second and third neighbors along the chain, as well as a bending and torsional angle potentials were determined by IBI to reproduce the pair distances and angular distributions of an isolated FM chain in the scaling limit, mapped onto four blobs DAdamoG:2011p2477 (). Moreover, an intermolecular Gaussian potential between any pair of blobs of different tetramers was assumed and optimized to reproduce the radial distribution function between the centers of mass of two isolated chains in the scaling limit Pelissetto:2005p296 (). The model perfectly reproduces the structure of the chains and provides accurate results for the thermodynamics of the solution up to . The non-uniform (universal) angular distributions of the tetramer are induced by the specific mapping of subchains onto their centers of mass, a procedure which introduces an explicit angular correlation along the chain even for ideal chains Laso:1991p1417 (). Failure in reproducing the correct angular distributions is at the basis of the inconsistent behavior observed for simpler multiblob models Pierleoni:2007p193 (); Pelissetto:2009p287 (); DAdamoG:2011p2477 (). Note that while central potentials are two-body interactions, bending and torsional angle potentials represent genuine three- and four-body interactions, respectively. Interestingly DAdamoG:2011p2477 (), a hierarchical order in the intensity of the various potentials is observed, the strongest and most relevant being the central first-neighbor interaction (bonding), followed by the second- and third-neighbor central interactions, and finally by the bending and the torsional angle potentials. Therefore, correlations along the chain decrease with the chemical distance and many-body effects are smaller than two-body ones, although of different nature. This observation suggests to use tetramers as building blocks of a resolution-invariant multiblob model.

Iii Tetramer-based multiblob models and results

iii.1 Zero-density transferability

A first tetramer-based multiblob model (4MB-1) is obtained by using the tetramer potentials of Ref. DAdamoG:2011p2477 () reported in the Appendix. In model 4MB-1 one considers six different intramolecular potentials, each of them depending on a single scalar variable. We use four bonding pair potentials: blobs and interact with a central potential with . Symmetry imposes and so that there are only four different central potentials. Then we use a three-body potential acting on the two equivalent bending angles and a four-body potential acting on the torsional angle, where and are defined as


where is the bond vector between blobs ad . The intermolecular potential is instead represented by a unique blob-blob pair interaction acting between any pair of blobs of different tetramers.

When using this potential field to build chains of blobs, the blob radius of gyration is taken as the length scale of the central potentials. In practice this means that when transferring the potentials from the model with 4 blobs to the model with blobs, the following identification is adopted


where, according to Eq. (6), is defined by


The angular potentials are instead invariant under this scale transformation.

This model, 4MB-1, is not found to be fully satisfactory, as illustrated in table 1, where we report results for the zero-density chain radius of gyration and the second dimensionless virial coefficient for increasing for polymer chains in the scaling limit (FM) footnote (), the multiblob model (4MB-1), and the modified multiblob model (4MB-2) still to be discussed. In particular, in model 4MB-1, becomes greater than 1 for increasing , a fact clearly incompatible with the multiblob procedure [see Eq. (5)]. Also a systematic increase with in the second virial coefficient of 4MB-1 is observed which witnesses the inability of such model to be fully consistent at all resolutions.

FM-CGR 4MB-1 4MB-2 (4MB-1) (4MB-2)
4 0.8921(1) 0.8933(3) 0.8903(2) 5.597(1) 5.596(3)
10 0.9642(1) 0.9812(3) 0.9638(7) 5.676(2) 5.619(2)
20 0.9842(1) 1.0060(6) 0.983(3) 5.714(2) 5.614(3)
30 0.9902(1) 1.0131(8) 0.990(3) 5.734(3) 5.618(6)
Table 1: and universal ratio ( is the second virial coefficient) as a function of for polymers obtained by full-monomer simulations (FM-CGR), and for the original and modified 4MB models. The universal value for polymers of the second virial coefficient combination is Caracciolo:2006p587 ().
Figure 1: Probability distribution of in the original tetramer model (4MB-1) and in the modified tetramer model (4MB-2). The distribution for a CGR of a full-monomer chain with four blobs is also reported for comparison (FM-CGR). The inset displays the relative difference of the two models with the FM-CGR prediction.

The origin of such inconsistency can be ascribed to a residual inaccuracy of the original tetramer model to reproduce four-body correlations. In Fig. 1 we show, for , the probability distribution of . Data of the 4MB-1 model are compared with FM-CGR results. Despite the explicit presence of the angular potential which reproduces very accurately the torsional angle distribution DAdamoG:2011p2477 (), is not well reproduced in the 4MB-1 model—more elongated configurations are apparently enhanced. This inaccuracy, irrelevant for the properties of the multiblob model at the tetramer level (), accumulates when transferring the tetramer potentials to longer chains, producing the observed increase of and with . A second, modified tetramer model which reproduces the FM-CGR behavior can be obtained by adding a new potential on . As for the other potentials, the IBI procedure is applied to extract the new optimal potential. At the same time the dihedral angle potential is re-optimized to keep the level of accuracy of the original tetramer model. Explicit expressions for the potentials are provided in the Appendix. The good accuracy of the new model is illustrated in figure 1. The modified tetramer model can now be safely used as the building block of a fully consistent model (4MB-2) which puts in action the ideas of the multiblob approach. As illustrated in table 1, at zero density the agreement with the FM-CGR predictions for the rescaled radius of gyration is excellent for all values investigated. Also the virial coefficient is independent on and in good agreement with the FM prediction Caracciolo:2006p587 (). This is the first requirement of transferability with at zero density. Note that the third virial coefficient appears to be much less sensible to model inaccuracies: we obtain both for 4MB-1 and 4MB-2 for all values investigated to be compared with the FM predictionCaracciolo:2006p587 (): .

To check that our parametrization exhausts all four-body terms we compare in Fig. 2 the cross correlation between the cosine of the two bending angles of the tetramer as obtained by the FM-CGR and by the 4MB-2 model. This quantity is defined as


In the absence of correlation we have . We observe the presence of a region of depleted correlation near , and two regions of enhanced correlation around and , respectively. The 4MB-2 model reproduces fairly well the results of the FM-CGR (see the right panel), except in a very narrow region around where statistical noise is also larger, indicating the accuracy of our parametrization.

Figure 2: Cross-correlation function . Comparison between the 4MB-2 model (left panel) FM-CGR chains (central panel) for at . In the right panel, we show the relative deviations of the 4MB-2 results from the FM-CGR ones, .

To further investigate if of the new model 4MB-2 is able to reproduce FM-CGR correlations, we computed genuine five-body intramolecular correlations for varying number of blobs . Note that two-, three- and four-body correlations at short chemical distance are explicitly built into the model, and five-body correlations are the first missing terms. A useful tool to characterize these correlations is the torsional angle correlation function between two subsequent dihedral angles along the backbone. If is the -th torsional angle along the chain, we adopt the following definition for the function :


where is the average over single-chain configurations, and the subscript emphasizes the dependence of on the number of blobs. It is possible to define similar functions for any pair of torsional angles along the backbone, not necessarily nearest neighbors. Since we expect the correlation between such pairs to decay rather quickly with the chemical distance we limit our analysis to nearest neighbors. Note that with this definition, in the absence of correlation. In Fig. 3 we report the function for in both 4MB-2 (left panels) and FM (central panels) representation and the relative deviations (right panels). The 4MB-2 model predictions of are quite close to the polymer results, the deviations observed are at most in well localized regions of the plane, and are nearly independent on , confirming the good transferability of the potential, as observed for the other structural properties. This observation supports our implicit assumption that irreducible -body terms in the intramolecular potential for are negligible.

Figure 3: . Comparison between the 4MB-2 model (left column) and FM-CGR chains (central column) for (upper row), (central row) and (lower row) at . On the right column, for each value of , we show the deviations of the 4MB-2 results from the polymer ones,

It might appear surprising that we spent so much effort in modeling the intramolecular potential while the intermolecular interaction is treated in a very simple way. The effective interaction among chains, of blobs each, is a general function of variables. According to the multiblob ideas, at any finite reduced chain density , we can adopt the small density expansion and neglect many-chain interaction terms for large enough. However the effective intermolecular potential between two chains of blobs remains a function of variables. In our model we adopted several simplifying hypotheses: we replace the intermolecular chain-chain interaction by a sum of pair-wise potentials between blobs of different chains, we assume that all these pair potentials can be adequately approximated by a single central potential, and finally we approximate the potential with a single Gaussian function. It would not be too difficult to avoid the last two assumptions. It is instead crucial, from a practical point of view, to assume pairwise blob-blob interactions. This unique intermolecular blob-blob potential is optimized in such a way to reproduce the polymer center-of-mass radial distribution function (RDF) determined from full monomer simulations DAdamoG:2011p2477 (). This criterion ensures the correct thermodynamics for the CG model since the isothermal compressibility is related to the integral of the RDF by the compressibility equationHansenMcDonald1987 ()


where is the osmotic pressure of the solution and . An equally good choice as thermodynamics is concerned, would have been to reproduce the RDF between any pair of blobs of different chains since in Eq. (14), can be replaced by any other intermolecular RDF Akkermans:2001p1711 (). With our simplifying hypothesis of a unique blob-blob intermolecular interaction we are of course unable to perfectly reproduce RDFs between center-of-mass and between any pair of blobs of the two chains simultaneously. In figure 4 we show some of those RDFs predicted by the 4MB-2 model and we compare them with the corresponding quantities from FM simulation. We observe perfect agreement for the RDF between the centers of mass but, as expected, some small deviations occur for the other RDFs at distances , a range irrelevant to thermodynamics. These small deviations are the price to pay in order to have a much simplified model, still accurate as thermodynamics is concerned.

Figure 4: Radial distribution functions (RDF) between a pair of isolated chains (zero density) as a function of . Comparison between FM-CGR and 4MB-2 model predictions for (left panel) and (right panel). RDF between centers of mass (CM-CM), between end monomers (E-E), between central monomers (C-C) and between end and central monomers (E-C) are shown.
Figure 5: Intramolecular radial distribution function versus the reduced blob distance for model 4MB-2 with and for the FM-CGR chain with 20 blobs. Four values of the reduced density are reported: . For sake of clarity, results at different densities are shifted upward according to the legend.

iii.2 Model transferability with density

The second and more difficult requirement for a transferable MB model is the ability to reproduce FM predictions at any finite reduced chain density by increasing the number of blobs. To illustrate the ability of our model to match this second requirement we compare single-chain properties and the Equation of State (EOS) with FM predictions in the reduced density range . Note that for the present MB model, which is resolution invariant at zero density, the definition of the reduced polymer density is given in terms of an -independent characteristic length, the FM radius of gyration . This implies that at given values of , MB models with different are all at the same density of the underlying polymer model. Conversely, when adopting non-resolution invariant models Pierleoni:2007p193 (); DAdamoG:2011p2477 (), comparing systems with different at the same value of implies comparing systems at different values of the absolute polymer density which is inconsistent.

In Fig. 5 we compare FM-CGR and 4MB-2 predictions for the intramolecular radial distribution function , which is universal in terms of the reduced intra-blob distance , for a chain of blobs and for increasing density. The observed agreement is remarkable in particular at the highest density investigated, , deep inside the semidilute regime. The same kind of agreement is observed for the distribution of the radius of gyration shown in Fig. 6. The agreement worsens at the highest densities for shorter chains, but it remains very good for . We would like to emphasize the nontrivial nature of the observed agreement. Indeed, is a weighted average of the length distribution over all pairs of blobs, so the agreement for this quantity demonstrates the accuracy of the 4MB-2 model to reproduce pair distances well beyond four neighbors. includes information on pair distances but also on all possible cross-terms between different pairs, so the observed agreement for this quantity is even less obvious.

Figure 6: Chain radius-of-gyration probability distribution function ( is the radius of gyration of a single chain and its zero-density average) at for a chain of blobs. Comparison between the 4MB-2 model and the FM-CGR predictions.

Finally, we have computed the EOS, using the standard molecular virial route to the pressure for several values of Akkermans204 (). In Table 2 we report the compressibility factor , where is the osmotic pressure of the solution.

0.54 1.875(1) 1.874(2) 1.870(2) 1.85(1)
1.09 2.9760(4) 2.9784(6) 2.9819(1) 2.96(1)
2.18 5.6004(7) 5.708(1) 5.729(3) 5.63(1)
4.36 11.4541(5) 12.058(1) 12.269(3) 12.3(1)
6.54 17.532(2) 18.931(4) 19.515(5) 20.0(1)
8.72 23.751(1) 25.592(4) 26.989(8) 28.7(1)
Table 2: Compressibility factor for the 4MB-2 model for several values of and for polymers in the scaling limit Pelissetto:2008p1683 (), at various reduced densities from the dilute to the semidilute regime.

Up to density , results for all values of studied differ by less than 1% from the asymptotic polymer value. When density increases longer chains are systematically more accurate. At the relative deviation in is for , but for both and . At , the deviation is for , and decreases to for and for . This is expected from the multiblob heuristic argument since from Eq. (6) the reduced blob density,


decreases, at fixed , for increasing as long as . At , and which explains why increasingly large -values are necessary in order to neglect three- (or more)-blobs contacts as in the present model.

It is of some interest to investigate the relation between the blob radius of gyration of the coarse grained representation and the density-dependent correlation length characteristic of semidilute polymer solutions Doi (), also called correlation blob size. The latter quantity can be inferred from the single chain structure factor, , as the characteristic length at which a crossover from the ideal behavior ) to the excluded-volume behavior () is observed Farnoux78 (); Paul91 (); Muller00 (). In figure 7 we show at the Kratky plot, against , for the full monomer Domb-Joyce model (FM), the FM-CGR with and the 4MB-2 model with the same number of blobs. In the FM-CGR case is defined as above with monomer positions replaced by the blob center-of-mass positions .

Figure 7: (color online) Kratky plot at as a function of . Full monomer Domb-Joyce model (FM, red squares), the same chain model mapped onto a 30-blobs chain (FM-CGR, blue line), the model 4MB-2 with blobs (open green circles), the Debye function (dotted line). Note that has been used for the FM-CGR and 4MB-2 models. For the difference between between and is already small and will become even smaller for larger values of .

The correlation blob size can be readily determined from the FM data to be , which corresponds to since at this reduced density. On the other hand from Eq. (6) we obtain , and for , and , respectively, which are considerably smaller than . This result is consistent with the multiblob argument: from Eq. (15), . On the other hand in the semidilute regime Doi () . Substituting the former relation in the latter we get: . Therefore, as far as , implies in agreement with our findings.

In Fig. 7 we note that the CGR at such resolution reproduces FM data up to . If one wishes to have access to more local properties he should take a CGR with more blobs. In order to follow the crossover from ideal to excluded volume behavior in the CGR a quite larger number of blobs should be considered. Note that the large- behavior of the CGR () is universal since our blobs contain a sufficiently large number of monomers to be representative of the scaling limit. Therefore this large- behavior is inherent to the CGR and not to the particular MB model we considered, namely the 4MB-2 model. Finally, note the that perfect agreement between the 4MB-2 prediction and the FM-CGR data for the single chain structure factor is fully compatible with previous results on (see Fig. 5).

Iv Conclusions

In conclusion, we have developed a fully consistent, scale preserving, multiblob model for linear polymers in good solvent. The model is built by transferring to smaller length scales a tetramer model parametrized to reproduce a number of scalar correlations of the FM chain at zero density. We have shown that this model is fully consistent when varying the number of blobs at zero density, and more relevant, it is able to reproduce the universal EOS for athermal semidilute polymer solutions at high chain concentrations if a sufficiently large numbers of blobs per chain is chosen in such a way to always remain in the dilute blob regime. In particular, the present multiblob model with only 30 blobs provides the compressibility factor at with an accuracy of , a level of accuracy which would require the use of the order of thousand monomers even with the most efficient lattice model, namely the Domb-Joyce model tuned in such a way to cancel the leading-order corrections to scaling Pelissetto:2008p1683 (). The present CG strategy may be extended to more complex systems and thus opens the way to a quantitative study of semidilute polymer solutions in situations where local structure is important like, for instance, colloid-polymer mixtures, polymer-decorated colloidal systems and polymer brushes at relatively low grafting density. It would also be extremely interesting to develop an analogous multiblob strategy to describe the cross-over between the good solvent and the solvent conditions and finally to extend this approach to diblock copolymer solutions. Work in these directions is in progress.

As a final remark it is interesting to observe that our model represents, to a very good approximation, the first example of a truly ”fixed-point” model in the renormalization-group language, the elusive model which reproduces the scaling behavior at any level of coarse-graining. Previous work focused on defining models in which the leading scaling corrections were zero within errors. Here, instead, we obtain a model in which all scaling corrections are apparently small, and which is thus able to predict the scaling behavior for any . Of course, there is a price to pay: for any , we do not have access to all possible observables, but only to the large-scale properties that can be modelled by the chosen CG representation.

V Acknowledgement

We thank J.-P. Hansen for inspiring discussions. CP is supported by IIT under the SEED project grant n 259 SIMBEDD.

Vi Appendix: Explicit expressions of the tetramer potentials

In this appendix we provide full details on the potentials of the two models used in the present work, namely 4MB-1 and 4MB-2. Distance in the central potentials is .

The bonding potentials and have been parametrized for as


where is the value of the potential at the origin, is a constant which has been fixed so that the potential vanishes at the minimum, and we only consider the range .

Figure 8: Potentials for the tetramer 4MB-2 and 4MB-1. In the left column we show the pair potentials between first-neighbors along the chain (top) and between next-to-nearest neighbors (bottom), as a function of . In the central column we report the bending-angle potential, the potential acting on the as a function of (top), and the torsion-angle potential as a function of (in radians) for both the models 4MB-1,4MB-2 (bottom). In the right column we report the intermolecular pair potential (top) and the potential between the first and the last atom of the tetramer (bottom) as a function of . The symbols represent the numerical results obtained through IBI procedure and the lines are the corresponding interpolations.
Model Potential
4MB-1,4MB-2 0.8808 6.8935 1.8693 0.4287 5.3453 4.4274 0.0202
4MB-1,4MB-2 1.0392 3.9339 1.54603 0.7602 6.8442 3.3646 0.0194
4MB-1,4MB-2 3.9518 2.9224 1.1555 0.1321
4MB-1,4MB-2 0.2238 0.2284 0.0769 0.0600
4MB-1 0.1425 1.0698 4.5475 0.0391 0.14345
4MB-2 0.0797 0.970 0.990 -0.0251 0.0785
4MB-2 0.1593 0.0371 0.0158

Table 3: Numerical coefficients for the tetramer potentials () for 4MB-1 and 4MB-2 models

The potentials between next-to-nearest neighbors has been parametrized as ()


The bending potential is parametrized as


(the minimum of the potential corresponds to ), while the torsion potential is given by


The potential acting on the has been parametrized as:


Numerical values for the coefficients for are reported in Table 3 and the potential are displayed in figure 8.




  • (1) G. Voth, ed., Coarse-Graining of Condensed Phases and Biomolecular Systems (CRC Press, Boca Raton, 2009).
  • (2) Themed issue, Phys. Chem. Chem. Phys. 11, 1853 (2009).
  • (3) Themed issue, Soft Matter 5, 4341 (2009).
  • (4) Themed issue, Faraday Discussion 144, 1 (2010).
  • (5) P. Carbone, H. Varzaneh,X. Chen, F. J. Müller-Plathe, J. Chem. Phys. 128, 064904 (2008).
  • (6) D. Fritz, V. A. Harmandaris, K. Kremer, and N. F. A. van der Vegt, Macromolecules 42, 7579 (2009).
  • (7) C. Peter and K. Kremer, Soft Matter 5, 4357 (2009).
  • (8) W. Tschöp, K. Kremer, J. Batoulis, T. Bürger and O. Hahn, Acta Polymer. 49, 61 (1998).
  • (9) A. Lyubartsev, A. Mirzoev, L. Chen, and A. Laaksonen, Faraday Discuss. 144, 43 (2010).
  • (10) J.-P. Hansen, C. Addison, and A. A. Louis, J. Phys.: Condens. Matter 17, S3185 (2005).
  • (11) A. J. Clark and M. G. Guenza, J. Chem. Phys. 132, 044902 (2010).
  • (12) P. J. Flory and W. R. Krigbaum, J. Chem. Phys 18, 1086 (1950).
  • (13) A. Y. Grosberg, P. G. Khalatur, and A. R. Khokhlov, Makromol. Chem., Rapid Commun. 3, 709 (1982).
  • (14) J. Dautenhahn and C. K. Hall, Macromolecules 27, 5399 (1994).
  • (15) M. Murat and K. Kremer, J. Chem. Phys. 108, 4340 (1998).
  • (16) P. G. Bolhuis, A. A. Louis, J. P. Hansen, and E. J. Meijer, J. Chem. Phys. 114, 4296 (2001).
  • (17) A. Pelissetto and J.-P. Hansen, J. Chem. Phys. 122, 134904 (2005).
  • (18) G. D’Adamo, A. Pelissetto and C. Pierleoni, “Polymers as soft compressible spheres”, J. Chem. Phys. in press (2012),
  • (19) I. Coluzza and J.-P. Hansen, Phys. Rev. Lett. 100, 016104 (2008).
  • (20) C. Pierleoni, C. I. Addison, J. P. Hansen, and V. Krakoviack, Phys. Rev. Lett. 96, 128302 (2006).
  • (21) J.-P. Hansen and C. Pearson, Mol. Phys. 104, 3389 (2006).
  • (22) E. Sambriski, M. Guenza, Phys. Rev. E 76, 051801 (2007).
  • (23) B. Capone, C. Pierleoni, J.-P. Hansen, and V. Krakoviack, J. Phys. Chem. B 113 , 3629 (2009).
  • (24) C. Gross,W. Paul, Soft Matter 6, 3273 (2010).
  • (25) I. Coluzza, B. Capone and J.-P. Hansen, Soft Matter 7, 5255 (2011).
  • (26) B. Capone, J.-P. Hansen and I. Coluzza, J. of Phys.: Cond. Mat. 23, 194102 (2011).
  • (27) M. Doi, Introduction to Polymer Physics (Clarendon Press, Oxford, 1992).
  • (28) C. Pierleoni, B. Capone and J. P. Hansen, J. Chem. Phys. 127, 171102 (2007).
  • (29) A. Pelissetto, J. Phys.: Condens. Matter 21, 115108 (2009).
  • (30) T. Vettorel, G. Besold, and K. Kremer, Soft Matter 6, 2282 (2010).
  • (31) G. D’Adamo, A. Pelissetto, and C. Pierleoni, Soft Matter 8, 5151 (2012).
  • (32) N. Clisby, Phys. Rev. Lett. 104, 55702 (2010).
  • (33) M. Laso, H. Öttinger, and U. Suter, J. Chem. Phys. 95, 2178 (1991).
  • (34) To study polymer chains in the scaling limit, we adopted the Domb-Joyce lattice model with , the value for which the coefficient of the leading-order corrections to scaling approximatively vanishes Pelissetto:2008p1683 (). As discussed in the appendix of Ref. DAdamoG:2011p2477 (), we considered chains of lengths and 2400 to extrapolate at the infinite chain length limit.
  • (35) A. Pelissetto, J. Chem. Phys. 129, 044901 (2008).
  • (36) S. Caracciolo, B. M. Mognetti, and A. Pelissetto, J. Chem. Phys. 125, 094903 (2006).
  • (37) J. P. Hansen and I. McDonald, Theory of Simple Liquids, 3rd ed. (Academic Press, Amsterdam, 2006)
  • (38) R. L. C. Akkermans and W. J. Briels, J. Chem. Phys. 114, 1020 (2001).
  • (39) R. L. C. Akkermans and G. Ciccotti, J. Phys. Chem. B 108, 6866 (2004).
  • (40) B. Farnoux, F. Boue, J. P. Cotton, M. Daoud, G. Jannink, M. Nierlich and P. G. de Gennes, J. Phys. (Paris) 39, 77 (1978).
  • (41) W. Paul, K. Binder, D.W. Hermann and K. Kremer, J. Phys. II (Paris) ÿ1, 37 (1991).
  • (42) M. Müller, K. Binder and L. Schäfer, Macromolecules 33, 4568 (2000).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description