Molecular dynamics simulation of diffusion coefficients and structural properties of some alkylbenzenes in supercritical carbon dioxide at infinite dilution

The binary infinite dilute diffusion coefficients, D12∞, of some alkylbenzenes (Ph-Cn, from Ph-H to Ph-C12) from 313 K to 333 K at 15 MPa in supercritical carbon dioxide (scCO2) have been studied by molecular dynamics (MD) simulation. The MD values agree well with the experimental ones, which indicate MD simulation technique is a powerful way to predict and obtain diffusion coefficients of solutes in supercritical fluids. Besides, the local structures of Ph-Cn/CO2 fluids are further investigated by calculating radial distribution functions and coordination numbers. It qualitatively convinces that the first solvation shell of Ph-Cn in scCO2 is significantly influenced by the structure of Ph-Cn solute. Meanwhile, the mean end-to-end distance, the mean radius of gyration and dihedral angle distribution are calculated to gain an insight into the structural properties of Ph-Cn in scCO2. The abnormal trends of radial distribution functions and coordination numbers can be reasonably explained in term of molecular flexibility. Moreover, the computed results of dihedral angle clarify that flexibility of long-chain Ph-Cn is the result of internal rotation of C-C single bond (σc-c) in alkyl chain. It is interesting that compared with n-alkane, because of the existence of benzene ring, the flexibility of alkyl chain in Ph-Cn with same carbon atom number is significantly reduced, as a result, the carbon chain dependence of diffusion behaviors for long-chain n-alkane (n ≥ 5) and long-chain Ph-Cn (n ≥ 4) in scCO2 are different