Coupled-decoupled configurational-bias Monte Carlo (CD-CBMC) simulations in the Gibbs ensemble were carried out to determine the vapor-liquid coexistence curves for n-triacontane and 2,6,10,15,19,23-hexamethyltetracosane (squalane). The transferable potentials for phase equilibria-united atom (TraPPE-UA) force field was used for these simulations. The simulated systems consisted of 200 molecules and the production periods extended to 100,000 Monte Carlo cycles, a system size about twice as large and a simulation length about one order of magnitude longer than used in previous simulations. The simulation results are in satisfactory agreement with the available experimental data. Examination of the coexistence curves in reduced units for the two triacontane isomers and for n-octane and 2,5-dimethylhexane shows that both molecular weight and branching can lead to deviations from the principle of corresponding states. Analysis of the molecular structures in the vapor and liquid phases points to a partial collapse (self-solvation) of the triacontane isomers as the likely origin of the deviations from the principle of corresponding states.
Bibliographical noteFunding Information:
Financial support from the National Science Foundation (CTS-9813601 and CHE-9816328), a Lando/NSF-REU (CHE-942420301) Summer Research Fellowship (NDZ), and a Department of Energy Computational Science Graduate Fellowship (MGM) is gratefully acknowledged. Part of the computer resources were provided by the Minnesota Supercomputing Institute.
- Critical properties
- Isomeric alkanes
- Molecular simulation
- Principle of corresponding states
- Vapor-liquid equilibria