Aggregation in dilute aqueous tert-butyl alcohol solutions: insights from large-scale simulations

Molecular dynamics simulations employing up to 64,000 particles are used to investigate aggregation and microheterogeneity in aqueous tert-butyl alcohol (TBA) solutions for TBA mole fractions X(t) ≤ 0.1. Four different force fields are considered. It is shown that the results obtained can be strongly dependent on the particular force field employed, and can be significantly influenced by system size. Two of the force fields considered show TBA aggregation in the concentration range X(t) ≈ 0.03 - 0.06. For these models, systems of 64,000 particles are minimally sufficient to accommodate the TBA aggregates. The structures resulting from TBA aggregation do not have a well-defined size and shape, as one might find in micellar systems, but are better described as TBA-rich and water-rich regions. All pair correlation functions exhibit long-range oscillatory behavior with wavelengths that are much larger than molecular length scales. The oscillations are not strongly damped and the correlations can easily exceed the size of the simulation cell, even for the low TBA concentrations considered here. We note that these long-range correlations pose a serious problem if one wishes to obtain certain physical properties such as Kirkwood-Buff integrals from simulation results. In contrast, two other force fields that we consider show little sign of aggregation for X(t) ≲ 0.08. In our 64,000 particle simulations all four models considered show demixing-like behavior for X(t) ≳ 0.1, although such behavior is not evident in smaller systems of 2000 particles. The meaning of the demixing-like behavior is unclear. Since real TBA-water solutions do not demix, it might be an indication that all four models we consider poorly represent the real system. Alternatively, it might be an artifact of finite system size. Possibly, the apparent demixing indicates that for X(t) ≳ 0.1, the stable TBA aggregates are simply too large to fit into the simulation cell. Our results provide a view of the possible nature of microheterogeneity in dilute TBA-water solutions, and of the associated long correlation lengths. It is clear that system size can be a very important factor in simulations of these solutions, and must be taken into account in the evaluation and development of TBA-water force fields.