Abstract

We investigate driving forces of the liquid–liquid phase separation of N-isopropylpropionamide (NiPPA) aqueous solutions above the lower critical solution temperature using molecular dynamics simulations. Spontaneous phase separations of the model aqueous solution with a modified OPLS-AA force field are observed above the experimentally determined cloud point. The destabilization toward the phase separation is confirmed by temperature dependence of the long-wavelength limit of the concentration-concentration structure factor, the dominant component of which is found to be an increasing effective attraction between NiPPA molecules. At varying temperatures, the potentials of mean force (PMFs) between a pair of NiPPA molecules at infinite dilution are obtained and decomposed into the nonpolar and Coulombic contributions. The nonpolar contribution, arising essentially from molecular volume, promotes association of NiPPA molecules with increasing temperature while the Coulombic one antagonizes the association. Thus, our analysis leads to a conclusion that the driving force of thermally induced aggregation of NiPPA molecules is the temperature dependence of the nonpolar contribution in PMF between NiPPA molecules, not the temperature dependence of the number or strength of hydrogen bonds between NiPPA and water molecules.

Highlights

  • The intermolecular interaction of water is of the TIP4P/2005 model[52], which provides the most realistic description of the bulk liquid density and the excess chemical potential of simple molecules[53]

  • The model of NiPPA is based on the all-atom optimized parameters for liquid simulations (OPLS-AA) force field[54], but with all the original electrostatic charges on the interaction sites in NiPPA increased by a factor of χ = 1.3149

  • The Lennard–Jones parameters of the cross-interaction are computed by εij = εiiεjj and σij = σiiσ jj

Read more

Summary

Introduction

MD simulations in the canonical (NVT) and the isothermal–isobaric (NPT) ensembles are performed using the GROMACS 4.6.6. The time step of the MD simulations is 1.0 fs. The pressure of 0.1 MPa and the temperatures ranging from 280 K to 360 K are controlled by the Nośe–Hoover thermostat[57,58] and the Parrinello–Rahman barostat[59], respectively, whereas the Berendsen algorithm[60] is used for equilibration. Periodic boundary conditions are applied for all directions of the simulation box. The lengths of equilibration and production runs depend on the system and the kind of analysis; details are given below and in Results and Discussion

Methods
Results
Conclusion
Full Text
Paper version not known

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.