324434 The Importance of Intra- and Inter-Molecular Interactions in the Development of Coarse-Grained Models for Chain Fluids Using SAFT-γ Mie in Molecular Simulations

Friday, November 8, 2013: 9:42 AM
Union Square 16 (Hilton)
Sadia Rahman1, Carlos Braga1, Olga Lobanova1, Vasileios Raptis1, Amparo Galindo2, George Jackson2 and Erich A. Muller1, (1)Chemical Engineering, Imperial College London, London, United Kingdom, (2)Centre for Process Systems Engineering, Department of Chemical Engineering, Imperial College London, London, United Kingdom


In common coarse-graining (CG) approaches one starts from an atomistic, bottom-up, description of molecules and proceed to integrate degrees of freedom, simplifying the force field in order to speed up the computer simulation of complex molecules and mixtures. [1-4] By contrast, the CG intermolecular potentials can also be obtained with a top-down approach by reproducing the macroscopic configurational properties of the system with a molecular-based equation of state. The latter approach is considered here, in particular, the use of the SAFT-γ Mie version of the Statistical Associating Fluid Theory (SAFT) which is based on interactions of the Mie form and has successfully been applied to simple organic molecules [5-7]. In our current work we explore the role of inter- and intra-molecular interactions in SAFT-γ Mie force fields of the linear alkane homologous series.

Alkane Study

Despite its success in correctly reproducing the thermodynamic properties of complex fluids, SAFT based approaches do not take into account the structure and flexibility of molecular systems explicitly, apart from the assumption that chains of m segments are formed [8]. We have addressed this issue by combining the SAFT-γ Mie top-down approach with a parameterization of the CG intramolecular (bonding) potentials on the basis of atomistic simulations.

As prototypical molecules we first perform the analysis for long chain n-alkanes. In order to obtain a CG model, three consecutive backbone carbon atoms and their bonded hydrogen's are all grouped together to form a "super-atom", or CG moiety.  The inter-molecular parameterization of the CG beads with the SAFT-γ Mie methodology then leads to two possibilities, depending on the recognition (or not) of the distinct identity of the terminal end CH3 group. In the first model (homonuclear) the beads representing the end and intervening CG alkyl groups are taken as identical; in the second model (heteronuclear) a distinction is made between the interactions of the two types of CG beads. (Fig 1)

Description: Macintosh HD:Users:sr1508:Desktop:Untitled.gif

Figure 1 Schematic of the CG strategy for n-nonane, lumping three backbone carbon atoms into single beads. For the homonuclear model (left) all beads are treated on an equal footing, and for the heteronuclear model (right) the CH2CH2CH3 beads that contain the CH3 end group and the CH2CH2CH2 beads that make up the intervening groups are distinguished.

To address the importance of chain flexibility, we have performed single component molecular dynamics simulations of an equivalent united-atom (UA) model using the NERD force field [9]. From the set of atomistic trajectories the CG intermolecular force field was computed using the potential of mean force of the corresponding CG degrees of freedom, e.g., CG bond distances and angles. The CG bond stretching and bond angle bending interactions were then superimposed to the CG intermolecular parameters. The full model includes both inter- and intra-molecular contributions.

An assessment of the different contributions to the full intermolecular potential is made by comparing the results obtained from the molecular dynamics simulations of fluid phase equilibria as well as the interfacial, transport and structural properties. An example of how the structural properties of nonane are affected when using different intra- and inter-molecular parameters can be seen in Figure 2. The g(r) results obtained for the homonuclear model was the same as the heteronuclear, so has not been included to keep the figure simple.

Description: gofr_nonane

Figure 2: Centre-to-centre radial distribution function for n-nonane comparing the heteronuclear SAFT-γ Mie CG descriptions (cf. Figure 1) with completely flexible (pearl-necklace) model or the restrained bonding angle model.

3. Conclusion

The inclusion of intra- molecular interactions CG models allows for an improved description of thermodynamic and structural properties without a detriment to the speed of the calculations. An additional advantage of the heteronuclear description is that they provide fully transferable CG models for linear alkanes.


[1] K. A. Maerzke and J. I. Siepmann, J. Phys. Chem. B, 115, 3452 (2011).

[2] S. J. Marrink, H. J. Risselada, S. Yefimov, D. P. Tieleman, and A. H. de Vries, J. Phys. Chem. B, 111, 7812 (2007).

[3] S. O. Nielsen, C. F. Lopez, G. Srinivas, and M. L. Klein, J. Chem. Phys., 119, 7043 (2003).

[4] W. Shinoda, R. DeVane, and M. L. Klein, Mol. Simul. 33, 27 (2007).

[5] C. Avendano, T. Lafitte,  A. Galindo, C. S. Adjiman, G. Jackson and E. A. Muller, J. Phys. Chem. B, 115, 11154 (2011).

[6] T. Lafitte, C. Avendano, V. Papaioannou, A. Galindo, C. S. Adjiman ,G. Jackson and E. A. Muller, Mol. Phys, 110, 1189 (2012).

[7] C. Avendano, T. Lafitte, C. S. Adjiman, A. Galindo, E. A. Muller and G. Jackson, J. Phys. Chem. B, 117, 2717 (2013).

[8] E.A. Muller and K.E. Gubbins, Mol. Phys., 80, 957 (1993).

[9] S.K.Nath, B.J. Banaszak & J.J. de Pablo, J. Chem. Phys, 114, 3612 (2001)                                                                                                                                                                                                                                                                                                                                                                                                                  

Extended Abstract: File Uploaded
See more of this Session: Development of Intermolecular Potential Models
See more of this Group/Topical: Engineering Sciences and Fundamentals